Dynamic genome evolution in a model fern.

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

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

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

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.

Cryptic diversity in the model fern genus Ceratopteris (Pteridaceae).

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

How Many Tree Species of Birch Are in Alaska? Implications for Wetland Designations.

Wetland areas are critical habitats, especially in northern regions of North America. Wetland classifications are based on several factors, including the presence of certain plant species and assemblages of species, of which trees play a significant role. Here we examined wetland species of birch (Betula) in North America, with a focus on Alaska, and the use of birche tree species in wetland delineation. We sampled over 200 trees from sites, including Alaska, Alberta, Minnesota, and New Hampshire. We used genetic data from over 3000 loci detected by restriction site associated DNA analysis. We used an indirect estimate of ploidy based on allelic ratios and we also examined population genetic structure. We find that inferred ploidy is strongly associated with genetic groupings. We find two main distinct groups; one found throughout most of Alaska, extending into Alberta. This group is probably attributable to Betula kenaica, Betula neoalaskana, or both. This group has a diploid genetic pattern although this could easily be a function of allopolyploidy. The second major genetic group appears to extend from Eastern North America into parts of southeastern Alaska. This group represents Betula papyrifera, and is not diploid based on allelic ratios. Published chromosome counts indicate pentaploidy. Because B. papyrifera is the only one of the above species that is distinctly associated with wetland habitats, our findings indicate that tree species of birch found in most parts of Alaska are not reliable indicators of wetland habitats. These results help to support stronger wetland ratings assigned to the tree species of birch for delineation purposes.

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

Ferns are notorious for possessing large genomes and numerous chromosomes. Despite decades of speculation, the processes underlying the expansive genomes of ferns are unclear, largely due to the absence of a sequenced homosporous fern genome. The lack of this crucial resource has not only hindered investigations of evolutionary processes responsible for the unusual genome characteristics of homosporous ferns, but also impeded synthesis of genome evolution across land plants. Here, we used the model fern species Ceratopteris richardii to address the processes (e.g., polyploidy, spread of repeat elements) by which the large genomes and high chromosome numbers typical of homosporous ferns may have evolved and have been maintained. We directly compared repeat compositions in species spanning the green plant tree of life and a diversity of genome sizes, as well as both short- and long-read-based assemblies of Ceratopteris. We found evidence consistent with a single ancient polyploidy event in the evolutionary history of Ceratopteris based on both genomic and cytogenetic data, and on repeat proportions similar to those found in large flowering plant genomes. This study provides a major stepping-stone in the understanding of land plant evolutionary genomics by providing the first homosporous fern reference genome, as well as insights into the processes underlying the formation of these massive genomes.

Worldwide relationships in the fern genus Pteridium (bracken) based on nuclear genome markers.

PREMISE: Spore-bearing plants are capable of dispersing very long distances. However, it is not known if gene flow can prevent genetic divergence in widely distributed taxa. Here we address this issue, and examine systematic relationships at a global geographic scale for the fern genus Pteridium. METHODS: We sampled plants from 100 localities worldwide, and generated nucleotide data from four nuclear genes and two plastid regions. We also examined 2801 single nucleotide polymorphisms detected by a restriction site-associated DNA approach. RESULTS: We found evidence for two distinct diploid species and two allotetraploids between them. The "northern" species (Pteridium aquilinum) has distinct groups at the continental scale (Europe, Asia, Africa, and North America). The northern European subspecies pinetorum appears to involve admixture among all of these. A sample from the Hawaiian Islands contained elements of both North American and Asian P. aquilinum. The "southern" species, P. esculentum, shows little genetic differentiation between South American and Australian samples. Components of African genotypes are detected on all continents. CONCLUSIONS: We find evidence of distinct continental-scale genetic differentiation in Pteridium. However, on top of this is a clear signal of recent hybridization. Thus, spore-bearing plants are clearly capable of extensive long-distance gene flow; yet appear to have differentiated genetically at the continental scale. Either gene flow in the past was at a reduced level, or vicariance is possible even in the face of long-distance gene flow.

ReFernment: An R package for annotating RNA editing in plastid genomes.

PREMISE OF THE STUDY: In the absence of cDNA, the annotation of RNA editing in plastomes must be done manually, representing a significant time cost to those studying the organellar genomes of ferns and hornworts. METHODS AND RESULTS: We developed an R package to automatically annotate apparent nonsense mutations in plastid genomes. The software successfully annotates such sites and results in no false positives for data with no sequencing or assembly errors. CONCLUSIONS: Compared to manual annotation, ReFernment offers greater speed and accuracy for annotating RNA editing sites. This software should be especially useful for researchers generating large numbers of plastome sequences for taxa with high levels of RNA editing.

Mobile Elements Shape Plastome Evolution in Ferns.

Plastid genomes display remarkable organizational stability over evolutionary time. From green algae to angiosperms, most plastid genomes are largely collinear, with only a few cases of inversion, gene loss, or, in extremely rare cases, gene addition. These plastome insertions are mostly clade-specific and are typically of nuclear or mitochondrial origin. Here, we expand on these findings and present the first family-level survey of plastome evolution in ferns, revealing a novel suite of dynamic mobile elements. Comparative plastome analyses of the Pteridaceae expose several mobile open reading frames that vary in sequence length, insertion site, and configuration among sampled taxa. Even between close relatives, the presence and location of these elements is widely variable when viewed in a phylogenetic context. We characterize these elements and refer to them collectively as Mobile Open Reading Frames in Fern Organelles (MORFFO). We further note that the presence of MORFFO is not restricted to Pteridaceae, but is found across ferns and other plant clades. MORFFO elements are regularly associated with inversions, intergenic expansions, and changes to the inverted repeats. They likewise appear to be present in mitochondrial and nuclear genomes of ferns, indicating that they can move between genomic compartments with relative ease. The origins and functions of these mobile elements are unknown, but MORFFO appears to be a major driver of structural genome evolution in the plastomes of ferns, and possibly other groups of plants.

Target sequence capture of nuclear-encoded genes for phylogenetic analysis in ferns.

PREMISE OF THE STUDY: Until recently, most phylogenetic studies of ferns were based on chloroplast genes. Evolutionary inferences based on these data can be incomplete because the characters are from a single linkage group and are uniparentally inherited. These limitations are particularly acute in studies of hybridization, which is prevalent in ferns; fern hybrids are common and ferns are able to hybridize across highly diverged lineages, up to 60 million years since divergence in one documented case. However, it not yet clear what effect such hybridization has on fern evolution, in part due to a paucity of available biparentally inherited (nuclear-encoded) markers. METHODS: We designed oligonucleotide baits to capture 25 targeted, low-copy nuclear markers from a sample of 24 species spanning extant fern diversity. RESULTS: Most loci were successfully sequenced from most accessions. Although the baits were designed from exon (transcript) data, we successfully captured intron sequences that should be useful for more focused phylogenetic studies. We present phylogenetic analyses of the new target sequence capture data and integrate these into a previous transcript-based data set. DISCUSSION: We make our bait sequences available to the community as a resource for further studies of fern phylogeny.

Fern genomes elucidate land plant evolution and cyanobacterial symbioses.

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade.

BACKGROUND: Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale. RESULTS: In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species. CONCLUSIONS: We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.

Evolutionary Conservation of ABA Signaling for Stomatal Closure.

Abscisic acid (ABA)-driven stomatal regulation reportedly evolved after the divergence of ferns, during the early evolution of seed plants approximately 360 million years ago. This hypothesis is based on the observation that the stomata of certain fern species are unresponsive to ABA, but exhibit passive hydraulic control. However, ABA-induced stomatal closure was detected in some mosses and lycophytes. Here, we observed that a number of ABA signaling and membrane transporter protein families diversified over the evolutionary history of land plants. The aquatic ferns Azolla filiculoides and Salvinia cucullata have representatives of 23 families of proteins orthologous to those of Arabidopsis (Arabidopsis thaliana) and all other land plant species studied. Phylogenetic analysis of the key ABA signaling proteins indicates an evolutionarily conserved stomatal response to ABA. Moreover, comparative transcriptomic analysis has identified a suite of ABA-responsive genes that differentially expressed in a terrestrial fern species, Polystichum proliferum These genes encode proteins associated with ABA biosynthesis, transport, reception, transcription, signaling, and ion and sugar transport, which fit the general ABA signaling pathway constructed from Arabidopsis and Hordeum vulgare The retention of these key ABA-responsive genes could have had a profound effect on the adaptation of ferns to dry conditions. Furthermore, stomatal assays have shown the primary evidence for ABA-induced closure of stomata in two terrestrial fern species Pproliferum and Nephrolepis exaltata In summary, we report, to our knowledge, new molecular and physiological evidence for the presence of active stomatal control in ferns.

Origins and diversity of a cosmopolitan fern genus on an island archipelago.

Isolated oceanic islands are characterized by patterns of biological diversity different from that on nearby continental mainlands. Isolation can provide the opportunity for evolutionary divergence, but also set the stage for hybridization between related taxa arriving from different sources. Ferns disperse by haploid spores, which are produced in large numbers and can travel long distances in air currents, enabling these plants to become established on most oceanic islands. Here, we examine the origins and patterns of diversity of the cosmopolitan fern genus Pteridium (Dennstaedtiaceae; bracken) on the Galapagos Islands. We use nucleotide sequences from two plastid genes, and two nuclear gene markers, to examine phylogeography of Pteridium on the Galapagos Islands. We incorporate data from a previous study to provide a worldwide context. We also sampled new specimens from South and Central America. We used flow cytometry to estimate genome size of some accessions. We found that both plastid and nuclear haplotypes fall into two distinct clades, consistent with a two-diploid-species taxonomy of P. aquilinum and P. esculentum. As predicted, the allotetraploid P. caudatum possesses nuclear haplotypes from both diploid species. Samples from the Galapagos include P. esculentum subsp. arachnoideum, P. caudatum and possible hybrids between them. Multiple Pteridium taxa were also observed growing together at some sites. We find evidence for multiple origins of Pteridium on the Galapagos Islands and multiple origins of tetraploid P. caudatum throughout its range in Central and South America. We also posit that P. caudatum may include recent diploid hybrids, backcrosses to P. esculentum, as well as allotetraploid plants. The Galapagos Islands are positioned close to the equator where they can receive dispersing propagules from both hemispheres. This may partly explain the high levels of diversity found for this cosmopolitan fern on these islands.

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.

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.

Genotyping-by-sequencing for Populus population genomics: an assessment of genome sampling patterns and filtering approaches.

Continuing advances in nucleotide sequencing technology are inspiring a suite of genomic approaches in studies of natural populations. Researchers are faced with data management and analytical scales that are increasing by orders of magnitude. With such dramatic advances comes a need to understand biases and error rates, which can be propagated and magnified in large-scale data acquisition and processing. Here we assess genomic sampling biases and the effects of various population-level data filtering strategies in a genotyping-by-sequencing (GBS) protocol. We focus on data from two species of Populus, because this genus has a relatively small genome and is emerging as a target for population genomic studies. We estimate the proportions and patterns of genomic sampling by examining the Populus trichocarpa genome (Nisqually-1), and demonstrate a pronounced bias towards coding regions when using the methylation-sensitive ApeKI restriction enzyme in this species. Using population-level data from a closely related species (P. tremuloides), we also investigate various approaches for filtering GBS data to retain high-depth, informative SNPs that can be used for population genetic analyses. We find a data filter that includes the designation of ambiguous alleles resulted in metrics of population structure and Hardy-Weinberg equilibrium that were most consistent with previous studies of the same populations based on other genetic markers. Analyses of the filtered data (27,910 SNPs) also resulted in patterns of heterozygosity and population structure similar to a previous study using microsatellites. Our application demonstrates that technically and analytically simple approaches can readily be developed for population genomics of natural populations.

Widespread triploidy in Western North American aspen (Populus tremuloides).

We document high rates of triploidy in aspen (Populus tremuloides) across the western USA (up to 69% of genets), and ask whether the incidence of triploidy across the species range corresponds with latitude, glacial history (as has been documented in other species), climate, or regional variance in clone size. Using a combination of microsatellite genotyping, flow cytometry, and cytology, we demonstrate that triploidy is highest in unglaciated, drought-prone regions of North America, where the largest clone sizes have been reported for this species. While we cannot completely rule out a low incidence of undetected aneuploidy, tetraploidy or duplicated loci, our evidence suggests that these phenomena are unlikely to be significant contributors to our observed patterns. We suggest that the distribution of triploid aspen is due to a positive synergy between triploidy and ecological factors driving clonality. Although triploids are expected to have low fertility, they are hypothesized to be an evolutionary link to sexual tetraploidy. Thus, interactions between clonality and polyploidy may be a broadly important component of geographic speciation patterns in perennial plants. Further, cytotypes are expected to show physiological and structural differences which may influence susceptibility to ecological factors such as drought, and we suggest that cytotype may be a significant and previously overlooked factor in recent patterns of high aspen mortality in the southwestern portion of the species range. Finally, triploidy should be carefully considered as a source of variance in genomic and ecological studies of aspen, particularly in western U.S. landscapes.

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.

The evolution of chloroplast genes and genomes in ferns.

Most of the publicly available data on chloroplast (plastid) genes and genomes come from seed plants, with relatively little information from their sister group, the ferns. Here we describe several broad evolutionary patterns and processes in fern plastid genomes (plastomes), and we include some new plastome sequence data. We review what we know about the evolutionary history of plastome structure across the fern phylogeny and we compare plastome organization and patterns of evolution in ferns to those in seed plants. A large clade of ferns is characterized by a plastome that has been reorganized with respect to the ancestral gene order (a similar order that is ancestral in seed plants). We review the sequence of inversions that gave rise to this organization. We also explore global nucleotide substitution patterns in ferns versus those found in seed plants across plastid genes, and we review the high levels of RNA editing observed in fern plastomes.