Hybridization and a loss of sex shape genome-wide diversity and the origin of species in the evening primroses (Oenothera, Onagraceae).

Hybridization is thought to promote speciation in at least two ways - by fixation of heterozygosity from diploid progenitors in allopolyploids, and by generation of transgressive phenotypes and shifting fitness optima during homoploid hybrid speciation. While recent studies support a hybrid origin for a growing number of species, the extent to which hybrid origins shape patterns of diversity in asexual species remains underexplored. Here we employed transcriptome sequencing and population genomic analysis to describe patterns of genomic variation in the 13 species belonging to Oenothera subsection Oenothera. Eight of these species are functionally asexual and arose by hybrid speciation from parents spanning a range of phylogenetic divergence. We showed that genomic divergence between parents has been retained as heterozygosity in functionally asexual species, and that genome-wide levels of heterozygosity in these asexuals largely reflects the divergence of parental haplotypes coupled with a breakdown in recombination and segregation across the genome. These results show that divergence between parental species and loss of sex in hybrids shape patterns of whole-genome diversity and the origin of asexual species.

Borrowed alleles and convergence in serpentine adaptation.

Serpentine barrens represent extreme hazards for plant colonists. These sites are characterized by high porosity leading to drought, lack of essential mineral nutrients, and phytotoxic levels of metals. Nevertheless, nature forged populations adapted to these challenges. Here, we use a population-based evolutionary genomic approach coupled with elemental profiling to assess how autotetraploid Arabidopsis arenosa adapted to a multichallenge serpentine habitat in the Austrian Alps. We first demonstrate that serpentine-adapted plants exhibit dramatically altered elemental accumulation levels in common conditions, and then resequence 24 autotetraploid individuals from three populations to perform a genome scan. We find evidence for highly localized selective sweeps that point to a polygenic, multitrait basis for serpentine adaptation. Comparing our results to a previous study of independent serpentine colonizations in the closely related diploid Arabidopsis lyrata in the United Kingdom and United States, we find the highest levels of differentiation in 11 of the same loci, providing candidate alleles for mediating convergent evolution. This overlap between independent colonizations in different species suggests that a limited number of evolutionary strategies are suited to overcome the multiple challenges of serpentine adaptation. Interestingly, we detect footprints of selection in A. arenosa in the context of substantial gene flow from nearby off-serpentine populations of A. arenosa, as well as from A. lyrata In several cases, quantitative tests of introgression indicate that some alleles exhibiting strong selective sweep signatures appear to have been introgressed from A. lyrata This finding suggests that migrant alleles may have facilitated adaptation of A. arenosa to this multihazard environment.

Genetic degeneration of old and young Y chromosomes in the flowering plant Rumex hastatulus.

Heteromorphic sex chromosomes have originated independently in many species, and a common feature of their evolution is the degeneration of the Y chromosome, characterized by a loss of gene content and function. Despite being of broad significance to our understanding of sex chromosome evolution, the genetic changes that occur during the early stages of Y-chromosome degeneration are poorly understood, especially in plants. Here, we investigate sex chromosome evolution in the dioecious plant Rumex hastatulus, in which X and Y chromosomes have evolved relatively recently and occur in two distinct systems: an ancestral XX/XY system and a derived XX/XY1Y2 system. This polymorphism provides a unique opportunity to investigate the effect of sex chromosome age on patterns of divergence and gene degeneration within a species. Despite recent suppression of recombination and low X-Y divergence in both systems, we find evidence that Y-linked genes have started to undergo gene loss, causing ∼ 28% and ∼ 8% hemizygosity of the ancestral and derived X chromosomes, respectively. Furthermore, genes remaining on Y chromosomes have accumulated more amino acid replacements, contain more unpreferred changes in codon use, and exhibit significantly reduced gene expression compared with their X-linked alleles, with the magnitude of these effects being greatest for older sex-linked genes. Our results provide evidence for reduced selection efficiency and ongoing Y-chromosome degeneration in a flowering plant, and indicate that Y degeneration can occur soon after recombination suppression between sex chromosomes.

Recurrent loss of sex is associated with accumulation of deleterious mutations in Oenothera.

Sexual reproduction is nearly universal among eukaryotes. Theory predicts that the rarity of asexual eukaryotic species is in part caused by accumulation of deleterious mutations and heightened extinction risk associated with suppressed recombination and segregation in asexual species. We tested this prediction with a large data set of 62 transcriptomes from 29 species in the plant genus Oenothera, spanning ten independent transitions between sexual and a functionally asexual genetic system called permanent translocation heterozygosity. Illumina short-read sequencing and de novo transcript assembly yielded an average of 16.4 Mb of sequence per individual. Here, we show that functionally asexual species accumulate more deleterious mutations than sexual species using both population genomic and phylogenetic analysis. At an individual level, asexual species exhibited 1.8 × higher heterozygosity than sexual species. Within species, we detected a higher proportion of nonsynonymous polymorphism relative to synonymous variation within asexual compared with sexual species, indicating reduced efficacy of purifying selection. Asexual species also exhibited a greater proportion of transcripts with premature stop codons. The increased proportion of nonsynonymous mutations was also positively correlated with divergence time between sexual and asexual species, consistent with Muller's ratchet. Between species, we detected repeated increases in the ratio of nonsynonymous to synonymous divergence in asexual species compared with sexually reproducing sister taxa, indicating increased accumulation of deleterious mutations. These results confirm that an important advantage of sex is that it facilitates selection against deleterious alleles, which might help to explain the dearth of extant asexual species.

Polyploidy: adaptation to the genomic environment.

Genomic evidence of ancestral whole genome duplication (WGD) and polyploidy is widespread among eukaryotic species, and especially among plants. WGD is thought to provide the raw material for adaptation in the form of duplicated genes, and polyploids are thought to benefit from both physiological and genetic buffering. Comparatively little attention has focused on the genomic challenge of polyploidy, however, although much evidence exists that polyploidy severely perturbs important cellular functions. Here, I review recent progress in the study of the re-establishment of stable meiosis in recently evolved polyploids, focusing on four plant species. This work has yielded an insight into the mechanisms underlying stabilization of genome transmission in polyploids, and is revealing remarkable parallels among diverse taxa. Importantly, these studies also provide a road map for investigating how polyploids respond to the challenge of WGD.

Genomic variation in Arabidopsis: tools and insights from next-generation sequencing.

The release of a reference genome for Arabidopsis thaliana in 2000 has been an enormous boon for the study of plant genetics. Less than a decade later, however, a revolution in sequencing technology had enabled rapid and inexpensive re-sequencing of whole A. thaliana genomes. Large-scale efforts to characterize natural genomic variation in A. thaliana have revealed remarkable intra-specific variation in this species, ranging from single-nucleotide differences to large structural rearrangements. The partitioning of this variation by geography and local adaptation has been described using powerful new methods and tools. Simultaneously, an ambitious research agenda has emerged to sequence 1001 A. thaliana lines from around the world, while sequencing of related species is enabling powerful evolutionary genomic analyses. In this review, I summarize recent progress in genomic analysis of natural variation in A. thaliana and its close relatives. This progress has set the stage for the emergence of Arabidopsis as a model genus for evolutionary and functional genomics.

Genetic adaptation associated with genome-doubling in autotetraploid Arabidopsis arenosa.

Genome duplication, which results in polyploidy, is disruptive to fundamental biological processes. Genome duplications occur spontaneously in a range of taxa and problems such as sterility, aneuploidy, and gene expression aberrations are common in newly formed polyploids. In mammals, genome duplication is associated with cancer and spontaneous abortion of embryos. Nevertheless, stable polyploid species occur in both plants and animals. Understanding how natural selection enabled these species to overcome early challenges can provide important insights into the mechanisms by which core cellular functions can adapt to perturbations of the genomic environment. Arabidopsis arenosa includes stable tetraploid populations and is related to well-characterized diploids A. lyrata and A. thaliana. It thus provides a rare opportunity to leverage genomic tools to investigate the genetic basis of polyploid stabilization. We sequenced the genomes of twelve A. arenosa individuals and found signatures suggestive of recent and ongoing selective sweeps throughout the genome. Many of these are at genes implicated in genome maintenance functions, including chromosome cohesion and segregation, DNA repair, homologous recombination, transcriptional regulation, and chromatin structure. Numerous encoded proteins are predicted to interact with one another. For a critical meiosis gene, ASYNAPSIS1, we identified a non-synonymous mutation that is highly differentiated by cytotype, but present as a rare variant in diploid A. arenosa, indicating selection may have acted on standing variation already present in the diploid. Several genes we identified that are implicated in sister chromatid cohesion and segregation are homologous to genes identified in a yeast mutant screen as necessary for survival of polyploid cells, and also implicated in genome instability in human diseases including cancer. This points to commonalities across kingdoms and supports the hypothesis that selection has acted on genes controlling genome integrity in A. arenosa as an adaptive response to genome doubling.

Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata.

Transposable elements (TEs) are often the primary determinant of genome size differences among eukaryotes. In plants, the proliferation of TEs is countered through epigenetic silencing mechanisms that prevent mobility. Recent studies using the model plant Arabidopsis thaliana have revealed that methylated TE insertions are often associated with reduced expression of nearby genes, and these insertions may be subject to purifying selection due to this effect. Less is known about the genome-wide patterns of epigenetic silencing of TEs in other plant species. Here, we compare the 24-nt siRNA complement from A. thaliana and a closely related congener with a two- to threefold higher TE copy number, Arabidopsis lyrata. We show that TEs--particularly siRNA-targeted TEs--are associated with reduced gene expression within both species and also with gene expression differences between orthologs. In addition, A. lyrata TEs are targeted by a lower fraction of uniquely matching siRNAs, which are associated with more effective silencing of TE expression. Our results suggest that the efficacy of RNA-directed DNA methylation silencing is lower in A. lyrata, a finding that may shed light on the causes of differential TE proliferation among species.

Local-scale patterns of genetic variability, outcrossing, and spatial structure in natural stands of Arabidopsis thaliana.

As Arabidopsis thaliana is increasingly employed in evolutionary and ecological studies, it is essential to understand patterns of natural genetic variation and the forces that shape them. Previous work focusing mostly on global and regional scales has demonstrated the importance of historical events such as long-distance migration and colonization. Far less is known about the role of contemporary factors or environmental heterogeneity in generating diversity patterns at local scales. We sampled 1,005 individuals from 77 closely spaced stands in diverse settings around Tübingen, Germany. A set of 436 SNP markers was used to characterize genome-wide patterns of relatedness and recombination. Neighboring genotypes often shared mosaic blocks of alternating marker identity and divergence. We detected recent outcrossing as well as stretches of residual heterozygosity in largely homozygous recombinants. As has been observed for several other selfing species, there was considerable heterogeneity among sites in diversity and outcrossing, with rural stands exhibiting greater diversity and heterozygosity than urban stands. Fine-scale spatial structure was evident as well. Within stands, spatial structure correlated negatively with observed heterozygosity, suggesting that the high homozygosity of natural A. thaliana may be partially attributable to nearest-neighbor mating of related individuals. The large number of markers and extensive local sampling employed here afforded unusual power to characterize local genetic patterns. Contemporary processes such as ongoing outcrossing play an important role in determining distribution of genetic diversity at this scale. Local "outcrossing hotspots" appear to reshuffle genetic information at surprising rates, while other stands contribute comparatively little. Our findings have important implications for sampling and interpreting diversity among A. thaliana accessions.

Indel-associated mutation rate varies with mating system in flowering plants.

A recently proposed mutational mechanism, indel-associated mutation (IDAM), posits that heterozygous insertions/deletions (indels) increase the point mutation rate at nearby nucleotides due to errors during meiosis. This mechanism could have especially dynamic consequences for the evolution of plant genomes, because the high degree of variation in the rate of self-fertilization among plant species causes differences in the heterozygosity of alleles, including indel alleles, segregating in plant species. In this study, we investigated the consequences of IDAM for species differing in mating system using both forward population genetic simulations and genomewide DNA resequencing data from Arabidopsis thaliana, Oryza sativa, and Oryza rufipogon. Simulations of different levels of selfing suggest that the effect of IDAM on surrounding nucleotide diversity should decrease with increasing selfing rate. Further simulations incorporating selfing rates and the time of onset of selfing suggest that the time since the switch to selfing also affects patterns of nucleotide diversity due to IDAM. Population genetic analyses of A. thaliana and Oryza DNA sequence data sets empirically confirmed our simulation results, revealing the strongest effect of IDAM in the outcrossing O. rufipogon, a weaker effect in the recently evolved selfer O. sativa, and the weakest effect in the relatively ancient selfer A. thaliana. These results support the novel idea that differences in life history, such as the level of selfing, can affect the per-individual mutation rate among species.

Meiotic adaptation to genome duplication in Arabidopsis arenosa.

Whole genome duplication (WGD) is a major factor in the evolution of multicellular eukaryotes, yet by doubling the number of homologs, WGD severely challenges reliable chromosome segregation, a process conserved across kingdoms. Despite this, numerous genome-duplicated (polyploid) species persist in nature, indicating early problems can be overcome. Little is known about which genes are involved--only one has been molecularly characterized. To gain new insights into the molecular basis of adaptation to polyploidy, we investigated genome-wide patterns of differentiation between natural diploids and tetraploids of Arabidopsis arenosa, an outcrossing relative of A. thaliana. We first show that diploids are not preadapted to polyploid meiosis. We then use a genome scanning approach to show that although polymorphism is extensively shared across ploidy levels, there is strong ploidy-specific differentiation in 39 regions spanning 44 genes. These are discrete, mostly single-gene peaks of sharply elevated differentiation. Among these peaks are eight meiosis genes whose encoded proteins coordinate a specific subset of early meiotic functions, suggesting these genes comprise a polygenic solution to WGD-associated chromosome segregation challenges. Our findings indicate that even conserved meiotic processes can be capable of nimble evolutionary shifts when required.

Epigenetic inheritance: what news for evolution?

Whether epigenetic variation is important in adaptive evolution has been contentious. Two recent studies in Arabidopsis thaliana significantly add to our understanding of genome-wide variation and stability of an epigenetic mark, and thus help pave the path for realistically incorporating epigenetics into evolutionary theory.

The Arabidopsis lyrata genome sequence and the basis of rapid genome size change.

We report the 207-Mb genome sequence of the North American Arabidopsis lyrata strain MN47 based on 8.3× dideoxy sequence coverage. We predict 32,670 genes in this outcrossing species compared to the 27,025 genes in the selfing species Arabidopsis thaliana. The much smaller 125-Mb genome of A. thaliana, which diverged from A. lyrata 10 million years ago, likely constitutes the derived state for the family. We found evidence for DNA loss from large-scale rearrangements, but most of the difference in genome size can be attributed to hundreds of thousands of small deletions, mostly in noncoding DNA and transposons. Analysis of deletions and insertions still segregating in A. thaliana indicates that the process of DNA loss is ongoing, suggesting pervasive selection for a smaller genome. The high-quality reference genome sequence for A. lyrata will be an important resource for functional, evolutionary and ecological studies in the genus Arabidopsis.

A triptych of the evolution of plant transposable elements.

Transposable elements (TEs) constitute the majority of angiosperm DNA, but the processes that govern their accumulation remain mysterious. Here we discuss the three major forces that govern the accumulation of TEs, corresponding to the three panels of a triptych. The first force, transposition, creates new copies of TEs, but is regulated by both host- and TE-specific mechanisms. The second force, deletion of TE DNA, is capable of removing vast swaths of genomic regions via recombinational processes, but we still have very little insight into how deletion varies across species and even among TE types. Finally, we focus on the often-ignored third panel of our triptych - the population processes that determine the ultimate evolutionary fate of TE insertions.

Epigenetic silencing of transposable elements: a trade-off between reduced transposition and deleterious effects on neighboring gene expression.

Transposable elements (TEs) are ubiquitous genomic parasites. The deleterious consequences of the presence and activity of TEs have fueled debate about the evolutionary forces countering their expansion. Purifying selection is thought to purge TE insertions from the genome, and TE sequences are targeted by hosts for epigenetic silencing. However, the interplay between epigenetic and evolutionary forces countering TE expansion remains unexplored. Here we analyze genomic, epigenetic, and population genetic data from Arabidopsis thaliana to yield three observations. First, gene expression is negatively correlated with the density of methylated TEs. Second, the signature of purifying selection is detectable for methylated TEs near genes but not for unmethylated TEs or for TEs far from genes. Third, TE insertions are distributed by age and methylation status, such that older, methylated TEs are farther from genes. Based on these observations, we present a model in which host silencing of TEs near genes has deleterious effects on neighboring gene expression, resulting in the preferential loss of methylated TEs from gene-rich chromosomal regions. This mechanism implies an evolutionary tradeoff in which the benefit of TE silencing imposes a fitness cost via deleterious effects on the expression of nearby genes.

Population and evolutionary dynamics of Helitron transposable elements in Arabidopsis thaliana.

Helitrons, a recently discovered superfamily of DNA transposons that capture host gene fragments, constitute up to 2% of the Arabidopsis thaliana genome. In this study, we identified 565 insertions of a family of nonautonomous Helitrons, known as Basho elements. We aligned subsets of these elements, estimated their phylogenetic relationships, and used branch lengths to yield insight into the age of each Basho insertion. The age distribution suggests that 87% of Bashos inserted within 5 Myr, subsequent to the divergence between A. thaliana and its sister species Arabidopsis lyrata. We screened 278 of these insertions for their presence or absence in a sample of 47 A. thaliana accessions. With both phylogenetic and population frequency data, we investigated the effects of gene density, recombination rate, and element length on Basho persistence. Our analyses suggested that longer Basho copies are less likely to persist in the genome, consistent with selection against the deleterious effects of ectopic recombination between Basho elements. Furthermore, we determined that 39% of Basho elements contain fragments of expressed protein-coding genes, but all of these fragments were explained by only 5 gene-capture events. Overall, the picture of A. thaliana Helitron evolution is one of rapid expansion, relatively few gene-capture events, and weak selection correlated with element length.