Copy number variation in transcriptionally active regions of sexual and apomictic Boechera demonstrates independently derived apomictic lineages

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits-mode of genomic inheritance-influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S-5.8S-28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA-histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction-with potential functional consequences related to the loss of sexual reproduction.

Gene Function Rather than Reproductive Mode Drives the Evolution of RNA Helicases in Sexual and Apomictic Boechera

In higher plants, sexual and asexual reproductions through seeds (apomixis) have evolved as alternative strategies. Evolutionary advantages leading to coexistence of both reproductive modes are currently not well understood. It is expected that accumulation of deleterious mutations leads to a rapid elimination of apomictic lineages from populations. In this line, apomixis originated repeatedly, likely from deregulation of the sexual pathway, leading to alterations in the development of reproductive lineages (germlines) in apomicts as compared with sexual plants. This potentially involves mutations in genes controlling reproduction.

Increasing evidence suggests that RNA helicases are crucial regulators of germline development. To gain insights into the evolution of 58 members of this diverse gene family in sexual and apomictic plants, we applied target enrichment combined with next-generation sequencing to identify allelic variants from 24 accessions of the genus Boechera, comprising sexual, facultative, and obligate apomicts. Interestingly, allelic variants from apomicts did not show consistently increased mutation frequency. Either sequences were highly conserved in any accession, or allelic variants preferentially harbored mutations in evolutionary less conserved C- and N-terminal domains, or presented high mutation load independent of the reproductive mode. Only for a few genes allelic variants harboring deleterious mutations were only identified in apomicts. To test if high sequence conservation correlates with roles in fundamental cellular or developmental processes, we analyzed Arabidopsis thaliana mutant lines in VASA-LIKE (VASL), and identified pleiotropic defects during ovule and reproductive development. This indicates that also in apomicts mechanisms of selection are in place based on gene function.

Evolutionary history of the buildup and breakdown of the heterostylous syndrome in Plumbaginaceae

The evolutionary pathways leading to the heterostylous syndrome are not well understood, and models concerning the origins of distyly differ in the order in which reciprocal herkogamy and self-incompatibility evolve. We investigated the evolution and breakdown of distyly in Plumbaginaceae, a family with considerable diversity of floral traits and reproductive systems. Using Bayesian Markov chain Monte Carlo analyses and stochastic character mapping, we examined the evolutionary assembly and breakdown of the heterostylous syndrome based on a well-resolved phylogeny of 121 species of Plumbaginaceae and six outgroup taxa using five nuclear and plastid gene regions. We used the distribution of reproductive traits and reconstructed ancestral characters across phylogenies to evaluate competing models for the evolution of distyly. The most likely common ancestor of Plumbaginaceae was self-incompatible and monomorphic for sex-organ arrangement and pollen-stigma characters. Character state reconstructions indicated that reciprocal herkogamy evolved at least three times and that shifts to selfing and apomixis occurred on multiple occasions. Our results provide comparative support for the early ideas of H. G. Baker on evolutionary pathways in Plumbaginaceae, and the more recent selfing avoidance model by D. & B. Charlesworth in which distyly evolves from self-incompatible ancestors.

Differential activity of F-box genes and E3 ligases distinguishes sexual versus apomictic germline specification in Boechera

Germline specification is the first step during sexual and apomictic plant reproduction, and takes place in the nucellus of the ovule, a specialized domain of the reproductive flower tissues. In each case, a sporophytic cell is determined to form the sexual megaspore mother cell (MMC) or an apomictic initial cell (AIC). These differ in their developmental fates: while the MMC undergoes meiosis, the AIC modifies or omits meiosis to form the female gametophyte. Despite great interest in these distinct developmental processes, little is known about their gene regulatory basis. To elucidate the gene regulatory networks underlying germline specification, we conducted tissue-specific transcriptional profiling using laser-assisted microdissection and RNA sequencing to compare the transcriptomes of nucellar tissues between different sexual and apomictic Boechera accessions representing four species and two ploidy levels. This allowed us to distinguish between expression differences caused by genetic background or reproductive mode. Statistical data analysis revealed 45 genes that were significantly differentially expressed, and which potentially play a role for determination of the reproductive mode. Based on annotations, these included F-box genes and E3 ligases that most likely relate to genes previously described as regulators important for germline development. Our findings provide novel insights into the transcriptional basis of sexual and apomictic reproduction.

Apospory and Diplospory in Diploid Boechera (Brassicaceae) May Facilitate Speciation by Recombination-Driven Apomixis-to-Sex Reversals

Apomixis (asexual seed formation) in angiosperms occurs either sporophytically, through adventitious embryony, or gametophytically, where an unreduced female gametophyte (embryo sac) forms and produces an unreduced egg that develops into an embryo parthenogenetically. Multiple types of gametophytic apomixis occur, and these are differentiated based on where and when the unreduced gametophyte forms, a process referred to as apomeiosis. Apomeiotic gametophytes form directly from ameiotic megasporocytes, as in Antennaria-type diplospory, from unreduced spores derived from 1st division meiotic restitutions, as in Taraxacum-type diplospory, or from cells of the ovule wall, as in Hieracium-type apospory. Multiple types of apomeiosis occasionally occur in the same plant, which suggests that the different types occur in response to temporal and/or spatial shifts in termination of sexual processes and onset timing of apomeiosis processes. To better understand the origins and evolutionary implications of apomixis in Boechera (Brassicaceae), we determined apomeiosis type for 64 accessions representing 44 taxonomic units. Plants expressing apospory and diplospory were equally common, and these generally produced reduced and unreduced pollen, respectively. Apospory and diplospory occurred simultaneously in individual plants of seven taxa. In Boechera, apomixis perpetuates otherwise sterile or semisterile interspecific hybrids (allodiploids) through multiple generations. Accordingly, ample time, in these multigenerational clones, is available for rare meioses to produce haploid, intergenomically recombined male and female gametes. The fusion of such gametes could then produce segmentally autoploidized progeny. If sex re-emerges among such progeny, then new and genomically unique sexual species could evolve. Herein, we present evidence that such apomixis-facilitated speciation is occurring in Boechera, and we hypothesize that it might also be occurring in facultatively apomictic allodiploids of other angiospermous taxa.

The Boechera Genus as a Resource for Apomixis Research

The genera Boechera (A. Löve et D. Löve) and Arabidopsis, the latter containing the model plant Arabidopsis thaliana, belong to the same clade within the Brassicaceae family. Boechera is the only among the more than 370 genera in the Brassicaceae where apomixis is well documented. Apomixis refers to the asexual reproduction through seed, and a better understanding of the underlying mechanisms has great potential for applications in agriculture. The Boechera genus currently includes 110 species (of which 38 are reported to be triploid and thus apomictic), which are distributed mostly in the North America. The apomictic lineages of Boechera occur at both the diploid and triploid level and show signs of a hybridogenic origin, resulting in a modification of their chromosome structure, as reflected by alloploidy, aneuploidy, substitutions of homeologous chromosomes, and the presence of aberrant chromosomes. In this review, we discuss the advantages of the Boechera genus to study apomixis, consider its modes of reproduction as well as the inheritance and possible mechanisms controlling apomixis. We also consider population genetic aspects and a possible role of hybridization at the origin of apomixis in Boechera. The molecular tools available to study Boechera, such as transformation techniques, laser capture microdissection, analysis of transcriptomes etc. are also discussed. We survey available genome assemblies of Boechera spp. and point out the challenges to assemble the highly heterozygous genomes of apomictic species. Due to these challenges, we argue for the application of an alternative reference-free method for the comparative analysis of such genomes, provide an overview of genomic sequencing data in the genus Boechera suitable for such analysis, and provide examples of its application.

Ecological differentiation facilitates fine-scale coexistence of sexual and asexual Boechera

PREMISE OF THE STUDY

Ecological differentiation (ED) between sexual and asexual organisms may permit the maintenance of reproductive polymorphism. Several studies of sexual/asexual ED in plants have shown that the geographic ranges of asexuals extend beyond those of sexuals, often in areas of higher latitude or elevation. But very little is known about ED at fine scales, wherein coexistence of sexuals and asexuals may be permitted by differential niche occupation.

METHODS

We used 149 populations of sexual and apomictic lineages in the genus Boechera (rock cress) collected across a portion of this mustard's vast range. We characterized reproductive mode, ploidy, and species identity or hybrid parentage of each individual, and then used a multipronged statistical approach to (1) identify ED between sexuals and asexuals; (2) investigate the impacts of two confounding factors, polyploidy and hybridization, on ED; and (3) determine the environmental variables underlying ED.

KEY RESULTS

We found that sexuals and asexuals are significantly ecologically differentiated across the landscape, despite fine-scale interdigitation of these two reproductive forms. Asexual reproduction was strongly associated with greater disturbance, reduced slope, and greater environmental variability. Although ploidy had little effect on the patterns observed, hybridization has a unique impact on the relationships between asexual reproduction and specific environmental variables.

CONCLUSIONS

Ecological differentiation along the axes of disturbance, slope, and climatic variability, as well as the effects of heterozygosity, may contribute to the maintenance of sexuality and asexuality across the landscape, ultimately impacting the establishment and spread of asexual lineages.

A reexamination of the North American Crepis agamic complex and comparison with the findings of Babcock and Stebbins' classic biosystematic monograph

PREMISE OF THE STUDY

Babcock and Stebbins coined the term agamic complex in their 1938 monograph of the North American Crepis agamic complex. Despite the historical role that this complex holds in the evolutionary literature, it has not been reexamined in over 75 years. We present a thorough reevaluation of the complex to test hypotheses proposed by Babcock and Stebbins about its origins and spread, the relationships of diploids, and the nature and origins of polyploids.

METHODS

We used flow cytometry to infer ploidy of roughly 600 samples spanning the morphological and taxonomic diversity of the complex and a phylogenetic analysis of plastid DNA variation to infer maternal relationships among diploids and to infer maternal origins of polyploids.

KEY RESULTS

We identified populations of all seven recognized diploids plus one new lineage. Phylogenetic analysis of plastid DNA variation in diploids revealed a well-resolved, but moderately supported phylogeny, with evidence for monophyly of the North America Crepis agamic complex and no evidence of widespread homoploid hybridization. Polyploids showed evidence of multiple origins and a pattern of frequent local co-occurrence consistent with repeated colonization of suitable sites.

CONCLUSIONS

Our findings agree broadly with the distribution and variation of ploidy within and among species described by Babcock and Stebbins. One key difference is finding support for monophyly of North American species, and refuting their hypothesis of polyphyly. Our results provide an explicit phylogenetic framework for further study of this classic agamic complex.

Laser-assisted Microdissection (LAM) as a Tool for Transcriptional Profiling of Individual Cell Types

The understanding of developmental processes at the molecular level requires insights into transcriptional regulation, and thus the transcriptome, at the level of individual cell types. While the methods described here are generally applicable to a wide range of species and cell types, our research focuses on plant reproduction. Plant cultivation and seed production is of crucial importance for human and animal nutrition. A detailed understanding of the regulatory networks that govern the formation of the reproductive lineage (germline) and ultimately of seeds is a precondition for the targeted manipulation of plant reproduction. In particular, the engineering of apomixis (asexual reproduction through seeds) into crop plants promises great improvements, as it leads to the formation of clonal seeds that are genetically identical to the mother plant. Consequently, the cell types of the female germline are of major importance for the understanding and engineering of apomixis. However, as the corresponding cells are deeply embedded within the floral tissues, they are very difficult to access for experimental analyses, including cell-type specific transcriptomics. To overcome this limitation, sections of individual cells can be isolated by laser-assisted microdissection (LAM). While LAM in combination with transcriptional profiling allows the identification of genes and pathways active in any cell type with high specificity, establishing a suitable protocol can be challenging. Specifically, the quality of RNA obtained after LAM can be compromised, especially when small, single cells are targeted. To circumvent this problem, we have established a workflow for LAM that reproducibly results in high RNA quality that is well suitable for transcriptomics, as exemplified here by the isolation of cells of the female germline in apomictic Boechera. In this protocol, procedures are described for tissue preparation and LAM, also with regard to RNA extraction and quality control.

Heritable gene expression differences between apomictic clone members in Taraxacum officinale: Insights into early stages of evolutionary divergence in asexual plants

BACKGROUND

Asexual reproduction has the potential to enhance deleterious mutation accumulation and to constrain adaptive evolution. One source of mutations that can be especially relevant in recent asexuals is activity of transposable elements (TEs), which may have experienced selection for high transposition rates in sexual ancestor populations. Predictions of genomic divergence under asexual reproduction therefore likely include a large contribution of transposable elements but limited adaptive divergence. For plants empirical insight into genome divergence under asexual reproduction remains limited. Here, we characterize expression divergence between clone members of a single apomictic lineage of the common dandelion (Taraxacum officinale) to contribute to our knowledge of genome evolution under asexuality.

RESULTS

Using RNA-Seq, we show that about one third of heritable divergence within the apomictic lineage is driven by TEs and TE-related gene activity. In addition, we identify non-random transcriptional differences in pathways related to acyl-lipid and abscisic acid metabolisms which might reflect functional divergence within the apomictic lineage. We analyze SNPs in the transcriptome to assess genetic divergence between the apomictic clone members and reveal that heritable expression differences between the accessions are not explained simply by genome-wide genetic divergence.

CONCLUSION

The present study depicts a first effort towards a more complete understanding of apomictic plant genome evolution. We identify abundant TE activity and ecologically relevant functional genes and pathways affecting heritable within-lineage expression divergence. These findings offer valuable resources for future work looking at epigenetic silencing and Cis-regulation of gene expression with particular emphasis on the effects of TE activity on asexual species' genome.

Copy number variation of a gene cluster encoding endopolygalacturonase mediates flesh texture and stone adhesion in peach

Texture is an important attribute affecting consumer perception of fruit quality. Peach melting flesh and flesh adhesion to stone (endocarp) are simply inherited and controlled by the F-M locus on linkage group (LG) 4. Here, we report that two genes encoding endopolygalacturonase (endoPG) in the F-M locus, designated PpendoPGF and PpendoPGM, are associated with the melting flesh and stone adhesion traits. PpendoPGM controls melting flesh while PpendoPGF has pleiotropic effects on both melting flesh and stone adhesion. The F-M locus has three allelic copy number variants of endoPG, H1 (PpendoPGF and PpendoPGM), H2 (PpendoPGM), and H3 (null). The H2 haplotype represents the ancestral one while the H1 and H3 haplotypes are two variants due to duplication and deletion of PpendoPGM, respectively. Accessions with H1H1, H1H2, or H1H3 genotypes show the freestone or semi-freestone and melting flesh phenotype, while both H2H2 and H2H3 accessions have the clingstone and melting flesh phenotype. The H3H3 accessions have the clingstone and non-melting flesh phenotype. Our study not only demonstrates a driving role of gene copy number variations in flesh texture diversification in fruit trees, but also provides a useful diagnostic tool for early seedling selection in peach breeding programmes.

Hybrid apomicts trapped in the ecological niches of their sexual ancestors

Asexual reproduction is expected to reduce the adaptive potential to novel or changing environmental conditions, restricting or altering the ecological niche of asexual lineages. Asexual lineages of plants and animals are typically polyploid, an attribute that may influence their genetic variation, plasticity, adaptive potential, and niche breadth. The genus Boechera (Brassicaceae) represents an ideal model to test the relative ecological and biogeographic impacts of reproductive mode and ploidy because it is composed of diploid sexual and both diploid and polyploid asexual (i.e., apomictic) lineages. Here, we demonstrate a strong association between a transcriptionally conserved allele and apomictic seed formation. We then use this allele as a proxy apomixis marker in 1,649 accessions to demonstrate that apomixis is likely to be a common feature across the Boechera phylogeny. Phylogeographic analyses of these data demonstrate (i) species-specific niche differentiation in sexuals, (ii) extensive niche conservation between differing reproductive modes of the same species, (iii) ploidy-specific niche differentiation within and among species, and (iv) occasional niche drift between apomicts and their sexual ancestors. We conclude that ploidy is a substantially stronger and more common driver of niche divergence within and across Boechera species although variation in both traits may not necessarily lead to niche evolution on the species scale.

A conserved apomixis-specific polymorphism is correlated with exclusive exonuclease expression in premeiotic ovules of apomictic boechera species

Apomixis (asexual seed production) is characterized by meiotically unreduced egg cell production (apomeiosis) followed by its parthenogenetic development into offspring that are genetic clones of the mother plant. Fertilization (i.e. pseudogamy) of the central cell is important for the production of a functional endosperm with a balanced 2:1 maternal:paternal genome ratio. Here, we present the APOLLO (for apomixis-linked locus) gene, an Aspartate Glutamate Aspartate Aspartate histidine exonuclease whose transcripts are down-regulated in sexual ovules entering meiosis while being up-regulated in apomeiotic ovules at the same stage of development in plants of the genus Boechera. APOLLO has both "apoalleles," which are characterized by a set of linked apomixis-specific polymorphisms, and "sexalleles." All apomictic Boechera spp. accessions proved to be heterozygous for the APOLLO gene (having at least one apoallele and one sexallele), while all sexual genotypes were homozygous for sexalleles. Apoalleles contained a 20-nucleotide polymorphism present in the 5' untranslated region that contains specific transcription factor-binding sites for ARABIDOPSIS THALIANA HOMEOBOX PROTEIN5, LIM1 (for LINEAGE ABNORMAL11, INSULIN1, MECHANOSENSORY PROTEIN3), SORLIP1AT (for SEQUENCES OVERREPRESENTED IN LIGHT-INDUCED PROMOTERS IN ARABIDOPSIS THALIANA1), SORLIP2AT, and POLYA SIGNAL1. In the same region, sexalleles contain transcription factor-binding sites for DNA BINDING WITH ONE FINGER2, DNA BINDING WITH ONE FINGER3, and PROLAMIN BOX-BINDING FACTOR. Our results suggest that the expression of a single deregulated allele could induce the cascade of events leading to asexual female gamete formation in an apomictic plant.

The conserved chimeric transcript UPGRADE2 is associated with unreduced pollen formation and is exclusively found in apomictic Boechera species

In apomictic Boechera spp., meiotic diplospory leads to the circumvention of meiosis and the suppression of recombination to produce unreduced male and female gametes (i.e. apomeiosis). Here, we have established an early flower developmental staging system and have performed microarray-based comparative gene expression analyses of the pollen mother cell stage in seven diploid sexual and seven diploid apomictic genotypes to identify candidate factors for unreduced pollen formation. We identified a transcript unique to apomictic Boechera spp. called UPGRADE2 (BspUPG2), which is highly up-regulated in their pollen mother cells. BspUPG2 is highly conserved among apomictic Boechera spp. genotypes but has no homolog in sexual Boechera spp. or in any other taxa. BspUPG2 undergoes posttranscriptional processing but lacks a prominent open reading frame. Together with the potential of stably forming microRNA-like secondary structures, we hypothesize that BspUPG2 functions as a long regulatory noncoding messenger RNA-like RNA. BspUPG2 has apparently arisen through a three-step process initiated by ancestral gene duplication of the original BspUPG1 locus, followed by sequential insertions of segmentally duplicated gene fragments, with final exonization of its sequence structure. Its genesis reflects the hybridization history that characterizes the genus Boechera.