Unfurling Fern Biology in the Genomics Age

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.

Dynamic evolution of MADS-box genes in extant ferns via large-scale phylogenomic analysis

Introduction

Several studies of MADS-box transcription factors in flowering plants have been conducted, and these studies have indicated that they have conserved functions in floral organ development; MIKC-type MADS-box genes has been proved to be expanded in ferns, however, few systematic studies of these transcription factors have been conducted in non-seed plants. Although ferns and seed plants are sister groups, they exhibit substantial morphological differences.

Methods

Here, we clarified the evolution of MADS-box genes across 71 extant fern species using available transcriptome, genome, and gene expression data.

Results

We obtained a total of 2,512 MADS-box sequences, ranging from 9 to 89 per species. The most recent common ancestor (MRCA) of ferns contained approximately three type I genes and at least 5-6 type II MADS-box genes. The domains, motifs, expression of type I and type II proteins, and the structure of the both type genes were conserved in ferns as to other land plants. Within type II genes, MIKC*-type proteins are involved in gametophyte development in ferns; MIKCC-type proteins have broader expression patterns in ferns than in seed plants, and these protein sequences are likely conserved in extant seed plants and ferns because of their diverse roles in diploid sporophyte development. More than 90% of MADS-box genes are type II genes, and MIKCC genes, especially CRM1 and CRM6-like genes, have undergone a large expansion in leptosporangiate ferns; the diverse expression patterns of these genes might be related to the fuctional diversification and increased complexity of the plant body plan. Tandem duplication of CRM1 and CRM6-like genes has contributed to the expansion of MIKCC genes.

Conclusion or Discussion

This study provides new insights into the diversity, evolution, and functions of MADS-box genes in extant ferns.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Proteome and Interactome Linked to Metabolism, Genetic Information Processing, and Abiotic Stress in Gametophytes of Two Woodferns

Ferns and lycophytes have received scant molecular attention in comparison to angiosperms. The advent of high-throughput technologies allowed an advance towards a greater knowledge of their elusive genomes. In this work, proteomic analyses of heart-shaped gametophytes of two ferns were performed: the apomictic Dryopteris affinis ssp. affinis and its sexual relative Dryopteris oreades. In total, a set of 218 proteins shared by these two gametophytes were analyzed using the STRING database, and their proteome associated with metabolism, genetic information processing, and responses to abiotic stress is discussed. Specifically, we report proteins involved in the metabolism of carbohydrates, lipids, and nucleotides, the biosynthesis of amino acids and secondary compounds, energy, oxide-reduction, transcription, translation, protein folding, sorting and degradation, and responses to abiotic stresses. The interactome of this set of proteins represents a total network composed of 218 nodes and 1792 interactions, obtained mostly from databases and text mining. The interactions among the identified proteins of the ferns D. affinis and D. oreades, together with the description of their biological functions, might contribute to a better understanding of the function and development of ferns as well as fill knowledge gaps in plant evolution.

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.

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.

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.

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.

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.

Sexual and Apogamous Species of Woodferns Show Different Protein and Phytohormone Profiles

The gametophyte of ferns reproduces either by sexual or asexual means. In the latter, apogamy represents a peculiar case of apomixis, in which an embryo is formed from somatic cells. A proteomic and physiological approach was applied to the apogamous fern Dryopteris affinis ssp. affinis and its sexual relative D. oreades. The proteomic analysis compared apogamous vs. female gametophytes, whereas the phytohormone study included, in addition to females, three apogamous stages (filamentous, spatulate, and cordate). The proteomic profiles revealed a total of 879 proteins and, after annotation, different regulation was found in 206 proteins of D. affinis and 166 of its sexual counterpart. The proteins upregulated in D. affinis are mostly associated to protein metabolism (including folding, transport, and proteolysis), ribosome biogenesis, gene expression and translation, while in the sexual counterpart, they account largely for starch and sucrose metabolism, generation of energy and photosynthesis. Likewise, ultra-performance liquid chromatography-tandem spectrometry (UHPLC-MS/MS) was used to assess the levels of indol-3-acetic acid (IAA); the cytokinins: 6-benzylaminopurine (BA), trans-Zeatine (Z), trans-Zeatin riboside (ZR), dyhidrozeatine (DHZ), dyhidrozeatin riboside (DHZR), isopentenyl adenine (iP), isopentenyl adenosine (iPR), abscisic acid (ABA), the gibberellins GA₃ and GA₄, salicylic acid (SA), and the brassinosteroids: brassinolide (BL) and castasterone (CS). IAA, the cytokinins Z, ZR, iPR, the gibberellin GA₄, the brassinosteoids castasterone, and ABA accumulated more in the sexual gametophyte than in the apogamous one. When comparing the three apogamous stages, BA and SA peaked in filamentous, GA₃ and BL in spatulate and DHRZ in cordate gametophytes. The results point to the existence of large metabolic differences between apogamous and sexual gametophytes, and invite to consider the fern gametophyte as a good experimental system to deepen our understanding of plant reproduction.

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

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

Animal chromosome counts reveal a similar range of chromosome numbers but with less polyploidy in animals compared to flowering plants

Understanding the mechanisms that underlie chromosome evolution could provide insights into the processes underpinning the origin, persistence and evolutionary tempo of lineages. Here, we present the first database of chromosome counts for animals (the Animal Chromosome Count database, ACC) summarizing chromosome numbers for ~15,000 species. We found remarkable a similarity in the distribution of chromosome counts between animals and flowering plants. Nevertheless, the similarity in the distribution of chromosome numbers between animals and plants is likely to be explained by different drivers. For instance, we found that while animals and flowering plants exhibit similar frequencies of speciation-related changes in chromosome number, plant speciation is more often related to changes in ploidy. By leveraging the largest data set of chromosome counts for animals, we describe a previously undocumented pattern across the Tree of Life-animals and flowering plants show remarkably similar distributions of haploid chromosome numbers.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

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.

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

Ferns, the second largest group of vascular plants, originated ~400 million years ago (Mya). They became dominant in the ancient Earth landscape before the angiosperms and are still important in current ecosystems. Many ferns have exceptionally high chromosome numbers, possibly resulting from whole-genome duplications (WGDs). However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates (Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%-97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strong evidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous (~145-66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.

Polyploidy in gymnosperms - Insights into the genomic and evolutionary consequences of polyploidy in Ephedra

While polyploidization is recognized as a major evolutionary driver for ferns and angiosperms, little is known about its impact in gymnosperms, where polyploidy is much less frequent. We explore Ephedra to evaluate (i) the extent of genome size diversity in the genus and the influence polyploidy has had on the evolution of nuclear DNA contents, and (ii) identify where shifts in genome size and polyploidy have occurred both temporally and spatially. A phylogenetic framework of all Ephedra species together with genome sizes and karyotypes for 87% and 67% of them respectively, were used to explore ploidy evolution and its global distribution patterns. Polyploidy was shown to be extremely common, with 41 species (83%) being polyploid (up to 8×) or having polyploid cytotypes - the highest frequency and level reported for any gymnosperm. Genome size was also diverse, with values ranging ~5-fold (8.09-38.34 pg/1C) - the largest range for any gymnosperm family - and increasing in proportion to ploidy level (i.e. no genome downsizing). Our findings provide novel data which support the view that gymnosperms have a more conserved mode of genomic evolution compared with angiosperms.

Differential gene expression profiling of one- and two-dimensional apogamous gametophytes of the fern Dryopteris affinis ssp. affinis

Apomixis was originally defined as the replacement of sexual reproduction by an asexual process that does not involve fertilization but, in angiosperms, it is often used in the more restricted sense of asexual reproduction through seeds. In ferns, apomixis combines the production of unreduced spores (diplospory) and the formation of sporophytes from somatic cells of the prothallium (apogamy). The genes that control the onset of apogamy in ferns are largely unknown. In this study, we describe the gametophyte transcriptome of the apogamous fern Dryopteris affinis ssp. affinis using an RNA-Seq approach to compare the gene expression profiles of one- and two-dimensional gametophytes, the latter containing apogamic centers. After collapsing highly similar de novo transcripts, we obtained 166,191 unigenes, of which 30% could be annotated using public databases. Multiple quality metrics indicate a good quality of the de novo transcriptome with a low level of fragmentation. Our data show a total of 10,679 genes (6% of all genes) to be differentially expressed between gametophytes of filamentous (one-dimensional) and prothallial (two-dimensional) architecture. 6,110 genes were up-regulated in two-dimensional relative to one-dimensional gametophytes, some of which are implicated in the regulation of meristem growth, auxin signaling, reproduction, and sucrose metabolism. 4,570 genes were down-regulated in two-dimensional versus one-dimensional gametophytes, which are enriched in stimulus and defense genes, as well as genes involved in epigenetic gene regulation and ubiquitin degradation. Our results provide insights into free-living gametophyte development, focusing on the filamentous-to-prothallus growth transition, and provide a useful resource for further investigations of asexual reproduction.

The evolutionary trajectory of root stem cells

Root stem cells are crucial for the establishment of roots and are therefore a major evolutionary innovation that enabled land plants to spread on land. Despite their importance, not too much is known about the origin and the molecular players installing and maintaining them. Although still fragmentary, the recent availability of new data for early land plants can be used to identify and analyze the conservation of key regulators of root meristems. In this review, we evaluate the possible conservation of important root stem cell regulators to suggest pathways that might have been important at the origin of roots.

One thousand plant transcriptomes and the phylogenomics of green plants

Green plants (Viridiplantae) include around 450,000-500,000 species1,2 of great diversity and have important roles in terrestrial and aquatic ecosystems. Here, as part of the One Thousand Plant Transcriptomes Initiative, we sequenced the vegetative transcriptomes of 1,124 species that span the diversity of plants in a broad sense (Archaeplastida), including green plants (Viridiplantae), glaucophytes (Glaucophyta) and red algae (Rhodophyta). Our analysis provides a robust phylogenomic framework for examining the evolution of green plants. Most inferred species relationships are well supported across multiple species tree and supermatrix analyses, but discordance among plastid and nuclear gene trees at a few important nodes highlights the complexity of plant genome evolution, including polyploidy, periods of rapid speciation, and extinction. Incomplete sorting of ancestral variation, polyploidization and massive expansions of gene families punctuate the evolutionary history of green plants. Notably, we find that large expansions of gene families preceded the origins of green plants, land plants and vascular plants, whereas whole-genome duplications are inferred to have occurred repeatedly throughout the evolution of flowering plants and ferns. The increasing availability of high-quality plant genome sequences and advances in functional genomics are enabling research on genome evolution across the green tree of life.

Dating Whole Genome Duplication in Ceratopteris thalictroides and Potential Adaptive Values of Retained Gene Duplicates

Whole-genome duplications (WGDs) are widespread in plants and frequently coincide with global climatic change events, such as the Cretaceous–Tertiary (KT) extinction event approximately 65 million years ago (mya). Ferns have larger genomes and higher chromosome numbers than seed plants, which likely resulted from multiple rounds of polyploidy. Here, we use diploid and triploid material from a model fern species, Ceratopteris thalictroides, for the detection of WGDs. High-quality RNA-seq data was used to infer the number of synonymous substitutions per synonymous site (Ks) between paralogs; Ks age distribution and absolute dating approach were used to determine the age of WGD events. Evidence of an ancient WGD event with a Ks peak value of approximately 1.2 was obtained for both samples; however, the Ks frequency distributions varied significantly. Importantly, we dated the WGD event at 51–53 mya, which coincides with the Paleocene-Eocene Thermal Maximum (PETM), when the Earth became warmer and wetter than any other period during the Cenozoic. Duplicate genes were preferentially retained for specific functions, such as environment response, further support that the duplicates may have promoted quick adaption to environmental changes and potentially resulted in evolutionary success, especially for pantropical species, such as C. thalictroides, which exhibits higher temperature tolerance.

Overlapping Patterns of Gene Expression Between Gametophyte and Sporophyte Phases in the Fern Polypodium amorphum (Polypodiales)

Ferns are unique among land plants in having sporophyte and gametophyte phases that are both free living and fully independent. Here, we examine patterns of sporophytic and gametophytic gene expression in the fern Polypodium amorphum, a member of the homosporous polypod lineage that comprises 80% of extant fern diversity, to assess how expression of a common genome is partitioned between two morphologically, ecologically, and nutritionally independent phases. Using RNA-sequencing, we generated transcriptome profiles for three replicates of paired samples of sporophyte leaf tissue and whole gametophytes to identify genes with significant differences in expression between the two phases. We found a nearly 90% overlap in the identity and expression levels of the genes expressed in both sporophytes and gametophytes, with less than 3% of genes uniquely expressed in either phase. We compare our results to those from similar studies to establish how phase-specific gene expression varies among major land plant lineages. Notably, despite having greater similarity in the identity of gene families shared between P. amorphum and angiosperms, P. amorphum has phase-specific gene expression profiles that are more like bryophytes and lycophytes than seed plants. Our findings suggest that shared patterns of phase-specific gene expression among seed-free plants likely reflect having relatively large, photosynthetic gametophytes (compared to the gametophytes of seed plants that are highly reduced). Phylogenetic analyses were used to further investigate the evolution of phase-specific expression for the phototropin, terpene synthase, and MADS-box gene families.

Proteogenomic Analysis Greatly Expands the Identification of Proteins Related to Reproduction in the Apogamous Fern Dryopteris affinis ssp. affinis

Performing proteomic studies on non-model organisms with little or no genomic information is still difficult. However, many specific processes and biochemical pathways occur only in species that are poorly characterized at the genomic level. For example, many plants can reproduce both sexually and asexually, the first one allowing the generation of new genotypes and the latter their fixation. Thus, both modes of reproduction are of great agronomic value. However, the molecular basis of asexual reproduction is not well understood in any plant. In ferns, it combines the production of unreduced spores (diplospory) and the formation of sporophytes from somatic cells (apogamy). To set the basis to study these processes, we performed transcriptomics by next-generation sequencing (NGS) and shotgun proteomics by tandem mass spectrometry in the apogamous fern D. affinis ssp. affinis. For protein identification we used the public viridiplantae database (VPDB) to identify orthologous proteins from other plant species and new transcriptomics data to generate a "species-specific transcriptome database" (SSTDB). In total 1,397 protein clusters with 5,865 unique peptide sequences were identified (13 decoy proteins out of 1,410, protFDR 0.93% on protein cluster level). We show that using the SSTDB for protein identification increases the number of identified peptides almost four times compared to using only the publically available VPDB. We identified homologs of proteins involved in reproduction of higher plants, including proteins with a potential role in apogamy. With the increasing availability of genomic data from non-model species, similar proteogenomics approaches will improve the sensitivity in protein identification for species only distantly related to models.

Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum

The mangrove fern genus Acrostichum grows in the extremely unstable marine intertidal zone under harsh conditions, such as high salt concentrations, tidal rhythms and long-term climate changes. To explore the phylogenetic relationships and molecular mechanisms underlying adaptations in this genus, we sequenced the transcriptomes of two species of Acrostichum, A. aureum and A. speciosum, as well as a species in the sister genus, Ceratopteris thalictroides. We obtained 47,517, 36,420 and 60,823 unigenes for the three ferns, of which 24.39-45.63% were annotated using public databases. The estimated divergence time revealed that Acrostichum adapted to the coastal region during the late Cretaceous, whereas the two mangrove ferns from the Indo West-Pacific (IWP) area diverged more recently. Two methods (the modified branch-site model and the K_h method) were used to identify several positively selected genes, which may contribute to differential adaptation of the two Acrostichum species to different light and salt conditions. Our study provides abundant transcriptome data and new insights into the evolution and adaptations of mangrove ferns in the inhospitable intertidal zone.

Transcriptome sequencing reveals genome-wide variation in molecular evolutionary rate among ferns

BACKGROUND

Transcriptomics in non-model plant systems has recently reached a point where the examination of nuclear genome-wide patterns in understudied groups is an achievable reality. This progress is especially notable in evolutionary studies of ferns, for which molecular resources to date have been derived primarily from the plastid genome. Here, we utilize transcriptome data in the first genome-wide comparative study of molecular evolutionary rate in ferns. We focus on the ecologically diverse family Pteridaceae, which comprises about 10 % of fern diversity and includes the enigmatic vittarioid ferns-an epiphytic, tropical lineage known for dramatically reduced morphologies and radically elongated phylogenetic branch lengths. Using expressed sequence data for 2091 loci, we perform pairwise comparisons of molecular evolutionary rate among 12 species spanning the three largest clades in the family and ask whether previously documented heterogeneity in plastid substitution rates is reflected in their nuclear genomes. We then inquire whether variation in evolutionary rate is being shaped by genes belonging to specific functional categories and test for differential patterns of selection.

RESULTS

We find significant, genome-wide differences in evolutionary rate for vittarioid ferns relative to all other lineages within the Pteridaceae, but we recover few significant correlations between faster/slower vittarioid loci and known functional gene categories. We demonstrate that the faster rates characteristic of the vittarioid ferns are likely not driven by positive selection, nor are they unique to any particular type of nucleotide substitution.

CONCLUSIONS

Our results reinforce recently reviewed mechanisms hypothesized to shape molecular evolutionary rates in vittarioid ferns and provide novel insight into substitution rate variation both within and among fern nuclear genomes.

The Small Nuclear Genomes of Selaginella Are Associated with a Low Rate of Genome Size Evolution

The haploid nuclear genome size (1C DNA) of vascular land plants varies over several orders of magnitude. Much of this observed diversity in genome size is due to the proliferation and deletion of transposable elements. To date, all vascular land plant lineages with extremely small nuclear genomes represent recently derived states, having ancestors with much larger genome sizes. The Selaginellaceae represent an ancient lineage with extremely small genomes. It is unclear how small nuclear genomes evolved in Selaginella We compared the rates of nuclear genome size evolution in Selaginella and major vascular plant clades in a comparative phylogenetic framework. For the analyses, we collected 29 new flow cytometry estimates of haploid genome size in Selaginella to augment publicly available data. Selaginella possess some of the smallest known haploid nuclear genome sizes, as well as the lowest rate of genome size evolution observed across all vascular land plants included in our analyses. Additionally, our analyses provide strong support for a history of haploid nuclear genome size stasis in Selaginella Our results indicate that Selaginella, similar to other early diverging lineages of vascular land plants, has relatively low rates of genome size evolution. Further, our analyses highlight that a rapid transition to a small genome size is only one route to an extremely small genome.

Genome evolution of ferns: evidence for relative stasis of genome size across the fern phylogeny

The genome evolution of ferns has been considered to be relatively static compared with angiosperms. In this study, we analyse genome size data and chromosome numbers in a phylogenetic framework to explore three hypotheses: the correlation of genome size and chromosome number, the origin of modern ferns from ancestors with high chromosome numbers, and the occurrence of several whole-genome duplications during the evolution of ferns. To achieve this, we generated new genome size data, increasing the percentage of fern species with genome sizes estimated to 2.8% of extant diversity, and ensuring a comprehensive phylogenetic coverage including at least three species from each fern order. Genome size was correlated with chromosome number across all ferns despite some substantial variation in both traits. We observed a trend towards conservation of the amount of DNA per chromosome, although Osmundaceae and Psilotaceae have substantially larger chromosomes. Reconstruction of the ancestral genome traits suggested that the earliest ferns were already characterized by possessing high chromosome numbers and that the earliest divergences in ferns were correlated with substantial karyological changes. Evidence for repeated whole-genome duplications was found across the phylogeny. Fern genomes tend to evolve slowly, albeit genome rearrangements occur in some clades.

Polyploidy and the proteome

Although major advances have been made during the past 20 years in our understanding of the genetic and genomic consequences of polyploidy, our knowledge of polyploidy and the proteome is in its infancy. One of our goals is to stimulate additional study, particularly broad-scale proteomic analyses of polyploids and their progenitors. Although it may be too early to generalize regarding the extent to which transcriptomic data are predictive of the proteome of polyploids, it is clear that the proteome does not always reflect the transcriptome. Despite limited data, important observations on the proteomes of polyploids are emerging. In some cases, proteomic profiles show qualitatively and/or quantitatively non-additive patterns, and proteomic novelty has been observed. Allopolyploids generally combine the parental contributions, but there is evidence of parental dominance of one contributing genome in some allopolyploids. Autopolyploids are typically qualitatively identical to but quantitatively different from their parents. There is also evidence of parental legacy at the proteomic level. Proteomes clearly provide insights into the consequences of genomic merger and doubling beyond what is obtained from genomic and/or transcriptomic data. Translating proteomic changes in polyploids to differences in morphology and physiology remains the holy grail of polyploidy--this daunting task of linking genotype to proteome to phenotype should emerge as a focus of polyploidy research in the next decade. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.

Mitosis and Interphase of the Highly Polyploid Palm Voanioalagerardii (2n = 606 ± 3)

The endemic, highly polyploid, monotypic Madagascan palm genus Voanioala (2n ≈ 606) was studied with regard to mitotic stages and interphase. Features of the cell cycle, morphology and sizes of metaphase chromosomes, fluorochrome banding patterns, and silver staining of NORs of such an extremely high polyploid organism are reported for the first time. On a whole, karyokinesis appears to be stable and efficient. A comparison with closely related palm taxa reveals that V. gerardii is 38-ploid, and comparison with the closely related genera Butia, Cocos (coconut) and Jubaea shows that Voanioala has lost ∼ 35% of its DNA amount subsequent to polyploidization and has suppressed between 74 and 88% of the original nucleolar organizers. About 10 active NORs are present in the nuclei. An auto- or allopolyploid origin of Voanioala is discussed with respect to currently available nuclear gene data. The biogeographic relations to Jubaeopsis, a closely related, monotypic, apparently likewise relict palm genus from eastern mainland South Africa are discussed. From a cytogenetic point of view, a common polyploid ancestor of both genera is most likely, but the available molecular phylogenetic data are not univocal.

An Exploration into Fern Genome Space

Ferns are one of the few remaining major clades of land plants for which a complete genome sequence is lacking. Knowledge of genome space in ferns will enable broad-scale comparative analyses of land plant genes and genomes, provide insights into genome evolution across green plants, and shed light on genetic and genomic features that characterize ferns, such as their high chromosome numbers and large genome sizes. As part of an initial exploration into fern genome space, we used a whole genome shotgun sequencing approach to obtain low-density coverage (∼0.4X to 2X) for six fern species from the Polypodiales (Ceratopteris, Pteridium, Polypodium, Cystopteris), Cyatheales (Plagiogyria), and Gleicheniales (Dipteris). We explore these data to characterize the proportion of the nuclear genome represented by repetitive sequences (including DNA transposons, retrotransposons, ribosomal DNA, and simple repeats) and protein-coding genes, and to extract chloroplast and mitochondrial genome sequences. Such initial sweeps of fern genomes can provide information useful for selecting a promising candidate fern species for whole genome sequencing. We also describe variation of genomic traits across our sample and highlight some differences and similarities in repeat structure between ferns and seed plants.

Cytological and Proteomic Analyses of Osmunda cinnamomea Germinating Spores Reveal Characteristics of Fern Spore Germination and Rhizoid Tip Growth

Fern spore is a good single-cell model for studying the sophisticated molecular networks in asymmetric cell division, differentiation, and polar growth. Osmunda cinnamomea L. var. asiatica is one of the oldest fern species with typical separate-growing trophophyll and sporophyll. The chlorophyllous spores generated from sporophyll can germinate without dormancy. In this study, the spore ultrastructure, antioxidant enzyme activities, as well as protein and gene expression patterns were analyzed in the course of spore germination at five typical stages (i.e. mature spores, rehydrated spores, double-celled spores, germinated spores, and spores with protonemal cells). Proteomic analysis revealed 113 differentially expressed proteins, which were mainly involved in photosynthesis, reserve mobilization, energy supplying, protein synthesis and turnover, reactive oxygen species scavenging, signaling, and cell structure modulation. The presence of multiple proteoforms of 25 differentially expressed proteins implies that post-translational modification may play important roles in spore germination. The dynamic patterns of proteins and their encoding genes exhibited specific characteristics in the processes of cell division and rhizoid tip growth, which include heterotrophic and autotrophic metabolisms, de novo protein synthesis and active protein turnover, reactive oxygen species and hormone (brassinosteroid and ethylene) signaling, and vesicle trafficking and cytoskeleton dynamic. In addition, the function skew of proteins in fern spores highlights the unique and common mechanisms when compared with evolutionarily divergent spermatophyte pollen. These findings provide an improved understanding of the typical single-celled asymmetric division and polar growth during fern spore germination.

Genome size evolution in Ontario ferns (Polypodiidae): evolutionary correlations with cell size, spore size, and habitat type and an absence of genome downsizing

Genome size is known to correlate with a number of traits in angiosperms, but less is known about the phenotypic correlates of genome size in ferns. We explored genome size variation in relation to a suite of morphological and ecological traits in ferns. Thirty-six fern taxa were collected from wild populations in Ontario, Canada. 2C DNA content was measured using flow cytometry. We tested for genome downsizing following polyploidy using a phylogenetic comparative analysis to explore the correlation between 1Cx DNA content and ploidy. There was no compelling evidence for the occurrence of widespread genome downsizing during the evolution of Ontario ferns. The relationship between genome size and 11 morphological and ecological traits was explored using a phylogenetic principal component regression analysis. Genome size was found to be significantly associated with cell size, spore size, spore type, and habitat type. These results are timely as past and recent studies have found conflicting support for the association between ploidy/genome size and spore size in fern polyploid complexes; this study represents the first comparative analysis of the trend across a broad taxonomic group of ferns.

Ferns: the missing link in shoot evolution and development

Shoot development in land plants is a remarkably complex process that gives rise to an extreme diversity of forms. Our current understanding of shoot developmental mechanisms comes almost entirely from studies of angiosperms (flowering plants), the most recently diverged plant lineage. Shoot development in angiosperms is based around a layered multicellular apical meristem that produces lateral organs and/or secondary meristems from populations of founder cells at its periphery. In contrast, non-seed plant shoots develop from either single apical initials or from a small population of morphologically distinct apical cells. Although developmental and molecular information is becoming available for non-flowering plants, such as the model moss Physcomitrella patens, making valid comparisons between highly divergent lineages is extremely challenging. As sister group to the seed plants, the monilophytes (ferns and relatives) represent an excellent phylogenetic midpoint of comparison for unlocking the evolution of shoot developmental mechanisms, and recent technical advances have finally made transgenic analysis possible in the emerging model fern Ceratopteris richardii. This review compares and contrasts our current understanding of shoot development in different land plant lineages with the aim of highlighting the potential role that the fern C. richardii could play in shedding light on the evolution of underlying genetic regulatory mechanisms.

Accommodating the load: The transposable element content of very large genomes

Very large genomes, that is, those above 20 Gb, are rare but widely distributed throughout the eukaryotes. They are found within the diatoms, dinoflagellates, metazoans and green plants, but so far have not been found in the excavates. There is a known positive correlation between genome size and the proportion of the genome composed of transposable elements (TEs). Very large genomes may therefore be expected to be almost entirely composed of TEs. Of the large genomes examined, in the angiosperms, gymnosperms and the dinoflagellates only a small portion of the genome was identified as TEs, most of these genomes were unidentified and may be novel or diverse TEs. In the salamanders and lungfish, 25 to 47% of the genome were identifiable retrotransposons, that is, TEs that copy themselves before insertion. However, the predominant class of TEs found in the lungfish was not the same as that found in the salamanders. The little data we have at the moment suggests therefore that the diversity and abundance of TEs is variable between taxa with large genomes, similar to patterns found in taxa with smaller genomes. Based on results from the human genome, we suggest that the ‘missing’ portion of the lungfish and salamander genomes are old, highly divergent, and therefore inactive copies of TEs. The data available indicate that, unlike plants with large genomes, neither the lungfish nor the salamanders show an increased risk of extinction. Based on a slow rate of DNA loss in salamanders it has been suggested that the large salamander genome is the result of run-away genome expansion involving genome size increases via TE proliferation associated with reduced recombination rate. We know of no studies on DNA loss or recombination rates in lungfish genomes, however a similar scenario could describe the process of genome expansion in the lungfish. A series of waves of TE transposition and sequence decay would describe the pattern of TE content seen in both the lungfish and the salamanders. The lungfish and salamanders, therefore, may accommodate their large load of TEs because these TEs have accumulated gradually over a long period of time and have been subject to inactivation and decay.

Between two fern genomes

Ferns are the only major lineage of vascular plants not represented by a sequenced nuclear genome. This lack of genome sequence information significantly impedes our ability to understand and reconstruct genome evolution not only in ferns, but across all land plants. Azolla and Ceratopteris are ideal and complementary candidates to be the first ferns to have their nuclear genomes sequenced. They differ dramatically in genome size, life history, and habit, and thus represent the immense diversity of extant ferns. Together, this pair of genomes will facilitate myriad large-scale comparative analyses across ferns and all land plants. Here we review the unique biological characteristics of ferns and describe a number of outstanding questions in plant biology that will benefit from the addition of ferns to the set of taxa with sequenced nuclear genomes. We explain why the fern clade is pivotal for understanding genome evolution across land plants, and we provide a rationale for how knowledge of fern genomes will enable progress in research beyond the ferns themselves.

The evolution of the actin binding NET superfamily

The Arabidopsis Networked (NET) superfamily are plant-specific actin binding proteins which specifically label different membrane compartments and identify specialized sites of interaction between actin and membranes unique to plants. There are 13 members of the superfamily in Arabidopsis, which group into four distinct clades or families. NET homologs are absent from the genomes of metazoa and fungi; furthermore, in plantae, NET sequences are also absent from the genome of mosses and more ancient extant plant clades. A single family of the NET proteins is found encoded in the club moss genome, an extant species of the earliest vascular plants. Gymnosperms have examples from families 4 and 3, with a hybrid form of NET1 and 2 which shows characteristics of both NET1 and NET2. In addition to NET3 and 4 families, the NET1 and pollen-expressed NET2 families are found only as independent sequences in Angiosperms. This is consistent with the divergence of reproductive actin. The four families are conserved across Monocots and Eudicots, with the numbers of members of each clade expanding at this point, due, in part, to regions of genome duplication. Since the emergence of the NET superfamily at the dawn of vascular plants, they have continued to develop and diversify in a manner which has mirrored the divergence and increasing complexity of land-plant species.

Reconstructing relative genome size of vascular plants through geological time

The strong positive relationship evident between cell and genome size in both animals and plants forms the basis of using the size of stomatal guard cells as a proxy to track changes in plant genome size through geological time. We report for the first time a taxonomic fine-scale investigation into changes in stomatal guard-cell length and use these data to infer changes in genome size through the evolutionary history of land plants. Our data suggest that many of the earliest land plants had exceptionally large genome sizes and that a predicted overall trend of increasing genome size within individual lineages through geological time is not supported. However, maximum genome size steadily increases from the Mississippian (c. 360 million yr ago (Ma)) to the present. We hypothesise that the functional relationship between stomatal size, genome size and atmospheric CO2 may contribute to the dichotomy reported between preferential extinction of neopolyploids and the prevalence of palaeopolyploidy observed in DNA sequence data of extant vascular plants.

Genome size expansion and the relationship between nuclear DNA content and spore size in the Asplenium monanthes fern complex (Aspleniaceae)

BACKGROUND

Homosporous ferns are distinctive amongst the land plant lineages for their high chromosome numbers and enigmatic genomes. Genome size measurements are an under exploited tool in homosporous ferns and show great potential to provide an overview of the mechanisms that define genome evolution in these ferns. The aim of this study is to investigate the evolution of genome size and the relationship between genome size and spore size within the apomictic Asplenium monanthes fern complex and related lineages.

RESULTS

Comparative analyses to test for a relationship between spore size and genome size show that they are not correlated. The data do however provide evidence for marked genome size variation between species in this group. These results indicate that Asplenium monanthes has undergone a two-fold expansion in genome size.

CONCLUSIONS

Our findings challenge the widely held assumption that spore size can be used to infer ploidy levels within apomictic fern complexes. We argue that the observed genome size variation is likely to have arisen via increases in both chromosome number due to polyploidy and chromosome size due to amplification of repetitive DNA (e.g. transposable elements, especially retrotransposons). However, to date the latter has not been considered to be an important process of genome evolution within homosporous ferns. We infer that genome evolution, at least in some homosporous fern lineages, is a more dynamic process than existing studies would suggest.

Transcriptome-mining for single-copy nuclear markers in ferns

BACKGROUND

Molecular phylogenetic investigations have revolutionized our understanding of the evolutionary history of ferns-the second-most species-rich major group of vascular plants, and the sister clade to seed plants. The general absence of genomic resources available for this important group of plants, however, has resulted in the strong dependence of these studies on plastid data; nuclear or mitochondrial data have been rarely used. In this study, we utilize transcriptome data to design primers for nuclear markers for use in studies of fern evolutionary biology, and demonstrate the utility of these markers across the largest order of ferns, the Polypodiales.

PRINCIPAL FINDINGS

We present 20 novel single-copy nuclear regions, across 10 distinct protein-coding genes: ApPEFP_C, cryptochrome 2, cryptochrome 4, DET1, gapCpSh, IBR3, pgiC, SQD1, TPLATE, and transducin. These loci, individually and in combination, show strong resolving power across the Polypodiales phylogeny, and are readily amplified and sequenced from our genomic DNA test set (from 15 diploid Polypodiales species). For each region, we also present transcriptome alignments of the focal locus and related paralogs-curated broadly across ferns-that will allow researchers to develop their own primer sets for fern taxa outside of the Polypodiales. Analyses of sequence data generated from our genomic DNA test set reveal strong effects of partitioning schemes on support levels and, to a much lesser extent, on topology. A model partitioned by codon position is strongly favored, and analyses of the combined data yield a Polypodiales phylogeny that is well-supported and consistent with earlier studies of this group.

CONCLUSIONS

The 20 single-copy regions presented here more than triple the single-copy nuclear regions available for use in ferns. They provide a much-needed opportunity to assess plastid-derived hypotheses of relationships within the ferns, and increase our capacity to explore aspects of fern evolution previously unavailable to scientific investigation.

Ceratopteris richardii (C-fern): a model for investigating adaptive modification of vascular plant cell walls

Plant cell walls are essential for most aspects of plant growth, development, and survival, including cell division, expansive cell growth, cell-cell communication, biomechanical properties, and stress responses. Therefore, characterizing cell wall diversity contributes to our overall understanding of plant evolution and development. Recent biochemical analyses, concomitantly with whole genome sequencing of plants located at pivotal points in plant phylogeny, have helped distinguish between homologous characters and those which might be more derived. Most plant lineages now have at least one fully sequenced representative and although genome sequences for fern species are in progress they are not yet available for this group. Ferns offer key advantages for the study of developmental processes leading to vascularisation and complex organs as well as the specific differences between diploid sporophyte tissues and haploid gametophyte tissues and the interplay between them. Ceratopteris richardii has been well investigated building a body of knowledge which combined with the genomic and biochemical information available for other plants will progress our understanding of wall diversity and its impact on evolution and development.

RNA-seq analysis identifies potential modulators of gravity response in spores of Ceratopteris (Parkeriaceae): evidence for modulation by calcium pumps and apyrase activity

PREMISE OF THE STUDY

Gravity regulates the magnitude and direction of a trans-cell calcium current in germinating spores of Ceratopteris richardii. Blocking this current with nifedipine blocks the spore's downward polarity alignment, a polarization that is fixed by gravity ∼10 h after light induces the spores to germinate. RNA-seq analysis at 10 h was used to identify genes potentially important for the gravity response. The data set will be valuable for other developmental and phylogenetic studies.

METHODS

De novo Newbler assembly of 958 527 reads from Roche 454 sequencing was executed. The sequences were identified and analyzed using in silico methods. The roles of endomembrane Ca(2+)-ATPase pumps and apyrases in the gravity response were further tested using pharmacological agents.

KEY RESULTS

Transcripts related to calcium signaling and ethylene biosynthesis were identified as notable constituents of the transcriptome. Inhibiting the activity of endomembrane Ca(2+)-ATPase pumps with 2,5-di-(t-butyl)-1,4-hydroquinone diminished the trans-cell current, but increased the orientation of the polar axis to gravity. The effects of applied nucleotides and purinoceptor antagonists gave novel evidence implicating extracellular nucleotides as regulators of the gravity response in these fern spores.

CONCLUSIONS

In addition to revealing general features of the transcriptome of germinating spores, the results highlight a number of calcium-responsive and light-receptive transcripts. Pharmacologic assays indicate endomembrane Ca(2+)-ATPases and extracellular nucleotides may play regulatory roles in the gravity response of Ceratopteris spores.

Ecological and genetic factors linked to contrasting genome dynamics in seed plants

The large-scale replacement of gymnosperms by angiosperms in many ecological niches over time and the huge disparity in species numbers have led scientists to explore factors (e.g. polyploidy, developmental systems, floral evolution) that may have contributed to the astonishing rise of angiosperm diversity. Here, we explore genomic and ecological factors influencing seed plant genomes. This is timely given the recent surge in genomic data. We compare and contrast the genomic structure and evolution of angiosperms and gymnosperms and find that angiosperm genomes are more dynamic and diverse, particularly amongst the herbaceous species. Gymnosperms typically have reduced frequencies of a number of processes (e.g. polyploidy) that have shaped the genomes of other vascular plants and have alternative mechanisms to suppress genome dynamism (e.g. epigenetics and activity of transposable elements). Furthermore, the presence of several characters in angiosperms (e.g. herbaceous habit, short minimum generation time) has enabled them to exploit new niches and to be viable with small population sizes, where the power of genetic drift can outweigh that of selection. Together these processes have led to increased rates of genetic divergence and faster fixation times of variation in many angiosperms compared with gymnosperms.

Fertility and precocity of Osmunda x intermedia offspring in culture

The feasibility of later-generation hybrid production in ferns has not been previously studied, although it is a significant factor in relation to reproductive isolation. Osmunda x intermedia, a hybrid between O. japonica and O. lancea, is semifertile and has moderate spore germination rates. Under the artificial conditions of this study, F2 and F3 offspring were formed. Some of the F2 offspring showed precocity, and some of the F3 offspring also showed precocity. This fertility suggests that introgressive hybridization might be ongoing in nature. This also indicates a currently unknown genetic control over the timing of fertile frond production in Osmunda.

De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum

Background

Because of their phylogenetic position and unique characteristics of their biology and life cycle, ferns represent an important lineage for studying the evolution of land plants. Large and complex genomes in ferns combined with the absence of economically important species have been a barrier to the development of genomic resources. However, high throughput sequencing technologies are now being widely applied to non-model species. We leveraged the Roche 454 GS-FLX Titanium pyrosequencing platform in sequencing the gametophyte transcriptome of bracken fern (Pteridium aquilinum) to develop genomic resources for evolutionary studies.

Results

681,722 quality and adapter trimmed reads totaling 254 Mbp were assembled de novo into 56,256 unique sequences (i.e. unigenes) with a mean length of 547.2 bp and a total assembly size of 30.8 Mbp with an average read-depth coverage of 7.0×. We estimate that 87% of the complete transcriptome has been sequenced and that all transcripts have been tagged. 61.8% of the unigenes had blastx hits in the NCBI nr protein database, representing 22,596 unique best hits. The longest open reading frame in 52.2% of the unigenes had positive domain matches in InterProScan searches. We assigned 46.2% of the unigenes with a GO functional annotation and 16.0% with an enzyme code annotation. Enzyme codes were used to retrieve and color KEGG pathway maps. A comparative genomics approach revealed a substantial proportion of genes expressed in bracken gametophytes to be shared across the genomes of Arabidopsis, Selaginella and Physcomitrella, and identified a substantial number of potentially novel fern genes. By comparing the list of Arabidopsis genes identified by blast with a list of gametophyte-specific Arabidopsis genes taken from the literature, we identified a set of potentially conserved gametophyte specific genes. We screened unigenes for repetitive sequences to identify 548 potentially-amplifiable simple sequence repeat loci and 689 expressed transposable elements.

Conclusions

This study is the first comprehensive transcriptome analysis for a fern and represents an important scientific resource for comparative evolutionary and functional genomics studies in land plants. We demonstrate the utility of high-throughput sequencing of a normalized cDNA library for de novo transcriptome characterization and gene discovery in a non-model plant.