Recent speciation associated with the evolution of selfing in Capsella

The Role of Small RNA-Based Epigenetic Silencing for Purifying Selection on Transposable Elements in Capsella grandiflora

To avoid negative effects of transposable element (TE) proliferation, plants epigenetically silence TEs using a number of mechanisms, including RNA-directed DNA methylation. These epigenetic modifications can extend outside the boundaries of TE insertions and lead to silencing of nearby genes, resulting in a trade-off between TE silencing and interference with nearby gene regulation. Therefore, purifying selection is expected to remove silenced TE insertions near genes more efficiently and prevent their accumulation within a population. To explore how effects of TE silencing on gene regulation shapes purifying selection on TEs, we analyzed whole genome sequencing data from 166 individuals of a large population of the outcrossing species Capsella grandiflora. We found that most TEs are rare, and in chromosome arms, silenced TEs are exposed to stronger purifying selection than those that are not silenced by 24-nucleotide small RNAs, especially with increasing proximity to genes. An age-of-allele test of neutrality on a subset of TEs supports our inference of purifying selection on silenced TEs, suggesting that our results are robust to varying transposition rates. Our results provide new insights into the processes affecting the accumulation of TEs in an outcrossing species and support the view that epigenetic silencing of TEs results in a trade-off between preventing TE proliferation and interference with nearby gene regulation. We also suggest that in the centromeric and pericentromeric regions, the negative aspects of epigenetic TE silencing are missing.

Parallel Evolution of Common Allelic Variants Confers Flowering Diversity in Capsella rubella

Flowering time is an adaptive life history trait. Capsella rubella, a close relative of Arabidopsis thaliana and a young species, displays extensive variation for flowering time but low standing genetic variation due to an extreme bottleneck event, providing an excellent opportunity to understand how phenotypic diversity can occur with a limited initial gene pool. Here, we demonstrate that common allelic variation and parallel evolution at the FLC locus confer variation in flowering time in C. rubella. We show that two overlapping deletions in the 5' untranslated region (UTR) of C. rubella FLC, which are associated with local changes in chromatin conformation and histone modifications, reduce its expression levels and promote flowering. We further show that these two pervasive variants originated independently in natural C. rubella populations after speciation and spread to an intermediate frequency, suggesting a role of this parallel cis-regulatory change in adaptive evolution. Our results provide an example of how parallel mutations in the same 5' UTR region can shape phenotypic evolution in plants.

Evolution of interspecies unilateral incompatibility in the relatives of Arabidopsis thaliana

The evolutionary concurrence of intraspecies self-incompatibility (SI) and explosive angiosperm radiation in the Cretaceous have led to the hypothesis that SI was one of the predominant drivers of rapid speciation in angiosperms. Interspecies unilateral incompatibility (UI) usually occurs when pollen from a self-compatible (SC) species is rejected by the pistils of a SI species, while the reciprocal pollination is compatible (UC). Although this SI × SC type UI is most prevalent and viewed as a prezygotic isolation barrier to promote incipient speciation of angiosperms, comparative evidence to support such a role is lacking. We show that SI × SI type UI in SI species pairs is also common in the well-characterized accessions representing the four major lineages of the Arabidopsis genus and is developmentally regulated. This allowed us to reveal a strong correlation between UI strength and species divergence in these representative accessions. In addition, analyses of a SC accession and the pseudo-self-compatible (psc) spontaneous mutant of Arabidopsis lyrata indicate that UI shares, at least, common pollen rejection pathway with SI. Furthermore, genetic and genomic analyses of SI × SI type UI in A. lyrata × A. arenosa species pair showed that two major-effect quantitative trait loci are the stigma and pollen-side determinant of UI, respectively, which could be involved in heterospecies pollen discrimination. By revealing a close link between UI and SI pathway, particularly between UI and species divergence in these representative accessions, our findings establish a connection between SI and speciation. Thus, the pre-existence of SI system would have facilitated the evolution of UI and accordingly promote speciation.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

Competitive ability of Capsella species with different mating systems and ploidy levels

BACKGROUND AND AIMS

Capsella is a model genus for studying the transition from outcrossing to selfing, with or without change in ploidy levels. The genomic consequences and changes in reproductive traits (selfing syndrome) associated with these shifts have been studied in depth. However, potential ecological divergence among species of the genus has not been determined. Among ecological traits, competitive ability could be relevant for selfing evolution, as selfing has been shown to be statistically associated with reduced competitiveness in a recent meta-analysis.

METHODS

We assessed the effect of competition on three Capsella species differing in their mating system and ploidy level. We used an experimental design where fitness related traits were measured in focal individuals with and without competitors.

KEY RESULTS

The diploid selfer (C. rubella) was most sensitive to competition, whereas the tetraploid selfer (C. bursa-pastoris) performed the best, with the diploid outcrosser (C. grandiflora) being intermediate.

CONCLUSIONS

These results add to the detailed characterization of Capsella species and highlight the possible roles of ecological context and ploidy in the evolutionary trajectories of selfing species.

Chloroplastic and nuclear diversity of wild beets at a large geographical scale: Insights into the evolutionary history of the Beta section

Historical demographic processes and mating systems are believed to be major factors in the shaping of the intraspecies genetic diversity of plants. Among Caryophyllales, the Beta section of the genus Beta, within the Amaranthaceae/Chenopodiaceae alliance, is an interesting study model with species and subspecies (Beta macrocarpa, Beta patula, Beta vulgaris maritima and B.v. adanensis) differing in geographical distribution and mating system. In addition, one of the species, B. macrocarpa, mainly diploid, varies in its level of ploidy with a tetraploid cytotype described in the Canary Islands and in Portugal. In this study, we analyzed the nucleotide diversity of chloroplastic and nuclear sequences on a representative sampling of species and subspecies of the Beta section (except B. patula). Our objectives were (1) to assess their genetic relationships through phylogenetic and multivariate analyses, (2) relate their genetic diversity to their mating system, and (3) reconsider the ploidy status and the origin of the Canarian Beta macrocarpa.

Consequences of multiple mating-system shifts for population and range-wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating-system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e ) caused by selfing, small-flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large-flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage-wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating-system differentiation observed across the range of this species.

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.

The Evolutionary Dynamics of Orthologs That Shift in Gene Body Methylation between Arabidopsis Species

DNA methylation labels a specific subset of genes in plant genomes. Recent work has shown that this gene-body methylation (gbM) is a conserved feature of orthologs, because highly methylated genes in one species tend to be highly methylated in another. In this study, we examined the exceptions to that rule by identifying genes that differ in gbM status between two plant species-Arabidopsis thaliana and Arabidopsis lyrata. Using Capsella grandiflora as an outgroup, we polarized the loss and gain of gbM for orthologs in the Arabidopsis lineage. Our survey identified a few hundred genes that differed between ingroup species, out of ∼9,000 orthologs. The estimated rate of gbM gain was ∼2 × 10-9 per gene per year for both ingroup taxa and was similar to the loss rate in A. lyrata. In contrast, A. thaliana had a ∼3-fold higher estimated rate of gbM loss per gene, consistent with a recent diminishment of genome size. As in previous studies, we found that body-methylated genes were expressed broadly across tissues and were also longer than other genic sets. Genes that differed in gbM status exhibited higher variance in expression between species than genes that were body-methylated in both species. Moreover, the gain of gbM in one lineage tended to be associated with an increase of expression in that lineage. The genes that varied in gbM status between species varied more significantly in length between species than other sets of genes; we hypothesize that length is a key feature in the transition between body-methylated and nonmethylated genes.

The potential role of hybridization in diversification and speciation in an insular plant lineage: insights from synthetic interspecific hybrids

Hybridization is recognized as an important process in plant evolution, and this may be particularly true for island plants where several biotic and abiotic factors facilitate interspecific hybridization. Although rarely done, experimental studies could provide insights into the potential of natural hybridization to generate diversity when species come into contact in the dynamic island setting. The potential of hybridization to generate morphological variation was analysed within and among 12 families (inbred lines) of an F4 hybrid generation between two species of Tolpis endemic to the Canary Islands. Combinations of characters not seen in the parents were present in hybrids. Several floral and vegetative characters were transgressive relative to their parents. Morphometric studies of floral, vegetative and fruit characters revealed that several F4 families were phenotypically distinct from other families, and from their parents. The study demonstrates that morphologically distinct pollen-fertile lines, potentially worthy of taxonomic recognition if occurring in nature, can be generated in four generations. The ability of the hybrid lines to set self-seed would reduce gene flow among the lines, and among the hybrids and their parental species. Selfing would also facilitate the fixation of characters within each of the lines. Overall, the results show the considerable potential of hybridization for generating diversity and distinct phenotypes in island lineages.

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.

On the young age of intraspecific herbaceous taxa

Dated phylogenies rarely include the divergence times of sister intraspecific taxa, and when they do little is said about this subject. We show that over 90% of the intraspecific plant taxa found in a literature search are estimated to be 5 million yr old or younger, with only 4% of taxa estimated to be over 10 million yr old or older. A Bayesian analysis of intraspecific taxon ages indicates that indeed these taxa are expected to be < 10 million yr old. This result for the young age of intraspecific taxa is consistent with the earlier observation that post-pollination reproductive barriers develop between 5 and 10 million yr after lineage splitting, thus leading to species formation. If lineages have not graduated to the species level of divergence by 10 million yr or so, they are likely to have gone extinct by that time as a result of narrow geographical distributions, narrow niche breadths, and relatively small numbers across populations.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Maintenance of mixed mating after the loss of self-incompatibility in a long-lived perennial herb

BACKGROUND AND AIMS

Many hermaphroditic plants avoid self-fertilization by means of a molecular self-incompatibility (SI) system, a complex trait that is difficult to evolve but relatively easy to lose. Loss of SI is a prerequisite for an evolutionary transition from obligate outcrossing to self-fertilization, which may bring about rapid changes in the genetic diversity and structure of populations. Loss of SI is also often followed by the evolution of a 'selfing syndrome', with plants having small flowers, little nectar and few pollen grains per ovule. Here, we document the loss of SI in the long-lived Spanish toadflax Linaria cavanillesii, which has led to mixed mating rather than a transition to a high rate of selfing and in which an outcrossing syndrome has been maintained.

METHODS

We performed crosses within and among six populations of L. cavanillesii in the glasshouse, measured floral traits in a common-garden experiment, performed a pollen-limitation experiment in the field and conducted population genetic analyses using microsatellites markers.

KEY RESULTS

Controlled crosses revealed variation in SI from fully SI through intermediate SI to fully self-compatible (SC). Flowers of SC individuals showed no evidence of a selfing syndrome. Although the SC population of L. cavanillesii had lower within-population genetic diversity than SI populations, as expected, population differentiation among all populations was extreme and represents an F_ST outlier in the distribution for both selfing and outcrossing species of flowering plants.

CONCLUSIONS

Together, our results suggest that the transition to SC in L. cavanillesii has probably been very recent, and may have been aided by selection during or following a colonization bottleneck rather than in the absence of pollinators. We find little indication that the transition to SC has been driven by selection for reproductive assurance under conditions currently prevailing in natural populations.

Shifts in outcrossing rates and changes to floral traits are associated with the evolution of herbicide resistance in the common morning glory

Human-mediated selection can strongly influence the evolutionary response of natural organisms within ecological timescales. But what traits allow for, or even facilitate, adaptation to the strong selection humans impose on natural systems? Using a combination of laboratory and greenhouse studies of 32 natural populations of the common agricultural weed, Ipomoea purpurea, we show that herbicide-resistant populations self-fertilise more than susceptible populations. We likewise show that anther-stigma distance, a floral trait associated with self-fertilisation in this species, exhibits a nonlinear relationship with resistance such that the most and least resistant populations exhibit lower anther-stigma separation compared to populations with moderate levels of resistance. Overall, our results extend the general finding that plant mating can be impacted by human-mediated agents of selection to that of the extreme selection of the agricultural system. This work highlights the influence of human-mediated selection on rapid responses of natural populations that can lead to unexpected long-term evolutionary consequences.

Repeated Inactivation of the First Committed Enzyme Underlies the Loss of Benzaldehyde Emission after the Selfing Transition in Capsella

The enormous species richness of flowering plants is at least partly due to floral diversification driven by interactions between plants and their animal pollinators [1, 2]. Specific pollinator attraction relies on visual and olfactory floral cues [3-5]; floral scent can not only attract pollinators but also attract or repel herbivorous insects [6-8]. However, despite its central role for plant-animal interactions, the genetic control of floral scent production and its evolutionary modification remain incompletely understood [9-13]. Benzenoids are an important class of floral scent compounds that are generated from phenylalanine via several enzymatic pathways [14-17]. Here we address the genetic basis of the loss of floral scent associated with the transition from outbreeding to selfing in the genus Capsella. While the outbreeding C. grandiflora emits benzaldehyde as a major constituent of its floral scent, this has been lost in the selfing C. rubella. We identify the Capsella CNL1 gene encoding cinnamate:CoA ligase as responsible for this variation. Population genetic analysis indicates that CNL1 has been inactivated twice independently in C. rubella via different novel mutations to its coding sequence. Together with a recent study in Petunia [18], this identifies cinnamate:CoA ligase as an evolutionary hotspot for mutations causing the loss of benzenoid scent compounds in association with a shift in the reproductive strategy of Capsella from pollination by insects to self-fertilization.

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.

Rapid Evolution of Genomic Imprinting in Two Species of the Brassicaceae

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.

Transposable element evolution in the allotetraploid Capsella bursa-pastoris

PREMISE OF THE STUDY

Shifts in ploidy affect the evolutionary dynamics of genomes in a myriad of ways. Population genetic theory predicts that transposable element (TE) proliferation may follow because the genomewide efficacy of selection should be reduced and the increase in gene copies may mask the deleterious effects of TE insertions. Moreover, in allopolyploids, TEs may further accumulate because of hybrid breakdown of TE silencing. However, to date the evidence of TE proliferation following an increase in ploidy is mixed, and the relative importance of relaxed selection vs. silencing breakdown remains unclear.

METHODS

We used high-coverage whole-genome sequence data to evaluate the abundance, genomic distribution, and population frequencies of TEs in the self-fertilizing recent allotetraploid Capsella bursa-pastoris (Brassicaceae). We then compared the C. bursa-pastoris TE profile with that of its two parental diploid species, outcrossing C. grandiflora and self-fertilizing C. orientalis.

KEY RESULTS

We found no evidence that C. bursa-pastoris has experienced a large genomewide proliferation of TEs relative to its parental species. However, when centromeric regions are excluded, we found evidence of significantly higher abundance of retrotransposons in C. bursa-pastoris along the gene-rich chromosome arms compared with C. grandiflora and C. orientalis.

CONCLUSIONS

The lack of a genomewide effect of allopolyploidy on TE abundance, combined with the increases TE abundance in gene-rich regions, suggests that relaxed selection rather than hybrid breakdown of host silencing explains the TE accumulation in C. bursa-pastoris.

Yeast two-hybrid interactions between Arabidopsis lyrata S Receptor Kinase and the ARC1 E3 ligase

Here we describe protein-protein interactions between signaling components in the conserved self-incompatibility pathway from Brassica spp. and Arabidopsis lyrata. Previously, we had demonstrated that ARC1 is necessary in A. lyrata for the rejection of self-pollen by the self-incompatibility pathway. The results described here demonstrate that A. lyrata ARC1 interacts with A. lyrata S Receptor Kinase (SRK1) in the yeast 2-hybrid system. A. lyrata ARC1 also interacted with B. napus SRK910 illustrating that interactions in this pathway are conserved across species. Finally, we discuss how the more widely occurring interactions between SRK and ARC1-related family members may be modulated in vivo by expression and subcellular localization patterns resulting in a particular response.

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.

Population and phylogenomic decomposition via genotyping-by-sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic-scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.

R-gene variation across Arabidopsis lyrata subspecies: effects of population structure, selection and mating system

Background

Examining allelic variation of R-genes in closely related perennial species of Arabidopsis thaliana is critical to understanding how population structure and ecology interact with selection to shape the evolution of innate immunity in plants. We finely sampled natural populations of Arabidopsis lyrata from the Great Lakes region of North America (A. l. lyrata) and broadly sampled six European countries (A. l. petraea) to investigate allelic variation of two R-genes (RPM1 and WRR4) and neutral genetic markers (Restriction Associated DNA sequences and microsatellites) in relation to mating system, phylogeographic structure and subspecies divergence.

Results

Fine-scale sampling of populations revealed strong effects of mating system and population structure on patterns of polymorphism for both neutral loci and R-genes, with no strong evidence for selection. Broad geographic sampling revealed evidence of balancing selection maintaining polymorphism in R-genes, with elevated heterozygosity and diversity compared to neutral expectations and sharing of alleles among diverged subspecies. Codon-based tests detected both positive and purifying selection for both R-genes, as commonly found for animal immune genes.

Conclusions

Our results highlight that combining fine and broad-scale sampling strategies can reveal the multiple factors influencing polymorphism and divergence at potentially adaptive genes such as R-genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0665-5) contains supplementary material, which is available to authorized users.

Divergence and isolation of cryptic sympatric taxa within the annual legume Amphicarpaea bracteata

The amphicarpic annual legume Amphicarpaea bracteata is unusual in producing aerial and subterranean cleistogamous flowers that always self‐fertilize and, less commonly, aerial chasmogamous flowers that outcross. Although both morphologic and genetic variants are known in this highly selfing species, debate continues over whether this variation is continuous, reflecting the segregation of standing genetic variation, or discontinuous, reflecting distinct taxa that rarely intercross. We characterized SNP variation in 128 individuals in southern Wisconsin to assess within‐ and among‐population variation at 3928 SNPs. We also assessed genotype and leaf morphology in an additional 76 individuals to connect phenotypic variation with genetic variation. Genetic variation maps onto three strongly divergent and highly inbred genetic groups showing little relation to site location. Each group has a distinct phenotype, but the divergence of these groups differs from the varietal divisions previously identified based on morphological characters. Like previous authors, we argue that the taxonomy of this species should be revised. Despite extensive sympatry, estimates of among‐group migration rates are low, and hybrid individuals were at low frequency (<2%) in our dataset. Restricted gene flow likely results from high selfing rates and partial reproductive incompatibility as evidenced by the U‐shaped distribution of pairwise F_ST values reflecting “islands” of genomic divergence. These islands may be associated with hybrid incompatibility loci that arose in allopatry. The coexistence of lineages within sites may reflect density‐dependent attack by species‐specific strains of pathogenic fungi and/or root‐nodulating bacteria specializing on distinct genotypes.

Transcriptome sequencing of Crucihimalaya himalaica (Brassicaceae) reveals how Arabidopsis close relative adapt to the Qinghai-Tibet Plateau

The extreme environment of the Qinghai-Tibet Plateau (QTP) provides an ideal natural laboratory for studies on adaptive evolution. Few genome/transcriptome based studies have been conducted on how plants adapt to the environments of QTP compared to numerous studies on vertebrates. Crucihimalaya himalaica is a close relative of Arabidopsis with typical QTP distribution, and is hoped to be a new model system to study speciation and ecological adaptation in extreme environment. In this study, we de novo generated a transcriptome sequence of C. himalaica, with a total of 49,438 unigenes. Compared to five relatives, 10,487 orthogroups were shared by all six species, and 4,286 orthogroups contain putative single copy gene. Further analysis identified 487 extremely significantly positively selected genes (PSGs) in C. himalaica transcriptome. Theses PSGs were enriched in functions related to specific adaptation traits, such as response to radiation, DNA repair, nitrogen metabolism, and stabilization of membrane. These functions are responsible for the adaptation of C. himalaica to the high radiation, soil depletion and low temperature environments on QTP. Our findings indicate that C. himalaica has evolved complex strategies for adapting to the extreme environments on QTP and provide novel insights into genetic mechanisms of highland adaptation in plants.

Mating system divergence affects the distribution of sequence diversity within and among populations of recently diverged subspecies of Clarkia xantiana (Onagraceae)

PREMISE OF THE STUDY

The population biology of outcrossing and self-fertilizing taxa is thought to differ because of the advantage that selfers have in colonizing unoccupied sites where mates and pollinators may be limiting (Baker's Law). This reduced tendency for outcrossers to colonize new sites, along with their greater dependence on pollinators to disperse pollen, has the potential to differently influence the genetic diversity and structure of outcrossing and selfing populations.

METHODS

We conducted a comparative population genetic study of two sister outcrossing and selfing subspecies of Clarkia xantiana that have very recently diverged. We used DNA sequence variation (>40 kb from eight nuclear loci) from large samples of individuals from 14 populations to assess geographic patterns of genetic diversity and make inferences about the demographic and colonization histories of each subspecies.

KEY RESULTS

We show that sequence variation is strongly reduced across all selfing populations. The demographic history of selfing populations exhibits recent colonization bottlenecks, whereas such bottlenecks are rarely observed for the outcrosser. The greater effect of genetic drift in the selfer has resulted in strong population genetic structure, but with no pattern of isolation by distance. By contrast, the stronger effect of gene flow in the outcrosser has resulted in considerably less structure, but a significant pattern of isolation by distance.

CONCLUSIONS

Taken together, our results suggest that selfing populations are not at migration-drift equilibrium, are affected by strong episodes of genetic drift during colonization, and experience little or no subsequent gene flow from other populations after those founder events.

No association between plant mating system and geographic range overlap

PREMISE OF THE STUDY

Automatic self-fertilization may influence the geography of speciation, promote reproductive isolation between incipient species, and lead to ecological differentiation. As such, selfing taxa are predicted to co-occur more often with their closest relatives than are outcrossing taxa. Despite suggestions that this pattern may be general, the extent to which mating system influences range overlap in close relatives has not been tested formally across a diverse group of plant species pairs.

METHODS

We tested for a difference in range overlap between species pairs for which zero, one, or both species are selfers, using data from 98 sister species pairs in 20 genera across 15 flowering plant families. We also used divergence time estimates from time-calibrated phylogenies to ask how range overlap changes with divergence time and whether this effect depends on mating system.

KEY RESULTS

We found no evidence that automatic self-fertilization influenced range overlap of closely related plant species. Sister pairs with more recent divergence times had modestly greater range overlap, but this effect did not depend on mating system.

CONCLUSIONS

The absence of a strong influence of mating system on range overlap suggests that mating system plays a minor or inconsistent role compared with many other mechanisms potentially influencing the co-occurrence of close relatives.

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.

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.

Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

Divergent sorting of a balanced ancestral polymorphism underlies the establishment of gene-flow barriers in Capsella

In the Bateson–Dobzhansky–Muller model of genetic incompatibilities post-zygotic gene-flow barriers arise by fixation of novel alleles at interacting loci in separated populations. Many such incompatibilities are polymorphic in plants, implying an important role for genetic drift or balancing selection in their origin and evolution. Here we show that NPR1 and RPP5 loci cause a genetic incompatibility between the incipient species Capsella grandiflora and C. rubella, and the more distantly related C. rubella and C. orientalis. The incompatible RPP5 allele results from a mutation in C. rubella, while the incompatible NPR1 allele is frequent in the ancestral C. grandiflora. Compatible and incompatible NPR1 haplotypes are maintained by balancing selection in C. grandiflora, and were divergently sorted into the derived C. rubella and C. orientalis. Thus, by maintaining differentiated alleles at high frequencies, balancing selection on ancestral polymorphisms can facilitate establishing gene-flow barriers between derived populations through lineage sorting of the alternative alleles.

A hybrid incompatibility between Capsella plant species is due to an interaction between two immune regulators. Here, the authors show that highly divergent haplotypes result from balancing selection in the ancestral lineage and their sorting into derived lineages facilitated the evolution of the incompatibility.

Evolution of the Selfing Syndrome in Arabis alpina (Brassicaceae)

Introduction

The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.

Methods

Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.

Results and Discussion

Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.

Conclusion

We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.

Non-reciprocal Interspecies Hybridization Barriers in the Capsella Genus Are Established in the Endosperm

The transition to selfing in Capsella rubella accompanies its recent divergence from the ancestral outcrossing C. grandiflora species about 100,000 years ago. Whether the change in mating system was accompanied by the evolution of additional reproductive barriers that enforced species divergence remained unknown. Here, we show that C. rubella and C. grandiflora are reproductively separated by an endosperm-based, non-reciprocal postzygotic hybridization barrier. While hybridizations of C. rubella maternal plants with C. grandiflora pollen donors resulted in complete seed abortion caused by endosperm cellularization failure, the reciprocal hybridization resulted in the formation of small seeds with precociously cellularized endosperm. Strikingly, the transcriptomic response of both hybridizations mimicked respectively the response of paternal and maternal excess hybridizations in Arabidopsis thaliana, suggesting unbalanced genome strength causes hybridization failure in both species. These results provide strong support for the theory that crosses between plants of different mating systems will be unbalanced, with the outcrosser behaving like a plant of increased ploidy, evoking a response that resembles an interploidy-type seed failure. Seed incompatilibity of C. rubella pollinated by C. grandiflora followed the Bateson-Dobzhansky-Muller model, involving negative genetic interaction of multiple paternal C. grandiflora loci with at least one maternal C. rubella locus. Given that both species only recently diverged, our data suggest that a fast evolving mechanism underlies the post-zygotic hybridization barrier(s) separating both species.

Full-genome evolutionary histories of selfing, splitting, and selection in Caenorhabditis

The nematode Caenorhabditis briggsae is a model for comparative developmental evolution with C. elegans. Worldwide collections of C. briggsae have implicated an intriguing history of divergence among genetic groups separated by latitude, or by restricted geography, that is being exploited to dissect the genetic basis to adaptive evolution and reproductive incompatibility; yet, the genomic scope and timing of population divergence is unclear. We performed high-coverage whole-genome sequencing of 37 wild isolates of the nematode C. briggsae and applied a pairwise sequentially Markovian coalescent (PSMC) model to 703 combinations of genomic haplotypes to draw inferences about population history, the genomic scope of natural selection, and to compare with 40 wild isolates of C. elegans. We estimate that a diaspora of at least six distinct C. briggsae lineages separated from one another approximately 200,000 generations ago, including the “Temperate” and “Tropical” phylogeographic groups that dominate most samples worldwide. Moreover, an ancient population split in its history approximately 2 million generations ago, coupled with only rare gene flow among lineage groups, validates this system as a model for incipient speciation. Low versus high recombination regions of the genome give distinct signatures of population size change through time, indicative of widespread effects of selection on highly linked portions of the genome owing to extreme inbreeding by self-fertilization. Analysis of functional mutations indicates that genomic context, owing to selection that acts on long linkage blocks, is a more important driver of population variation than are the functional attributes of the individually encoded genes.

Frequent introgressions from diploid species contribute to the adaptation of the tetraploid Shepherd's purse (Capsella bursa-pastoris)

Adaptation is the most important ability for organisms to survive in diverse habitats. Animals have the option to escape from stressful environments, but plants do not. In plants, polyploids consist of about 30%-70% angiosperms and 95% ferns, of which some are important crops such as cotton and wheat. How polyploid plants adapt to various habitats has been a fundamental question remained largely unanswered. The tetraploid Shepherd's purse (Capsella bursa-pastoris) is one of the most successful plants on earth and has been distributed across the world, thus being an ideal model system for studying the adaptation of polyploids. We found that there are frequent introgressions from congeneric diploids to Shepherd's purse. Ecological niche modeling suggests that ecological differentiation is evident between the introgressed and non-introgressed C. bursa-pastoris, and the introgressions are a source of adaptation. This result links an evolutionary process to the adaptation of polyploids, and sheds light on the breeding strategy of polyploids as well. We conclude that frequent introgressions from congeneric diploids contributed to the acquisition of adequate genetic variations, thereby allowing C. bursa-pastoris to adapt to various habitats across the world. Our results highlight how a polyploid could have successfully established after it originated.

Hybrid origins and the earliest stages of diploidization in the highly successful recent polyploid Capsella bursa-pastoris

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.

The genetic architecture of traits associated with the evolution of self-pollination in Mimulus

Quantitative trait locus (QTL) mapping is a first step toward understanding the genetic basis of adaptive evolution and may also reveal reproductive incompatibilities unique to hybrids. In plants, the shift from outcrossing to self-pollination is common, providing the opportunity for comparisons of QTL architecture among parallel evolutionary transitions. We used QTL mapping in hybrids between the bee-pollinated monkeyflower Mimulus lewisii and the closely related selfer Mimulus parishii to determine the genetic basis of divergence in floral traits and flowering time associated with mating-system evolution, and to characterize hybrid anther sterility. We found a moderately polygenic and highly directional basis for floral size evolution, suggesting adaptation from standing variation or in pursuit of a moving optimum, whereas only a few major loci accounted for substantial flowering-time divergence. Cytonuclear incompatibilities caused hybrid anther sterility, confounding estimation of reproductive organ QTLs. The genetic architecture of floral traits associated with selfing in M. parishii was primarily polygenic, as in other QTL studies of this transition, but in contrast to the previously characterized oligogenic basis of a pollinator shift in close relatives. Hybrid anther sterility appeared parallel at the molecular level to previously characterized incompatibilities, but also raised new questions about cytonuclear co-evolution in plants.

Secondary extinctions of biodiversity

Extinctions beget further extinctions when species lose obligate mutualists, predators, prey, or hosts. Here, we develop a conceptual model of species and community attributes affecting secondary extinction likelihood, incorporating mechanisms that buffer organisms against partner loss. Specialized interactors, including 'cryptic specialists' with diverse but nonredundant partner assemblages, incur elevated risk. Risk is also higher for species that cannot either evolve new traits following partner loss or obtain novel partners in communities reorganizing under changing environmental conditions. Partner loss occurs alongside other anthropogenic impacts; multiple stressors can circumvent ecological buffers, enhancing secondary extinction risk. Stressors can also offset each other, reducing secondary extinction risk, a hitherto unappreciated phenomenon. This synthesis suggests improved conservation planning tactics and critical directions for research on secondary extinctions.

Repeated evolutionary changes of leaf morphology caused by mutations to a homeobox gene

Elucidating the genetic basis of morphological changes in evolution remains a major challenge in biology. Repeated independent trait changes are of particular interest because they can indicate adaptation in different lineages or genetic and developmental constraints on generating morphological variation. In animals, changes to "hot spot" genes with minimal pleiotropy and large phenotypic effects underlie many cases of repeated morphological transitions. By contrast, only few such genes have been identified from plants, limiting cross-kingdom comparisons of the principles of morphological evolution. Here, we demonstrate that the REDUCED COMPLEXITY (RCO) locus underlies more than one naturally evolved change in leaf shape in the Brassicaceae. We show that the difference in leaf margin dissection between the sister species Capsella rubella and Capsella grandiflora is caused by cis-regulatory variation in the homeobox gene RCO-A, which alters its activity in the developing lobes of the leaf. Population genetic analyses in the ancestral C. grandiflora indicate that the more-active C. rubella haplotype is derived from a now rare or lost C. grandiflora haplotype via additional mutations. In Arabidopsis thaliana, the deletion of the RCO-A and RCO-B genes has contributed to its evolutionarily derived smooth leaf margin, suggesting the RCO locus as a candidate for an evolutionary hot spot. We also find that temperature-responsive expression of RCO-A can explain the phenotypic plasticity of leaf shape to ambient temperature in Capsella, suggesting a molecular basis for the well-known negative correlation between temperature and leaf margin dissection.