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

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

Order-level fern plastome phylogenomics: new insights from Hymenophyllales

PREMISE OF THE STUDY

Filmy ferns (Hymenophyllales) are a highly specialized lineage, having mesophyll one-cell layer thick and inhabiting particularly shaded and humid environments. The phylogenetic placement of Hymenophyllales has been inconclusive, and while over 87 whole fern plastomes have been published, none was from Hymenophyllales. To better understand the evolutionary history of filmy ferns, we sequenced the first complete plastome for this order.

METHODS

We compiled a phylogenomic plastome data set encompassing all 11 fern orders, and reconstructed phylogenies using different data types (nucleotides, codons, and amino acids) and partition schemes (codon positions and loci). To infer the evolution of fern plastome organization, we coded plastome features, including inversions, inverted repeat boundary shifts, gene losses, and tRNA anticodon sequences as characters, and reconstructed the ancestral states for these characters.

KEY RESULTS

We discovered a suite of novel, Hymenophyllales-specific plastome structures that likely resulted from repeated expansions and contractions of the inverted repeat regions. Our phylogenetic analyses reveal that Hymenophyllales is highly supported as either sister to Gleicheniales or to Gleicheniales + the remaining non-Osmundales leptosporangiates, depending on the data type and partition scheme.

CONCLUSIONS

Although our analyses could not confidently resolve the phylogenetic position of Hymenophyalles, the results here highlight the danger of drawing conclusions from "all-in" phylogenomic data set without exploring potential inconsistencies in the data. Finally, our first order-level reconstruction of fern plastome structural evolution provides a useful framework for future plastome research.

Origin of Equisetum: Evolution of horsetails (Equisetales) within the major euphyllophyte clade Sphenopsida

PREMISE OF THE STUDY

Equisetum is the sole living representative of Sphenopsida, a clade with impressive species richness, a long fossil history dating back to the Devonian, and obscure relationships with other living pteridophytes. Based on molecular data, the crown group age of Equisetum is mid-Paleogene, although fossils with possible crown synapomorphies appear in the Triassic. The most widely circulated hypothesis states that the lineage of Equisetum derives from calamitaceans, but no comprehensive phylogenetic studies support the claim. Using a combined approach, we provide a comprehensive phylogenetic analysis of Equisetales, with special emphasis on the origin of genus Equisetum.

METHODS

We performed parsimony phylogenetic analyses to address relationships of 43 equisetalean species (15 extant, 28 extinct) using a combination of morphological and molecular characters.

KEY RESULTS

We recovered Equisetaceae + Neocalamites as sister to Calamitaceae + a clade of Angaran and Gondwanan horsetails, with the four groups forming a clade that is sister to Archaeocalamitaceae. The estimated age for the Equisetum crown group is mid-Mesozoic.

CONCLUSIONS

Modern horsetails are not nested within calamitaceans; instead, both groups have explored independent evolutionary trajectories since the Carboniferous. Diverse fossil taxon sampling helps to shed light on the position and relationships of equisetalean lineages, of which only a tiny remnant is present within the extant flora. Understanding these relationships and early character configurations of ancient plant clades as Equisetales provide useful tests of hypotheses about overall phylogenetic relationships of euphyllophytes and foundations for future tests of molecular dates with paleontological data.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.

Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralog in rpl2

Intron patterns in plant mitochondrial genomes differ significantly between the major land plant clades. We here report on a new, clade-specific group II intron in the rps1 gene of monilophytes (ferns). This intron, rps1i25g2, is strikingly similar to rpl2i846g2 previously identified in the mitochondrial rpl2 gene of seed plants, ferns, and the lycophyte Phlegmariurus squarrosus Although mitochondrial ribosomal protein genes are frequently subject to endosymbiotic gene transfer among plants, we could retrieve the mitochondrial rps1 gene in a taxonomically wide sampling of 44 monilophyte taxa including basal lineages such as the Ophioglossales, Psilotales, and Marattiales with the only exception being the Equisetales (horsetails). Introns rps1i25g2 and rpl2i846g2 were likewise consistently present with only two exceptions: Intron rps1i25g2 is lost in the genus Ophioglossum and intron rpl2i846g2 is lost in Equisetum bogotense Both intron sequences are moderately affected by RNA editing. The unprecedented primary and secondary structure similarity of rps1i25g2 and rpl2i846g2 suggests an ancient retrotransposition event copying rpl2i846g2 into rps1, for which we suggest a model. Our phylogenetic analysis adding the new rps1 locus to a previous data set is fully congruent with recent insights on monilophyte phylogeny and further supports a sister relationship of Gleicheniales and Hymenophyllales.

Convergent Evolution of Fern-Specific Mitochondrial Group II Intron atp1i361g2 and Its Ancient Source Paralogue rps3i249g2 and Independent Losses of Intron and RNA Editing among Pteridaceae

Mitochondrial intron patterns are highly divergent between the major land plant clades. An intron in the atp1 gene, atp1i361g2, is an example for a group II intron specific to monilophytes (ferns). Here, we report that atp1i361g2 is lost independently at least 4 times in the fern family Pteridaceae. Such plant organelle intron losses have previously been found to be accompanied by loss of RNA editing sites in the flanking exon regions as a consequence of genomic recombination of mature cDNA. Instead, we now observe that RNA editing events in both directions of pyrimidine exchange (C-to-U and U-to-C) are retained in atp1 exons after loss of the intron in Pteris argyraea/biaurita and in Actiniopteris and Onychium We find that atp1i361g2 has significant similarity with intron rps3i249g2 present in lycophytes and gymnosperms, which we now also find highly conserved in ferns. We conclude that atp1i361g2 may have originated from the more ancestral rps3i249g2 paralogue by a reverse splicing copy event early in the evolution of monilophytes. Secondary structure elements of the two introns, most characteristically their domains III, show strikingly convergent evolution in the monilophytes. Moreover, the intron paralogue rps3i249g2 reveals relaxed evolution in taxa where the atp1i361g2 paralogue is lost. Our findings may reflect convergent evolution of the two related mitochondrial introns exerted by co-evolution with an intron-binding protein simultaneously acting on the two paralogues.

Reverse U-to-C editing exceeds C-to-U RNA editing in some ferns - a monilophyte-wide comparison of chloroplast and mitochondrial RNA editing suggests independent evolution of the two processes in both organelles

BACKGROUND

RNA editing by C-to-U conversions is nearly omnipresent in land plant chloroplasts and mitochondria, where it mainly serves to reconstitute conserved codon identities in the organelle mRNAs. Reverse U-to-C RNA editing in contrast appears to be restricted to hornworts, some lycophytes, and ferns (monilophytes). A well-resolved monilophyte phylogeny has recently emerged and now allows to trace the side-by-side evolution of both types of pyrimidine exchange editing in the two endosymbiotic organelles.

RESULTS

Our study of RNA editing in four selected mitochondrial genes show a wide spectrum of divergent RNA editing frequencies including a dominance of U-to-C over the canonical C-to-U editing in some taxa like the order Schizaeales. We find that silent RNA editing leaving encoded amino acids unchanged is highly biased with more than ten-fold amounts of silent C-to-U over U-to-C edits. In full contrast to flowering plants, RNA editing frequencies are low in early-branching monilophyte lineages but increase in later emerging clades. Moreover, while editing rates in the two organelles are usually correlated, we observe uncoupled evolution of editing frequencies in fern mitochondria and chloroplasts. Most mitochondrial RNA editing sites are shared between the recently emerging fern orders whereas chloroplast editing sites are mostly clade-specific. Finally, we observe that chloroplast RNA editing appears to be completely absent in horsetails (Equisetales), the sister clade of all other monilophytes.

CONCLUSIONS

C-to-U and U-to-C RNA editing in fern chloroplasts and mitochondria follow disinct evolutionary pathways that are surprisingly different from what has previously been found in flowering plants. The results call for careful differentiation of the two types of RNA editing in the two endosymbiotic organelles in comparative evolutionary studies.

The evolutionary history of ferns inferred from 25 low-copy nuclear genes

•

PREMISE OF THE STUDY

Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.•

METHODS

Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.•

KEY RESULTS

Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.•

CONCLUSIONS

Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.

Stasis and convergence characterize morphological evolution in eupolypod II ferns

BACKGROUND AND AIMS

Patterns of morphological evolution at levels above family rank remain underexplored in the ferns. The present study seeks to address this gap through analysis of 79 morphological characters for 81 taxa, including representatives of all ten families of eupolypod II ferns. Recent molecular phylogenetic studies demonstrate that the evolution of the large eupolypod II clade (which includes nearly one-third of extant fern species) features unexpected patterns. The traditional 'athyrioid' ferns are scattered across the phylogeny despite their apparent morphological cohesiveness, and mixed among these seemingly conservative taxa are morphologically dissimilar groups that lack any obvious features uniting them with their relatives. Maximum-likelihood and maximum-parsimony character optimizations are used to determine characters that unite the seemingly disparate groups, and to test whether the polyphyly of the traditional athyrioid ferns is due to evolutionary stasis (symplesiomorphy) or convergent evolution. The major events in eupolypod II character evolution are reviewed, and character and character state concepts are reappraised, as a basis for further inquiries into fern morphology.

METHODS

Characters were scored from the literature, live plants and herbarium specimens, and optimized using maximum-parsimony and maximum-likelihood, onto a highly supported topology derived from maximum-likelihood and Bayesian analysis of molecular data. Phylogenetic signal of characters were tested for using randomization methods and fitdiscrete.

KEY RESULTS

The majority of character state changes within the eupolypod II phylogeny occur at the family level or above. Relative branch lengths for the morphological data resemble those from molecular data and fit an ancient rapid radiation model (long branches subtended by very short backbone internodes), with few characters uniting the morphologically disparate clades. The traditional athyrioid ferns were circumscribed based upon a combination of symplesiomorphic and homoplastic characters. Petiole vasculature consisting of two bundles is ancestral for eupolypods II and a synapomorphy for eupolypods II under deltran optimization. Sori restricted to one side of the vein defines the recently recognized clade comprising Rhachidosoraceae through Aspleniaceae, and sori present on both sides of the vein is a synapomorphy for the Athyriaceae sensu stricto. The results indicate that a chromosome base number of x =41 is synapomorphic for all eupolypods, a clade that includes over two-thirds of extant fern species.

CONCLUSIONS

The integrated approach synthesizes morphological studies with current phylogenetic hypotheses and provides explicit statements of character evolution in the eupolypod II fern families. Strong character support is found for previously recognized clades, whereas few characters support previously unrecognized clades. Sorus position appears to be less complicated than previously hypothesized, and linear sori restricted to one side of the vein support the clade comprising Aspleniaceae, Diplaziopsidaceae, Hemidictyaceae and Rachidosoraceae - a lineage only recently identified. Despite x =41 being a frequent number among extant species, to our knowledge it has not previously been demonstrated as the ancestral state. This is the first synapomorphy proposed for the eupolypod clade, a lineage comprising 67 % of extant fern species. This study provides some of the first hypotheses of character evolution at the family level and above in light of recent phylogenetic results, and promotes further study in an area that remains open for original observation.

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.

Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes

BACKGROUND

Plastid genome structure and content is remarkably conserved in land plants. This widespread conservation has facilitated taxon-rich phylogenetic analyses that have resolved organismal relationships among many land plant groups. However, the relationships among major fern lineages, especially the placement of Equisetales, remain enigmatic.

RESULTS

In order to understand the evolution of plastid genomes and to establish phylogenetic relationships among ferns, we sequenced the plastid genomes from three early diverging species: Equisetum hyemale (Equisetales), Ophioglossum californicum (Ophioglossales), and Psilotum nudum (Psilotales). A comparison of fern plastid genomes showed that some lineages have retained inverted repeat (IR) boundaries originating from the common ancestor of land plants, while other lineages have experienced multiple IR changes including expansions and inversions. Genome content has remained stable throughout ferns, except for a few lineage-specific losses of genes and introns. Notably, the losses of the rps16 gene and the rps12i346 intron are shared among Psilotales, Ophioglossales, and Equisetales, while the gain of a mitochondrial atp1 intron is shared between Marattiales and Polypodiopsida. These genomic structural changes support the placement of Equisetales as sister to Ophioglossales + Psilotales and Marattiales as sister to Polypodiopsida. This result is augmented by some molecular phylogenetic analyses that recover the same relationships, whereas others suggest a relationship between Equisetales and Polypodiopsida.

CONCLUSIONS

Although molecular analyses were inconsistent with respect to the position of Marattiales and Equisetales, several genomic structural changes have for the first time provided a clear placement of these lineages within the ferns. These results further demonstrate the power of using rare genomic structural changes in cases where molecular data fail to provide strong phylogenetic resolution.

Evolution of the rpoB-psbZ region in fern plastid genomes: notable structural rearrangements and highly variable intergenic spacers

BACKGROUND

The rpoB-psbZ (BZ) region of some fern plastid genomes (plastomes) has been noted to go through considerable genomic changes. Unraveling its evolutionary dynamics across all fern lineages will lead to clarify the fundamental process shaping fern plastome structure and organization.

RESULTS

A total of 24 fern BZ sequences were investigated with taxon sampling covering all the extant fern orders. We found that: (i) a tree fern Plagiogyria japonica contained a novel gene order that can be generated from either the ancestral Angiopteris type or the derived Adiantum type via a single inversion; (ii) the trnY-trnE intergenic spacer (IGS) of the filmy fern Vandenboschia radicans was expanded 3-fold due to the tandem 27-bp repeats which showed strong sequence similarity with the anticodon domain of trnY; (iii) the trnY-trnE IGSs of two horsetail ferns Equisetum ramosissimum and E. arvense underwent an unprecedented 5-kb long expansion, more than a quarter of which was consisted of a single type of direct repeats also relevant to the trnY anticodon domain; and (iv) ycf66 has independently lost at least four times in ferns.

CONCLUSIONS

Our results provided fresh insights into the evolutionary process of fern BZ regions. The intermediate BZ gene order was not detected, supporting that the Adiantum type was generated by two inversions occurring in pairs. The occurrence of Vandenboschia 27-bp repeats represents the first evidence of partial tRNA gene duplication in fern plastomes. Repeats potentially forming a stem-loop structure play major roles in the expansion of the trnY-trnE IGS.

Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages

Background

Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats.

Results

We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes.

Conclusions

Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Utility of a large, multigene plastid data set in inferring higher-order relationships in ferns and relatives (monilophytes)

PREMISE OF THE STUDY

The monilophytes (ferns and relatives)-the third largest group of land plants-exhibit a diverse array of vegetative and reproductive morphologies. Investigations into their early ecological and life-history diversification require accurate, well-corroborated phylogenetic estimates. We examined the utility of a large plastid-based data set in inferring backbone relationships for monilophytes. •

METHODS

We recovered 17 plastid genes for exemplar taxa using published and new primers. We compared results from maximum-likelihood and parsimony analyses, assessed the effects of removing rapidly evolving characters, and examined the extent to which our data corroborate or contradict the results of other studies, or resolve current ambiguities. •

KEY RESULTS

Considering multifamily clades, we found bootstrap support comparable to or better than that in published studies that used fewer genes from fewer or more taxa. We firmly establish filmy ferns (Hymenophyllales) as the sister group of all leptosporangiates except Osmundaceae, resolving the second deepest split in leptosporangiate-fern phylogeny. A clade comprising Ophioglossaceae and Psilotaceae is currently accepted as the sister group of other monilophytes, but we recover Equisetum in this position. We also recover marattioid and leptosporangiate ferns as sister groups. Maximum-likelihood rate-class estimates are somewhat skewed when a long-branch lineage (Selaginella) is included, negatively affecting bootstrap support for early branches. •

CONCLUSIONS

Our findings support the utility of this gene set in corroborating relationships found in previous studies, improving support, and resolving uncertainties in monilophyte phylogeny. Despite these advances, our results also underline the need for continued work on resolving the very earliest splits in monilophyte phylogeny.

The horsetail Equisetum arvense mitochondria share two group I introns with the liverwort Marchantia, acquired a novel group II intron but lost intron-encoded ORFs

We studied the genomic structure and RNA editing of mitochondrial cox1, cox2, cob and atp9 from the horsetail Equisetum arvense, a representative of an old fern lineage. Editing of cox1, cob and atp9 mRNAs occur only by C-to-U transitions. No changes were found in cox2 transcripts constituting one of the rare examples of unedited mitochondrial mRNA in land plants. From three intervening sequences in cox1, cox1i395 and cox1i624 are group IB introns homologous to the Marchantia polymorpha cox1 introns, and cox1i747 is a group IIA intron different to other introns found in plant mtDNA. The group II intron cox2i373 is very similar to other introns found in cox2 from vascular plants. While cob and atp9 have no introns and display the gene structure found in seed plants, various nucleotide substitutions abolish the only potential ORF, a LAGLIDADG endonuclease present in cox1i395. Thus, E. arvense mitochondria conserve two group I introns from non-vascular plants, probably inherited from a common ancestor with liverworts. Analogous to seed plants, E. arvense has no potential mitochondrial splicing factors encoded in these introns. This is the first report concerning the presence of vertically inherited group I introns in vascular plant mitochondria.

Phylogeny of marattioid ferns (Marattiaceae): inferring a root in the absence of a closely related outgroup

Closely related outgroups are optimal for rooting phylogenetic trees; however, such ideal outgroups are not always available. A phylogeny of the marattioid ferns (Marattiaceae), an ancient lineage with no close relatives, was reconstructed using nucleotide sequences of multiple chloroplast regions (rps4 + rps4-trnS spacer, trnS-trnG spacer + trnG intron, rbcL, atpB), from 88 collections, selected to cover the broadest possible range of morphologies and geographic distributions within the extant taxa. Because marattioid ferns are phylogenetically isolated from other lineages, and internal branches are relatively short, rooting was problematic. Root placement was strongly affected by long-branch attraction under maximum parsimony and by model choice under maximum likelihood. A multifaceted approach to rooting was employed to isolate the sources of bias and produce a consensus root position. In a statistical comparison of all possible root positions with three different outgroups, most root positions were not significantly less optimal than the maximum likelihood root position, including the consensus root position. This phylogeny has several important taxonomic implications for marattioid ferns: Marattia in the broad sense is paraphyletic; the Hawaiian endemic Marattia douglasii is most closely related to tropical American taxa; and Angiopteris is monophyletic only if Archangiopteris and Macroglossum are included.

Molecular phylogeny of horsetails (Equisetum) including chloroplast atpB sequences

Equisetum is a genus of 15 extant species that are the sole surviving representatives of the class Sphenopsida. The generally accepted taxonomy of Equisetum recognizes two subgenera: Equisetum and Hippochaete. Two recent phylogenetical studies have independently questioned the monophyly of subgenus Equisetum. Here, I use original (atpB) and published (rbcL, trnL-trnF, rps4) sequence data to investigate the phylogeny of the genus. Analyses of atpB sequences give an unusual topology, with E. bogotense branching within Hippochaete. A Bayesian analysis based on all available sequences yields a tree with increased resolution, favoring the sister relationships of E. bogotense with subgenus Hippochaete.

A molecular phylogeny of scaly tree ferns (Cyatheaceae)

Tree ferns recently were identified as the closest sister group to the hyperdiverse clade of ferns, the polypods. Although most of the 600 species of tree ferns are arborescent, the group encompasses a wide range of morphological variability, from diminutive members to the giant scaly tree ferns, Cyatheaceae. This well-known family comprises most of the tree fern diversity (∼500 species) and is widespread in tropical, subtropical, and south temperate regions of the world. Here we investigate the phylogenetic relationships of scaly tree ferns based on DNA sequence data from five plastid regions (rbcL, rbcL-accD IGS, rbcL-atpB IGS, trnG-trnR, and trnL-trnF). A basal dichotomy resolves Sphaeropteris as sister to all other taxa and scale features support these two clades: Sphaeropteris has conform scales, whereas all other taxa have marginate scales. The marginate-scaled clade consists of a basal trichotomy, with the three groups here termed (1) Cyathea (including Cnemidaria, Hymenophyllopsis, Trichipteris), (2) Alsophila sensu stricto, and (3) Gymnosphaera (previously recognized as a section within Alsophila) + A. capensis. Scaly tree ferns display a wide range of indusial structures, and although indusium shape is homoplastic it does contain useful phylogenetic information that supports some of the larger clades recognised.

Reconciling Extreme Branch Length Differences: Decoupling Time and Rate through the Evolutionary History of Filmy Ferns

The rate of molecular evolution is not constant across the Tree of Life. Characterizing rate discrepancies and evaluating the relative roles of time and rate along branches through the past are both critical to a full understanding of evolutionary history. In this study, we explore the interactions of time and rate in filmy ferns (Hymenophyllaceae), a lineage with extreme branch length differences between the two major clades. We test for the presence of significant rate discrepancies within and between these clades, and we separate time and rate across the filmy fern phylogeny to simultaneously yield an evolutionary time scale of filmy fern diversification and reconstructions of ancestral rates of molecular evolution. Our results indicate that the branch length disparity observed between the major lineages of filmy ferns is indeed due to a significant difference in molecular evolutionary rate. The estimation of divergence times reveals that the timing of crown group diversification was not concurrent for the two lineages, and the reconstruction of ancestral rates of molecular evolution points to a substantial rate deceleration in one of the clades. Further analysis suggests that this may be due to a genome-wide deceleration in the rate of nucleotide substitution.

Tree ferns: monophyletic groups and their relationships as revealed by four protein-coding plastid loci

Tree ferns are a well-established clade within leptosporangiate ferns. Most of the 600 species (in seven families and 13 genera) are arborescent, but considerable morphological variability exists, spanning the giant scaly tree ferns (Cyatheaceae), the low, erect plants (Plagiogyriaceae), and the diminutive endemics of the Guayana Highlands (Hymenophyllopsidaceae). In this study, we investigate phylogenetic relationships within tree ferns based on analyses of four protein-coding, plastid loci (atpA, atpB, rbcL, and rps4). Our results reveal four well-supported clades, with genera of Dicksoniaceae (sensu ) interspersed among them: (A) (Loxomataceae, (Culcita, Plagiogyriaceae)), (B) (Calochlaena, (Dicksonia, Lophosoriaceae)), (C) Cibotium, and (D) Cyatheaceae, with Hymenophyllopsidaceae nested within. How these four groups are related to one other, to Thyrsopteris, or to Metaxyaceae is weakly supported. Our results show that Dicksoniaceae and Cyatheaceae, as currently recognised, are not monophyletic and new circumscriptions for these families are needed.

Multiple major increases and decreases in mitochondrial substitution rates in the plant family Geraniaceae

BACKGROUND

Rates of synonymous nucleotide substitutions are, in general, exceptionally low in plant mitochondrial genomes, several times lower than in chloroplast genomes, 10-20 times lower than in plant nuclear genomes, and 50-100 times lower than in many animal mitochondrial genomes. Several cases of moderate variation in mitochondrial substitution rates have been reported in plants, but these mostly involve correlated changes in chloroplast and/or nuclear substitution rates and are therefore thought to reflect whole-organism forces rather than ones impinging directly on the mitochondrial mutation rate. Only a single case of extensive, mitochondrial-specific rate changes has been described, in the angiosperm genus Plantago.

RESULTS

We explored a second potential case of highly accelerated mitochondrial sequence evolution in plants. This case was first suggested by relatively poor hybridization of mitochondrial gene probes to DNA of Pelargonium hortorum (the common geranium). We found that all eight mitochondrial genes sequenced from P. hortorum are exceptionally divergent, whereas chloroplast and nuclear divergence is unexceptional in P. hortorum. Two mitochondrial genes were sequenced from a broad range of taxa of variable relatedness to P. hortorum, and absolute rates of mitochondrial synonymous substitutions were calculated on each branch of a phylogenetic tree of these taxa. We infer one major, approximately 10-fold increase in the mitochondrial synonymous substitution rate at the base of the Pelargonium family Geraniaceae, and a subsequent approximately 10-fold rate increase early in the evolution of Pelargonium. We also infer several moderate to major rate decreases following these initial rate increases, such that the mitochondrial substitution rate has returned to normally low levels in many members of the Geraniaceae. Finally, we find unusually little RNA editing of Geraniaceae mitochondrial genes, suggesting high levels of retroprocessing in their history.

CONCLUSION

The existence of major, mitochondrial-specific changes in rates of synonymous substitutions in the Geraniaceae implies major and reversible underlying changes in the mitochondrial mutation rate in this family. Together with the recent report of a similar pattern of rate heterogeneity in Plantago, these findings indicate that the mitochondrial mutation rate is a more plastic character in plants than previously realized. Many molecular factors could be responsible for these dramatic changes in the mitochondrial mutation rate, including nuclear gene mutations affecting the fidelity and efficacy of mitochondrial DNA replication and/or repair and--consistent with the lack of RNA editing--exceptionally high levels of "mutagenic" retroprocessing. That the mitochondrial mutation rate has returned to normally low levels in many Geraniaceae raises the possibility that, akin to the ephemerality of mutator strains in bacteria, selection favors a low mutation rate in plant mitochondria.