The evolution of recombination in a heterogeneous environment

Local adaptation, recombination, and the fate of neopolyploids

Polyploidy is widely recognized as an important speciation mechanism because it isolates tetraploids from their diploid progenitors. Polyploidy also provides new genetic material that may facilitate adaptive evolution. However, new mutations are more likely to arise after a neopolyploid has already successfully invaded a population. Thus, the role of adaptive forces in establishing a polyploid remains unclear. One solution to this apparent paradox may lie in the capacity of polyploids to suppress recombination among preexisting locally adapted alleles. The local adaptation mechanism requires that spatially heterogeneous selection acts on multiple loci and that gene flow introduces maladapted alleles to the population where the polyploid forms. The mechanism requires neither strong genetic drift nor any intrinsic benefit of genome doubling and can accommodate any mode of gene action. A unique prediction of the mechanism is that adaptive alleles should predate polyploidization, a pattern consistent with observations from a few well-studied polyploids. The mechanism is also consistent with the coexistence of both diploid and tetraploid cytotypes, fitness heterogeneity among independently derived polyploids, and the prevalence of outcrossing among older polyploids. The local adaptation mechanism also makes novel predictions about circumstances favoring polyploid invasions that can be tested using molecular genetic or comparative approaches.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

Outcrossing in Caenorhabditis elegans increases in response to food limitation

Theory predicts that organisms should diversify their offspring when faced with a stressful environment. This prediction has received empirical support across diverse groups of organisms and stressors. For example, when encountered by Caenorhabditis elegans during early development, food limitation (a common environmental stressor) induces the nematodes to arrest in a developmental stage called dauer and to increase their propensity to outcross when they are subsequently provided with food and enabled to develop to maturity. Here we tested whether food limitation first encountered during late development/early adulthood can also induce increased outcrossing propensity in C. elegans. Previously well-fed C. elegans increased their propensity to outcross when challenged with food limitation during the final larval stage of development and into early adulthood, relative to continuously well-fed (control) nematodes. Our results thus support previous research demonstrating that the stress of food limitation can induce increased outcrossing propensity in C. elegans. Furthermore, our results expand on previous work by showing that food limitation can still increase outcrossing propensity even when it is not encountered until late development, and this can occur independently of the developmental and gene expression changes associated with dauer.

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

The genomic consequences of selection across development

Understanding how natural selection drives diversification in nature has been at the forefront of biological research for over a century. The main idea is simple: natural selection favours individuals best suited to pass on their genes. However, the journey from birth to reproduction is complex as organisms experience multiple developmental stages, each influenced by genetic and environmental factors (Orr, 2009). These complexities compound even further as each stage of development might be governed by a unique underlying set of alleles and genes. In this issue of Molecular Ecology, Goebl et al. (2022) examine the role of natural selection in driving ecotypic divergence across different life history stages of the prairie sunflower Helianthus petiolaris. The authors used reciprocal transplant experiments, demographic models, and genomic sequencing to explore fitness variation across developmental stages. They show how natural selection impacts population divergence across multiple life history stages and evaluate the resulting allele frequency changes. Goebl et al. link these results to the role of chromosomal inversions, thus furthering our understanding of how ecological divergence proceeds in the face of gene flow. Below, we explore these results in detail and complement their interpretation by considering the evolution of genetic correlations amongst traits governing fitness.

The role of recombination landscape in species hybridisation and speciation

It is now well recognised that closely related species can hybridize and exchange genetic material, which may promote or oppose adaptation and speciation. In some cases, interspecific hybridisation is very common, making it surprising that species identity is preserved despite active gene exchange. The genomes of most eukaryotic species are highly heterogeneous with regard to gene density, abundance of repetitive DNA, chromatin compactisation etc, which can make certain genomic regions more prone or more resistant to introgression of genetic material from other species. Heterogeneity in local recombination rate underpins many of the observed patterns across the genome (e.g. actively recombining regions are typically gene rich and depleted for repetitive DNA) and it can strongly affect the permeability of genomic regions to interspecific introgression. The larger the region lacking recombination, the higher the chance for the presence of species incompatibility gene(s) in that region, making the entire non- or rarely recombining block impermeable to interspecific introgression. Large plant genomes tend to have highly heterogeneous recombination landscape, with recombination frequently occurring at the ends of the chromosomes and central regions lacking recombination. In this paper we review the relationship between recombination and introgression in plants and argue that large rarely recombining regions likely play a major role in preserving species identity in actively hybridising plant species.

Quantitative trait loci mapping reveals an oligogenic architecture of a rapidly adapting trait during the European invasion of common ragweed

Biological invasions offer a unique opportunity to investigate evolution over contemporary timescales. Rapid adaptation to local climates during range expansion can be a major determinant of invasion success, yet fundamental questions remain about its genetic basis. This study sought to investigate the genetic basis of climate adaptation in invasive common ragweed (Ambrosia artemisiifolia). Flowering time adaptation is key to this annual species' invasion success, so much so that it has evolved repeated latitudinal clines in size and phenology across its native and introduced ranges despite high gene flow among populations. Here, we produced a high-density linkage map (4493 SNPs) and paired this with phenotypic data from an F2 mapping population (n = 336) to identify one major and two minor quantitative trait loci (QTL) underlying flowering time and height differentiation in this species. Within each QTL interval, several candidate flowering time genes were also identified. Notably, the major flowering time QTL detected in this study was found to overlap with a previously identified haploblock (putative inversion). Multiple genetic maps of this region identified evidence of suppressed recombination in specific genotypes, consistent with inversions. These discoveries support the expectation that a concentrated genetic architecture with fewer, larger, and more tightly linked alleles should underlie rapid local adaptation during invasion, particularly when divergently adapting populations experience high levels of gene flow.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Chasing genetic correlation breakers to stimulate population resilience to climate change

Global climate change introduces new combinations of environmental conditions, which is expected to increase stress on plants. This could affect many traits in multiple ways that are as yet unknown but will likely require the modification of existing genetic relationships among functional traits potentially involved in local adaptation. Theoretical evolutionary studies have determined that it is an advantage to have an excess of recombination events under heterogeneous environmental conditions. Our study, conducted on a population of radiata pine (Pinus radiata D. Don), was able to identify individuals that show high genetic recombination at genomic regions, which potentially include pleiotropic or collocating QTLs responsible for the studied traits, reaching a prediction accuracy of 0.80 in random cross-validation and 0.72 when whole family was removed from the training population and predicted. To identify these highly recombined individuals, a training population was constructed from correlation breakers, created through tandem selection of parents in the previous generation and their consequent mating. Although the correlation breakers showed lower observed heterogeneity possibly due to direct selection in both studied traits, the genomic regions with statistically significant differences in the linkage disequilibrium pattern showed higher level of heretozygosity, which has the effect of decomposing unfavourable genetic correlation. We propose undertaking selection of correlation breakers under current environmental conditions and using genomic predictions to increase the frequency of these ’recombined’ individuals in future plantations, ensuring the resilience of planted forests to changing climates. The increased frequency of such individuals will decrease the strength of the population-level genetic correlations among traits, increasing the opportunity for new trait combinations to be developed in the future.

Pathways to de novo domestication of crop wild relatives

Growing knowledge about crop domestication, combined with increasingly powerful gene-editing toolkits, sets the stage for the continual domestication of crop wild relatives and other lesser-known plant species.

The role of phenotypic plasticity in the establishment of range margins

It has been argued that adaptive phenotypic plasticity may facilitate range expansions over spatially and temporally variable environments. However, plasticity may induce fitness costs. This may hinder the evolution of plasticity. Earlier modelling studies examined the role of plasticity during range expansions of populations with fixed genetic variance. However, genetic variance evolves in natural populations. This may critically alter model outcomes. We ask: how does the capacity for plasticity in populations with evolving genetic variance alter range margins that populations without the capacity for plasticity are expected to attain? We answered this question using computer simulations and analytical approximations. We found a critical plasticity cost above which the capacity for plasticity has no impact on the expected range of the population. Below the critical cost, by contrast, plasticity facilitates range expansion, extending the range in comparison to that expected for populations without plasticity. We further found that populations may evolve plasticity to buffer temporal environmental fluctuations, but only when the plasticity cost is below the critical cost. Thus, the cost of plasticity is a key factor involved in range expansions of populations with the potential to express plastic response in the adaptive trait. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait's genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor "concentrated" genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.

Chromosomal fusions facilitate adaptation to divergent environments in threespine stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.

Negative heterosis for meiotic recombination rate in spermatocytes of the domestic chicken Gallus gallus

Benef its and costs of meiotic recombination are a matter of discussion. Because recombination breaks allele combinations already tested by natural selection and generates new ones of unpredictable f itness, a high recombination rate is generally benef icial for the populations living in a f luctuating or a rapidly changing environment and costly in a stable environment. Besides genetic benef its and costs, there are cytological effects of recombination, both positive and negative. Recombination is necessary for chromosome synapsis and segregation. However, it involves a massive generation of double-strand DNA breaks, erroneous repair of which may lead to germ cell death or various mutations and chromosome rearrangements. Thus, the benef its of recombination (generation of new allele combinations) would prevail over its costs (occurrence of deleterious mutations) as long as the population remains suff iciently heterogeneous. Using immunolocalization of MLH1, a mismatch repair protein, at the synaptonemal complexes, we examined the number and distribution of recombination nodules in spermatocytes of two chicken breeds with high (Pervomai) and low (Russian Crested) recombination rates and their F1 hybrids and backcrosses. We detected negative heterosis for recombination rate in the F1 hybrids. Backcrosses to the Pervomai breed were rather homogenous and showed an intermediate recombination rate. The differences in overall recombination rate between the breeds, hybrids and backcrosses were mainly determined by the differences in the crossing over number in the seven largest macrochromosomes. The decrease in recombination rate in F1 is probably determined by diff iculties in homology matching between the DNA sequences of genetically divergent breeds. The suppression of recombination in the hybrids may impede gene f low between parapatric populations and therefore accelerate their genetic divergence.

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.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Understanding the evolution of multiple drug resistance in structured populations

The evolution of multidrug resistance (MDR) is a pressing public health concern. Yet many aspects, such as the role played by population structure, remain poorly understood. Here, we argue that studying MDR evolution by focusing upon the dynamical equations for linkage disequilibrium (LD) can greatly simplify the calculations, generate more insight, and provide a unified framework for understanding the role of population structure. We demonstrate how a general epidemiological model of MDR evolution can be recast in terms of the LD equations. These equations reveal how the different forces generating and propagating LD operate in a dynamical setting at both the population and metapopulation levels. We then apply these insights to show how the LD perspective: (i) explains equilibrium patterns of MDR, (ii) provides a simple interpretative framework for transient evolutionary dynamics, and (iii) can be used to assess the consequences of different drug prescription strategies for MDR evolution.

The evolution of sex along an environmental gradient

Although temporally changing environments generally favor sex and recombination, the effects of spatial environmental heterogeneity have been less explored. In this article, we use a classical model of adaptation along with an environmental gradient to study the selective forces acting on reproductive mode evolution in the central and marginal parts of the distribution range of a species. The model considers a polygenic trait under stabilizing selection (the optimal trait value changing across space) and includes a demographic component imposing range limits. The results show that in the central part of the range (where populations are well adapted), recombination tends to increase the mean fitness of offspring in regimes where drift is sufficiently strong (generating a benefit for sex), while it has the opposite effect when the effect of drift stays negligible. However, these effects remain weak and are easily overwhelmed by slight intrinsic fitness differences between sexuals and asexuals. In agreement with previous results, asexuality may be favored in marginal populations, as it can preserve adaptation to extreme conditions. However, a substantial advantage of asexuality is possible only in conditions maintaining a strong maladaptation of sexuals at range limits (high effective environmental gradient, weak selection at loci coding for the trait).

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Genome expansion in early eukaryotes drove the transition from lateral gene transfer to meiotic sex

Prokaryotes acquire genes from the environment via lateral gene transfer (LGT). Recombination of environmental DNA can prevent the accumulation of deleterious mutations, but LGT was abandoned by the first eukaryotes in favour of sexual reproduction. Here we develop a theoretical model of a haploid population undergoing LGT which includes two new parameters, genome size and recombination length, neglected by previous theoretical models. The greater complexity of eukaryotes is linked with larger genomes and we demonstrate that the benefit of LGT declines rapidly with genome size. The degeneration of larger genomes can only be resisted by increases in recombination length, to the same order as genome size - as occurs in meiosis. Our results can explain the strong selective pressure towards the evolution of sexual cell fusion and reciprocal recombination during early eukaryotic evolution - the origin of meiotic sex.

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.

Evolutionary potential for genomic islands of sexual divergence on recombining sex chromosomes

Differentiated sex chromosomes are thought to develop through the accumulation of polymorphisms at loci subject to opposing selection between males and females, and/or between haploids and diploids. As sex chromosomes differentiate, reduced recombination becomes favored between selected loci and the sex-determining region, strengthening genetic associations between alleles favored in a sex and the corresponding sex chromosome. Here a model is analyzed to explore whether polymorphism at one sexually or ploidally antagonistic locus facilitates the spread of rare alleles at other loci experiencing antagonistic selection, promoting further differentiation of the sex chromosomes. It is found that antagonistic polymorphisms can spread and capture other such loci, building 'genomic islands' of differentiation on sex chromosomes, but the conditions are very restrictive, requiring the loci to be strongly selected, tightly linked and distant from the sex-determining region. Epistatic interactions can facilitate the promotion of polymorphism among selected loci, but only if preferentially favoring heterozygotes. Although these results apply to any taxa, plants provide a fertile ground for testing these and related theories given the recurrent evolutionary transitions to dioecy, which provide multiple opportunities to track the early evolution of sex chromosomes.

Mating system impacts the genetic architecture of adaptation to heterogeneous environments

Self-fertilisation has consequences for variation across the genome as it reduces effective population size, effect recombination rates and pollen flow, with implications for local adaptation. We conducted simulations of divergent stabilising selection on a quantitative trait with drift, pollen flow, mutation, recombination and different outcrossing rates. We quantified trait divergence and the genetic architecture of adaptation. We conducted an F_ST outlier analysis to identify candidate loci and quantified the impact of mating system on detectability. Selfing promoted trait divergence mainly through reductions in pollen flow. Moreover, trait architecture became more diffuse with selfing. Average effect size of trait loci was lower, while the number of loci, and their clustering distance increased. The genetic architecture of selfers was also more diffuse than outcrossers for equivalent migration rates. However, when deleterious alleles were included, architectures became more concentrated in selfers, likely to be because of reductions in population size caused by mutational meltdown and impacts of background selection on N_e . Our simulations demonstrate that mating system has important impacts on adaptive divergence of traits and the genetic landscape underlying that divergence. Selfing has a significant effect on detectability of regions of the genome important for adaptation because of neutral divergence and diffuse trait architecture.

Genomic Patterns of Local Adaptation under Gene Flow in Arabidopsis lyrata

Short-scale local adaptation is a complex process involving selection, migration, and drift. The expected effects on the genome are well grounded in theory but examining these on an empirical level has proven difficult, as it requires information about local selection, demographic history, and recombination rate variation. Here, we use locally adapted and phenotypically differentiated Arabidopsis lyrata populations from two altitudinal gradients in Norway to test these expectations at the whole-genome level. Demography modeling indicates that populations within the gradients diverged <2 kya and that the sites are connected by gene flow. The gene flow estimates are, however, highly asymmetric with migration from high to low altitudes being several times more frequent than vice versa. To detect signatures of selection for local adaptation, we estimate patterns of lineage-specific differentiation among these populations. Theory predicts that gene flow leads to concentration of adaptive loci in areas of low recombination; a pattern we observe in both lowland-alpine comparisons. Although most selected loci display patterns of conditional neutrality, we found indications of genetic trade-offs, with one locus particularly showing high differentiation and signs of selection in both populations. Our results further suggest that resistance to solar radiation is an important adaptation to alpine environments, while vegetative growth and bacterial defense are indicated as selected traits in the lowland habitats. These results provide insights into genetic architectures and evolutionary processes driving local adaptation under gene flow. We also contribute to understanding of traits and biological processes underlying alpine adaptation in northern latitudes.

Intercontinental karyotype-environment parallelism supports a role for a chromosomal inversion in local adaptation in a seaweed fly

Large chromosomal rearrangements are thought to facilitate adaptation to heterogeneous environments by limiting genomic recombination. Indeed, inversions have been implicated in adaptation along environmental clines and in ecotype specialization. Here, we combine classical ecological studies and population genetics to investigate an inversion polymorphism previously documented in Europe among natural populations of the seaweed fly Coelopa frigida along a latitudinal cline in North America. We test if the inversion is present in North America and polymorphic, assess which environmental conditions modulate the inversion karyotype frequencies, and document the relationship between inversion karyotype and adult size. We sampled nearly 2000 flies from 20 populations along several environmental gradients to quantify associations of inversion frequencies to heterogeneous environmental variables. Genotyping and phenotyping showed a widespread and conserved inversion polymorphism between Europe and America. Variation in inversion frequency was significantly associated with environmental factors, with parallel patterns between continents, indicating that the inversion may play a role in local adaptation. The three karyotypes of the inversion are differently favoured across micro-habitats and represent life-history strategies likely to be maintained by the collective action of several mechanisms of balancing selection. Our study adds to the mounting evidence that inversions are facilitators of adaptation and enhance within-species diversity.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

Modularity of genes involved in local adaptation to climate despite physical linkage

BACKGROUND

Linkage among genes experiencing different selection pressures can make natural selection less efficient. Theory predicts that when local adaptation is driven by complex and non-covarying stresses, increased linkage is favored for alleles with similar pleiotropic effects, with increased recombination favored among alleles with contrasting pleiotropic effects. Here, we introduce a framework to test these predictions with a co-association network analysis, which clusters loci based on differing associations. We use this framework to study the genetic architecture of local adaptation to climate in lodgepole pine, Pinus contorta, based on associations with environments.

RESULTS

We identify many clusters of candidate genes and SNPs associated with distinct environments, including aspects of aridity and freezing, and discover low recombination rates among some candidate genes in different clusters. Only a few genes contain SNPs with effects on more than one distinct aspect of climate. There is limited correspondence between co-association networks and gene regulatory networks. We further show how associations with environmental principal components can lead to misinterpretation. Finally, simulations illustrate both benefits and caveats of co-association networks.

CONCLUSIONS

Our results support the prediction that different selection pressures favor the evolution of distinct groups of genes, each associating with a different aspect of climate. But our results went against the prediction that loci experiencing different sources of selection would have high recombination among them. These results give new insight into evolutionary debates about the extent of modularity, pleiotropy, and linkage in the evolution of genetic architectures.

Turnover in local parasite populations temporarily favors host outcrossing over self-fertilization during experimental evolution

The ubiquity of outcrossing in plants and animals is difficult to explain given its costs relative to self-fertilization. Despite these costs, exposure to changing environmental conditions can temporarily favor outcrossing over selfing. Therefore, recurring episodes of environmental change are predicted to favor the maintenance of outcrossing. Studies of host-parasite coevolution have provided strong support for this hypothesis. However, it is unclear whether multiple exposures to novel parasite genotypes in the absence of coevolution are sufficient to favor outcrossing. Using the nematode Caenorhabditis elegans and the bacterial parasite Serratia marcescens, we studied host responses to parasite turnover. We passaged several replicates of a host population that was well-adapted to the S. marcescens strain Sm2170 with either Sm2170 or one of three novel S. marcescens strains, each derived from Sm2170, for 18 generations. We found that hosts exposed to novel parasites maintained higher outcrossing rates than hosts exposed to Sm2170. Nonetheless, host outcrossing rates declined over time against all but the most virulent novel parasite strain. Hosts exposed to the most virulent novel strain exhibited increased outcrossing rates for approximately 12 generations, but did not maintain elevated levels of outcrossing throughout the experiment. Thus, parasite turnover can transiently increase host outcrossing. These results suggest that recurring episodes of parasite turnover have the potential to favor the maintenance of host outcrossing. However, such maintenance may require frequent exposure to novel virulent parasites, rapid rates of parasite turnover, and substantial host gene flow.

Varied Genomic Responses to Maladaptive Gene Flow and Their Evidence

Adaptation to a local environment often occurs in the face of maladaptive gene flow. In this perspective, I discuss several ideas on how a genome may respond to maladaptive gene flow during adaptation. On the one hand, selection can build clusters of locally adaptive alleles at fortuitously co-localized loci within a genome, thereby facilitating local adaptation with gene flow ('allele-only clustering'). On the other hand, the selective pressure to link adaptive alleles may drive co-localization of the actual loci relevant for local adaptation within a genome through structural genome changes or an evolving intra-genomic crossover rate ('locus clustering'). While the expected outcome is, in both cases, a higher frequency of locally adaptive alleles in some genome regions than others, the molecular units evolving in response to gene flow differ (i.e., alleles versus loci). I argue that, although making this distinction is important, we commonly lack the critical empirical evidence to do so. This is mainly because many current approaches are biased towards detecting local adaptation in genome regions with low crossover rates. The importance of low-crossover genome regions for adaptation with gene flow, such as in co-localizing relevant loci within a genome, thus remains unclear. Future empirical investigations should address these questions by making use of comparative genomics, where multiple de novo genome assemblies from species evolved under different degrees of genetic exchange are compared. This research promises to advance our understanding of how a genome adapts to maladaptive gene flow, thereby promoting adaptive divergence and reproductive isolation.

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

The evolution of recombination rates in finite populations during ecological speciation

Recombination can impede ecological speciation with gene flow by mixing locally adapted genotypes with maladapted migrant genotypes from a divergent population. In such a scenario, suppression of recombination can be selectively favoured. However, in finite populations evolving under the influence of random genetic drift, recombination can also facilitate adaptation by reducing Hill-Robertson interference between loci under selection. In this case, increased recombination rates can be favoured. Although these two major effects on recombination have been studied individually, their joint effect on ecological speciation with gene flow remains unexplored. Using a mathematical model, we investigated the evolution of recombination rates in two finite populations that exchange migrants while adapting to contrasting environments. Our results indicate a two-step dynamic where increased recombination is first favoured (in response to the Hill-Robertson effect), and then disfavoured, as the cost of recombining locally with maladapted migrant genotypes increases over time (the maladaptive gene flow effect). In larger populations, a stronger initial benefit for recombination was observed, whereas high migration rates intensify the long-term cost of recombination. These dynamics may have important implications for our understanding of the conditions that facilitate incipient speciation with gene flow and the evolution of recombination in finite populations.

Contribution of gene flow to the evolution of recombination suppression in sex chromosomes

Polymorphism of alleles that benefit one sex but harm the other (sexually antagonistic alleles) generates selective pressures for reduced recombination between themselves and sex-determination loci. Such polymorphism can be maintained within a population when selection coefficients are sufficiently balanced between males and females. However, if regulatory mutations restrict gene expression only to one sex, these alleles become neutral in the other sex and easily fixed within a population, removing the selective pressures for recombination suppression in sex chromosomes. When there is spatial variation in selection regimes, however, alleles that are deleterious in one sex and neutral in the other can be maintained in other neighboring populations and gene flow may continuously supply deleterious alleles. We hypothesized that this maintenance of genetic variation may promote the establishment of recombination suppression in sex chromosomes even in cases where selection is limited to one sex. Using individual-based simulations, we show that spatial variation in male-limited selection and gene flow can promote the establishment of Y-autosome fusions, a special case of recombination suppression in sex chromosomes. This can be explained by the fact that fused Y-chromosomes that capture alleles that are beneficial for local males have a higher mean fitness compared to unfused Y chromosomes in the presence of deleterious gene flow. We also simulated the case of sex-concordant selection and found that gene flow of alleles that are deleterious in both sexes did not substantially increase the establishment rates of Y-autosome fusions across the parameter space examined. This can be accounted for by the fact that foreign alleles that are deleterious in both sexes can be efficiently removed from the population compared to alleles that are neutral in females. These results indicate that how gene flow affects the establishment rates of Y-autosome fusions depends largely on selection regimes. Spatial variation in sex-specific selection and gene flow should be appreciated as a factor affecting sex chromosome evolution.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Habitat heterogeneity favors asexual reproduction in natural populations of grassthrips

Explaining the overwhelming success of sex among eukaryotes is difficult given the obvious costs of sex relative to asexuality. Different studies have shown that sex can provide benefits in spatially heterogeneous environments under specific conditions, but whether spatial heterogeneity commonly contributes to the maintenance of sex in natural populations remains unknown. We experimentally manipulated habitat heterogeneity for sexual and asexual thrips lineages in natural populations and under seminatural mesocosm conditions by varying the number of hostplants available to these herbivorous insects. Asexual lineages rapidly replaced the sexual ones, independently of the level of habitat heterogeneity in mesocosms. In natural populations, the success of sexual thrips decreased with increasing habitat heterogeneity, with sexual thrips apparently only persisting in certain types of hostplant communities. Our results illustrate how genetic diversity-based mechanisms can favor asexuality instead of sex when sexual lineages co-occur with genetically variable asexual lineages.

Recombination Rate Evolution and the Origin of Species

A recipe for dissolving incipient species into a continuum of phenotypes is to recombine their genetic material. Therefore, students of speciation have become increasingly interested in the mechanisms by which recombination between locally adapted lineages is reduced. Evidence abounds that chromosomal rearrangements, via their suppression of recombination during meiosis in hybrids, play a major role in adaptation and speciation. By contrast, genic modifiers of recombination rates have been largely ignored in studies of speciation. We show how both types of reduction in recombination rates facilitate divergence in the face of gene flow, including the early stages of adaptive divergence, the persistence of species after secondary contact, and reinforcement.

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

BACKGROUND

Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each).

RESULTS

For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination.

CONCLUSIONS

Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

Recombining without Hotspots: A Comprehensive Evolutionary Portrait of Recombination in Two Closely Related Species of Drosophila

Meiotic recombination rate varies across the genome within and between individuals, populations, and species in virtually all taxa studied. In almost every species, this variation takes the form of discrete recombination hotspots, determined in some mammals by a protein called PRDM9. Hotspots and their determinants have a profound effect on the genomic landscape, and share certain features that extend across the tree of life. Drosophila, in contrast, are anomalous in their absence of hotspots, PRDM9, and other species-specific differences in the determination of recombination. To better understand the evolution of meiosis and general patterns of recombination across diverse taxa, we present a truly comprehensive portrait of recombination across time, combining recently published cross-based contemporary recombination estimates from each of two sister species with newly obtained linkage-disequilibrium-based historic estimates of recombination from both of these species. Using Drosophila pseudoobscura and Drosophila miranda as a model system, we compare recombination rate between species at multiple scales, and we suggest that Drosophila replicate the pattern seen in human-chimpanzee in which recombination rate is conserved at broad scales. We also find evidence of a species-wide recombination modifier(s), resulting in both a present and historic genome-wide elevation of recombination rates in D. miranda, and identify broad scale effects on recombination from the presence of an inversion. Finally, we reveal an unprecedented view of the distribution of recombination in D. pseudoobscura, illustrating patterns of linked selection and where recombination is taking place. Overall, by combining these estimation approaches, we highlight key similarities and differences in recombination between Drosophila and other organisms.

Contribution of recombination to the evolutionary history of HIV

PURPOSE OF REVIEW

An improved understanding of how recombination affects the evolutionary history of HIV is crucial to understand its current and future evolution. The present review aims to disentangle the manifold effects of recombination on HIV by discussing its effects on the evolutionary history and the adaptive potential of HIV in the context of concepts from evolutionary genetics and genomics.

RECENT FINDINGS

The increasing occurrence of secondary contacts between divergent subtype populations (during coinfection) results in increased observations of recombinants worldwide. Recombination is heterogeneous along the HIV genome. Consequences of recombination of HIV evolution are, in combination with other demographic processes, expected to either homogenize the genetic composition of HIV populations (homogenization) or provide the potential for novel adaptations (diversification). New methods in population genomics allow deep characterization of recombinant genome (the segment composition and origin) and their evolutionary trajectories.

SUMMARY

HIV recombinants increase worldwide and invade geographical regions where pure subtypes were previously predominant. This trend is expected to continue in the future, as ease to travel worldwide increases opportunities for recombination between divergent HIV strains. While the effects of recombination in HIV are much researched, more effort is required to characterize current HIV recombinant composition and dynamics. This can be achieved with new population genetic and genomic methods.

The evolution of obligate sex: the roles of sexual selection and recombination

The evolution of sex is one of the greatest mysteries in evolutionary biology. An even greater mystery is the evolution of obligate sex, particularly when competing with facultative sex and not with complete asexuality. Here, we develop a stochastic simulation of an obligate allele invading a facultative population, where males are subject to sexual selection. We identify a range of parameters where sexual selection can contribute to the evolution of obligate sex: Especially when the cost of sex is low, mutation rate is high, and the facultative individuals do not reproduce sexually very often. The advantage of obligate sex becomes larger in the absence of recombination. Surprisingly, obligate sex can take over even when the population has a lower mean fitness as a result. We show that this is due to the high success of obligate males that can compensate the cost of sex.