Exploring the phylogeny of the marattialean ferns

The eusporangiate marattialean ferns represent an ancient radiation with a rich fossil record but limited modern diversity in the tropics. The long evolutionary history without close extant relatives has confounded studies of the phylogenetic origin, rooting and timing of marattialean ferns. Here we present new complete plastid genomes of six marattialean species and compiled a plastid genome dataset representing all of the currently accepted marattialean genera. We further supplemented this dataset by compiling a large dataset of mitochondrial genes and a phenotypic data matrix covering both extant and extinct representatives of the lineage. Our phylogenomic and total-evidence analyses corroborated the postulated position of marattialean ferns as the sister to leptosporangiate ferns, and the position of Danaea as the sister to the remaining extant marattialean genera. However, our results provide new evidence that Christensenia is sister to Marattia and that M. cicutifolia actually belongs to Eupodium. The apparently highly reduced rate of molecular evolution in marattialean ferns provides a challenge for dating the key phylogenetic events with molecular clock approaches. We instead applied a parsimony-based total-evidence dating approach, which suggested a Triassic age for the extant crown group. The modern distribution can best be explained as mainly resulting from vicariance following the breakup of Pangaea and Gondwana. We resolved the fossil genera Marattiopsis, Danaeopsis and Qasimia as members of the monophyletic family Marattiaceae, and the Carboniferous genera Sydneia and Radstockia as the monophyletic sister of all other marattialean ferns.     

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions

Background

Ferns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes.

Results

We characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN).

Conclusions

IRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

     

Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails

Using DNA sequence data from multiple genes (often from more than one genome compartment) to reconstruct phylogenetic relationships has become routine. Augmenting this approach with genomic structural characters (e.g., intron gain and loss, changes in gene order) as these data become available from comparative studies already has provided critical insight into some long-standing questions about the evolution of land plants. Here we report on the presence of a group II intron located in the mitochondrial atp1 gene of leptosporangiate and marattioid ferns. Primary sequence data for the atp1 gene are newly reported for 27 taxa, and results are presented from maximum likelihood-based phylogenetic analyses using Bayesian inference for 34 land plants in three data sets: (1) single-gene mitochondrial atp1 (exon+intron sequences); (2) five combined genes (mitochondrial atp1 [exon only]; plastid rbcL, atpB, rps4; nuclear SSU rDNA); and (3) same five combined genes plus morphology. All our phylogenetic analyses corroborate results from previous fern studies that used plastid and nuclear sequence data: the monophyly of euphyllophytes, as well as of monilophytes; whisk ferns (Psilotidae) sister to ophioglossoid ferns (Ophioglossidae); horsetails (Equisetopsida) sister to marattioid ferns (Marattiidae), which together are sister to the monophyletic leptosporangiate ferns. In contrast to the results from the primary sequence data, the genomic structural data (atp1 intron distribution pattern) would seem to suggest that leptosporangiate and marattioid ferns are monophyletic, and together they are the sister group to horsetails--a topology that is rarely reconstructed using primary sequence data.     

Evolutionary relationships within the Neotropical, eusporangiate fern genus Danaea (Marattiaceae)

Genera within the eusporangiate fern family Marattiaceae have long been neglected in taxonomic and systematic studies. Here we present the first phylogenetic hypothesis of relationships within the exclusively Neotropical genus Danaea based on a sampling of 60 specimens representing 31 species from various Neotropical sites. We used DNA sequence data from three plastid regions (atpB, rbcL, and trnL-F), morphological characters from both herbarium specimens and live plants observed in the field, and geographical and ecological information to examine evolutionary patterns. Eleven representatives of five other marattioid genera (Angiopteris, Archangiopteris, Christensenia, Macroglossum, and Marattia) were used to root the topology. We identified three well-supported clades within Danaea that are consistent with morphological characters: the "leprieurii" clade (containing species traditionally associated with the name D. elliptica), the "nodosa" clade (containing all species traditionally associated with the name D. nodosa), and the "alata" clade (containing all other species). All three clades are geographically and ecologically widely distributed, but subclades within them show various distribution patterns. Our phylogenetic hypothesis provides a robust framework within which broad questions related to the morphology, taxonomy, biogeography, evolution, and ecology of these ferns can be addressed.     

Complete chloroplast genome sequence of a tree fern Alsophila spinulosa : insights into evolutionary changes in fern chloroplast genomes

Background  

Ferns have generally been neglected in studies of chloroplast genomics. Before this study, only one polypod and two basal ferns had their complete chloroplast (cp) genome reported. Tree ferns represent an ancient fern lineage that first occurred in the Late Triassic. In recent phylogenetic analyses, tree ferns were shown to be the sister group of polypods, the most diverse group of living ferns. Availability of cp genome sequence from a tree fern will facilitate interpretation of the evolutionary changes of fern cp genomes. Here we have sequenced the complete cp genome of a scaly tree fern Alsophila spinulosa (Cyatheaceae).     

Phylogenetics and classification of the pantropical fern family Lindsaeaceae

The classification and generic definition in the tropical–subtropical fern family Lindsaeaceae have been uncertain and have so far been based on morphological characters only. We have now studied the evolutionary history of the Lindsaeaceae by simultaneously optimizing 55 morphological characters, two plastid genes (rpoC1 and rps4) and three non‐coding plastid intergenic spacers (trnL‐F, rps4‐trnS and trnH‐psbA). Our data set included all genera associated with Lindsaeaceae, except Xyropteris, and c. 73% of the currently accepted species. The phylogenetic relationships of the lindsaeoid ferns with two enigmatic genera that have recently been included in the Lindsaeaceae, Cystodium and Lonchitis, remain ambiguous. Within the monophyletic lindsaeoids, we found six well‐supported and diagnostic clades that can be recognized as genera: Sphenomeris, Odontosoria, Osmolindsaea, Nesolindsaea, Tapeinidium and Lindsaea. Sphenomeris was shown to be monotypic; most taxa formerly placed in that genus belong to the Odontosoria clade. Ormoloma is embedded within Lindsaea and therefore does not merit recognition as a genus. Tapeinidium is sister to a clade with some species formerly placed in Lindsaea that are morphologically distinct from that genus and are transferred to Osmolindsaea and Nesolindsaea, proposed here as two new genera. We do not maintain the current subgeneric classification of Lindsaea itself, because neither of the two generally accepted subgenera (Lindsaea and Odontoloma) is monophyletic, and most of the sections also appear unnatural. Nesolindsaea shows an ancient biogeographical link between Sri Lanka and the Seychelles and many of the main clades within Lindsaea have geographically disjunct distributions. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 163 , 305–359.     

A global plastid phylogeny of the fern genus Asplenium (Aspleniaceae)

The infrageneric relationships and taxonomy of the largest fern genus, Asplenium (Aspleniaceae), have remained poorly understood. Previous studies have focused mainly on specific species complexes involving a few or dozens of species only, or have achieved a large taxon sampling but only one plastid marker was used. In the present study, DNA sequences from six plastid markers (atpB, rbcL, rps4, rps4-trnS, trnL and trnL-F) of 1030 accessions (616 of them newly sequenced here) representing c. 420 species of Asplenium (60% of estimated species diversity), 16 species of Hymenasplenium, three Diplaziopsidaceae, and four Rhachidosoraceae were used to produce the largest genus-level phylogeny yet for ferns. Our major results include: (i) Asplenium as broadly circumscribed is monophyletic based on our inclusion of representatives of 32 of 38 named segregate genera; (ii) 11 major clades in Asplenium are identified, and their relationships are mostly well-resolved and strongly supported; (iii) numerous species, unsampled in previous studies, suggest new relationships and numerous cryptic species and species complexes in Asplenium; and (iv) the accrued molecular evidence provides an essential foundation for further investigations of complex patterns of geographical diversification, speciation and reticulate evolution in this family.     

First insights into fern matK phylogeny

MatK, the only maturase gene in the land plant plastid genome, is a very popular phylogenetic marker that has been extensively applied in reconstructing angiosperm phylogeny. However, the use of matK in fern phylogeny is largely unknown, due to difficulties with amplification: ferns have lost the flanking trnK exons, typically the region used for designing stable priming sites. We developed primers that are either universal or lineage-specific that successfully amplify matK across all fern families. To evaluate whether matK is as powerful a phylogenetic marker in ferns as in angiosperms, we compared its sequence characteristics and phylogenetic performance to those of rbcL and atpA. Among these three genes, matK has the highest variability and substitution evenness, yet shows the least homoplasy. Most importantly, applying matK in fern phylogenetics better resolved relationships among families, especially within eupolypods I and II. Here we demonstrate the power of matK for fern phylogenetic reconstruction, as well as provide primers and extensive sequence data that will greatly facilitate future evolutionary studies of ferns.     

Plastid genome sequencing,comparative genomics,and phylogenomics: Current status and prospects

Abstract More than 190 plastid genomes have been completely sequenced during the past two decades due to advances in DNA sequencing technologies. Based on this unprecedented abundance of data, extensive genomic changes have been revealed in the plastid genomes. Inversion is the most common mechanism that leads to gene order changes. Several inversion events have been recognized as informative phylogenetic markers, such as a 30‐kb inversion found in all living vascular plants minus lycopsids and two short inversions putatively shared by all ferns. Gene loss is a common event throughout plastid genome evolution. Many genes were independently lost or transferred to the nuclear genome in multiple plant lineages. The trnR‐CCG gene was lost in some clades of lycophytes, ferns, and seed plants, and all the ndh genes were absent in parasitic plants, gnetophytes, Pinaceae, and the Taiwan moth orchid. Certain parasitic plants have, in particular, lost plastid genes related to photosynthesis because of the relaxation of functional constraint. The dramatic growth of plastid genome sequences has also promoted the use of whole plastid sequences and genomic features to solve phylogenetic problems. Chloroplast phylogenomics has provided additional evidence for deep‐level phylogenetic relationships as well as increased phylogenetic resolutions at low taxonomic levels. However, chloroplast phylogenomics is still in its infant stage and rigorous analysis methodology has yet to be developed.     

Comparative Analysis of the Complete Plastid Genome Sequence of the Red Alga Gracilaria tenuistipitata var. liui Provides Insights into the Evolution of Rhodoplasts and Their Relationship to Other Plastids

We sequenced to completion the circular plastid genome of the red alga Gracilaria tenuistipitata var. liui. This is the first plastid genome sequence from the subclass Florideophycidae (Rhodophyta). The genome is composed of 183,883 bp and contains 238 predicted genes, including a single copy of the ribosomal RNA operon. Comparisons with the plastid genome of Porphyra pupurea reveal strong conservation of gene content and order, but we found major genomic rearrangements and the presence of coding regions that are specific to Gracilaria. Phylogenetic analysis of a data set of 41 concatenated proteins from 23 plastid and two cyanobacterial genomes support red algal plastid monophyly and a specific evolutionary relationship between the Florideophycidae and the Bangiales. Gracilaria maintains a surprisingly ancient gene content in its plastid genome and, together with other Rhodophyta, contains the most complete repertoire of plastid genes known in photosynthetic eukaryotes.Supplementary material () is available for this article.[Reviewing Editor: Dr. W. Ford Doolittle]     

Volatile organic compounds in the strongly fragrant fern genus Melpomene (Polypodiaceae)

Volatile organic compounds (VOCs) are common among plants, both as attractants for pollinators and as defence against herbivores. While much studied among flowering plants, the prevalence and function of VOCs among ferns is little known. Using headspace sorption and gas chromatography, we analysed the VOCs of dried specimens of six species of grammitid fern (Polypodiaceae), including two species of the genus Melpomene, which is characterised by a distinctive sweet smell. We identified 38 VOCs, including 22 not previously recorded among ferns. The two species of Melpomene had distinct VOC cocktails, including 12 substances not found in the other four studied genera, mainly involving fatty acid derivatives (FADs) and aromatics. We propose that these VOCs have, at least in part, a function in herbivore defence, but note that the VOC bouquet of Melpomene is distinct from that typically found in angiosperms.     

Spatial genetic analyses reveal strong genetic structure in two populations of the outcrossing tree fern Alsophila firma (Cyatheaceae)

The development of spatial genetic structure (SGS) in seed plants has been linked to several biological attributes of species, such as breeding system and life form. However, little is known about SGS in ferns, which together with lycopods are unique among land plants in having two free‐living life stages. We combined spatial aggregation statistics and spatial genetic autocorrelation analyses using five plastid microsatellites and one nuclear gene to investigate SGS in two populations of the outcrossing tree fern Alsophila firma (Cyatheaceae). We assessed how the observed patterns compare with those estimated for other ferns and seed plants. Populations of A. firma exhibited strong SGS, spatial clustering of individuals, substantial clonal diversity and no inbreeding. SGS in ferns appears to be higher than in most seed plants analysed to date. Contrary to our expectations, an outcrossing breeding system, wind dispersal and an arborescent life form did not translate into weak or no SGS. In ferns, SGS is probably being affected by the life cycle with two free‐living life stages. The reproductive biology of ferns appears to be more complex than previously thought. This implies that SGS in ferns is affected by some factors that cannot be inferred from the study of flowering plants. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 439–449.     

Gene-rich plastid genomes of two parasitic red algal species,Laurencia australis and L. verruciformis (Rhodomelaceae,Ceramiales), and a taxonomic revision of Janczewskia

Parasitic red algae are an interesting system for investigating the genetic changes that occur in parasites. These parasites have evolved independently multiple times within the red algae. The functional loss of plastid genomes can be investigated in these multiple independent examples, and fine-scale patterns may be discerned. The only plastid genomes from red algal parasites known so far are highly reduced and missing almost all photosynthetic genes. Our study assembled and annotated plastid genomes from the parasites Janczewskia tasmanica and its two Laurencia host species (Laurencia elata and one unidentified Laurencia sp. A25) from Australia and Janczewskia verruciformis, its host species (Laurencia catarinensis), and the closest known free-living relative (Laurencia obtusa) from the Canary Islands (Spain). For the first time we show parasitic red algal plastid genomes that are similar in size and gene content to free-living host species without any gene loss or genome reduction. The only exception was two pseudogenes (moeB and ycf46) found in the plastid genome of both isolates of J. tasmanica, indicating potential for future loss of these genes. Further comparative analyses with the three highly reduced plastid genomes showed possible gene loss patterns, in which photosynthetic gene categories were lost followed by other gene categories. Phylogenetic analyses did not confirm monophyly of Janczewskia, and the genus was subsumed into Laurencia. Further investigations will determine if any convergent small-scale patterns of gene loss exist in parasitic red algae and how these are applicable to other parasitic systems.     

Molecular phylogeny and species delimitation of Stachyuraceae: Advocating a herbarium specimen‐based phylogenomic approach in resolving species boundaries

Species concept and delimitation are fundamental to taxonomic and evolutionary studies. Both inadequate informative sites in the molecular data and limited taxon sampling have often led to poor phylogenetic resolution and incorrect species delineation. Recently, the whole chloroplast genome sequences from extensive herbarium specimen samples have been shown to be effective to amend the problem. Stachyuraceae are a small family consisting of only one genus Stachyurus of six to 16 species. However, species delimitation in Stachyurus has been highly controversial because of few and generally unstable morphological characters used for classification. In this study, we sampled 69 individuals of seven species (each with at least three individuals) covering the entire taxonomic diversity, geographic range, and morphological variation of Stachyurus from herbarium specimens for genome‐wide plastid gene sequencing to address species delineation in the genus. We obtained high‐quality DNAs from specimens using a recently developed DNA reconstruction technique. We first assembled four whole chloroplast genome sequences. Based on the chloroplast genome and one nuclear ribosomal DNA sequence of Stachyurus, we designed primers for multiplex polymerase chain reaction and high throughput sequencing of 44 plastid loci for species of Stachyurus. Data of these chloroplast DNA and nuclear ribosomal DNA internal transcribed spacer sequences were used for phylogenetic analyses. The phylogenetic results showed that the Japanese species Stachyurus praecox Siebold & Zucc. was sister to the rest in mainland China, which indicated a typical Sino‐Japanese distribution pattern. Based on diagnostic morphological characters, distinct distributional range, and monophyly of each clade, we redefined seven species for Stachyurus following an integrative species concept, and revised the taxonomy of the family based on previous reports and specimens, in particular the type specimens. Furthermore, our divergence time estimation results suggested that Stachyuraceae split from its sister group Crossosomataceae from the New World at ca. 54.29 Mya, but extant species of Stachyuraceae started their diversification only recently at ca. 6.85 Mya. Diversification time of Stachyurus in mainland China was estimated to be ca. 4.45 Mya. This research has provided an example of using the herbarium specimen‐based phylogenomic approach in resolving species boundaries in a taxonomically difficult genus.     

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron

Standard plant DNA barcodes based on 2–3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3–4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.     

Complete DNA sequences of the plastid genomes of two parasitic flowering plant species, <Emphasis Type="Italic">Cuscuta reflexa</Emphasis> and <Emphasis Type="Italic">Cuscuta gronovii</Emphasis>

Background  

The holoparasitic plant genus Cuscuta comprises species with photosynthetic capacity and functional chloroplasts as well as achlorophyllous and intermediate forms with restricted photosynthetic activity and degenerated chloroplasts. Previous data indicated significant differences with respect to the plastid genome coding capacity in different Cuscuta species that could correlate with their photosynthetic activity. In order to shed light on the molecular changes accompanying the parasitic lifestyle, we sequenced the plastid chromosomes of the two species Cuscuta reflexa and Cuscuta gronovii. Both species are capable of performing photosynthesis, albeit with varying efficiencies. Together with the plastid genome of Epifagus virginiana, an achlorophyllous parasitic plant whose plastid genome has been sequenced, these species represent a series of progression towards total dependency on the host plant, ranging from reduced levels of photosynthesis in C. reflexa to a restricted photosynthetic activity and degenerated chloroplasts in C. gronovii to an achlorophyllous state in E. virginiana.     

Complete Plastid Genome Sequencing of Four Tilia Species (Malvaceae): A Comparative Analysis and Phylogenetic Implications

Tilia is an ecologically and economically important genus in the family Malvaceae. However, there is no complete plastid genome of Tilia sequenced to date, and the taxonomy of Tilia is difficult owing to frequent hybridization and polyploidization. A well-supported interspecific relationships of this genus is not available due to limited informative sites from the commonly used molecular markers. We report here the complete plastid genome sequences of four Tilia species determined by the Illumina technology. The Tilia plastid genome is 162,653 bp to 162,796 bp in length, encoding 113 unique genes and a total number of 130 genes. The gene order and organization of the Tilia plastid genome exhibits the general structure of angiosperms and is very similar to other published plastid genomes of Malvaceae. As other long-lived tree genera, the sequence divergence among the four Tilia plastid genomes is very low. And we analyzed the nucleotide substitution patterns and the evolution of insertions and deletions in the Tilia plastid genomes. Finally, we build a phylogeny of the four sampled Tilia species with high supports using plastid phylogenomics, suggesting that it is an efficient way to resolve the phylogenetic relationships of this genus.