Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

Variation in fine-scale recombination rate in temperature-evolved Drosophila melanogaster populations in response to selection

Meiotic recombination plays a critical evolutionary role in maintaining fitness in response to selective pressures due to changing environments. Variation in recombination rate has been observed amongst and between species and populations and within genomes across numerous taxa. Studies have demonstrated a link between changes in recombination rate and selection, but the extent to which fine-scale recombination rate varies between evolved populations during the evolutionary period in response to selection is under active research. Here, we utilize a set of 3 temperature-evolved Drosophila melanogaster populations that were shown to have diverged in several phenotypes, including recombination rate, based on the temperature regime in which they evolved. Using whole-genome sequencing data from these populations, we generated linkage disequilibrium-based fine-scale recombination maps for each population. With these maps, we compare recombination rates and patterns among the 3 populations and show that they have diverged at fine scales but are conserved at broader scales. We further demonstrate a correlation between recombination rates and genomic variation in the 3 populations. Lastly, we show variation in localized regions of enhanced recombination rates, termed warm spots, between the populations with these warm spots and associated genes overlapping areas previously shown to have diverged in the 3 populations due to selection. These data support the existence of recombination modifiers in these populations which are subject to selection during evolutionary change.

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.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

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.

Variation in Recombination Rate Is Shaped by Domestication and Environmental Conditions in Barley

Meiotic recombination generates genetic diversity upon which selection can act. Recombination rates are highly variable between species, populations, individuals, sexes, chromosomes, and chromosomal regions. The underlying mechanisms are controlled at the genetic and epigenetic level and show plasticity toward the environment. Environmental plasticity may be divided into short- and long-term responses. We estimated recombination rates in natural populations of wild barley and domesticated landraces using a population genetics approach. We analyzed recombination landscapes in wild barley and domesticated landraces at high resolution. In wild barley, high recombination rates are found in more interstitial chromosome regions in contrast to distal chromosome regions in domesticated barley. Among subpopulations of wild barley, natural variation in effective recombination rate is correlated with temperature, isothermality, and solar radiation in a nonlinear manner. A positive linear correlation was found between effective recombination rate and annual precipitation. We discuss our findings with respect to how the environment might shape effective recombination rates in natural populations. Higher recombination rates in wild barley populations subjected to specific environmental conditions could be a means to maintain fitness in a strictly inbreeding species.

Crossover Interference: Shedding Light on the Evolution of Recombination

Through recombination, genes are freed to evolve more independently of one another, unleashing genetic variance hidden in the linkage disequilibrium that accumulates through selection combined with drift. Yet crossover numbers are evolutionarily constrained, with at least one and not many more than one crossover per bivalent in most taxa. Crossover interference, whereby a crossover reduces the probability of a neighboring crossover, contributes to this homogeneity. The mechanisms by which interference is achieved and crossovers are regulated are a major current subject of inquiry, facilitated by novel methods to visualize crossovers and to pinpoint recombination events. Here, we review patterns of crossover interference and the models built to describe this process. We then discuss the selective forces that have likely shaped interference and the regulation of crossover numbers.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

The Genomic Landscape of Crossover Interference in the Desert Tree Populus euphratica

Crossover (CO) interference is a universal phenomenon by which the occurrence of one CO event inhibits the simultaneous occurrence of other COs along a chromosome. Because of its critical role in the evolution of genome structure and organization, the cytological and molecular mechanisms underlying CO interference have been extensively investigated. However, the genome-wide distribution of CO interference and its interplay with sex-, stress-, and age-induced differentiation remain poorly understood. Multi-point linkage analysis has proven to be a powerful tool for landscaping CO interference, especially within species for which CO mutants are rarely available. We implemented four-point linkage analysis to landscape a detailed picture of how CO interference is distributed through the entire genome of Populus euphratica, the only forest tree that can survive and grow in saline desert. We identified an extensive occurrence of CO interference, and found that its strength depends on the length of chromosomes and the genomic locations within the chromosome. We detected high-order CO interference, possibly suggesting a highly complex mechanism crucial for P. euphratica to grow, reproduce, and evolve in its harsh environment.

Role of Cis, Trans, and Inbreeding Effects on Meiotic Recombination in Saccharomyces cerevisiae

Meiotic recombination is a major driver of genome evolution by creating new genetic combinations. To probe the factors driving variability of meiotic recombination, we used a high-throughput method to measure recombination rates in hybrids between SK1 and a total of 26 Saccharomyces cerevisiae strains from different geographic origins and habitats. Fourteen intervals were monitored for each strain, covering chromosomes VI and XI entirely, and part of chromosome I. We found an average number of crossovers per chromosome ranging between 1.0 and 9.5 across strains ("domesticated" or not), which is higher than the average between 0.5 and 1.5 found in most organisms. In the different intervals analyzed, recombination showed up to ninefold variation across strains but global recombination landscapes along chromosomes varied less. We also built an incomplete diallel experiment to measure recombination rates in one region of chromosome XI in 10 different crosses involving five parental strains. Our overall results indicate that recombination rate is increasingly positively correlated with sequence similarity between homologs (i) in DNA double-strand-break-rich regions within intervals, (ii) in entire intervals, and (iii) at the whole genome scale. Therefore, these correlations cannot be explained by cis effects only. We also estimated that cis and trans effects explained 38 and 17%, respectively, of the variance of recombination rate. In addition, by using a quantitative genetics analysis, we identified an inbreeding effect that reduces recombination rate in homozygous genotypes, while other interaction effects (specific combining ability) or additive effects (general combining ability) are found to be weak. Finally, we measured significant crossover interference in some strains, and interference intensity was positively correlated with crossover number.

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'.

Evolutionary mysteries in meiosis

Meiosis is a key event of sexual life cycles in eukaryotes. Its mechanistic details have been uncovered in several model organisms, and most of its essential features have received various and often contradictory evolutionary interpretations. In this perspective, we present an overview of these often 'weird' features. We discuss the origin of meiosis (origin of ploidy reduction and recombination, two-step meiosis), its secondary modifications (in polyploids or asexuals, inverted meiosis), its importance in punctuating life cycles (meiotic arrests, epigenetic resetting, meiotic asymmetry, meiotic fairness) and features associated with recombination (disjunction constraints, heterochiasmy, crossover interference and hotspots). We present the various evolutionary scenarios and selective pressures that have been proposed to account for these features, and we highlight that their evolutionary significance often remains largely mysterious. Resolving these mysteries will likely provide decisive steps towards understanding why sex and recombination are found in the majority of eukaryotes.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Landscaping Crossover Interference Across a Genome

The evolutionary success of eukaryotic organisms crucially depends on the capacity to produce genetic diversity through reciprocal exchanges of each chromosome pair, or crossovers (COs), during meiosis. It has been recognized that COs arise more evenly across a given chromosome than at random. This phenomenon, termed CO interference, occurs pervasively in eukaryotes and may confer a selective advantage. We describe here a multipoint linkage analysis procedure for segregating families to quantify the strength of CO interference over the genome, and extend this procedure to illustrate the landscape of CO interference in natural populations. We further discuss the crucial role of CO interference in amplifying and maintaining genetic diversity through sex-, stress-, and age-induced differentiation.

Building Ultra-High-Density Linkage Maps Based on Efficient Filtering of Trustable Markers

The study is focused on addressing the problem of building genetic maps in the presence of ∼10³-10⁴ of markers per chromosome. We consider a spectrum of situations with intrachromosomal heterogeneity of recombination rate, different level of genotyping errors, and missing data. In the ideal scenario of the absence of errors and missing data, the majority of markers should appear as groups of cosegregating markers ("twins") representing no challenge for map construction. The central aspect of the proposed approach is to take into account the structure of the marker space, where each twin group (TG) and singleton markers are represented as points of this space. The confounding effect of genotyping errors and missing data leads to reduction of TG size, but upon a low level of these effects surviving TGs can still be used as a source of reliable skeletal markers. Increase in the level of confounding effects results in a considerable decrease in the number or even disappearance of usable TGs and, correspondingly, of skeletal markers. Here, we show that the paucity of informative markers can be compensated by detecting kernels of markers in the marker space using a clustering procedure, and demonstrate the utility of this approach for high-density genetic map construction on simulated and experimentally obtained genotyping datasets.

Hill-Robertson interference maintained by Red Queen dynamics favours the evolution of sex

Although it is well established theoretically that selective interference among mutations (Hill-Robertson interference) favours meiotic recombination, genomewide mean rates of mutation and strengths of selection appear too low to support this as the mechanism favouring recombination in nature. A possible solution to this discrepancy between theory and observation is that selection is at least intermittently very strong due to the antagonistic coevolution between a host and its parasites. The Red Queen theory posits that such coevolution generates fitness epistasis among loci, which generates negative linkage disequilibrium among beneficial mutations, which in turn favours recombination. This theory has received only limited support. However, Red Queen dynamics without epistasis may provide the ecological conditions that maintain strong and frequent selective interference in finite populations that indirectly selects for recombination. This hypothesis is developed here through the simulation of Red Queen dynamics. This approach required the development of a method to calculate the exact frequencies of multilocus haplotypes after recombination. Simulations show that recombination is favoured by the moderately weak selection of many loci involved in the interaction between a host and its parasites, which results in substitution rates that are compatible with empirical estimates. The model also reproduces the previously reported rapid increase in the rate of outcrossing in Caenorhabditis elegans coevolving with a bacterial pathogen.

Selection Signatures in Four Lignin Genes from Switchgrass Populations Divergently Selected for In Vitro Dry Matter Digestibility

Switchgrass is undergoing development as a dedicated cellulosic bioenergy crop. Fermentation of lignocellulosic biomass to ethanol in a bioenergy system or to volatile fatty acids in a livestock production system is strongly and negatively influenced by lignification of cell walls. This study detects specific loci that exhibit selection signatures across switchgrass breeding populations that differ in in vitro dry matter digestibility (IVDMD), ethanol yield, and lignin concentration. Allele frequency changes in candidate genes were used to detect loci under selection. Out of the 183 polymorphisms identified in the four candidate genes, twenty-five loci in the intron regions and four loci in coding regions were found to display a selection signature. All loci in the coding regions are synonymous substitutions. Selection in both directions were observed on polymorphisms that appeared to be under selection. Genetic diversity and linkage disequilibrium within the candidate genes were low. The recurrent divergent selection caused excessive moderate allele frequencies in the cycle 3 reduced lignin population as compared to the base population. This study provides valuable insight on genetic changes occurring in short-term selection in the polyploid populations, and discovered potential markers for breeding switchgrass with improved biomass quality.

Effect of sex, age and genetics on crossover interference in cattle

Crossovers generated by homologous recombination ensure proper chromosome segregation during meiosis. Crossover interference results in chiasmata being more evenly distributed along chromosomes, but the mechanism underlying crossover interference remains elusive. Based on large pedigrees of Holstein and Jersey cattle with genotype data, we extracted three-generation families, including 147,327 male and 71,687 female meioses in Holstein, and 108,163 male and 37,008 female meioses in Jersey, respectively. We identified crossovers in these meioses and fitted the Housworth-Stahl "interference-escape" model to study crossover interference patterns in the cattle genome. Our result reveals that the degree of crossover interference is stronger in females than in males. We found evidence for inter-chromosomal variation in the level of crossover interference, with smaller chromosomes exhibiting stronger interference. In addition, crossover interference levels decreased with maternal age. Finally, sex-specific GWAS analyses identified one locus near the NEK9 gene on chromosome 10 to have a significant effect on crossover interference levels. This locus has been previously associated with recombination rate in cattle. Collectively, this large-scale analysis provided a comprehensive description of crossover interference across chromosome, sex and age groups, identified associated candidate genes, and produced useful insights into the mechanism of crossover interference.

Rapid genomic changes in Drosophila melanogaster adapting to desiccation stress in an experimental evolution system

BACKGROUND

Experimental evolution studies, coupled with whole genome resequencing and advances in bioinformatics, have become a powerful tool for exploring how populations respond to selection at the genome-wide level, complementary to genome-wide association studies (GWASs) and linkage mapping experiments as strategies to connect genotype and phenotype. In this experiment, we analyzed genomes of Drosophila melanogaster from lines evolving under long-term directional selection for increased desiccation resistance in comparison with control (no-selection) lines.

RESULTS

We demonstrate that adaptive responses to desiccation stress have exerted extensive footprints on the genomes, manifested through a high degree of fixation of alleles in surrounding neighborhoods of eroded heterozygosity. These patterns were highly convergent across replicates, consistent with signatures of 'soft' selective sweeps, where multiple alleles present as standing genetic variation become beneficial and sweep through the replicate populations at the same time. Albeit much less frequent, we also observed line-unique sweep regions with zero or near-zero heterozygosity, consistent with classic, or 'hard', sweeps, where novel rather than pre-existing adaptive mutations may have been driven to fixation. Genes responsible for cuticle and protein deubiquitination seemed to be central to these selective sweeps. High divergence within coding sequences between selected and control lines was also reflected by significant results of the McDonald-Kreitman and Ka/Ks tests, showing that as many as 347 genes may have been under positive selection.

CONCLUSIONS

Desiccation stress, a common challenge to many organisms inhabiting dry environments, proves to be a very potent selecting factor having a big impact on genome diversity.