Liverwort oil bodies: diversity, biochemistry, and molecular cell biology of the earliest secretory structure of land plants.

Liverworts are known for their large chemical diversity. Much of this diversity is synthesized and enclosed within oil bodies (OBs), a synapomorphy of the lineage. OBs contain the enzymes to biosynthesize and store large quantities of sesquiterpenoids and other compounds while limiting their cytotoxicity. Recent important biochemical and molecular discoveries related to OB formation, diversity, and biochemistry allow comparison with other secretory structures of land plants from an evo-devo perspective. This review addresses and discusses the most recent advances in OB origin, development, and function towards understanding the importance of these organelles in liverwort physiology and adaptation to changing environments. Our mapping of OB types and chemical compounds to the current liverwort phylogeny suggests that OBs were present in the most recent common ancestor of liverworts, supporting that OBs evolved as the first secretory structures in land plants. Yet, we require better sampling to define the macroevolutionary pattern governing the ancestral type of OB. We conclude that current efforts to find molecular mechanisms responsible for the morphological and chemical diversity of secretory structures will help understand the evolution of each major group of land plants, and open new avenues in biochemical research on bioactive compounds in bryophytes and vascular plants.

A fossil record of land plant origins from charophyte algae.

Molecular time trees indicating that embryophytes originated around 500 million years ago (Ma) during the Cambrian are at odds with the record of fossil plants, which first appear in the mid-Silurian almost 80 million years later. This time gap has been attributed to a missing fossil plant record, but that attribution belies the case for fossil spores. Here, we describe a Tremadocian (Early Ordovician, about 480 Ma) assemblage with elements of both Cambrian and younger embryophyte spores that provides a new level of evolutionary continuity between embryophytes and their algal ancestors. This finding suggests that the molecular phylogenetic signal retains a latent evolutionary history of the acquisition of the embryophytic developmental genome, a history that perhaps began during Ediacaran-Cambrian time but was not completed until the mid-Silurian (about 430 Ma).

Genome-wide analysis of the AREB/ABF gene lineage in land plants and functional analysis of TaABF3 in Arabidopsis.

BACKGROUND: Previous studies have shown that ABFs (abscisic acid-responsive transcription factors) are important ABA-signaling components that participate in abiotic stress response. However, little is known about the function of ABFs in Triticum aestivum. In addition, although various ABFs have been identified in other species, the phylogenetic relationship between ABF transcription factors has not been systemically investigated in land plants. RESULTS: In this study, we systemically collected ABFs from land plants and analyzed the phylogenetic relationship of these ABF genes. The ABF genes are present in all the land plants we investigated, including moss, lycophyte, monocots, and eudicots. Furthermore, these ABF genes are phylogenetically divided into seven subgroups, differentiations that are supported by variation in the gene structure, protein properties, and motif patterns. We further demonstrated that the expression of ABF genes varies among different tissues and developmental stages, and are induced by one or more environmental stresses. Furthermore, we found that three wheat ABFs (TaABF1, TaABF2, and TaABF3) were significantly induced by drought stress. Compared with wild-type (WT) plants, transgenic Arabidopsis plants overexpressing TaABF3 displayed enhanced drought tolerance. CONCLUSIONS: These results provide important ground work for understanding the phylogenetic relationships between plant ABF genes. Our results also indicate that TaABFs may participate in regulating plant response to abiotic stresses.

Three-dimensional growth: a developmental innovation that facilitated plant terrestrialization.

One of the most transformative events in the history of life on earth was the transition of plants from water to land approximately 470 million years ago. Within the Charophyte green algae, the closest living relatives of land plants, body plans have evolved from those that comprise simple unicells to those that are morphologically complex, large and multicellular. The Charophytes developed these broad ranging body plans by exploiting a range of one-dimensional and two-dimensional growth strategies to produce filaments, mats and branches. When plants were confronted with harsh conditions on land, they were required to make significant changes to the way they shaped their body plans. One of the fundamental developmental transitions that occurred was the evolution of three-dimensional growth and the acquisition of apical cells with three or more cutting faces. Plants subsequently developed a range of morphological adaptations (e.g. vasculature, roots, flowers, seeds) that enabled them to colonise progressively drier environments. 3D apical growth also evolved convergently in the brown algae, completely independently of the green lineage. This review summarises the evolving developmental complexities observed in the early divergent Charophytes all the way through to the earliest conquerors of land, and investigates 3D apical growth in the brown algae.

The Expansion and Diversification of Pentatricopeptide Repeat RNA-Editing Factors in Plants.

The RNA-binding pentatricopeptide repeat (PPR) family comprises hundreds to thousands of genes in most plants, but only a few dozen in algae, indicating massive gene expansions during land plant evolution. The nature and timing of these expansions has not been well defined due to the sparse sequence data available from early-diverging land plant lineages. In this study, we exploit the comprehensive OneKP datasets of over 1000 transcriptomes from diverse plants and algae toward establishing a clear picture of the evolution of this massive gene family, focusing on the proteins typically associated with RNA editing, which show the most spectacular variation in numbers and domain composition across the plant kingdom. We characterize over 2 250 000 PPR motifs in over 400 000 proteins. In lycophytes, polypod ferns, and hornworts, nearly 10% of expressed protein-coding genes encode putative PPR editing factors, whereas they are absent from algae and complex-thalloid liverworts. We show that rather than a single expansion, most land plant lineages with high numbers of editing factors have continued to generate novel sequence diversity. We identify sequence variations that imply functional differences between PPR proteins in seed plants versus non-seed plants and variations we propose to be linked to seed-plant-specific editing co-factors. Finally, using the sequence variations across the datasets, we develop a structural model of the catalytic DYW domain associated with C-to-U editing and identify a clade of unique DYW variants that are strong candidates as U-to-C RNA-editing factors, given their phylogenetic distribution and sequence characteristics.

Mitochondrial genomes of the early land plant lineage liverworts (Marchantiophyta): conserved genome structure, and ongoing low frequency recombination.

BACKGROUND: In contrast to the highly labile mitochondrial (mt) genomes of vascular plants, the architecture and composition of mt genomes within the main lineages of bryophytes appear stable and invariant. The available mt genomes of 18 liverwort accessions representing nine genera and five orders are syntenous except for Gymnomitrion concinnatum whose genome is characterized by two rearrangements. Here, we expanded the number of assembled liverwort mt genomes to 47, broadening the sampling to 31 genera and 10 orders spanning much of the phylogenetic breadth of liverworts to further test whether the evolution of the liverwort mitogenome is overall static. RESULTS: Liverwort mt genomes range in size from 147 Kb in Jungermanniales (clade B) to 185 Kb in Marchantiopsida, mainly due to the size variation of intergenic spacers and number of introns. All newly assembled liverwort mt genomes hold a conserved set of genes, but vary considerably in their intron content. The loss of introns in liverwort mt genomes might be explained by localized retroprocessing events. Liverwort mt genomes are strictly syntenous in genome structure with no structural variant detected in our newly assembled mt genomes. However, by screening the paired-end reads, we do find rare cases of recombination, which means multiple concurrent genome structures may exist in the vegetative tissues of liverworts. Our phylogenetic analyses of the nuclear encoded double stand break repair protein families revealed liverwort-specific subfamilies expansions. CONCLUSIONS: The low repeat recombination level, selection, along with the intensified nuclear surveillance, might together shape the structural evolution of liverwort mt genomes.

The commonness of rarity: Global and future distribution of rarity across land plants.

A key feature of life's diversity is that some species are common but many more are rare. Nonetheless, at global scales, we do not know what fraction of biodiversity consists of rare species. Here, we present the largest compilation of global plant diversity to quantify the fraction of Earth's plant biodiversity that are rare. A large fraction, ~36.5% of Earth's ~435,000 plant species, are exceedingly rare. Sampling biases and prominent models, such as neutral theory and the k-niche model, cannot account for the observed prevalence of rarity. Our results indicate that (i) climatically more stable regions have harbored rare species and hence a large fraction of Earth's plant species via reduced extinction risk but that (ii) climate change and human land use are now disproportionately impacting rare species. Estimates of global species abundance distributions have important implications for risk assessments and conservation planning in this era of rapid global change.

Long-Reads Reveal That the Chloroplast Genome Exists in Two Distinct Versions in Most Plants.

The chloroplast genome usually has a quadripartite structure consisting of a large single copy region and a small single copy region separated by two long inverted repeats. It has been known for some time that a single cell may contain at least two structural haplotypes of this structure, which differ in the relative orientation of the single copy regions. However, the methods required to detect and measure the abundance of the structural haplotypes are labor-intensive, and this phenomenon remains understudied. Here, we develop a new method, Cp-hap, to detect all possible structural haplotypes of chloroplast genomes of quadripartite structure using long-read sequencing data. We use this method to conduct a systematic analysis and quantification of chloroplast structural haplotypes in 61 land plant species across 19 orders of Angiosperms, Gymnosperms, and Pteridophytes. Our results show that there are two chloroplast structural haplotypes which occur with equal frequency in most land plant individuals. Nevertheless, species whose chloroplast genomes lack inverted repeats or have short inverted repeats have just a single structural haplotype. We also show that the relative abundance of the two structural haplotypes remains constant across multiple samples from a single individual plant, suggesting that the process which maintains equal frequency of the two haplotypes operates rapidly, consistent with the hypothesis that flip-flop recombination mediates chloroplast structural heteroplasmy. Our results suggest that previous claims of differences in chloroplast genome structure between species may need to be revisited.

Strigolactone synthesis is ancestral in land plants, but canonical strigolactone signalling is a flowering plant innovation.

BACKGROUND: Strigolactones (SLs) are an important class of carotenoid-derived signalling molecule in plants, which function both as exogenous signals in the rhizosphere and as endogenous plant hormones. In flowering plants, SLs are synthesized by a core pathway of four enzymes and are perceived by the DWARF14 (D14) receptor, leading to degradation of SMAX1-LIKE7 (SMXL7) target proteins in a manner dependent on the SCFMAX2 ubiquitin ligase. The evolutionary history of SLs is poorly understood, and it is not clear whether SL synthesis and signalling are present in all land plant lineages, nor when these traits evolved. RESULTS: We have utilized recently-generated genomic and transcriptomic sequences from across the land plant clade to resolve the origin of each known component of SL synthesis and signalling. We show that all enzymes in the core SL synthesis pathway originated at or before the base of land plants, consistent with the previously observed distribution of SLs themselves in land plant lineages. We also show that the late-acting enzyme LATERAL BRANCHING OXIDOREDUCTASE (LBO) may be considerably more ancient than previously thought. We perform a detailed phylogenetic analysis of SMXL proteins and show that specific SL target proteins only arose in flowering plants. We also assess diversity and protein structure in the SMXL family, identifying several previously unknown clades. CONCLUSIONS: Overall, our results suggest that SL synthesis is much more ancient than canonical SL signalling, consistent with the idea that SLs first evolved as rhizosphere signals and were only recruited much later as hormonal signals.

Recent gene duplications dominate evolutionary dynamics of adaptor protein complex subunits in embryophytes.

Adaptor protein complexes and the related complexes COPI and TSET function in packaging vesicles for transport among endomembrane compartments in eukaryotic cells. Differences in the complement of these complexes in lineages such as yeast and mammals as well as apicomplexan and kinetoplastid parasites via loss or duplication of subunits appears to reflect specialization in their respective trafficking systems. The model plant Arabidopsis thaliana possesses multiple paralogues for adaptor protein complex subunits, raising questions as to the timing and extent of these duplications in embryophytes (land plants). However, adaptor protein complex evolution in embryophytes is unexplored. Therefore, we analyzed genomes of diverse embryophytes and closely related green algae using extensive homology searches and phylogenetic analysis of 35 complex subunit proteins. The results reveal numerous paralogues, the vast majority of which, approximately 97%, arose from recent duplication events. This suggests that specialization of these protein complexes may occur frequently but independently in embryophytes.

Evolutionary and Comparative Expression Analyses of TCP Transcription Factor Gene Family in Land Plants.

The plant-specific Teosinte-branched 1/Cycloidea/Proliferating (TCP) transcription factor genes are involved in plants' development, hormonal pathways, and stress response but their evolutionary history is uncertain. The genome-wide analysis performed here for 47 plant species revealed 535 TCP candidates in terrestrial plants and none in aquatic plants, and that TCP family genes originated early in the history of land plants. Phylogenetic analysis divided the candidate genes into Classes I and II, and Class II was further divided into CYCLOIDEA (CYC) and CINCINNATA (CIN) clades; CYC is more recent and originated from CIN in angiosperms. Protein architecture, intron pattern, and sequence characteristics were conserved in each class or clade supporting this classification. The two classes significantly expanded through whole-genome duplication during evolution. Expression analysis revealed the conserved expression of TCP genes from lower to higher plants. The expression patterns of Class I and CIN genes in different stages of the same tissue revealed their function in plant development and their opposite effects in the same biological process. Interaction network analysis showed that TCP proteins tend to form protein complexes, and their interaction networks were conserved during evolution. These results contribute to further functional studies on TCP family genes.

Potential of Transcript Editing Across Mitogenomes of Early Land Plants Shows Novel and Familiar Trends.

RNA editing alters the identity of nucleotides in an RNA sequence so that the mature transcript differs from the template defined in the genome. This process has been observed in chloroplasts and mitochondria of both seed and early land plants. However, the frequency of RNA editing in plant mitochondria ranges from zero to thousands of editing sites. To date, analyses of RNA editing in mitochondria of early land plants have been conducted on a small number of genes or mitochondrial genomes of a single species. This study provides an overview of the mitogenomic RNA editing potential of the main lineages of these two groups of early land plants by predicting the RNA editing sites of 33 mitochondrial genes of 37 species of liverworts and mosses. For the purpose of the research, we newly assembled seven mitochondrial genomes of liverworts. The total number of liverwort genera with known complete mitogenome sequences has doubled and, as a result, the available complete mitogenome sequences now span almost all orders of liverworts. The RNA editing site predictions revealed that C-to-U RNA editing in liverworts and mosses is group-specific. This is especially evident in the case of liverwort lineages. The average level of C-to-U RNA editing appears to be over three times higher in liverworts than in mosses, while the C-to-U editing frequency of the majority of genes seems to be consistent for each gene across bryophytes.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants.

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

No mass extinction for land plants at the Permian-Triassic transition.

The most severe mass extinction among animals took place in the latest Permian (ca. 252 million years ago). Due to scarce and impoverished fossil floras from the earliest Triassic, the common perception has been that land plants likewise suffered a mass extinction, but doubts remained. Here we use global occurrence data of both plant macro- and microfossils to analyse plant biodiversity development across the Permian-Triassic boundary. We show that the plant fossil record is strongly biased and that evidence for a mass extinction among plants in the latest Permian is not robust. The taxonomic diversities of gymnosperm macrofossils and of the pollen produced by this group are particularly incongruent. Our results indicate that gymnosperm macrofossils are considerably undersampled for the Early Triassic, which creates the impression of increased gymnosperm extinction in the latest Permian.

Mechanistic basis for the evolution of chalcone synthase catalytic cysteine reactivity in land plants.

Flavonoids are important polyphenolic natural products, ubiquitous in land plants, that play diverse functions in plants' survival in their ecological niches, including UV protection, pigmentation for attracting pollinators, symbiotic nitrogen fixation, and defense against herbivores. Chalcone synthase (CHS) catalyzes the first committed step in plant flavonoid biosynthesis and is highly conserved in all land plants. In several previously reported crystal structures of CHSs from flowering plants, the catalytic cysteine is oxidized to sulfinic acid, indicating enhanced nucleophilicity in this residue associated with its increased susceptibility to oxidation. In this study, we report a set of new crystal structures of CHSs representing all five major lineages of land plants (bryophytes, lycophytes, monilophytes, gymnosperms, and angiosperms), spanning 500 million years of evolution. We reveal that the structures of CHS from a lycophyte and a moss species preserve the catalytic cysteine in a reduced state, in contrast to the cysteine sulfinic acid seen in all euphyllophyte CHS structures. In vivo complementation, in vitro biochemical and mutagenesis analyses, and molecular dynamics simulations identified a set of residues that differ between basal-plant and euphyllophyte CHSs and modulate catalytic cysteine reactivity. We propose that the CHS active-site environment has evolved in euphyllophytes to further enhance the nucleophilicity of the catalytic cysteine since the divergence of euphyllophytes from other vascular plant lineages 400 million years ago. These changes in CHS could have contributed to the diversification of flavonoid biosynthesis in euphyllophytes, which in turn contributed to their dominance in terrestrial ecosystems.

Community phylogeny of the globally critically imperiled pine rockland ecosystem.

PREMISE OF THE STUDY: Community phylogenetic methods incorporate information on evolutionary relationships into studies of organismal assemblages. We used a community phylogenetic framework to investigate relationships and biogeographic affinities and to calculate phylogenetic signal of endemism and invasiveness for the flora of the pine rocklands-a globally critically imperiled ecosystem with a significant portion of its distribution in South Florida, United States. METHODS: We reconstructed phylogenetic relationships of 538 vascular plant taxa, which represent 92.28% of the vascular flora of the pine rocklands. We estimated phylogenetic signal for endemism and invasiveness using phylogenetic generalized linear mixed models. We determined the native range for each species in the data set and calculated the total number of species sourced from each region and all possible combinations of these regions. KEY RESULTS: The pine rockland flora includes representatives of all major vascular plant lineages, and most species have native ranges in the New World. There was strong phylogenetic signal for endemism, but not for invasiveness. CONCLUSIONS: Community phylogenetics has high potential value for conservation planning, particularly for fragmented and endangered ecosystems like the pine rockland. Strong phylogenetic signal for endemic species in our data set, which also tend to be threatened or endangered, can help to identify species at risk, as well as fragments where those species occur, highlighting conservation priorities. Our results indicate, at least in the pine rockland ecosystem, no phylogenetic signal for invasive species, and thus other information must be used to predict the potential for invasiveness.

Genome-wide identification and analysis of MICU genes in land plants and their potential role in calcium stress.

Mitochondrial calcium uptake (MICU) plays a vital role in the regulation of mitochondrial calcium homeostasis, and, consequently, influences calcium signaling transduction. Although genes involved in mitochondrial calcium uptake have been well studied in animals, less is known about their ubiquity and function in plants. In this study, we identified 96 MICU genes in land plants. On the basis of phylogenetic analysis of MICU proteins, they were classified into three clades: MICU from eudicots (Clade I), from monocots (Clade II), and from a basal angiosperm, a bryophyte, and a lycophyte (Clade III). Pairwise identity analysis across all MICU proteins showed that they are highly conserved among land plants at the protein level. Conserved motif analysis showed that most MICU proteins contained three EF-hands, and an additional EF-hand motif first identified in the MICU of Arabidopsis thaliana but not mammals was found in all 96 putative MICU proteins. This suggests that a cellular pathway of calcium uptake and signaling that requires three EF-hand motifs is evolutionarily conserved in plants. In addition, we discovered that MICU-defective mutants of Arabidopsis thaliana exhibited longer roots than wild-type under high calcium stress. Concurrently, the mRNA transcription levels of MICU were decreased under high calcium conditions. These results suggest that loss-of-function mutations of MICU may have potential roles in helping plants resist high calcium stress. This study provides clues to the possible role of plant MICU in mitochondrial calcium uptake, as well as useful information to support further studies on MICU function in plants.

Discovery of the World's Smallest Terrestrial Pteridophyte.

Ophioglossum L. commonly known as "adder's tongue fern", has been of great interest due to the highest number of chromosomes in any organism so far known in biological world. Here, a new species of adder's tongue fern has been discovered and reported from Western Ghats of India. It is prominently distinct from the other known taxa in Ophioglossaceae family. Phylogenetic analysis of three chloroplast DNA (cpDNA) regions (trnL-F, rbcL and psbA-trnH) unambiguously designate this adder's tongue fern as the distinct lineage and is sister to the clade containing O. parvifolium and O. nudicaule. Azolla caroliniana - an aquatic fern (average size, 0.5-1.5 cm), is the smallest fern on the earth. Our discovery discloses a new species of adder's tongue fern and ranking it among the smallest terrestrial fern in the world, attaining an average size of only 1-1.2 cm.

Contrasting spatial, temporal and environmental patterns in observation and specimen based species occurrence data.

Species occurrence data records the location and time of an encounter with a species, and is valuable for many aspects of ecological and evolutionary analyses. A key distinction within species occurrence data is between (1) collected and preserved specimens that can be taxonomically validated (i.e., natural history collections), and (2) observations, which are more error prone but richer in terms of number and spread of observations. In this study we analyse the distribution in temporal, spatial, taxonomic and environmental coverage of specimen- and observation based species occurrence data for land plants in Norway, a region with strong climatic and human population density gradients. Of 4.8 million species occurrence records, the majority (78%) were observations. However, there was a greater species richness in the specimen record (N = 4691) than in the observation record (N = 3193) and most species were recorded more as specimens than observations. Specimen data was on average older, and collected later during the year. Both record types were highly influenced by a small number of prolific contributors. The species most highly represented in the observation data set were widespread or invasive, while in the specimen records, taxonomically challenging species were overrepresented. Species occurrence records were unevenly spatially distributed. Both specimen and observation records were concentrated in regions of Norway with high human population density and with high temperatures and precipitation, but in different regions within Norway. Observation and specimen records thus differ in taxonomic, temporal, spatial and environmental coverage for a well-sampled group and study region, potentially influencing the ecological inferences made from studies utilizing species occurrence data. The distribution of observation data dominates the dataset, so inferences of species diversity and distributions do not correspond to the evolutionary or physiological knowledge of species, which is based on specimen data. We make recommendations for users of biodiversity data, and collectors to better exploit the complementary strengths of these distinct biodiversity data types.

Accounting for Calibration Uncertainty: Bayesian Molecular Dating as a "Doubly Intractable" Problem.

This study introduces a new Bayesian technique for molecular dating that explicitly accommodates for uncertainty in the phylogenetic position of calibrated nodes derived from the analysis of fossil data. The proposed approach thus defines an adequate framework for incorporating expert knowledge and/or prior information about the way fossils were collected in the inference of node ages. Although it belongs to the class of "node-dating" approaches, this method shares interesting properties with "tip-dating" techniques. Yet, it alleviates some of the computational and modeling difficulties that hamper tip-dating approaches. The influence of fossil data on the probabilistic distribution of trees is the crux of the matter considered here. More specifically, among all the phylogenies that a tree model (e.g., the birth-death process) generates, only a fraction of them "agree" with the fossil data. Bayesian inference under the new model requires taking this fraction into account. However, evaluating this quantity is difficult in practice. A generic solution to this issue is presented here. The proposed approach relies on a recent statistical technique, the so-called exchange algorithm, dedicated to drawing samples from "doubly intractable" distributions. A small example illustrates the problem of interest and the impact of uncertainty in the placement of calibration constraints in the phylogeny given fossil data. An analysis of land plant sequences and multiple fossils further highlights the pertinence of the proposed approach.