3 1 l e t t e r sThe shift from outcrossing to selfing is common in flowering plants 1,2 , but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.

Whole-genome duplication (WGD) events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection, and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole-genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole-genome data provide strong support for recent hybrid origins of the tetraploid species within the past 100,000-300,000 y from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide due to the combined effects of demography, selfing, and genome redundancy from WGD. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that nonfunctionalization can receive a "head start" through a legacy of deleterious variants and differential expression originating in parental diploid populations.polyploidy | population genomics | speciation | gene loss

Whole-genome duplication (WGD) events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection, and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole-genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole-genome data provide strong support for recent hybrid origins of the tetraploid species within the past 100,000-300,000 y from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide due to the combined effects of demography, selfing, and genome redundancy from WGD. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that nonfunctionalization can receive a "head start" through a legacy of deleterious variants and differential expression originating in parental diploid populations.polyploidy | population genomics | speciation | gene loss

Genomic imprinting is an epigenetic phenomenon occurring in mammals and flowering plants that causes genes to adopt a parent-of-origin-specific mode of expression. While the imprinting status of genes is well conserved in mammals, clear estimates for the degree of conservation were lacking in plants. We therefore analyzed the genome-wide imprinting status of Capsella rubella, which shared a common recent ancestor with Arabidopsis thaliana ∼10 to 14 million years ago. However, only ∼14% of maternally expressed genes (MEGs) and ∼29% of paternally expressed genes (PEGs) in C. rubella were commonly imprinted in both species, revealing that genomic imprinting is a rapidly evolving phenomenon in plants. Nevertheless, conserved PEGs exhibited signs of selection, suggesting that a subset of imprinted genes play an important functional role and are therefore maintained in plants. Like in Arabidopsis, PEGs in C. rubella are frequently associated with the presence of transposable elements that preferentially belong to helitron and MuDR families. Our data further reveal that MEGs and PEGs differ in their targeting by 24-nucleotide small RNAs and asymmetric DNA methylation, suggesting different mechanisms establishing DNA methylation at MEGs and PEGs.

Hybrid seed lethality is a widespread type of reproductive barrier among angiosperm taxa that contributes to species divergence by preventing gene flow between natural populations. Besides its ecological importance, it is an important obstacle to plant breeding strategies . Hybrid seed lethality is mostly due to a failure of the nourishing endosperm tissue, resulting in embryo arrest. The cause of this failure is a parental dosage imbalance in the endosperm that can be a consequence of either differences in parental ploidy levels or differences in the 'effective ploidy', also known as the endosperm balance number (EBN). Hybrid seed defects exhibit a parent-of-origin pattern, suggesting that differences in number or expression strength of parent-of-origin-specific imprinted genes underpin, as the primary or the secondary cause, the molecular basis of the EBN. Here, we have tested this concept in the genus Capsella and show that the effective ploidy of three Capsella species correlates with the number and expression level of paternally expressed genes (PEGs). Importantly, the number of PEGs and the effective ploidy decrease with the selfing history of a species: the obligate outbreeder Capsella grandiflora had the highest effective ploidy, followed by the recent selfer Capsella rubella and the ancient selfer Capsella orientalis. PEGs were associated with the presence of transposable elements and their silencing mark, DNA methylation in CHH context (where H denotes any base except C). This suggests that transposable elements have driven the imprintome divergence between Capsella species. Together, we propose that variation in transposable element insertions, the resulting differences in PEG number and divergence in their expression level form one component of the effective ploidy variation between species of different breeding system histories, and, as a consequence, allow the establishment of endosperm-based hybridization barriers.

Plant mating systems have profound effects on levels and structuring of genetic variation and can affect the impact of natural selection. Although theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population-genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles, whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, mating system and demography shape the impact of purifying selection on genomic variation in These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.

Plant mating systems have profound effects on levels and structuring of genetic variation, and can affect the impact of natural selection. While theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer Arabis alpina. We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian A. alpina show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, both mating system and demography shape the impact of purifying selection on genomic variation in A. alpina. These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.SignificanceIntermediate outcrossing rates are theoretically predicted to maintain effective selection against harmful alleles, but few studies have empirically tested this prediction using genomic data. We used whole-genome resequencing data from alpine rock-cress to study how genetic variation and purifying selection vary with mating system. We find that populations with intermediate outcrossing rates have similar levels of genetic diversity as outcrossing populations, and that purifying selection against harmful alleles is efficient in mixed-mating populations. In contrast, self-fertilizing populations from Scandinavia have strongly reduced genetic diversity, and accumulate harmful mutations, likely as a result of demographic effects of postglacial colonization. Our results suggest that mixed-mating populations can avoid the negative evolutionary consequences of high self-fertilization rates.

Summary A crucial step in the transition from outcrossing to self‐fertilization is the loss of genetic self‐incompatibility (SI). In the Brassicaceae, SI involves the interaction of female and male specificity components, encoded by the genes SRK and SCR at the self‐incompatibility locus (S‐locus). Theory predicts that S‐linked mutations, and especially dominant mutations in SCR, are likely to contribute to loss of SI. However, few studies have investigated the contribution of dominant mutations to loss of SI in wild plant species. Here, we investigate the genetic basis of loss of SI in the self‐fertilizing crucifer species Capsella orientalis, by combining genetic mapping, long‐read sequencing of complete S‐haplotypes, gene expression analyses and controlled crosses. We show that loss of SI in C. orientalis occurred < 2.6 Mya and maps as a dominant trait to the S‐locus. We identify a fixed frameshift deletion in the male specificity gene SCR and confirm loss of male SI specificity. We further identify an S‐linked small RNA that is predicted to cause dominance of self‐compatibility. Our results agree with predictions on the contribution of dominant S‐linked mutations to loss of SI, and thus provide new insights into the molecular basis of mating system transitions.