Rapid evolution of recombination landscapes during the divergence of cichlid ecotypes in Lake Masoko

Variation of recombination rate along the genome is of crucial importance to rapid adaptation and organismal diversification. Many unknowns remain regarding how and why recombination landscapes evolve in nature. Here, we reconstruct recombination maps based on linkage disequilibrium and use subsampling and simulations to derive a new measure of recombination landscape evolution: the Population Recombination Divergence Index (PRDI). Using PRDI, we show that fine-scale recombination landscapes differ substantially between two cichlid fish ecotypes of Astatotilapia calliptera that diverged only ~2,500 generations ago. Perhaps surprisingly, recombination landscape differences are not driven by divergence in terms of allele frequency (FST) and nucleotide diversity (Δ(π)): although there is some association, we observe positive PRDI in regions where FST and Δ(π) are zero. We found a stronger association between the evolution of recombination and 47 large haplotype blocks that are polymorphic in Lake Masoko, cover 21% of the genome, and appear to include multiple inversions. Among haplotype blocks, there is a strong and clear association between the degree of recombination divergence and differences between ecotypes in heterozygosity, consistent with recombination suppression in heterozygotes. Overall, our work provides a holistic view of changes in population recombination landscapes during the early stages of speciation with gene flow.