Across the tree of life, molecular phylogenetic studies often reveal surprising relationships between taxa with radically different morphologies that have long obscured their close affiliations. A spectacular botanical example is Rafflesia, a holoparasite that produces the largest flowers in the world, but that evolved from tiny-flowered ancestors within the Euphorbiaceae. Outside of parasitic lineages, such abrupt transformations are rarely seen. One exception involves the “maidenhair ferns” (Adiantum), which are quintessential ferns: beautifully dissected, terrestrial, and shade loving. The closely related “shoestring ferns” (vittarioids), in contrast, have an extremely simplified morphology, are canopy-dwelling epiphytes, and exhibit greatly accelerated rates of molecular evolution. While Adiantum and the vittarioids together have been shown to form a robust monophyletic group (adiantoids), there remain unanswered questions regarding the monophyly of Adiantum and the evolutionary history of the vittarioids. Here we review recent phylogenetic evidence suggesting support for the monophyly of Adiantum, and analyze new plastid data to confirm this result. We find that Adiantum is monophyletic and sister to the vittarioids. With this robust phylogenetic framework established for the broadest relationships in the adiantoid clade, we can now focus on understanding the evolutionary processes associated with the extreme morphological, ecological, and genetic transitions that took place within this lineage.
Early molecular phylogenetic analyses of ferns (Hasebe et al. 1994, 1995) inferred several unexpected associations that had not previously been suspected. Most of these newly recognized relationships, which subsequently drew considerable attention, have stood the test of time. A prominent example is the monophyly of the heterosporous water ferns in the Marsileaceae and Salviniaceae (Hasebe et al. 1994; Rothwell and Stockey 1994; Pryer et al. 1995; Pryer 1999). Another surprise emerged within the Pteridaceae, grouping Adiantum (maidenhair ferns) together with the vittarioids (shoestring ferns) in a well-supported clade now referred to as the adiantoids (Schuettpelz et al. 2007).
Adiantum and the vittarioids could not be more morphologically or ecologically disparate. In coarse morphology, their conspicuous sporophytes look nothing like one another. The leaves of Adiantum are typically broad and finely divided, whereas those of vittarioids are almost always simple and strap-like (Fig. 1A; Tryon and Tryon 1982; Kramer 1990). Fertile Adiantum leaves are uniquely distinguished by their sporangia borne on, and limited to, false indusia, whereas vittarioid sori occur on the laminae (Crane et al. 1995). These groups also display major differences in the morphology of their gametophytes, although these are less obvious to the naked eye (Fig. 1B). The gametophytes of Adiantum, like those of most ferns, are determinate and heart-shaped, with a distinct midrib and broad wings (Nayar and Kaur 1971). They are generally ephemeral to short-lived (months) and are incapable of vegetative reproduction. Vittarioid gametophytes, on the other hand, are indeterminate and ribbon-like. They can be exceptionally long-lived (years) and can also reproduce asexually via propagules called gemmae (Atkinson and Stokey 1964; Farrar 1974, 1985). In addition, Adiantum and vittarioids occupy two dramatically different niches. Whereas the cosmopolitan genus Adiantum usually occurs on shady forest floors, vittarioids generally grow as epiphytes, colonizing tree trunks and canopies of tropical rain forests. The differences between these two groups also extend to their genomes. Although most Adiantum species have diploid chromosome numbers of n = 29 or n = 30 (Löve et al. 1977), nearly all vittarioids studied are n = 60 (Löve et al. 1977), suggesting that at least one genome duplication event occurred early in the evolutionary history of this lineage.
Because of their distinctive simplified morphology, vittarioids have until quite recently been regarded as a distinct family, the Vittariaceae (Tryon and Tryon 1982; Kramer 1990). However, in recent years, studies have not only identified a close relationship with Adiantum, but also suggested that perhaps vittarioids may even be nested within this genus (Prado et al. 2007; Schuettpelz and Pryer 2007). These analyses of plastid data have further revealed yet another dissimilarity between Adiantum and the vittarioids: a striking difference in branch lengths. Vittarioid branches are extraordinarily long relative to Adiantum, or to any other fern lineage, for that matter. As a consequence, branch support across the vittarioid topology is consistently robust, whereas there is mostly weak support among species of Adiantum, especially for the backbone nodes that lie far deeper than the tips (Schuettpelz and Pryer 2007; Schuettpelz et al. 2007).
While it is clear that Adiantum and the vittarioids together compose a robust clade, there are still unanswered questions regarding the monophyly of Adiantum and the evolutionary history of this genus, as well as that of the vittarioids. For example, are there correlates in the morphology, ecology, and life history of vittarioid ferns that may be contributing to their faster rate of molecular evolution? Here we review recent studies and analyze new plastid data, and find strong support for the monophyly of Adiantum and for other deep divergences within adiantoids. This robust phylogenetic framework will permit us, in future studies, to explore the evolutionary processes that resulted in this extreme morphological, ecological, and genomic makeover in ferns.
Materials and Methods
Taxon Sampling—Our taxon sample consists of 16 species of Adiantum, eight vittarioid species, eight species of cheilanthoid ferns (the sister group to the adiantoids), and two outgroup species (Cryptogramma crispa [cryptogrammoid] and Pityrogramma austroamericana [pteridoid]) that were selected based on Schuettpelz et al. (2007) and Rothfels and Schuettpelz (2014); see Appendix 1. While 16 of the approximately 200 species of Adiantum may appear to be sparse coverage, our sampling was carefully informed by our ongoing phylogenetic study focused on the bulk of Adiantum species. As a consequence, our study is the only one to date to capture the deepest divergences within the genus, and includes representatives from all known major clades of Adiantum.
DNA Extraction, Amplification, and Sequencing—DNA was isolated from either silica-dried or herbarium material. Protocols for DNA extraction, amplification, and sequencing followed Schuettpelz and Pryer (2007), Cochran et al. (2014), and Rothfels and Schuettpelz (2014). Sequences were obtained for six plastid loci (atpA, atpB, chlN, rbcL, rpoA, rps4). Primers for atpA, atpB, and rbcL were identical to those used in Schuettpelz et al. (2007) and Rothfels and Schuettpelz (2014); primers used for chlN, rpoA, and rps4 are listed in Table 1.
Sequence Alignment and Data Sets—DNA sequence chromatograms were manually edited and assembled using Sequencher 4.5 (Gene Codes Corporation, Ann Arbor, Michigan). Each plastid region was aligned with AliView (Larsson 2014), which integrates MUSCLE (Edgar, 2004) as the default alignment program. Each alignment was manually inspected and edited. Although alignment was straightforward for the protein-coding loci, there were some indels in the non-protein-coding regions that rendered the alignment ambiguous; these were excluded prior to subsequent analyses. Unsequenced portions of plastid regions were coded as missing data. Six individual data sets were compiled, one for each of the plastid regions. Data set alignments and phylogenetic trees are deposited in the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.4m6s6. Eighty-one newly obtained DNA sequences were deposited in GenBank (Appendix 1).
Phylogenetic Analyses—Separate phylogenetic analyses were conducted for each of the six plastid data sets using maximum likelihood (ML) in PAUP* 4.0a136 (Swofford 2002). We first inferred a maximum parsimony tree and used the AutoModel function in PAUP* to perform model selection under the AICc. With the best-fit model selected, and the model parameters estimated for each data set, ML analyses were conducted in PAUP* with tree bisection and reconnection branch swapping and 100 random-addition-sequence replicates. Maximum likelihood trees for each individual plastid region were visually inspected for conflicts supported by bootstrap values ≥ 70% (Mason-Gamer and Kellogg 1996). Because no instances of mutually well-supported incongruence were detected when comparing phylogenies across different plastid regions, the six separate data sets were concatenated and analyzed together using ML and Bayesian inference (BI).
For the concatenated data set, we used PartitionFinder (Lanfear et al. 2012) to determine the optimal data-partitioning scheme and substitution models according to the AICc (Table 2). ML tree searches and ML bootstrap (MLBS) analyses (1,000 replicates) were carried out from eight independent random-addition-starting trees in Garli 2.0 (Zwickl 2006) with “genthreshfortopoterm” set to 1,000,000 and 100,000, respectively. MrBayes 3.2 (Ronquist et al. 2012) was used to conduct BI analyses. Because it is not possible to implement some of the best-fitting models in MrBayes, another PartitionFinder analysis was run to choose more applicable models. Two independent Markov chain Monte Carlo (MCMC) runs were carried out, each with four chains (one heated and three cold) running for 20 million generations. Priors followed the default settings with a flat Dirichlet distribution for both the stationary state frequencies and the substitution rates, and trees were sampled every 1,000 generations. The substitution parameters were unlinked, and the rate prior was set to allow variation among the subsets. After the MCMC runs, the output parameters were inspected in Tracer v1.5 (Rambaut and Drummond 2009) to ensure convergence and proper mixing. The first 25% of the sample was discarded as burn-in and the remainder was used to calculate a 50% majority-rule consensus tree.
Seven aligned data matrices (six single-region matrices and one concatenated) were analyzed for this study; a summary of sequence characteristics, best-fit models of sequence evolution, and tree statistics appears in Table 3. Trees resulting from the maximum likelihood and Bayesian analyses were identical in topology. The best ML tree (InL = -38,392.964) from the analysis of the concatenated six-plastid loci (atpA, atpB, chlN, rbcL, rpoA, and rps4) is shown in Fig. 2. Nearly all of its internal nodes (24 out of 31) are highly supported, with ML bootstrap support ≥ 70% and Bayesian posterior probability (PP) support ≥ 0.99 (Fig. 2). The monophyly of both Adiantum (82% MLBS, 1.0 PP) and the vittarioid clade (100% MLBS, 1.0 PP) is robustly supported, as is their sister relationship to one another (100% MLBS, 1.0 PP; Fig. 2). There is strong support for all relationships across vittarioids, whereas within Adiantum, the entire backbone of the clade is weakly supported (<70% MLBS, <0.97 PP).
Phylogenetic analyses from across the tree of life have revealed that the rate of molecular evolution in closely related lineages can be very similar, or can vary dramatically (Lanfear et al. 2010). Within ferns, notable molecular rate heterogeneity has been reported from horsetails (Des Marais et al. 2003) and filmy ferns (Schuettpelz and Pryer 2006) to vittarioids (Rothfels and Schuettpelz 2014) and polygrammoid ferns (Schneider et al. 2004). Sometimes, but not always, this rate variation appears to correlate strongly with certain other aspects of biology. For example, Soltis et al. (2002) and Korall et al. (2010) found that an abrupt rate deceleration coincided with the evolution of the long-lived, tree-like habit at the base of the tree fern clade (i.e. tree ferns, with longer generation times, consistently have slower rates of molecular evolution).
Primers used for DNA amplification and sequencing of plastid chlN, rpoA, rps4 in this study.
Optimal data-partitioning scheme and substitution models for the concatenated data set of six plastid loci; determined according to the AICc using PartitionFinder (Lanfear et al. 2012).
One of the challenges posed by molecular evolutionary rate heterogeneity has been the associated difficulty of recovering phylogenetic topologies that reflect accurate relationships (Schuettpelz and Pryer 2006; Rothfels et al. 2012). Although all molecular analyses to date bring Adiantum together with the vittarioids in a strongly supported clade (adiantoids), they have also repeatedly struggled to find support for the monophyly of Adiantum with respect to the morphologically and ecologically highly dissimilar vittarioids. Three conflicting topologies have been recovered for relationships within adiantoids. The first places Adiantum sister to vittarioids, but with weak support for the monophyly of Adiantum (Fig. 3A; Schuettpelz et al. 2007, their Fig. 1): 67 taxa (16 adiantoids), atpA/atpB/rbcL; Lu et al. 2012 (their Fig. 2): 98 taxa (74 adiantoids), atpA/atpB/rbcL). A second topology suggests a paraphyletic Adiantum, with A. raddianum sister to a weakly supported clade of vittarioids + the rest of Adiantum (Fig. 3B; Schuettpelz et al. 2007, their Fig. 3): 147 taxa (36 adiantoids), rbcL). The third topological option is again for a paraphyletic Adiantum, but this time with A. raddianum sister to vittarioids, and the rest of Adiantum sister to that clade (Fig. 3C; Schuettpelz and Pryer 2007, their Fig. 1B): 400 taxa (15 adiantoids), atpA/atpB/rbcL). The last two topologies both suggest that the vittarioids may actually be derived from within Adiantum (Fig. 3B, C), a hypothesis that is difficult to reconcile with the extraordinary degree of morphological conservatism within Adiantum compared to other large fern genera.
The most significant phylogenetic result from our study of six-plastid loci (atpA/atpB/chlN/rbcL/rpoA/rps4) from 34 taxa (24 adiantoids) is strong support for a monophyletic Adiantum, with 82% ML bootstrap and 1.0 posterior probability support (Figs. 2, 3D). The only other study to date to convincingly demonstrate the monophyly of Adiantum was a six-locus (from across all three genome compartments) data set for 26 taxa (16 adiantoids, Rothfels and Schuettpelz 2014, their Fig. 1C: atpA/atpB/rbcL [plastid] + gapCp [nuc] + atp1/ nad5 [mt]). Also noteworthy—although only three Adiantum species (but including A. raddianum) and two vittarioids out of 73 ferns were included in their study—the fern phylogeny inferred by Rothfels et al. (2015) from 25 low-copy nuclear genes robustly refutes the hypothesis that vittarioids are nested within Adiantum.
In their combined three-plastid marker (atpA/atpB/rbcL) analysis of 98 taxa (74 adiantoids), Lu et al. (2012, see their Fig. 2) found strong maximum parsimony bootstrap support (87%) for the monophyly of Adiantum, but significant Bayesian support at that node was lacking (PP = 0.81). While bootstrap support and posterior probabilities measure different things and are thus not directly comparable, posterior probability support values almost universally exceed support values from parsimony and likelihood bootstrap support (except, of course, when they are equal; Hillis and Bull 1993; Alfaro et al. 2003; Rothfels et al. 2012). The fact that the posterior probability for this particular node is so low in Lu et al. (2012), especially given that this is a situation (long branches subtended by very short internodes) that is particularly prone to long-branch attraction issues (and thus unreliable MP inference), undermines any confidence in their result. Put simply, the two methods (MP and BI) have contrasting support for a node where MP might reasonably be expected to fail. The node in question is in fact the only node with a Bayesian posterior probability lower than its maximum parsimony bootstrap in the entire Lu et al. (2012) study. Thus, we were not convinced that their study demonstrated the monophyly of Adiantum, which prompted our study using an expanded plastid sequence data set.
Sequence characteristics for six plastid loci, best-fit sequence evolution models, and resulting tree statistics. Missing data does not include indels.
We believe that our increased sampling of plastid markers (including the two new loci chlN and rpoA) together with model-based approaches, such as maximum likelihood and Bayesian analyses (Swofford et al. 1996; Ronquist and Huelsenbeck 2003), as well as informed taxon sampling from our parallel large-scale analysis of Adiantum, allowed us to succeed in achieving credible support, comparable across both measures, for the monophyly of Adiantum. This level of confidence is requisite to proceeding further with our ongoing within-clade studies; we are thus now moving forward with an appropriate sampling of vittarioids as outgroup for a nearly complete phylogenetic study based on four plastid markers (rbcL, atpA, chlN, rpoA) focused on the bulk of Adiantum species (Huiet et al. unpubl.), and a comparable study for the vittarioid clade (Schuettpelz et al. unpubl.).
Most other relationships within the adiantoids (Fig. 2) are in agreement with earlier studies (Schuettpelz et al. 2007; Lu et al. 2012), including strong support for a Rheopteris + Vaginularia clade (Ruhfel et al. 2008; Rothfels and Schuettpelz 2014). Broad relationships within the cheilanthoid ferns also mirror those found in recent studies (Prado et al. 2007; Schuettpelz et al. 2007; Rothfels et al. 2008; Windham et al. 2009; Eiserhardt et al. 2011); however, strong support for Calciphilopteris ludens as a member of this clade, which was recovered by Rothfels and Schuettpelz (2014) using data from all three genomic compartments, was not achieved here. To date, analyses using only plastid data typically resolve C. ludens (without strong support) as either sister to the rest of the cheilanthoids (Schuettpelz et al. 2007) or the adiantoids (this study).
An extreme makeover occurred during the evolutionary history of the vittarioid ferns, involving morphological transformations in both the sporophyte and gametophyte phases of their life cycle, the evolution of epiphytism, genome duplication(s), and molecular rate acceleration. With a robust phylogenetic framework now in place, we aim to identify, in future studies, the underlying causal mechanisms that may have contributed to this rather spectacular transformation.
The authors thank the herbarium curators and staff of A, COLO, DUKE, GOET, P, TUR, UC, and UTC for access to specimens needed for this study. This research was funded by National Science Foundation grants DEB-1145614 to K. M. P. and L. H. and DEB-1405181 to E. S. as well as NSF DDIG awards DEB-1110652 to K. M. P. and C. J. R. and DEB-1407158 to K. M. P. and F.-W. L. Assistance from T.-T. Kao is gratefully acknowledged, as are the thoughtful comments of M. D. Windham.
Appendix 1. Vouchers and GenBank accession numbers for taxa used in our molecular phylogenetic analysis. Taxon, collection locality; voucher specimen collector and no. (herbarium acronym), Fern DNA database number (fernlab.biology.duke.edu), GenBank accession, citation for previously published data for atpA; atpB; chlN; rbcL; rpoA; and rps4 (in that order). A dash (—) indicates not applicable; asterisk (*) indicates only 71 bp of sequence data were obtained for the chlN gene of Adiantum hispidulum, which is too short for GenBank to accession, therefore, the sequence is provided here: CGGCAACTACATTAATGCGTCGAAGGAAATGCCAGTTAGTTGGAGCACCTTTCCCAATTGGTCCAGATGGG