Target Enrichment and Extensive Population Sampling Help Untangle the Recent, Rapid Radiation of Oenothera Sect. Calylophus.

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking noncoding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic, and morphological diversity in the group. We found that the combination of exons and flanking noncoding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation. [Gypsum endemism; Oenothera sect. Calylophus; Onagraceae; phylogenomics; pollinator shift; recent radiation; target enrichment.].

The evolution of multi-gene families and metabolic pathways in the evening primroses (Oenothera: Onagraceae): A comparative transcriptomics approach.

The plant genus Oenothera has played an important role in the study of plant evolution of genomes and plant defense and reproduction. Here, we build on the 1kp transcriptomic dataset by creating 44 new transcriptomes and analyzing a total of 63 transcriptomes to present a large-scale comparative study across 29 Oenothera species. Our dataset included 30.4 million reads per individual and 2.3 million transcripts on average. We used this transcriptome resource to examine genome-wide evolutionary patterns and functional diversification by searching for orthologous genes and performed gene family evolution analysis. We found wide heterogeneity in gene family evolution across the genus, with section Oenothera exhibiting the most pronounced evolutionary changes. Overall, more significant gene family expansions occurred than contractions. We also analyzed the molecular evolution of phenolic metabolism by retrieving proteins annotated for phenolic enzymatic complexes. We identified 1,568 phenolic genes arranged into 83 multigene families that varied widely across the genus. All taxa experienced rapid phenolic evolution (fast rate of genomic turnover) involving 33 gene families, which exhibited large expansions, gaining about 2-fold more genes than they lost. Upstream enzymes phenylalanine ammonia-lyase (PAL) and 4-coumaroyl: CoA ligase (4CL) accounted for most of the significant expansions and contractions. Our results suggest that adaptive and neutral evolutionary processes have contributed to Oenothera diversification and rapid gene family evolution.

Correction of frameshift mutations in the atpB gene by translational recoding in chloroplasts of Oenothera and tobacco.

Translational recoding, also known as ribosomal frameshifting, is a process that causes ribosome slippage along the messenger RNA, thereby changing the amino acid sequence of the synthesized protein. Whether the chloroplast employs recoding is unknown. I-iota, a plastome mutant of Oenothera (evening primrose), carries a single adenine insertion in an oligoA stretch [11A] of the atpB coding region (encoding the ß-subunit of the ATP synthase). The mutation is expected to cause synthesis of a truncated, nonfunctional protein. We report that a full-length AtpB protein is detectable in I-iota leaves, suggesting operation of a recoding mechanism. To characterize the phenomenon, we generated transplastomic tobacco lines in which the atpB reading frame was altered by insertions or deletions in the oligoA motif. We observed that insertion of two adenines was more efficiently corrected than insertion of a single adenine, or deletion of one or two adenines. We further show that homopolymeric composition of the oligoA stretch is essential for recoding, as an additional replacement of AAA lysine codon by AAG resulted in an albino phenotype. Our work provides evidence for the operation of translational recoding in chloroplasts. Recoding enables correction of frameshift mutations and can restore photoautotrophic growth in the presence of a mutation that otherwise would be lethal.

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.

Why So Variable: Can Genetic Variance in Flowering Thresholds Be Maintained by Fluctuating Selection?

We use integral projection models (IPMs) and individual-based simulations to study the evolution of genetic variance in two monocarpic plant systems. Previous approaches combining IPMs with an adaptive dynamics-style invasion analysis predicted that genetic variability in the size threshold for flowering will not be maintained, which conflicts with empirical evidence. We ask whether this discrepancy can be resolved by making more realistic assumptions about the underlying genetic architecture, assuming a multilocus quantitative trait in an outcrossing diploid species. To do this, we embed the infinitesimal model of quantitative genetics into an IPM for a size-structured cosexual plant species. The resulting IPM describes the joint dynamics of individual size and breeding value of the evolving trait. We apply this general framework to the monocarpic perennials Oenothera glazioviana and Carlina vulgaris. The evolution of heritable variation in threshold size is explored in both individual-based models (IBMs) and IPMs, using a mutation rate modifier approach. In the Oenothera model, where the environment is constant, there is selection against producing genetically variable offspring. In the Carlina model, where the environment varies between years, genetically variable offspring provide a selective advantage, allowing the maintenance of genetic variability. The contrasting predictions of adaptive dynamics and quantitative genetics models for the same system suggest that fluctuating selection may be more effective at maintaining genetic variation than previously thought.

Proteomic analysis on roots of Oenothera glazioviana under copper-stress conditions.

Proteomic studies were performed to identify proteins involved in the response of Oenothera glazioviana seedlings under Cu stress. Exposure of 28-d-old seedlings to 50 µM CuSO4 for 3 d led to inhibition of shoot and root growth as well as a considerable increase in the level of lipid peroxidation in the roots. Cu absorbed by O. glazioviana accumulated more easily in the root than in the shoot. Label-free proteomic analysis indicated 58 differentially abundant proteins (DAPs) of the total 3,149 proteins in the roots of O. glazioviana seedlings, of which 36 were upregulated and 22 were downregulated under Cu stress conditions. Gene Ontology analysis showed that most of the identified proteins could be annotated to signal transduction, detoxification, stress defence, carbohydrate, energy, and protein metabolism, development, and oxidoreduction. We also retrieved 13 proteins from the enriched Kyoto Encyclopaedia of Genes and Genomes and the protein-protein interaction databases related to various pathways, including the citric acid (CA) cycle. Application of exogenous CA to O. glazioviana seedlings exposed to Cu alleviated the stress symptoms. Overall, this study provided new insights into the molecular mechanisms of plant response to Cu at the protein level in relation to soil properties.

Spontaneous Chloroplast Mutants Mostly Occur by Replication Slippage and Show a Biased Pattern in the Plastome of Oenothera.

Spontaneous plastome mutants have been used as a research tool since the beginning of genetics. However, technical restrictions have severely limited their contributions to research in physiology and molecular biology. Here, we used full plastome sequencing to systematically characterize a collection of 51 spontaneous chloroplast mutants in Oenothera (evening primrose). Most mutants carry only a single mutation. Unexpectedly, the vast majority of mutations do not represent single nucleotide polymorphisms but are insertions/deletions originating from DNA replication slippage events. Only very few mutations appear to be caused by imprecise double-strand break repair, nucleotide misincorporation during replication, or incorrect nucleotide excision repair following oxidative damage. U-turn inversions were not detected. Replication slippage is induced at repetitive sequences that can be very small and tend to have high A/T content. Interestingly, the mutations are not distributed randomly in the genome. The underrepresentation of mutations caused by faulty double-strand break repair might explain the high structural conservation of seed plant plastomes throughout evolution. In addition to providing a fully characterized mutant collection for future research on plastid genetics, gene expression, and photosynthesis, our work identified the spectrum of spontaneous mutations in plastids and reveals that this spectrum is very different from that in the nucleus.

Geographic variation in floral allometry suggests repeated transitions between selfing and outcrossing in a mixed mating plant.

PREMISE OF THE STUDY: Isometric and allometric scaling of a conserved floral plan could provide a parsimonious mechanism for rapid and reversible transitions between breeding systems. This scaling may occur during transitions between predominant autogamy and xenogamy, contributing to the maintenance of a stable mixed mating system. METHODS: We compared nine disjunct populations of the polytypic, mixed mating species Oenothera flava (Onagraceae) to two parapatric relatives, the obligately xenogamous species O. acutissima and the mixed mating species O. triloba. We compared floral morphology of all taxa using principal component analysis (PCA) and developmental trajectories of floral organs using ANCOVA homogeneity of slopes. KEY RESULTS: The PCA revealed both isometric and allometric scaling of a conserved floral plan. Three principal components (PCs) explained 92.5% of the variation in the three species. PC1 predominantly loaded on measures of floral size and accounts for 36% of the variation. PC2 accounted for 35% of the variation, predominantly in traits that influence pollinator handling. PC3 accounted for 22% of the variation, primarily in anther-stigma distance (herkogamy). During O. flava subsp. taraxacoides development, style elongation was accelerated relative to anthers, resulting in positive herkogamy. During O. flava subsp. flava development, style elongation was decelerated, resulting in zero or negative herkogamy. Of the two populations with intermediate morphology, style elongation was accelerated in one population and decelerated in the other. CONCLUSIONS: Isometric and allometric scaling of floral organs in North American Oenothera section Lavauxia drive variation in breeding system. Multiple developmental paths to intermediate phenotypes support the likelihood of multiple mating system transitions.

No evidence that sex and transposable elements drive genome size variation in evening primroses.

Genome size varies dramatically across species, but despite an abundance of attention there is little agreement on the relative contributions of selective and neutral processes in governing this variation. The rate of sex can potentially play an important role in genome size evolution because of its effect on the efficacy of selection and transmission of transposable elements (TEs). Here, we used a phylogenetic comparative approach and whole genome sequencing to investigate the contribution of sex and TE content to genome size variation in the evening primrose (Oenothera) genus. We determined genome size using flow cytometry for 30 species that vary in genetic system and find that variation in sexual/asexual reproduction cannot explain the almost twofold variation in genome size. Moreover, using whole genome sequences of three species of varying genome sizes and reproductive system, we found that genome size was not associated with TE abundance; instead the larger genomes had a higher abundance of simple sequence repeats. Although it has long been clear that sexual reproduction may affect various aspects of genome evolution in general and TE evolution in particular, it does not appear to have played a major role in genome size evolution in the evening primroses.

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.

Local topography shapes fine-scale spatial genetic structure in the Arkansas Valley evening primrose, Oenothera harringtonii (Onagraceae).

Identifying factors that shape the spatial distribution of genetic variation is crucial to understanding many population- and landscape-level processes. In this study, we explore fine-scale spatial genetic structure in Oenothera harringtonii (Onagraceae), an insect-pollinated, gravity-dispersed herb endemic to the grasslands of south-central and southeastern Colorado, USA. We genotyped 315 individuals with 11 microsatellite markers and utilized a combination of spatial autocorrelation analyses and landscape genetic models to relate life history traits and landscape features to dispersal processes. Spatial genetic structure was consistent with theoretical expectations of isolation by distance, but this pattern was weak (Sp = 0.00374). Anisotropic analyses indicated that spatial genetic structure was markedly directional, in this case consistent with increased dispersal along prominent slopes. Landscape genetic models subsequently confirmed that spatial genetic variation was significantly influenced by local topographic heterogeneity, specifically that geographic distance, elevation and aspect were important predictors of spatial genetic structure. Among these variables, geographic distance was ~68% more important than elevation in describing spatial genetic variation, and elevation was ~42% more important than aspect after removing the effect of geographic distance. From these results, we infer a mechanism of hydrochorous seed dispersal along major drainages aided by seasonal monsoon rains. Our findings suggest that landscape features may shape microevolutionary processes at much finer spatial scales than typically considered, and stress the importance of considering how particular dispersal vectors are influenced by their environmental context.

Macroevolution of plant defenses against herbivores in the evening primroses.

Plant species vary greatly in defenses against herbivores, but existing theory has struggled to explain this variation. Here, we test how phylogenetic relatedness, tradeoffs, trait syndromes, and sexual reproduction affect the macroevolution of defense. To examine the macroevolution of defenses, we studied 26 Oenothera (Onagraceae) species, combining chemistry, comparative phylogenetics and experimental assays of resistance against generalist and specialist herbivores. We detected dozens of phenolic metabolites within leaves, including ellagitannins (ETs), flavonoids, and caffeic acid derivatives (CAs). The concentration and composition of phenolics exhibited low to moderate phylogenetic signal. There were clear negative correlations between multiple traits, supporting the prediction of allocation tradeoffs. There were also positively covarying suites of traits, but these suites did not strongly predict resistance to herbivores and thus did not act as defensive syndromes. By contrast, specific metabolites did correlate with the performance of generalist and specialist herbivores. Finally, that repeated losses of sex in Oenothera was associated with the evolution of increased flavonoid diversity and altered phenolic composition. These results show that secondary chemistry has evolved rapidly during the diversification of Oenothera. This evolution has been marked by allocation tradeoffs between traits, some of which are related to herbivore performance. The repeated loss of sex appears also to have constrained the evolution of plant secondary chemistry, which may help to explain variation in defense among plants.

Microsatellite primers in Oenothera harringtonii (Onagraceae), an annual endemic to the shortgrass prairie of Colorado.

PREMISE OF THE STUDY: Microsatellite markers were developed in the annual herb, Oenothera harringtonii, to investigate patterns of genetic diversity, gene flow, and parentage within and among populations of this Colorado endemic. METHODS AND RESULTS: Ten polymorphic loci were identified in O. harringtonii and tested in four populations sampled across the range of the species. These loci contained trinucleotide repeats with 7-29 alleles per locus. Nine of the 10 loci also amplified in O. caespitosa subsp. macroglottis, O. caespitosa subsp. marginata, and O. caespitosa subsp. navajoensis. In addition, we optimized three markers developed for O. biennis and provide reports of their effectiveness in all four taxa. CONCLUSIONS: These results indicate the utility of these markers in O. harringtonii for future studies of genetic structure, gene flow, and parentage as well as their applicability in other members of the O. caespitosa species complex.

Social dilemmas among supergenes: intragenomic sexual conflict and a selfing solution in Oenothera.

Recombination is a powerful policing mechanism to control intragenomic cheats. The "parliament of the genes" can often rapidly block driving genes from cheating during meiosis. But what if the genome parliament is reduced to only two members, or supergenes? Using a series of simple game-theoretic models inspired by the peculiar genetics of Oenothera sp., we illustrate that a two supergene genome (α and ß) can produce a number of surprising evolutionary dynamics, including increases in lineage longevity following a transition from sexuality (outcrossing) to asexuality (clonal self-fertilization). We end by interpreting the model in the broader context of the evolution of mutualism, which highlights that greater α, ß cooperation in the self-fertilizing model can be viewed as an example of partner fidelity driving multilineage cooperation.

Loss of sexual recombination and segregation is associated with increased diversification in evening primroses.

The loss of sexual recombination and segregation in asexual organisms has been portrayed as an irreversible process that commits asexually reproducing lineages to reduced diversification. We test this hypothesis by estimating rates of speciation, extinction, and transition between sexuality and functional asexuality in the evening primroses. Specifically, we estimate these rates using the recently developed BiSSE (Binary State Speciation and Extinction) phylogenetic comparative method, which employs maximum likelihood and Bayesian techniques. We infer that net diversification rates (speciation minus extinction) in functionally asexual evening primrose lineages are roughly eight times faster than diversification rates in sexual lineages, largely due to higher speciation rates in asexual lineages. We further reject the hypothesis that a loss of recombination and segregation is irreversible because the transition rate from functional asexuality to sexuality is significantly greater than zero and in fact exceeded the reverse rate. These results provide the first empirical evidence in support of the alternative theoretical prediction that asexual populations should instead diversify more rapidly than sexual populations because they are free from the homogenizing effects of sexual recombination and segregation. Although asexual reproduction may often constrain adaptive evolution, our results show that the loss of recombination and segregation need not be an evolutionary dead end in terms of diversification of lineages.

Uncoupling of sexual reproduction from homologous recombination in homozygous Oenothera species.

Salient features of the first meiotic division are independent segregation of chromosomes and homologous recombination (HR). In non-sexually reproducing, homozygous species studied to date HR is absent. In this study, we constructed the first linkage maps of homozygous, bivalent-forming Oenothera species and provide evidence that HR was exclusively confined to the chromosome ends of all linkage groups in our population. Co-segregation of complementary DNA-based markers with the major group of AFLP markers indicates that HR has only a minor role in generating genetic diversity of this taxon despite its efficient adaptation capability. Uneven chromosome condensation during meiosis in Oenothera may account for restriction of HR. The use of plants with ancient chromosomal arm arrangement demonstrates that limitation of HR occurred before and independent from species hybridizations and reciprocal translocations of chromosome arms-a phenomenon, which is widespread in the genus. We propose that consecutive loss of HR favored the evolution of reciprocal translocations, beneficial superlinkage groups and ultimately permanent translocation heterozygosity.

Climate, niche evolution, and diversification of the "bird-cage" evening primroses (Oenothera, sections Anogra and Kleinia).

We integrate climatic niche models and dated phylogenies to characterize the evolution of climatic niches in Oenothera sections Anogra and Kleinia (Onagraceae), and from that we make inferences on diversification in relation to climate. The evolution of climatic tolerances in Anogra + Kleinia has been heterogeneous, across phylogenetic groups and across different dimensions of climate. All the extant taxa occur in semiarid to arid conditions (annual precipitation of 10.1-49.1 cm and high temperatures in the warmest month of 28.5 degrees-40.1 degrees C), but there is striking variation among taxa in their climatic tolerances, especially temperature (minimum temperatures in the coldest month of -14.0 degrees to 5.3 degrees C) and summer versus winter precipitation (precipitation in the warmest quarter of 0.6-19.4 cm). Climatic disparity is especially pronounced in two subclades (californica, deltoides) that radiated in the southwestern United States and California, apparently including both divergent and convergent evolution of climatic tolerances. This niche evolution is remarkable, given the probable timescale of the radiation (approximately 1 million years). We suggest that the spatiotemporal climatic heterogeneity of western North America has served as a driver of diversification. Our data are also consistent with Axelrod's hypothesis that the spread of arid conditions in western North America stimulated diversification of arid-adapted lineages.

Meiotic events in Oenothera - a non-standard pattern of chromosome behaviour.

The genus Oenothera shows an intriguing extent of permanent translocation heterozygosity. Reciprocal translocations of chromosome arms in species or populations result in various kinds of chromosome multivalents in diakinesis. Early meiotic events conditioning such chromosome behaviour are poorly understood. We found a surprising uniformity of the leptotene-diplotene period, regardless of the chromosome configuration at diakinesis (ring of 14, 7 bivalents, mixture of bivalents and multivalents). It appears that the earliest chromosome interactions at Oenothera meiosis are untypical, since they involve pericentromeric regions. During early leptotene, proximal chromosome parts cluster and form a highly polarized Rabl configuration. Telomeres associated in pairs were seen at zygotene. The high degree of polarization of meiotic nuclei continues for an exceptionally long period, i.e., during zygotene-pachytene into the diplotene contraction stage. The Rabl-polarized meiotic architecture and clustering of pericentromeres suggest a high complexity of karyotypes, not only in structural heterozygotes but also in bivalent-forming homozygous species.

Molecular marker systems for Oenothera genetics.

The genus Oenothera has an outstanding scientific tradition. It has been a model for studying aspects of chromosome evolution and speciation, including the impact of plastid nuclear co-evolution. A large collection of strains analyzed during a century of experimental work and unique genetic possibilities allow the exchange of genetically definable plastids, individual or multiple chromosomes, and/or entire haploid genomes (Renner complexes) between species. However, molecular genetic approaches for the genus are largely lacking. In this study, we describe the development of efficient PCR-based marker systems for both the nuclear genome and the plastome. They allow distinguishing individual chromosomes, Renner complexes, plastomes, and subplastomes. We demonstrate their application by monitoring interspecific exchanges of genomes, chromosome pairs, and/or plastids during crossing programs, e.g., to produce plastome-genome incompatible hybrids. Using an appropriate partial permanent translocation heterozygous hybrid, linkage group 7 of the molecular map could be assigned to chromosome 9.8 of the classical Oenothera map. Finally, we provide the first direct molecular evidence that homologous recombination and free segregation of chromosomes in permanent translocation heterozygous strains is suppressed.

The complete nucleotide sequences of the 5 genetically distinct plastid genomes of Oenothera, subsection Oenothera: II. A microevolutionary view using bioinformatics and formal genetic data.

A unique combination of genetic features and a rich stock of information make the flowering plant genus Oenothera an appealing model to explore the molecular basis of speciation processes including nucleus-organelle coevolution. From representative species, we have recently reported complete nucleotide sequences of the 5 basic and genetically distinguishable plastid chromosomes of subsection Oenothera (I-V). In nature, Oenothera plastid genomes are associated with 6 distinct, either homozygous or heterozygous, diploid nuclear genotypes of the 3 basic genomes A, B, or C. Artificially produced plastome-genome combinations that do not occur naturally often display interspecific plastome-genome incompatibility (PGI). In this study, we compare formal genetic data available from all 30 plastome-genome combinations with sequence differences between the plastomes to uncover potential determinants for interspecific PGI. Consistent with an active role in speciation, a remarkable number of genes have high Ka/Ks ratios. Different from the Solanacean cybrid model Atropa/tobacco, RNA editing seems not to be relevant for PGIs in Oenothera. However, predominantly sequence polymorphisms in intergenic segments are proposed as possible sources for PGI. A single locus, the bidirectional promoter region between psbB and clpP, is suggested to contribute to compartmental PGI in the interspecific AB hybrid containing plastome I (AB-I), consistent with its perturbed photosystem II activity.