Impact of mating systems on patterns of sequence polymorphism in flowering plants

A fundamental challenge in population genetics and molecular evolution is to understand the forces shaping the patterns of genetic diversity within and among species. Among them, mating systems are thought to have important influences on molecular diversity and genome evolution. Selfing is expected to reduce effective population size, Ne, and effective recombination rates, directly leading to reduced polymorphism and increased linkage disequilibrium compared with outcrossing. Increased isolation between populations also results directly from selfing or indirectly from evolutionary changes, such as small flowers and low pollen output, leading to greater differentiation of molecular markers than under outcrossing. The lower effective recombination rate increases the likelihood of hitch-hiking, further reducing within-deme diversity of selfers and thus increasing their genetic differentiation. There are also indirect effects on molecular evolutionary processes. Low Ne reduces the efficacy of selection; in selfers, selection should thus be less efficient in removing deleterious mutations. The rarity of heterozygous sites in selfers leads to infrequent action of biased conversion towards GC, which tends to increase sequences' GC content in the most highly recombining genome regions of outcrossers. To test these predictions in plants, we used a newly developed sequence polymorphism database to investigate the effects of mating system differences on sequence polymorphism and genome evolution in a wide set of plant species. We also took into account other life-history traits, including life form (whether annual or perennial herbs, and woody perennial) and the modes of pollination and seed dispersal, which are known to affect enzyme and DNA marker polymorphism. We show that among various life-history traits, mating systems have the greatest influence on patterns of polymorphism.     

AN INTEGRATIVE TEST OF THE DEAD‐END HYPOTHESIS OF SELFING EVOLUTION IN TRITICEAE (POACEAE)

Self‐fertilization is hypothesized to be an evolutionary dead end because reversion to outcrossing can rarely happen, and selfing lineages are thought to rapidly become extinct because of limited potential for adaptation and/or accumulation of deleterious mutations. We tested these two assumptions by combining morphological characters and molecular‐evolution analyses in a tribe of hermaphroditic grasses (Triticeae). First, we determined the mating system of the 19 studied species. Then, we sequenced 27 protein‐coding loci and compared base composition and substitution patterns between selfers and outcrossers. We found that the evolution of the mating system is best described by a model including outcrossing‐to‐selfing transitions only. At the molecular level, we showed that regions of low recombination exhibit signatures of relaxed selection. However, we did not detect any evidence of accumulation of nonsynonymous substitutions in selfers compared to outcrossers. Additionally, we tested for the potential deleterious effects of GC‐biased gene conversion in outcrossing species. We found that recombination and not the mating system affected substitution patterns and base composition. We suggest that, in Triticeae, although recombination patterns have remained stable, selfing lineages are of recent origin and inbreeding may have persisted for insufficient time for differences between the two mating systems to evolve.     

Mating systems and the efficacy of selection at the molecular level

Mating systems are thought to play a key role in molecular evolution through their effects on effective population size (N(e)) and effective recombination rate. Because of reduced N(e), selection in self-fertilizing species is supposed to be less efficient, allowing fixation of weakly deleterious alleles or lowering adaptation, which may jeopardize their long-term evolution. Relaxed selection pressures in selfers should be detectable at the molecular level through the analyses of the ratio of nonsynonymous and synonymous divergence, D(n)/D(s), or the ratio of nonsynonymous and synonymous polymorphism, pi(n)/pi(s). On the other hand, selfing reveals recessive alleles to selection (homozygosity effect), which may counterbalance the reduction in N(e). Through population genetics models, this study investigates which process may prevail in natural populations and which conditions are necessary to detect evidence for relaxed selection signature at the molecular level in selfers. Under a wide range of plausible population and mutation parameters, relaxed selection against deleterious mutations should be detectable, but the differences between the two mating systems can be weak. At equilibrium, differences between outcrossers and selfers should be more pronounced using divergence measures (D(n)/D(s) ratio) than using polymorphism data (pi(n)/pi(s) ratio). The difference in adaptive substitution rates between outcrossers and selfers is much less predictable because it critically depends on the dominance levels of new advantageous mutations, which are poorly known. Different ways of testing these predictions are suggested, and implications of these results for the evolution of self-fertilizing species are also discussed.     

Biased gene conversion: implications for genome and sex evolution

Classical genetic studies show that gene conversion can favour some alleles over others. Molecular experiments suggest that gene conversion could favour GC over AT basepairs, leading to the concept of biased gene conversion towards GC (BGC(GC)). The expected consequence of such a process is the GC-enrichment of DNA sequences under gene conversion. Recent genomic work suggests that BGC(GC) affects the base composition of yeast, invertebrate and mammalian genomes. Hypotheses for the mechanisms and evolutionary origin of such a strange phenomenon have been proposed. Most BGC(GC) events probably occur during meiosis, which has implications for our understanding of the evolution of sex and recombination.     

The effects of Hill-Robertson interference between weakly selected mutations on patterns of molecular evolution and variation

Associations between selected alleles and the genetic backgrounds on which they are found can reduce the efficacy of selection. We consider the extent to which such interference, known as the Hill-Robertson effect, acting between weakly selected alleles, can restrict molecular adaptation and affect patterns of polymorphism and divergence. In particular, we focus on synonymous-site mutations, considering the fate of novel variants in a two-locus model and the equilibrium effects of interference with multiple loci and reversible mutation. We find that weak selection Hill-Robertson (wsHR) interference can considerably reduce adaptation, e.g., codon bias, and, to a lesser extent, levels of polymorphism, particularly in regions of low recombination. Interference causes the frequency distribution of segregating sites to resemble that expected from more weakly selected mutations and also generates specific patterns of linkage disequilibrium. While the selection coefficients involved are small, the fitness consequences of wsHR interference across the genome can be considerable. We suggest that wsHR interference is an important force in the evolution of nonrecombining genomes and may explain the unexpected constancy of codon bias across species of very different census population sizes, as well as several unusual features of codon usage in Drosophila.     

Species and recombination effects on DNA variability in the tomato genus

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.     

Microsatellite variation and evolution in the Mimulus guttatus species complex with contrasting mating systems

Mutational variability at microsatellite loci is shaped by both population history and the mating system. In turn, alternate mating systems in flowering plants can resolve aspects of microsatellite loci evolution. Five species of yellow monkeyflowers (Mimulus sect. Simiolis) differing for historical rates of inbreeding were surveyed for variation at six microsatellite loci. High levels of diversity at these loci were found in both outcrossing and selfing taxa. In line with allozyme studies, inbreeders showed more partitioning of diversity among populations, and diversity in selfing taxa was lower than expected from reductions in effective population size due to selfing alone, suggesting the presence of either population bottlenecks or background selection in selfers. Evaluation of the stepwise mutation model (a model of DNA replication slippage) suggests that these loci evolve in a stepwise fashion. Inferred coalescent times of microsatellite alleles indicate that past bottlenecks of population size or colonization events were important in reducing diversity in the inbreeding taxon.      

Direct and indirect consequences of meiotic recombination: implications for genome evolution

There is considerable variation within eukaryotic genomes in the local rate of crossing over. Why is this and what effect does it have on genome evolution? On the genome scale, it is known that by shuffling alleles, recombination increases the efficacy of selection. By contrast, the extent to which differences in the recombination rate modulate the efficacy of selection between genomic regions is unclear. Recombination also has direct consequences on the origin and fate of mutations: biased gene conversion and other forms of meiotic drive promote the fixation of mutations in a similar way to selection, and recombination itself may be mutagenic. Consideration of both the direct and indirect effects of recombination is necessary to understand why its rate is so variable and for correct interpretation of patterns of genome evolution.     

GC-biased gene conversion and selection affect GC content in the Oryza genus (rice)

Base composition varies among and within eukaryote genomes. Although mutational bias and selection have initially been invoked, more recently GC-biased gene conversion (gBGC) has been proposed to play a central role in shaping nucleotide landscapes, especially in yeast, mammals, and birds. gBGC is a kind of meiotic drive in favor of G and C alleles, associated with recombination. Previous studies have also suggested that gBGC could be at work in grass genomes. However, these studies were carried on third codon positions that can undergo selection on codon usage. As most preferred codons end in G or C in grasses, gBGC and selection can be confounded. Here we investigated further the forces that might drive GC content evolution in the rice genus using both coding and noncoding sequences. We found that recombination rates correlate positively with equilibrium GC content and that selfing species (Oryza sativa and O. glaberrima) have significantly lower equilibrium GC content compared with more outcrossing species. As recombination is less efficient in selfing species, these results suggest that recombination drives GC content. We also detected a positive relationship between expression levels and GC content in third codon positions, suggesting that selection favors codons ending with G or C bases. However, the correlation between GC content and recombination cannot be explained by selection on codon usage alone as it was also observed in noncoding positions. Finally, analyses of polymorphism data ruled out the hypothesis that genomic variation in GC content is due to mutational processes. Our results suggest that both gBGC and selection on codon usage affect GC content in the Oryza genus and likely in other grass species.     

Integrating genomics,bioinformatics, and classical genetics to study the effects of recombination on genome evolution

This study presents compelling evidence that recombination significantly increases the silent GC content of a genome in a selectively neutral manner, resulting in a highly significant positive correlation between recombination and "GC3s" in the yeast Saccharomyces cerevisiae. Neither selection nor mutation can explain this relationship. A highly significant GC-biased mismatch repair system is documented for the first time in any member of the Kingdom Fungi. Much of the variation in the GC3s within yeast appears to result from GC-biased gene conversion. Evidence suggests that GC-biased mismatch repair exists in numerous organisms spanning six kingdoms. This transkingdom GC mismatch repair bias may have evolved in response to a ubiquitous AT mutational bias. A significant positive correlation between recombination and GC content is found in many of these same organisms, suggesting that the processes influencing the evolution of the yeast genome may be a general phenomenon. Nonrecombining regions of the genome and nonrecombining genomes would not be subject to this type of molecular drive. It is suggested that the low GC content characteristic of many nonrecombining genomes may be the result of three processes (1) a prevailing AT mutational bias, (2) random fixation of the most common types of mutation, and (3) the absence of the GC-biased gene conversion which, in recombining organisms, permits the reversal of the most common types of mutation. A model is proposed to explain the observation that introns, intergenic regions, and pseudogenes typically have lower GC content than the silent sites of corresponding open reading frames. This model is based on the observation that the greater the heterology between two sequences, the less likely it is that recombination will occur between them. According to this "Constraint" hypothesis, the formation and propagation of heteroduplex DNA is expected to occur, on average, more frequently within conserved coding and regulatory regions of the genome. In organisms possessing GC-biased mismatch repair, this would enhance the GC content of these regions through biased gene conversion. These findings have a number of important implications for the way we view genome evolution and suggest a new model for the evolution of sex.     

植物交酸系统的进化、资源分配对策与遗传多样性

影响植物自交率进化的选择力量主要体现在两个方面：当外来花粉量不足时,自交可以提高植物的结实率,即雌性适合度（繁殖保障）;而如果进行自交的花粉比异交花粉更易获得使胚珠受精的机会,那么自交也可以提高植物的雄性适合度（自动选择优势）。但是,鉴别什么时候是繁殖保障、什么时候是自动选择优势导致了自交的进化却是极其困难的。花粉贴现降低了自交植物通过异交花粉途径获得的适合度,即减弱了自动选择优势,而近交衰退既减少了自动选择优势也减少了繁残给自交者带来的利益。具有不同交配系统的植物种群将具有不同的资源分配对策。理论研究已经说明,自交率增加将减少植物对雄性功能的资源分配比例,但将使繁殖分配加大,而且在一定条件下交配系统在改变甚至可以导致植物生活史发生剧烈变化,即从多年生变为一年生。文献中支持自交减少植物雄性投入的证据有很多,但是对繁殖分配与自交率的关系目前还没有系统的研究,资源分配理论可以解释植物繁育系统的多样性,尤其是能够3说明为什么大多数植物都是雌雄同体的,自交对植物种群遗传结构的影响是减少种群内的遗传变异,增加种群间的遗传分化,长期以来人们一直猜测,自交者可能会丢掉一些长期进化的潜能,目前这个假说得到了一些支持。     

Emergence and Maintenance of Sex among Diploid Organisms Aided by Assortative Mating

Using computer simulations I studied the simultaneous effect of variable environments, mutation rates, ploidy, number of loci subject to evolution and random and assortative mating on various reproductive systems. The simulations showed that mutants for sex and recombination are evolutionarily stable, displacing alleles for monosexuality in diploid populations mating assortatively under variable selection pressure. Assortative mating reduced excessive allelic variance induced by recombination and sex, especially among diploids. Results suggest a novel adaptive value for sex and recombination. They show that the adaptive value of diploidy and that of the segregation of sexes is different to that of sex and recombination. The results suggest that the emergence of sex had to be preceded by the emergence of diploid monosexual organisms and provide an explanation for the emergence and maintenance of sex among diploids and for the scarcity of sex among haploid organisms.     

The advantages of segregation and the evolution of sex

In diploids, sexual reproduction promotes both the segregation of alleles at the same locus and the recombination of alleles at different loci. This article is the first to investigate the possibility that sex might have evolved and been maintained to promote segregation, using a model that incorporates both a general selection regime and modifier alleles that alter an individual's allocation to sexual vs. asexual reproduction. The fate of different modifier alleles was found to depend strongly on the strength of selection at fitness loci and on the presence of inbreeding among individuals undergoing sexual reproduction. When selection is weak and mating occurs randomly among sexually produced gametes, reductions in the occurrence of sex are favored, but the genome-wide strength of selection is extremely small. In contrast, when selection is weak and some inbreeding occurs among gametes, increased allocation to sexual reproduction is expected as long as deleterious mutations are partially recessive and/or beneficial mutations are partially dominant. Under strong selection, the conditions under which increased allocation to sex evolves are reversed. Because deleterious mutations are typically considered to be partially recessive and weakly selected and because most populations exhibit some degree of inbreeding, this model predicts that higher frequencies of sex would evolve and be maintained as a consequence of the effects of segregation. Even with low levels of inbreeding, selection is stronger on a modifier that promotes segregation than on a modifier that promotes recombination, suggesting that the benefits of segregation are more likely than the benefits of recombination to have driven the evolution of sexual reproduction in diploids.     

植物交配系统的进化、资源分配对策与遗传多样性

影响植物自交率进化的选择力量主要体现在两个方面：当外来花粉量不足时,自交可以提高植物的结实率,即雌性适合度（繁殖保障）;而如果进行自交的花粉比异交花粉更易获得使胚珠受精的机会,那么自交也可以提高植物的雄性适合度（自动选择优势）。但是,鉴别什么时候是繁殖保障、什么时候是自动选择优势导致了自交的进化却是极其困难的。花粉贴现降低了自交植物通过异交花粉途径获得的适合度,即减弱了自动选择优势,而近交衰退既减少了自动选择优势也减少了繁殖保障给自交者带来的利益。具有不同交配系统的植物种群将具有不同的资源分配对策。理论研究已经说明,自交率增加将减少植物对雄性功能的资源分配比例,但将使繁殖分配加大,而且在一定条件下交配系统的改变甚至可以导致植物生活史发生剧烈变化,即从多年生变为一年生。文献中支持自交减少植物雄性投入的证据有很多,但是对繁殖分配与自交率的关系目前还没有系统的研究。资源分配理论可以解释植物繁育系统的多样性,尤其是能够说明为什么大多数植物都是雌雄同体的。自交对植物种群遗传结构的影响是减少种群内的遗传变异,增加种群间的遗传分化。长期以来人们一直猜测,自交者可能会丢掉一些长期进化的潜能,目前这个假说得到了一些支持。     

Population,evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.     

Assortative mating for fitness and the evolution of recombination

To understand selection on recombination, we need to consider how linkage disequilibria develop and how recombination alters these disequilibria. Any factor that affects the development of disequilibria, including nonrandom mating, can potentially change selection on recombination. Assortative mating is known to affect linkage disequilibria but its effects on the evolution of recombination have not been previously studied. Given that assortative mating for fitness can arise indirectly via a number of biologically realistic scenarios, it is plausible that weak assortative mating occurs across a diverse set of taxa. Using a modifier model, we examine how assortative mating for fitness affects the evolution of recombination under two evolutionary scenarios: selective sweeps and mutation-selection balance. We find there is no net effect of assortative mating during a selective sweep. In contrast, assortative mating could have a large effect on recombination when deleterious alleles are maintained at mutation-selection balance but only if assortative mating is sufficiently strong. Upon considering reasonable values for the number of loci affecting fitness components, the strength of selection, and the mutation rate, we conclude that the correlation in fitness between mates is unlikely to be sufficiently high for assortative mating to affect the evolution of recombination in most species.     

Epistasis,inbreeding depression,and the evolution of self-fertilization

Inbreeding depression resulting from partially recessive deleterious alleles is thought to be the main genetic factor preventing self-fertilizing mutants from spreading in outcrossing hermaphroditic populations. However, deleterious alleles may also generate an advantage to selfers in terms of more efficient purging, while the effects of epistasis among those alleles on inbreeding depression and mating system evolution remain little explored. In this article, we use a general model of selection to disentangle the effects of different forms of epistasis (additive-by-additive, additive-by-dominance, and dominance-by-dominance) on inbreeding depression and on the strength of selection for selfing. Models with fixed epistasis across loci, and models of stabilizing selection acting on quantitative traits (generating distributions of epistasis) are considered as special cases. Besides its effects on inbreeding depression, epistasis may increase the purging advantage associated with selfing (when it is negative on average), while the variance in epistasis favors selfing through the generation of linkage disequilibria that increase mean fitness. Approximations for the strengths of these effects are derived, and compared with individual-based simulation results.     

The demography and population genomics of evolutionary transitions to self-fertilization in plants

The evolution of self-fertilization from outcrossing has occurred on numerous occasions in flowering plants. This shift in mating system profoundly influences the morphology, ecology, genetics and evolution of selfing lineages. As a result, there has been sustained interest in understanding the mechanisms driving the evolution of selfing and its environmental context. Recently, patterns of molecular variation have been used to make inferences about the selective mechanisms associated with mating system transitions. However, these inferences can be complicated by the action of linked selection following the transition. Here, using multilocus simulations and comparative molecular data from related selfers and outcrossers, we demonstrate that there is little evidence for strong bottlenecks associated with initial transitions to selfing, and our simulation results cast doubt on whether it is possible to infer the role of bottlenecks associated with reproductive assurance in the evolution of selfing. They indicate that the effects of background selection on the loss of diversity and efficacy of selection occur rapidly following the shift to high selfing. Future comparative studies that integrate explicit ecological and genomic details are necessary for quantifying the independent and joint effects of selection and demography on transitions to selfing and the loss of genetic diversity.     

Natural selection affects multiple aspects of genetic variation at putatively neutral sites across the human genome

A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.     

Tempo of Degeneration Across Independently Evolved Nonrecombining Regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.