Connecting theory and data to understand recombination rate evolution

Natural selection on apical membrane antigen 1 (AMA1) of an emerging zoonotic malaria parasite Plasmodium inui

Apical membrane antigen 1 (AMA1) of malaria parasites plays an important role in host cell invasion. Antibodies to AMA1 can inhibit malaria merozoite invasion of erythrocytes while vaccine-induced specific cytotoxic T cell responses to this protein are associated with clinical protection. Polymorphisms in AMA1 of Plasmodium falciparum (PfAMA1) and P. vivax (PvAMA1) are of concern for vaccine development. To date, little is known about sequence diversity in ama1 of P. inui (Piama1), an emerging zoonotic malaria parasite. In this study, 80 complete Piama1 coding sequences were obtained from 57 macaques in Thailand that defined 60 haplotypes clustering in two phylogenetic lineages. In total, 74 nucleotide substitutions were identified and distributed unevenly across the gene. Blockwise analysis of the rates of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions did not show a significant deviation from neutrality among Thai isolates. However, significantly negative Tajima's D values were detected in domain I and the loop region of domain II, implying purifying selection. Codon-based analysis of dN/dS has identified 12 and 14 codons under positive and negative selections, respectively. Meanwhile, 85 amino acid substitutions were identified among 80 Thai and 11 non-Thai PiAMA1 sequences. Of these, 48 substituted residues had a significant alteration in physicochemical properties, suggesting positive selection. More than half of these positively selected amino acids (32 of 48) corresponded to the predicted B-cell or T-cell epitopes, suggesting that selective pressure could be mediated by host immunity. Importantly, 14 amino acid substitutions were singletons and predicted to be deleterious that could be subject to ongoing purifying selection or elimination. Besides genetic drift and natural selection, intragenic recombination identified in domain II could generate sequence variation in Piama1. It is likely that malarial ama1 exhibits interspecies differences in evolutionary histories. Knowledge of the sequence diversity of the Piama1 locus further provides an evolutionary perspective of this important malaria vaccine candidate.

Evolution of unexpected diversity in a putative mating type locus and its correlation with genome variability reveals likely asexuality in the model mycorrhizal fungus Rhizophagus irregularis

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) form mutualistic partnerships with approximately 80% of plant species. AMF, and their diversity, play a fundamental role in plant growth, driving plant diversity, and global carbon cycles. Knowing whether AMF are sexual or asexual has fundamental consequences for how they can be used in agricultural applications. Evidence for and against sexuality in the model AMF, Rhizophagus irregularis, has been proposed. The discovery of a putative mating-type locus (MAT locus) in R. irregularis, and the previously suggested recombination among nuclei of a dikaryon R. irregularis isolate, potentially suggested sexuality. Unless undergoing frequent sexual reproduction, evolution of MAT-locus diversity is expected to be very low. Additionally, in sexual species, MAT-locus evolution is decoupled from the evolution of arbitrary genome-wide loci.

RESULTS

We studied MAT-locus diversity of R. irregularis. This was then compared to diversification in a phosphate transporter gene (PTG), that is not involved in sex, and to genome-wide divergence, defined by 47,378 single nucleotide polymorphisms. Strikingly, we found unexpectedly high MAT-locus diversity indicating that either it is not involved in sex, or that AMF are highly active in sex. However, a strongly congruent evolutionary history of the MAT-locus, PTG and genome-wide arbitrary loci allows us to reject both the hypothesis that the MAT-locus is involved in mating and that the R. irregularis lineage is sexual.

CONCLUSION

Our finding shapes the approach to developing more effective AMF strains and is highly informative as it suggests that introduced strains applied in agriculture will not exchange DNA with native populations.

Creating Meiotic Recombination-Regulating DNA Sites by SpEDIT in Fission Yeast Reveals Inefficiencies, Target-Site Duplications, and Ectopic Insertions

Recombination hotspot-activating DNA sites (e.g., M26, CCAAT, Oligo-C) and their binding proteins (e.g., Atf1-Pcr1 heterodimer; Php2-Php3-Php5 complex, Rst2, Prdm9) regulate the distribution of Spo11 (Rec12)-initiated meiotic recombination. We sought to create 14 different candidate regulatory DNA sites via bp substitutions in the ade6 gene of Schizosaccharomyces pombe. We used a fission yeast-optimized CRISPR-Cas9 system (SpEDIT) and 196 bp-long dsDNA templates with centrally located bp substitutions designed to ablate the genomic PAM site, create specific 15 bp-long DNA sequences, and introduce a stop codon. After co-transformation with a plasmid that encoded both the guide RNA and Cas9 enzyme, about one-third of colonies had a phenotype diagnostic for DNA sequence changes at ade6. PCR diagnostics and DNA sequencing revealed a diverse collection of alterations at the target locus, including: (A) complete or (B) partial template-directed substitutions; (C) non-homologous end joinings; (D) duplications; (E) bp mutations, and (F) insertions of ectopic DNA. We concluded that SpEDIT can be used successfully to generate a diverse collection of DNA sequence elements within a reporter gene of interest. However, its utility is complicated by low efficiency, incomplete template-directed repair events, and undesired alterations to the target locus.

Divergence time shapes gene reuse during repeated adaptation

When diverse lineages repeatedly adapt to similar environmental challenges, the extent to which the same genes are involved (gene reuse) varies across systems. We propose that divergence time among lineages is a key factor driving this variability: as lineages diverge, the extent of gene reuse should decrease due to reductions in allele sharing, functional differentiation among genes, and restructuring of genome architecture. Indeed, we show that many genomic studies of repeated adaptation find that more recently diverged lineages exhibit higher gene reuse during repeated adaptation, but the relationship becomes less clear at older divergence time scales. Thus, future research should explore the factors shaping gene reuse and their interplay across broad divergence time scales for a deeper understanding of evolutionary repeatability.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

rec-1 loss of function increases recombination in the central gene clusters at the expense of autosomal pairing centers

Meiotic control of crossover (CO) number and position is critical for homologous chromosome segregation and organismal fertility, recombination of parental genotypes, and the generation of novel genetic combinations. We here characterize the recombination rate landscape of a rec-1 loss of function modifier of CO position in Caenorhabditis elegans, one of the first ever modifiers discovered. By averaging CO position across hermaphrodite and male meioses and by genotyping 203 single-nucleotide variants covering about 95% of the genome, we find that the characteristic chromosomal arm-center recombination rate domain structure is lost in the loss of function rec-1 mutant. The rec-1 loss of function mutant smooths the recombination rate landscape but is insufficient to eliminate the nonuniform position of CO. Lower recombination rates in the rec-1 mutant are particularly found in the autosomal arm domains containing the pairing centers. We further find that the rec-1 mutant is of little consequence for organismal fertility and egg viability and thus for rates of autosomal nondisjunction. It nonetheless increases X chromosome nondisjunction rates and thus male appearance. Our findings question the maintenance of recombination rate heritability and genetic diversity among C. elegans natural populations, and they further suggest that manipulating genetic modifiers of CO position will help find quantitative trait loci located in low-recombining genomic regions normally refractory to discovery.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Modeling Recombination Rate as a Quantitative Trait Reveals New Insight into Selection in Humans

Meiotic recombination is both a fundamental biological process required for proper chromosomal segregation during meiosis and an important genomic parameter that shapes major features of the genomic landscape. However, despite the central importance of this phenotype, we lack a clear understanding of the selective pressures that shape its variation in natural populations, including humans. While there is strong evidence of fitness costs of low rates of recombination, the possible fitness costs of high rates of recombination are less defined. To determine whether a single lower fitness bound can explain the variation in recombination rates observed in human populations, we simulated the evolution of recombination rates as a sexually dimorphic quantitative trait. Under each scenario, we statistically compared the resulting trait distribution with the observed distribution of recombination rates from a published study of the Icelandic population. To capture the genetic architecture of recombination rates in humans, we modeled it as a moderately complex trait with modest heritability. For our fitness function, we implemented a hyperbolic tangent curve with several flexible parameters to capture a wide range of existing hypotheses. We found that costs of low rates of recombination alone are likely insufficient to explain the current variation in recombination rates in both males and females, supporting the existence of fitness costs of high rates of recombination in humans. With simulations using both upper and lower fitness boundaries, we describe a parameter space for the costs of high recombination rates that produces results consistent with empirical observations.

The fine-scale recombination rate variation and associations with genomic features in a butterfly

Recombination is a key molecular mechanism that has profound implications on both micro- and macroevolutionary processes. However, the determinants of recombination rate variation in holocentric organisms are poorly understood, in particular in Lepidoptera (moths and butterflies). The wood white butterfly (Leptidea sinapis) shows considerable intraspecific variation in chromosome numbers and is a suitable system for studying regional recombination rate variation and its potential molecular underpinnings. Here, we developed a large whole-genome resequencing data set from a population of wood whites to obtain high-resolution recombination maps using linkage disequilibrium information. The analyses revealed that larger chromosomes had a bimodal recombination landscape, potentially caused by interference between simultaneous chiasmata. The recombination rate was significantly lower in subtelomeric regions, with exceptions associated with segregating chromosome rearrangements, showing that fissions and fusions can have considerable effects on the recombination landscape. There was no association between the inferred recombination rate and base composition, supporting a limited influence of GC-biased gene conversion in butterflies. We found significant but variable associations between the recombination rate and the density of different classes of transposable elements, most notably a significant enrichment of short interspersed nucleotide elements in genomic regions with higher recombination rate. Finally, the analyses unveiled significant enrichment of genes involved in farnesyltranstransferase activity in recombination coldspots, potentially indicating that expression of transferases can inhibit formation of chiasmata during meiotic division. Our results provide novel information about recombination rate variation in holocentric organisms and have particular implications for forthcoming research in population genetics, molecular/genome evolution, and speciation.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

The effect of DNA polymorphisms and natural variation on crossover hotspot activity in Arabidopsis hybrids

In hybrid organisms, genetically divergent homologous chromosomes pair and recombine during meiosis; however, the effect of specific types of polymorphisms on crossover is poorly understood. Here, to analyze this in Arabidopsis, we develop the seed-typing method that enables the massively parallel fine-mapping of crossovers by sequencing. We show that structural variants, observed in one of the generated intervals, do not change crossover frequency unless they are located directly within crossover hotspots. Both natural and Cas9-induced deletions result in lower hotspot activity but are not compensated by increases in immediately adjacent hotspots. To examine the effect of single nucleotide polymorphisms on crossover formation, we analyze hotspot activity in mismatch detection-deficient msh2 mutants. Surprisingly, polymorphic hotspots show reduced activity in msh2. In lines where only the hotspot-containing interval is heterozygous, crossover numbers increase above those in the inbred (homozygous). We conclude that MSH2 shapes crossover distribution by stimulating hotspot activity at polymorphic regions.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Gene flow biases population genetic inference of recombination rate

Accurate estimates of the rate of recombination are key to understanding a host of evolutionary processes as well as the evolution of the recombination rate itself. Model-based population genetic methods that infer recombination rates from patterns of linkage disequilibrium in the genome have become a popular method to estimate rates of recombination. However, these linkage disequilibrium-based methods make a variety of simplifying assumptions about the populations of interest that are often not met in natural populations. One such assumption is the absence of gene flow from other populations. Here, we use forward-time population genetic simulations of isolation-with-migration scenarios to explore how gene flow affects the accuracy of linkage disequilibrium-based estimators of recombination rate. We find that moderate levels of gene flow can result in either the overestimation or underestimation of recombination rates by up to 20–50% depending on the timing of divergence. We also find that these biases can affect the detection of interpopulation differences in recombination rate, causing both false positives and false negatives depending on the scenario. We discuss future possibilities for mitigating these biases and recommend that investigators exercise caution and confirm that their study populations meet assumptions before deploying these methods.

On the prospect of achieving accurate joint estimation of selection with population history

As both natural selection and population history can affect genome-wide patterns of variation, disentangling the contributions of each has remained as a major challenge in population genetics. We here discuss historical and recent progress towards this goal-highlighting theoretical and computational challenges that remain to be addressed, as well as inherent difficulties in dealing with model complexity and model violations-and offer thoughts on potentially fruitful next steps.

Adaptive Divergence of Meiotic Recombination Rate in Ecological Speciation

Theories predict that directional selection during adaptation to a novel habitat results in elevated meiotic recombination rate. Yet the lack of population-level recombination rate data leaves this hypothesis untested in natural populations. Here, we examine the population-level recombination rate variation in two incipient ecological species, the microcrustacean Daphnia pulex (an ephemeral-pond species) and Daphnia pulicaria (a permanent-lake species). The divergence of D. pulicaria from D. pulex involved habitat shifts from pond to lake habitats as well as strong local adaptation due to directional selection. Using a novel single-sperm genotyping approach, we estimated the male-specific recombination rate of two linkage groups in multiple populations of each species in common garden experiments and identified a significantly elevated recombination rate in D. pulicaria. Most importantly, population genetic analyses show that the divergence in recombination rate between these two species is most likely due to divergent selection in distinct ecological habitats rather than neutral evolution.

Sex-specific variation in the genome-wide recombination rate

In most species that reproduce sexually, successful gametogenesis requires recombination during meiosis. The number and placement of crossovers (COs) vary among individuals, with females and males often presenting the most striking contrasts. Despite the recognition that the sexes recombine at different rates (heterochiasmy), existing data fail to answer the question of whether patterns of genetic variation in recombination rate are similar in the two sexes. To fill this gap, we measured the genome-wide recombination rate in both sexes from a panel of wild-derived inbred strains from multiple subspecies of house mice (Mus musculus) and from a few additional species of Mus. To directly compare recombination rates in females and males from the same genetic backgrounds, we applied established methods based on immunolocalization of recombination proteins to inbred strains. Our results reveal discordant patterns of genetic variation in the two sexes. Whereas male genome-wide recombination rates vary substantially among strains, female recombination rates measured in the same strains are more static. The direction of heterochiasmy varies within two subspecies, Mus musculus molossinus and Mus musculus musculus. The direction of sex differences in the length of the synaptonemal complex and CO positions is consistent across strains and does not track sex differences in genome-wide recombination rate. In males, contrasts between strains with high recombination rate and strains with low recombination rate suggest more recombination is associated with stronger CO interference and more double-strand breaks. The sex-specific patterns of genetic variation we report underscore the importance of incorporating sex differences into recombination research.

Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C. obscurior that could contribute to the species' adaptive potential.

Evolution and Plasticity of Genome-Wide Meiotic Recombination Rates

Sex, as well as meiotic recombination between homologous chromosomes, is nearly ubiquitous among eukaryotes. In those species that use it, recombination is important for chromosome segregation during gamete production, and thus for fertility. Strikingly, although in most species only one crossover event per chromosome is required to ensure proper segregation, recombination rates vary considerably above this minimum and show variation within and among species. However, whether this variation in recombination is adaptive or neutral and what might shape it remain unclear. Empirical studies and theory support the idea that recombination is generally beneficial but can also have costs. Here, we review variation in genome-wide recombination rates, explore what might cause this, and discuss what is known about its mechanistic basis. We end by discussing the environmental sensitivity of meiosis and recombination rates, how these features may relate to adaptation, and their implications for a broader understanding of recombination rate evolution. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Male meiotic recombination rate varies with seasonal temperature fluctuations in wild populations of autotetraploid Arabidopsis arenosa

Meiosis, the cell division by which eukaryotes produce haploid gametes, is essential for fertility in sexually reproducing species. This process is sensitive to temperature, and can fail outright at temperature extremes. At less extreme values, temperature affects the genome‐wide rate of homologous recombination, which has important implications for evolution and population genetics. Numerous studies in laboratory conditions have shown that recombination rate plasticity is common, perhaps nearly universal, among eukaryotes. These studies have also shown that variation in the length or timing of stresses can strongly affect results, raising the important question whether these findings translate to more variable field conditions. Moreover, lower or higher recombination rate could cause certain kinds of meiotic aberrations, especially in polyploid species—raising the additional question whether temperature fluctuations in field conditions cause problems. Here, we tested whether (1) recombination rate varies across a season in the wild in two natural populations of autotetraploid Arabidopsis arenosa, (2) whether recombination rate correlates with temperature fluctuations in nature, and (3) whether natural temperature fluctuations might cause meiotic aberrations. We found that plants in two genetically distinct populations showed a similar plastic response with recombination rate increases correlated with both high and low temperatures. In addition, increased recombination rate correlated with increased multivalent formation, especially at lower temperature, hinting that polyploids in particular may suffer meiotic problems in conditions they encounter in nature. Our results show that studies of recombination rate plasticity done in laboratory settings inform our understanding of what happens in nature.

Genetic variation in recombination rate in the pig

Background

Meiotic recombination results in the exchange of genetic material between homologous chromosomes. Recombination rate varies between different parts of the genome, between individuals, and is influenced by genetics. In this paper, we assessed the genetic variation in recombination rate along the genome and between individuals in the pig using multilocus iterative peeling on 150,000 individuals across nine genotyped pedigrees. We used these data to estimate the heritability of recombination and perform a genome-wide association study of recombination in the pig.

Results

Our results confirmed known features of the recombination landscape of the pig genome, including differences in genetic length of chromosomes and marked sex differences. The recombination landscape was repeatable between lines, but at the same time, there were differences in average autosome-wide recombination rate between lines. The heritability of autosome-wide recombination rate was low but not zero (on average 0.07 for females and 0.05 for males). We found six genomic regions that are associated with recombination rate, among which five harbour known candidate genes involved in recombination: RNF212, SHOC1, SYCP2, MSH4 and HFM1.

Conclusions

Our results on the variation in recombination rate in the pig genome agree with those reported for other vertebrates, with a low but nonzero heritability, and the identification of a major quantitative trait locus for recombination rate that is homologous to that detected in several other species. This work also highlights the utility of using large-scale livestock data to understand biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00643-0.

Meiosis and beyond - understanding the mechanistic and evolutionary processes shaping the germline genome

The separation of germ cell populations from the soma is part of the evolutionary transition to multicellularity. Only genetic information present in the germ cells will be inherited by future generations, and any molecular processes affecting the germline genome are therefore likely to be passed on. Despite its prevalence across taxonomic kingdoms, we are only starting to understand details of the underlying micro-evolutionary processes occurring at the germline genome level. These include segregation, recombination, mutation and selection and can occur at any stage during germline differentiation and mitotic germline proliferation to meiosis and post-meiotic gamete maturation. Selection acting on germ cells at any stage from the diploid germ cell to the haploid gametes may cause significant deviations from Mendelian inheritance and may be more widespread than previously assumed. The mechanisms that affect and potentially alter the genomic sequence and allele frequencies in the germline are pivotal to our understanding of heritability. With the rise of new sequencing technologies, we are now able to address some of these unanswered questions. In this review, we comment on the most recent developments in this field and identify current gaps in our knowledge.

A simple expression for the strength of selection on recombination generated by interference among mutations

One of the most widely cited hypotheses to explain the evolutionary maintenance of genetic recombination states that the reshuffling of genotypes at meiosis increases the efficiency of natural selection by reducing interference among selected loci. However, and despite several decades of theoretical work, a quantitative estimation of the possible selective advantage of a mutant allele increasing chromosomal map length (the average number of cross-overs at meiosis) remains difficult. This article derives a simple expression for the strength of selection acting on a modifier gene affecting the genetic map length of a whole chromosome or genome undergoing recurrent mutation. In particular, it shows that indirect selection for recombination caused by interference among mutations is proportional to [Formula: see text], where [Formula: see text] is the effective population size, U is the deleterious mutation rate per chromosome, and R is the chromosome map length. Indirect selection is relatively insensitive to the fitness effects of deleterious alleles, epistasis, or the genetic architecture of recombination rate variation and may compensate for substantial costs associated with recombination when linkage is tight. However, its effect generally stays weak in large, highly recombining populations.

Reduced meiotic recombination in rhesus macaques and the origin of the human recombination landscape

Characterizing meiotic recombination rates across the genomes of nonhuman primates is important for understanding the genetics of primate populations, performing genetic analyses of phenotypic variation and reconstructing the evolution of human recombination. Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primates in biomedical research. We constructed a high-resolution genetic map of the rhesus genome based on whole genome sequence data from Indian-origin rhesus macaques. The genetic markers used were approximately 18 million SNPs, with marker density 6.93 per kb across the autosomes. We report that the genome-wide recombination rate in rhesus macaques is significantly lower than rates observed in apes or humans, while the distribution of recombination across the macaque genome is more uniform. These observations provide new comparative information regarding the evolution of recombination in primates.

Population-Specific Recombination Maps from Segments of Identity by Descent

Recombination rates vary significantly across the genome, and estimates of recombination rates are needed for downstream analyses such as haplotype phasing and genotype imputation. Existing methods for recombination rate estimation are limited by insufficient amounts of informative genetic data or by high computational cost. We present a method and software, called IBDrecomb, for using segments of identity by descent to infer recombination rates. IBDrecomb can be applied to sequenced population cohorts to obtain high-resolution, population-specific recombination maps. In simulated admixed data, IBDrecomb obtains higher accuracy than admixture-based estimation of recombination rates. When applied to 2,500 simulated individuals, IBDrecomb obtains similar accuracy to a linkage-disequilibrium (LD)-based method applied to 96 individuals (the largest number for which computation is tractable). Compared to LD-based maps, our IBD-based maps have the advantage of estimating recombination rates in the recent past rather than the distant past. We used IBDrecomb to generate new recombination maps for European Americans and for African Americans from TOPMed sequence data from the Framingham Heart Study (1,626 unrelated individuals) and the Jackson Heart Study (2,046 unrelated individuals), and we compare them to LD-based, admixture-based, and family-based maps.

Mechanistic Insight into Crossing over during Mouse Meiosis

Crossover recombination is critical for meiotic chromosome segregation, but how mammalian crossing over is accomplished is poorly understood. Here, we illuminate how strands exchange during meiotic recombination in male mice by analyzing patterns of heteroduplex DNA in recombinant molecules preserved by the mismatch correction deficiency of Msh2^-/- mutants. Surprisingly, MSH2-dependent recombination suppression was not evident. However, a substantial fraction of crossover products retained heteroduplex DNA, and some provided evidence of MSH2-independent correction. Biased crossover resolution was observed, consistent with asymmetry between DNA ends in earlier intermediates. Many crossover products yielded no heteroduplex DNA, suggesting dismantling by D-loop migration. Unlike the complexity of crossovers in yeast, these simple modifications of the original double-strand break repair model-asymmetry in recombination intermediates and D-loop migration-may be sufficient to explain most meiotic crossing over in mice while also addressing long-standing questions related to Holliday junction resolution.

Natural Selection Shapes Variation in Genome-wide Recombination Rate in Drosophila pseudoobscura

While recombination is widely recognized to be a key modulator of numerous evolutionary phenomena, we have a poor understanding of how recombination rate itself varies and evolves within a species. Here, we performed a comprehensive study of recombination rate (rate of meiotic crossing over) in two natural populations of Drosophila pseudoobscura from Utah and Arizona, USA. We used an amplicon sequencing approach to obtain high-quality genotypes in approximately 8,000 individual backcrossed offspring (17 mapping populations with roughly 530 individuals each), for which we then quantified crossovers. Interestingly, variation in recombination rate within and between populations largely manifested as differences in genome-wide recombination rate rather than remodeling of the local recombination landscape. Comparing populations, we discovered individuals from the Utah population displayed on average 8% higher crossover rates than the Arizona population, a statistically significant difference. Using a Q_ST-F_ST analysis, we found that this difference in crossover rate was dramatically higher than expected under neutrality, indicating that this difference may have been driven by natural selection. Finally, using a combination of short- and long-read whole-genome sequencing, we found no significant association between crossover rate and structural variation at the 200-400 kb scale. Our results demonstrate that (1) there is abundant variation in genome-wide crossover rate in natural populations, (2) at the 200-400 kb scale, recombination rate appears to vary largely genome-wide, rather than in specific intervals, and (3) interpopulation differences in recombination rate may be the result of local adaptation.

Weak Correlation between Nucleotide Variation and Recombination Rate across the House Mouse Genome

Positive selection and purifying selection reduce levels of variation at linked neutral loci. One consequence of these processes is that the amount of neutral diversity and the meiotic recombination rate are predicted to be positively correlated across the genome-a prediction met in some species but not others. To better document the prevalence of selection at linked sites, we used new and published whole-genome sequences to survey nucleotide variation in population samples of the western European house mouse (Mus musculus domesticus) from Germany, France, and Gough Island, a remote volcanic island in the south Atlantic. Correlations between sequence variation and recombination rates estimated independently from dense linkage maps were consistently very weak (ρ ≤ 0.06), though they exceeded conventional significance thresholds. This pattern persisted in comparisons between genomic regions with the highest and lowest recombination rates, as well as in models incorporating the density of transcribed sites, the density of CpG dinucleotides, and divergence between mouse and rat as covariates. We conclude that natural selection affects linked neutral variation in a restricted manner in the western European house mouse.

Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.

Mapping Recombination Rate on the Autosomal Chromosomes Based on the Persistency of Linkage Disequilibrium Phase Among Autochthonous Beef Cattle Populations in Spain

In organisms with sexual reproduction, genetic diversity, and genome evolution are governed by meiotic recombination caused by crossing-over, which is known to vary within the genome. In this study, we propose a simple method to estimate the recombination rate that makes use of the persistency of linkage disequilibrium (LD) phase among closely related populations. The biological material comprised 171 triplets (sire/dam/offspring) from seven populations of autochthonous beef cattle in Spain (Asturiana de los Valles, Avileña-Negra Ibérica, Bruna dels Pirineus, Morucha, Pirenaica, Retinta, and Rubia Gallega), which were genotyped for 777,962 SNPs with the BovineHD BeadChip. After standard quality filtering, we reconstructed the haplotype phases in the parental individuals and calculated the LD by the correlation -r- between each pair of markers that had a genetic distance < 1 Mb. Subsequently, these correlations were used to calculate the persistency of LD phase between each pair of populations along the autosomal genome. Therefore, the distribution of the recombination rate along the genome can be inferred since the effect of the number of generations of divergence should be equivalent throughout the genome. In our study, the recombination rate was highest in the largest chromosomes and at the distal portion of the chromosomes. In addition, the persistency of LD phase was highly heterogeneous throughout the genome, with a ratio of 25.4 times between the estimates of the recombination rates from the genomic regions that had the highest (BTA18-7.1 Mb) and the lowest (BTA12-42.4 Mb) estimates. Finally, an overrepresentation enrichment analysis (ORA) showed differences in the enriched gene ontology (GO) terms between the genes located in the genomic regions with estimates of the recombination rate over (or below) the 95^th (or 5^th) percentile throughout the autosomal genome.

Molecular evolution of the meiotic recombination pathway in mammals

Meiotic recombination shapes evolution and helps to ensure proper chromosome segregation in most species that reproduce sexually. Recombination itself evolves, with species showing considerable divergence in the rate of crossing-over. However, the genetic basis of this divergence is poorly understood. Recombination events are produced via a complicated, but increasingly well-described, cellular pathway. We apply a phylogenetic comparative approach to a carefully selected panel of genes involved in the processes leading to crossovers-spanning double-strand break formation, strand invasion, the crossover/non-crossover decision, and resolution-to reconstruct the evolution of the recombination pathway in eutherian mammals and identify components of the pathway likely to contribute to divergence between species. Eleven recombination genes, predominantly involved in the stabilization of homologous pairing and the crossover/non-crossover decision, show evidence of rapid evolution and positive selection across mammals. We highlight TEX11 and associated genes involved in the synaptonemal complex and the early stages of the crossover/non-crossover decision as candidates for the evolution of recombination rate. Evolutionary comparisons to MLH1 count, a surrogate for the number of crossovers, reveal a positive correlation between genome-wide recombination rate and the rate of evolution at TEX11 across the mammalian phylogeny. Our results illustrate the power of viewing the evolution of recombination from a pathway perspective.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

A first genetic portrait of synaptonemal complex variation

The synaptonemal complex (SC) is a proteinaceous scaffold required for synapsis and recombination between homologous chromosomes during meiosis. Although the SC has been linked to differences in genome-wide crossover rates, the genetic basis of standing variation in SC structure remains unknown. To investigate the possibility that recombination evolves through changes to the SC, we characterized the genetic architecture of SC divergence on two evolutionary timescales. Applying a novel digital image analysis technique to spermatocyte spreads, we measured total SC length in 9,532 spermatocytes from recombinant offspring of wild-derived mouse strains with differences in this fundamental meiotic trait. Using this large dataset, we identified the first known genomic regions involved in the evolution of SC length. Distinct loci affect total SC length divergence between and within subspecies, with the X chromosome contributing to both. Joint genetic analysis of MLH1 foci—immunofluorescent markers of crossovers—from the same spermatocytes revealed that two of the identified loci also confer differences in the genome-wide recombination rate. Causal mediation analysis suggested that one pleiotropic locus acts early in meiosis to designate crossovers prior to SC assembly, whereas a second locus primarily shapes crossover number through its effect on SC length. One genomic interval shapes the relationship between SC length and recombination rate, likely modulating the strength of crossover interference. Our findings pinpoint SC formation as a key step in the evolution of recombination and demonstrate the power of genetic mapping on standing variation in the context of the recombination pathway.

During the first stages of meiosis, the chromosome axes are organized along a protein scaffold in preparation for recombination and their subsequent segregation. This scaffold, known as the synaptonemal complex (SC), is critical for the regular progression of recombination. A complex relationship exists between the organization of the SC, the frequency of recombination, and the likelihood of improper chromosome segregation. In this study, we investigate the genetics of synaptonemal complex variation in the house mouse and connect it with variation in the rate of recombination. We found five loci and several compelling candidate genes responsible for the evolution of synaptonemal complex length within and between mouse subspecies. Several of these loci also affect recombination rate, and our joint analyses of the phenotypes suggest an order by which their effects manifest within the recombination pathway. Our results show that evolution of SC length is crucial to recombination rate divergence. Our work here also demonstrates that genetic analysis of additional meiotic phenotypes can help explain the evolution of recombination, a fundamental evolutionary force.

Conservation of the genome-wide recombination rate in white-footed mice

Despite being linked to the fundamental processes of chromosome segregation and offspring diversification, meiotic recombination rates vary within and between species. Recent years have seen progress in quantifying recombination rate evolution across multiple temporal and genomic scales. Nevertheless, the level of variation in recombination rate within wild populations-a key determinant of evolution in this trait-remains poorly documented on the genomic scale. To address this notable gap, we used immunofluorescent cytology to quantify genome-wide recombination rates in males from a wild population of the white-footed mouse, Peromyscus leucopus. For comparison, we measured recombination rates in a second population of male P. leucopus raised in the laboratory and in male deer mice from the subspecies Peromyscus maniculatus bairdii. Although we found differences between individuals in the genome-wide recombination rate, levels of variation were low-within populations, between populations, and between species. Quantification of synaptonemal complex length and crossover positions along chromosome 1 using a novel automated approach also revealed conservation in broad-scale crossover patterning, including strong crossover interference. We propose stabilizing selection targeting recombination or correlated processes as the explanation for these patterns.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

A rigorous measure of genome-wide genetic shuffling that takes into account crossover positions and Mendel's second law

Comparative studies in evolutionary genetics rely critically on evaluation of the total amount of genetic shuffling that occurs during gamete production. Such studies have been hampered by the absence of a direct measure of this quantity. Existing measures consider crossing-over by simply counting the average number of crossovers per meiosis. This is qualitatively inadequate, because the positions of crossovers along a chromosome are also critical: a crossover toward the middle of a chromosome causes more shuffling than a crossover toward the tip. Moreover, traditional measures fail to consider shuffling from independent assortment of homologous chromosomes (Mendel's second law). Here, we present a rigorous measure of genome-wide shuffling that does not suffer from these limitations. We define the parameter [Formula: see text] as the probability that the alleles at two randomly chosen loci are shuffled during gamete production. This measure can be decomposed into separate contributions from crossover number and position and from independent assortment. Intrinsic implications of this metric include the fact that [Formula: see text] is larger when crossovers are more evenly spaced, which suggests a selective advantage of crossover interference. Utilization of [Formula: see text] is enabled by powerful emergent methods for determining crossover positions either cytologically or by DNA sequencing. Application of our analysis to such data from human male and female reveals that (i) [Formula: see text] in humans is close to its maximum possible value of 1/2 and that (ii) this high level of shuffling is due almost entirely to independent assortment, the contribution of which is ∼30 times greater than that of crossovers.

Variation in recombination frequency and distribution across eukaryotes: patterns and processes

Recombination, the exchange of DNA between maternal and paternal chromosomes during meiosis, is an essential feature of sexual reproduction in nearly all multicellular organisms. While the role of recombination in the evolution of sex has received theoretical and empirical attention, less is known about how recombination rate itself evolves and what influence this has on evolutionary processes within sexually reproducing organisms. Here, we explore the patterns of, and processes governing recombination in eukaryotes. We summarize patterns of variation, integrating current knowledge with an analysis of linkage map data in 353 organisms. We then discuss proximate and ultimate processes governing recombination rate variation and consider how these influence evolutionary processes. Genome-wide recombination rates (cM/Mb) can vary more than tenfold across eukaryotes, and there is large variation in the distribution of recombination events across closely related taxa, populations and individuals. We discuss how variation in rate and distribution relates to genome architecture, genetic and epigenetic mechanisms, sex, environmental perturbations and variable selective pressures. There has been great progress in determining the molecular mechanisms governing recombination, and with the continued development of new modelling and empirical approaches, there is now also great opportunity to further our understanding of how and why recombination rate varies.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Recombination: the good, the bad and the variable

Recombination, the process by which DNA strands are broken and repaired, producing new combinations of alleles, occurs in nearly all multicellular organisms and has important implications for many evolutionary processes. The effects of recombination can be good, as it can facilitate adaptation, but also bad when it breaks apart beneficial combinations of alleles, and recombination is highly variable between taxa, species, individuals and across the genome. Understanding how and why recombination rate varies is a major challenge in biology. Most theoretical and empirical work has been devoted to understanding the role of recombination in the evolution of sex—comparing between sexual and asexual species or populations. How recombination rate evolves and what impact this has on evolutionary processes within sexually reproducing organisms has received much less attention. This Theme Issue focusses on how and why recombination rate varies in sexual species, and aims to coalesce knowledge of the molecular mechanisms governing recombination with our understanding of the evolutionary processes driving variation in recombination within and between species. By integrating these fields, we can identify important knowledge gaps and areas for future research, and pave the way for a more comprehensive understanding of how and why recombination rate varies.

Molecular Evolution at a Meiosis Gene Mediates Species Differences in the Rate and Patterning of Recombination

Crossing over between homologous chromosomes during meiosis repairs programmed DNA double-strand breaks, ensures proper segregation at meiosis I [1], shapes the genomic distribution of nucleotide variability in populations, and enhances the efficacy of natural selection among genetically linked sites [2]. Between closely related Drosophila species, large differences exist in the rate and chromosomal distribution of crossing over. Little, however, is known about the molecular genetic changes or population genetic forces that mediate evolved differences in recombination between species [3, 4]. Here, we show that a meiosis gene with a history of rapid evolution acts as a trans-acting modifier of species differences in crossing over. In transgenic flies, the dicistronic gene, mei-217/mei-218, recapitulates a large part of the species differences in the rate and chromosomal distribution of crossing over. These phenotypic differences appear to result from changes in protein sequence not gene expression. Our population genetics analyses show that the protein-coding sequence of mei-218, but not mei-217, has a history of recurrent positive natural selection. By modulating the intensity of centromeric and telomeric suppression of crossing over, evolution at mei-217/-218 has incidentally shaped gross differences in the chromosomal distribution of nucleotide variability between species. We speculate that recurrent bouts of adaptive evolution at mei-217/-218 might reflect a history of coevolution with selfish genetic elements.

Low recombination rates in sexual species and sex-asex transitions

In most sexual, diploid eukaryotes, at least one crossover occurs between each pair of homologous chromosomes during meiosis, presumably in order to ensure proper segregation. Well-known exceptions to this rule are species in which one sex does not recombine and specific chromosomes lacking crossover. We review other possible exceptions, including species with chromosome maps of less than 50 cM in one or both sexes. We discuss the idea that low recombination rates may favour sex-asex transitions, or, alternatively may be a consequence of it. We then show that a yet undescribed species of brine shrimp Artemia from Kazakhstan (A sp. Kazakhstan), the closest known relative of the asexual Artemia parthenogenetica, has one of the shortest genetic linkage maps known. Based on a family of 42 individuals and 589 RAD markers, we find that many linkage groups are considerably shorter than 50 cM, suggesting either no obligate crossover or crossovers concentrated at terminal positions with little effect on recombination. We contrast these findings with the published map of the more distantly related sexual congener, A. franciscana, and conclude that the study of recombination in non-model systems is important to understand the evolutionary causes and consequences of recombination.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.