Phylotranscriptomic analysis and genome evolution of the Cypripedioideae (Orchidaceae)

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

Statistical Genomics Analysis of Simple Sequence Repeats from the Paphiopedilum Malipoense Transcriptome Reveals Control Knob Motifs Modulating Gene Expression

Simple sequence repeats (SSRs) are found in nonrandom distributions in genomes and are thought to impact gene expression. The distribution patterns of 48 295 SSRs of Paphiopedilum malipoense are mined and characterized based on the first full-length transcriptome and comprehensive transcriptome dataset from 12 organs. Statistical genomics analyses are used to investigate how SSRs in transcripts affect gene expression. The results demonstrate the correlations between SSR distributions, characteristics, and expression level. Nine expression-modulating motifs (expMotifs) are identified and a model is proposed to explain the effect of their key features, potency, and gene function on an intra-transcribed region scale. The expMotif-transcribed region combination is the most predominant contributor to the expression-modulating effect of SSRs, and some intra-transcribed regions are critical for this effect. Genes containing the same type of expMotif-SSR elements in the same transcribed region are likely linked in function, regulation, or evolution aspects. This study offers novel evidence to understand how SSRs regulate gene expression and provides potential regulatory elements for plant genetic engineering.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Integrated evolutionary pattern analyses reveal multiple origins of steroidal saponins in plants

Steroidal saponins are a class of specialized metabolites essential for plant's response to biotic and abiotic stresses. They are also important raw materials for the industrial production of steroid drugs. Steroidal saponins are present in some monocots, such as Dioscorea and Paris, but their distribution, origin, and evolution in plants remain poorly understood. By reconstructing the evolutionary history of the steroidal saponin-associated module (SSAM) in plants, we reveal that the steroidal saponin pathway has its origin in Asparagus and Dioscorea. Through evaluating the distribution and evolutionary pattern of steroidal saponins in angiosperms, we further show that steroidal saponins originated multiple times in angiosperms, and exist in early diverged lineages of certain monocot lineages including Asparagales, Dioscoreales, and Liliales. In these lineages, steroidal saponins are synthesized through the high copy and/or high expression mechanisms of key genes in SSAM. Together with shifts in gene evolutionary rates and amino acid usage, these molecular mechanisms shape the current distribution and diversity of steroidal saponins in plants. Consequently, our results provide new insights into the distribution, diversity and evolutionary history of steroidal saponins in plants, and enhance our understanding of plants' resistance to abiotic and biotic stresses. Additionally, fundamental understanding of the steroidal saponin biosynthesis will facilitate their industrial production and pharmacological applications.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Advances and prospects of orchid research and industrialization

Orchidaceae is one of the largest, most diverse families in angiosperms with significant ecological and economical values. Orchids have long fascinated scientists by their complex life histories, exquisite floral morphology and pollination syndromes that exhibit exclusive specializations, more than any other plants on Earth. These intrinsic factors together with human influences also make it a keystone group in biodiversity conservation. The advent of sequencing technologies and transgenic techniques represents a quantum leap in orchid research, enabling molecular approaches to be employed to resolve the historically interesting puzzles in orchid basic and applied biology. To date, 16 different orchid genomes covering four subfamilies (Apostasioideae, Vanilloideae, Epidendroideae, and Orchidoideae) have been released. These genome projects have given rise to massive data that greatly empowers the studies pertaining to key innovations and evolutionary mechanisms for the breadth of orchid species. The extensive exploration of transcriptomics, comparative genomics, and recent advances in gene engineering have linked important traits of orchids with a multiplicity of gene families and their regulating networks, providing great potential for genetic enhancement and improvement. In this review, we summarize the progress and achievement in fundamental research and industrialized application of orchids with a particular focus on molecular tools, and make future prospects of orchid molecular breeding and post-genomic research, providing a comprehensive assemblage of state of the art knowledge in orchid research and industrialization.

Mesoamerican Cypripedium: Mycorrhizal Contributions to Promote Their Conservation as Critically Endangered Species

In the valuable orchid genus Cypripedium, the section Irapeana consists of a distinctive group of Mesoamerican species that is formed by Cypripedium dickinsonianum Hágsater, C. irapeanum Lex., and C. molle Lindl. All lady slipper orchids exhibit different distributions and abundances. Data analysis that used herbarium accessions and field investigations indicated that the habitats of these three species have been dramatically reduced. Prospecting for suitable habitats based on climatic, vegetation, and soil parameters allows us to predict potential distributions. Conservation strategies, such as ex situ propagation by asymbiotic and symbiotic approaches, have indicated that the culture media used are a determining factor for seedling development. Mycorrhizal isolates play a main role in the compatibility and further development of germinated seeds. The fungi isolated from adult plants belong to two different families, which makes it possible that widely distributed C. irapeanum populations will be fungal-specific as well as restricted for C. molle. Root mycorrhization patterns occur high on the secondary roots. In contrast with other species of the genus, in situ germination can occur over a short period of two months, but we have documented periods as long as ten years. Cypripedium is a highly problematic genus for ex situ conservation because the germination requirements and cultures are poorly documented, and there is great urgency for in situ conservation to develop strategies for identifying hotspot habitats and actualize the protection status to avoid extinction of this genus.

Orchid Phylotranscriptomics: The Prospects of Repurposing Multi-Tissue Transcriptomes for Phylogenetic Analysis and Beyond

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed re-assembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer.

Target sequence capture in orchids: Developing a kit to sequence hundreds of single-copy loci

PREMISE

Understanding relationships among orchid species and populations is of critical importance for orchid conservation. Target sequence capture has become a standard method for extracting hundreds of orthologous loci for phylogenomics. Up-front cost and time associated with design of bait sets makes this method prohibitively expensive for many researchers. Therefore, we designed a target capture kit to reliably sequence hundreds of orthologous loci across orchid lineages.

METHODS

We designed an Orchidaceae target capture bait set for 963 single-copy genes identified in published orchid genome sequences. The bait set was tested on 28 orchid species, with representatives of the subfamilies Cypripedioideae, Orchidoideae, and Epidendroideae.

RESULTS

Between 1,518,041 and 87,946,590 paired-end 150-base reads were generated for target-enriched genomic libraries. We assembled an average of 812 genes per library for Epidendroideae species and a mean of 501 genes for species in the subfamilies Orchidoideae and Cypripedioideae. Furthermore, libraries had on average 107 of the 254 genes that are included in the Angiosperms353 bait set, allowing for direct comparison of studies using either bait set.

DISCUSSION

The Orchidaceae963 kit will enable greater accessibility and utility of next-generation sequencing for orchid systematics, population genetics, and identification in the illegal orchid trade.

Hundreds of nuclear and plastid loci yield novel insights into orchid relationships

PREMISE

The inference of evolutionary relationships in the species-rich family Orchidaceae has hitherto relied heavily on plastid DNA sequences and limited taxon sampling. Previous studies have provided a robust plastid phylogenetic framework, which was used to classify orchids and investigate the drivers of orchid diversification. However, the extent to which phylogenetic inference based on the plastid genome is congruent with the nuclear genome has been only poorly assessed.

METHODS

We inferred higher-level phylogenetic relationships of orchids based on likelihood and ASTRAL analyses of 294 low-copy nuclear genes sequenced using the Angiosperms353 universal probe set for 75 species (representing 69 genera, 16 tribes, 24 subtribes) and a concatenated analysis of 78 plastid genes for 264 species (117 genera, 18 tribes, 28 subtribes). We compared phylogenetic informativeness and support for the nuclear and plastid phylogenetic hypotheses.

RESULTS

Phylogenetic inference using nuclear data sets provides well-supported orchid relationships that are highly congruent between analyses. Comparisons of nuclear gene trees and a plastid supermatrix tree showed that the trees are mostly congruent, but revealed instances of strongly supported phylogenetic incongruence in both shallow and deep time. The phylogenetic informativeness of individual Angiosperms353 genes is in general better than that of most plastid genes.

CONCLUSIONS

Our study provides the first robust nuclear phylogenomic framework for Orchidaceae and an assessment of intragenomic nuclear discordance, plastid-nuclear tree incongruence, and phylogenetic informativeness across the family. Our results also demonstrate what has long been known but rarely thoroughly documented: nuclear and plastid phylogenetic trees can contain strongly supported discordances, and this incongruence must be reconciled prior to interpretation in evolutionary studies, such as taxonomy, biogeography, and character evolution.

A multitiered sequence capture strategy spanning broad evolutionary scales: Application for phylogenetic and phylogeographic studies of orchids

With over 25,000 species, the drivers of diversity in the Orchidaceae remain to be fully understood. Here, we outline a multitiered sequence capture strategy aimed at capturing hundreds of loci to enable phylogenetic resolution from subtribe to subspecific levels in orchids of the tribe Diurideae. For the probe design, we mined subsets of 18 transcriptomes, to give five target sequence sets aimed at the tribe (Sets 1 & 2), subtribe (Set 3), and within subtribe levels (Sets 4 & 5). Analysis included alternative de novo and reference-guided assembly, before target sequence extraction, annotation and alignment, and application of a homology-aware k-mer block phylogenomic approach, prior to maximum likelihood and coalescence-based phylogenetic inference. Our evaluation considered 87 taxa in two test data sets: 67 samples spanning the tribe, and 72 samples involving 24 closely related Caladenia species. The tiered design achieved high target loci recovery (>89%), with the median number of recovered loci in Sets 1-5 as follows: 212, 219, 816, 1024, and 1009, respectively. Interestingly, as a first test of the homologous k-mer approach for targeted sequence capture data, our study revealed its potential for enabling robust phylogenetic species tree inferences. Specifically, we found matching, and in one case improved phylogenetic resolution within species complexes, compared to conventional phylogenetic analysis involving target gene extraction. Our findings indicate that a customized multitiered sequence capture strategy, in combination with promising yet underutilized phylogenomic approaches, will be effective for groups where interspecific divergence is recent, but information on deeper phylogenetic relationships is also required.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

Plastome structure and adaptive evolution of Calanthe s.l. species

Calanthe s.l. is the most diverse group in the tribe Collabieae (Orchidaceae), which are pantropical in distribution. Illumina sequencing followed by de novo assembly was used in this study, and the plastid genetic information of Calanthe s.l. was used to investigate the adaptive evolution of this taxon. Herein, the complete plastome of five Calanthe s.l. species (Calanthe davidii, Styloglossum lyroglossa, Preptanthe rubens, Cephalantheropsis obcordata, and Phaius tankervilliae) were determined, and the two other published plastome sequences of Calanthe s.l. were added for comparative analyses to examine the evolutionary pattern of the plastome in the alliance. The seven plastomes ranged from 150,181 bp (C. delavayi) to 159,014 bp (C. davidii) in length and were all mapped as circular structures. Except for the three ndh genes (ndhC, ndhF, and ndhK) lost in C. delavayi, the remaining six species contain identical gene orders and numbers (115 gene). Nucleotide diversity was detected across the plastomes, and we screened 14 mutational hotspot regions, including 12 non-coding regions and two gene regions. For the adaptive evolution investigation, three species showed positive selected genes compared with others, C. obcordata (cemA), S. lyroglossa (infA, ycf1 and ycf2) and C. delavayi (nad6 and ndhB). Six genes were under site-specific positive selection in Calanthe s.l., namely, accD, ndhB, ndhD, rpoC2, ycf1, and ycf2, most of which are involved in photosynthesis. These results, including the new plastomes, provide resources for the comparative plastome, breeding, and plastid genetic engineering of orchids and flowering plants.

Phylogeny and multiple independent whole-genome duplication events in the Brassicales

PREMISE

Whole-genome duplications (WGDs) are prevalent throughout the evolutionary history of plants. For example, dozens of WGDs have been phylogenetically localized across the order Brassicales, specifically, within the family Brassicaceae. A WGD event has also been identified in the Cleomaceae, the sister family to Brassicaceae, yet its placement, as well as that of WGDs in other families in the order, remains unclear.

METHODS

Phylo-transcriptomic data were generated and used to infer a nuclear phylogeny for 74 Brassicales taxa. Genome survey sequencing was also performed on 66 of those taxa to infer a chloroplast phylogeny. These phylogenies were used to assess and confirm relationships among the major families of the Brassicales and within Brassicaceae. Multiple WGD inference methods were then used to assess the placement of WGDs on the nuclear phylogeny.

RESULTS

Well-supported chloroplast and nuclear phylogenies for the Brassicales and the putative placement of the Cleomaceae-specific WGD event Th-ɑ are presented. This work also provides evidence for previously hypothesized WGDs, including a well-supported event shared by at least two members of the Resedaceae family, and a possible event within the Capparaceae.

CONCLUSIONS

Phylogenetics and the placement of WGDs within highly polyploid lineages continues to be a major challenge. This study adds to the conversation on WGD inference difficulties by demonstrating that sampling is especially important for WGD identification and phylogenetic placement. Given its economic importance and genomic resources, the Brassicales continues to be an ideal group for assessing WGD inference methods.

Gene and genome duplications in the evolution of chemodiversity: perspectives from studies of Lamiaceae

Plants are reservoirs of extreme chemical diversity, yet biosynthetic pathways remain underexplored in the majority of taxa. Access to improved, inexpensive genomic and computational technologies has recently enhanced our understanding of plant specialized metabolism at the biochemical and evolutionary levels including the elucidation of pathways leading to key metabolites. Furthermore, these approaches have provided insights into the mechanisms of chemical evolution, including neofunctionalization and subfunctionalization, structural variation, and modulation of gene expression. The broader utilization of genomic tools across the plant tree of life, and an expansion of genomic resources from multiple accessions within species or populations, will improve our overall understanding of chemodiversity. These data and knowledge will also lead to greater insight into the selective pressures contributing to and maintaining this diversity, which in turn will enable the development of more accurate predictive models of specialized metabolism in 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.

Ancient Polyploidy and Genome Evolution in Palms

Mechanisms of genome evolution are fundamental to our understanding of adaptation and the generation and maintenance of biodiversity, yet genome dynamics are still poorly characterized in many clades. Strong correlations between variation in genomic attributes and species diversity across the plant tree of life suggest that polyploidy or other mechanisms of genome size change confer selective advantages due to the introduction of genomic novelty. Palms (order Arecales, family Arecaceae) are diverse, widespread, and dominant in tropical ecosystems, yet little is known about genome evolution in this ecologically and economically important clade. Here, we take a phylogenetic comparative approach to investigate palm genome dynamics using genomic and transcriptomic data in combination with a recent, densely sampled, phylogenetic tree. We find conclusive evidence of a paleopolyploid event shared by the ancestor of palms but not with the sister clade, Dasypogonales. We find evidence of incremental chromosome number change in the palms as opposed to one of recurrent polyploidy. We find strong phylogenetic signal in chromosome number, but no signal in genome size, and further no correlation between the two when correcting for phylogenetic relationships. Palms thus add to a growing number of diverse, ecologically successful clades with evidence of whole-genome duplication, sister to a species-poor clade with no evidence of such an event. Disentangling the causes of genome size variation in palms moves us closer to understanding the genomic conditions facilitating adaptive radiation and ecological dominance in an evolutionarily successful, emblematic tropical clade.