Search

Background and Aims: The size and shape of reproductive structures is especially relevant in evolution because these characters are directly related to the capacity for pollination and seed dispersal, a process that plays a basic role in evolutionary patterns. The evolutionary trajectories of reproductive phenotypes in gymnosperms have received special attention in terms of pollination and innovations related to the emergence of the spermatophytes. However, variability of reproductive structures, evolutionary trends and the role of environment in the evolution of cycad species have not been well documented and explored. This study considered this topic under an explicitly phylogenetic and evolutionary approach that included a broad sampling of reproductive structures in the genus Ceratozamia., Methods: We sampled 1400 individuals of 36 Ceratozamia species to explore the evolutionary pattern and identify and evaluate factors that potentially drove their evolution. We analysed characters for both pollen and ovulate strobili within a phylogenetic framework using different methods and characters (i.e. molecular and both quantitative and qualitative morphological) to infer phylogenetic relationships. Using this phylogenetic framework, evolutionary models of trait evolution for strobilar size were evaluated. In addition, quantitative morphological variation and its relation to environmental variables across species were analysed., Key Results: We found contrasting phylogenetic signals between characters of pollen and ovulate strobili. These structures exhibited high morphological disparity in several characters related to size. Results of analyses of evolutionary trajectories suggested a stabilizing selection model. With regard to phenotype-environment, the analysis produced mixed results and differences for groups in the vegetation type where the species occur; however, a positive relationship with climatic variables was found., Conclusions: The integrated approach synthesized reproductive phenotypic variation with current phylogenetic hypotheses and provided explicit statements of character evolution. The characters of volume for ovulate strobili were the most informative, and could provide a reference for further study of the evolutionary complexity in Ceratozamia. Finally, heterogeneous environments, which are under changing weather conditions, promote variability of reproductive structures., (© The Author(s) 2024. Published by Oxford University Press on behalf of the Annals of Botany Company.)

Plants synthesize natural products via lineage-specific offshoots of their core metabolic pathways, including fatty acid synthesis. Recent studies have shed light on new fatty acid-derived natural products and their biosynthetic pathways in disparate plant species. Inspired by this progress, we set out to develop tools for exploring the evolution of fatty-acid derived products. We sampled multiple species from all major clades of euphyllophytes, including ferns, gymnosperms, and angiosperms (monocots and eudicots), and we show that the compositional profiles (though not necessarily the total amounts) of fatty-acid derived surface waxes from preserved plant specimens are consistent with those obtained from freshly collected tissue in a semi-quantitative and sometimes quantitative manner. We then sampled herbarium specimens representing 57 monocot species to assess the phylogenetic distribution and evolution, of two fatty acid-derived natural products found in that clade: beta-diketones and alkylresorcinols. These chemical data, combined with analyses of 26 monocot genomes, revealed a co-occurrence (though not necessarily a causal relationship) between whole genome duplication and the evolution of diketone synthases from an ancestral alkylresorcinol synthase-like polyketide synthase. Limitations of using herbarium specimen wax profiles as proxies for those of fresh tissue seem likely to include effects from loss of epicuticular wax crystals, effects from preservation techniques, and variation in wax chemical profiles due to genotype or environment. Nevertheless, this work reinforces the widespread utility of herbarium specimens for studying leaf surface waxes (and possibly other chemical classes) and reveals some of the evolutionary history of fatty acid-derived natural products within monocots., (© 2024 The Author(s). The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: Paulownia, an ecologically and economically valuable plant species native to China, is notable for its excellent timber quality and strong adaptability. Among them, Paulownia catalpifolia displays the ability to survive in cold climate, a trait associated with northern China. Yet, the molecular information for its cold-tolerance has not been explored. This study was to investigate the changes in physiological indices and transcript levels of P. catalpifolia following cold exposure, which could provide evidence for revealing whether there were differences in the genetic basis of inducing physiological perturbations between moderate low temperature (MLT) and extreme low temperature (ELT)., Results: The detection of physiological indices under diverse degrees of chilling stress showed similar patterns of alteration. Enhanced accumulation of osmoregulatory substances, such as soluble sugar and soluble protein, were more conducive under ELT compared to MLT in P. catalpifolia. Moreover, we observed leaf wilting symptoms distinctly after exposure to ELT for 48 h, while this effect was not obvious after MLT exposure for 48 h. Comparative transcriptomic analysis between MLT and ELT demonstrated 13,688 differentially expressed genes (DEGs), most of them appeared after 12 h and 48 h of treatment. GO and KEGG analyses elucidated prominent enrichment in aromatic-L-amino-acid decarboxylase activity term and carbohydrate metabolism pathways. Therefore, it was speculated that the DEGs involved in the above processes might be related to the difference in the contents of soluble protein and soluble sugar between MLT and ELT. Time series clustering analyses further highlighted several key genes engaged in the 'Glycosyltransferases', 'Galactose metabolism' and 'Starch and sucrose metabolism' pathways as well as the 'tyrosine decarboxylase activity' term. For instance, cellulose synthase-like A (CLSA2/9), raffinose synthase (RafS2), β-amylase (BAM1) and tyrosine/DOPA decarboxylase (TYDC1/2/5) genes, diverging in their expression trends between MLT and ELT, might significantly affect the soluble sugar and soluble protein abundance within P. catalpifolia., Conclusion: Between MLT and ELT treatments, partial overlaps in response pathways of P. catalpifolia were identified, while several genes regulating the accumulation of osmotic adjustment substances had disparate expression patterns. These findings could provide a novel physiological and molecular perspective for P. catalpifolia to adapt to complex low temperature habitats., (© 2024. The Author(s).)

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies., (© 2024 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: Gymnospermium kiangnanense is the only species distributed in the subtropical region within the spring ephemeral genus Gymnospermium . Extensive human exploitation and habitat destruction have resulted in a rapid shrink of G. kiangnanense populations. This study utilizes microsatellite markers to analyze the genetic diversity and structure and to deduce historical population events of extant populations of G. kiangnanense ., Methods: A total of 143 individuals from eight extant populations of G. kiangnanense , including two populations from Anhui Province and six populations from Zhejiang Province, were analyzed with using 21 pairs of microsatellite markers. Genetic diversity indices were calculated using Cervus, GENEPOP, GenALEX. Population structure was assessed using genetic distance (UPGMA), principal coordinate analysis (PCoA), Bayesian clustering method (STRUCTURE), and molecular variation analysis of variance (AMOVA). Population history events were inferred using DIYABC., Results: The studied populations of G. kiangnanense exhibited a low level of genetic diversity ( He = 0.179, I = 0.286), but a high degree of genetic differentiation ( F ST = 0.521). The mean value of gene flow ( N m ) among populations was 1.082, indicating prevalent gene exchange via pollen dispersal. Phylogeographic analyses suggested that the populations of G. kiangnanense were divided into two lineages, Zhejiang (ZJ) and Anhui (AH). These two lineages were separated by the Huangshan-Tianmu Mountain Range. AMOVA analysis revealed that 36.59% of total genetic variation occurred between the two groups. The ZJ lineage was further divided into the Hangzhou (ZJH) and Zhuji (ZJZ) lineages, separated by the Longmen Mountain and Fuchun River. DIYABC analyses suggested that the ZJ and AH lineages were separated at 5.592 ka, likely due to the impact of Holocene climate change and human activities. Subsequently, the ZJZ lineage diverged from the ZJH lineage around 2.112 ka. Given the limited distribution of G. kiangnanense and the significant genetic differentiation among its lineages, both in-situ and ex-situ conservation strategies should be implemented to protect the germplasm resources of G. kiangnanense ., Competing Interests: The authors declare that they have no competing interests., (© 2024 Liu et al.)

Main Conclusion: Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Gymnosperms are long-lived, cone-bearing seed plants that include some of the most ancient extant plant species. These relict land plants have evolved to survive in habitats marked by chronic or episodic stress. Their ability to thrive in these environments is partly due to their phenotypic flexibility, and epigenetic regulation likely plays a crucial part in this plasticity. We review the current knowledge on abiotic and biotic stress memory in gymnosperms and the possible epigenetic mechanisms underlying long-term phenotypic adaptations. We also discuss recent technological improvements and new experimental possibilities that likely will advance our understanding of epigenetic regulation in these ancient and hard-to-study plants., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Gymnosperms are mostly dioecious, and their staminate strobili undergo a longer developmental period than those of angiosperms. However, the underlying molecular mechanisms remain unclear. This study aimed to identify key genes and pathways involved in staminate strobilus development and dehiscence in Torreya grandis. Through weighted gene co-expression network analysis (WGCNA), we identified fast elongation-related genes enriched in carbon metabolism and auxin signal transduction, whereas dehiscence-related genes were abundant in alpha-linolenic acid metabolism and the phenylpropanoid pathway. Based on WGCNA, we also identified PHYTOCHROME-INTERACTING FACTOR4 (TgPIF4) as a potential regulator for fast elongation of staminate strobilus and 2 WRKY proteins (TgWRKY3 and TgWRKY31) as potential regulators for staminate strobilus dehiscence. Multiple protein-DNA interaction analyses showed that TgPIF4 directly activates the expression of TRANSPORT INHIBITOR RESPONSE2 (TgTIR2) and NADP-MALIC ENZYME (TgNADP-ME). Overexpression of TgPIF4 significantly promoted staminate strobilus elongation by elevating auxin signal transduction and pyruvate content. TgWRKY3 and TgWRKY31 bind to the promoters of the lignin biosynthesis gene PHENYLALANINE AMMONIA-LYASE (TgPAL) and jasmonic acid metabolism gene JASMONATE O-METHYLTRANSFERASE (TgJMT), respectively, and directly activate their transcription. Overexpression of TgWRKY3 and TgWRKY31 in the staminate strobilus led to early dehiscence, accompanied by increased lignin and methyl jasmonate levels, respectively. Collectively, our findings offer a perspective for understanding the growth of staminate strobili in gymnosperms., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus., (© 2024. The Author(s).)

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba , Taxus chinensis , and Pinus tabuliformis . In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Mutations and subsequent repair processes are known to be strongly context-dependent in the flowering-plant chloroplast genome. At least six flanking bases, three on each side, can have an influence on the relative rates of different types of mutation at any given site. In this analysis, examine context and substitution at noncoding and fourfold degenerate coding sites in gymnosperm DNA. The sequences are analyzed in sets of three, allowing the inference of the substitution direction and the generation of context-dependent rate matrices. The size of the dataset limits the analysis to the tetranucleotide context of the sites, but the evidence shows that there are significant contextual effects, with patterns that are similar to those observed in angiosperms. These effects most likely represent an influence on the underlying mutation/repair dynamics. The data extend the plastome lineages that feature very complex patterns of mutation, which can have significant effects on the evolutionary dynamics of the chloroplast genome.

Forests, comprising 31% of the Earth's surface, play pivotal roles in regulating the carbon, water, and energy cycles. Despite being far less diverse than angiosperms, gymnosperms account for over 50% of the global woody biomass production. To sustain growth and development, gymnosperms have evolved the capacity to sense and respond to cyclical environmental signals, such as changes in photoperiod and seasonal temperature, which initiate growth (spring and summer) and dormancy (fall and winter). Cambium, the lateral meristem responsible for wood formation, is reactivated through a complex interplay among hormonal, genetic, and epigenetic factors. Temperature signals perceived in early spring induce the synthesis of several phytohormones, including auxins, cytokinins, and gibberellins, which in turn reactivate cambium cells. Additionally, microRNA-mediated genetic and epigenetic pathways modulate cambial function. As a result, the cambium becomes active during the summer, resulting in active secondary xylem (i.e., wood) production, and starts to become inactive in autumn. This review summarizes and discusses recent findings regarding the climatic, hormonal, genetic, and epigenetic regulation of wood formation in gymnosperm trees (i.e., conifers) in response to seasonal changes.

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C 18 Δ 9 -elongase and a C 20 Δ 5 -desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ 5 -desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants., (© 2023. The Author(s).)

During the course of evolution, organisms have developed genetic mechanisms in response to various environmental stresses including wounding from mechanical damage or herbivory-caused injury. A previous study of wounding response in the plant tobacco identified a unique wound-induced gene, aptly named KED due to its coding for a protein that has an unusually high content of amino acids lysine (K), glutamic acid (E) and aspartic acid (D). However, by far little is known about this intriguing gene. In this study, we investigated the evolutionary aspects of the KED-rich coding genes. We found that a consistent pattern of wound-induced KED gene expression is maintained across representative species of angiosperm and gymnosperm. KED genes can be identified in species from all groups of land plants (Embryophyta). All the KED proteins from vascular plants (Tracheophyta) including angiosperm, gymnosperm, fern and lycophyte share a conserved 19-amino acid domain near the C-terminus, whereas bryophytes (moss, liverwort and hornwort) possess KED-rich, multi-direct-repeat sequences that are distinct from the vascular plant KEDs. We detected KED-rich sequences in Charophyta species but not in Chlorophyta wherever genome sequences are available. Our studies suggest diverse and complex evolution pathways for land plant KED genes. Vascular plant KEDs exhibit high evolutionary conservation, implicating their shared function in response to wounding stress. The extraordinary enrichment of amino acids K, E and D in these groups of distinct and widely distributed proteins may reflect the structural and functional requirement for these three residues during some 600 million years of land plant evolution., Competing Interests: All authors have declared that no competing interests exist., (Copyright: © 2023 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: The complex physical structure and abundant repeat sequences make it difficult to assemble the mitogenomes of seed plants, especially gymnosperms. Only approximately 33 mitogenomes of gymnosperms have been reported. However, as the most widely distributed and the second largest family among gymnosperms, Cupressaceae has only six assembled mitogenomes, including five draft mitogenomes and one complete mitogenome, which has greatly hindered the understanding of mitogenome evolution within this large family, even gymnosperms., Results: In this study, we assembled and validated the complete mitogenome of Thuja sutchuenensis, with a size of 2.4 Mb. Multiple sequence units constituted its complex structure, which can be reduced to three linear contigs and one small circular contig. The analysis of repeat sequences indicated that the numbers of simple sequence repeats increased during the evolutionary history of gymnosperms, and the mitogenome of Thuja sutchuenensis harboured abundant extra-long repeats (more than 5 kb). Additionally, the longest repeat sequence identified in these seven gymnosperms also came from the mitogenome of Thuja sutchuenensis, with a length of up to 47 kb. The analysis of colinear blocks and gene clusters both revealed that the orders of mitochondrial genes within gymnosperms was not conserved. The comparative analysis showed that only four tRNAs were shared by seven gymnosperms, namely, trnD-GUC, trnE-UUC, trnI-CAU and trnY-GUA. Furthermore, four genes have undergone potential positive selection in most gymnosperm species, namely, atp8, ccmB, mttB and sdh4., Conclusion: We successfully assembled the second complete mitogenome within Cupressaceae and verified that it consisted of multiple sequence units. Our study also indicated that abundant long repeats may contribute to the generation of the complex conformation of the mitogenome of Thuja sutchuenensis. The investigation of Thuja sutchuenensis's mitogenome in our study provides new insight into further understanding the complex mitogenome architecture within gymnosperms., (© 2023. The Author(s).)

The angiosperms, or flowering plants, arose at least 135 million years ago (Ma) and rapidly diversified to form over 300,000 species alive today. This group appears, however, to have separated from its closest living relatives, the extant gymnosperms, much earlier: over 300 Ma. Representatives of basally-diverging angiosperm lineages are of key importance to studies aimed at reconstructing the most recent common ancestor of living angiosperms, including its morphological, anatomical, eco-physiological and molecular aspects. Furthermore, evo-devo comparisons of angiosperms with living gymnosperms may help to determine how the many novel aspects of angiosperms, including those of the flower, first came about. This chapter reviews literature on the origin of angiosperms and focusses on basally-diverging angiosperms and gymnosperms that show advantages as potential experimental models, reviewing information and protocols for the use of these species in an evo-devo context. The final section suggests a means by which data from living and fossil groups could be integrated to better elucidate evolutionary events that took place on the long stem-lineage that apparently preceded the radiation of living angiosperms., (© 2023. Springer Science+Business Media, LLC, part of Springer Nature.)

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads., (© 2022 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported., (© 2022 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5-15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants., (© The Author(s) 2022. Published by Oxford University Press GigaScience.)

Premise: Anatomically preserved evidence for a novel clade of gymnosperms emphasizes diversity of seed plants immediately prior to the appearance of angiosperm fossils in the paleontological record., Methods: Cupulate seeds from the Early Cretaceous Apple Bay locality (Vancouver Island) are described from serial cellulose acetate peels and three-dimensional reconstruction. Phylogenetic context is assessed through the comparative analysis of gymnosperm seed producing fructifications and maximum parsimony analysis of a revised morphological data set for seed plant phylogeny., Results: Xadzigacalix quatsinoensis gen. et sp. nov. is characterized by an orthotropous ovule with an elongated micropyle and complex integument, enclosed within a radial cupule. The micropylar canal is elongated; and the nucellus extends into the micropyle to seal the post pollination ovule. Except at the apex of the micropyle, the seed is completely enclosed by a parenchymatous cupule with ca. 20 axially elongated secretory ducts. The cupulate seed is produced upon a triangular woody stele, consisting of a parenchymatous pith surrounded by radially aligned tracheids. The stele produces three short terete traces that terminate within the base of the cupule as transfusion tissue at the seed chalaza., Conclusions: Organography, vascularization, nature of the integument and nucellus, and configuration of the micropylar canal distinguish Xadzigacalix quatsinoensis from all other gymnosperm clades. Cladistic analyses suggest the new plant may have affinities with gnetophytes or angiosperms. These results are complemented with a critical re-evaluation of ovulate structures for Mesozoic gymnosperms, providing new insight into plant diversity immediately antecedent to the explosive diversification of flowering plants., (© 2022 The Authors. American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

Glutamine synthetase (GS) is a key enzyme responsible for the incorporation of inorganic nitrogen in the form of ammonium into the amino acid glutamine. In plants, two groups of functional GS enzymes are found: eubacterial GSIIb (GLN2) and eukaryotic GSIIe (GLN1/GS). Only GLN1/GS genes are found in vascular plants, which suggests that they are involved in the final adaptation of plants to terrestrial life. The present phylogenetic study reclassifies the different GS genes of seed plants into three clusters: GS1a, GS1b and GS2. The presence of genes encoding GS2 has been expanded to Cycadopsida gymnosperms, which suggests the origin of this gene in a common ancestor of Cycadopsida, Ginkgoopsida and angiosperms. GS1a genes have been identified in all gymnosperms, basal angiosperms and some Magnoliidae species. Previous studies in conifers and the gene expression profiles obtained in ginkgo and magnolia in the present work could explain the absence of GS1a in more recent angiosperm species (e.g. monocots and eudicots) as a result of the redundant roles of GS1a and GS2 in photosynthetic cells. Altogether, the results provide a better understanding of the evolution of plant GS isoenzymes and their physiological roles, which is valuable for improving crop nitrogen use efficiency and productivity. This new view of GS evolution in plants, including a new cytosolic GS group (GS1a), has important functional implications in the context of plant metabolism adaptation to global changes., (© 2022 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists., (© 2022. The Author(s).)

Background: Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections., Results: All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny., Conclusions: Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection., (© 2022. The Author(s).)

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution., (© 2021 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Background: Soil salinization is causing ecosystem degradation and crop yield reduction worldwide, and elucidation of the mechanism of salt-tolerant plants to improve crop yield is highly significant. Podocarpus macrophyllus is an ancient gymnosperm species with a unique environmental adaptation strategy that may be attributed to its lengthy evolutionary process. The present study investigated the physiological and molecular responses of P. macrophyllus plants to salt stress by analyzing its photosynthetic system and antioxidant enzyme activity. We also analyzed the differentially expressed genes (DEGs) in P. macrophyllus under salt stress using RNA sequencing and de novo transcriptome assembly., Results: Salt treatment significantly affected the photosynthetic system in P. macrophyllus seedlings, which decreased chlorophyll content, altered chloroplast ultrastructure, and reduced photosynthesis. The activities of antioxidant enzymes increased significantly following salt stress treatment. Transcriptome analysis showed that salt stress induced a large number of genes involved in multiple metabolic and biological regulation processes. The transcription levels of genes that mediate phytohormone transport or signaling were altered. K + and Ca 2+ transporter-encoding genes and the MYB transcription factor were upregulated under salt stress. However, the genes involved in cell wall biosynthesis and secondary metabolism were downregulated., Conclusion: Our research identified some important pathways and putative genes involved in salt tolerance in P. macrophyllus and provided clues for elucidating the mechanism of salt tolerance and the utilization of the salt tolerance genes of P. macrophyllus for crop improvement., (© 2021. The Author(s).)

Background: Chloroplast transfer RNAs (tRNAs) can participate in various vital processes. Gymnosperms have important ecological and economic value, and they are the dominant species in forest ecosystems in the Northern Hemisphere. However, the evolution and structural changes in chloroplast tRNAs in gymnosperms remain largely unclear., Results: In this study, we determined the nucleotide evolution, phylogenetic relationships, and structural variations in 1779 chloroplast tRNAs in gymnosperms. The numbers and types of tRNA genes present in the chloroplast genomes of different gymnosperms did not differ greatly, where the average number of tRNAs was 33 and the frequencies of occurrence for various types of tRNAs were generally consistent. Nearly half of the anticodons were absent. Molecular sequence variation analysis identified the conserved secondary structures of tRNAs. About a quarter of the tRNA genes were found to contain precoded 3' CCA tails. A few tRNAs have undergone novel structural changes that are closely related to their minimum free energy, and these structural changes affect the stability of the tRNAs. Phylogenetic analysis showed that tRNAs have evolved from multiple common ancestors. The transition rate was higher than the transversion rate in gymnosperm chloroplast tRNAs. More loss events than duplication events have occurred in gymnosperm chloroplast tRNAs during their evolutionary process., Conclusions: These findings provide novel insights into the molecular evolution and biological characteristics of chloroplast tRNAs in gymnosperms., (© 2021. The Author(s).)

Morphological complexity is a notable feature of multicellular life, although whether it evolves gradually or in early bursts is unclear. Vascular plant reproductive structures, such as flowers, are familiar examples of complex morphology. In this study, we use a simple approach based on the number of part types to analyze changes in complexity over time. We find that reproductive complexity increased in two pulses separated by ~250 million years of stasis, including an initial rise in the Devonian with the radiation of vascular plants and a pronounced increase in the Late Cretaceous that reflects flowering plant diversification. These pulses are associated with innovations that increased functional diversity, suggesting that shifts in complexity are linked to changes in function regardless of whether they occur early or late in the history of vascular plants.

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Background: Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes., Results: Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants., Conclusions: This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants., (© 2021. The Author(s).)

The gymnosperm Welwitschia mirabilis belongs to the ancient, enigmatic gnetophyte lineage. It is a unique desert plant with extreme longevity and two ever-elongating leaves. We present a chromosome-level assembly of its genome (6.8 Gb/1 C) together with methylome and transcriptome data to explore its astonishing biology. We also present a refined, high-quality assembly of Gnetum montanum to enhance our understanding of gnetophyte genome evolution. The Welwitschia genome has been shaped by a lineage-specific ancient, whole genome duplication (~86 million years ago) and more recently (1-2 million years) by bursts of retrotransposon activity. High levels of cytosine methylation (particularly at CHH motifs) are associated with retrotransposons, whilst long-term deamination has resulted in an exceptionally GC-poor genome. Changes in copy number and/or expression of gene families and transcription factors (e.g. R2R3MYB, SAUR) controlling cell growth, differentiation and metabolism underpin the plant's longevity and tolerance to temperature, nutrient and water stress., (© 2021. The Author(s).)

Gymnosperms are a unique lineage of plants that currently lack a high-quality reference genome due to their large genome size and high repetitive sequence content. Here, we report a nearly complete genome assembly for Ginkgo biloba with a genome size of 9.87 Gb, an N50 contig size of 1.58 Mb and an N50 scaffold size of 775 Mb. We were able to accurately annotate 27,832 protein-coding genes in total, superseding the inaccurate annotation of 41,840 genes in a previous draft genome assembly. We found that expansion of the G. biloba genome, accompanied by the notable extension of introns, was mainly caused by the insertion of long terminal repeats rather than the recent occurrence of whole-genome duplication events, in contrast to the findings of a previous report. We also identified candidate genes in the central pair, intraflagellar transport and dynein protein families that are associated with the formation of the spermatophore flagellum, which has been lost in all seed plants except ginkgo and cycads. The newly obtained Ginkgo genome provides new insights into the evolution of the gymnosperm genome.

Established model systems in the flowering plants have greatly advanced our understanding of plant developmental biology, facilitating in turn its investigation across diverse land plants. The reliance on a limited number of model organisms, however, constitutes a barrier for future progress in evolutionary developmental biology (evo-devo). In particular, a more thorough understanding of seed plant character evolution and of its genetic and developmental basis has been hampered in part by a lack of gymnosperm model systems, since most are trees with decades-long generation times. Guided by the premise that future model organisms should be selected based on their character diversity, rather than simply phylogenetic "position," we highlight biological questions of potential interest that can be addressed via comparative studies in Ephedra (Gnetales). In addition to having relatively small genomes and shorter generation times in comparison to most other gymnosperms, Ephedra are amenable to investigations on the evolution of the key reproductive seed plant innovations of pollination and seed dispersal, as well as on polyploidy, and adaptation to extreme environments., (© 2021 Wiley Periodicals LLC.)

Gymnosperm taxa are incredibly diverse in many traits, including taxa with evergreen or deciduous leaves; broad, compound, needle-like or scale-like leaves; trees, shrubs, and lianas; and taxa with seed cones that range from fleshy to woody. Although less appreciated, xylem conduits are also diverse in structure among extant gymnosperm taxa. Within the xylem of gymnosperms, axial transport occurs predominantly via tracheids, although 10-40% of gymnosperm taxa, particularly within the Gnetophyta and Cycadophyta, also contain vessels. Gymnosperm taxa vary greatly in their inter-conduit pit structure, with different types of pit membranes and pitting arrangements. While some taxa display torus-margo bordered pits (60%), many others do not contain tori (40%), and at least some taxa without tori occur within each of the four extant phyla (Coniferophyta, Cycadophyta, Ginkgophyta, and Gnetophyta). Pit membrane types vary within families but appear relatively conserved within genera. Woody species with torus-bearing pit membranes occur in colder environments (lower mean annual temperature) compared to those without tori; but occurrence does not differ with mean annual precipitation. Detailed descriptions of pit membrane types are lacking for many species and genera, indicating a need for increased anatomical study. Increased knowledge of these traits could provide a unique experimental context in which to study the evolution of conduit networks, the hydraulic implications of conduit and pit structure, and the diverse structural and functional strategies utilized by gymnosperms. There are myriad potential study questions and research opportunities within this unique and diverse group of plants., (© 2021 Botanical Society of America.)

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Competing Interests: The author declares no competing interest.

Background: Gnetum is an economically important tropical and subtropical gymnosperm genus with various dietary, industrial and medicinal uses. Many carbohydrates, proteins and fibers accumulate during the ripening of Gnetum seeds. However, the molecular mechanisms related to this process remain unknown., Results: We therefore assembled a full-length transcriptome from immature and mature G. luofuense seeds using PacBio sequencing reads. We identified a total of 5726 novel genes, 9061 alternative splicing events, 3551 lncRNAs, 2160 transcription factors, and we found that 8512 genes possessed at least one poly(A) site. In addition, gene expression comparisons of six transcriptomes generated by Illumina sequencing showed that 14,323 genes were differentially expressed from an immature stage to a mature stage with 7891 genes upregulated and 6432 genes downregulated. The expression of 14 differentially expressed transcription factors from the MADS-box, Aux/IAA and bHLH families was validated by qRT-PCR, suggesting that they may have important roles in seed ripening of G. luofuense., Conclusions: These findings provide a valuable molecular resource for understanding seed development of gymnosperms.

Given the 2,400-fold range of genome sizes (0.06-148.9 Gbp (gigabase pair)) of seed plants (angiosperms and gymnosperms) with a broadly similar gene content (amounting to approximately 0.03 Gbp), the repeat-sequence content of the genome might be expected to increase with genome size, resulting in the largest genomes consisting almost entirely of repetitive sequences. Here we test this prediction, using the same bioinformatic approach for 101 species to ensure consistency in what constitutes a repeat. We reveal a fundamental change in repeat turnover in genomes above around 10 Gbp, such that species with the largest genomes are only about 55% repetitive. Given that genome size influences many plant traits, habits and life strategies, this fundamental shift in repeat dynamics is likely to affect the evolutionary trajectory of species lineages.

The term rheophyte describes a biological group of flood-tolerant plants that are confined to the beds of swift-running streams and rivers in nature and grow up to flood level, but not beyond the reach of regularly occurring flash floods. Although over 35 yr have passed since the first global census of rheophytes, no updates have been recorded regarding the number of taxa in this biological group in seed plants. Therefore, the present work aimed to (1) review the main topics associated with rheophytism (e.g., morphological characteristics, genetic studies, geographic distribution, conservation, and evolutionary aspects); (2) provide an updated checklist of rheophytes distributed around the world considering the two main groups in seed plants (gymnosperms and angiosperms); (3) demonstrate the distribution of rheophytism in the angiosperm phylogeny; and (4) estimate the geographical distribution and richness of selected taxa on the world map for the first time. All data compiled for the present study originated from a search of peer-reviewed articles, secondary literature (theses, dissertations, reports, books, and floras), and electronic facilities. We compiled a data set composed of four taxa in gymnosperms (Podocarpaceae family) and 1,368 taxa (including obligate, facultative, and unclassified rheophytes) distributed in 114 families and 508 genera in angiosperms. Most of the studied taxa belong to eudicotyledons (72.81%), while 1.46% belong to magnoliids, and 25.73% belong to monocotyledons. The families with the highest number of taxa in descending order are Podostemaceae, Araceae, Myrtaceae, Rubiaceae, Asteraceae, Apocynaceae, Arecaceae, Fabaceae, Phyllanthaceae, and Poaceae. Of the 114 families plotted in angiosperm phylogeny, at least 80 harbor obligate rheophytes. The geographical distribution of rheophytes in angiosperms, as expected based on the first census of this biological group, is mainly in the tropical and subtropical regions. The high richness of rheophytic taxa was mostly found in southern Mexico, southern China, Borneo, and northern and eastern Australia. In contrast, the geographical distribution of rheophytes in gymnosperms is restricted to New Caledonia and Tasmania. The present study will help to advance knowledge regarding the diversity of rheophytes in angiosperms and gymnosperms while drawing attention to this biological group, which has often been overlooked. There are no copyright restrictions. Please cite this data paper when the data are used in publications and teaching events., (© 2020 by the Ecological Society of America.)

Hundreds of plant mitogenomes have been sequenced from angiosperms, but relatively few mitogenomes are available from its sister lineage, gymnosperms. To examine mitogenomic diversity among extant gymnosperms, we generated draft mitogenomes from 11 diverse species and compared them with four previously published mitogenomes. Examined mitogenomes from Pinaceae and cycads retained all 41 protein genes and 26 introns present in the common ancestor of seed plants, whereas gnetophyte and cupressophyte mitogenomes experienced extensive gene and intron loss. In Pinaceae and cupressophyte mitogenomes, an unprecedented number of exons are distantly dispersed, requiring trans-splicing of 50-70% of mitochondrial introns to generate mature transcripts. RNAseq data confirm trans-splicing of these dispersed exons in Pinus. The prevalence of trans-splicing in vascular plant lineages with recombinogenic mitogenomes suggests that genomic rearrangement is the primary cause of shifts from cis- to trans-splicing in plant mitochondria., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

DNA barcoding is a valuable tool to support species identification with broad applications from traditional taxonomy, ecology, forensics, food analysis, and environmental science. We introduce Microfluidic Enrichment Barcoding (MEBarcoding) for plant DNA Barcoding, a cost-effective method for high-throughput DNA barcoding. MEBarcoding uses the Fluidigm Access Array to simultaneously amplify targeted regions for 48 DNA samples and hundreds of PCR primer pairs (producing up to 23,040 PCR products) during a single thermal cycling protocol. As a proof of concept, we developed a microfluidic PCR workflow using the Fluidigm Access Array and Illumina MiSeq. We tested 96 samples for each of the four primary DNA barcode loci in plants: rbcL, matK, trnH-psbA, and ITS. This workflow was used to build a reference library for 78 families and 96 genera from all major plant lineages - many currently lacking in public databases. Our results show that this technique is an efficient alternative to traditional PCR and Sanger sequencing to generate large amounts of plant DNA barcodes and build more comprehensive barcode databases.

Seasonal phenomena in plants are primarily affected by day length and temperature. The shoot transcriptomes of trees grown in the field and a controlled-environment chamber were compared to characterize genes that control annual rhythms and the effects of day length- and temperature-regulated genes in the gymnosperm Japanese cedar (Cryptomeria japonica D. Don), which exhibits seasonally indeterminate growth. Annual transcriptome dynamics were clearly demonstrated by principal component analysis using microarray data obtained under field-grown conditions. Analysis of microarray data from trees grown in a controlled chamber identified 2,314 targets exhibiting significantly different expression patterns under short-day (SD) and long-day conditions, and 2,045 targets exhibited significantly different expression patterns at 15°C (LT; low temperature) versus 25°C. Interestingly, although growth was suppressed under both SD and LT conditions, approximately 80% of the SD- and LT-regulated targets differed, suggesting that each factor plays a unique role in the annual cycle. The top 1,000 up-regulated targets in the growth/dormant period in the field coincided with more than 50% of the SD- and LT-regulated targets, and gene co-expression network analysis of the annual transcriptome indicated a close relationship between the SD- and LT-regulated targets. These results indicate that the respective effects of day length and temperature interact to control annual transcriptome dynamics. Well-known upstream genes of signaling pathways responsive to environmental conditions, such as the core clock (LHY/CjLHYb and CCA1/CjLHYa) and PEBP family (MFT) genes, exhibited unique expression patterns in Japanese cedar compared with previous reports in other species, suggesting that these genes control differences in seasonal regulation mechanisms between species. The results of this study provide new insights into seasonal regulation of transcription in Japanese cedar., Competing Interests: The authors have declared that no competing interests exist.

Angiosperms that are biotically pollinated typically produce flowers with bright and contrasting colours that help to attract pollinators and hence contribute to the reproductive success of the species. This colourful array contrasts with the much less multicoloured reproductive structures of the four living gymnosperm lineages, which are mostly wind pollinated, though cycads and Gnetales are predominantly pollinated by insects that feed on surface fluids from the pollination drops. This review examines the possible evolutionary pathways and cryptic clues for flower colour in both living and fossil seed plants. It investigates how the ancestral flowering plants could have overcome the inevitable trade-off that exists between attracting pollinators and minimizing herbivory, and explores the possible evolutionary and biological inferences from the colours that occur in some living gymnosperms. The red colours present in the seed-cone bracts of some living conifers result from accumulation of anthocyanin pigments; their likely primary function is to help protect the growing plant tissues under particular environmental conditions. Thus, the visual cue provided by colour in flower petals could have first evolved as a secondary effect, probably post-dating the evolution of bee colour vision but occurring before the subsequent functional accumulation of a range of different flower pigments., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Experimental Biology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms., Results: Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species., Conclusions: We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Plant mitogenomes can be difficult to assemble because they are structurally dynamic and prone to intergenomic DNA transfers, leading to the unusual situation where an organelle genome is far outnumbered by its nuclear counterparts. As a result, comparative mitogenome studies are in their infancy and some key aspects of genome evolution are still known mainly from pregenomic, qualitative methods. To help address these limitations, we combined machine learning and in silico enrichment of mitochondrial-like long reads to assemble the bacterial-sized mitogenome of Norway spruce (Pinaceae: Picea abies). We conducted comparative analyses of repeat abundance, intergenomic transfers, substitution and rearrangement rates, and estimated repeat-by-repeat homologous recombination rates. Prompted by our discovery of highly recombinogenic small repeats in P. abies, we assessed the genomic support for the prevailing hypothesis that intramolecular recombination is predominantly driven by repeat length, with larger repeats facilitating DNA exchange more readily. Overall, we found mixed support for this view: Recombination dynamics were heterogeneous across vascular plants and highly active small repeats (ca. 200 bp) were present in about one-third of studied mitogenomes. As in previous studies, we did not observe any robust relationships among commonly studied genome attributes, but we identify variation in recombination rates as a underinvestigated source of plant mitogenome diversity., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Caffeoyl shikimate esterase (CSE) has been reported to be involved in lignin biosynthesis; however, studies of CSE in gymnosperms are lacking. In this study, CSE was successfully cloned from Larix kaempferi ( LkCSE ) based on Larix laricina transcriptome screening. LkCSE was likely to have catalytic activity based on homologous sequence alignment and phylogenetic analyses of CSEs from different species. In vitro assays with the recombinant enzyme validated the catalytic activity of LkCSE, indicating its function in converting caffeoyl shikimate into caffeate and shikimate. Additionally, the optimum reaction pH and temperature of LkCSE were determined to be 6.0 and 30 °C, respectively. The values of K m and V max of CSE for caffeoyl shikimate were 98.11 μM and 14.44 nM min -1 , respectively. Moreover, LkCSE was observed to have tissue expression specificity and was abundantly expressed in stems and leaves, especially stems, which was 50 times higher than the expression levels of roots. Lastly, translational fusion assays using LkCSE fused with green fluorescent proteins (GFP) in tobacco leaves indicated that LkCSE was localized in the plasma membrane and endoplasmic reticulum (ER). These results revealed that CSE clearly functions in gymnosperms and it is possible for LkCSE to interact with other ER-resident proteins and regulate mass flux in the monolignol biosynthesis pathway.

The adaptation of herbivorous insects to new host plants is key to their evolutionary success in diverse environments. Many insects are associated with mutualistic gut bacteria that contribute to the host's nutrition and can thereby facilitate dietary switching in polyphagous insects. However, how gut microbial communities differ between populations of the same species that feed on different host plants remains poorly understood. Most species of Pyrrhocoridae (Hemiptera: Heteroptera) are specialist seed-feeders on plants in the family Malvaceae, although populations of one species, Probergrothius angolensis, have switched to the very distantly related Welwitschia mirabilis plant in the Namib Desert. We first compared the development and survival of laboratory populations of Pr. angolensis with two other pyrrhocorids on seeds of Welwitschia and found only Pr. angolensis was capable of successfully completing its development. We then collected Pr. angolensis in Namibia from Malvaceae and Welwitschia host plants, respectively, to assess their bacterial and fungal community profiles using high-throughput amplicon sequencing. Comparison with long-term laboratory-reared insects indicated stable associations of Pr. angolensis with core bacteria (Commensalibacter, Enterococcus, Bartonella and Klebsiella), but not with fungi or yeasts. Phylogenetic analyses of core bacteria revealed relationships to other insect-associated bacteria, but also found new taxa indicating potential host-specialized nutritional roles. Importantly, the microbial community profiles of bugs feeding on Welwitschia versus Malvaceae revealed stark and consistent differences in the relative abundance of core bacterial taxa that correlate with the host-plant switch; we were able to reproduce this result through feeding experiments. Thus, a dynamic gut microbiota may provide a means for insect adaptation to new host plants in new environments when food plants are extremely divergent., (© 2019 John Wiley & Sons Ltd.)

Taiwania cryptomerioides is a monotypic gymnosperm species, valued for the high decay resistance of its wood. This durability has been attributed to the abundance of terpenoids, especially the major diterpenoid metabolite ferruginol, with antifungal and antitermite activity. Specialized diterpenoid metabolism in gymnosperms primarily recruits bifunctional class-I/II diterpene synthases (diTPSs), whereas monofunctional class-II and class-I enzymes operate in angiosperms. In this study, we identified a previously unrecognized group of monofunctional diTPSs in T. cryptomerioides, which suggests a distinct evolutionary divergence of the diTPS family in this species. Specifically, five monofunctional diTPS functions not previously observed in gymnosperms were characterized, including monofunctional class-II enzymes forming labda-13-en-8-ol diphosphate (LPP, TcCPS2) and (+)-copalyl diphosphate (CPP, TcCPS4), and three class-I diTPSs producing biformene (TcKSL1), levopimaradiene (TcKSL3) and phyllocladanol (TcKSL5), respectively. Methyl jasmonate (MeJA) elicited the accumulation of levopimaradiene and the corresponding biosynthetic diTPS genes, TcCPS4 and TcKSL3, is consistent with a possible role in plant defense. Furthermore, TcCPS4 and TcKSL3 are likely to contribute to abietatriene biosynthesis via levopimaradiene as an intermediate in ferruginol biosynthesis in Taiwania. In conclusion, this study provides deeper insight into the functional landscape and molecular evolution of specialized diterpenoid metabolism in gymnosperms as a basis to better understand the role of these metabolites in tree chemical defense., (© 2019 The Authors The Plant Journal © 2019 John Wiley & Sons Ltd.)

Species assemble into communities through ecological and evolutionary processes. Phylogenetic niche conservatism-the tendency of species to retain ancestral ecological distributions-is thought to influence which species from a regional species pool can persist in a particular environment. We analyzed data for seed plants in China to test hypotheses about the distribution of species within regional floras. Of 16 environmental variables, actual evapotranspiration, minimum temperature of the coldest month, and annual precipitation most strongly influenced regional species richness, phylogenetic dispersion, and phylogenetic diversity for both gymnosperms (cone-bearing plants) and angiosperms (flowering plants). For most evolutionary clades at, and above, the family level, the relationships between metrics of phylogenetic dispersion (i.e., average phylogenetic distance among species), or phylogenetic diversity, and the 3 environmental variables were consistent with the tropical niche conservatism hypothesis, which predicts closer phylogenetic relatedness and reduced phylogenetic diversity with increasing environmental stress. The slopes of the relationships between phylogenetic relatedness and the 3 environmental drivers identified in this analysis were steeper for primarily tropical clades, implying greater niche conservatism, than for primarily temperate clades. These observations suggest that the distributions of seed plants across large-scale environmental gradients in China are constrained by conserved adaptations to the physical environment, i.e., phylogenetic niche conservatism., Competing Interests: The authors declare no competing interest.