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

Comparative transcriptomics in ferns reveals key innovations and divergent evolution of the secondary cell walls

Ferns are essential for understanding plant evolution; however, their large and intricate genomes have kept their genetic landscape largely unexplored, with only a few genomes sequenced and limited transcriptomic data available. To bridge this gap, we generated extensive RNA-sequencing data across various organs from 22 representative fern species, resulting in high-quality transcriptome assemblies. These data enabled us to construct a time-calibrated phylogeny for ferns, encompassing all major clades, which revealed numerous instances of whole-genome duplication. We highlighted the distinctiveness of fern genetics, discovering that half of the identified gene families are unique to ferns. Our exploration of fern cell walls through biochemical and immunological analyses uncovered the presence of the lignin syringyl unit, along with evidence of its independent evolution in ferns. Additionally, the identification of an unusual sugar in fern cell walls suggests a divergent evolutionary trajectory in cell wall biochemistry, probably influenced by gene duplication and sub-functionalization. To facilitate further research, we have developed an online database that includes preloaded genomic and transcriptomic data for ferns and other land plants. We used this database to demonstrate the independent evolution of lignocellulosic gene modules in ferns. Our findings provide a comprehensive framework illustrating the unique evolutionary journey ferns have undertaken since diverging from the last common ancestor of euphyllophytes more than 360 million years ago.

Ancient large-scale gene duplications and diversification in bryophytes illuminate the plant terrestrialization

Large-scale gene duplications (LSGDs) are crucial for evolutionary adaptation and recurrent in vascular plants. However, the role of ancient LSGDs in the terrestrialization and diversification of bryophytes, the second most species-rich group of land plants, remains largely elusive due to limited sampling in bryophytes. Employing the most extensive nuclear gene dataset in bryophytes to date, we reconstructed a time-calibrated phylogenetic tree from 209 species, covering virtually all key bryophyte lineages, for phylogenomic analyses of LSGDs and diversification. We newly identified two ancient LSGDs: one in the most recent common ancestor (MRCA) of extant bryophytes and another in the MRCA of the majority of Jungermanniales s. lato. Duplicated genes from these two LSGDs show significant enrichment in photosynthesis-related processes and structures. Rhizoid-responsive ROOTHAIR DEFECTIVE SIX-LIKE (RSL) genes from ancient LSGDs are present in rhizoidless bryophytes, challenging assumptions about rhizoid absence mechanisms. We highlighted four major diversification rate upshifts, two of which slightly postdated LSGDs, potentially linked to the flourishing of gymnosperms and angiosperms and explaining over 80% of bryophyte diversity. Our findings, supported by extensive bryophyte sampling, highlight the significance of LSGDs in the early terrestrialization and diversification of bryophytes, offering new insights into land plant evolution.

Genomic correlates of vascular plant reproductive complexity and the uniqueness of angiosperms

Whole genome duplication (WGD) likely plays an important role in plant macroevolution, and has been implicated in diversification rate shifts, structural innovations, and increased disparity. But the general effects of WGD are challenging to evaluate, in part due to the difficulty of directly comparing morphological patterns across disparate clades. We explored relationships between WGD and the evolution of reproductive complexity across vascular plants using a metric based on the number of reproductive part types. We used multiple regression models to evaluate the relative importance of inferred WGD events, genome size, and a suite of additional variables relating to growth habit and reproductive biology in explaining part type complexity. WGD was a consistent predictor of reproductive complexity only among angiosperms. Across vascular plants generally, reproductive biology, clade identity, and the presence of bisexual strobili (those that produce microsporangiate and megasporangiate organs) were better predictors of complexity. Angiosperms are unique among vascular plants in combining frequent polyploidy with high-reproductive complexity. Whether WGD is mechanistically linked to floral complexity is unclear, but we suggest widespread polyploidy and increased complexity were ultimately facilitated by the evolution of herbaceous growth habits in early angiosperms.

Evolutionary dynamics of nitrate uptake, assimilation, and signalling in plants: adapting to a changing environment

Nitrogen (N) is a crucial macronutrient for plant growth, with nitrate as a primary inorganic N source for most plants. Beyond its role as a nutrient, nitrate also functions as a signalling molecule, influencing plant morphogenetic development. While nitrate utilization and signalling mechanisms have been extensively studied in model plants, the origin, evolution, and diversification of core components in nitrate uptake, assimilation, and signalling remain largely unexplored. In our investigation, we discovered that deep sea algae living in low nitrate conditions developed a high-affinity transport system (HATS) for nitrate uptake and a pathway of nitrate primary assimilation (NR-NiR-GS-GOGAT). In contrast, low-affinity transport systems (LATS) and the plastid GS originated from the ancestors of land and seed plants, respectively. These adaptations facilitated amino acid acquisition as plants conquered terrestrial environments. Furthermore, the intricate nitrate signalling, relying on NRT1.1 and NLP7, evolved stepwise, potentially establishing systematic regulation in bryophytes for self-regulation under complex terrestrial nitrate environments. As plants underwent terrestrialization, they underwent adaptive changes to thrive in dynamic nitrate environments, continually enhancing their nitrate uptake, assimilation, and signal transduction abilities.

Phylogenomics resolves the backbone of Poales and identifies signals of hybridization and polyploidy

Poales, as one of the largest orders of angiosperm, holds crucial economic and ecological importance. Nevertheless, achieving a consensus topology has been challenging in previous studies due to limited molecular data and sparse taxon sampling. The uneven distribution of species diversity among families and the factors leading to elevated species richness in certain lineages have also been subjects of ongoing discussion and investigation. In this study, we conducted a comprehensive sampling, including representatives from all 14 families and 85 taxa of Poales, along with five additional outgroups. To reconstruct the phylogeny of Poales, we employed a combination of coalescent and concatenation methods on three nuclear gene sets (1093, 491, 143) and one plastid gene set (53), which were inferenced from genomic data. We also conducted phylogenetic hypothesis analyses to evaluate two major conflicting nodes detected in phylogenetic analyses. As a result, we successfully resolved the backbone of Poales and provided a timeline for its evolutionary history. We recovered the sister relationship between Typhaceae and Bromeliaceae as the earliest diverging families within Poales. The clade consisting of Ecdeiocoleaceae and Joinvilleaceae was recovered as the sister group of Poaceae. Within the xyrid clade, Mayacaceae and Erioaculaceae + Xyridaceae successively diverged along the backbone of Poales. The topology of [Aristidoideae, ((Micrairoideae, Panicoideae), (Arundinoideae, (Chloridoideae, Danthonioideae)))] within the PACMAD clade has received strong support from multiple findings. We also delved into the underlying biological factors that contributed to the conflicting nodes observed in the phylogenetic analysis. Apart from the uncertainty regarding the sister group of Poaceae caused by cytonuclear discordance, frequent hybridization and polyploidy may have contributed to other conflicting nodes. We identified 26 putative whole-genome duplication (WGD) events within Poales. However, apart from the σ-WGD and the ρ-WGD, we did not observe any potential polyploid events that could be directly linked to the species diversification in specific lineages. Furthermore, there was a significant increase in the net diversification rate of Poales following the K-Pg boundary.

All tangled up: Unraveling phylogenetics and reticulate evolution in the vining ferns, Lygodium (Schizaeales)

PREMISE

Reticulate evolution, often accompanied by polyploidy, is prevalent in plants, and particularly in the ferns. Resolving the resulting non-bifurcating histories remains a major challenge for plant phylogenetics. Here, we present a phylogenomic investigation into the complex evolutionary history of the vining ferns, Lygodium (Lygodiaceae, Schizaeales).

METHODS

Using a targeted enrichment approach with the GoFlag 408 flagellate land plant probe set, we generated large nuclear and plastid sequence datasets for nearly all taxa in the genus and constructed the most comprehensive phylogeny of the family to date using concatenated maximum likelihood and coalescence approaches. We integrated this phylogeny with cytological and spore data to explore karyotype evolution and generate hypotheses about the origins of putative polyploids and hybrids.

RESULTS

Our data and analyses support the origins of several putative allopolyploids (e.g., L. cubense, L. heterodoxum) and hybrids (e.g., L. ×fayae) and also highlight the potential prevalence of autopolyploidy in this clade (e.g., L. articulatum, L. flexuosum, and L. longifolium).

CONCLUSIONS

Our robust phylogenetic framework provides valuable insights into dynamic reticulate evolution in this clade and demonstrates the utility of target-capture data for resolving these complex relationships.

The Possible Earliest Allopolyploidization in Tracheophytes Revealed by Phylotranscriptomics and Morphology of Selaginellaceae

Selaginellaceae, originated in the Carboniferous and survived the Permian-Triassic mass extinction, is the largest family of lycophyte, which is sister to other tracheophytes. It stands out from tracheophytes by exhibiting extraordinary habitat diversity and lacking polyploidization. The organelle genome-based phylogenies confirmed the monophyly of Selaginella, with six or seven subgenera grouped into two superclades, but the phylogenetic positions of the enigmatic Selaginella sanguinolenta clade remained problematic. Here, we conducted a phylogenomic study on Selaginellaceae utilizing large-scale nuclear gene data from RNA-seq to elucidate the phylogeny and explore the causes of the phylogenetic incongruence of the S. sanguinolenta clade. Our phylogenetic analyses resolved three different positions of the S. sanguinolenta clade, which were supported by the sorted three nuclear gene sets, respectively. The results from the gene flow test, species network inference, and plastome-based phylogeny congruently suggested a probable hybrid origin of the S. sanguinolenta clade involving each common ancestor of the two superclades in Selaginellaceae. The hybrid hypothesis is corroborated by the evidence from rhizophore morphology and spore micromorphology. The chromosome observation and Ks distributions further suggested hybridization accompanied by polyploidization. Divergence time estimation based on independent datasets from nuclear gene sets and plastid genome data congruently inferred that allopolyploidization occurred in the Early Triassic. To our best knowledge, the allopolyploidization in the Mesozoic reported here represents the earliest record of tracheophytes. Our study revealed a unique triad of phylogenetic positions for a hybrid-originated group with comprehensive evidence and proposed a hypothesis for retaining both parental alleles through gene conversion.

Unfurling an improved method for visualizing mitotic chromosomes in ferns

Premise

Cytotaxonomy employs chromosome visualization to study organismal relationships and evolution. Despite the critical value of cytogenetic data, cytotypes are lacking for many plant groups. Here, we present an improved approach for visualizing mitotic chromosomes in ferns, a key lineage of land plants, using the dividing cells of unfurling croziers (fiddleheads).

Methods and Results

Our modified mitotic chromosome preparation incorporates a brief pectinase-cellulase pretreatment, as well as colchicine fixation and the Feulgen reaction to improve the staining and separation of mitotic chromosomes. To demonstrate this easy and efficient assessment, we determined the sporophytic (2n) chromosome number for three fern species: Cheilanthes mollis (2n = 60), Cheilanthes hypoleuca (2n = 120), and Nephrolepis cordifolia (2n = 82).

Conclusions

The new method presented here improves visualizations of mitotic chromosomes from the dividing nuclei of young fern croziers. Fiddleheads are widely accessible in nature and in living collections worldwide, and this modified approach increases their suitability for fern cytotaxonomic studies.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Organellar phylogenomics of Ophioglossaceae fern genera

Previous phylogenies showed conflicting relationships among the subfamilies and genera within the fern family Ophioglossaceae. However, their classification remains unsettled where contrasting classifications recognize four to 15 genera. Since these treatments are mostly based on phylogenetic evidence using limited, plastid-only loci, a phylogenomic understanding is actually necessary to provide conclusive insight into the systematics of the genera. In this study, we have therefore compiled datasets with the broadest sampling of Ophioglossaceae genera to date, including all fifteen currently recognized genera, especially for the first time the South African endemic genus Rhizoglossum. Notably, our comprehensive phylogenomic matrix is based on both plastome and mitogenome genes. Inferred from the coding sequences of 83 plastid and 37 mitochondrial genes, a strongly supported topology for these subfamilies is presented, and is established by analyses using different partitioning approaches and substitution models. At the generic level, most relationships are well resolved except for few within the subfamily Ophioglossoideae. With this new phylogenomic scheme, key morphological and genomic changes were further identified along this backbone. In addition, we confirmed numerous horizontally transferred (HGT) genes in the genera Botrypus, Helminthostachys, Mankyua, Sahashia, and Sceptridium. These HGT genes are most likely located in mitogenomes and are predominately donated from angiosperm Santalales or non-Ophioglossaceae ferns. By our in-depth searches of the organellar genomes, we also provided phylogenetic overviews for the plastid and mitochondrial MORFFO genes found in these Ophioglossaceae ferns.

Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

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.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Phylogenomic evolutionary insights in the fern family Gleicheniaceae

The pantropical fern family Gleicheniaceae comprises approximately 157 species. Seven genera are currently recognized in the family, although their monophyly is still uncertain due to low sampling in phylogenetic studies. We examined the monophyly of the genera through extended sampling, using the first phylogenomic inference of the family including data from both nuclear and plastid genomes. Seventy-six samples were sequenced (70 Gleicheniaceae species and six outgroups) using high throughput sequencing, including all seven currently recognized genera. Plastid and nuclear data were recovered and assembled; the nuclear data was phased to reduce paralogy as well as hybrid noise in the final recovered topology. Maximum likelihood trees were built for each locus, and a concatenated dataset was built for both datasets. A species tree based on a multispecies coalescent model was generated, and divergence time analyses performed. We here present the first genomic phylogenetic inferences concerning Gleicheniaceae, confirming the monophyly of most genera except Sticherus, which we recovered as paraphyletic. Although most of the extant genera of Gleicheniaceae originated during the Mesozoic, several genera show Neogene and even Quaternary diversifications, and our results suggest that reticulation and polyploidy may have played significant roles during this diversification. However, some genera, such as Rouxopteris and Stromatopteris, appear to represent evolutionary relicts.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

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.

SHINE clade of ERF transcription factors: A significant player in abiotic and biotic stress tolerance in plants

SHINE (SHN) clade transcription factors (TFs) represents a subfamily of APETALA2/ethylene-responsive factor (AP2/ERF) proteins. The latter, is characterized by its responsiveness to the phytohormone ethylene and the presence of AP2 DNA-binding domain. They are involved in many biological processes and in responses to different environmental constraints. SHN TFs were among the first identified regulators of cuticle formation. Cuticle plays crucial role in plant tolerance to drought, salinity and high temperature as well as in defense against pathogens. In addition, SHN were shown to be involved in the regulation of stomatal development which influences resistance to drought and diseases. Interestingly, recent studies have also shown that SHN TFs are involved in mediating the beneficial effects of arbuscular mycorrhizal fungi (AMF) as well as disease resistance conferred by nanoparticles. To fulfill their roles, SHN TFs are controlled upstream by other TFs and they control, in their turn, different downstream genes. In this review, we highlight the role of SHN TFs in different abiotic and biotic stresses through their involvement in cuticle biosynthesis, stomatal development and molecular regulation of biochemical and physiological traits. In addition, we discuss the regulation of SHN TFs by plant hormones and their influence on hormone biosynthesis and signaling pathways. Knowledge of this complex regulation can be put into contribution to increase multiple abiotic stress tolerances through transgenesis, gene editing and classical breeding.

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.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Unlocking the Hidden Genetic Diversity of Varicosaviruses, the Neglected Plant Rhabdoviruses

The genus Varicosavirus is one of six genera of plant-infecting rhabdoviruses. Varicosaviruses have non-enveloped, flexuous, rod-shaped virions and a negative-sense, single-stranded RNA genome. A distinguishing feature of varicosaviruses, which is shared with dichorhaviruses, is a bi-segmented genome. Before 2017, a sole varicosavirus was known and characterized, and then two more varicosaviruses were identified through high-throughput sequencing in 2017 and 2018. More recently, the number of known varicosaviruses has substantially increased in concert with the extensive use of high-throughput sequencing platforms and data mining approaches. The novel varicosaviruses have revealed not only sequence diversity, but also plasticity in terms of genome architecture, including a virus with a tentatively unsegmented genome. Here, we report the discovery of 45 novel varicosavirus genomes which were identified in publicly available metatranscriptomic data. The identification, assembly, and curation of the raw Sequence Read Archive reads has resulted in 39 viral genome sequences with full-length coding regions and 6 with nearly complete coding regions. The highlights of the obtained sequences include eight varicosaviruses with unsegmented genomes, which are linked to a phylogenetic clade associated with gymnosperms. These findings have resulted in the most complete phylogeny of varicosaviruses to date and shed new light on the phylogenetic relationships and evolutionary landscape of this group of plant rhabdoviruses. Thus, the extensive use of sequence data mining for virus discovery has allowed us to unlock of the hidden genetic diversity of varicosaviruses, the largely neglected plant rhabdoviruses.

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

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

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.

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.

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

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

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.

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.

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.

Lycophyte transcriptomes reveal two whole-genome duplications in Lycopodiaceae: Insights into the polyploidization of Phlegmariurus

Lycophytes are an ancient clade of the non-flowering vascular plants with chromosome numbers that vary from tens to hundreds. They are an excellent study system for examining whole-genome duplications (WGDs), or polyploidization, in spore-dispersed vascular plants. However, a lack of genome sequence data limits the reliable detection of very ancient WGDs, small-scale duplications (SSDs), and recent WGDs. Here, we integrated phylogenomic analysis and the distribution of synonymous substitutions per synonymous sites (Ks) of the transcriptomes of 13 species of lycophytes to identify, locate, and date multiple WGDs in the lycophyte family Lycopodiaceae. Additionally, we examined the genus Phlegmariurus for signs of genetic discordance, which can provide valuable insight into the underlying causes of such conflict (e.g., hybridization, incomplete lineage sorting, or horizontal gene transfer).We found strong evidence that two WGD events occurred along the phylogenetic backbone of Lycopodiaceae, with one occurring in the common ancestor of extant Phlegmariurus (Lycopodiaceae) approximately 22-23 million years ago (Mya) and the other occurring in the common ancestor of Lycopodiaceae around 206-214 Mya. Interestingly, we found significant genetic discordance in the genus Phlegmariurus, indicating that the genus has a complex evolutionary history. This study provides molecular evidence for multiple WGDs in Lycopodiaceae and offers phylogenetic clues to the evolutionary history of Lycopodiaceae.

Analysis of Comparative Transcriptome and Positively Selected Genes Reveal Adaptive Evolution in Leaf-Less and Root-Less Whisk Ferns

While roots and leaves have evolved independently in lycophytes, ferns and seed plants, there is still confusion regarding the morphological evolution of ferns, especially in whisk ferns, which lack true leaves and roots and instead only exhibit leaf-like appendages and absorptive rhizoids. In this study, analyses of comparative transcriptomics on positively selected genes were performed to provide insights into the adaptive evolution of whisk fern morphologies. Significantly clustered gene families specific to whisk ferns were mainly enriched in Gene Ontology (GO) terms “binding proteins” and “transmembrane transporter activity”, and positive selection was detected in genes involved in transmembrane transporter activities and stress response (e.g., sodium/hydrogen exchanger and heat shock proteins), which could be related to the adaptive evolution of tolerance to epiphytic environments. The analysis of TF/TR gene family sizes indicated that some rapidly evolving gene families (e.g., the GRF and the MADS-MIKC families) related to the development of morphological organs were commonly reduced in whisk ferns and ophioglossoid ferns. Furthermore, the WUS homeobox-containing (WOX) gene family and the knotted1-like homeobox (KNOX) gene family, both associated with root and leaf development, were phylogenetically conserved in whisk ferns and ophioglossoid ferns. In general, our results suggested that adaptive evolution to epiphytic environments might have occurred in whisk ferns. We propose that the simplified and reduced leaf and root system in whisk ferns is the result of reduction from the common ancestor of whisk ferns and ophioglossoid ferns, rather than an independent origin.

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Evolution and expression of LEAFY genes in ferns and lycophytes

BACKGROUND

The LEAFY (LFY) transcription factors are present in algae and across land plants. The available expression and functional data of these genes in embryophytes suggest that LFY genes control a plethora of processes including the first zygotic cell division in bryophytes, shoot cell divisions of the gametophyte and sporophyte in ferns, cone differentiation in gymnosperms and floral meristem identity in flowering plants. However, their putative plesiomorphic role in plant reproductive transition in vascular plants remains untested.

RESULTS

We perform Maximum Likelihood (ML) phylogenetic analyses for the LFY gene lineage in embryophytes with expanded sampling in lycophytes and ferns. We recover the previously identified seed plant duplication that results in LEAFY and NEEDLY paralogs. In addition, we recover multiple species-specific duplications in ferns and lycophytes and large-scale duplications possibly correlated with the occurrence of whole genome duplication (WGD) events in Equisetales and Salviniales. To test putative roles in diverse ferns and lycophytes we perform LFY expression analyses in Adiantum raddianum, Equisetum giganteum and Selaginella moellendorffii. Our results show that LFY genes are active in vegetative and reproductive tissues, with higher expression in early fertile developmental stages and during sporangia differentiation.

CONCLUSIONS

Our data point to previously unrecognized roles of LFY genes in sporangia differentiation in lycophytes and ferns and suggests that functions linked to reproductive structure development are not exclusive to seed plant LFY homologs.

Inferring the Potential Geographic Distribution and Reasons for the Endangered Status of the Tree Fern, Sphaeropteris lepifera, in Lingnan, China Using a Small Sample Size

In this study, we investigated suitable habitats for the endangered tree fern, Sphaeropteris lepifera (J. Sm. ex Hook.) R.M. Tryon, based on fieldwork, ecological niche modeling, and regression approaches. We combined these data with the characterization of spore germination and gametophytic development in the laboratory to assess the reasons why S. lepifera is endangered and to propose a conservation strategy that focuses on suitable sites for reintroduction and accounts for the ecology and biphasic life cycle of the species. Our methods represent an integration of process- and correlation-based approaches to understanding the distributional patterns of this species, and this combined approach, while uncommonly applied, is a more robust strategy than either approach used in isolation. Our ecological niche models indicated that cold temperature extremes, temperature stability over long- and short-terms, and the seasonality of precipitation were among the most important abiotic environmental factors affecting the distribution of S. lepifera among the variables that we measured. Moreover, distribution of this fern species is also strongly influenced by the timing of development of male and female gametes. Additionally, we observed that slope aspect, specifically south-facing slopes, facilitates more incoming sunlight for mature trees, and simultaneously, provides greater, much-needed shade for fiddleheads on account of the canopy being denser. We believe that our study can provide important guidance on the restoration of S. lepifera in the wild. Specifically, potential restoration areas can be screened for the specific environmental factors that we infer to have a critical impact on the survival of the species.

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.

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.

A drought-driven model for the evolution of obligate apomixis in ferns: evidence from pellaeids (Pteridaceae)

PREMISE

Xeric environments impose major constraints on the fern life cycle, yet many lineages overcome these limitations by evolving apomixis. Here, we synthesize studies of apomixis in ferns and present an evidence-based model for the evolution and establishment of this reproductive strategy, focusing on genetic and environmental factors associated with its two defining traits: the production of "unreduced" spores (n = 2n) and the initiation of sporophytes from gametophyte tissue (i.e., diplospory and apogamy, respectively).

METHODS

We evaluated existing literature in light of the hypothesis that abiotic characteristics of desert environments (e.g., extreme diurnal temperature fluctuations, high light intensity, and water limitation) drive the evolution of obligate apomixis. Pellaeid ferns (Cheilanthoideae: Pteridaceae) were examined in detail, as an illustrative example. We reconstructed a plastid (rbcL, trnG-trnR, atpA) phylogeny for the clade and mapped reproductive mode (sexual versus apomictic) and ploidy across the resulting tree.

RESULTS

Our six-stage model for the evolution of obligate apomixis in ferns emphasizes the role played by drought and associated abiotic conditions in the establishment of this reproductive approach. Furthermore, our updated phylogeny of pellaeid ferns reveals repeated origins of obligate apomixis and shows an increase in the frequency of apomixis, and rarity of sexual reproduction, among taxa inhabiting increasingly dry North American deserts.

CONCLUSIONS

Our findings reinforce aspects of other evolutionary, physiological, developmental, and omics-based studies, indicating a strong association between abiotic factors and the establishment of obligate apomixis in ferns. Water limitation, in particular, appears critical to establishment of this reproductive mode.