A phylotranscriptomic analysis of gene family expansion and evolution in the largest order of pleurocarpous mosses (Hypnales, Bryophyta)

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Fatty Acid Profiles of Some Siberian Bryophytes and Prospects of Their Use in Chemotaxonomy

The composition of fatty acids (FAs) in gametophyte samples of 20 Siberian bryophyte species from four orders of mosses and four orders of liverworts collected in relatively cold months (April and/or October) was examined. FA profiles were obtained using gas chromatography. Thirty-seven FAs were found, from 12:0 to 26:0; they included mono-, polyunsaturated (PUFAs) and rare FAs, such as 22:5n-3 and two acetylenic FAs, 6a,9,12-18:3 and 6a,9,12,15-18:4 (dicranin). Acetylenic FAs were found in all examined species of the Bryales and Dicranales orders, dicranin being the predominant FA. The role of particular PUFAs in mosses and liverworts is discussed. Multivariate discriminant analysis (MDA) was performed to determine whether FAs can be used in the chemotaxonomy of bryophytes. Based on the MDA results, FA composition is related to the taxonomic status of species. Thus, several individual FAs were identified as chemotaxonomic markers at the level of bryophyte orders. These were 18:3n-3; 18:4n-3; 6a,9,12-18:3; 6a,9,12,15-18:4; 20:4n-3 and EPA in mosses and 16:3n-3; 16:2n-6; 18:2n-6; 18:3n-3 and EPA in liverworts. These findings indicate that further research into bryophyte FA profiles can shed light on phylogenetic relationships within this group of plants and the evolution of their metabolic pathways.

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

On the potential of Angiosperms353 for population genomic studies

PREMISE

The successful application of universal targeted sequencing markers, such as those developed for the Angiosperms353 probe set, within populations could reduce or eliminate the need for specific marker development, while retaining the benefits of full-gene sequences in population-level analyses. However, whether the Angiosperms353 markers provide sufficient variation within species to calculate demographic parameters is untested.

METHODS

Using herbarium specimens from a 50-year-old floristic survey in Texas, we sequenced 95 samples from 24 species using the Angiosperms353 probe set. Our data workflow calls variants within species and prepares data for population genetic analysis using standard metrics. In our case study, gene recovery was affected by genomic library concentration only at low concentrations and displayed limited phylogenetic bias.

RESULTS

We identified over 1000 segregating variants with zero missing data for 92% of species and demonstrate that Angiosperms353 markers contain sufficient variation to estimate pairwise nucleotide diversity (π)-typically between 0.002 and 0.010, with most variation found in flanking non-coding regions. In a subset of variants that were filtered to reduce linkage, we uncovered high heterozygosity in many species, suggesting that denser sampling within species should permit estimation of gene flow and population dynamics.

DISCUSSION

Angiosperms353 should benefit conservation genetic studies by providing universal repeatable markers, low missing data, and haplotype information, while permitting inclusion of decades-old herbarium specimens.

Gene-rich UV sex chromosomes harbor conserved regulators of sexual development

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

Charting the genomic landscape of seed-free plants

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae

A new paradigm has slowly emerged regarding the diversification of bryophytes, with inferences from molecular data highlighting a dynamic evolution of their genome. However, comparative studies of expressed genes among closely related taxa is so far missing. Here we contrast the dimensions of the vegetative transcriptome of Funaria hygrometrica and Physcomitrium pyriforme against the genome of their relative, Physcomitrium (Physcomitrella) patens. These three species of Funariaceae share highly conserved vegetative bodies, and are partially sympatric, growing on mineral soil in mostly temperate regions. We analyzed the vegetative gametophytic transcriptome of F. hygrometrica and P. pyriforme and mapped short reads, transcripts, and proteins to the genome and gene space of P. patens. Only about half of the transcripts of F. hygrometrica map to their ortholog in P. patens, whereas at least 90% of those of P. pyriforme align to loci in P. patens. Such divergence is unexpected given the high morphological similarity of the gametophyte but reflects the estimated times of divergence of F. hygrometrica and P. pyriforme from P. patens, namely 55 and 20 mya, respectively. The newly sampled transcriptomes bear signatures of at least one, rather ancient, whole genome duplication (WGD), which may be shared with one reported for P. patens. The transcriptomes of F. hygrometrica and P. pyriforme reveal significant contractions or expansions of different gene families. While transcriptomes offer only an incomplete estimate of the gene space, the high number of transcripts obtained suggest a significant divergence in gene sequences, and gene number among the three species, indicative of a rather strong, dynamic genome evolution, shaped in part by whole, partial or localized genome duplication. The gene ontology of their specific and rapidly-evolving protein families, suggests that the evolution of the Funariaceae may have been driven by the diversification of metabolic genes that may optimize the adaptations to environmental conditions, a hypothesis well in line with ecological patterns in the genetic diversity and structure in seed plants.

The Ceratodon purpureus genome uncovers structurally complex, gene rich sex chromosomes.. [PMID: 0 DOI: 10.1101/2020.07.03.163634]

Non-recombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration1,2. The correlation between suppressed recombination and degeneration is clear in animal XY systems1,2, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions3,4. Here we test for degeneration in the bryophyte UV sex chromosome system through genomic comparisons with new female and male chromosome-scale reference genomes of the moss Ceratodon purpureus. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes show signs of weaker purifying selection than autosomes, we find suppressed recombination alone is insufficient to drive gene loss on sex-specific chromosomes. Instead, the U and V sex chromosomes harbor thousands of broadly-expressed genes, including numerous key regulators of sexual development across land plants.

Phylotranscriptomic Analyses Reveal Asymmetrical Gene Duplication Dynamics and Signatures of Ancient Polyploidy in Mints

Ancient duplication events and retained gene duplicates have contributed to the evolution of many novel plant traits and, consequently, to the diversity and complexity within and across plant lineages. Although mounting evidence highlights the importance of whole-genome duplication (WGD; polyploidy) and its key role as an evolutionary driver, gene duplication dynamics and mechanisms, both of which are fundamental to our understanding of evolutionary process and patterns of plant diversity, remain poorly characterized in many clades. We use newly available transcriptomic data and a robust phylogeny to investigate the prevalence, occurrence, and timing of gene duplications in Lamiaceae (mints), a species-rich and chemically diverse clade with many ecologically, economically, and culturally important species. We also infer putative WGDs—an extreme mechanism of gene duplication—using large-scale data sets from synonymous divergence (K_S), phylotranscriptomic, and divergence time analyses. We find evidence for widespread but asymmetrical levels of gene duplication and ancient polyploidy in Lamiaceae that correlate with species richness, including pronounced levels of gene duplication and putative ancient WGDs (7–18 events) within the large subclade Nepetoideae and up to 10 additional WGD events in other subclades. Our results help disentangle WGD-derived gene duplicates from those produced by other mechanisms and illustrate the nonuniformity of duplication dynamics in mints, setting the stage for future investigations that explore their impacts on trait diversity and species diversification. Our results also provide a practical context for evaluating the benefits and limitations of transcriptome-based approaches to inferring WGD, and we offer recommendations for researchers interested in investigating ancient WGDs in other plant groups.

Genome size and endopolyploidy evolution across the moss phylogeny

BACKGROUND AND AIMS

Compared with other plant lineages, bryophytes have very small genomes with little variation across species, and high levels of endopolyploid nuclei. This study is the first analysis of moss genome evolution over a broad taxonomic sampling using phylogenetic comparative methods. We aim to determine whether genome size evolution is unidirectional as well as examine whether genome size and endopolyploidy are correlated in mosses.

METHODS

Genome size and endoreduplication index (EI) estimates were newly generated using flow cytometry from moss samples collected in Canada. Phylogenetic relationships between moss species were reconstructed using GenBank sequence data and maximum likelihood methods. Additional 1C-values were compiled from the literature and genome size and EI were mapped onto the phylogeny to reconstruct ancestral character states, test for phylogenetic signal and perform phylogenetic independent contrasts.

KEY RESULTS

Genome size and EI were obtained for over 50 moss taxa. New genome size estimates are reported for 33 moss species and new EIs are reported for 20 species. In combination with data from the literature, genome sizes were mapped onto a phylogeny for 173 moss species with this analysis, indicating that genome size evolution in mosses does not appear to be unidirectional. Significant phylogenetic signal was detected for genome size when evaluated across the phylogeny, whereas phylogenetic signal was not detected for EI. Genome size and EI were not found to be significantly correlated when using phylogenetically corrected values.

CONCLUSIONS

Significant phylogenetic signal indicates closely related mosses have more similar genome sizes and EI values. This study supports that DNA content in mosses is defined by small genomes that are highly endopolyploid, suggesting strong selective pressure to maintain these features. Further research is needed to understand the functional significance of DNA content evolution in mosses.

GingerRoot: A Novel DNA Transposon Encoding Integrase-Related Transposase in Plants and Animals

Transposable elements represent the largest components of many eukaryotic genomes and different genomes harbor different combinations of elements. Here, we discovered a novel DNA transposon in the genome of the clubmoss Selaginella lepidophylla. Further searching for related sequences to the conserved DDE region uncovered the presence of this superfamily of elements in fish, coral, sea anemone, and other animal species. However, this element appears restricted to Bryophytes and Lycophytes in plants. This transposon, named GingerRoot, is associated with a 6 bp (base pair) target site duplication, and 100-150 bp terminal inverted repeats. Analysis of transposase sequences identified the DDE motif, a catalytic domain, which shows similarity to the integrase of Gypsy-like long terminal repeat retrotransposons, the most abundant component in plant genomes. A total of 77 intact and several hundred truncated copies of GingerRoot elements were identified in S. lepidophylla. Like Gypsy retrotransposons, GingerRoots show a lack of insertion preference near genes, which contrasts to the compact genome size of about 100 Mb. Nevertheless, a considerable portion of GingerRoot elements was found to carry gene fragments, suggesting the capacity of duplicating gene sequences is unlikely attributed to the proximity to genes. Elements carrying gene fragments appear to be less methylated, more diverged, and more distal to genes than those without gene fragments, indicating they are preferentially retained in gene-poor regions. This study has identified a broadly dispersed, novel DNA transposon, and the first plant DNA transposon with an integrase-related transposase, suggesting the possibility of de novo formation of Gypsy-like elements in plants.

Genome Sequencing of Pleurozium schreberi: The Assembled and Annotated Draft Genome of a Pleurocarpous Feather Moss

The pleurocarpous feather moss Pleurozium schreberi is a ubiquitous moss species which plays a fundamental role in many terrestrial ecosystems, for instance within the boreal forest, the Earth's largest terrestrial biome, this species plays a significant role in driving ecosystem nitrogen and carbon inputs and fluxes. By hosting dinitrogen (N₂)-fixing cyanobacteria, the moss-cyanobacteria symbiosis constitutes the main nitrogen input into the ecosystem and by the high productivity and the low decomposability of the moss litter, P schreberi contributes significantly to build-up soil organic matter, and therefore long-term C sequestration. Knowledge on P. schreberi genome will facilitate the development of 'omics' and system's biology approaches to gain a more complete understanding of the physiology and ecological adaptation of the moss and the mechanisms underpinning the establishment of the symbiosis. Here we present the de novo assembly and annotation of P. schreberi genome that will help investigating these questions. The sequencing was performed using the HiSeq X platform with Illumina paired-end and mate-pair libraries prepared with CTAB extracted DNA. In total, the assembled genome was approximately 318 Mb, while repetitive elements account for 28.42% of the genome and 15,992 protein-coding genes were predicted from the genome, of which 84.23% have been functionally annotated. We anticipate that the genomic data generated will constitute a significant resource to study ecological and evolutionary genomics of P. schreberi, and will be valuable for evo-devo investigations as well as our understanding of the evolution of land plants by providing the genome of a pleurocarpous moss.

Assessing the Performance of Ks Plots for Detecting Ancient Whole Genome Duplications

Genomic data have provided evidence of previously unknown ancient whole genome duplications (WGDs) and highlighted the role of WGDs in the evolution of many eukaryotic lineages. Ancient WGDs often are detected by examining distributions of synonymous substitutions per site (Ks) within a genome, or “Ks plots.” For example, WGDs can be detected from Ks plots by using univariate mixture models to identify peaks in Ks distributions. We performed gene family simulation experiments to evaluate the effects of different Ks estimation methods and mixture models on our ability to detect ancient WGDs from Ks plots. The simulation experiments, which accounted for variation in substitution rates and gene duplication and loss rates across gene families, tested the effects of WGD age and gene retention rates following WGD on inferring WGDs from Ks plots. Our simulations reveal limitations of Ks plot analyses. Strict interpretations of mixture model analyses often overestimate the number of WGD events, and Ks plot analyses typically fail to detect WGDs when ≤10% of the duplicated genes are retained following the WGD. However, WGDs can accurately be characterized over an intermediate range of Ks. The simulation results are supported by empirical analyses of transcriptomic data, which also suggest that biases in gene retention likely affect our ability to detect ancient WGDs. Although our results indicate mixture model results should be interpreted with great caution, using node-averaged Ks estimates and applying more appropriate mixture models can improve the accuracy of detecting WGDs.

Phylogenomics reveals an extensive history of genome duplication in diatoms (Bacillariophyta)

PREMISE OF THE STUDY

Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well.

METHODS

We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species.

KEY RESULTS

Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades.

CONCLUSIONS

Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity.

Evolutionary dynamism in bryophytes: Phylogenomic inferences confirm rapid radiation in the moss family Funariaceae

Rapid diversifications of plants are primarily documented and studied in angiosperms, which are perceived as evolutionarily dynamic. Recent studies have, however, revealed that bryophytes have also undergone periods of rapid radiation. The speciose family Funariaceae, including the model taxon Physcomitrella patens, is one such lineage. Here, we infer relationships among major lineages within the Entosthodon-Physcomitrium complex from virtually complete organellar exomes (i.e., 123 genes) obtained through high throughput sequencing of genomic libraries enriched in these loci via targeted locus capture. Based on these extensive exonic data we (1) reconstructed a robust backbone topology of the Funariaceae, (2) confirmed the monophyly of Funaria and the polyphyly of Entosthodon, Physcomitrella, and Physcomitrium, and (3) argue for the occurrence of a rapid radiation within the Entosthodon-Physcomitrium complex that began 28 mya and gave rise more than half of the species diversity of the family. This diversification may have been triggered by a whole genome duplication and coincides with global Eocene cooling that continued through the Oligocene and Miocene. The Funariaceae join a growing list of bryophyte lineages whose history is marked by at least one burst of diversification, and our study thereby strengthens the view that bryophytes are evolutionarily dynamic lineages and that patterns and processes characterizing the evolution of angiosperms may be universal among land plants.

Multiple innovations underpinned branching form diversification in mosses

Broad-scale evolutionary comparisons have shown that branching forms arose by convergence in vascular plants and bryophytes, but the trajectory of branching form diversification in bryophytes is unclear. Mosses are the most species-rich bryophyte lineage and two sub-groups are circumscribed by alternative reproductive organ placements. In one, reproductive organs form apically, terminating growth of the primary shoot (gametophore) axis. In the other, reproductive organs develop on very short lateral branches. A switch from apical to lateral reproductive organ development is proposed to have primed branching form diversification. Moss gametophores have modular development and each module develops from a single apical cell. Here we define the architectures of 175 mosses by the number of module classes, branching patterns and the pattern in which similar modules repeat. Using ancestral character state reconstruction we identify two stages of architectural diversification. During a first stage there were sequential changes in the module repetition pattern, reproductive organ position, branching pattern and the number of module classes. During a second stage, vegetative changes occurred independently of reproductive fate. The results pinpoint the nature of developmental change priming branching form diversification in mosses and provide a framework for mechanistic studies of architectural diversification.