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

A 160 Gbp fork fern genome shatters size record for eukaryotes

Vascular plants are exceptional among eukaryotes due to their outstanding genome size diversity which ranges ∼2,400-fold, including the largest genome so far recorded in the angiosperm Paris japonica (148.89 Gbp/1C). Despite available data showing that giant genomes are restricted across the Tree of Life, the biological limits to genome size expansion remain to be established. Here, we report the discovery of an even larger eukaryotic genome in Tmesipteris oblanceolata, a New Caledonian fork fern. At 160.45 Gbp/1C, this record-breaking genome challenges current understanding and opens new avenues to explore the evolutionary dynamics of genomic gigantism.

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.

Modified CTAB protocols for high-molecular-weight DNA extractions from ferns

Premise

Efficient protocols for extracting high-molecular-weight (HMW) DNA from ferns facilitate the long-read sequencing of their large and complex genomes. Here, we perform two cetyltrimethylammonium bromide (CTAB)-based protocols to extract HMW DNA and evaluate their applicability in diverse fern taxa for the first time.

Methods and Results

We describe two modified CTAB protocols, with key adjustments to minimize mechanical disruption during lysis to prevent DNA shearing. One of these protocols uses a small amount of fresh tissue but yields a considerable quantity of HMW DNA with high efficiency. The other accommodates a large amount of input tissue, adopts an initial step of nuclei isolation, and thus ensures a high yield in a short period of time. Both methods were proven to be robust and effective in obtaining HMW DNA from diverse fern lineages, including 33 species in 19 families. The DNA extractions mostly had high DNA integrity, with mean sizes larger than 50 kbp, as well as high purity (A₂₆₀/A₂₃₀ and A₂₆₀/A₂₈₀ > 1.8).

Conclusions

This study provides HMW DNA extraction protocols for ferns in the hope of facilitating further attempts to sequence their genomes, which will bridge our genomic understanding of land plant diversity.

Polyploidy: its consequences and enabling role in plant diversification and evolution

BACKGROUND

Most, if not all, green plant (Virdiplantae) species including angiosperms and ferns are polyploids themselves or have ancient polyploid or whole genome duplication signatures in their genomes. Polyploids are not only restricted to our major crop species such as wheat, maize, potato and the brassicas, but also occur frequently in wild species and natural habitats. Polyploidy has thus been viewed as a major driver in evolution, and its influence on genome and chromosome evolution has been at the centre of many investigations. Mechanistic models of the newly structured genomes are being developed that incorporate aspects of sequence evolution or turnover (low-copy genes and regulatory sequences, as well as repetitive DNAs), modification of gene functions, the re-establishment of control of genes with multiple copies, and often meiotic chromosome pairing, recombination and restoration of fertility.

SCOPE

World-wide interest in how green plants have evolved under different conditions - whether in small, isolated populations, or globally - suggests that gaining further insight into the contribution of polyploidy to plant speciation and adaptation to environmental changes is greatly needed. Forward-looking research and modelling, based on cytogenetics, expression studies, and genomics or genome sequencing analyses, discussed in this Special Issue of the Annals of Botany, consider how new polyploids behave and the pathways available for genome evolution. They address fundamental questions about the advantages and disadvantages of polyploidy, the consequences for evolution and speciation, and applied questions regarding the spread of polyploids in the environment and challenges in breeding and exploitation of wild relatives through introgression or resynthesis of polyploids.

CONCLUSION

Chromosome number, genome size, repetitive DNA sequences, genes and regulatory sequences and their expression evolve following polyploidy - generating diversity and possible novel traits and enabling species diversification. There is the potential for ever more polyploids in natural, managed and disturbed environments under changing climates and new stresses.

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.

Fern Spores-"Ready-to-Use" Standards for Plant Genome Size Estimation Using a Flow Cytometric Approach

Spores and pollen of plants were used as flow cytometric materials to efficiently infer genome sizes. Given this advantage, they hold great potential for various flow cytometric applications, particularly as plant genome size standards. To develop such novel standards, we investigated conditions of pretreatment (bead vortex), buffer, and reliable genome sizes of three fern spore collections-Cibotium taiwanense "Kuo4395", Sphaeropteris lepifera "Tang0001", and Alsophila metteniana "Lee s.n.". Additionally, up to 30 year-old spore collections were obtained from herbarium specimens and from samples stored at 4 °C; their spore nuclei were extracted, and the quality and quantity of these nucleus extractions through storage ages were examined. Nuclear extractions with a longer bead vortex duration or lower spore/bead ratio generally resulted in a higher recovered quantity but a lower quality or purity. For each spore standard, the protocol optimization was determined by their performance in bead vortex conditions, and a 1C genome size was further inferred by linear regression (C. taiwanense "Kuo4395" = 5.058 pg; S. lepifera "Tang0001" = 7.117 pg; and A. metteniana "Lee s.n." = 19.379 pg). Spore nucleus quality and quantity are significantly negatively correlated with storage ages. Nuclear extractions of 10-year-old refrigerated spores remained qualified as a genome size standard; however, none of the herbarium spore collections fit such criteria. Our study is the first to develop and apply dried and refrigerated spores for genome size standards. These standards are ready to use, easy to manipulate, and feature long-term storage in comparison with traditionally used standards of fresh leaves.

Kinetochore size scales with chromosome size in bimodal karyotypes of Agavoideae

BACKGROUND AND AIMS

In eukaryotes, the total kinetochore size (defined as a chromosomal region containing CENH3-positive nucleosomes) per nucleus strongly correlates with genome size, a relationship that has been hypothesized to stem from general intracellular scaling principles. However, if larger chromosomes within a karyotype required larger kinetochores to move properly, it could also be derived from the mechanics of cell division.

METHODS

We selected seven species of the plant subfamily Agavoideae whose karyotypes are characterized by the presence of small and very large chromosomes. We visualized the kinetochore regions and chromosomes by immunolabelling with an anti-CENH3 antibody and DAPI (6'-diamidino-2-phenylindole) staining. We then employed 2D widefield and 3D super-resolution microscopy to measure chromosome and kinetochore areas and volumes, respectively. To assess the scaling relationship of kinetochore size to chromosome size inside a karyotype, we log-transformed the data and analysed them with linear mixed models which allowed us to control for the inherent hierarchical structure of the dataset (metaphases within slides and species).

KEY RESULTS

We found a positive intra-karyotype relationship between kinetochore and chromosome size. The slope of the regression line of the observed relationship (0.277 for areas, 0.247 for volumes) was very close to the theoretical slope of 0.25 for chromosome width based on the expected physics of chromosome passage through the cytoplasm during cell division. We obtained similar results by reanalysing available data from human and maize.

CONCLUSIONS

Our findings suggest that the total kinetochore size to genome size scaling observed across eukaryotes may also originate from the mechanics of cell division. Moreover, the potential causal link between kinetochore and chromosome size indicates that evolutionary mechanisms capable of leading kinetochore size changes to fixation, such as centromere drive, could promote the size evolution of entire chromosomes and genomes.

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.

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.

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.

A dormant resource for genome size estimation in ferns: C-value inference of the Ophioglossaceae using herbarium specimen spores

PREMISE

The great variation of genome size (C-value) across land plants is linked to various adaptative features. Flow cytometry (FCM), the standard approach to estimating C-values, relies mostly on fresh materials, performing poorly when used with herbarium materials. No fern C-value reports have been derived from herbarium specimens; however, the herbarium spores of some ferns remain highly viable for decades and are thus promising for further investigation. To explore this possibility, we evaluated herbarium spore collections of Ophioglossaceae ferns using FCM.

METHODS

Flow cytometry was conducted on 24 spore samples, representing eight of the 12 genera of the Ophioglossaceae, using specimens ranging in age from 2.6 to 111 years obtained from five herbaria.

RESULTS

Regardless of the genus or the source herbarium, high-quality C-value data were generated from 17 samples, with the oldest being 26 years old. Estimates of the C-values from sporophytic tissues of known ploidy did not reveal any evidence of apomixis for the species surveyed here. We also detected a pronounced genome downsizing in Sceptridium polyploids.

DISCUSSION

The recent success of FCM for C-value estimation using spores provides a much more convenient method of utilizing "dry" refrigerated materials. We demonstrate here that herbarium spores of some ferns are also promising for this use, even for older specimens.