The Genetic Architecture of Natural Variation in Recombination Rate in Drosophila melanogaster

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Genome-wide association in Drosophila identifies a role for Piezo and Proc-R in sleep latency

Sleep latency, the amount of time that it takes an individual to fall asleep, is a key indicator of sleep need. Sleep latency varies considerably both among and within species and is heritable, but lacks a comprehensive description of its underlying genetic network. Here we conduct a genome-wide association study of sleep latency. Using previously collected sleep and activity data on a wild-derived population of flies, we calculate sleep latency, confirming significant, heritable genetic variation for this complex trait. We identify 520 polymorphisms in 248 genes contributing to variability in sleep latency. Tests of mutations in 23 candidate genes and additional putative pan-neuronal knockdown of 9 of them implicated CG44153, Piezo, Proc-R and Rbp6 in sleep latency. Two large-effect mutations in the genes Proc-R and Piezo were further confirmed via genetic rescue. This work greatly enhances our understanding of the genetic factors that influence variation in sleep latency.

Genetic architecture of individual meiotic crossover rate and distribution in Atlantic Salmon

Meiotic recombination through chromosomal crossovers ensures proper segregation of homologous chromosomes during meiosis, while also breaking down linkage disequilibrium and shuffling alleles at loci located on the same chromosome. Rates of recombination can vary between species, but also between and within individuals, sex and chromosomes within species. Indeed, the Atlantic salmon genome is known to have clear sex differences in recombination with female biased heterochiasmy and markedly different landscapes of crossovers between males and females. In male meiosis, crossovers occur strictly in the telomeric regions, whereas in female meiosis crossovers tend to occur closer to the centromeres. However, little is known about the genetic control of these patterns and how this differs at the individual level. Here, we investigate genetic variation in individual measures of recombination in > 5000 large full-sib families of a Norwegian Atlantic salmon breeding population with high-density SNP genotypes. We show that females had 1.6 × higher crossover counts (CC) than males, with autosomal linkage maps spanning a total of 2174 cM in females and 1483 cM in males. However, because of the extreme telomeric bias of male crossovers, female recombination is much more important for generation of new haplotypes with 8 × higher intra-chromosomal genetic shuffling than males. CC was heritable in females (h2 = 0.11) and males (h2 = 0.10), and shuffling was also heritable in both sex but with a lower heritability in females (h2 = 0.06) than in males (h2 = 0.11). Inter-sex genetic correlations for both traits were close to zero, suggesting that rates and distribution of crossovers are genetically distinct traits in males and females, and that there is a potential for independent genetic change in both sexes in the Atlantic Salmon. Together, these findings give novel insights into the genetic architecture of recombination in salmonids and contribute to a better understanding of how rates and distribution of recombination may evolve in eukaryotes more broadly.

Measuring the Impact of Genetic Heterogeneity and Chromosomal Inversions on the Efficacy of CRISPR-Cas9 Gene Drives in Different Strains of Anopheles gambiae

The human malaria vector Anopheles gambiae is becoming increasingly resistant to insecticides, spurring the development of genetic control strategies. CRISPR-Cas9 gene drives can modify a population by creating double-stranded breaks at highly specific targets, triggering copying of the gene drive into the cut site ("homing"), ensuring its inheritance. The DNA repair mechanism responsible requires homology between the donor and recipient chromosomes, presenting challenges for the invasion of laboratory-developed gene drives into wild populations of target species An. gambiae species complex, which show high levels of genome variation. Two gene drives (vas2-5958 and zpg-7280) were introduced into three An. gambiae strains collected across Africa with 5.3-6.6% variation around the target sites, and the effect of this variation on homing was measured. Gene drive homing across different karyotypes of the 2La chromosomal inversion was also assessed. No decrease in gene drive homing was seen despite target site heterology, demonstrating the applicability of gene drives to wild populations.

Modeling the evolution of recombination plasticity: A prospective review

Meiotic recombination is one of the main sources of genetic variation, a fundamental factor in the evolutionary adaptation of sexual eukaryotes. Yet, the role of variation in recombination rate and other recombination features remains underexplored. In this review, we focus on the sensitivity of recombination rates to different extrinsic and intrinsic factors. We briefly present the empirical evidence for recombination plasticity in response to environmental perturbations and/or poor genetic background and discuss theoretical models developed to explain how such plasticity could have evolved and how it can affect important population characteristics. We highlight a gap between the evidence, which comes mostly from experiments with diploids, and theory, which typically assumes haploid selection. Finally, we formulate open questions whose solving would help to outline conditions favoring recombination plasticity. This will contribute to answering the long-standing question of why sexual recombination exists despite its costs, since plastic recombination may be evolutionary advantageous even in selection regimes rejecting any non-zero constant recombination.

Ovarian transcriptional response to Wolbachia infection in D. melanogaster in the context of between-genotype variation in gene expression

Wolbachia is a maternally transmitted endosymbiotic bacteria that infects a wide variety of arthropod and nematode hosts. The effects of Wolbachia on host biology are far-reaching and include changes in host gene expression. However, previous work on the host transcriptional response has generally been investigated in the context of a single host genotype. Thus, the relative effect of Wolbachia infection versus vs. host genotype on gene expression is unknown. Here, we explicitly test the relative roles of Wolbachia infection and host genotype on host gene expression by comparing the ovarian transcriptomes of 4 strains of Drosophila melanogaster (D. melanogaster) infected and uninfected with Wolbachia. Our data suggest that infection explains a small amount of transcriptional variation, particularly in comparison to variation in gene expression among strains. However, infection specifically affects genes related to cell cycle, translation, and metabolism. We also find enrichment of cell division and recombination processes among genes with infection-associated differential expression. Broadly, the transcriptomic changes identified in this study provide novel understanding of the relative magnitude of the effect of Wolbachia infection on gene expression in the context of host genetic variation and also point to genes that are consistently differentially expressed in response to infection among multiple genotypes.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

Variation and genetic control of individual recombination rates in Norwegian Red dairy cattle

Meiotic recombination is an important evolutionary mechanism that breaks up linkages between loci and creates novel haplotypes for selection to act upon. Understanding the genetic control of variation in recombination rates is therefore of great interest in both natural and domestic breeding populations. In this study, we used pedigree information and medium-density (∼50K) genotyped data in a large cattle (Bos taurus) breeding population in Norway (Norwegian Red cattle) to investigate recombination rate variation between sexes and individual animals. Sex-specific linkage mapping showed higher rates in males than in females (total genetic length of autosomes = 2,492.9 cM in males and 2,308.9 cM in females). However, distribution of recombination along the genome showed little variation between males and females compared with that in other species. The heritability of autosomal crossover count was low but significant in both sexes (h² = 0.04 and 0.09 in males and females, respectively). We identified 2 loci associated with variation in individual crossover counts in female, one close to the candidate gene CEP55 and one close to both MLH3 and NEK9. All 3 genes have been associated with recombination rates in other cattle breeds. Our study contributes to the understanding of how recombination rates are controlled and how they may vary between closely related breeds as well as between species.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

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.

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.

Genomic prediction of the recombination rate variation in barley - A route to highly recombinogenic genotypes

Meiotic recombination is not only fundamental to the adaptation of sexually reproducing eukaryotes in nature but increased recombination rates facilitate the combination of favourable alleles into a single haplotype in breeding programmes. The main objectives of this study were to (i) assess the extent and distribution of the recombination rate variation in cultivated barley (Hordeum vulgare L.), (ii) quantify the importance of the general and specific recombination effects, and (iii) evaluate a genomic selection approach's ability to predict the recombination rate variation. Genetic maps were created for the 45 segregating populations that were derived from crosses among 23 spring barley inbreds with origins across the world. The genome-wide recombination rate among populations ranged from 0.31 to 0.73 cM/Mbp. The crossing design used in this study allowed to separate the general recombination effects (GRE) of individual parental inbreds from the specific recombination effects (SRE) caused by the combinations of parental inbreds. The variance of the genome-wide GRE was found to be about eight times the variance of the SRE. This finding indicated that parental inbreds differ in the efficiency of their recombination machinery. The ability to predict the chromosome or genome-wide recombination rate of an inbred ranged from 0.80 to 0.85. These results suggest that a reliable screening of large genetic materials for their potential to cause a high extent of genetic recombination in their progeny is possible, allowing to systematically manipulate the recombination rate using natural variation.

Genome-wide recombination map construction from single sperm sequencing in cattle

BACKGROUND

Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle.

RESULTS

To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle.

CONCLUSIONS

To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Meiotic, genomic and evolutionary properties of crossover distribution in Drosophila yakuba

The number and location of crossovers across genomes are highly regulated during meiosis, yet the key components controlling them are fast evolving, hindering our understanding of the mechanistic causes and evolutionary consequences of changes in crossover rates. Drosophila melanogaster has been a model species to study meiosis for more than a century, with an available high-resolution crossover map that is, nonetheless, missing for closely related species, thus preventing evolutionary context. Here, we applied a novel and highly efficient approach to generate whole-genome high-resolution crossover maps in D. yakuba to tackle multiple questions that benefit from being addressed collectively within an appropriate phylogenetic framework, in our case the D. melanogaster species subgroup. The genotyping of more than 1,600 individual meiotic events allowed us to identify several key distinct properties relative to D. melanogaster. We show that D. yakuba, in addition to higher crossover rates than D. melanogaster, has a stronger centromere effect and crossover assurance than any Drosophila species analyzed to date. We also report the presence of an active crossover-associated meiotic drive mechanism for the X chromosome that results in the preferential inclusion in oocytes of chromatids with crossovers. Our evolutionary and genomic analyses suggest that the genome-wide landscape of crossover rates in D. yakuba has been fairly stable and captures a significant signal of the ancestral crossover landscape for the whole D. melanogaster subgroup, even informative for the D. melanogaster lineage. Contemporary crossover rates in D. melanogaster, on the other hand, do not recapitulate ancestral crossovers landscapes. As a result, the temporal stability of crossover landscapes observed in D. yakuba makes this species an ideal system for applying population genetic models of selection and linkage, given that these models assume temporal constancy in linkage effects. Our studies emphasize the importance of generating multiple high-resolution crossover rate maps within a coherent phylogenetic context to broaden our understanding of crossover control during meiosis and to improve studies on the evolutionary consequences of variable crossover rates across genomes and time.

Soft selective sweeps: Addressing new definitions, evaluating competing models, and interpreting empirical outliers

The ability to accurately identify and quantify genetic signatures associated with soft selective sweeps based on patterns of nucleotide variation has remained controversial. We here provide counter viewpoints to recent publications in PLOS Genetics that have argued not only for the statistical identifiability of soft selective sweeps, but also for their pervasive evolutionary role in both Drosophila and HIV populations. We present evidence that these claims owe to a lack of consideration of competing evolutionary models, unjustified interpretations of empirical outliers, as well as to new definitions of the processes themselves. Our results highlight the dangers of fitting evolutionary models based on hypothesized and episodic processes without properly first considering common processes and, more generally, of the tendency in certain research areas to view pervasive positive selection as a foregone conclusion.

Diet-induced changes in titer support a discrete response of Wolbachia-associated plastic recombination in Drosophila melanogaster

Plastic recombination in Drosophila melanogaster has been associated with a variety of extrinsic and intrinsic factors such as temperature, starvation, and parasite infection. The bacterial endosymbiont Wolbachia pipientis has also been associated with plastic recombination in D. melanogaster. Wolbachia infection is pervasive in arthropods and this infection induces a variety of phenotypes in its hosts, the strength of which can depend on bacterial titer. Here, we test the hypothesis that the magnitude of Wolbachia-associated plastic recombination in D. melanogaster depends on titer. To manipulate titer, we raised Wolbachia-infected and uninfected flies on diets that have previously been shown to increase or decrease Wolbachia titer relative to controls. We measured recombination in treated and control individuals using a standard backcrossing scheme with two X-linked visible markers. Our results recapitulate previous findings that Wolbachia infection is associated with increased recombination rate across the yellow-vermillion interval of the X chromosome. Our data show no significant effect of diet or diet by Wolbachia interactions on recombination, suggesting that diet-induced changes in Wolbachia titer have no effect on the magnitude of plastic recombination. These findings represent one of the first steps toward investigating Wolbachia-associated plastic recombination and demonstrate that the phenotype is a discrete response rather than a continuous one.

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs' deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.

Recombination mapping of the Brazilian stingless bee Frieseomelitta varia confirms high recombination rates in social hymenoptera

Background

Meiotic recombination is a fundamental genetic process that shuffles allele combinations and promotes accurate segregation of chromosomes. Analyses of the ubiquitous variation of recombination rates within and across species suggest that recombination is evolving adaptively. All studied insects with advanced eusociality have shown exceptionally high recombination rates, which may represent a prominent case of adaptive evolution of recombination. However, our understanding of the relationship between social evolution and recombination rates is incomplete, partly due to lacking empirical data. Here, we present a linkage map of the monandrous, advanced eusocial Brazilian stingless bee, Frieseomelitta varia, providing the first recombination analysis in the diverse Meliponini (Hymenoptera, Apidae).

Results

Our linkage map includes 1417 markers in 19 linkage groups. This map spans approximately 2580 centimorgans, and comparisons to the physical genome assembly indicate that it covers more than 75 % of the 275 Megabasepairs (Mbp) F. varia genome. Thus, our study results in a genome-wide recombination rate estimate of 9.3–12.5 centimorgan per Mbp. This value is higher than estimates from nonsocial insects and comparable to other highly social species, although it does not support our prediction that monandry and strong queen-worker caste divergence of F. varia lead to even higher recombination rates than other advanced eusocial species.

Conclusions

Our study expands the association between elevated recombination and sociality in the order Hymenoptera and strengthens the support for the hypothesis that advanced social evolution in hymenopteran insects invariably selects for high genomic recombination rates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07987-3.

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

Fitness dependence preserves selection for recombination across diverse mixed mating strategies

Meiotic recombination and the factors affecting its rate and fate in nature have inspired many studies in theoretical evolutionary biology. Classical theoretical models have inferred that recombination can be favored under a rather restricted parameter range. Thus, the ubiquity of recombination in nature remains an open question. However, these models assumed constant recombination with an equal rate across all individuals within the population, whereas empirical evidence suggests that recombination may display certain sensitivity to ecological stressors and/or genotype fitness. Models assuming condition-dependent recombination show that such a strategy can often be favored over constant recombination. Moreover, in our recent model with panmictic populations subjected to purifying selection, fitness-dependent recombination was quite often favored even when any constant recombination was rejected. By using numerical modeling, we test whether such a 'recombination-rescuing potential' of fitness dependence holds also beyond panmixia, given the recognized effect of mating strategy on the evolution of recombination. We show that deviations from panmixia generally increase the recombination-rescuing potential of fitness dependence, with the strongest effect under intermediate selfing or high clonality. We find that under partial clonality, the evolutionary advantage of fitness-dependent recombination is determined mostly by selection against heterozygotes and additive-by-additive epistasis, while under partial selfing, additive-by-dominance epistasis is also a driver.

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.

Genomics of Recombination Rate Variation in Temperature-Evolved Drosophila melanogaster Populations

Meiotic recombination is a critical process that ensures proper segregation of chromosome homologs through DNA double-strand break repair mechanisms. Rates of recombination are highly variable among various taxa, within species, and within genomes with far-reaching evolutionary and genomic consequences. The genetic basis of recombination rate variation is therefore crucial in the study of evolutionary biology but remains poorly understood. In this study, we took advantage of a set of experimental temperature-evolved populations of Drosophila melanogaster with heritable differences in recombination rates depending on the temperature regime in which they evolved. We performed whole-genome sequencing and identified several chromosomal regions that appear to be divergent depending on temperature regime. In addition, we identify a set of single-nucleotide polymorphisms and associated genes with significant differences in allele frequency when the different temperature populations are compared. Further refinement of these gene candidates emphasizing those expressed in the ovary and associated with DNA binding reveals numerous potential candidate genes such as Hr38, EcR, and mamo responsible for observed differences in recombination rates in these experimental evolution lines thus providing insight into the genetic basis of recombination rate variation.

The Intracellular Symbiont Wolbachia pipientis Enhances Recombination in a Dose-Dependent Manner

Wolbachia pipientis is an intracellular alphaproteobacterium that infects 40%-60% of insect species and is well known for host reproductive manipulations. Although Wolbachia are primarily maternally transmitted, evidence of horizontal transmission can be found in incongruent host-symbiont phylogenies and recent acquisitions of the same Wolbachia strain by distantly related species. Parasitoids and predator-prey interactions may indeed facilitate the transfer of Wolbachia between insect lineages, but it is likely that Wolbachia are acquired via introgression in many cases. Many hypotheses exist to explain Wolbachia prevalence and penetrance, such as nutritional supplementation, protection from parasites, protection from viruses, or direct reproductive parasitism. Using classical genetics, we show that Wolbachia increase recombination in infected lineages across two genomic intervals. This increase in recombination is titer-dependent as the wMelPop variant, which infects at higher load in Drosophila melanogaster, increases recombination 5% more than the wMel variant. In addition, we also show that Spiroplasma poulsonii, another bacterial intracellular symbiont of D. melanogaster, does not induce an increase in recombination. Our results suggest that Wolbachia infection specifically alters its host's recombination landscape in a dose-dependent manner.

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.

Local Variation in Recombination Rates of the Honey Bee (Apis mellifera) Genome among Samples from Six Disparate Populations

Meiotic recombination is an essential component of eukaryotic sexual reproduction but its frequency varies within and between genomes. Although it is well-established that honey bees have a high recombination rate with about 20 cM/Mbp, the proximate and ultimate causes of this exceptional rate are poorly understood. Here, we describe six linkage maps of the Western Honey Bee Apis mellifera that were produced with consistent methodology from samples from distinct parts of the species' near global distribution. We compared the genome-wide rates and distribution of meiotic crossovers among the six maps and found considerable differences. Overall similarity of local recombination rates among our samples was unrelated to geographic or phylogenetic distance of the populations that our samples were derived from. However, the limited sampling constrains the interpretation of our results because it is unclear how representative these samples are. In contrast to previous studies, we found only in two datasets a significant relation between local recombination rate and GC content. Focusing on regions of particularly increased or decreased recombination in specific maps, we identified several enriched gene ontologies in these regions and speculate about their local adaptive relevance. These data are contributing to an increasing comparative effort to gain an understanding of the intra-specific variability of recombination rates and their evolutionary role in honey bees and other social insects.

The evolutionary advantage of fitness-dependent recombination in diploids: A deterministic mutation-selection balance model

Recombination's omnipresence in nature is one of the most intriguing problems in evolutionary biology. The question of why recombination exhibits certain general features is no less interesting than that of why it exists at all. One such feature is recombination's fitness dependence (FD). The so far developed population genetics models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower- versus higher-fitness genotypes, which can be manifested in selected systems with at least three loci. We also show that if the modifier is linked to such selected system, it can additionally benefit from modifying this linkage in a fitness-dependent manner. The revealed evolutionary advantage of FD recombination appeared robust to crossover interference within the selected system, either positive or negative. Remarkably, FD recombination was often favored in situations where any constant nonzero recombination was evolutionarily disfavored, implying a relaxation of the rather strict constraints on major parameters (e.g., selection intensity and epistasis) required for the evolutionary advantage of nonzero recombination formulated by classical models.

An ultra-high-density genetic map provides insights into genome synteny, recombination landscape and taproot skin colour in radish (Raphanus sativus L.)

High-density genetic map is a valuable tool for exploring novel genomic information, quantitative trait locus (QTL) mapping and gene discovery of economically agronomic traits in plant species. However, high-resolution genetic map applied to tag QTLs associated with important traits and to investigate genomic features underlying recombination landscape in radish (Raphanus sativus) remains largely unexplored. In this study, an ultra-high-density genetic map with 378 738 SNPs covering 1306.8 cM in nine radish linkage groups (LGs) was developed by a whole-genome sequencing-based approach. A total of 18 QTLs for 11 horticulture traits were detected. The map-based cloning data indicated that the R2R3-MYB transcription factor RsMYB90 was a crucial candidate gene determining the taproot skin colour. Comparative genomics analysis among radish, Brassica rapa and B. oleracea genome revealed several genomic rearrangements existed in the radish genome. The highly uneven distribution of recombination was observed across the nine radish chromosomes. Totally, 504 recombination hot regions (RHRs) were enriched near gene promoters and terminators. The recombination rate in RHRs was positively correlated with the density of SNPs and gene, and GC content, respectively. Functional annotation indicated that genes within RHRs were mainly involved in metabolic process and binding. Three QTLs for three traits were found in the RHRs. The results provide novel insights into the radish genome evolution and recombination landscape, and facilitate the development of effective strategies for molecular breeding by targeting and dissecting important traits in radish.

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.

Positive Selection and Functional Divergence at Meiosis Genes That Mediate Crossing Over Across the Drosophila Phylogeny

Meiotic crossing over ensures proper segregation of homologous chromosomes and generates genotypic diversity. Despite these functions, little is known about the genetic factors and population genetic forces involved in the evolution of recombination rate differences among species. The dicistronic meiosis gene, mei-217/mei-218, mediates most of the species differences in crossover rate and patterning during female meiosis between the closely related fruitfly species, Drosophila melanogaster and D. mauritiana The MEI-218 protein is one of several meiosis-specific mini-chromosome maintenance (mei-MCM) proteins that form a multi-protein complex essential to crossover formation, whereas the BLM helicase acts as an anti-crossover protein. Here we study the molecular evolution of five genes- mei-218, the other three known members of the mei-MCM complex, and Blm- over the phylogenies of three Drosophila species groups- melanogaster, obscura, and virilis We then use transgenic assays in D. melanogaster to test if molecular evolution at mei-218 has functional consequences for crossing over using alleles from the distantly related species D. pseudoobscura and D. virilis Our molecular evolutionary analyses reveal recurrent positive selection at two mei-MCM genes. Our transgenic assays show that sequence divergence among mei-218 alleles from D. melanogaster, D. pseudoobscura, and D. virilis has functional consequences for crossing over. In a D. melanogaster genetic background, the D. pseudoobscura mei-218 allele nearly rescues wildtype crossover rates but alters crossover patterning, whereas the D. virilis mei-218 allele conversely rescues wildtype crossover patterning but not crossover rates. These experiments demonstrate functional divergence at mei-218 and suggest that crossover rate and patterning are separable functions.

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.

Substantial Heritable Variation in Recombination Rate on Multiple Scales in Honeybees and Bumblebees

Meiotic recombination shuffles genetic variation and promotes correct segregation of chromosomes. Rates of recombination vary on several scales, both within genomes and between individuals, and this variation is affected by both genetic and environmental factors. Social insects have extremely high rates of recombination, although the evolutionary causes of this are not known. Here, we estimate rates of crossovers and gene conversions in 22 colonies of the honeybee, Apis mellifera, and 9 colonies of the bumblebee, Bombus terrestris, using direct sequencing of 299 haploid drone offspring. We confirm that both species have extremely elevated crossover rates, with higher rates measured in the highly eusocial honeybee than the primitively social bumblebee. There are also significant differences in recombination rate between subspecies of honeybee. There is substantial variation in genome-wide recombination rate between individuals of both A. mellifera and B. terrestris and the distribution of these rates overlap between species. A large proportion of interindividual variation in recombination rate is heritable, which indicates the presence of variation in trans-acting factors that influence recombination genome-wide. We infer that levels of crossover interference are significantly lower in honeybees compared to bumblebees, which may be one mechanism that contributes to higher recombination rates in honeybees. We also find a significant increase in recombination rate with distance from the centromere, mirrored by methylation differences. We detect a strong transmission bias due to GC-biased gene conversion associated with noncrossover gene conversions. Our results shed light on the mechanistic causes of extreme rates of recombination in social insects and the genetic architecture of recombination rate variation.

DNA Motifs Are Not General Predictors of Recombination in Two Drosophila Sister Species

Meiotic recombination is crucial for chromosomal segregation and facilitates the spread of beneficial and removal of deleterious mutations. Recombination rates frequently vary along chromosomes and Drosophila melanogaster exhibits a remarkable pattern. Recombination rates gradually decrease toward centromeres and telomeres, with a dramatic impact on levels of variation in natural populations. Two close sister species, Drosophila simulans and Drosophila mauritiana do not only have higher recombination rates but also exhibit a much more homogeneous recombination rate that only drops sharply very close to centromeres and telomeres. Because certain sequence motifs are associated with recombination rate variation in D. melanogaster, we tested whether the difference in recombination landscape between D. melanogaster and D. simulans can be explained by the genomic distribution of recombination rate–associated sequence motifs. We constructed the first high-resolution recombination map for D. simulans based on 189 haplotypes from a natural D. simulans population and searched for short sequence motifs linked with higher than average recombination in both sister species. We identified five consensus motifs significantly associated with higher than average chromosome-wide recombination rates in at least one species and present in both. Testing fine resolution associations between motif density and recombination, we found strong and positive associations genome-wide over a range of scales in D. melanogaster, while the results were equivocal in D. simulans. Despite the strong association in D. melanogaster, we did not find a decreasing density of these short-repeat motifs toward centromeres and telomeres. We conclude that the density of recombination-associated repeat motifs cannot explain the large-scale recombination landscape in D. melanogaster, nor the differences to D. simulans. The strong association seen for the sequence motifs in D. melanogaster likely reflects their impact influencing local differences in recombination rates along the genome.

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.

What can genetic association panels tell us about evolutionary processes in insects?

If we are to fully comprehend the evolution of insect diversity at a genomic level we need to understand how natural selection can alter genetically encoded characters within populations. Genetic association panels have the potential to be standard bearers in this endeavour. They enable the mapping of phenotypes to genotypes at unprecedented resolution while simultaneously providing population genomic samples that can be interrogated for the tell-tale signs of selection. Analyses of these panels promise to elucidate the entanglement of gene ontologies, pathways, developmental processes and evolutionary constraints, and inform how these are shaped by adaptation.

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.

Wolbachia Infection Associated with Increased Recombination in Drosophila

Wolbachia is a maternally-transmitted endosymbiotic bacteria that infects a large diversity of arthropod and nematode hosts. Some strains of Wolbachia are parasitic, manipulating host reproduction to benefit themselves, while other strains of Wolbachia exhibit obligate or facultative mutualisms with their host. The effects of Wolbachia on its host are many, though primarily relate to host immune and reproductive function. Here we test the hypothesis that Wolbachia infection alters the frequency of homologous recombination during meiosis. We use D. melanogaster as a model system, and survey recombination in eight wild-derived Wolbachia-infected (strain wMel) and Wolbachia-uninfected strains, controlling for genotype. We measure recombination in two intervals of the genome. Our results indicate that Wolbachia infection is associated with increased recombination in one genomic interval and not the other. The effect of Wolbachia infection on recombination is thus heterogenous across the genome. Our data also indicate a reproductive benefit of Wolbachia infection; infected females show higher fecundity than their uninfected genotypic controls. Given the prevalence of Wolbachia infection in natural populations, our findings suggest that Wolbachia infection is likely to contribute to recombination rate and fecundity variation among individuals in nature.

Genome-Wide Association Study of Circadian Behavior in Drosophila melanogaster

Circadian rhythms influence physiological processes from sleep–wake cycles to body temperature and are controlled by highly conserved cycling molecules. Although the mechanistic basis of the circadian clock has been known for decades, the extent to which circadian rhythms vary in nature and the underlying genetic basis for that variation is not well understood. We measured circadian period (Ʈ) and rhythmicity index in the Drosophila Genetic Reference Panel (DGRP) and observed extensive genetic variation in both. Seven DGRP lines had sexually dimorphic arrhythmicity and one line had an exceptionally long Ʈ. Genome-wide analyses identified 584 polymorphisms in 268 genes. We observed differences among transcripts for nine genes predicted to interact among themselves and canonical clock genes in the long period line and a control. Mutations/RNAi knockdown targeting these genes also affected circadian behavior. Our observations reveal that complex genetic interactions influence high levels of variation in circadian phenotypes.

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.

The Sleep Inbred Panel, a Collection of Inbred Drosophila melanogaster with Extreme Long and Short Sleep Duration

Understanding how genomic variation causes differences in observable phenotypes remains a major challenge in biology. It is difficult to trace the sequence of events originating from genomic variants to changes in transcriptional responses or protein modifications. Ideally, one would conduct experiments with individuals that are at either extreme of the trait of interest, but such resources are often not available. Further, advances in genome editing will enable testing of candidate polymorphisms individually and in combination. Here we have created a resource for the study of sleep with 39 inbred lines of Drosophila-the Sleep Inbred Panel (SIP). SIP lines have stable long- and short-sleeping phenotypes developed from naturally occurring polymorphisms. These lines are fully sequenced, enabling more accurate targeting for genome editing and transgenic constructs. This panel facilitates the study of intermediate transcriptional and proteomic correlates of sleep, and supports genome editing studies to verify polymorphisms associated with sleep duration.

Local Inversion Heterozygosity Alters Recombination throughout the Genome

Crossovers (COs) are formed during meiosis by the repair of programmed DNA double-strand breaks (DSBs) and are required for the proper segregation of chromosomes. More DSBs are made than COs, and the remaining DSBs are repaired as noncrossovers (NCOs). The distribution of recombination events along a chromosome occurs in a stereotyped pattern that is shaped by CO-promoting and CO-suppressing forces, collectively referred to as crossover patterning mechanisms. Chromosome inversions are structural aberrations that, when heterozygous, disrupt the recombination landscape by suppressing crossing over. In Drosophila species, the local suppression of COs by heterozygous inversions triggers an increase in crossing over on freely recombining chromosomes termed the interchromosomal (IC) effect [1, 2]. The molecular mechanism(s) by which heterozygous inversions suppress COs, whether noncrossover gene conversions (NCOGCs) are similarly affected, and what mediates the increase in COs in the rest of the genome remain open questions. By sequencing whole genomes of individual offspring from mothers containing heterozygous inversions, we show that, although COs are suppressed by inversions, NCOGCs occur throughout inversions at higher than wild-type frequencies. We confirm that CO frequency increases on the freely recombining chromosomes, yet CO interference remains intact. Intriguingly, NCOGCs do not increase in frequency on the freely recombining chromosomes and the total number of DSBs is approximately the same per genome. Together, our data show that heterozygous inversions change the recombination landscape by altering the relative proportions of COs and NCOGCs and suggest that DSB fate may be plastic until a CO assurance checkpoint has been satisfied.

Connecting theory and data to understand recombination rate evolution

Meiotic recombination is necessary for successful gametogenesis in most sexually reproducing organisms and is a fundamental genomic parameter, influencing the efficacy of selection and the fate of new mutations. The molecular and evolutionary functions of recombination should impose strong selective constraints on the range of recombination rates. Yet, variation in recombination rate is observed on a variety of genomic and evolutionary scales. In the past decade, empirical studies have described variation in recombination rate within genomes, between individuals, between sexes, between populations and between species. At the same time, theoretical work has provided an increasingly detailed picture of the evolutionary advantages to recombination. Perhaps surprisingly, the causes of natural variation in recombination rate remain poorly understood. We argue that empirical and theoretical approaches to understand the evolution of recombination have proceeded largely independently of each other. Most models that address the evolution of recombination rate were created to explain the evolutionary advantage of recombination rather than quantitative differences in rate among individuals. Conversely, most empirical studies aim to describe variation in recombination rate, rather than to test evolutionary hypotheses. In this Perspective, we argue that efforts to integrate the rich bodies of empirical and theoretical work on recombination rate are crucial to moving this field forward. We provide new directions for the development of theory and the production of data that will jointly close this gap.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

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

What drives the evolution of condition-dependent recombination in diploids? Some insights from simulation modelling

While the evolutionary advantages of non-zero recombination rates have prompted diverse theoretical explanations, the evolution of essential recombination features remains underexplored. We focused on one such feature, the condition dependence of recombination, viewed as the variation in within-generation sensitivity of recombination to external (environment) and/or internal (genotype) conditions. Limited empirical evidence for its existence comes mainly from diploids, whereas theoretical models show that it only easily evolves in haploids. The evolution of condition-dependent recombination can be explained by its advantage for the selected system (indirect effect), or by benefits to modifier alleles, ensuring this strategy regardless of effects on the selected system (direct effect). We considered infinite panmictic populations of diploids exposed to a cyclical two-state environment. Each organism had three selected loci. Examining allele dynamics at a fourth, selectively neutral recombination modifier locus, we frequently observed that a modifier allele conferring condition-dependent recombination between the selected loci displaced the allele conferring the optimal constant recombination rate. Our simulations also confirm the results of theoretical studies showing that condition-dependent recombination cannot evolve in diploids on the basis of direct fitness-dependent effects alone. Therefore, the evolution of condition-dependent recombination in diploids can be driven by indirect effects alone, i.e. by modifier effects on the selected system.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Meiotic Recombination: Mixing It Up in Plants

Meiosis halves diploid chromosome numbers to haploid levels that are essential for sexual reproduction in most eukaryotes. Meiotic recombination ensures the formation of bivalents between homologous chromosomes (homologs) and their subsequent proper segregation. It also results in genetic diversity among progeny that influences evolutionary responses to selection. Moreover, crop breeding depends upon the action of meiotic recombination to rearrange elite traits between parental chromosomes. An understanding of the molecular mechanisms that drive meiotic recombination is important for both fundamental research and practical applications. This review emphasizes advances made during the past 5 years, primarily in Arabidopsis and rice, by summarizing newly characterized genes and proteins and examining the regulatory mechanisms that modulate their action.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

A Genomic Reference Panel for Drosophila serrata

Here we describe a collection of re-sequenced inbred lines of Drosophila serrata, sampled from a natural population situated deep within the species endemic distribution in Brisbane, Australia. D. serrata is a member of the speciose montium group whose members inhabit much of south east Asia and has been well studied for aspects of climatic adaptation, sexual selection, sexual dimorphism, and mate recognition. We sequenced 110 lines that were inbred via 17-20 generations of full-sib mating at an average coverage of 23.5x with paired-end Illumina reads. 15,228,692 biallelic SNPs passed quality control after being called using the Joint Genotyper for Inbred Lines (JGIL). Inbreeding was highly effective and the average levels of residual heterozygosity (0.86%) were well below theoretical expectations. As expected, linkage disequilibrium decayed rapidly, with r² dropping below 0.1 within 100 base pairs. With the exception of four closely related pairs of lines which may have been due to technical errors, there was no statistical support for population substructure. Consistent with other endemic populations of other Drosophila species, preliminary population genetic analyses revealed high nucleotide diversity and, on average, negative Tajima’s D values. A preliminary GWAS was performed on a cuticular hydrocarbon trait, 2-Me-C₂₈ revealing 4 SNPs passing Bonferroni significance residing in or near genes. One gene Cht9 may be involved in the transport of CHCs from the site of production (oenocytes) to the cuticle. Our panel will facilitate broader population genomic and quantitative genetic studies of this species and serve as an important complement to existing D. melanogaster panels that can be used to test for the conservation of genetic architectures across the Drosophila genus.