Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

Genome structure and the benefit of sex

We examine the behavior of sexual and asexual populations in modular multipeaked fitness landscapes and show that sexuals can systematically reach different, higher fitness adaptive peaks than asexuals. Whereas asexuals must move against selection to escape local optima, sexuals reach higher fitness peaks reliably because they create specific genetic variants that "skip over" fitness valleys, moving from peak to peak in the fitness landscape. This occurs because recombination can supply combinations of mutations in functional composites or "modules," that may include individually deleterious mutations. Thus when a beneficial module is substituted for another less-fit module by sexual recombination it provides a genetic variant that would require either several specific simultaneous mutations in an asexual population or a sequence of individual mutations some of which would be selected against. This effect requires modular genomes, such that subsets of strongly epistatic mutations are tightly physically linked. We argue that such a structure is provided simply by virtue of the fact that genomes contain many genes each containing many strongly epistatic nucleotides. We briefly discuss the connections with "building blocks" in the evolutionary computation literature. We conclude that there are conditions in which sexuals can systematically evolve high-fitness genotypes that are essentially unevolvable for asexuals.     

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp ) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon.     

Recombination enhances protein adaptation in Drosophila melanogaster

Evolutionary theory predicts that the rate and level of adaptation will be enhanced in sexual relative to asexual genomes because sexual recombination facilitates the elimination of deleterious mutations and the fixation of beneficial ones by natural selection. To date, the most compelling evidence for this prediction comes from experimental evolution studies and from loci completely lacking recombination, such as those on Y chromosomes, which often show reduced adaptation and even degeneration. Here, by analyzing replacement and silent DNA polymorphism and divergence at 98 loci, I show that recombination increases the efficacy of protein adaptation throughout the genome of the fruit fly Drosophila melanogaster. Genes residing in genomic regions with reduced recombination rates suffer a greater load of segregating, mildly deleterious mutations and fix fewer beneficial mutations than genes residing in regions with higher recombination rates. These findings suggest that the capacity to respond to natural selection varies with recombination rate across the genome, consistent with theory on the evolutionary advantages of sex and recombination.     

Mutation Load in Sunflower Inversions Is Negatively Correlated with Inversion Heterozygosity

Recombination is critical both for accelerating adaptation and purging deleterious mutations. Chromosomal inversions can act as recombination modifiers that suppress local recombination in heterozygotes and thus, under some conditions, are predicted to accumulate such mutations. In this study, we investigated patterns of recombination, transposable element abundance, and coding sequence evolution across the genomes of 1,445 individuals from three sunflower species, as well as within nine inversions segregating within species. We also analyzed the effects of inversion genotypes on 87 phenotypic traits to test for overdominance. We found significant negative correlations of long terminal repeat retrotransposon abundance and deleterious mutations with recombination rates across the genome in all three species. However, we failed to detect an increase in these features in the inversions, except for a modest increase in the proportion of stop codon mutations in several very large or rare inversions. Consistent with this finding, there was little evidence of overdominance of inversions in phenotypes that may relate to fitness. On the other hand, significantly greater load was observed for inversions in populations polymorphic for a given inversion compared to populations monomorphic for one of the arrangements, suggesting that the local state of inversion polymorphism affects deleterious load. These seemingly contradictory results can be explained by the low frequency of inversion heterozygotes in wild sunflower populations, apparently due to divergent selection and associated geographic structure. Inversions contributing to local adaptation represent ideal recombination modifiers, acting to facilitate adaptive divergence with gene flow, while largely escaping the accumulation of deleterious mutations.     

Balancing selection maintains cryptic colour morphs

Animals display incredibly diverse colour patterns, a testament to evolution's endless innovation in shaping life. In many species, the interplay between males and females in the pursuit of mates has driven the evolution of a myriad of colour forms, from the flashy peacock tail feathers to the tiniest colour markings in damselflies. In others, colour provides crypsis by allowing to blend into the background and to escape the eyes of predators. While the obvious benefits of this dazzling diversity for reproduction and survival seem straightforward, its maintenance is not. Theory predicts that genetic drift and various forms of selection reduce variation over time, making the persistence of colour variants over generations a puzzle. In this issue of Molecular Ecology, Lindtke et al. ( 2017 ) study the cryptic colour morphs of Timema cristinae walking sticks to shed light on the genetic architecture and mechanisms that allow colour polymorphism maintenance over long timescales. By combining genome‐wide data with phenotyping information from natural populations, they were able to map the green and melanistic colour to one genomic region with highly reduced effective recombination rate between two main chromosomal variants, consistent with an inversion polymorphism. These two main chromosomal variants showed geographically widespread heterozygote excess, and genomic signatures consistent with long‐term balancing selection. A younger chromosomal variant was detected for the third morph, the green‐striped colour morphs, in the same genomic regions as the melanistic and the green‐unstriped morphs. Together, these results suggest that the genetic architecture of cryptic T. cristinae morphs is caused by nonrecombining genomic blocks that have been maintained over extended time periods by balancing selection making this study one of the few available empirical examples documenting that balancing selection of various forms may play an important role in maintaining adaptive genetic variation in nature.     

Genetic linkage and natural selection

The prevalence of recombination in eukaryotes poses one of the most puzzling questions in biology. The most compelling general explanation is that recombination facilitates selection by breaking down the negative associations generated by random drift (i.e. Hill–Robertson interference, HRI). I classify the effects of HRI owing to: deleterious mutation, balancing selection and selective sweeps on: neutral diversity, rates of adaptation and the mutation load. These effects are mediated primarily by the density of deleterious mutations and of selective sweeps. Sequence polymorphism and divergence suggest that these rates may be high enough to cause significant interference even in genomic regions of high recombination. However, neither seems able to generate enough variance in fitness to select strongly for high rates of recombination. It is plausible that spatial and temporal fluctuations in selection generate much more fitness variance, and hence selection for recombination, than can be explained by uniformly deleterious mutations or species-wide selective sweeps.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Molecular genetics of natural populations

Significant progress in evolutionary genetics has been made by studying, on the one hand, patterns of DNA sequence polymorphism and, on the other, genetic architecture of complex adaptive traits. However, connections between nucleotide variants under selection and adaptively relevant phenotypes are missing. Such connections can be established using precise gene replacement. We review the recent successful introduction of this technique to the analysis of two evolutionarily interesting loci--Odysseus and desaturase2. Both genes have subtle phenotypes that nevertheless could be identified using gene replacement, demonstrating that effects of naturally occurring alleles can be measured in the laboratory. This is an important first step in connecting statistical signatures of selection with adaptation in nature. More candidate genes involved in adaptation, for example, through cloning of genes responsible for reproductive isolation, now need to be identified. Molecular genetic manipulation, DNA polymorphism analysis, and field studies then have to be integrated to provide fresh insights into the mechanisms of evolutionary change.     

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI) , PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7) , and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein–protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.     

Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata

Genetic diversity is unusually high at loci in the S-locus region of the self-incompatible species of the flowering plant, Arabidopsis lyrata, not just in the S loci themselves, but also at two nearby loci. In a previous study of a single natural population from Iceland, we attributed this elevated polymorphism to linkage disequilibrium (LD) between variants at loci close to the S locus and the S alleles, which are maintained in the population by balancing selection. With the four S-flanking loci whose diversity we previously studied, we could not determine the extent of the region linked to the S loci in which neutral sites are affected. We also could not exclude the possibility of a population bottleneck, or of admixture, as causes of the LD. We have now studied four more distant loci flanking the S-locus region, and more populations, and we analyze the results using a theoretical model of the effect of balancing selection on diversity at linked neutral sites within and between different functional S-allelic classes. In the model, diversity is a function of the number of selectively maintained alleles and the recombination distances from the selectively maintained sites. We use the model to estimate the number of different functional S alleles, their turnover rate, and recombination rates between the S-locus region and other loci. Our estimates suggest that there is a small region of very low recombination surrounding the S-locus region.     

Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

The long-term evolution of multilocus traits under frequency-dependent disruptive selection

Frequency-dependent disruptive selection is widely recognized as an important source of genetic variation. Its evolutionary consequences have been extensively studied using phenotypic evolutionary models, based on quantitative genetics, game theory, or adaptive dynamics. However, the genetic assumptions underlying these approaches are highly idealized and, even worse, predict different consequences of frequency-dependent disruptive selection. Population genetic models, by contrast, enable genotypic evolutionary models, but traditionally assume constant fitness values. Only a minority of these models thus addresses frequency-dependent selection, and only a few of these do so in a multilocus context. An inherent limitation of these remaining studies is that they only investigate the short-term maintenance of genetic variation. Consequently, the long-term evolution of multilocus characters under frequency-dependent disruptive selection remains poorly understood. We aim to bridge this gap between phenotypic and genotypic models by studying a multilocus version of Levene's soft-selection model. Individual-based simulations and deterministic approximations based on adaptive dynamics theory provide insights into the underlying evolutionary dynamics. Our analysis uncovers a general pattern of polymorphism formation and collapse, likely to apply to a wide variety of genetic systems: after convergence to a fitness minimum and the subsequent establishment of genetic polymorphism at multiple loci, genetic variation becomes increasingly concentrated on a few loci, until eventually only a single polymorphic locus remains. This evolutionary process combines features observed in quantitative genetics and adaptive dynamics models, and it can be explained as a consequence of changes in the selection regime that are inherent to frequency-dependent disruptive selection. Our findings demonstrate that the potential of frequency-dependent disruptive selection to maintain polygenic variation is considerably smaller than previously expected.     

Occasional recombination of a selfish X‐chromosome may permit its persistence at high frequencies in the wild

The sex‐ratio X‐chromosome (SR) is a selfish chromosome that promotes its own transmission to the next generation by destroying Y‐bearing sperm in the testes of carrier males. In some natural populations of the fly Drosophila neotestacea, up to 30% of the X‐chromosomes are SR chromosomes. To investigate the molecular evolutionary history and consequences of SR, we sequenced SR and standard (ST) males at 11 X‐linked loci that span the ST X‐chromosome and at seven arbitrarily chosen autosomal loci from a sample of D. neotestacea males from throughout the species range. We found that the evolutionary relationship between ST and SR varies among individual markers, but genetic differentiation between SR and ST is chromosome‐wide and likely due to large chromosomal inversions that suppress recombination. However, SR does not consist of a single multilocus haplotype: we find evidence for gene flow between ST and SR at every locus assayed. Furthermore, we do not find long‐distance linkage disequilibrium within SR chromosomes, suggesting that recombination occurs in females homozygous for SR. Finally, polymorphism on SR is reduced compared to that on ST, and loci displaying signatures of selection on ST do not show similar patterns on SR. Thus, even if selection is less effective on SR, our results suggest that gene flow with ST and recombination between SR chromosomes may prevent the accumulation of deleterious mutations and allow its long‐term persistence at relatively high frequencies.     

Selection against accumulating mutations in niche-preference genes can drive speciation

Our current understanding of sympatric speciation is that it occurs primarily through disruptive selection on ecological genes driven by competition, followed by reproductive isolation through reinforcement-like selection against inferior intermediates/heterozygotes. Our evolutionary model of selection on resource recognition and preference traits suggests a new mechanism for sympatric speciation. We find speciation can occur in three phases. First a polymorphism of functionally different phenotypes is established through evolution of specialization. On the gene level, regulatory functions have evolved in which some alleles are conditionally switched off (i.e. are silent). These alleles accumulate harmful mutations that potentially may be expressed in offspring through recombination. Second mating associated with resource preference invades because harmful mutations in parents are not expressed in the offspring when mating assortatively, thereby dividing the population into two pre-zygotically isolated resource-specialist lineages. Third, silent alleles that evolved in phase one now accumulate deleterious mutations over the following generations in a Bateson-Dobzhansky-Muller fashion, establishing a post-zygotic barrier to hybridization.     

Self‐fertilization and inbreeding limit the scope for sexually antagonistic polymorphism

Sexual antagonism occurs when there is a positive intersexual genetic correlation in trait expression but opposite fitness effects of the trait(s) in males and females. As such, it constrains the evolution of sexual dimorphism and may therefore have implications for adaptive evolution. There is currently considerable evidence for the existence of sexually antagonistic genetic variation in laboratory and natural populations, but how sexual antagonism interacts with other evolutionary phenomena is still poorly understood in many cases. Here, we explore how self‐fertilization and inbreeding affect the maintenance of polymorphism for sexually antagonistic loci. We expected a priori that selfing should reduce the region of polymorphism, as inbreeding reduces the frequency of heterozygotes and speeds fixation. This expectation was supported, but although previous results suggest that the more an allele that is deleterious to one sex is dominant in that sex, the smaller the region of parameter space that will admit polymorphism, we found that this effect is weakened by self‐fertilization. However, the effect of inbreeding is not strong enough to completely cancel out the effect of dominance: For a given frequency of inbreeding, it will still be the case that the more dominant the alleles are in their deleterious context, the smaller the region of parameter space in which they can exist at polymorphism.     

Microsatellite variation and recombination rate in the human genome

Background (purifying) selection on deleterious mutations is expected to remove linked neutral mutations from a population, resulting in a positive correlation between recombination rate and levels of neutral genetic variation, even for markers with high mutation rates. We tested this prediction of the background selection model by comparing recombination rate and levels of microsatellite polymorphism in humans. Published data for 28 unrelated Europeans were used to estimate microsatellite polymorphism (number of alleles, heterozygosity, and variance in allele size) for loci throughout the genome. Recombination rates were estimated from comparisons of genetic and physical maps. First, we analyzed 61 loci from chromosome 22, using the complete sequence of this chromosome to provide exact physical locations. These 61 microsatellites showed no correlation between levels of variation and recombination rate. We then used radiation-hybrid and cytogenetic maps to calculate recombination rates throughout the genome. Recombination rates varied by more than one order of magnitude, and most chromosomes showed significant suppression of recombination near the centromere. Genome-wide analyses provided no evidence for a strong positive correlation between recombination rate and polymorphism, although analyses of loci with at least 20 repeats suggested a weak positive correlation. Comparisons of microsatellites in lowest-recombination and highest-recombination regions also revealed no difference in levels of polymorphism. Together, these results indicate that background selection is not a major determinant of microsatellite variation in humans.     

Adaptive introgression in animals: examples and comparison to new mutation and standing variation as sources of adaptive variation

Adaptive genetic variation has been thought to originate primarily from either new mutation or standing variation. Another potential source of adaptive variation is adaptive variants from other (donor) species that are introgressed into the (recipient) species, termed adaptive introgression. Here, the various attributes of these three potential sources of adaptive variation are compared. For example, the rate of adaptive change is generally thought to be faster from standing variation, slower from mutation and potentially intermediate from adaptive introgression. Additionally, the higher initial frequency of adaptive variation from standing variation and lower initial frequency from mutation might result in a higher probability of fixation of the adaptive variants for standing variation. Adaptive variation from introgression might have higher initial frequency than new adaptive mutations but lower than that from standing variation, again making the impact of adaptive introgression variation potentially intermediate. Adaptive introgressive variants might have multiple changes within a gene and affect multiple loci, an advantage also potentially found for adaptive standing variation but not for new adaptive mutants. The processes that might produce a common variant in two taxa, convergence, trans‐species polymorphism from incomplete lineage sorting or from balancing selection and adaptive introgression, are also compared. Finally, potential examples of adaptive introgression in animals, including balancing selection for multiple alleles for major histocompatibility complex (MHC), S and csd genes, pesticide resistance in mice, black colour in wolves and white colour in coyotes, Neanderthal or Denisovan ancestry in humans, mimicry genes in Heliconius butterflies, beak traits in Darwin's finches, yellow skin in chickens and non‐native ancestry in an endangered native salamander, are examined.     

Genetic Diversity,Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex,Daphnia pulicaria,and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.     

Regions of lower crossing over harbor more rare variants in African populations of Drosophila melanogaster

A correlation between diversity levels and rates of recombination is predicted both by models of positive selection, such as hitchhiking associated with the rapid fixation of advantageous mutations, and by models of purifying selection against strongly deleterious mutations (commonly referred to as "background selection"). With parameter values appropriate for Drosophila populations, only the first class of models predicts a marked skew in the frequency spectrum of linked neutral variants, relative to a neutral model. Here, we consider 29 loci scattered throughout the Drosophila melanogaster genome. We show that, in African populations, a summary of the frequency spectrum of polymorphic mutations is positively correlated with the meiotic rate of crossing over. This pattern is demonstrated to be unlikely under a model of background selection. Models of weakly deleterious selection are not expected to produce both the observed correlation and the extent to which nucleotide diversity is reduced in regions of low (but nonzero) recombination. Thus, of existing models, hitchhiking due to the recurrent fixation of advantageous variants is the most plausible explanation for the data.     

Background Selection as Baseline for Nucleotide Variation across the Drosophila Genome

The constant removal of deleterious mutations by natural selection causes a reduction in neutral diversity and efficacy of selection at genetically linked sites (a process called Background Selection, BGS). Population genetic studies, however, often ignore BGS effects when investigating demographic events or the presence of other types of selection. To obtain a more realistic evolutionary expectation that incorporates the unavoidable consequences of deleterious mutations, we generated high-resolution landscapes of variation across the Drosophila melanogaster genome under a BGS scenario independent of polymorphism data. We find that BGS plays a significant role in shaping levels of variation across the entire genome, including long introns and intergenic regions distant from annotated genes. We also find that a very large percentage of the observed variation in diversity across autosomes can be explained by BGS alone, up to 70% across individual chromosome arms at 100-kb scale, thus indicating that BGS predictions can be used as baseline to infer additional types of selection and demographic events. This approach allows detecting several outlier regions with signal of recent adaptive events and selective sweeps. The use of a BGS baseline, however, is particularly appropriate to investigate the presence of balancing selection and our study exposes numerous genomic regions with the predicted signature of higher polymorphism than expected when a BGS context is taken into account. Importantly, we show that these conclusions are robust to the mutation and selection parameters of the BGS model. Finally, analyses of protein evolution together with previous comparisons of genetic maps between Drosophila species, suggest temporally variable recombination landscapes and, thus, local BGS effects that may differ between extant and past phases. Because genome-wide BGS and temporal changes in linkage effects can skew approaches to estimate demographic and selective events, future analyses should incorporate BGS predictions and capture local recombination variation across genomes and along lineages.