Analyses of transcriptome sequences reveal multiple ancient large-scale duplication events in the ancestor of Sphagnopsida (Bryophyta)

Genome-Wide Characterization and Expression Profiling of ABA Biosynthesis Genes in a Desert Moss Syntrichia caninervis

Syntrichia caninervis can survive under 80-90% protoplasmic water losses, and it is a model plant in desiccation tolerance research. A previous study has revealed that S. caninervis would accumulate ABA under dehydration stress, while the ABA biosynthesis genes in S. caninervis are still unknown. This study identified one ScABA1, two ScABA4s, five ScNCEDs, twenty-nine ScABA2s, one ScABA3, and four ScAAOs genes, indicating that the ABA biosynthesis genes were complete in S. caninervis. Gene location analysis showed that the ABA biosynthesis genes were evenly distributed in chromosomes but were not allocated to sex chromosomes. Collinear analysis revealed that ScABA1, ScNCED, and ScABA2 had homologous genes in Physcomitrella patens. RT-qPCR detection found that all of the ABA biosynthesis genes responded to abiotic stress; it further indicated that ABA plays an important role in S. caninervis. Moreover, the ABA biosynthesis genes in 19 representative plants were compared to study their phylogenetic and conserved motifs; the results suggested that the ABA biosynthesis genes were closely associated with plant taxa, but these genes had the same conserved domain in each plant. In contrast, there is a huge variation in the exon number between different plant taxa; it revealed that ABA biosynthesis gene structures are closely related to plant taxa. Above all, this study provides strong evidence demonstrating that ABA biosynthesis genes were conserved in the plant kingdom and deepens our understanding of the evolution of the phytohormone ABA.

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.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

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.

Evolution of an intermediate C4 photosynthesis in the non-foliar tissues of the Poaceae

Carbon concentrating mechanisms (CCMs) in plants are abaptive features that have evolved to sustain plant growth in unfavorable environments, especially at low atmospheric carbon levels and high temperatures. Uptake of CO₂ and its storage in the aerenchyma tissues of Lycopsids and diurnal acidity fluctuation in aquatic plants during the Palaeozoic era (ca. 300 Ma.) would represent the earliest evolution of a CCM. The CCM parts of the dark reactions of photosynthesis have evolved many times, while the light reactions are conserved across plant lineages. A C₄ type CCM, leaf C₄ photosynthesis is evolved in the PACMAD clade of the Poaceae family. The evolution of C₄ photosynthesis from C₃ photosynthesis was an abaptation. Photosynthesis in reproductive tissues of sorghum and maize (PACMAD clade) has been shown to be of a weaker C₄ type (high CO₂ compensation point, low carbon isotope discrimination, and lack of Rubisco compartmentalization, when compared to the normal C₄ types) than that in the leaves (normal C₄ type). However, this does not fit well with the character polarity concept from an evolutionary perspective. In a recent model proposed for CCM evolution, the development of a rudimentary CCM prior to the evolution of a more efficient CCM (features contrasting to a weaker C₄ type, leading to greater biomass production rate) has been suggested. An intermediate crassulacean acid metabolism (CAM) type of CCM (rudimentary) was reported in the genera, Brassia, Coryanthes, Eriopsis, Peristeria, of the orchids (well-known group of plants that display the CAM pathway). Similarly, we propose here the evolution of a rudimentary CCM (C₄-like type pathway) in the non-foliar tissues of the Poaceae, prior to the evolution of the C₄ pathway as identified in the leaves of the C₄ species of the PACMAD clade.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Moss transcription factors regulating development and defense responses to stress

Transcription factors control gene expression, leading to regulation of biological processes that determine plant development and adaptation to the environment. Land colonization by plants occurred 450-470 million years ago and was accompanied by an increase in the complexity of transcriptional regulation associated to transcription factor gene expansions. AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY transcription factor families increased in land plants compared with algae. In angiosperms, they play crucial roles in regulating plant growth and responses to environmental stressors. However, less information is available in bryophytes and only in a few cases is the functional role of moss transcription factors in stress mechanisms known. In this review, we discuss current knowledge of the transcription factor families involved in development and defense responses to stress in mosses and other bryophytes. By exploring and analysing the Physcomitrium patens public database and published transcriptional profiles, we show that a high number of AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY genes are differentially expressed in response to abiotic stresses and during biotic interactions. Expression profiles together with a comprehensive analysis provide insights into relevant transcription factors involved in moss defenses, and hint at distinct and conserved biological roles between bryophytes and angiosperms.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

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.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

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.

Phylogenomic synteny network analyses reveal ancestral transpositions of auxin response factor genes in plants

BACKGROUND

Auxin response factors (ARFs) have long been a research focus and represent a class of key regulators of plant growth and development. Integrated phylogenomic synteny network analyses were able to provide novel insights into the evolution of the ARF gene family.

RESULTS

Here, more than 3500 ARFs collected from plant genomes and transcriptomes covering major streptophyte lineages were used to reconstruct the broad-scale family phylogeny, where the early origin and diversification of ARF in charophytes was delineated. Based on the family phylogeny, we proposed a unified six-group classification system for angiosperm ARFs. Phylogenomic synteny network analyses revealed the deeply conserved genomic syntenies within each of the six ARF groups and the interlocking syntenic relationships connecting distinct groups. Recurrent duplication events, such as those that occurred in seed plants, angiosperms, core eudicots and grasses contributed to the expansion of ARF genes which facilitated functional diversification. Ancestral transposition activities in important plant families, including crucifers, legumes and grasses, were unveiled by synteny network analyses. Ancestral gene duplications along with transpositions have profound evolutionary significance which may have accelerated the functional diversification process of paralogues.

CONCLUSIONS

The broad-scale family phylogeny in combination with the state-of-art phylogenomic synteny network analyses not only allowed us to infer the evolutionary trajectory of ARF genes across distinct plant lineages, but also facilitated to generate a more robust classification regime for this transcription factor family. Our study provides insights into the evolution of ARFs which will enhance our current understanding of this important transcription factor family.

Recurrent genome duplication events likely contributed to both the ancient and recent rise of ferns

Ferns, the second largest group of vascular plants, originated ~400 million years ago (Mya). They became dominant in the ancient Earth landscape before the angiosperms and are still important in current ecosystems. Many ferns have exceptionally high chromosome numbers, possibly resulting from whole-genome duplications (WGDs). However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates (Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%-97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strong evidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous (~145-66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.

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.

Genome-wide identification of the CLAVATA3/EMBRYO SURROUNDING REGION (CLE) family in grape (Vitis vinifera L.)

Background

CLE genes play various biological roles in plant growth and development, as well as in responses to environmental stimuli.

Results

In the present study, we identified nine CLE genes in the grape genome using an effective identification method. We analyzed the expression profiles of grape CLE genes in different tissues and under environmental different stimuli. VvCLE3 was expressed in shoot apical meristem (SAM) enriched regions, and VvCLE6 was expressed in shoot tissue without SAM. When grapes were infected with bois noir, VvCLE2 was up-regulated. Under ABA treatment, VvCLE3 was down-regulated. VvCLE6 was up-regulated under high temperature stress. We found that VvCLE6 and VvCLE1 were highly expressed in root tissue. In addition, we compared the characteristics of CLEs from grape and other plant species. The CLE family in Sphagnum fallax underwent positive selection, while the CLE family in grape underwent purifying selection. The frequency of optimal codons and codon adaptation index of rice and grape CLE family members were positively correlated with GC content at the third site of synonymous codons, indicating that the dominant evolutionary pressure acting on rice and grape CLE genes was mutation pressure. We also found that closely related species had higher levels of similarity in relative synonymous codon usage in CLE genes. The rice CLE family was biased toward C and G nucleotides at third codon positions. Gene duplication and loss events were also found in grape CLE genes.

Conclusion

These results demonstrate an effective identification method for CLE motifs and increase the understanding of grape CLEs. Future research on CLE genes may have applications for grape breeding and cultivation to better understand root and nodulation development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5944-2) contains supplementary material, which is available to authorized users.

Ancient duplications and grass-specific transposition influenced the evolution of LEAFY transcription factor genes

The LFY transcription factor gene family are important in the promotion of cell proliferation and floral development. Understanding their evolution offers an insight into floral development in plant evolution. Though a promiscuous transition intermediate and a gene duplication event within the LFY family had been identified previously, the early evolutionary path of this family remained elusive. Here, we reconstructed the LFY family phylogeny using maximum-likelihood and Bayesian inference methods incorporating LFY genes from all major lineages of streptophytes. The well-resolved phylogeny unveiled a high-confidence duplication event before the functional divergence of types I and II LFY genes in the ancestry of liverworts, mosses and tracheophytes, supporting sub-functionalization of an ancestral promiscuous gene. The identification of promiscuous genes in Osmunda suggested promiscuous LFY genes experienced an ancient transient duplication. Genomic synteny comparisons demonstrated a deep genomic positional conservation of LFY genes and an ancestral lineage-specific transposition activity in grasses.

Evolutionary history and structure of nuclear matrix constituent proteins, the plant analogues of lamins

Nuclear matrix constituent proteins (NMCPs), the structural components of the plant lamina, are considered to be the analogues of lamins in plants based on numerous structural and functional similarities. Current phylogenetic knowledge suggests that, in contrast to lamins, which are widely distributed in eukaryotes, NMCPs are taxonomically restricted to Streptophyta. At present, most information about NMCPs comes from angiosperms, and virtually no data are available from more ancestral groups. In angiosperms, the NMCP family comprises two phylogenetic groups, NMCP1 and NMCP2, which evolved from the NMCP1 and NMCP2 progenitor genes. Based on sequence conservation and the presence of NMCP-specific domains, we determined the structure and number of NMCP genes present in different Streptophyta clades. We analysed 91 species of embryophytes and report additional NMCP sequences from mosses, liverworts, clubmosses, horsetail, ferns, gymnosperms, and Charophyta algae. Our results confirm an origin of NMCPs in Charophyta (the earliest diverging group of Streptophyta), resolve the number and structure of NMCPs in the different clades, and propose the emergence of additional NMCP homologues by whole-genome duplication events. Immunofluorescence microscopy demonstrated localization of a basal NMCP from the moss Physcomitrella patens at the nuclear envelope, suggesting a functional conservation for basal and more evolved NMCPs.

Evolution of TOR-SnRK dynamics in green plants and its integration with phytohormone signaling networks

The target of rapamycin (TOR)-sucrose non-fermenting 1 (SNF1)-related protein kinase 1 (SnRK1) signaling is an ancient regulatory mechanism that originated in eukaryotes to regulate nutrient-dependent growth. Although the TOR-SnRK1 signaling cascade shows highly conserved functions among eukaryotes, studies in the past two decades have identified many important plant-specific innovations in this pathway. Plants also possess SnRK2 and SnRK3 kinases, which originated from the ancient SnRK1-related kinases and have specialized roles in controlling growth, stress responses and nutrient homeostasis in plants. Recently, an integrative picture has started to emerge in which different SnRKs and TOR kinase are highly interconnected to control nutrient and stress responses of plants. Further, these kinases are intimately involved with phytohormone signaling networks that originated at different stages of plant evolution. In this review, we highlight the evolution and divergence of TOR-SnRK signaling components in plants and their communication with each other as well as phytohormone signaling to fine-tune growth and stress responses in plants.

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.

Evolution by duplication: paleopolyploidy events in plants reconstructed by deciphering the evolutionary history of VOZ transcription factors

BACKGROUND

Facilitated by the rapid progress of sequencing technology, comparative genomic studies in plants have unveiled recurrent whole genome duplication (i.e. polyploidization) events throughout plant evolution. The evolutionary past of plant genes should be analyzed in a background of recurrent polyploidy events in distinctive plant lineages. The Vascular Plant One Zinc-finger (VOZ) gene family encode transcription factors associated with a number of important traits including control of flowering time and photoperiodic pathways, but the evolutionary trajectory of this gene family remains uncharacterized.

RESULTS

In this study, we deciphered the evolutionary history of the VOZ gene family by analyses of 107 VOZ genes in 46 plant genomes using integrated methods: phylogenic reconstruction, Ks-based age estimation and genomic synteny comparisons. By scrutinizing the VOZ gene family phylogeny the core eudicot γ event was well circumscribed, and relics of the precommelinid τ duplication event were detected by incorporating genes from oil palm and banana. The more recent T and ρ polyploidy events, closely coincident with the species diversification in Solanaceae and Poaceae, respectively, were also identified. Other important polyploidy events captured included the "salicoid" event in poplar and willow, the "early legume" and "soybean specific" events in soybean, as well as the recent polyploidy event in Physcomitrella patens. Although a small transcription factor gene family, the evolutionary history of VOZ genes provided an outstanding record of polyploidy events in plants. The evolutionary past of VOZ gene family demonstrated a close correlation with critical plant polyploidy events which generated species diversification and provided answer to Darwin's "abominable mystery".

CONCLUSIONS

We deciphered the evolutionary history of VOZ transcription factor family in plants and ancestral polyploidy events in plants were recapitulated simultaneously. This analysis allowed for the generation of an idealized plant gene tree demonstrating distinctive retention and fractionation patterns following polyploidy events.

Whole-Genome Duplication and Plant Macroevolution

Whole-genome duplication (WGD) is characteristic of almost all fundamental lineages of land plants. Unfortunately, the timings of WGD events are loosely constrained and hypotheses of evolutionary consequence are poorly formulated, making them difficult to test. Using examples from across the plant kingdom, we show that estimates of timing can be improved through the application of molecular clock methodology to multigene datasets. Further, we show that phenotypic change can be quantified in morphospaces and that relative phenotypic disparity can be compared in the light of WGD. Together, these approaches facilitate tests of hypotheses on the role of WGD in plant evolution, underscoring the potential of plants as a model system for investigating the role WGD in macroevolution.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

Constraining the timing of whole genome duplication in plant evolutionary history

Whole genome duplication (WGD) has occurred in many lineages within the tree of life and is invariably invoked as causal to evolutionary innovation, increased diversity, and extinction resistance. Testing such hypotheses is problematic, not least since the timing of WGD events has proven hard to constrain. Here we show that WGD events can be dated through molecular clock analysis of concatenated gene families, calibrated using fossil evidence for the ages of species divergences that bracket WGD events. We apply this approach to dating the two major genome duplication events shared by all seed plants (ζ) and flowering plants (ɛ), estimating the seed plant WGD event at 399-381 Ma, and the angiosperm WGD event at 319-297 Ma. These events thus took place early in the stem of both lineages, precluding hypotheses of WGD conferring extinction resistance, driving dramatic increases in innovation and diversity, but corroborating and qualifying the more permissive hypothesis of a 'lag-time' in realizing the effects of WGD in plant evolution.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

The Sphagnome Project: enabling ecological and evolutionary insights through a genus-level sequencing project

Considerable progress has been made in ecological and evolutionary genetics with studies demonstrating how genes underlying plant and microbial traits can influence adaptation and even 'extend' to influence community structure and ecosystem level processes. Progress in this area is limited to model systems with deep genetic and genomic resources that often have negligible ecological impact or interest. Thus, important linkages between genetic adaptations and their consequences at organismal and ecological scales are often lacking. Here we introduce the Sphagnome Project, which incorporates genomics into a long-running history of Sphagnum research that has documented unparalleled contributions to peatland ecology, carbon sequestration, biogeochemistry, microbiome research, niche construction, and ecosystem engineering. The Sphagnome Project encompasses a genus-level sequencing effort that represents a new type of model system driven not only by genetic tractability, but by ecologically relevant questions and hypotheses.

Co-expression and Transcriptome Analysis of Marchantia polymorpha Transcription Factors Supports Class C ARFs as Independent Actors of an Ancient Auxin Regulatory Module

We performed differential gene expression (DGE) and co-expression analyses with genes encoding components of hormonal signaling pathways and the ∼400 annotated transcription factors (TFs) of M. polymorpha across multiple developmental stages of the life cycle. We identify a putative auxin-related co-expression module that has significant overlap with transcripts induced in auxin-treated tissues. Consistent with phylogenetic and functional studies, the class C ARF, MpARF3, is not part of the auxin-related co-expression module and instead is associated with transcripts enriched in gamete-producing gametangiophores. We analyze the Mparf3 and MpmiR160 mutant transcriptomes in the context of coexpression to suggest that MpARF3 may antagonize the reproductive transition via activating the MpMIR11671 and MpMIR529c precursors whose encoded microRNAs target SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcripts of MpSPL1 and MpSPL2. Both MpSPL genes are part of the MpARF3 co-expression group corroborating their functional significance. We provide evidence of the independence of MpARF3 from the auxin-signaling module and provide new testable hypotheses on the role of auxin-related genes in patterning meristems and differentiation events in liverworts.

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.

Transcriptome-Wide Identification and Expression Analysis of DIVARICATA- and RADIALIS-Like Genes of the Mediterranean Orchid Orchis italica

Bilateral symmetry of flowers is a relevant novelty that has occurred many times throughout the evolution of flowering plants. In Antirrhinum majus, establishment of flower dorso-ventral asymmetry is mainly due to interaction of TCP (CYC and DICH) and MYB (DIV, RAD, and DRIF) transcription factors. In the present study, we characterized 8 DIV-, 4 RAD-, and 2 DRIF-like genes from the transcriptome of Orchis italica, an orchid species with bilaterally symmetric and resupinate flowers. We found a similar number of DIV- and RAD-like genes within the genomes of Phalaenopsis equestris and Dendrobium catenatum orchids. Orchid DIV- and RAD-like proteins share conserved motifs whose distribution reflects their phylogeny and analysis of the genomic organization revealed a single intron containing many traces of transposable elements. Evolutionary analysis has shown that purifying selection acts on the DIV- and RAD-like coding regions in orchids, with relaxation of selective constraints in a branch of the DIV-like genes. Analysis of the expression patterns of DIV- and RAD-like genes in O. italica revealed possible redundant functions for some of them. In the perianth of O. italica, the ortholog of DIV and DRIF of A. majus are expressed in all tissues, whereas RAD is mainly expressed in the outer tepals and lip. These data allow for proposal of an evolutionary conserved model in which the expression of the orthologs of the DIV, RAD, and DRIF genes might be related to establishment of flower bilateral symmetry in the nonmodel orchid species O. italica.