The C-Fern (Ceratopteris richardii) genome: insights into plant genome evolution with the first partial homosporous fern genome assembly

Crosslinking of base modified RNAs by synthetic DYW-KP base editors implicates an enzymatic lysine as the nitrogen donor for U-to-C RNA editing

Base editing mechanisms are being investigated as potential therapeutic tools to alleviate genetic diseases. Sequence specific C-to-U and A-to-I base editing tools are capable of altering RNA and DNA sequences and utilize a hydrolytic deamination mechanism requiring an active site zinc ion and a glutamate residue. In plant organelles, DYW-PG domain containing enzymes catalyze C-to-U edits and likely use the canonical deamination mechanism. Proteins developed from consensus sequences from the DYW-KP domain family catalyze what initially appeared to be U-to-C edits leading to this investigation into mechanistic insights into U-to-C editing. The synthetic DYW-KP enzyme KP6 was found sufficient for C-to-U editing activity stimulated by the addition of carboxylic acids in vitro. Despite addition of 14 different putative amine/amide donors in vitro U-to-C editing could not be observed. C-to-U editing was found not to be concomitant with U-to-C editing thus discounting a pyrimidine transaminase mechanism. RNAs containing base modifications were covalently crosslinked to KP6, KP2, and KP3 proteins. Mass spectrometry of purified KP2 and KP6 proteins revealed an additional mass of 319 Da. A U-to-C crosslinking mechanism was projected to explain the link between crosslinking, RNA base changes, and the ∼319 Da mass. In this model an enzymatic lysine attacks C4 of uridine to form a Schiff base RNA-protein conjugate. Sequenced RT-PCR products from the fern Ceratopteris richardii indicate U-to-C base edits do not preserve proteinaceous crosslinks in planta. Hydrolysis of a protonated Schiff base conjugate releasing cytidine is hypothesized to explain the completed pathway in plants.

Conserved autophagy and diverse cell wall composition: unifying features of vascular tissues in evolutionarily distinct plants

BACKGROUND AND AIMS

The formation of multifunctional vascular tissues represents a significant advancement in plant evolution. Differentiation of conductive cells is specific, involving two main pathways, namely protoplast clearance and cell wall modification. In xylogenesis, autophagy is a crucial process for complete protoplast elimination in tracheary elements, whose cell wall also undergoes strong changes. Knowledge pertaining to living sieve elements, which lose most of their protoplast during phloemogenesis, remains limited. We hypothesized that autophagy plays a crucial role, not only in complete cytoplasmic clearance in xylem but also in partial degradation in phloem. Cell wall elaborations of mature sieve elements are not so extensive. These analyses performed on evolutionarily diverse model species potentially make it possible to understand phloemogenesis to an equal extent to xylogenesis.

METHODS

We investigated the distribution of ATG8 protein, which is an autophagy marker, and cell wall components in the roots of ferns, gymnosperms and angiosperms (monocots, dicot herbaceous plants and trees). Furthermore, we conducted a bioinformatic analysis of complete data on ATG8 isoforms for Ceratopteris richardii.

KEY RESULTS

The presence of ATG8 protein was confirmed in both tracheary elements and sieve elements; however, the composition of cell wall components varied considerably among vascular tissues in the selected plants. Arabinogalactan proteins and β-1,4-galactan were detected in the roots of all studied species, suggesting their potential importance in phloem formation or function. In contrast, no evolutionary pattern was observed for xyloglucan, arabinan or homogalacturonan.

CONCLUSIONS

Our findings indicate that the involvement of autophagy in plants is universal during the development of tracheary elements that are dead at maturity and sieve elements that remain alive. Given the conserved nature of autophagy and its function in protoplast degradation for uninterrupted flow, autophagy might have played a vital role in the development of increasingly complex biological organizations, including the formation of vascular tissues. However, different cell wall compositions of xylem and phloem in different species might indicate diverse functionality and potential for substance transport, which is crucial in plant evolution.

Differential expression and co-localization of transcriptional factors during callus transition to differentiation for shoot organogenesis in the water fern Ceratopteris richardii

BACKGROUND AND AIMS

In flowering plants, regeneration can be achieved by a variety of approaches, and different sets of transcriptional factors are involved in these processes. However, regeneration in taxa other than flowering plants remains a mystery. Ceratopteris richardii is a representative fern capable of both direct and indirect organogenesis, and we aimed to investigate the genetics underlying the transition from callus proliferation to differentiation.

METHODS

Morphological and histological analyses were used to determine the type of regeneration involved. RNA sequencing and differential gene expression were used to investigate how the callus switches from proliferation to differentiation. Phylogenetic analysis and RNA in situ hybridization were used to understand whether transcriptional factors are involved in this transition.

KEY RESULTS

The callus formed on nascent leaves and subsequently developed the shoot pro-meristem and shoot meristem, thus completing indirect de novo shoot organogenesis in C. richardii. Genes were differentially expressed during the callus transition from proliferation to differentiation, indicating a role for photosynthesis, stimulus response and transmembrane signalling in this transition and the involvement of almost all cell layers that make up the callus. Transcriptional factors were either downregulated or upregulated, which were generally in many-to-many orthology with genes known to be involved in callus development in flowering plants, suggesting that the genetics of fern callus development are both conserved and divergent. Among them, an STM-like, a PLT-like and an ethylene- and salt-inducible ERF gene3-like gene were expressed simultaneously in the vasculature but not in the other parts of the callus, indicating that the vasculature played a role in the callus transition from proliferation to differentiation.

CONCLUSIONS

Indirect de novo shoot organogenesis could occur in ferns, and the callus transition from proliferation to differentiation required physiological changes, differential expression of transcriptional factors and involvement of the vasculature.

Re-examining meristems through the lens of evo-devo

The concept of the meristem was introduced in 1858 to characterize multicellular, formative, and proliferative tissues that give rise to the entire plant body, based on observations of vascular plants. Although its original definition did not encompass bryophytes, this concept has been used and continuously refined over the past 165 years to describe the diverse apices of all land plants. Here, we re-examine this matter in light of recent evo-devo research and show that, despite displaying high anatomical diversity, land plant meristems are unified by shared genetic control. We also propose a modular view of meristem function and highlight multiple evolutionary mechanisms that are likely to have contributed to the assembly and diversification of the varied meristems during the course of plant evolution.

Recurrent evolutionary switches of mitochondrial cytochrome c maturation systems in Archaeplastida

Mitochondrial cytochrome c maturation (CCM) requires heme attachment via distinct pathways termed systems I and III. The mosaic distribution of these systems in Archaeplastida raises questions about the genetic mechanisms and evolutionary forces promoting repeated evolution. Here, we show a recurrent shift from ancestral system I to the eukaryotic-specific holocytochrome c synthase (HCCS) of system III in 11 archaeplastid lineages. Archaeplastid HCCS is sufficient to rescue mutants of yeast system III and Arabidopsis system I. Algal HCCS mutants exhibit impaired growth and respiration, and altered biochemical and metabolic profiles, likely resulting from deficient CCM and reduced cytochrome c-dependent respiratory activity. Our findings demonstrate that archaeplastid HCCS homologs function as system III components in the absence of system I. These results elucidate the evolutionary trajectory and functional divergence of CCM pathways in Archaeplastida, providing insight into the causes, mechanisms, and consequences of repeated cooption of an entire biological pathway.

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.

A positive feedback regulation of SnRK1 signaling by autophagy in plants

SnRK1, an evolutionarily conserved heterotrimeric kinase complex that acts as a key metabolic sensor in maintaining energy homeostasis in plants, is an important upstream activator of autophagy that serves as a cellular degradation mechanism for the healthy growth of plants. However, whether and how the autophagy pathway is involved in regulating SnRK1 activity remains unknown. In this study, we identified a clade of plant-specific and mitochondria-localized FCS-like zinc finger (FLZ) proteins as currently unknown ATG8-interacting partners that actively inhibit SnRK1 signaling by repressing the T-loop phosphorylation of the catalytic α subunits of SnRK1, thereby negatively modulating autophagy and plant tolerance to energy deprivation caused by long-term carbon starvation. Interestingly, these AtFLZs are transcriptionally repressed by low-energy stress, and AtFLZ proteins undergo a selective autophagy-dependent pathway to be delivered to the vacuole for degradation, thereby constituting a positive feedback regulation to relieve their repression of SnRK1 signaling. Bioinformatic analyses show that the ATG8-FLZ-SnRK1 regulatory axis first appears in gymnosperms and seems to be highly conserved during the evolution of seed plants. Consistent with this, depletion of ATG8-interacting ZmFLZ14 confers enhanced tolerance, whereas overexpression of ZmFLZ14 leads to reduced tolerance to energy deprivation in maize. Collectively, our study reveals a previously unknown mechanism by which autophagy contributes to the positive feedback regulation of SnRK1 signaling, thereby enabling plants to better adapt to stressful environments.

The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

DELLA proteins regulate spore germination and reproductive development in Physcomitrium patens

Proteins of the DELLA family integrate environmental signals to regulate growth and development throughout the plant kingdom. Plants expressing non-degradable DELLA proteins underpinned the development of high-yielding 'Green Revolution' dwarf crop varieties in the 1960s. In vascular plants, DELLAs are regulated by gibberellins, diterpenoid plant hormones. How DELLA protein function has changed during land plant evolution is not fully understood. We have examined the function and interactions of DELLA proteins in the moss Physcomitrium (Physcomitrella) patens, in the sister group of vascular plants (Bryophytes). PpDELLAs do not undergo the same regulation as flowering plant DELLAs. PpDELLAs are not degraded by diterpenes, do not interact with GID1 gibberellin receptor proteins and do not participate in responses to abiotic stress. PpDELLAs do share a function with vascular plant DELLAs during reproductive development. PpDELLAs also regulate spore germination. PpDELLAs interact with moss-specific photoreceptors although a function for PpDELLAs in light responses was not detected. PpDELLAs likely act as 'hubs' for transcriptional regulation similarly to their homologues across the plant kingdom. Taken together, these data demonstrate that PpDELLA proteins share some biological functions with DELLAs in flowering plants, but other DELLA functions and regulation evolved independently in both plant lineages.

Evolution of genome space occupation in ferns: linking genome diversity and species richness

BACKGROUND AND AIMS

The dynamics of genome evolution caused by whole genome duplications and other processes are hypothesized to shape the diversification of plants and thus contribute to the astonishing variation in species richness among the main lineages of land plants. Ferns, the second most species-rich lineage of land plants, are highly suitable to test this hypothesis because of several unique features that distinguish fern genomes from those of seed plants. In this study, we tested the hypothesis that genome diversity and disparity shape fern species diversity by recording several parameters related to genome size and chromosome number.

METHODS

We conducted de novo measurement of DNA C-values across the fern phylogeny to reconstruct the phylogenetic history of the genome space occupation in ferns by integrating genomic parameters such as genome size, chromosome number and average DNA amount per chromosome into a time-scaled phylogenetic framework. Using phylogenetic generalized least square methods, we determined correlations between chromosome number and genome size, species diversity and evolutionary rates of their transformation.

KEY RESULTS

The measurements of DNA C-values for 233 species more than doubled the taxon coverage from ~2.2 % in previous studies to 5.3 % of extant diversity. The dataset not only documented substantial differences in the accumulation of genomic diversity and disparity among the major lineages of ferns but also supported the predicted correlation between species diversity and the dynamics of genome evolution.

CONCLUSIONS

Our results demonstrated substantial genome disparity among different groups of ferns and supported the prediction that alterations of reproductive modes alter trends of genome evolution. Finally, we recovered evidence for a close link between the dynamics of genome evolution and species diversity in ferns for the first time.

Revisiting ancient polyploidy in leptosporangiate ferns

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building K_S -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.

Giant Fern Genomes Show Complex Evolution Patterns: A Comparative Analysis in Two Species of Tmesipteris (Psilotaceae)

Giant genomes are rare across the plant kingdom and their study has focused almost exclusively on angiosperms and gymnosperms. The scarce genetic data that are available for ferns, however, indicate differences in their genome organization and a lower dynamism compared to other plant groups. Tmesipteris is a small genus of mainly epiphytic ferns that occur in Oceania and several Pacific Islands. So far, only two species with giant genomes have been reported in the genus, T. tannensis (1C = 73.19 Gbp) and T. obliqua (1C = 147.29 Gbp). Low-coverage genome skimming sequence data were generated in these two species and analyzed using the RepeatExplorer2 pipeline to identify and quantify the repetitive DNA fraction of these genomes. We found that both species share a similar genomic composition, with high repeat diversity compared to taxa with small (1C < 10 Gbp) genomes. We also found that, in general, characterized repetitive elements have relatively high heterogeneity scores, indicating ancient diverging evolutionary trajectories. Our results suggest that a whole genome multiplication event, accumulation of repetitive elements, and recent activation of those repeats have all played a role in shaping these genomes. It will be informative to compare these data in the future with data from the giant genome of the angiosperm Paris japonica, to determine if the structures observed here are an emergent property of massive genomic inflation or derived from lineage specific processes.

Cell Division and Meristem Dynamics in Fern Gametophytes

One of the most important questions in all multicellular organisms is how to define and maintain different cell fates during continuous cell division and proliferation. Plant meristems provide a unique research system to address this fundamental question because meristems dynamically maintain themselves and sustain organogenesis through balancing cell division and cell differentiation. Different from the gametophytes of seed plants that depend on their sporophytes and lack meristems, the gametophytes of seed-free ferns develop different types of meristems (including apical cell-based meristems and multicellular apical and marginal meristems) to promote independent growth and proliferation during the sexual gametophyte phase. Recent studies combining confocal time-lapse imaging and computational image analysis reveal the cellular basis of the initiation and proliferation of different types of meristems in fern gametophytes, providing new insights into the evolution of meristems in land plants. In this review, we summarize the recent progress in understanding the cell growth dynamics in fern gametophytes and discuss both conserved and diversified mechanisms underlying meristem cell proliferation in seed-free vascular plants.

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

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.

Deep origin and gradual evolution of transporting tissues: Perspectives from across the land plants

The evolution of transporting tissues was an important innovation in terrestrial plants that allowed them to adapt to almost all nonaquatic environments. These tissues consist of water-conducting cells and food-conducting cells and bridge plant-soil and plant-air interfaces over long distances. The largest group of land plants, representing about 95% of all known plant species, is associated with morphologically complex transporting tissue in plants with a range of additional traits. Therefore, this entire clade was named tracheophytes, or vascular plants. However, some nonvascular plants possess conductive tissues that closely resemble vascular tissue in their organization, structure, and function. Recent molecular studies also point to a highly conserved toolbox of molecular regulators for transporting tissues. Here, we reflect on the distinguishing features of conductive and vascular tissues and their evolutionary history. Rather than sudden emergence of complex, vascular tissues, plant transporting tissues likely evolved gradually, building on pre-existing developmental mechanisms and genetic components. Improved knowledge of the intimate structure and developmental regulation of transporting tissues across the entire taxonomic breadth of extant plant lineages, combined with more comprehensive documentation of the fossil record of transporting tissues, is required for a full understanding of the evolutionary trajectory of transporting tissues.

Transcriptional analysis of Ceratopteris richardii young sporophyte reveals conservation of stem cell factors in the root apical meristem

Gene expression in roots has been assessed in different plant species in studies ranging from complete organs to specific cell layers, and more recently at the single cell level. While certain genes or functional categories are expressed in the root of all or most plant species, lineage-specific genes have also been discovered. An increasing amount of transcriptomic data is available for angiosperms, while a limited amount of data is available for ferns, and few studies have focused on fern roots. Here, we present a de novo transcriptome assembly from three different parts of the Ceratopteris richardii young sporophyte. Differential gene expression analysis of the root tip transcriptional program showed an enrichment of functional categories related to histogenesis and cell division, indicating an active apical meristem. Analysis of a diverse set of orthologous genes revealed conserved expression in the root meristem, suggesting a preserved role for different developmental roles in this tissue, including stem cell maintenance. The reconstruction of evolutionary trajectories for ground tissue specification genes suggests a high degree of conservation in vascular plants, but not for genes involved in root cap development, showing that certain genes are absent in Ceratopteris or have intricate evolutionary paths difficult to track. Overall, our results suggest different processes of conservation and divergence of genes involved in root development.

Genomic insights into genetic diploidization in the homosporous fern Adiantum nelumboides

Whole genome duplication has been recognized as a major process in speciation of land plants, especially in ferns. Whereas genome downsizing contributes greatly to the post-genome shock responses of polyploid flowering plants, diploidization of polyploid ferns diverges by maintaining most of the duplicated DNA and is thus expected to be dominated by genic processes. As a consequence, fern genomes provide excellent opportunities to study ecological speciation enforced by expansion of protein families via polyploidy. To test the key predictions of this hypothesis, we reported the de novo genome sequence of Adiantum nelumboides, a tetraploid homosporous fern. The obtained draft genome had a size of 6.27 Gb assembled into 11,767 scaffolds with the contig N50 of 1.37 Mb. Repetitive DNA sequences contributed with about 81.7%, a remarkably high proportion of the genome. With 69,568 the number of predicted protein-coding genes exceeded those reported in most other land plant genomes. Intragenomic synteny analyses recovered 443 blocks with the average block size of 1.29 Mb and the average gene content of 16 genes. The results are consistent with the hypothesis of high ancestral chromosome number, lack of substantial genome downsizing, and dominance of genic diploidization. As expected in the calciphilous plants, a notable number of detected genes were involved in calcium uptake and transport. In summary, the genome sequence of a tetraploid homosporous fern not only provides access to a genomic resource of a derived fern, but also supports the hypothesis of maintenance of high chromosome numbers and duplicated DNA in young polyploid ferns.

Structure, mechanism, and phylogeny of LysM-chitinase conjugates specifically found in fern plants

A unique GH18 chitinase containing two N-terminal lysin motifs (PrLysM1 and PrLysM2) was first found in fern, Pteris ryukyuensis (Onaga and Taira, Glycobiology, 18, 414-423, 2008). This type of LysM-chitinase conjugates is not usually found in plants but in fungi. Here, we produced a similar GH18 chitinase with one N-terminal LysM module (EaLysM) from the fern, Equisetum arvense (EaChiA, Inamine et al., Biosci. Biotechnol. Biochem., 79, 1296-1304, 2015), using an Escherichia coli expression system and characterized for its structure and mechanism of action. The crystal structure of EaLysM exhibited an almost identical fold (βααβ) to that of PrLysM2. From isothermal titration calorimetry and nuclear magnetic resonance, the binding mode and affinities of EaLysM for chitooligosaccharides (GlcNAc)_n (3, 4, 5, and 6) were found to be comparable to those of PrLysM2. The LysM module in EaChiA is likely to bind (GlcNAc)_n almost independently through CH-π stacking of a Tyr residue with the pyranose ring. The (GlcNAc)_n-binding mode of LysMs in the LysM-chitinase conjugates from fern plants appears to differ from that of plant LysMs acting in chitin- or Nod-signal perception, in which multiple LysMs cooperatively act on (GlcNAc)_n. Phylogenetic analysis suggested that LysM-GH18 conjugates of fern plants formed a monophyletic group and had been separated earlier than forming the clade of fungal chitinases with LysMs.

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.

Positional cues and cell division dynamics drive meristem development and archegonium formation in Ceratopteris gametophytes

Fern gametophytes are autotrophic and independent of sporophytes, and they develop pluripotent meristems that drive prothallus development and sexual reproduction. To reveal cellular dynamics during meristem development in fern gametophytes, we performed long-term time-lapse imaging and determined the real-time lineage, identity and division activity of each single cell from meristem initiation to establishment in gametophytes of the fern Ceratopteris richardii. Our results demonstrate that in Ceratopteris gametophytes, only a few cell lineages originated from the marginal layer contribute to meristem initiation and proliferation, and the meristem lacks a distinguishable central zone or apical cell with low division activity. Within the meristem, cell division is independent of cell lineages and cells at the marginal layer are more actively dividing than inner cells. Furthermore, the meristem triggers differentiation of adjacent cells into egg-producing archegonia in a position-dependent manner. These findings advance the understanding of diversified meristem and gametophyte development in land plants.

Time-lapse imaging of the fern Ceratopteris richardii during meristem initiation and proliferation provides insights into the lineage, identity and division activity of each cell throughout the growth of gametophytes.

Cell growth dynamics in two types of apical meristems in fern gametophytes

In contrast to seed plants, the gametophytes of seed-free plants develop pluripotent meristems, which promote and sustain their independent growth and development. To date, the cellular basis of meristem development in gametophytes of seed-free ferns remains largely unknown. In this study, we used Woodsia obtusa, the blunt-lobe cliff fern, to quantitatively determine cell growth dynamics in two different types of apical meristems - the apical initial centered meristem and the multicellular apical meristem in gametophytes. Through confocal time-lapse live imaging and computational image analysis and quantification, we determined unique patterns of cell division and growth that sustain or terminate apical initials, dictate the transition from apical initials to multicellular apical meristems, and drive proliferation of apical meristems in ferns. Quantitative results showed that small cells correlated to active cell division in fern gametophytes. The marginal cells of multicellular apical meristems in fern gametophytes undergo division in both anticlinal and periclinal orientations, not only increasing cell numbers but also playing a dominant role in increasing cell layers during gametophyte development. All these findings provide insights into the function and regulation of meristems in gametophytes of seed-free vascular plants, suggesting both conserved and diversified mechanisms underlying meristem cell proliferation across land plants.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence

To date, little is known about the evolution of fern genomes, with only two small genomes published from the heterosporous Salviniales. Here we assembled the genome of Alsophila spinulosa, known as the flying spider-monkey tree fern, onto 69 pseudochromosomes. The remarkable preservation of synteny, despite resulting from an ancient whole-genome duplication over 100 million years ago, is unprecedented in plants and probably speaks to the uniqueness of tree ferns. Our detailed investigations into stem anatomy and lignin biosynthesis shed new light on the evolution of stem formation in tree ferns. We identified a phenolic compound, alsophilin, that is abundant in xylem, and we provided the molecular basis for its biosynthesis. Finally, analysis of demographic history revealed two genetic bottlenecks, resulting in rapid demographic declines of A. spinulosa. The A. spinulosa genome fills a crucial gap in the plant genomic landscape and helps elucidate many unique aspects of tree fern biology.

The biology of C. richardii as a tool to understand plant evolution

The fern Ceratopteris richardii has been studied as a model organism for over 50 years because it is easy to grow and has a short life cycle. In particular, as the first homosporous vascular plant for which genomic resources were developed, C. richardii has been an important system for studying plant evolution. However, we know relatively little about the natural history of C. richardii. In this article, we summarize what is known about this aspect of C. richardii, and discuss how learning more about its natural history could greatly increase our understanding of the evolution of land plants.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

Why Do Heterosporous Plants Have So Few Chromosomes?

The mechanisms controlling chromosome number, size, and shape, and the relationship of these traits to genome size, remain some of the least understood aspects of genome evolution. Across vascular plants, there is a striking disparity in chromosome number between homosporous and heterosporous lineages. Homosporous plants (comprising most ferns and some lycophytes) have high chromosome numbers compared to heterosporous lineages (some ferns and lycophytes and all seed plants). Many studies have investigated why homosporous plants have so many chromosomes. However, homospory is the ancestral condition from which heterospory has been derived several times. Following this phylogenetic perspective, a more appropriate question to ask is why heterosporous plants have so few chromosomes. Here, we review life history differences between heterosporous and homosporous plants, previous work on chromosome number and genome size in each lineage, known mechanisms of genome downsizing and chromosomal rearrangements, and conclude with future prospects for comparative research.

First Multi-Organ Full-Length Transcriptome of Tree Fern Alsophila spinulosa Highlights the Stress-Resistant and Light-Adapted Genes

Alsophila spinulosa, a relict tree fern, is a valuable plant for investigating environmental adaptations. Its genetic resources, however, are scarce. We used the PacBio and Illumina platforms to sequence the polyadenylated RNA of A. spinulosa root, rachis, and pinna, yielding 125,758, 89,107, and 89,332 unigenes, respectively. Combining the unigenes from three organs yielded a non-redundant reference transcriptome with 278,357 unigenes and N50 of 4141 bp, which were further reconstructed into 38,470 UniTransModels. According to functional annotation, pentatricopeptide repeat genes and retrotransposon-encoded polyprotein genes are the most abundant unigenes. Clean reads mapping to the full-length transcriptome is used to assess the expression of unigenes. The stress-induced ASR genes are highly expressed in all three organs, which is validated by qRT-PCR. The organ-specific upregulated genes are enriched for pathways involved in stress response, secondary metabolites, and photosynthesis. Genes for five types of photoreceptors, CRY signaling pathway, ABA biosynthesis and transduction pathway, and stomatal movement-related ion channel/transporter are profiled using the high-quality unigenes. The gene expression pattern coincides with the previously identified stomatal characteristics of fern. This study is the first multi-organ full-length transcriptome report of a tree fern species, the abundant genetic resources and comprehensive analysis of A. spinulosa, which provides the groundwork for future tree fern research.

Green plant genomes: What we know in an era of rapidly expanding opportunities

Green plants play a fundamental role in ecosystems, human health, and agriculture. As de novo genomes are being generated for all known eukaryotic species as advocated by the Earth BioGenome Project, increasing genomic information on green land plants is essential. However, setting standards for the generation and storage of the complex set of genomes that characterize the green lineage of life is a major challenge for plant scientists. Such standards will need to accommodate the immense variation in green plant genome size, transposable element content, and structural complexity while enabling research into the molecular and evolutionary processes that have resulted in this enormous genomic variation. Here we provide an overview and assessment of the current state of knowledge of green plant genomes. To date fewer than 300 complete chromosome-scale genome assemblies representing fewer than 900 species have been generated across the estimated 450,000 to 500,000 species in the green plant clade. These genomes range in size from 12 Mb to 27.6 Gb and are biased toward agricultural crops with large branches of the green tree of life untouched by genomic-scale sequencing. Locating suitable tissue samples of most species of plants, especially those taxa from extreme environments, remains one of the biggest hurdles to increasing our genomic inventory. Furthermore, the annotation of plant genomes is at present undergoing intensive improvement. It is our hope that this fresh overview will help in the development of genomic quality standards for a cohesive and meaningful synthesis of green plant genomes as we scale up for the future.

Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions

Third-generation long-read sequencing is transforming plant genomics. Oxford Nanopore Technologies and Pacific Biosciences are offering competing long-read sequencing technologies and enable plant scientists to investigate even large and complex plant genomes. Sequencing projects can be conducted by single research groups and sequences of smaller plant genomes can be completed within days. This also resulted in an increased investigation of genomes from multiple species in large scale to address fundamental questions associated with the origin and evolution of land plants. Increased accessibility of sequencing devices and user-friendly software allows more researchers to get involved in genomics. Current challenges are accurately resolving diploid or polyploid genome sequences and better accounting for the intra-specific diversity by switching from the use of single reference genome sequences to a pangenome graph.

Quantitative live-imaging reveals the dynamics of apical cells during gametophyte development in ferns

Meristems in land plants share conserved functions but develop highly variable structures. Meristems in seed-free plants, including ferns, usually contain one or a few pyramid-/wedge-shaped apical cells (ACs) as initials, which are lacking in seed plants. It remained unclear how ACs promote cell proliferation in fern gametophytes and whether any persistent AC exists to sustain fern gametophyte development continuously. Here, we uncovered previously undefined ACs maintained even at late developmental stages in fern gametophytes. Through quantitative live-imaging, we determined division patterns and growth dynamics that maintain the persistent AC in Sphenomeris chinensis, a representative fern. The AC and its immediate progenies form a conserved cell packet, driving cell proliferation and prothallus expansion. At the apical centre of gametophytes, the AC and its adjacent progenies display small dimensions resulting from active cell division instead of reduced cell expansion. These findings provide insight into diversified meristem development in land plants.

The evolution of imprinting in plants: beyond the seed

Genomic imprinting results in the biased expression of alleles depending on if the allele was inherited from the mother or the father. Despite the prevalence of sexual reproduction across eukaryotes, imprinting is only found in placental mammals, flowering plants, and some insects, suggesting independent evolutionary origins. Numerous hypotheses have been proposed to explain the selective pressures that favour the innovation of imprinted gene expression and each differs in their experimental support and predictions. Due to the lack of investigation of imprinting in land plants, other than angiosperms with triploid endosperm, we do not know whether imprinting occurs in species lacking endosperm and with embryos developing on maternal plants. Here, we discuss the potential for uncovering additional examples of imprinting in land plants and how these observations may provide additional support for one or more existing imprinting hypotheses.

Timing of meristem initiation and maintenance determines the morphology of fern gametophytes

The alternation of generations in land plants occurs between the sporophyte phase and the gametophyte phase. The sporophytes of seed plants develop self-maintained, multicellular meristems, and these meristems determine plant architecture. The gametophytes of seed plants lack meristems and are heterotrophic. In contrast, the gametophytes of seed-free vascular plants, including ferns, are autotrophic and free-living, developing meristems to sustain their independent growth and proliferation. Compared with meristems in the sporophytes of seed plants, the cellular mechanisms underlying meristem development in fern gametophytes remain largely unknown. Here, using confocal time-lapse live imaging and computational segmentation and quantification, we determined different patterns of cell divisions associated with the initiation and proliferation of two distinct types of meristems in gametophytes of two closely related Pteridaceae ferns, Pteris vittata and Ceratopteris richardii. Our results reveal how the simple timing of a switch between two meristems has considerable consequences for the divergent gametophyte morphologies of the two ferns. They further provide evolutionary insight into the function and regulation of gametophyte meristems in seed-free vascular plants.

Conditional stomatal closure in a fern shares molecular features with flowering plant active stomatal responses

Stomata evolved as plants transitioned from water to land, enabling carbon dioxide uptake and water loss to be controlled. In flowering plants, the most recently divergent land plant lineage, stomatal pores actively close in response to drought. In this response, the phytohormone abscisic acid (ABA) triggers signaling cascades that lead to ion and water loss in the guard cells of the stomatal complex, causing a reduction in turgor and pore closure. Whether this stimulus-response coupling pathway acts in other major land plant lineages is unclear, with some investigations reporting that stomatal closure involves ABA but others concluding that closure is passive. Here, we show that in the model fern Ceratopteris richardii active stomatal closure is conditional on sensitization by pre-exposure to either low humidity or exogenous ABA and is promoted by ABA. RNA-seq analysis and de novo transcriptome assembly reconstructed the protein-coding complement of the C. richardii genome, with coverage comparable to other plant models, enabling transcriptional signatures of stomatal sensitization and closure to be inferred. In both cases, changes in abundance of homologs of ABA, Ca²⁺, and ROS-related signaling components were observed, suggesting that the closure-response pathway is conserved in ferns and flowering plants. These signatures further suggested that sensitization is achieved by lowering the threshold required for a subsequent closure-inducing signal to trigger a response. We conclude that the canonical signaling network for active stomatal closure functioned in at least a rudimentary form in the stomata of the last common ancestor of ferns and flowering plants.

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.

De Novo Sporophyte Transcriptome Assembly and Functional Annotation in the Endangered Fern Species Vandenboschia speciosa (Willd.) G. Kunkel

We sequenced the sporophyte transcriptome of Killarney fern (Vandenboschia speciosa (Willd.) G. Kunkel). In addition to being a rare endangered Macaronesian-European endemism, this species has a huge genome (10.52 Gb) as well as particular biological features and extreme ecological requirements. These characteristics, together with the systematic position of ferns among vascular plants, make it of high interest for evolutionary, conservation and functional genomics studies. The transcriptome was constructed de novo and contained 36,430 transcripts, of which 17,706 had valid BLAST hits. A total of 19,539 transcripts showed at least one of the 7362 GO terms assigned to the transcriptome, whereas 6547 transcripts showed at least one of the 1359 KEGG assigned terms. A prospective analysis of functional annotation results provided relevant insights on genes involved in important functions such as growth and development as well as physiological adaptations. In this context, a catalogue of genes involved in the genetic control of plant development, during the vegetative to reproductive transition, in stress response as well as genes coding for transcription factors is given. Altogether, this study provides a first step towards understanding the gene expression of a significant fern species and the in silico functional and comparative analyses reported here provide important data and insights for further comparative evolutionary studies in ferns and land plants in general.

Bioactive compounds and their libraries: An insight into prospective phytotherapeutics approach for oral mucocutaneous cancers

Oral mucocutaneous cancers (OMCs) are cancers that affect both the oral mucosa and perioral cutaneous structures. Common OMCs are squamous cell carcinoma (SCC), basal cell carcinoma (BCC) and malignant melanoma (MM). Anatomical similarities and conventions which categorizes these lesions blur the magnitude of OMCs in diverse populations. The burden of OMC is high in the sub-Saharan Africa and Indian subcontinents, and the cost of management is prohibitive in the resource-limited, developing world. Hence, there is a pressing demand for the use of cost-effective in silico approaches to identify diagnostic tools and treatment targets for diseases with high burdens in these regions. Due to their ubiquitousness and accessibility, the use of therapeutic efficacy of plant bioactive compounds in the management of OMC is both appropriate and plausible. Furthermore, screening known mechanistic disease targets with well annotated plant bioactive compound libraries is poised to improve the routine management of OMCs provided that the requisite access to database resources are available and accessible. Using natural products minimizes the side effects and morbidities associated with conventional therapies. The development of innovative treatments approaches would tremendously benefit the African and Indian populace and reduce the mortalities associated with OMCs in the developing world. Hence, we discuss herein, the potential benefits, opportunities and challenges of using bioactive compound libraries in the management of OMCs.

Genome-wide comparative analyses of GATA transcription factors among 19 Arabidopsis ecotype genomes: Intraspecific characteristics of GATA transcription factors

GATA transcription factors (TFs) are widespread eukaryotic regulators whose DNA-binding domain is a class IV zinc finger motif (CX2CX17-20CX2C) followed by a basic region. Due to the low cost of genome sequencing, multiple strains of specific species have been sequenced: e.g., number of plant genomes in the Plant Genome Database (http://www.plantgenome.info/) is 2,174 originated from 713 plant species. Thus, we investigated GATA TFs of 19 Arabidopsis thaliana genome-widely to understand intraspecific features of Arabidopsis GATA TFs with the pipeline of GATA database (http://gata.genefamily.info/). Numbers of GATA genes and GATA TFs of each A. thaliana genome range from 29 to 30 and from 39 to 42, respectively. Four cases of different pattern of alternative splicing forms of GATA genes among 19 A. thaliana genomes are identified. 22 of 2,195 amino acids (1.002%) from the alignment of GATA domain amino acid sequences display variations across 19 ecotype genomes. In addition, maximally four different amino acid sequences per each GATA domain identified in this study indicate that these position-specific amino acid variations may invoke intraspecific functional variations. Among 15 functionally characterized GATA genes, only five GATA genes display variations of amino acids across ecotypes of A. thaliana, implying variations of their biological roles across natural isolates of A. thaliana. PCA results from 28 characteristics of GATA genes display the four groups, same to those defined by the number of GATA genes. Topologies of bootstrapped phylogenetic trees of Arabidopsis chloroplasts and common GATA genes are mostly incongruent. Moreover, no relationship between geographical distribution and their phylogenetic relationships was found. Our results present that intraspecific variations of GATA TFs in A. thaliana are conserved and evolutionarily neutral along with 19 ecotypes, which is congruent to the fact that GATA TFs are one of the main regulators for controlling essential mechanisms, such as seed germination and hypocotyl elongation.

Charting the genomic landscape of seed-free plants

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

Polyploid phylogenetics

Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.

Plant evolution driven by interactions with symbiotic and pathogenic microbes

During 450 million years of diversification on land, plants and microbes have evolved together. This is reflected in today's continuum of associations, ranging from parasitism to mutualism. Through phylogenetics, cell biology, and reverse genetics extending beyond flowering plants into bryophytes, scientists have started to unravel the genetic basis and evolutionary trajectories of plant-microbe associations. Protection against pathogens and support of beneficial, symbiotic, microorganisms are sustained by a blend of conserved and clade-specific plant mechanisms evolving at different speeds. We propose that symbiosis consistently emerges from the co-option of protection mechanisms and general cell biology principles. Exploring and harnessing the diversity of molecular mechanisms used in nonflowering plant-microbe interactions may extend the possibilities for engineering symbiosis-competent and pathogen-resilient crops.

A De Novo Transcriptome Assembly of Ceratopteris richardii Provides Insights into the Evolutionary Dynamics of Complex Gene Families in Land Plants

As the closest extant sister group to seed plants, ferns are an important reference point to study the origin and evolution of plant genes and traits. One bottleneck to the use of ferns in phylogenetic and genetic studies is the fact that genome-level sequence information of this group is limited, due to the extreme genome sizes of most ferns. Ceratopteris richardii (hereafter Ceratopteris) has been widely used as a model system for ferns. In this study, we generated a transcriptome of Ceratopteris, through the de novo assembly of the RNA-seq data from 17 sequencing libraries that are derived from two sexual types of gametophytes and five different sporophyte tissues. The Ceratopteris transcriptome, together with 38 genomes and transcriptomes from other species across the Viridiplantae, were used to uncover the evolutionary dynamics of orthogroups (predicted gene families using OrthoFinder) within the euphyllophytes and identify proteins associated with the major shifts in plant morphology and physiology that occurred in the last common ancestors of euphyllophytes, ferns, and seed plants. Furthermore, this resource was used to identify and classify the GRAS domain transcriptional regulators of many developmental processes in plants. Through the phylogenetic analysis within each of the 15 GRAS orthogroups, we uncovered which GRAS family members are conserved or have diversified in ferns and seed plants. Taken together, the transcriptome database and analyses reported here provide an important platform for exploring the evolution of gene families in land plants and for studying gene function in seed-free vascular plants.

Origin and diversification of ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) genes in land plants and phylogenetic evidence that the ancestral CER1/3 gene resulted from the fusion of pre-existing domains

ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) are key genes in synthesis of alkanes, a major component of cuticular waxes in land plants. The genes share extensive similarity, including the N-terminal (ERG3/FAH) and C-terminal (WAX2) domains. This study, traces the origin, evolutionary history, phylogenetic relationships and variation in copy number of the two genes within and beyond the Viridiplantae (green plants). Protein homologs of both CER1 and CER3 were identified across most Embryophyta (land plants), a single homolog (CER1/3) in charophytes and prasinophytes, and none in the other green, red or brown algae. Ancestral state reconstructions in 100 sequenced Archaeplastida using presence/absence of CER1/3 family genes revealed that the CER1/3 gene probably originated in the common ancestor of Viridiplantae. Phylogenetic analysis of CER1 and CER3 protein sequences from 146 plant species strongly suggests that the two genes originated by duplication of CER1/3 in the ancestral embryophyte. The evolution of CER1 and CER3 genes involved differential divergence of the two domains. Outside Embryophyta, CER1/3 similar sequences identified in diatoms and a cryptophyte, were the closest relatives of the CER1/3 family proteins. Proteins harbouring WAX2-wxAR (WAX2 associated region) similar regions were identified in proteins of bacteria, Archaea, cryptophytes, dinoflagellates and Stramenopiles. The independent existence of both ERG3/FAH and WAX2-wxAR domains in diverse lineages strongly points to the origin of CER1/3 gene in green plants by the fusion of pre-existing domains.

HAM Gene Family and Shoot Meristem Development

Land plants develop highly diversified shoot architectures, all of which are derived from the pluripotent stem cells in shoot apical meristems (SAMs). As sustainable resources for continuous organ formation in the aboveground tissues, SAMs play an important role in determining plant yield and biomass production. In this review, we summarize recent advances in understanding one group of key regulators - the HAIRY MERISTEM (HAM) family GRAS domain proteins - in shoot meristems. We highlight the functions of HAM family members in dictating shoot stem cell initiation and proliferation, the signaling cascade that shapes HAM expression domains in shoot meristems, and the conservation and diversification of HAM family members in land plants. We also discuss future directions that potentially lead to a more comprehensive view of the HAM gene family and stem cell homeostasis in land plants.

Development and Cell Cycle Activity of the Root Apical Meristem in the Fern Ceratopteris richardii

Ferns are a representative clade in plant evolution although underestimated in the genomic era. Ceratopteris richardii is an emergent model for developmental processes in ferns, yet a complete scheme of the different growth stages is necessary. Here, we present a developmental analysis, at the tissue and cellular levels, of the first shoot-borne root of Ceratopteris. We followed early stages and emergence of the root meristem in sporelings. While assessing root growth, the first shoot-borne root ceases its elongation between the emergence of the fifth and sixth roots, suggesting Ceratopteris roots follow a determinate developmental program. We report cell division frequencies in the stem cell niche after detecting labeled nuclei in the root apical cell (RAC) and derivatives after 8 h of exposure. These results demonstrate the RAC has a continuous mitotic activity during root development. Detection of cell cycle activity in the RAC at early times suggests this cell acts as a non-quiescent organizing center. Overall, our results provide a framework to study root function and development in ferns and to better understand the evolutionary history of this organ.

Light- and hormone-mediated development in non-flowering plants: An overview

Light, hormones and their interaction regulate different aspects of development in non-flowering plants. They might have played a role in the evolution of different plant groups by conferring specific adaptive evolutionary changes. Plants are sessile organisms. Unlike animals, they lack the opportunity to abandon their habitat in unfavorable conditions. They respond to different environmental cues and adapt accordingly to control their growth and developmental pattern. While phytohormones are known to be internal regulators of plant development, light is a major environmental signal that shapes plant processes. It is plausible that light-hormone crosstalk might have played an important role in plant evolution. But how the crosstalk between light and phytohormone signaling pathways might have shaped the plant evolution is unclear. One of the possible reasons is that flowering plants have been studied extensively in context of plant development, which cannot serve the purpose of evolutionary comparisons. In order to elucidate the role of light, hormone and their crosstalk in the evolutionary adaptation in plant kingdom, one needs to understand various light- and hormone-mediated processes in diverse non-flowering plants. This review is an attempt to outline major light- and phytohormone-mediated responses in non-flowering plant groups such as algae, bryophytes, pteridophytes and gymnosperms.

Leaf dorsoventrality candidate gene CpARF4 has conserved expression pattern but divergent tasiR-ARF regulation in the water fern Ceratopteris pteridoides

PREMISE

Leaves are traditionally classified into microphylls and megaphylls, and recently have been regarded as independently originating in lycophytes, ferns, and seed plants. The developmental genetics of leaf dorsoventrality, a synapomorphy in vascular plants, has been extensively studied in flowering plants. AUXIN RESPONSE FACTOR4 (ARF4) genes are key to leaf abaxial identity in flowering plants, but whether they exist in ferns is still an open question.

METHODS

ARF4 genes from Ceratopteris pteridoides, Cyrtomium guizhouense, and Parathelypteris nipponica were mined from transcriptomes and investigated in terms of evolutionary phylogeny and sequence motifs, with a focus on the tasiR-ARF binding site. In situ hybridization was used to localize expression of CpARF4 in Ceratopteris pteridoides. 5'RNA ligase-mediated-RACE was employed to verify whether CpARF4 transcripts were sliced by tasiR-ARF.

RESULTS

ARF4 genes exist in ferns, and this lineage originates from a gene duplication in the common ancestor of ferns and seed plants. ARF4 genes are of a single copy in the ferns studied here, and they contain divergent and, at most, one tasiR-ARF binding site. CpARF4 is expressed in the abaxial but not the adaxial domain of leaf primordia at various developmental stages. Transcript slicing guided by tasiR-ARF is active in C. pteridoides, but CpARF4 probably has not been affected by it.

CONCLUSIONS

Fern ARF4 genes differ in copy number and tasiR-ARF regulation relative to flowering plants, though they can be similarly expressed in the abaxial domain of leaves, revealing a key role for ARF4 genes in the evolution of leaf dorsoventrality of vascular plants.

Delving into the evolutionary origin of steroid sensing in plants

Brassinosteroids (BRs) are steroid hormones that play crucial roles in plant growth, development and adaptation to shifting environmental conditions. Our current understanding of the origin, evolution and functional significance of BRs is influenced by a double-edged bias: most we know stems from studies on a single species and, on the flip side, dearth of information from a phylogenetically broad and significant array of land plants precludes well-grounded comparisons. Here, we provide an update on BR presence and sensing along land plant evolution. Furthermore, a comprehensive search in all major plant lineages reveals the widespread presence of BR-receptor related sequences, suggesting that steroid-related signals may have been functional early in the evolution of land plants.

What can lycophytes teach us about plant evolution and development? Modern perspectives on an ancient lineage

All Evo-Devo studies rely on representative sampling across the tree of interest to elucidate evolutionary trajectories through time. In land plants, genetic resources are well established in model species representing lineages including bryophytes (mosses, liverworts, and hornworts), monilophytes (ferns and allies), and seed plants (gymnosperms and flowering plants), but few resources are available for lycophytes (club mosses, spike mosses, and quillworts). Living lycophytes are a sister group to the euphyllophytes (the fern and seed plant clade), and have retained several ancestral morphological traits despite divergence from a common ancestor of vascular plants around 420 million years ago. This sister relationship offers a unique opportunity to study the conservation of traits such as sporophyte branching, vasculature, and indeterminacy, as well as the convergent evolution of traits such as leaves and roots which have evolved independently in each vascular plant lineage. To elucidate the evolution of vascular development and leaf formation, molecular studies using RNA Seq, quantitative reverse transcription polymerase chain reaction, in situ hybridisation and phylogenetics have revealed the diversification and expression patterns of KNOX, ARP, HD-ZIP, KANADI, and WOX gene families in lycophytes. However, the molecular basis of further trait evolution is not known. Here we describe morphological traits of living lycophytes and their extinct relatives, consider the molecular underpinnings of trait evolution and discuss future research required in lycophytes to understand the key evolutionary innovations enabling the growth and development of all vascular plants.

Cryptic diversity in the model fern genus Ceratopteris (Pteridaceae)

Cryptic species are present throughout the tree of life. They are especially prevalent in ferns, because of processes such hybridization, polyploidy, and reticulate evolution. In addition, the simple morphology of ferns limits phenotypic variation and makes it difficult to detect cryptic species. The model fern genus Ceratopteris has long been suspected to harbor cryptic diversity, in particular within the highly polymorphic C. thalictroides. Yet no studies have included samples from throughout its pan-tropical range or utilized genomic sequencing, making it difficult to assess the full extent of cryptic variation within this genus. Here, we present the first multilocus phylogeny of the genus using reduced representation genomic sequencing (RADseq) and examine population structure, phylogenetic relationships, and ploidy level variation. We recover similar species relationships found in previous studies, find support for the cryptic species C. gaudichaudii as genetically distinct, and identify novel genomic variation within two of the mostly broadly distributed species in the genus, C. thalictroides and C. cornuta. Finally, we detail the utility of our approach for working on cryptic, reticulate groups of ferns. Specifically, it does not require a reference genome, of which there are very few available for ferns. RADseq is a cost-effective way to work with study groups lacking genomic resources, and to obtain the thousands of nuclear markers needed to untangle species complexes.