Polygenic selection to a changing optimum under self-fertilisation.

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

Purging due to self-fertilization does not prevent accumulation of expansion load.

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly.

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa.

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

A Similar Genetic Architecture Underlies the Convergent Evolution of the Selfing Syndrome in Capsella.

Whether, and to what extent, phenotypic evolution follows predictable genetic paths remains an important question in evolutionary biology. Convergent evolution of similar characters provides a unique opportunity to address this question. The transition to selfing and the associated changes in flower morphology are among the most prominent examples of repeated evolution in plants. In this study, we take advantage of the independent transitions to self-fertilization in the genus Capsella to compare the similarities between parallel modifications of floral traits and test for genetic and developmental constraints imposed on flower evolution in the context of the selfing syndrome. Capsella rubella and Capsella orientalis emerged independently but evolved almost identical flower characters. Not only is the evolutionary outcome identical but the same developmental strategies underlie the convergent reduction of flower size. This has been associated with convergent evolution of gene expression changes. The transcriptomic changes common to both selfing lineages are enriched in genes with low network connectivity and with organ-specific expression patterns. Comparative genetic mapping also suggests that, at least in the case of petal size evolution, these similarities have a similar genetic basis. Based on these results, we hypothesize that the limited availability of low-pleiotropy paths predetermines closely related species to similar evolutionary outcomes.

Genetic and demographic signatures accompanying the evolution of the selfing syndrome in Daphne kiusiana, an evergreen shrub.

BACKGROUND AND AIMS: The evolution of mating systems from outcrossing to self-fertilization is a common transition in flowering plants. This shift is often associated with the 'selfing syndrome', which is characterized by less visible flowers with functional changes to control outcrossing. In most cases, the evolutionary history and demographic dynamics underlying the evolution of the selfing syndrome remain poorly understood. METHODS: Here, we characterize differences in the demographic genetic consequences and associated floral-specific traits between two distinct geographical groups of a wild shrub, Daphne kiusiana, endemic to East Asia; plants in the eastern region (southeastern Korea and Kyushu, Japan) exhibit smaller and fewer flowers compared to those of plants in the western region (southwestern Korea). Genetic analyses were conducted using nuclear microsatellites and chloroplast DNA (multiplexed phylogenetic marker sequencing) datasets. KEY RESULTS: A high selfing rate with significantly increased homozygosity characterized the eastern lineage, associated with lower levels of visibility and herkogamy in the floral traits. The two lineages harboured independent phylogeographical histories. In contrast to the western lineage, the eastern lineage showed a gradual reduction in the effective population size with no signs of a severe bottleneck despite its extreme range contraction during the last glacial period. CONCLUSIONS: Our results suggest that the selfing-associated morphological changes in D. kiusiana are of relatively old origin (at least 100 000 years ago) and were driven by directional selection for efficient self-pollination. We provide evidence that the evolution of the selfing syndrome in D. kiusiana is not strongly associated with a severe population bottleneck.

Evolution of sex ratio through gene loss.

The maintenance of males at intermediate frequencies is an important evolutionary problem. Several species of Caenorhabditis nematodes have evolved a mating system in which selfing hermaphrodites and males coexist. While selfing produces XX hermaphrodites, cross-fertilization produces 50% XO male progeny. Thus, male mating success dictates the sex ratio. Here, we focus on the contribution of the male secreted short (mss) gene family to male mating success, sex ratio, and population growth. The mss family is essential for sperm competitiveness in gonochoristic species, but has been lost in parallel in androdioecious species. Using a transgene to restore mss function to the androdioecious Caenorhabditis briggsae, we examined how mating system and population subdivision influence the fitness of the mss+ genotype. Consistent with theoretical expectations, when mss+ and mss-null (i.e., wild type) genotypes compete, mss+ is positively selected in both mixed-mating and strictly outcrossing situations, though more strongly in the latter. Thus, while sexual mode alone affects the fitness of mss+, it is insufficient to explain its parallel loss. However, in genetically homogenous androdioecious populations, mss+ both increases male frequency and depresses population growth. We propose that the lack of inbreeding depression and the strong subdivision that characterize natural Caenorhabditis populations impose selection on sex ratio that makes loss of mss adaptive after self-fertility evolves.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude.

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Genomic characterization of selfing in the cyclic parthenogen Daphnia magna.

Inbreeding refers to the fusion of related individuals' gametes, with self-fertilization (selfing) being an extreme form of inbreeding-involving gametes produced by the same individual. Selfing is expected to reduce heterozygosity by an average of 50% in one generation; however, little is known about the empirical variation on a genome level surrounding this figure and the factors that affect variation. We selfed genotypes of the cyclic parthenogen Daphnia magna and analysed whole genomes of mothers and selfed offspring, observing the predicted 50% heterozygosity reduction on average. We also saw substantial variation around this value and significant differences among mother-offspring pairs. Crossover analysis confirmed the known trend of recombination occurring more often towards the telomeres. This effect was shown, through simulations, to increase the variance of heterozygosity reduction compared to when a uniform distribution of crossovers was used. Similarly, we simulated inbred line production after several generations of selfing and we observed higher variance in achieved homozygosity when we consider a higher recombination rate towards the telomeres. Our empirical and simulation study highlights that the expected mean values of heterozygosity reduction show remarkable variation, which can help understand, for example, differences among inbred individuals.

Beyond the thale: comparative genomics and genetics of Arabidopsis relatives.

For decades a small number of model species have rightly occupied a privileged position in laboratory experiments, but it is becoming increasingly clear that our knowledge of biology is greatly improved when informed by a broader diversity of species and evolutionary context. Arabidopsis thaliana has been the primary model organism for plants, benefiting from a high-quality reference genome sequence and resources for reverse genetics. However, recent studies have made a group of species also in the Brassicaceae family and closely related to A. thaliana a focal point for comparative molecular, genomic, phenotypic and evolutionary studies. In this Review, we emphasize how such studies complement continued study of the model plant itself, provide an evolutionary perspective and summarize our current understanding of genetic and phenotypic diversity in plants.

The Evolutionary Interplay between Adaptation and Self-Fertilization.

Genome-wide surveys of nucleotide polymorphisms, obtained from next-generation sequencing, have uncovered numerous examples of adaptation in self-fertilizing organisms, especially regarding changes to climate, geography, and reproductive systems. Yet existing models for inferring attributes of adaptive mutations often assume idealized outcrossing populations, which risks mischaracterizing properties of these variants. Recent theoretical work is emphasizing how various aspects of self-fertilization affects adaptation, yet empirical data on these properties are lacking. We review theoretical and empirical studies demonstrating how self-fertilization alters the process of adaptation, illustrated using examples from current sequencing projects. We propose ideas for how future research can more accurately quantify aspects of adaptation in self-fertilizers, including incorporating the effects of standing variation, demographic history, and polygenic adaptation.

Phylogenomic analysis reveals multiple evolutionary origins of selfing from outcrossing in a lineage of heterostylous plants.

Evolutionary transitions from outcrossing to selfing often occur in heterostylous plants. Selfing homostyles originate within distylous populations and frequently evolve to become reproductively isolated species. We investigated this process in 10 species of Primula section Obconicolisteri using phylogenomic approaches and inferred how often homostyly originated from distyly and its consequences for population genetic diversity and floral trait evolution. We estimated phylogenetic relationships and reconstructed character evolution using the whole plastome comprised of 76 protein-coding genes. To investigate mating patterns and genetic diversity we screened 15 microsatellite loci in 40 populations. We compared floral traits among distylous and homostylous populations to determine how phenotypically differentiated homostyles were from their distylous ancestors. Section Obconicolisteri was monophyletic and we estimated multiple independent transitions from distyly to homostyly. High selfing rates characterised homostylous populations and this was associated with reduced genetic diversity. Flower size and pollen production were reduced in homostylous populations, but pollen size was significantly larger in some homostyles than in distylous morphs. Repeated transitions to selfing in section Obconicolisteri are likely to have been fostered by the complex montane environments that species occupy. Unsatisfactory pollinator service is likely to have promoted reproductive assurance in homostyles leading to subsequent population divergence through isolation.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata.

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.

Sibling competition does not magnify inbreeding depression in North American Arabidopsis lyrata.

About half of all angiosperms have some form of molecular self-incompatibility to promote outcrossing. If self-incompatibility breaks down, inbreeding depression (δ) is the main barrier to the evolution of self-fertilisation (selfing). If inbreeding depression is lower than 50% (δ < 0.5), the inherent transmission advantage of selfers should theoretically drive the evolution of selfing. However, this does not always happen in practice. For example, despite frequent breakdowns of self-incompatibility in North American Arabidopsis lyrata, selfing has only evolved in few populations. This is surprising given that previous inbreeding-depression estimates were well below the 0.5 threshold. Here, we test whether this could be due to underestimation of true inbreeding depression in competition-free environments. Specifically, we tested whether direct competition between crossed and selfed siblings magnified inbreeding-depression estimates in A. lyrata. We found that this was neither the case for belowground nor for aboveground biomass. For reproductive traits, there was hardly any significant inbreeding depression regardless of competition. Combined with previous findings that drought stress and inducing defence also did not magnify inbreeding depression, our results suggest that the relatively low estimates of inbreeding depression for biomass are indeed realistic estimates of the true inbreeding depression in North American A. lyrata.

Does the evolution of self-fertilization rescue populations or increase the risk of extinction?

Background and Aims: As a major evolutionary transition in seed plants, the evolution of plant mating systems has been much debated in evolutionary ecology. Over the last 10 years, well-established patterns of evolution have emerged. On the one hand, experimental studies have shown that self-fertilization is likely to evolve in a few generations (microevolution) as a response to rapid environmental change (e.g. pollinator decline), eventually rescuing a population. On the other, phylogenetic studies have demonstrated that repeated evolution towards self-fertilization (macroevolution) leads to a higher risk of lineage extinction and is thus likely to be disadvantageous in the long term. Scope: In either case - the short-term or long-term evolution of self-fertilization (selfing) - these findings indicate that a mating system is not neutral with respect to population or lineage persistence. They also suggest that selfing can have contrasting effects depending on time scale. This raises the question of whether mating system evolution can rescue populations facing environmental change. In this review, empirical and theoretical evidence of the direct and indirect effects of mating systems on population demography and lineage persistence were analysed. A simple theoretical evolutionary rescue model was also developed to investigate the potential for evolutionary rescue through selfing. Key Findings: Demographic studies consistently show a short-term advantage of selfing provided by reproductive assurance, but a long-term disadvantage for selfing lineages, suggesting indirect genomic consequences of selfing (e.g. mutation load and lower adaptability). However, our theoretical evolutionary rescue model found that even in the short term, while mating system evolution can lead to evolutionary rescue, it can also lead to evolutionary suicide, due to the inherent frequency-dependent selection of mating system traits. Conclusions: These findings point to the importance of analysing the demographic consequences of self-fertilization in order to predict the effect of selfing on population persistence as well as take into account the indirect genomic consequences of selfing. The pace at which processes such as inbreeding depression, purging, reproductive assurance and genomic rearrangements occur after the selfing transition is the key to clarifying whether or not selfing will result in evolutionary rescue.

Twenty-four-nucleotide siRNAs produce heritable trans-chromosomal methylation in F1 Arabidopsis hybrids.

Hybrid Arabidopsis plants undergo epigenetic reprogramming producing decreased levels of 24-nt siRNAs and altered patterns of DNA methylation that can affect gene expression. Driving the changes in methylation are the processes trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM). In TCM/TCdM the methylation state of one allele is altered to resemble the other allele. We show that Pol IV-dependent sRNAs are required to establish TCM events. The changes in DNA methylation and the associated changes in sRNA levels in the F1 hybrid can be maintained in subsequent generations and affect hundreds of regions in the F2 epigenome. The inheritance of these altered epigenetic states varies in F2 individuals, resulting in individuals with genetically identical loci displaying different epigenetic states and gene expression profiles. The change in methylation at these regions is associated with the presence of sRNAs. Loci without any sRNA activity can have altered methylation states, suggesting that a sRNA-independent mechanism may also contribute to the altered methylation state of the F1 and F2 generations.

Standing genetic variation in a tissue-specific enhancer underlies selfing-syndrome evolution in Capsella.

Mating system shifts recurrently drive specific changes in organ dimensions. The shift in mating system from out-breeding to selfing is one of the most frequent evolutionary transitions in flowering plants and is often associated with an organ-specific reduction in flower size. However, the evolutionary paths along which polygenic traits, such as size, evolve are poorly understood. In particular, it is unclear how natural selection can specifically modulate the size of one organ despite the pleiotropic action of most known growth regulators. Here, we demonstrate that allelic variation in the intron of a general growth regulator contributed to the specific reduction of petal size after the transition to selfing in the genus Capsella Variation within this intron affects an organ-specific enhancer that regulates the level of STERILE APETALA (SAP) protein in the developing petals. The resulting decrease in SAP activity leads to a shortening of the cell proliferation period and reduced number of petal cells. The absence of private polymorphisms at the causal region in the selfing species suggests that the small-petal allele was captured from standing genetic variation in the ancestral out-crossing population. Petal-size variation in the current out-crossing population indicates that several small-effect mutations have contributed to reduce petal-size. These data demonstrate how tissue-specific regulatory elements in pleiotropic genes contribute to organ-specific evolution. In addition, they provide a plausible evolutionary explanation for the rapid evolution of flower size after the out-breeding-to-selfing transition based on additive effects of segregating alleles.

Patterns of Polymorphism at the Self-Incompatibility Locus in 1,083 Arabidopsis thaliana Genomes.

Although the transition to selfing in the model plant Arabidopsis thaliana involved the loss of the self-incompatibility (SI) system, it clearly did not occur due to the fixation of a single inactivating mutation at the locus determining the specificities of SI (the S-locus). At least three groups of divergent haplotypes (haplogroups), corresponding to ancient functional S-alleles, have been maintained at this locus, and extensive functional studies have shown that all three carry distinct inactivating mutations. However, the historical process of loss of SI is not well understood, in particular its relation with the last glaciation. Here, we took advantage of recently published genomic resequencing data in 1,083 Arabidopsis thaliana accessions that we combined with BAC sequencing to obtain polymorphism information for the whole S-locus region at a species-wide scale. The accessions differed by several major rearrangements including large deletions and interhaplogroup recombinations, forming a set of haplogroups that are widely distributed throughout the native range and largely overlap geographically. "Relict" A. thaliana accessions that directly derive from glacial refugia are polymorphic at the S-locus, suggesting that the three haplogroups were already present when glacial refugia from the last Ice Age became isolated. Interhaplogroup recombinant haplotypes were highly frequent, and detailed analysis of recombination breakpoints suggested multiple independent origins. These findings suggest that the complete loss of SI in A. thaliana involved independent self-compatible mutants that arose prior to the last Ice Age, and experienced further rearrangements during postglacial colonization.

Genome Sequencing Reveals the Origin of the Allotetraploid Arabidopsis suecica.

Polyploidy is an example of instantaneous speciation when it involves the formation of a new cytotype that is incompatible with the parental species. Because new polyploid individuals are likely to be rare, establishment of a new species is unlikely unless polyploids are able to reproduce through self-fertilization (selfing), or asexually. Conversely, selfing (or asexuality) makes it possible for polyploid species to originate from a single individual-a bona fide speciation event. The extent to which this happens is not known. Here, we consider the origin of Arabidopsis suecica, a selfing allopolyploid between Arabidopsis thaliana and Arabidopsis arenosa, which has hitherto been considered to be an example of a unique origin. Based on whole-genome re-sequencing of 15 natural A. suecica accessions, we identify ubiquitous shared polymorphism with the parental species, and hence conclusively reject a unique origin in favor of multiple founding individuals. We further estimate that the species originated after the last glacial maximum in Eastern Europe or central Eurasia (rather than Sweden, as the name might suggest). Finally, annotation of the self-incompatibility loci in A. suecica revealed that both loci carry non-functional alleles. The locus inherited from the selfing A. thaliana is fixed for an ancestral non-functional allele, whereas the locus inherited from the outcrossing A. arenosa is fixed for a novel loss-of-function allele. Furthermore, the allele inherited from A. thaliana is predicted to transcriptionally silence the allele inherited from A. arenosa, suggesting that loss of self-incompatibility may have been instantaneous.