The evolution of selfing in Arabidopsis thaliana

Fixation of Expression Divergences by Natural Selection in Arabidopsis Coding Genes

Functional divergences of coding genes can be caused by divergences in their coding sequences and expression. However, whether and how expression divergences and coding sequence divergences coevolve is not clear. Gene expression divergences in differentiated cells and tissues recapitulate developmental models within a species, while gene expression divergences between analogous cells and tissues resemble traditional phylogenies in different species, suggesting that gene expression divergences are molecular traits that can be used for evolutionary studies. Using transcriptomes and evolutionary proxies to study gene expression divergences among differentiated cells and tissues in Arabidopsis, expression divergences of coding genes are shown to be strongly anti-correlated with phylostrata (gene ages), indicators of selective constraint Ka/Ks (nonsynonymous replacement rate/synonymous substitution rate) and indicators of positive selection (frequency of loci with Ka/Ks > 1), but only weakly or not correlated with indicators of neutral selection (Ks). Our results thus suggest that expression divergences largely coevolve with coding sequence divergences, suggesting that expression divergences of coding genes are selectively fixed by natural selection but not neutral selection, which provides a molecular framework for trait diversification, functional adaptation and speciation. Our findings therefore support that positive selection rather than negative or neutral selection is a major driver for the origin and evolution of Arabidopsis genes, supporting the Darwinian theory at molecular levels.

Evaluating the Roles of Drift and Selection in Trait Loss along an Elevational Gradient

Traits that have lost function sometimes persist through evolutionary time. These traits may persist if there is not enough standing genetic variation for the trait to allow a response to selection, if selection against the trait is weak relative to drift, or if the trait has a residual function. To determine the evolutionary processes shaping whether nonfunctional traits are retained or lost, we investigated short stamens in 16 populations of Arabidopsis thaliana along an elevational cline in northeast Spain. We found a cline in short stamen number from retention of short stamens in high elevation populations to incomplete loss in low elevation populations. We did not find evidence that limited genetic variation constrains the loss of short stamens at high elevations, nor evidence for divergent selection on short stamens between high and low elevations. Finally, we identified loci associated with short stamens in northeast Spain that are different from loci associated with variation in short stamen number across latitudes from a previous study. Overall, we did not identify the evolutionary mechanisms contributing to an elevational cline in short stamen number but did identify different genetic loci underlying variation in short stamen along similar phenotypic clines.

Fine mapping of a major co-localized QTL associated with self-incompatibility identified in two F2 populations (broccoli × cauliflower and cauliflower × Chinese kale)

KEY MESSAGE

A major QTL responsible for self-incompatibility was stably identified in two F2 populations. Through fine mapping and qRT-PCR analysis, ARK3 emerged as the most promising candidate gene, playing a pivotal role in regulating self-incompatibility in Brassica oleracea. Self-incompatibility (SI) is a common phenomenon in Brassica oleracea species, which can maintain genetic diversity but will also limit seed production. Although the S locus has been extensively studied in Arabidopsis and some Brassicaceae crops, map-based cloning of self-incompatibility genes has not been conducted in Brassica oleracea, such as cauliflower and broccoli. In the present study, we identified a major co-localized QTL on chromosome C6 that control SI in two F2 populations derived from intervarietal crosses: broccoli × cauliflower (CL_F2) and cauliflower × Chinese kale (CJ_F2). Subsequently, this QTL was narrowed down to 168.5 Kb through fine mapping using 3,429 F2:3 progenies and 12 available KASP markers. Within this 168.5 Kb region, BolC6t39084H, a homologue of Arabidopsis ARK3, could be a candidate gene that plays a key role in regulating SI in B. oleracea species. This finding can pave the way for an in-depth understanding of the molecular mechanisms underlying SI, and will contribute to the seed production of B. oleracea vegetables.

Analysis of randomly mutated AlSRKb genes reveals that most loss-of-function mutations cause defects in plasma membrane localization

Only very limited information is available on why some nonsynonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on SRK which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 300 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed using error-prone PCR and the genotype and self-incompatibility phenotype of each transformant were determined. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. Our findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of AlSRKb function. In addition, we examined whether the RandomForest and Extreme Gradient Boosting methods could predict the self-incompatibility phenotypes of SRK mutants.

A pollen selection system links self and interspecific incompatibility in the Brassicaceae

Self-incompatibility and recurrent transitions to self-compatibility have shaped the extant mating systems underlying the nonrandom mating critical for speciation in angiosperms. Linkage between self-incompatibility and speciation is illustrated by the shared pollen rejection pathway between self-incompatibility and interspecific unilateral incompatibility (UI) in the Brassicaceae. However, the pollen discrimination system that activates this shared pathway for heterospecific pollen rejection remains unknown. Here we show that Stigma UI3.1, the genetically identified stigma determinant of UI in Arabidopsis lyrata × Arabidopsis arenosa crosses, encodes the S-locus-related glycoprotein 1 (SLR1). Heterologous expression of A. lyrata or Capsella grandiflora SLR1 confers on some Arabidopsis thaliana accessions the ability to discriminate against heterospecific pollen. Acquisition of this ability also requires a functional S-locus receptor kinase (SRK), whose ligand-induced dimerization activates the self-pollen rejection pathway in the stigma. SLR1 interacts with SRK and interferes with SRK homomer formation. We propose a pollen discrimination system based on competition between basal or ligand-induced SLR1-SRK and SRK-SRK complex formation. The resulting SRK homomer levels would be sensed by the common pollen rejection pathway, allowing discrimination among conspecific self- and cross-pollen as well as heterospecific pollen. Our results establish a mechanistic link at the pollen recognition phase between self-incompatibility and interspecific incompatibility.

Self-incompatibility phenotypes of SRK mutants can be predicted with high accuracy

Only very limited information is available on why some non-synonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on S-locus receptor kinase, SRK, which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 299 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed and analyzed to determine the genotype and self-incompatibility phenotype of each transformant. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. These mutations were likely to underlie the self-incompatibility defect as they caused significant changes to amino acid properties. Such findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of function. In addition, this study showed the RandomForest and Extreme Gradient Boosting methods could predict self-incompatibility phenotypes of SRK mutants with high accuracy.

Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. SI is controlled by a single S-locus with multiple haplotypes (S-haplotypes). When the pistil and pollen share the same S-haplotype, the pollen is recognized as self and rejected by the pistil. This review introduces our research on Brassicaceae and Solanaceae SI systems to identify the S-determinants encoded at the S-locus and uncover the mechanisms of self/nonself-discrimination and pollen rejection. The recognition mechanisms of SI systems differ between these families. A self-recognition system is adopted by Brassicaceae, whereas a collaborative nonself-recognition system is used by Solanaceae. Work by our group and subsequent studies indicate that plants have evolved diverse SI systems.

Balancing growth amidst salt stress - lifestyle perspectives from the extremophyte model Schrenkiella parvula

Schrenkiella parvula, a leading extremophyte model in Brassicaceae, can grow and complete its lifecycle under multiple environmental stresses, including high salinity. Yet, the key physiological and structural traits underlying its stress-adapted lifestyle are unknown along with trade-offs when surviving salt stress at the expense of growth and reproduction. We aimed to identify the influential adaptive trait responses that lead to stress-resilient and uncompromised growth across developmental stages when treated with salt at levels known to inhibit growth in Arabidopsis and most crops. Its resilient growth was promoted by traits that synergistically allowed primary root growth in seedlings, the expansion of xylem vessels across the root-shoot continuum, and a high capacity to maintain tissue water levels by developing thicker succulent leaves while enabling photosynthesis during salt stress. A successful transition from vegetative to reproductive phase was initiated by salt-induced early flowering, resulting in viable seeds. Self-fertilization in salt-induced early flowering was dependent upon filament elongation in flowers otherwise aborted in the absence of salt during comparable plant ages. The maintenance of leaf water status promoting growth, and early flowering to ensure reproductive success in a changing environment, were among the most influential traits that contributed to the extremophytic lifestyle of S. parvula.

Genomic insights into inter- and intraspecific mating system shifts in Primulina

The mating system shift from outcrossing to selfing is one of the most frequent evolutionary trends in flowering plants. However, the genomic consequences of this shift remain poorly understood. Specifically, the relative importance of the demographic and genetic processes causing changes in genetic variation and selection efficacy associated with the evolution of selfing is unclear. Here we sequenced the genomes of two Primulina species with contrasting mating systems: P. eburnea (outcrossing) versus P. tabacum (outcrossing, mixed-mating and selfing populations). Whole-genome resequencing data were used to investigate the genomic consequences of mating system shifts within and between species. We found that highly selfing populations of P. tabacum display loss of genetic diversity, increased deleterious mutations, higher genomic burden and fewer adaptive substitutions. However, compared with outcrossing populations, mixed-mating populations did not display loss of genetic diversity and accumulation of genetic load. We find no evidence of population bottlenecks associated with the shift to selfing, which suggests that the genetic effects of selfing on N_e and possibly linked selection, rather than demographic history, are the primary drivers of diversity reduction in highly selfing populations. Our results highlight the importance of distinguishing the relative contribution of mating system and demography on the genomic consequences associated with mating system evolution in plants.

The selfing syndrome and beyond: diverse evolutionary consequences of mating system transitions in plants

The shift from outcrossing to self-fertilization (selfing) is considered one of the most prevalent evolutionary transitions in flowering plants. Selfing species tend to share similar reproductive traits in morphology and function, and such a set of traits is called the 'selfing syndrome'. Although the genetic basis of the selfing syndrome has been of great interest to evolutionary biologists, knowledge of the causative genes or mutations was limited until recently. Thanks to advances in population genomic methodologies combined with high-throughput sequencing technologies, several studies have successfully unravelled the molecular and genetic basis for evolution of the selfing syndrome in Capsella, Arabidopsis, Solanum and other genera. Here we first introduce recent research examples that have explored the loci, genes and mutations responsible for the selfing syndrome traits, such as reductions in petal size or in pollen production, that are mainly relevant to pre-pollination processes. Second, we review the relationship between the evolution of selfing and interspecific pollen transfer, highlighting the findings of post-pollination reproductive barriers at the molecular level. We then discuss the emerging view of patterns in evolution of the selfing syndrome, such as the pervasive involvement of loss-of-function mutations and the relative importance of selection versus neutral degradation. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Generation of Arabidopsis thaliana transformants showing the self-recognition activity of Brassica rapa

Self-incompatibility in the Brassicaceae family is governed by SRK and SCR, which are two highly polymorphic genes located at the S-locus. Previously, the Arabidopsis lyrata SRK and SCR genes were introduced into Arabidopsis thaliana to generate self-incompatible lines. However, there are no reports showing that Brassica SRK and SCR genes confer self-incompatibility in A. thaliana. Doing so would further advance the mechanistic understanding of self-incompatibility in Brassicaceae. Therefore, we attempted to generate A. thaliana transformants showing the self-recognition activity of Brassica rapa by introducing BrSCR along with a chimeric BrSRK (BrSRK chimera, in which the kinase domain of BrSRK was replaced with that of AlSKR-b). We found that the BrSRK chimera and BrSCR of B. rapa S-9 and S-46 haplotypes, but not those of S-29, S-44, and S-60 haplotypes, conferred self-recognition activity in A. thaliana. Analyses of A. thaliana transformants expressing mutant variants of the BrSRK-9 chimera and BrSCR-9 revealed that mutations at the amino acid residues involved in BrSRK9-BrSCR9 interaction caused defects in the self-incompatibility response. The method developed in this study for generating self-incompatible A. thaliana transformants showing B. rapa self-recognition activity will be useful for analysis of self-recognition mechanisms in Brassicaceae.

Cell-cell signaling during the Brassicaceae self-incompatibility response

Self-incompatibility (SI) is a mechanism that many plant families employ to prevent self-fertilization. In the Brassicaceae, the S-haplotype-specific interaction of the pollen-borne ligand, and a stigma-specific receptor protein kinase triggers a signaling cascade that culminates in the rejection of self-pollen. While the upstream molecular components at the receptor level of the signaling pathway have been extensively studied, the intracellular responses beyond receptor activation were not as well understood. Recent research has uncovered several key molecules and signaling events that operate in concert for the manifestation of the self-incompatible responses in Brassicaceae stigmas. Here, we review the recent discoveries in both the compatible and self-incompatible pathways and provide new perspectives on the early stages of Brassicaceae pollen-pistil interactions.

Origin, loss, and regain of self-incompatibility in angiosperms

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

The rate and molecular spectrum of mutation are selectively maintained in yeast

What determines the rate (μ) and molecular spectrum of mutation is a fundamental question. The prevailing hypothesis asserts that natural selection against deleterious mutations has pushed μ to the minimum achievable in the presence of genetic drift, or the drift barrier. Here we show that, contrasting this hypothesis, μ substantially exceeds the drift barrier in diverse organisms. Random mutation accumulation (MA) in yeast frequently reduces μ, and deleting the newly discovered mutator gene PSP2 nearly halves μ. These results, along with a comparison between the MA and natural yeast strains, demonstrate that μ is maintained above the drift barrier by stabilizing selection. Similar comparisons show that the mutation spectrum such as the universal AT mutational bias is not intrinsic but has been selectively preserved. These findings blur the separation of mutation from selection as distinct evolutionary forces but open the door to alleviating mutagenesis in various organisms by genome editing.

How natural selection shapes the rate and molecular spectrum of mutations is debated. Yeast mutation accumulation experiments identify a gene promoting mutagenesis and show stabilizing selection maintaining the mutation rate above the drift barrier. Selection also preserves the mutation spectrum.

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.

Self-compatibility in peach [Prunus persica (L.) Batsch]: patterns of diversity surrounding the S-locus and analysis of SFB alleles

Self-incompatibility (SI) to self-compatibility (SC) transition is one of the most frequent and prevalent evolutionary shifts in flowering plants. Prunus L. (Rosaceae) is a genus of over 200 species most of which exhibit a Gametophytic SI system. Peach [Prunus persica (L.) Batsch; 2n = 16] is one of the few exceptions in the genus known to be a fully self-compatible species. However, the evolutionary process of the complete and irreversible loss of SI in peach is not well understood and, in order to fill that gap, in this study 24 peach accessions were analyzed. Pollen tube growth was controlled in self-pollinated flowers to verify their self-compatible phenotypes. The linkage disequilibrium association between alleles at the S-locus and linked markers at the end of the sixth linkage group was not significant (P > 0.05), except with the closest markers suggesting the absence of a signature of negative frequency dependent selection at the S-locus. Analysis of SFB1 and SFB2 protein sequences allowed identifying the absence of some variable and hypervariable domains and the presence of additional α-helices at the C-termini. Molecular and evolutionary analysis of SFB nucleotide sequences showed a signature of purifying selection in SFB2, while the SFB1 seemed to evolve neutrally. Thus, our results show that the SFB2 allele diversified after P. persica and P. dulcis (almond) divergence, a period which is characterized by an important bottleneck, while SFB1 diversified at a transition time between the bottleneck and population expansion.

Double-Locking Mechanism of Self-Compatibility in Arabidopsis thaliana: The Synergistic Effect of Transcriptional Depression and Disruption of Coding Region in the Male Specificity Gene

Self-compatibility in Arabidopsis thaliana represents the relatively recent disruption of ancestral obligate cross pollination, recognized as one of the prevalent evolutionary pathways in flowering plants, as noted by Darwin. Our previous study found that inversion of the male specificity gene (SP11/SCR) disrupted self-incompatibility, which was restored by overexpressing the SCR with the reversed inversion. However, SCR in A. thaliana has other mutations aside from the pivotal inversion, in both promoter and coding regions, with probable effects on transcriptional regulation. To examine the functional consequences of these mutations, we conducted reciprocal introductions of native promoters and downstream sequences from orthologous loci of self-compatible A. thaliana and self-incompatible A. halleri. Use of this inter-species pair enabled us to expand the scope of the analysis to transcriptional regulation and deletion in the intron, in addition to inversion in the native genomic background. Initial analysis revealed that A. thaliana has a significantly lower basal expression level of SCR transcripts in the critical reproductive stage compared to that of A. halleri, suggesting that the promoter was attenuated in inducing transcription in A. thaliana. However, in reciprocal transgenic experiments, this A. thaliana promoter was able to restore partial function if coupled with the functional A. halleri coding sequence, despite extensive alterations due to the self-compatible mode of reproduction in A. thaliana. This represents a synergistic effect of the promoter and the inversion resulting in fixation of self-compatibility, primarily enforced by disruption of SCR. Our findings elucidate the functional and evolutionary context of the historical transition in A. thaliana thus contributing to the understanding of the molecular events leading to development of self-compatibility.

Biased Gene Conversion Constrains Adaptation in Arabidopsis thaliana

Reduction of fitness due to deleterious mutations imposes a limit to adaptive evolution. By characterizing features that influence this genetic load we may better understand constraints on responses to both natural and human-mediated selection. Here, using whole-genome, transcriptome, and methylome data from >600 Arabidopsis thaliana individuals, we set out to identify important features influencing selective constraint. Our analyses reveal that multiple factors underlie the accumulation of maladaptive mutations, including gene expression level, gene network connectivity, and gene-body methylation. We then focus on a feature with major effect, nucleotide composition. The ancestral vs. derived status of segregating alleles suggests that GC-biased gene conversion, a recombination-associated process that increases the frequency of G and C nucleotides regardless of their fitness effects, shapes sequence patterns in A. thaliana Through estimation of mutational effects, we present evidence that biased gene conversion hinders the purging of deleterious mutations and contributes to a genome-wide signal of decreased efficacy of selection. By comparing these results to two outcrossing relatives, Arabidopsis lyrata and Capsella grandiflora, we find that protein evolution in A. thaliana is as strongly affected by biased gene conversion as in the outcrossing species. Last, we perform simulations to show that natural levels of outcrossing in A. thaliana are sufficient to facilitate biased gene conversion despite increased homozygosity due to selfing. Together, our results show that even predominantly selfing taxa are susceptible to biased gene conversion, suggesting that it may constitute an important constraint to adaptation among plant species.

Adaptive reduction of male gamete number in the selfing plant Arabidopsis thaliana

The number of male gametes is critical for reproductive success and varies between and within species. The evolutionary reduction of the number of pollen grains encompassing the male gametes is widespread in selfing plants. Here, we employ genome-wide association study (GWAS) to identify underlying loci and to assess the molecular signatures of selection on pollen number-associated loci in the predominantly selfing plant Arabidopsis thaliana. Regions of strong association with pollen number are enriched for signatures of selection, indicating polygenic selection. We isolate the gene REDUCED POLLEN NUMBER1 (RDP1) at the locus with the strongest association. We validate its effect using a quantitative complementation test with CRISPR/Cas9-generated null mutants in nonstandard wild accessions. In contrast to pleiotropic null mutants, only pollen numbers are significantly affected by natural allelic variants. These data support theoretical predictions that reduced investment in male gametes is advantageous in predominantly selfing species.

Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data

Several methods based on the Sequential Markovian coalescence (SMC) have been developed that make use of genome sequence data to uncover population demographic history, which is of interest in its own right and is a key requirement to generate a null model for selection tests. While these methods can be applied to all possible kind of species, the underlying assumptions are sexual reproduction in each generation and non-overlapping generations. However, in many plants, invertebrates, fungi and other taxa, those assumptions are often violated due to different ecological and life history traits, such as self-fertilization or long term dormant structures (seed or egg-banking). We develop a novel SMC-based method to infer 1) the rates/parameters of dormancy and of self-fertilization, and 2) the populations' past demographic history. Using simulated data sets, we demonstrate the accuracy of our method for a wide range of demographic scenarios and for sequence lengths from one to 30 Mb using four sampled genomes. Finally, we apply our method to a Swedish and a German population of Arabidopsis thaliana demonstrating a selfing rate of ca. 0.87 and the absence of any detectable seed-bank. In contrast, we show that the water flea Daphnia pulex exhibits a long lived egg-bank of three to 18 generations. In conclusion, we here present a novel method to infer accurate demographies and life-history traits for species with selfing and/or seed/egg-banks. Finally, we provide recommendations for the use of SMC-based methods for non-model organisms, highlighting the importance of the per site and the effective ratios of recombination over mutation.

Deleterious Mutation Accumulation in Arabidopsis thaliana Pollen Genes: A Role for a Recent Relaxation of Selection

In many studies, sex-related genes have been found to evolve rapidly. We therefore expect plant pollen genes to evolve faster than sporophytic genes. In addition, pollen genes are expressed as haploids which can itself facilitate rapid evolution because recessive advantageous and deleterious alleles are not masked by dominant alleles. However, this mechanism is less straightforward to apply in the model plant species Arabidopsis thaliana. For 1 Myr, A. thaliana has been self-compatible, a life history switch that has caused: a reduction in pollen competition, increased homozygosity, and a dilution of masking in diploid expressed, sporophytic genes. In this study, we have investigated the relative strength of selection on pollen genes compared with sporophytic genes in A. thaliana. We present two major findings: 1) before becoming self-compatible, positive selection was stronger on pollen genes than sporophytic genes for A. thaliana and 2) current polymorphism data indicate that selection is weaker on pollen genes compared with sporophytic genes. This weaker selection on pollen genes can in part be explained by their higher tissue specificity, which in outbreeding plants can be outweighed by the effects of haploid expression and pollen competition. These results indicate that since A. thaliana has become self-compatible, selection on pollen genes has become more relaxed. This has led to higher polymorphism levels and a higher build-up of deleterious mutations in pollen genes compared with sporophytic genes.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

Generation of Transgenic Self-Incompatible Arabidopsis thaliana Shows a Genus-Specific Preference for Self-Incompatibility Genes

Brassicaceae species employ both self-compatibility and self-incompatibility systems to regulate post-pollination events. Arabidopsis halleri is strictly self-incompatible, while the closely related Arabidopsis thaliana has transitioned to self-compatibility with the loss of functional S-locus genes during evolution. The downstream signaling protein, ARC1, is also required for the self-incompatibility response in some Arabidopsis and Brassica species, and its gene is deleted in the A. thaliana genome. In this study, we attempted to reconstitute the SCR-SRK-ARC1 signaling pathway to restore self-incompatibility in A. thaliana using genes from A. halleri and B. napus, respectively. Several of the transgenic A. thaliana lines expressing the A. halleriSCR₁₃-SRK₁₃-ARC1 transgenes displayed self-incompatibility, while all the transgenic A. thaliana lines expressing the B. napusSCR₁-SRK₁-ARC1 transgenes failed to show any self-pollen rejection. Furthermore, our results showed that the intensity of the self-incompatibility response in transgenic A. thaliana plants was not associated with the expression levels of the transgenes. Thus, this suggests that there are differences between the Arabidopsis and Brassica self-incompatibility signaling pathways, which perhaps points to the existence of other factors downstream of B. napusSRK that are absent in Arabidopsis species.

Family quarrels in seeds and rapid adaptive evolution in Arabidopsis

Evolutionary conflict can drive rapid adaptive evolution, sometimes called an arms race, because each party needs to respond continually to the adaptations of the other. Evidence for such arms races can sometimes be seen in morphology, in behavior, or in the genes underlying sexual interactions of host-pathogen interactions, but is rarely predicted a priori. Kin selection theory predicts that conflicts of interest should usually be reduced but not eliminated among genetic relatives, but there is little evidence as to whether conflict within families can drive rapid adaptation. Here we test multiple predictions about how conflict over the amount of resources an offspring receives from its parent would drive rapid molecular evolution in seed tissues of the flowering plant Arabidopsis As predicted, there is more adaptive evolution in genes expressed in Arabidopsis seeds than in other specialized organs, more in endosperms and maternal tissues than in embryos, and more in the specific subtissues involved in nutrient transfer. In the absence of credible alternative hypotheses, these results suggest that kin selection and conflict are important in plants, that the conflict includes not just the mother and offspring but also the triploid endosperm, and that, despite the conflict-reducing role of kinship, family members can engage in slow but steady tortoise-like arms races.

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.

Transposable elements drive rapid phenotypic variation in Capsella rubella

The mechanisms underlying rapid adaptation to changing environments in species with reduced genetic variation, referred to as the “genetic paradox of invasion,” are unknown. We report that transposable elements (TEs) are highly enriched in the gene promoter regions of Capsella rubella compared with its outcrossing sister species Capsella grandiflora. We also show that a number of polymorphic TEs in C. rubella are associated with changes in gene expression. Frequent TE insertions at FLOWERING LOCUS C of C. rubella affect flowering-time variation, an important life history trait correlated with fitness. These results indicate that TE insertions drive rapid phenotypic variation, which could potentially help adapting to novel environments in species with limited genetic variation.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3′ UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Genome Sequence of Jaltomata Addresses Rapid Reproductive Trait Evolution and Enhances Comparative Genomics in the Hyper-Diverse Solanaceae

Within the economically important plant family Solanaceae, Jaltomata is a rapidly evolving genus that has extensive diversity in flower size and shape, as well as fruit and nectar color, among its ∼80 species. Here, we report the whole-genome sequencing, assembly, and annotation, of one representative species (Jaltomata sinuosa) from this genus. Combining PacBio long reads (25×) and Illumina short reads (148×) achieved an assembly of ∼1.45 Gb, spanning ∼96% of the estimated genome. Ninety-six percent of curated single-copy orthologs in plants were detected in the assembly, supporting a high level of completeness of the genome. Similar to other Solanaceous species, repetitive elements made up a large fraction (∼80%) of the genome, with the most recently active element, Gypsy, expanding across the genome in the last 1–2 Myr. Computational gene prediction, in conjunction with a merged transcriptome data set from 11 tissues, identified 34,725 protein-coding genes. Comparative phylogenetic analyses with six other sequenced Solanaceae species determined that Jaltomata is most likely sister to Solanum, although a large fraction of gene trees supported a conflicting bipartition consistent with substantial introgression between Jaltomata and Capsicum after these species split. We also identified gene family dynamics specific to Jaltomata, including expansion of gene families potentially involved in novel reproductive trait development, and loss of gene families that accompanied the loss of self-incompatibility. This high-quality genome will facilitate studies of phenotypic diversification in this rapidly radiating group and provide a new point of comparison for broader analyses of genomic evolution across the Solanaceae.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Atypical Myrosinase as a Mediator of Glucosinolate Functions in Plants

Glucosinolates (GLSs) are a well-known class of specialized plant metabolites, distributed mostly in the order Brassicales. A vast research field in basic and applied sciences has grown up around GLSs owing to their presence in important agricultural crops and the model plant Arabidopsis thaliana, and their broad range of bioactivities beneficial to human health. The major purpose of GLSs in plants has been considered their function as a chemical defense against predators. GLSs are physically separated from a specialized class of beta-thioglucosidases called myrosinases, at the tissue level or at the single-cell level. They are brought together as a consequence of tissue damage, primarily triggered by herbivores, and their interaction results in the release of toxic volatile chemicals including isothiocyanates. In addition, recent studies have suggested that plants may adopt other strategies independent of tissue disruption for initiating GLS breakdown to cope with certain biotic/abiotic stresses. This hypothesis has been further supported by the discovery of an atypical class of GLS-hydrolyzing enzymes possessing features that are distinct from those of the classical myrosinases. Nevertheless, there is only little information on the physiological importance of atypical myrosinases. In this review, we focus on the broad diversity of the beta-glucosidase subclasses containing known atypical myrosinases in A. thaliana to discuss the hypothesis that numerous members of these subclasses can hydrolyze GLSs to regulate their diverse functions in plants. Also, the increasingly broadening functional repertoires of known atypical/classical myrosinases are described with reference to recent findings. Assessment of independent insights gained from A. thaliana with respect to (1) the phenotype of mutants lacking genes in the GLS metabolic/breakdown pathways, (2) fluctuation in GLS contents/metabolism under specific conditions, and (3) the response of plants to exogenous GLSs or their hydrolytic products, will enable us to reconsider the physiological importance of GLS breakdown in particular situations, which is likely to be regulated by specific beta-glucosidases.

Higher Rates of Protein Evolution in the Self-Fertilizing Plant Arabidopsis thaliana than in the Out-Crossers Arabidopsis lyrata and Arabidopsis halleri

The common transition from out-crossing to self-fertilization in plants decreases effective population size. This is expected to result in a reduced efficacy of natural selection and in increased rates of protein evolution in selfing plants compared with their outcrossing congeners. Prior analyses, based on a very limited number of genes, detected no differences between the rates of protein evolution in the selfing Arabidopsis thaliana compared with the out-crosser Arabidopsis lyrata. Here, we reevaluate this trend using the complete genomes of A. thaliana, A. lyrata, Arabidopsis halleri, and the outgroups Capsella rubella and Thellungiella parvula. Our analyses indicate slightly but measurably higher nonsynonymous divergences (dN), synonymous divergences (dS) and dN/dS ratios in A. thaliana compared with the other Arabidopsis species, indicating that purifying selection is indeed less efficacious in A. thaliana.

Archaic lineages broaden our view on the history of Arabidopsis thaliana

Contents Summary 1194 I. Introduction 1194 II. Origin of the A. thaliana species 1194 III. The classic model of the history of A. thaliana 1195 IV. New genomic data from outside Eurasia challenge our view of A. thaliana history 1195 V. Conclusions 1197 Acknowledgements 1197 References 1197 SUMMARY: Natural variation in Arabidopsis thaliana has contributed to discoveries in diverse areas of plant biology. While A. thaliana has typically been considered a weed associated primarily with human-mediated environments, including agricultural and urban sites and railways, it has recently been shown that it is also native in remote natural areas, including high altitude sites in Eurasia and Africa, from the Atlas mountains in Morocco to the afro-alpine regions in Eastern and South Africa to Yunnan in China, the Himalayas and the Tibetan Plateau. This finding suggests that while A. thaliana has been extensively studied in Europe and Western Asia there are still many open questions about its population history, genotype-phenotype relationships and mechanisms of adaptation.

The Molecular and Cellular Regulation of Brassicaceae Self-Incompatibility and Self-Pollen Rejection

In flowering plants, sexual reproduction is actively regulated by cell-cell communication between the male pollen and female pistil, and many species possess self-incompatibility systems for the selective rejection of self-pollen to maintain genetic diversity. The Brassicaceae self-incompatibility pathway acts early on when pollen grains have landed on the stigmatic papillae at the top of the pistil. Extensive studies have revealed that self-pollen rejection in the Brassicaceae is initiated by an S-haplotype-specific interaction between two polymorphic proteins: the pollen S-locus protein 11/S cysteine-rich (SP11/SCR) ligand and the stigma S receptor kinase (SRK). While the different S-haplotypes are typically codominant, there are several examples of dominant-recessive interactions, and a small RNA-based regulation of SP11/SCR expression has been uncovered as a mechanism behind these genetic interactions. Recent research has also added to our understanding of various cellular components in the pathway leading from the SP11/SCR-SRK interaction, including two signaling proteins, the M-locus protein kinase (MLPK) and the ARM-repeat containing 1 (ARC1) E3 ligase, as well as calcium fluxes and induction of autophagy in the stigmatic papillae. Finally, a better understanding of the compatible pollen responses that are targeted by the self-incompatibility pathway is starting to emerge, and this will allow us to more fully understand how the Brassicaceae self-incompatibility pathway causes self-pollen rejection. Here, we provide an overview of the field, highlighting recent contributions to our understanding of Brassicaceae self-incompatibility, and draw comparisons to a recently discovered unilateral incompatibility system.

Polyploid Arabidopsis species originated around recent glaciation maxima

Polyploidy may provide adaptive advantages and is considered to be important for evolution and speciation. Polyploidy events are found throughout the evolutionary history of plants, however they do not seem to be uniformly distributed along the time axis. For example, many of the detected ancient whole-genome duplications (WGDs) seem to cluster around the K/Pg boundary (∼66Mya), which corresponds to a drastic climate change event and a mass extinction. Here, we discuss more recent polyploidy events using Arabidopsis as the most developed plant model at the level of the entire genus. We review the history of the origin of allotetraploid species A. suecica and A. kamchatica, and tetraploid lineages of A. lyrata, A. arenosa and A. thaliana, and discuss potential adaptive advantages. Also, we highlight an association between recent glacial maxima and estimated times of origins of polyploidy in Arabidopsis. Such association might further support a link between polyploidy and environmental challenge, which has been observed now for different time-scales and for both ancient and recent polyploids.

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2-3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster- and species-specific marker analyses revealed that, apart from four early-generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1-3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.

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.

Long-term balancing selection contributes to adaptation in Arabidopsis and its relatives

BACKGROUND

In contrast to positive selection, which reduces genetic variation by fixing beneficial alleles, balancing selection maintains genetic variation within a population or species and plays crucial roles in adaptation in diverse organisms. However, which genes, genome-wide, are under balancing selection and the extent to which these genes are involved in adaptation are largely unknown.

RESULTS

We performed a genome-wide scan for genes under balancing selection across two plant species, Arabidopsis thaliana and its relative Capsella rubella, which diverged about 8 million generations ago. Among hundreds of genes with shared coding-region polymorphisms, we find evidence for long-term balancing selection in five genes: AT1G35220, AT2G16570, AT4G29360, AT5G38460, and AT5G44000. These genes are involved in the response to biotic and abiotic stress and other fundamental biochemical processes. More intriguingly, for these genes, we detected significant ecological diversification between the two haplotype groups, suggesting that balancing selection has been very important for adaptation.

CONCLUSIONS

Our results indicate that beyond the well-known S-locus genes and resistance genes, many loci are under balancing selection. These genes are mostly correlated with resistance to stress or other fundamental functions and likely play a more important role in adaptation to diverse habitats than previously thought.

Plant mating systems: self-incompatibility and evolutionary transitions to self-fertility in the mustard family

Flowering plants have evolved diverse mechanisms that promote outcrossing. The most widespread of these outbreeding devices are self-incompatibility systems, the highly selective prefertilization mating barriers that prevent self-fertilization by disrupting pollen-pistil interactions. Despite the advantages of outcrossing, loss of self-incompatibility has occurred repeatedly in many plant families. In the mustard family, the highly polymorphic receptors and ligands that mediate the recognition and inhibition of self-pollen in self-incompatibility have been characterized and the 3D structure of the receptor-ligand complex has been solved. Sequence analyses and empirical studies in self-incompatible and self-compatible species are elucidating the genetic basis of switches from the outcrossing to selfing modes of mating and beginning to provide clues to the diversification of the self recognition repertoire.

Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis Species

The generation and maintenance of functional variation in the pathogen defense system of plants is central to the constant evolutionary battle between hosts and parasites. If a species is susceptible to a given pathogen, hybridization and subsequent introgression of a resistance allele from a related species can potentially be an important source of new immunity and is therefore expected to be selected for in a process referred to as adaptive introgression. Here, we survey sequence variation in 10 resistance (R-) genes and compare them with 37 reference genes in natural populations of the two closely related and interfertile species: Arabidopsis lyrata and A. halleri. The R-genes are highly polymorphic in both species and show clear signs of trans-species polymorphisms. We show that A. lyrata and A. halleri have had a history of limited introgression for the reference genes. For the R-genes, the introgression rate has been significantly higher than for the reference genes, resulting in fewer fixed differences between species and a higher sharing of identical haplotypes. We conclude that R-genes likely cross the species boundaries at a higher rate than reference genes and therefore also that some of the increased diversity and trans-specific polymorphisms in R-genes is due to adaptive introgression.

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.

African genomes illuminate the early history and transition to selfing in Arabidopsis thaliana

Over the past 20 y, many studies have examined the history of the plant ecological and molecular model, Arabidopsis thaliana, in Europe and North America. Although these studies informed us about the recent history of the species, the early history has remained elusive. In a large-scale genomic analysis of African A. thaliana, we sequenced the genomes of 78 modern and herbarium samples from Africa and analyzed these together with over 1,000 previously sequenced Eurasian samples. In striking contrast to expectations, we find that all African individuals sampled are native to this continent, including those from sub-Saharan Africa. Moreover, we show that Africa harbors the greatest variation and represents the deepest history in the A. thaliana lineage. Our results also reveal evidence that selfing, a major defining characteristic of the species, evolved in a single geographic region, best represented today within Africa. Demographic inference supports a model in which the ancestral A. thaliana population began to split by 120-90 kya, during the last interglacial and Abbassia pluvial, and Eurasian populations subsequently separated from one another at around 40 kya. This bears striking similarities to the patterns observed for diverse species, including humans, implying a key role for climatic events during interglacial and pluvial periods in shaping the histories and current distributions of a wide range of species.

Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes

Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.

Recent selection for self-compatibility in a population of Leavenworthia alabamica

The evolution of self-compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self-incompatibility (SI) to self-fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S-locus. Despite permitting inbreeding, mutations at the S-locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self- and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S-locus, we collected sequence data from a ∼690 Kb region surrounding the S-locus, as well as from regions not linked to the S-locus. To test for recent positive selection acting at the S-locus, we examined polymorphism and the site-frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.

Identification and Evolution of Functional Alleles of the Previously Described Pollen Specific Myrosinase Pseudogene AtTGG6 in Arabidopsis thaliana

Myrosinases are β-thioglucoside glucohydrolases and serve as defense mechanisms against insect pests and pathogens by producing toxic compounds. AtTGG6 in Arabidopsis thaliana was previously reported to be a myrosinase pseudogene but specifically expressed in pollen. However, we found that AlTGG6, an ortholog to AtTGG6 in A. lyrata (an outcrossing relative of A. thaliana) was functional, suggesting that functional AtTGG6 alleles may still exist in A. thaliana. AtTGG6 alleles in 29 A. thaliana ecotypes were cloned and sequenced. Results indicate that ten alleles were functional and encoded Myr II type myrosinase of 512 amino acids, and myrosinase activity was confirmed by overexpressing AtTGG6 in Pichia pastoris. However, the 19 other ecotypes had disabled alleles with highly polymorphic frame-shift mutations and diversified sequences. Thirteen frame-shift mutation types were identified, which occurred independently many times in the evolutionary history within a few thousand years. The functional allele was expressed specifically in pollen similar to the disabled alleles but at a higher expression level, suggesting its role in defense of pollen against insect pests such as pollen beetles. However, the defense function may have become less critical after A. thaliana evolved to self-fertilization, and thus resulted in loss of function in most ecotypes.

Quantitative Genetics Identifies Cryptic Genetic Variation Involved in the Paternal Regulation of Seed Development

Embryonic development requires a correct balancing of maternal and paternal genetic information. This balance is mediated by genomic imprinting, an epigenetic mechanism that leads to parent-of-origin-dependent gene expression. The parental conflict (or kinship) theory proposes that imprinting can evolve due to a conflict between maternal and paternal alleles over resource allocation during seed development. One assumption of this theory is that paternal alleles can regulate seed growth; however, paternal effects on seed size are often very low or non-existent. We demonstrate that there is a pool of cryptic genetic variation in the paternal control of Arabidopsis thaliana seed development. Such cryptic variation can be exposed in seeds that maternally inherit a medea mutation, suggesting that MEA acts as a maternal buffer of paternal effects. Genetic mapping using recombinant inbred lines, and a novel method for the mapping of parent-of-origin effects using whole-genome sequencing of segregant bulks, indicate that there are at least six loci with small, paternal effects on seed development. Together, our analyses reveal the existence of a pool of hidden genetic variation on the paternal control of seed development that is likely shaped by parental conflict.

Evolution of Selfing: Recurrent Patterns in Molecular Adaptation

Selfing has evolved in animals, fungi, and plants, and since Darwin's pioneering study, it is considered one of the most frequent evolutionary trends in flowering plants. Generally, the evolution of selfing is characterized by a loss of self-incompatibility, the selfing syndrome, and changes in genome-wide polymorphism patterns. Recent interdisciplinary studies involving molecular functional experiments, genome-wide data, experimental evolution, and evolutionary ecology using Arabidopsis thaliana, Caenorhabditis elegans, and other species show that the evolution of selfing is not merely a degradation of outcrossing traits but a model for studying the recurrent patterns underlying adaptive molecular evolution. For example, in wild Arabidopsis relatives, self-compatibility evolved from mutations in the male specificity gene, S-LOCUS CYSTEINE-RICH PROTEIN/S-LOCUS PROTEIN 11 (SCR/SP11), rather than the female specificity gene, S-LOCUS RECEPTOR KINASE (SRK), supporting the theoretical prediction of sexual asymmetry. Prevalence of dominant self-compatible mutations is consistent with Haldane's sieve, which acts against recessive adaptive mutations. Time estimates based on genome-wide polymorphisms and self-incompatibility genes generally support the recent origin of selfing.

Population Genomics for Understanding Adaptation in Wild Plant Species

Darwin's theory of evolution by natural selection is the foundation of modern biology. However, it has proven remarkably difficult to demonstrate at the genetic, genomic, and population level exactly how wild species adapt to their natural environments. We discuss how one can use large sets of multiple genome sequences from wild populations to understand adaptation, with an emphasis on the small herbaceous plant Arabidopsis thaliana. We present motivation for such studies; summarize progress in describing whole-genome, species-wide sequence variation; and then discuss what insights have emerged from these resources, either based on sequence information alone or in combination with phenotypic data. We conclude with thoughts on opportunities with other plant species and the impact of expected progress in sequencing technology and genome engineering for studying adaptation in nature.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

BACKGROUND

The transition from outcrossing to selfing has long been portrayed as an 'evolutionary dead end' because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a 'rostellum') co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.

RESULTS

The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4-3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.

CONCLUSIONS

The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the 'rostellum') is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.