Horsetails and ferns are a monophyletic group and the closest living relatives to seed plants

Patterns in Genome-Wide Codon Usage Bias in Representative Species of Lycophytes and Ferns

The latest research shows that ferns and lycophytes have distinct evolutionary lineages. The codon usage patterns of lycophytes and ferns have not yet been documented. To investigate the gene expression profiles across various plant lineages with respect to codon usage, analyze the disparities and determinants of gene evolution in primitive plant species, and identify appropriate exogenous gene expression platforms, the whole-genome sequences of four distinct species were retrieved from the NCBI database. The findings indicated that Ceratopteris richardii, Adiantum capillus-veneris, and Selaginella moellendorffii exhibited an elevated A/U content in their codon base composition and a tendency to end with A/U. Additionally, S. capillus-veneris had more C/G in its codons and a tendency to end with C/G. The ENC values derived from both ENC-plot and ENC-ratio analyses deviated significantly from the standard curves, suggesting that the codon usage preferences of these four species were primarily influenced by genetic mutations and natural selection, with natural selection exerting a more prominent influence. This finding was further supported by PR2-Plot, neutrality plot analysis, and COA. A combination of RSCU and ENC values was used as a reference criterion to rank the codons and further identify the optimal codons. The study identified 24 high-frequency codons in C. richardii, A. capillus-veneris, and Diphasiastrum complanatum, with no shared optimal codons among the four species. Arabidopsis thaliana and Ginkgo biloba exhibited similar codon preferences to the three species, except for S. moellendorffii. This research offers a theoretical framework at the genomic codon level for investigating the phylogenetic relationships between lycophytes and ferns, shedding light on gene codon optimization and its implications for genetic engineering in breeding.

Transcriptomics and metabolomics analyses provide insights into resistance genes of tree ferns

As ancient organisms, tree ferns play a crucial role as an evolutionary bridge between lower and higher plant species, providing various utilitarian benefits. However, they face challenges such as overexploitation, climate change, adverse environmental conditions, and insect pests, resulting in conservation concerns. In this study, we provide an overview of metabolic and transcriptomic resources of leaves in two typical tree ferns, A. spinulosa and A. metteniana, and explore the resistance genes for the first time. The landscape of metabolome showed that the compound skimmin may hold medicinal significance. A total of 111 differentially accumulated metabolites (DAMs) were detected, with pathway enrichment analysis highlighting 14 significantly enriched pathways, including 2-oxocarboxylic acid metabolism possibly associated with environmental adaptations. A total of 14,639 differentially expressed genes (DEGs) were found, among which 606 were resistance (R) genes. We identified BAM1 as a significantly differentially expressed R gene, which is one of the core genes within the R gene interaction network. Both the maximum-likelihood phylogenetic tree and the PPI network revealed a close relationship between BAM1, FLS2, and TMK. Moreover, BAM1 showed a significant positive correlation with neochlorogenic acid and kaempferol-7-O-glucoside. These metabolites, known for their antioxidant and anti-inflammatory properties, likely play a crucial role in the defense response of tree ferns. This research provides valuable insights into the metabolic and transcriptomic differences between A. spinulosa and A. metteniana, enhancing our understanding of resistance genes in tree ferns.

Calamitean Cones and Their In Situ Spores from the Pennsylvanian Limnic Basins of the Czech Republic

This paper describes the in situ spores of the Calamospora type, macerated from sixty-one specimens of calamitean cones belonging to sixteen species of genera, such as the Palaeostachya, Macrostachya, Calamostachys, and Huttonia from the Pennsylvanian Czech Republic period, specifically from the Moscovian/Kasimovian ages (i.e., Duckmantian-Stephanian). The in situ spores were compared to twenty dispersed species of Calamospora. The majority of spores were microspores; however, some cones yielded both micro- and megaspores. Morphological variations of the in situ spores, including the diameter, labrum, contact area, ontogenetic stages, and secondary folds of the exine, are described, including their importance for the classification of calamospores. The relationships of Elaterites, Pteroretis, Vestispora, and some monopseudosaccate spores are discussed. All Paleozoic Calamospora-producing parent plants are summarized.

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

Indication-Specific Effect of a Phytotherapeutic Remedy on Human Fetal Osteoblastic Cells: An in vitro Analysis

BACKGROUND

Impaired fracture healing is a recurring interdisciplinary medical challenge. Alternative treatment concepts, apart from conventional medicine, are popular, but scientific evidence on their effects is still lacking. Plant-derived substances are widely assumed to support bone homeostasis. To clarify the effects on bone healing mechanisms, a commercially available, homeopathic-spagyric remedy, containing inter alia two herbal substances with assumed osteogenic potential, equisetum arvense and bellis perennis, was analyzed.

METHODS

Human fetal osteoblastic (hFOB) 1.19 cells were incubated with the test substance in serial dilutions from 10 to 0.00001%. Cell viability has been evaluated through ATP level (CTG assay) and MTT tetrazolium reduction. Cell proliferation was analyzed by BrdU incorporation and cell migration by wound healing assay (WHA) via image analysis. Additionally, determination of the expression of key genes via real-time PCR and proteins via proteome array for inflammation, cell proliferation, and angiogenesis were performed.

RESULTS

An incubation of hFOB 1.19 cells with the test substance for 24/72 h showed no reduction in cell number, viability, or proliferation. Cell migration was unimpaired. The test substance induced inflammatory genes and growth factors along with genes of osseous regeneration (ALP, Col1, IL-1α, IL-6, IL-8, IL-10, Osteocalcin, Osteonectin, RUMX2, TGF, VEGFA). Increased protein expression was found in multiple cytokines, chemokines, and acute phase proteins.

CONCLUSION

The test substance did not impair cell vitality parameters (MTT, CTG, BrdU, and WHA). A tendency to activate growth factors, bone regeneration genes, and proteins was shown for osteoblasts, indicating a possible positive effect on osteogenic processes.

Hintergrund

Störungen des komplexen Prozesses der Knochenheilung stellen auch heutzutage noch eine interdisziplinäre Herausforderung dar. Es existieren zahlreiche alternative Therapiekonzepte, deren Evidenz jedoch häufig nicht belegt ist. Es wird davon ausgegangen, dass pflanzliche Substanzen die Knochenheilung unterstützen können. Wir analysierten die Wirkung eines kommerziellen, homeopathisch-spagyrischen Heilmittels, welches unter anderen zwei Pflanzenstoffe enthält, denen ein osteogenes Potential zugeschrieben wird (Equisetum arvense und Bellis perennis).

Methoden

Es erfolgte eine Inkubation humaner fetaler Osteoblastenzellen (hFOB 1.19) mit der Testsubstanz in absteigender Verdünnung von 10 bis 0.00001%. Die Zellvitalität wurde anhand der Zellzahlbestimmung durch ATP-abhängige metabolische Aktivität mittels CellTiter-Glo® (CTG) Test sowie durch Tetrazolium Reduktion (MTT) evaluiert. Die Zellproliferation wurde durch Inkorporation von Bromdesoxyuridin (BrdU) in die DNA aktiver Zellen analysiert. Der Wound Healing Assay (WHA) diente der Quantifizierung der Zellmigration. Zusätzlich wurde die Expression bestimmter Schlüsselgene mittels real-time PCR und die Proteinexpression via proteom array für Inflammation, Zellproliferation und Angiogenese erhoben.

Ergebnisse

Die Inkubation von hFOB 1.19 mit der Testsubstanz für 24/72 Stunden führte zu keiner Reduktion von Zellzahl, -vitalität oder -proliferation. Auch die Zellmigration war unbeeinträchtigt. Es zeigte sich eine Induktion inflammatorischer Gene, Wachstumsfaktoren sowie Genen der knöchernen Regeneration (ALP, Col1, IL-1α, IL-6, IL-8, IL-10, Osteocalcin, Osteonectin, RUMX2, TGF, VEGFA). Verschiedene Zytokine, Chemokine und Akute Phase Proteine wurden vermehrt exprimiert.

Schlussfolgerung

Die Testsubstanz hatte keine negativen Auswirkungen auf die gemessenen Zellvitalitätsparameter (MTT, CTG, BrdU and WHA). Es zeigte sich eine Aktivierungstendenz für Wachstumsfaktoren, Gene und Proteine der Knochenregeneration, die auf einen möglichen positiven Effekt der Substanz auf den Prozess des Knochenheilung hinweisen.

Embracing uncertainty: The way forward in plant fossil phylogenetics

Although molecular phylogenetics remains the most widely used method of inferring the evolutionary history of living groups, the last decade has seen a renewed interest in morphological phylogenetics, mostly driven by the promises that integrating the fossil record in phylogenetic trees offers to our understanding of macroevolutionary processes and dynamics and the possibility that the inclusion of fossil taxa could lead to more accurate phylogenetic hypotheses. The plant fossil record presents some challenges to its integration in a phylogenetic framework. Phylogenies including plant fossils often retrieve uncertain relationships with low support, or lack of resolution. This low support is due to the pervasiveness of morphological convergence among plant organs and the fragmentary nature of many plant fossils, and it is often perceived as a fundamental weakness reducing the utility of plant fossils in phylogenetics. Here I discuss the importance of uncertainty in morphological phylogenetics and how we can identify important information from different patterns and types of uncertainty. I also review a set of methodologies that can allow us to understand the causes underpinning uncertainty and how these practices can help us to further our knowledge of plant fossils. I also propose that a new visual language, including the use of networks instead of trees, represents an improvement on the old visualization based on consensus trees and more adequately serves phylogeneticists working with plant fossils. This set of methods and visualization tools represents an important way forward in a fundamental field for our understanding of the evolutionary history of plants.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

The genome of homosporous maidenhair fern sheds light on the euphyllophyte evolution and defences

Euphyllophytes encompass almost all extant plants, including two sister clades, ferns and seed plants. Decoding genomes of ferns is the key to deep insight into the origin of euphyllophytes and the evolution of seed plants. Here we report a chromosome-level genome assembly of Adiantum capillus-veneris L., a model homosporous fern. This fern genome comprises 30 pseudochromosomes with a size of 4.8-gigabase and a contig N50 length of 16.22 Mb. Gene co-expression network analysis uncovered that homospore development in ferns has relatively high genetic similarities with that of the pollen in seed plants. Analysing fern defence response expands understanding of evolution and diversity in endogenous bioactive jasmonates in plants. Moreover, comparing fern genomes with those of other land plants reveals changes in gene families important for the evolutionary novelties within the euphyllophyte clade. These results lay a foundation for studies on fern genome evolution and function, as well as the origin and evolution of euphyllophytes.

Pteridofitas: helechos y licófitos de El mundo vegetal de los Andes peruanos en el contexto actual

Los helechos y licófitos mencionados en El mundo Vegetal de los Andes peruanos de Weberbauer brindan una breve información de su riqueza y sustento a las observaciones ambientales para cada categoría de su esquema fitogeográfico del Perú. Este estudio actualiza tanto el concepto y características de los pteridofitos mencionados en el capítulo Unidades Sistemáticas, como la nomenclatura de los 126 nombres de los taxones mencionados en la obra. Además, se ofrece una breve historia del estudio del grupo en el Perú, para proveer el contexto de la inclusión de estas plantas por Weberbauer. Y se mencionan las colecciones que Weberbauer realizara del grupo y, en particular, las citadas como ejemplos. Se ofrece un perfil de las tareas por realizar para completar la labor biogeográfica iniciada por Weberbauer.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

Comparative chloroplast genome and transcriptome analysis on the ancient genus Isoetes from China

Isoetes is a famous living fossil that plays a significant role in the evolutionary studies of the plant kingdom. To explore the adaptive evolution of the ancient genus Isoetes from China, we focused on Isoetes yunguiensis (Q.F. Wang and W.C. Taylor), I. shangrilaensis (X. Li, Y.Q. Huang, X.K. Dai & X. Liu), I. taiwanensis (DeVol), I. sinensis (T.C. Palmer), I. hypsophila_GHC (Handel-Mazzetti), and I. hypsophila_HZS in this study. We sequenced, assembled, and annotated six individuals' chloroplast genomes and transcriptomes, and performed a series of analyses to investigate their chloroplast genome structures, RNA editing events, and adaptive evolution. The six chloroplast genomes of Isoetes exhibited a typical quadripartite structure with conserved genome sequence and structure. Comparative analyses of Isoetes species demonstrated that the gene organization, genome size, and GC contents of the chloroplast genome are highly conserved across the genus. Besides, our positive selection analyses suggested that one positively selected gene was statistically supported in Isoetes chloroplast genomes using the likelihood ratio test (LRT) based on branch-site models. Moreover, we detected positive selection signals using transcriptome data, suggesting that nuclear-encoded genes involved in the adaption of Isoetes species to the extreme environment of the Qinghai-Tibetan Plateau (QTP). In addition, we identified 291-579 RNA editing sites in the chloroplast genomes of six Isoetes based on transcriptome data, well above the average of angiosperms. RNA editing in protein-coding transcripts results from amino acid changes to increase their hydrophobicity and conservation in Isoetes, which may help proteins form functional three-dimensional structure. Overall, the results of this study provide comprehensive transcriptome and chloroplast genome resources and contribute to a better understanding of adaptive evolutionary and molecular biology in Isoetes.

Phylotranscriptomics Illuminates the Placement of Whole Genome Duplications and Gene Retention in Ferns

Ferns are the second largest clade of vascular plants with over 10,000 species, yet the generation of genomic resources for the group has lagged behind other major clades of plants. Transcriptomic data have proven to be a powerful tool to assess phylogenetic relationships, using thousands of markers that are largely conserved across the genome, and without the need to sequence entire genomes. We assembled the largest nuclear phylogenetic dataset for ferns to date, including 2884 single-copy nuclear loci from 247 transcriptomes (242 ferns, five outgroups), and investigated phylogenetic relationships across the fern tree, the placement of whole genome duplications (WGDs), and gene retention patterns following WGDs. We generated a well-supported phylogeny of ferns and identified several regions of the fern phylogeny that demonstrate high levels of gene tree-species tree conflict, which largely correspond to areas of the phylogeny that have been difficult to resolve. Using a combination of approaches, we identified 27 WGDs across the phylogeny, including 18 large-scale events (involving more than one sampled taxon) and nine small-scale events (involving only one sampled taxon). Most inferred WGDs occur within single lineages (e.g., orders, families) rather than on the backbone of the phylogeny, although two inferred events are shared by leptosporangiate ferns (excluding Osmundales) and Polypodiales (excluding Lindsaeineae and Saccolomatineae), clades which correspond to the majority of fern diversity. We further examined how retained duplicates following WGDs compared across independent events and found that functions of retained genes were largely convergent, with processes involved in binding, responses to stimuli, and certain organelles over-represented in paralogs while processes involved in transport, organelles derived from endosymbiotic events, and signaling were under-represented. To date, our study is the most comprehensive investigation of the nuclear fern phylogeny, though several avenues for future research remain unexplored.

Club-mosses (Diphasiastrum, Lycopodiaceae) from the Far East - Introgression and possible cryptic speciation

Hybridization occurs often in the genus Diphasiastrum (Lycopodiaceae), which corroborates reports for the two other recognized lycophyte families, Isoëtaceae and Selaginellaceae. Here we investigate the case of D. alpinum and D. sitchense from the Russian Far East (Kamchatka). Their hybrid, D. × takedae, was morphologically recognizable in 16 out of 22 accessions showing molecular signatures of hybridization; the remaining accessions displayed the morphology of either D. alpinum (3) or D. sitchense (3). We sequenced markers for chloroplast microsatellites (cp, 175 accessions from Kamchatka) and for the two nuclear markers RPB and LFY (175 and 152 accessions). A selection of 42 accessions, including all hybrid accessions, was analysed via genotyping by sequencing (GBS). We found multiple, but apparently uniparental hybridization, clearly characterized by a deviating group of haplotypes for D. sitchense and all hybrids. All accessions showing molecular signatures of hybridization in nuclear markers revealed the parental haplotype of D. sitchense, however only the LFY marker differentiated between the parent species. GBS, including 69,819 quality-filtered single nucleotid polymorphisms, unambiguously identified the hybrids and revealed introgression to occur. Most of the hybrids were F₁, but three turned out to be backcrosses with D. alpinum (one) and with D. sitchense (two). These observations are in contrast to prior findings on three European species and their intermediates where all three hybrids turned out to be independent F₁ crosses without evidence of recent backcrossing. In this study, backcrossing was detected, which indicates a limited fertility of the hybrid taxon D. × takedae. A comparison of accessions of Kamchatkian D. alpinum with plants from Europe indicated possible cryptic speciation. Accessions from the Far East had (i) a lower DNA content (7.0 vs. 7.5 pg/2C), (ii) different prevailing cp haplotypes, and (iii) RPB genotypes, and (iv) a clearly different SNP pattern in GBS. Diphasiastrum sitchense and the similar D. nikoënse, for the latter additional accessions from Japan were investigated, appeared as forms of one diverse species, sharing genotypes in both nuclear markers, although chloroplast haplotypes and DNA content show slight variations.

Phylogenomic Analysis Reconstructed the Order Matoniales from Paleopolyploidy Veil

Phylogenetic conflicts limit our understanding of the evolution of terrestrial life under multiple whole genome duplication events, and the phylogeny of early terrestrial plants remains full of controversy. Although much incongruence has been solved with so-called robust topology based on single or lower copy genes, the evolutionary mechanisms behind phylogenetic conflicts such as polyploidization remain poorly understood. Here, through decreasing the effects of polyploidization and increasing the samples of species, which represent all four orders and eight families that comprise early leptosporangiate ferns, we have reconstructed a robust phylogenetic tree and network with 1125 1-to-1 orthologs based on both coalescent and concatenation methods. Our data consistently suggest that Matoniales, as a monophyletic lineage including Matoniaceae and Dipteridaceae, should be redefined as an ordinal rank. Furthermore, we have identified and located at least 11 whole-genome duplication events within the evolutionary history of four leptosporangiates lineages, and associated polyploidization with higher speciation rates and mass extinction events. We hypothesize that paleopolyploidization may have enabled leptosporangiate ferns to survive during mass extinction events at the end Permian period and then flourish throughout the Mesozoic era, which is supported by extensive fossil records. Our results highlight how ancient polyploidy can result in rapid species radiation, thus causing phylogenetic conflicts yet allowing plants to survive and thrive during mass extinction events.

Plastid Genome of Equisetum xylochaetum from the Atacama Desert, Chile and the Relationships of Equisetum Based on Frequently Used Plastid Genes and Network Analysis

The modern pteridophyte genus Equisetum is the only survivor of Sphenopsida, an ancient clade known from the Devonian. This genus, of nearly worldwide distribution, comprises approximately 15 extant species. However, genomic information is limited. In this study, we assembled the complete chloroplast genome of the giant species Equisetum xylochaetum from a metagenomic sequence and compared the plastid genome structure and protein-coding regions with information available for two other Equisetum species using network analysis. Equisetum chloroplast genomes showed conserved traits of quadripartite structure, gene content, and gene order. Phylogenetic analysis based on plastome protein-coding regions corroborated previous reports that Equisetum is monophyletic, and that E. xylochaetum is more closely related to E. hyemale than to E. arvense. Single-gene phylogenetic estimation and haplotype analysis showed that E. xylochaetum belonged to the subgenus Hippochaete. Single-gene haplotype analysis revealed that E. arvense, E. hyemale, E. myriochaetum, and E. variegatum resolved more than one haplotype per species, suggesting the presence of a high diversity or a high mutation rate of the corresponding nucleotide sequence. Sequences from E. bogotense appeared as a distinct group of haplotypes representing the subgenus Paramochaete that diverged from Hippochaete and Equisetum. In addition, the taxa that were frequently located at the joint region of the map were E. scirpoides and E. pratense, suggesting the presence of some plastome characters among the Equiseum subgenera.

Phylogeny of Merlin's grass (Isoetaceae): revealing an "Amborella syndrome" and the importance of geographic distribution for understanding current and historical diversity

Background

Merlin’s grass (Isoetes, Isoetaceae, Lycopsida), is the extant remnant of the isoetalean wood-producing lycopsids that originated during the Paleozoic, possibly in aquatic or boggy habitats. Modern day species are aquatic, semi-aquatic or terrestrial and occur almost worldwide. They display little morphological variation; the lobed corm has helically arranged leaves with internal air channels and basal sporangia. Genetic variation has also proven limited, which has hampered phylogenetic inference. We investigate evolutionary relationships in Isoetes, using molecular data and an extended sample of species compared to previous work, adding species that have never before been included in a phylogenetic study.

Results

Our results reveal an unexpected discovery of an “Amborella syndrome” in Isoetaceae: a single poorly known species is sister to the remaining family. The species, Isoetes wormaldii, is a rare endemic to the Eastern Cape of South Africa. Its leaves are flattened with a rounded point, which sharply contrasts with the awl-shaped leaves of most other species of Isoetes. The remaining species of Isoetes are resolved in five major clades, also indicated in previous work. While the phylogeny shows geographic structure, the patterns are complex. For example, tropical-southern African species occur in at least five clades, and Indian, Australian and Mediterranean species in at least three clades each.

Conclusion

The evolutionary and biogeographical history of Isoetes is not easily explained, and may conceivably include ample extinction and a mixture of ancient and more recent processes. Previously shown difficulties with node age estimation increase the problem. The here demonstrated sister-relationship between the phylogenetically, morphologically and genetically distinct Isoetes wormaldii and the remaining family appears to bridge the morphological gap between Isoetes and its extinct relatives, although further studies are needed. Moreover, it shortens the branch length to its living sister genus Selaginella, and may enhance node age estimation in future studies. Isoetes wormaldii is critically endangered, known only from one (to a few) minor populations. Immediate actions need to be taken if we want to prevent this unique species from going extinct.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01988-w.

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

State of the art in cytogenetics, insights into chromosome number evolution, and new C-value reports for the fern family Gleicheniaceae

Studies concerning the cytogenetics of Gleicheniaceae have been scarce, especially those employing evolutionary approaches. Two chromosome number evolutionary models have been hypothesized for Gleicheniaceae. One proposes that ancestral haploid numbers were small and that the chromosome numbers of extant species evolved through polyploidy. The other model proposes that, at the genus level, fern chromosome evolution occurred from ancestors with essentially the same high chromosome numbers seen in living lineages. Neither of those hypotheses has been tested based on phylogenetic frameworks. We sought to (i) present the state of the art of Gleicheniaceae chromosome numbers; (ii) test the two evolutionary models of chromosome numbers within a phylogenetic framework; (iii) test correlations between DNA contents and chromosome numbers in the family. We report here DNA C-values for five species, which increases the number of investigated taxa nearly twofold and report two new genera records. Ancestral state chromosome reconstruction corroborates the hypothesis that ancestral chromosome numbers in Gleicheniaceae were as high as those of extant lineages. Our results demonstrate the possible role of dysploidy in the evolutionary chromosome history of Gleicheniaceae at the genus level and suggest that the relationship between chromosome number and DNA content does not appear to be linear.

Transposable element landscapes illuminate past evolutionary events in the endangered fern Vandenboschia speciosa

Vandenboschia speciosa is an endangered tetraploid fern species with a large genome (10.5 Gb). Its geographical distribution is characterized by disjoined tertiary flora refuges, with relict populations that survived past climate crises. Here, we analyzed the transposable elements (TEs) and found that they comprise approximately 76% of the V. speciosa genome, thus being the most abundant type of DNA sequence in this gigantic genome. The V. speciosa genome is composed of 51% and 5.6% of Class I and Class II elements, respectively. LTR retrotransposons were the most abundant TEs in this species (at least 42% of the genome), followed by non-LTR retrotransposons, which constituted at least 8.7% of the genome of this species. We introduce an additional analysis to identify the nature of non-annotated elements (19% of the genome). A BLAST search of the non-annotated contigs against the V. speciosa TE database allowed for the identification of almost half of them, which were most likely diverged sequence variants of the annotated TEs. In general, the TE composition in V. speciosa resembles the TE composition in seed plants. In addition, repeat landscapes revealed three episodes of amplification for all TEs, most likely due to demographic changes associated with past climate crises.

Transcriptome-wide SNPs for Botrychium lunaria ferns enable fine-grained analysis of ploidy and population structure

Ferns are the second most diverse group of land plants after angiosperms. Extant species occupy a wide range of habitats and contribute significantly to ecosystem functioning. Despite the importance of ferns, most taxa are poorly covered by genomic resources and within‐species studies based on high‐resolution markers are entirely lacking. The genus Botrychium belongs to the family Ophioglossaceae, which includes species with very large genomes and chromosome numbers (e.g., Ophioglossum reticulatum 2n = 1520). The genus has a cosmopolitan distribution with 35 species, half of which are polyploids. Here, we establish a transcriptome for Botrychium lunaria (L.) Sw., a diploid species with an extremely large genome of about ~19.0–23.7 Gb. We assembled 25,677 high‐quality transcripts with an average length of 1,333 bp based on deep RNA‐sequencing of a single individual. We sequenced 11 additional transcriptomes of individuals from two populations in Switzerland, including the population of the reference individual. Based on read mapping to reference transcript sequences, we identified 374,463 single nucleotide polymorphisms (SNPs) segregating among individuals for an average density of 14 SNPs per kilobase. We found that all 12 transcriptomes were most likely from diploid individuals. The transcriptome‐wide markers provided unprecedented resolution of the population genetic structure, revealing substantial variation in heterozygosity among individuals. We also constructed a phylogenomic tree of 92 taxa representing all fern orders to ascertain the placement of the genus Botrychium. High‐quality transcriptomic resources and SNP sets constitute powerful population genomic resources to investigate the ecology, and evolution of fern populations.

Phylogenomic conflict coincides with rapid morphological innovation

Evolutionary biologists have long been fascinated with the episodes of rapid phenotypic innovation that underlie the emergence of major lineages. Although our understanding of the environmental and ecological contexts of such episodes has steadily increased, it has remained unclear how population processes contribute to emergent macroevolutionary patterns. One insight gleaned from phylogenomics is that gene-tree conflict, frequently caused by population-level processes, is often rampant during the origin of major lineages. With the understanding that phylogenomic conflict is often driven by complex population processes, we hypothesized that there may be a direct correspondence between instances of high conflict and elevated rates of phenotypic innovation if both patterns result from the same processes. We evaluated this hypothesis in six clades spanning vertebrates and plants. We found that the most conflict-rich regions of these six clades also tended to experience the highest rates of phenotypic innovation, suggesting that population processes shaping both phenotypic and genomic evolution may leave signatures at deep timescales. Closer examination of the biological significance of phylogenomic conflict may yield improved connections between micro- and macroevolution and increase our understanding of the processes that shape the origin of major lineages across the Tree of Life.

Biogeography and genome size evolution of the oldest extant vascular plant genus, Equisetum (Equisetaceae)

BACKGROUND AND AIMS

Extant plant groups with a long fossil history are key elements in understanding vascular plant evolution. Horsetails (Equisetum, Equisetaceae) have a nearly continuous fossil record dating back to the Carboniferous, but their phylogenetic and biogeographic patterns are still poorly understood. We use here the most extensive phylogenetic analysis to date as a framework to evaluate their age, biogeography and genome size evolution.

METHODS

DNA sequences of four plastid loci were used to estimate divergence times and investigate the biogeographic history of all extant species of Equisetum. Flow cytometry was used to study genome size evolution against the framework of phylogenetic relationships in Equisetum.

KEY RESULTS

On a well-supported phylogenetic tree including all extant Equisetum species, a molecular clock calibrated with multiple fossils places the node at which the outgroup and Equisetum diverged at 343 Mya (Early Carboniferous), with the first major split among extant species occurring 170 Mya (Middle Jurassic). These dates are older than those reported in some other recent molecular clock studies but are largely in agreement with a timeline established by fossil appearance in the geological record. Representatives of evergreen subgenus Hippochaete have much larger genome sizes than those of deciduous subgenus Equisetum, despite their shared conserved chromosome number. Subgenus Paramochaete has an intermediate genome size and maintains the same number of chromosomes.

CONCLUSIONS

The first divergences among extant members of the genus coincided with the break-up of Pangaea and the resulting more humid, warmer climate. Subsequent tectonic activity most likely involved vicariance events that led to species divergences combined with some more recent, long-distance dispersal events. We hypothesize that differences in genome size between subgenera may be related to the number of sperm flagellae.

Reconstructing development of the earliest seed integuments raises a new hypothesis for the evolution of ancestral seed-bearing structures

How plant seeds originated remains unresolved, in part due to disconnects between fossil intermediates and developmental genetics in extant species. The Carboniferous fossil Genomosperma is considered among the most primitive known seeds, with highly lobed integument and exposed nucellus. We have used this key fossil taxon to investigate the evolutionary origins of seed development. We examined sectioned Genomosperma specimens using modern digital 3D reconstruction techniques and established population-level measurements of Genomosperma ovules for quantitative analysis. Genomosperma ovules show significant variation in integumentary lobe fusion and curvature. Our analysis suggests that this variation represents a single species with significant variations in lobe number and fusion, reminiscent of floral development in extant species. We conclude that changes in lobe flexure occurred late in development, consistent with a previously hypothesized function in pollen guidance/retention. We also identify seeds of Genomosperma within cupules for the first time. The presence of a cupule adds evidence towards the plesiomorphy of cupules within seed plants. Together with the similarities identified between the Genomosperma lobed integument and floral organs, we propose that the cupule, integument and nucellus together developed in a shoot-like fashion, potentially ancestral to extant seed plant reproductive shoots.

Plant Taxonomy: A Historical Perspective, Current Challenges, and Perspectives

Taxonomy is the science that explores, describes, names, and classifies all organisms. In this introductory chapter, we highlight the major historical steps in the elaboration of this science, which provides baseline data for all fields of biology and plays a vital role for society but is also an independent, complex, and sound hypothesis-driven scientific discipline.In a first part, we underline that plant taxonomy is one of the earliest scientific disciplines that emerged thousands of years ago, even before the important contributions of the Greeks and Romans (e.g., Theophrastus, Pliny the Elder, and Dioscorides). In the fifteenth-sixteenth centuries, plant taxonomy benefited from the Great Navigations, the invention of the printing press, the creation of botanic gardens, and the use of the drying technique to preserve plant specimens. In parallel with the growing body of morpho-anatomical data, subsequent major steps in the history of plant taxonomy include the emergence of the concept of natural classification , the adoption of the binomial naming system (with the major role of Linnaeus) and other universal rules for the naming of plants, the formulation of the principle of subordination of characters, and the advent of the evolutionary thought. More recently, the cladistic theory (initiated by Hennig) and the rapid advances in DNA technologies allowed to infer phylogenies and to propose true natural, genealogy-based classifications.In a second part, we put the emphasis on the challenges that plant taxonomy faces nowadays. The still very incomplete taxonomic knowledge of the worldwide flora (the so-called taxonomic impediment) is seriously hampering conservation efforts that are especially crucial as biodiversity has entered its sixth extinction crisis. It appears mainly due to insufficient funding, lack of taxonomic expertise, and lack of communication and coordination. We then review recent initiatives to overcome these limitations and to anticipate how taxonomy should and could evolve. In particular, the use of molecular data has been era-splitting for taxonomy and may allow an accelerated pace of species discovery. We examine both strengths and limitations of such techniques in comparison to morphology-based investigations, we give broad recommendations on the use of molecular tools for plant taxonomy, and we highlight the need for an integrative taxonomy based on evidence from multiple sources.

Isolation and functional characterization of four microbial type terpene synthases from ferns

Typical plant terpene synthases (TPSs) are responsible for the production of terpenes, a major class of plant secondary metabolites. However, various nonseed plants also harbor genes encoding microbial terpene synthase-like (MTPSL) enzymes. Here, a scan of 31 ferns transcriptomes revealed 40 sequences putatively encoding MTPSLs. Two groups of sequences were recognized based on the key conserved motifs. Four representative genes were isolated from each of the four species Adiantum capillus-veneris, Cyclosorus parasiticus, Drynaria bonii and Microlepia platyphylla. Following their heterologous expression in E. coli, the recombinant proteins were tested for monoterpene synthase and sesquiterpene synthase activity. These enzymatic products were typical monoterpenes and sesquiterpenes that have been previous shown to be generated by classical plant TPSs when provided with GPP and FPP as substrates. Subcellular localization experiments in the leaf epidermis of Nicotiana benthamiana and onion (Allium cepa) inner epidermal cells indicated that AcMTPSL1 and DbMTPSL were deposited in both the cytoplasm and nucleus, whereas CpMTPSL1 and MpMTPSL were localized in the cytoplasm, chloroplasts and nucleus. AcMTPSL1 was up-regulated in plants exposed to methyl jasmonate treatment, suggesting a role for this gene in host defense. This study provides more information about the catalytic function of MTPSLs in nonseed plants and for the first time, the subcellular localization of MTPSLs was experimentally characterized.

Three highly acidic Equisetum XTHs differ from hetero-trans-β-glucanase in donor substrate specificity and are predominantly xyloglucan homo-transglucosylases

Transglycanases are enzymes that remodel the primary cell wall in plants, potentially loosening and/or strengthening it. Xyloglucan endotransglucosylase (XET; EC 2.4.1.207), ubiquitous in land plants, is a homo-transglucanase activity (donor, xyloglucan; acceptor, xyloglucan) exhibited by XTH (xyloglucan endotransglucosylase/hydrolase) proteins. By contrast, hetero-trans-β-glucanase (HTG) is the only known enzyme that is preferentially a hetero-transglucanase. Its two main hetero-transglucanase activities are MLG : xyloglucan endotransglucosylase (MXE) and cellulose : xyloglucan endotransglucosylase (CXE). HTG is highly acidic and found only in the evolutionarily isolated genus of fern-allies, Equisetum. We now report genes for three new highly acidic HTG-related XTHs in E. fluviatile (EfXTH-A, EfXTH-H and EfXTH-I). We expressed them heterologously in Pichia and tested the encoded proteins' enzymic activities to determine whether their acidity and/or their Equisetum-specific sequences might confer high hetero-transglucanase activity. Untransformed Pichia was found to secrete MLG-degrading enzyme(s), which had to be removed for reliable MXE assays. All three acidic EfXTHs exhibited very predominantly XET activity, although low but measurable hetero-transglucanase activities (MXE and CXE) were also detected in EfXTH-H and EfXTH-I. We conclude that the extremely high hetero-transglucanase activities of Equisetum HTG are not emulated by similarly acidic Equisetum XTHs that share up to 55.5% sequence identity with HTG.

Evolutionary relationships of the ancient fern lineage the adder's tongues (Ophioglossaceae) with description of Sahashia gen. nov

As an ancient lineage of ferns, Ophioglossaceae are evolutionarily among the most fascinating because they have the highest chromosome count of any known organism as well as the presence of sporophores, subterranean gametophytes, eusporangiate sporangia without annuli, and endophytic fungi. Previous studies have produced conflicting results, identifyingsome lineages with unresolved relationships, and have paid much attention to the subfamily Botrychioideae. But the other species-rich subfamily, Ophioglossoideae, has remained largely understudied and only up to 12 accessions of Ophioglossoideae have been sampled. In this study, DNA sequences of seven plastid markers of 149 accessions (75 in Ophioglossoideae) representing approximately 82 species (approximately 74% of estimated species diversity sensu J. Syst. Evol., 2016, 54, 563) in the family, and two Marattiaceae and two Psilotaceae, are used to infer a phylogeny. Our major results include: (1) Ophioglossaceae are resolved as monophyletic with strong support, and so are all four subfamilies and genera sensu PPG I except Botrypus and Ophioglossum; (2) a new genus Sahashia is segregated from Botrypus so that the monophyly of Botrypus can be retained; (3) the monophyly of Ophioglossum in its current circumscription is uncertain in spite of our large character sampling; (4) there is substantial cryptic speciation in Ophioderma detected by our molecular and morphological study; (5) the recognition of Holubiella is advocated based on its morphology and its sister relationship with Sceptridium; and (6) a novel sister relationship between Botrychium and the JHS clade (Japanobotrychium + (Holubiella + Sceptridium)) is discovered.

Keraphyton gen. nov., a new Late Devonian fern-like plant from Australia

The first plants related to the ferns are represented by several extinct groups that emerged during the Devonian. Among them, the iridopterids are closely allied to the sphenopsids, a group represented today by the genus Equisetum. They have been documented in Middle to early Late Devonian deposits of Laurussia and the Kazakhstan plate. Their Gondwanan record is poor, with occurrences limited to Venezuela and Morocco. Here we describe a new genus from a late Late Devonian locality of New South Wales. It is represented by a single anatomically preserved large stem characterized by a star-shaped vascular system with protoxylem strands located at rib tips, and by a lack of secondary tissues. Within the first fern-like plants, this stem shares the largest number of characters with iridopterid axes but differs by the pattern of its vascular system. Keraphyton mawsoniae gen. et sp. nov. adds a new record of early fern-like plants in eastern Gondwana. It provides new insights into the anatomical diversity within this key group of plants and supports the distinctiveness of the Australian flora in the latest Devonian.

Patterns and Rates of Plastid rps12 Gene Evolution Inferred in a Phylogenetic Context using Plastomic Data of Ferns

The trans-splicing rps12 gene of fern plastomes (plastid genomes) exhibits a unique structure owing to its variations in intragenic exon location and intron content, and thus, it provides an excellent model system for examining the effect of plastid gene structure on rates and patterns of molecular evolution. In this study, 16 complete fern plastome sequences were newly generated via the Illumina HiSeq sequencing platform. We reconstructed the phylogeny of ferns and inferred the patterns and rates of plastid rps12 gene evolution in a phylogenetic context by combining these plastome data with those of previously published fern species. We uncovered the diversity of fern plastome evolution by characterizing the structures of these genomes and obtained a highly supported phylogenetic framework for ferns. Furthermore, our results revealed molecular evolutionary patterns that were completely different from the patterns revealed in previous studies. There were significant differences in the patterns and rates of nucleotide substitutions in both intron-containing and intron-less rps12 alleles. Rate heterogeneity between single-copy (SC) and inverted repeat (IR) exons was evident. Unexpectedly, however, IR exons exhibited significantly higher synonymous substitution rates (dS) than SC exons, a pattern that contrasts the regional effect responsible for decreased rates of nucleotide substitutions in IRs. Our results reveal that structural changes in plastid genes have important effects on evolutionary rates, and we propose possible mechanisms to explain the variations in the nucleotide substitution rates of this unusual gene.

Taxonomic, Phylogenetic, and Functional Diversity of Ferns at Three Differently Disturbed Sites in Longnan County, China

Human disturbances are greatly threatening to the biodiversity of vascular plants. Compared to seed plants, the diversity patterns of ferns have been poorly studied along disturbance gradients, including aspects of their taxonomic, phylogenetic, and functional diversity. Longnan County, a biodiversity hotspot in the subtropical zone in South China, was selected to obtain a more thorough picture of the fern–disturbance relationship, in particular, the taxonomic, phylogenetic, and functional diversity of ferns at different levels of disturbance. In 90 sample plots of 5 × 5 m2 along roadsides at three sites, we recorded a total of 20 families, 50 genera, and 99 species of ferns, as well as 9759 individual ferns. The sample coverage curve indicated that the sampling effort was sufficient for biodiversity analysis. In general, the taxonomic, phylogenetic, and functional diversity measured by Hill numbers of order q = 0–3 indicated that the fern diversity in Longnan County was largely influenced by the level of human disturbance, which supports the ‘increasing disturbance hypothesis’. Many functional traits of ferns at the most disturbed site were adaptive to the disturbance. There were also some indicators of fern species responding to the different disturbance levels. Hence, ferns may be considered as a good indicator group for environmental stress.

Identification and evolutionary analysis of chalcone isomerase-fold proteins in ferns

The distribution of type I and II chalcone isomerases (CHIs) in plants is highly family specific. We have previously reported that ancient land plants, such as the liverworts and Selaginella moellendorffii, harbor type II CHIs. To better understand the function and evolution of CHI-fold proteins, transcriptomic data obtained from 52 pteridophyte species were subjected to sequence alignment and phylogenetic analysis. The residues determining type I/II CHI identity in the pteridophyte CHIs were identical to those of type I CHIs. The enzymatic characterization of a sample of 24 CHIs, representing all the key pteridophyte lineages, demonstrated that 19 of them were type I enzymes and that five exhibited some type II activity due to an amino acid mutation. Two pteridophyte chalcone synthases (CHSs) were also characterized, and a type IV CHI (CHIL) was demonstrated to interact physically with CHSs and CHI, and to increase CHS activity by decreasing derailment products, thus enhancing flavonoid production. These findings suggest that the emergence of type I CHIs may have coincided with the divergence of the pteridophytes. This study deepens our understanding of the molecular mechanism of CHIL as an enhancer in the flavonoid biosynthesis pathway.

Lycophytes of the Chapada das Mesas National Park, Cerrado, Maranhão, Brazil

Abstract: Despite an increase in studies involving lycophytes and ferns in recent years, Maranhão is still one of the northeastern states with a poorly known flora, mainly due to low sampling. For Chapada das Mesas National Park, a protected conservation unit in the Cerrado maranhense with phytophysiognomies that favor floristic diversity, there are no floristic studies about lycophytes. This study aimed to conduct the floristic inventory and taxonomic study of lycophytes in Chapada das Mesas National Park. Three collection expeditions, of three days each, were made between March 2017 and June 2018. Three families of lycophytes represented by three genera and six species were identified in the study area. The most diverse family was Selaginellaceae with four species: Selaginella conduplicata, S. erythropus (Mart.) Spring, S. radiata (Aubl.) Baker and S. simplex Baker. The other two families, Isoëtaceae and Lycopodiaceae, had one species each, Isoëtes sp. and Palhinhaea cernua (L.) Franco & Vasc. This study also reports two new record for the Brazilian Cerrado phytogeographic domain, S. conduplicate and S. radiata, of which the first is a new record for Maranhão State. Most specimens were found in riparian forest on rocky outcrops.

Importance of Sporopollenin Structure and Accessibility in the Sorption of Phenanthrene by Biota Spores and Pollens

Although spores/pollens are so abundant and ubiquitous in the environment, the role of these natural organic matter concerning fate and transport of organic pollutants in the environment is neglected. Lipid-free fractions and sporopollenins were isolated from seven spores/pollens collected from lower and higher biota species and were characterized by elemental analysis, CO₂ adsorption techniques, and advanced solid-state ¹³C nuclear magnetic resonance spectroscopy. Then, the sorption isotherms of phenanthrene (Phen) on all the samples were investigated by a batch technique. The sporopollenins were a highly cross-linked polymer including alkyl carbon, poly(methylene) carbon, and aromatic carbon as well as oxygen functionalities; additionally, their sorption capacities (K_oc) for Phen reached up to 1 170 000 mL/g, suggesting that some of the sporopollenins were good biopolymeric sorbents for the removal of hydrophobic organic contaminants in aquatic media. A highly significant and positive correlation between the sorption capacity of Phen and the aliphaticity of the sporopollenins suggested that their structure was critical to Phen sorption. Meanwhile, the (O + N)/C atomic ratios and polar groups were significantly and negatively correlated with the sorption capacity of Phen, indicating that accessibility also played a significant role in the sorption process. Moreover, variable correlations between the sorption capacities (K_oc) and the micropore volumes of the spore/pollen fractions were observed. This study sheds light on the importance of the polarity, microporosity, and structure of sporopollenins in the sorption process of Phen.

Virtual issue: Ecology and evolution of pteridophytes in the era of molecular genetics

The past quarter-century has witnessed a revolution in our understanding of the phylogenetics, systematics, and ecology of pteridophytes (ferns and lycophytes), particularly due to the rapid accumulation of plastid sequence data and a renewed interest in the ecology of the sexual phase of the life cycle. We here compile 19 papers recently published in the Journal of Plant Research dealing with the biology of pteridophytes, grouped into six categories: (1) breeding systems, (2) species complexes and polyploidization, (3) fossil taxa, (4) gametophyte ecology, (5) systematics, (6) biodiversity. We hope this collection of papers will be of value to researchers interested in this fascinating group of plants.

Plastid genome and composition analysis of two medical ferns: Dryopteris crassirhizoma Nakai and Osmunda japonica Thunb

Background

Dryopteris crassirhizoma Nakai and Osmunda japonica Thunb. are ferns that are popularly used for medicine, as recorded by the Chinese pharmacopoeia, and are distributed in different regions of China. However, O. japonica is not record in the Standards of Chinese Herbal Medicines in Hong Kong. Research on identification methods of D. crassirhizoma and O. japonica is necessary and the phylogenetic position of the two species should be identified. The plastid genome is structurally highly conserved, providing valuable sources of genetic markers for phylogenetic analyses and development of molecule makers for identification.

Methods

The plastid genome DNA was extracted from both fern species and then sequenced on the Illumina Hiseq 4000. Sequences were assembled into contigs by SOAPdenovo2.04, aligned to the reference genome using BLAST, and then manually corrected. Genome annotation was performed by the online DOGMA tool. General characteristics of the plastid genomes of the two species were analyzed and compared with closely related species. Additionally, phylogenetical trees were reconstructed by maximum likelihood methods. The content of dryocrassin of the two species were determined according to the Standards of Chinese Herbal Medicines in Hong Kong.

Results

The genome structures of D. crassirhizoma and O. japonica have different characteristics including the genome size, the size of each area, gene location, and types. Moreover, the (simple sequence repeats) SSRs of the plastid genomes were more similar to other species in the same genera. Compared with D. fragrans, D. crassirhizoma shows an inversion (approximately 1.6 kb), and O. japonica shows two inversions (1.9 kb and 216 bp). The nucleotide diversity (polymorphism information, Pi) analysis showed that the psbK gene and rpl14-rpl16 region have the highest Pi value in Dryopteris, and the ycf2-CDS3 and rpl14-rpl16 regions show the highest Pi vale in O. japonica. Phylogenetic analyses showed that the two species were grouped in two separate clades from each other, with both individually located with other members of their genus. The marker content of dryocrassin is not found in O. japonica.

Conclusions

The study is the first to identify plastid genome features of D. crassirhizoma and O. japonica. The results may provide a theoretical basis for the identification and the application of the two medically important fern species.

Electronic supplementary material

The online version of this article (10.1186/s13020-019-0230-4) contains supplementary material, which is available to authorized users.

Full plastome sequence of the fern Vandenboschia speciosa (Hymenophyllales): structural singularities and evolutionary insights

We provide here the first full chloroplast genome sequence, i.e., the plastome, for a species belonging to the fern order Hymenophyllales. The phylogenetic position of this order within leptosporangiate ferns, together with the general scarcity of information about fern plastomes, places this research as a valuable study on the analysis of the diversity of plastomes throughout fern evolution. Gene content of V. speciosa plastome was similar to that in most ferns, although there were some characteristic gene losses and lineage-specific differences. In addition, an important number of genes required U to C RNA editing for proper protein translation and two genes showed start codons alternative to the canonical AUG (AUA). Concerning gene order, V. speciosa shared the specific 30-kb inversion of euphyllophytes plastomes and the 3.3-kb inversion of fern plastomes, keeping the ancestral gene order shared by eusporangiate and early leptosporangiate ferns. Conversely, V. speciosa has expanded IR regions comprising the rps7, rps12, ndhB and trnL genes in addition to rRNA and other tRNA genes, a condition shared with several eusporangiate ferns, lycophytes and hornworts, as well as most seed plants.

Assessing Chemical Diversity in Psilotum nudum (L.) Beauv., a Pantropical Whisk Fern That Has Lost Many of Its Fern-Like Characters

Members of the Psilotales (whisk ferns) have a unique anatomy, with conducting tissues but lacking true leaves and roots. Based on recent phyogenies, these features appear to represent a reduction from a more typical modern fern plant rather than the persistence of ancestral features. In this study, extracts of several Psilotum organs and tissues were analyzed by Gas Chromatography - Mass Spectrometry (GC-MS) and High Performance Liquid Chromatography - Quadrupole Time of Flight - Mass Spectrometry (HPLC-QTOF-MS). Some arylpyrones and biflavonoids had previously been reported to occur in Psilotum and these metabolite classes were found to be prominent constituents in the present study. Some of these were enriched and further characterized by Nuclear Magnetic Resonance (NMR) spectroscopy. HPLC-QTOF-MS and NMR data were searched against an updated Spektraris database (expanded by incorporating over 300 new arylpyrone and biflavonoid spectral records) to aid significantly with peak annotation. Principal Component Analysis (PCA) with combined GC-MS and HPLC-QTOF-MS data sets obtained with several Psilotum organs and tissues indicated a clear separation of the sample types. The principal component scores for below-ground rhizome samples corresponded to the vectors for carbohydrate monomers and dimers and small organic acids. Above-ground rhizome samples had principal component scores closer to the direction of vectors for arylpyrone glycosides and sucrose (which had high concentrations in above-and below-ground rhizomes). The unique position of brown synangia in a PCA plot correlated with the vector for biflavonoid glycosides. Principal component scores for green and yellow synangia correlated with the direction of vectors for arylpyrone glycosides and biflavonoid aglycones. Localization studies with cross sections of above-ground rhizomes, using Matrix-Assisted Laser Desorption/Ionization - Mass Spectrometry (MALDI-MS), provided evidence for a preferential accumulation of arylpyrone glycosides and biflavonoid aglycones in cells of the chlorenchyma. Our results indicate a differential localization of metabolites with potentially tissue-specific functions in defenses against biotic and abiotic stresses. The data are also a foundation for follow-up work to better understand chemical diversity in the Psilotales and other members of the fern lineage.

Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution

We present reference-quality genome assembly and annotation for the stout camphor tree (Cinnamomum kanehirae (Laurales, Lauraceae)), the first sequenced member of the Magnoliidae comprising four orders (Laurales, Magnoliales, Canellales and Piperales) and over 9,000 species. Phylogenomic analysis of 13 representative seed plant genomes indicates that magnoliid and eudicot lineages share more recent common ancestry than monocots. Two whole-genome duplication events were inferred within the magnoliid lineage: one before divergence of Laurales and Magnoliales and the other within the Lauraceae. Small-scale segmental duplications and tandem duplications also contributed to innovation in the evolutionary history of Cinnamomum. For example, expansion of the terpenoid synthase gene subfamilies within the Laurales spawned the diversity of Cinnamomum monoterpenes and sesquiterpenes.

Large-scale phylogenomic analysis resolves a backbone phylogeny in ferns

Background

Ferns, originated about 360 million years ago, are the sister group of seed plants. Despite the remarkable progress in our understanding of fern phylogeny, with conflicting molecular evidence and different morphological interpretations, relationships among major fern lineages remain controversial.

Results

With the aim to obtain a robust fern phylogeny, we carried out a large-scale phylogenomic analysis using high-quality transcriptome sequencing data, which covered 69 fern species from 38 families and 11 orders. Both coalescent-based and concatenation-based methods were applied to both nucleotide and amino acid sequences in species tree estimation. The resulting topologies are largely congruent with each other, except for the placement of Angiopteris fokiensis, Cheiropleuria bicuspis, Diplaziopsis brunoniana, Matteuccia struthiopteris, Elaphoglossum mcclurei, and Tectaria subpedata.

Conclusions

Our result confirmed that Equisetales is sister to the rest of ferns, and Dennstaedtiaceae is sister to eupolypods. Moreover, our result strongly supported some relationships different from the current view of fern phylogeny, including that Marattiaceae may be sister to the monophyletic clade of Psilotaceae and Ophioglossaceae; that Gleicheniaceae and Hymenophyllaceae form a monophyletic clade sister to Dipteridaceae; and that Aspleniaceae is sister to the rest of the groups in eupolypods II. These results were interpreted with morphological traits, especially sporangia characters, and a new evolutionary route of sporangial annulus in ferns was suggested. This backbone phylogeny in ferns sets a foundation for further studies in biology and evolution in ferns, and therefore in plants.

Genetic diversity and hybrid formation in Central European club-mosses (Diphasiastrum, Lycopodiaceae) - New insights from cp microsatellites, two nuclear markers and AFLP

In Europe, the genus Diphasiastrum (Lycopodiophyta) forms a reticulate network of six diploid taxa, including three parent species (D. alpinum, D. complanatum and D. tristachyum) and three hybrids (D. × issleri, D. × oellgaardii and D. × zeilleri). It was not clear if the hybrids arose once or repeatedly, if they have reproductive competence and if backcrossing occurs. We addressed these questions by analysing 209 accessions for chloroplast microsatellites (cp), two nuclear markers (introns of the RPB and LFY genes) and AFLP. For D. complanatum we show a sexual life cycle with alternation of generations: the gametophytic DNA amount is half of that of the sporophyte. With the exception of a single accession all hybrids display one of the two parental cp haplotypes; their frequencies do not differ significantly from a 1:1 ratio. Genotypes of nuclear markers are species-specific, displaying 2/4/1 (RPB) and 1/8/1 alleles (LFY) for the three parents mentioned above; all hybrids have one allele from each parent. All three hybrid taxa apparently represent independent F1 crosses. Hybridisation occurs bidirectional; no evidence for recent backcrossing was found. Asexual reproduction via agamospory is at least rare, since AFLP showed all hybrid plants to be different.

Order-level fern plastome phylogenomics: new insights from Hymenophyllales

PREMISE OF THE STUDY

Filmy ferns (Hymenophyllales) are a highly specialized lineage, having mesophyll one-cell layer thick and inhabiting particularly shaded and humid environments. The phylogenetic placement of Hymenophyllales has been inconclusive, and while over 87 whole fern plastomes have been published, none was from Hymenophyllales. To better understand the evolutionary history of filmy ferns, we sequenced the first complete plastome for this order.

METHODS

We compiled a phylogenomic plastome data set encompassing all 11 fern orders, and reconstructed phylogenies using different data types (nucleotides, codons, and amino acids) and partition schemes (codon positions and loci). To infer the evolution of fern plastome organization, we coded plastome features, including inversions, inverted repeat boundary shifts, gene losses, and tRNA anticodon sequences as characters, and reconstructed the ancestral states for these characters.

KEY RESULTS

We discovered a suite of novel, Hymenophyllales-specific plastome structures that likely resulted from repeated expansions and contractions of the inverted repeat regions. Our phylogenetic analyses reveal that Hymenophyllales is highly supported as either sister to Gleicheniales or to Gleicheniales + the remaining non-Osmundales leptosporangiates, depending on the data type and partition scheme.

CONCLUSIONS

Although our analyses could not confidently resolve the phylogenetic position of Hymenophyalles, the results here highlight the danger of drawing conclusions from "all-in" phylogenomic data set without exploring potential inconsistencies in the data. Finally, our first order-level reconstruction of fern plastome structural evolution provides a useful framework for future plastome research.

Origin of Equisetum: Evolution of horsetails (Equisetales) within the major euphyllophyte clade Sphenopsida

PREMISE OF THE STUDY

Equisetum is the sole living representative of Sphenopsida, a clade with impressive species richness, a long fossil history dating back to the Devonian, and obscure relationships with other living pteridophytes. Based on molecular data, the crown group age of Equisetum is mid-Paleogene, although fossils with possible crown synapomorphies appear in the Triassic. The most widely circulated hypothesis states that the lineage of Equisetum derives from calamitaceans, but no comprehensive phylogenetic studies support the claim. Using a combined approach, we provide a comprehensive phylogenetic analysis of Equisetales, with special emphasis on the origin of genus Equisetum.

METHODS

We performed parsimony phylogenetic analyses to address relationships of 43 equisetalean species (15 extant, 28 extinct) using a combination of morphological and molecular characters.

KEY RESULTS

We recovered Equisetaceae + Neocalamites as sister to Calamitaceae + a clade of Angaran and Gondwanan horsetails, with the four groups forming a clade that is sister to Archaeocalamitaceae. The estimated age for the Equisetum crown group is mid-Mesozoic.

CONCLUSIONS

Modern horsetails are not nested within calamitaceans; instead, both groups have explored independent evolutionary trajectories since the Carboniferous. Diverse fossil taxon sampling helps to shed light on the position and relationships of equisetalean lineages, of which only a tiny remnant is present within the extant flora. Understanding these relationships and early character configurations of ancient plant clades as Equisetales provide useful tests of hypotheses about overall phylogenetic relationships of euphyllophytes and foundations for future tests of molecular dates with paleontological data.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.

Recruitment of IC-WOX Genes in Root Evolution

Root evolution has resulted in the extant bifurcating roots in lycophytes, adventitious/lateral roots in euphyllophytes (ferns and seed plants), and primary roots in seed plants. Here, we hypothesize a role for intermediate-clade-WUSCHEL-RELATED HOMEOBOX (IC-WOX) genes in root evolution. IC-WOX might not be specifically involved in lycophyte bifurcation rooting. In the fern Ceratopteris richardii, IC-WOX is expressed in adventitious/lateral root founder cells. In the seed plant Arabidopsis thaliana, there are two IC-WOX subclades, AtWOX11/12 and AtWOX8/9, in adventitious and primary root founder cells, respectively. Thus, IC-WOX was recruited in the common ancestor of ferns and seed plants for adventitious/lateral root organogenesis and evolved into two subclades in seed plants: one was retained in adventitious root organogenesis, while the other was recruited for primary root organogenesis.

The Huperzia selago Shoot Tip Transcriptome Sheds New Light on the Evolution of Leaves

Lycopodiophyta—consisting of three orders, Lycopodiales, Isoetales and Selaginellales, with different types of shoot apical meristems (SAMs)—form the earliest branch among the extant vascular plants. They represent a sister group to all other vascular plants, from which they differ in that their leaves are microphylls—that is, leaves with a single, unbranched vein, emerging from the protostele without a leaf gap—not megaphylls. All leaves represent determinate organs originating on the flanks of indeterminate SAMs. Thus, leaf formation requires the suppression of indeterminacy, that is, of KNOX transcription factors. In seed plants, this is mediated by different groups of transcription factors including ARP and YABBY.

We generated a shoot tip transcriptome of Huperzia selago (Lycopodiales) to examine the genes involved in leaf formation. Our H. selago transcriptome does not contain any ARP homolog, although transcriptomes of Selaginella spp. do. Surprisingly, we discovered a YABBY homolog, although these transcription factors were assumed to have evolved only in seed plants.

The existence of a YABBY homolog in H. selago suggests that YABBY evolved already in the common ancestor of the vascular plants, and subsequently was lost in some lineages like Selaginellales, whereas ARP may have been lost in Lycopodiales. The presence of YABBY in the common ancestor of vascular plants would also support the hypothesis that this common ancestor had a simplex SAM. Furthermore, a comparison of the expression patterns of ARP in shoot tips of Selaginella kraussiana (Harrison CJ, etal. 2005. Independent recruitment of a conserved developmental mechanism during leaf evolution. Nature 434(7032):509–514.) and YABBY in shoot tips of H. selago implies that the development of microphylls, unlike megaphylls, does not seem to depend on the combined activities of ARP and YABBY. Altogether, our data show that Lycopodiophyta are a diverse group; so, in order to understand the role of Lycopodiophyta in evolution, representatives of Lycopodiales, Selaginellales, as well as of Isoetales, have to be examined.