Am. J. Bot.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


(American Journal of Botany. 2008;95:626-641.)
doi: 10.3732/ajb.2007308
© 2008 Botanical Society of America, Inc.
  Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter
What's this?
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Murdock, A. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Murdock, A. G.
Agricola
Right arrow Articles by Murdock, A. G.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

Systematics and Phytogeography

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

Andrew G. Murdock2

Department of Integrative Biology, University of California, Berkeley, 1001 Valley Life Sciences Bldg., California 94720-2465 USA

Received for publication 1 January 2007. Accepted for publication 14 February 2008.

ABSTRACT

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.

Key Words: AngiopterisChristenseniaDanaea • long-branch attraction • Marattia • Marattiaceae • rooting

Recent phylogenetic studies have greatly improved our understanding of the evolutionary relationships of ferns. In particular, relationships of lineages within the diverse leptosporangiate ferns have been significantly clarified (e.g., Pryer et al., 2001Go, 2004Go; Schuettpelz et al., 2006Go; Schuettpelz and Pryer, 2007Go). The leptosporangiate ferns, Marattiales (marattioids), Ophioglossales (ophioglossoids), Psilotales (whisk ferns), Equisetales (horsetails), and Spermatophytes (seed plants) (classification following Smith et al., 2006Go) are all well-supported clades in phylogenetic studies; however the interrelationships among these lineages, that together constitute the Euphyllophytes, remain unresolved and controversial (Rothwell and Nixon, 2006Go; Schuettpelz et al., 2006Go). DNA sequence data from all compartments of the genome have generated an array of conflicting, but generally poorly supported topologies for the early branching events in the Euphyllophytes (Nickrent et al., 2000Go; Pryer et al., 2001Go, 2004Go; Wikström and Pryer, 2005Go; Qiu et al., 2006Go) and genomic structural characters have also proven inconclusive (Kelch et al., 2004Go; Wolf et al., 2005Go). Previous studies that have attempted to integrate fossils with molecular data in phylogenetic analyses have found that inclusion of fossils can alter tree topology, resolution, and support depending on the approach, in some cases causing basal relationships in Euphyllophytes to be more poorly or differently resolved than in analyses including only extant taxa (Rothwell and Nixon, 2006Go; Schneider, 2007Go).

The lack of a clear phylogenetic consensus at the base of the Euphyllophyte tree can be attributed to the loss of historical signal through time and through extinction-limited sampling, methodological difficulties arising from long branches separating the extant crown groups, and relatively high levels of evolutionary rate heterogeneity (Soltis et al. 2002Go), problems common to most deep, unresolved nodes of the tree of life (Mishler, 2000Go; Giribet, 2002Go; Burleigh and Mathews, 2004Go; Embley and Martin, 2006Go). All the extant members of the early-diverging Euphyllophyte lineages have been sampled sufficiently in previous studies to be certain that the long-branch problem cannot be mitigated by further sampling, with one possible exception: Christensenia, a rare and morphologically unique marattioid fern genus that was resolved as sister to all other marattioid ferns in a morphological cladistic analysis (Hill and Camus, 1986Go), has not been included in previous molecular analyses.

Of the extant early-branching lineages of the Euphyllophytes, the marattioid ferns (also "marattialean" or "marattiaceous" in some earlier literature) are perhaps the least studied. Marattioid ferns are an ancient lineage of eusporangiate ferns with no close relatives whose extant taxa are essentially restricted to the tropical regions of the world. Historically, marattioids were thought to be closely related to ophioglossoids based on a small number of comparable morphological characters (Campbell, 1911Go), but recent results have found marattioids to be phylogenetically isolated from other fern lineages (Smith, 1995Go), with most molecular analyses resolving Ophioglossales as sister to Psilotales and several suggesting a possible affinity between Marattiales and Equisetales (Pryer et al., 2001Go, 2004Go; Wikström and Pryer, 2005Go; Qiu et al., 2006Go). Despite the uncertainty concerning the interrelationships of early-branching fern lineages, the hypothesis that marattioid and leptosporangiate ferns are closely related lineages (if not sister lineages) has not been questioned.

Marattioids have one of the most extensive fossil records of any modern fern lineage, with numerous fossils of extinct lineages dating from the Carboniferous, scattered Mesozoic fossils that more closely resemble modern genera, but are essentially lacking from the Tertiary fossil record (Millay, 1979Go; Hill and Camus, 1986Go; Taylor and Taylor, 1993Go; Collinson, 2001Go). Extant marattioids can be diagnosed morphologically on the basis of multiple characters, some of which are apparent in putative fossil relatives: large, thick-walled sporangia that are fully or partially fused into synangia, with each sporangium producing large quantities of small spores; complex polycyclic stem vasculature; leaves basally flanked by paired, starchy, vascularized auricles (often called stipules); leaves possessing prominent fleshy swellings (pulvini) at nodes and/or bases of segments; extensive mucilage canal systems throughout the tissues of the plant; multicellular root hairs; cyclocytic multicyclic stomata (Campbell, 1911Go; Hill and Camus, 1986Go; Rolleri et al., 2003Go).

Taxonomy
Modern marattioids have most commonly been treated in a single family Marattiaceae Kaulf. (e.g., Rolleri et al., 2003Go; Smith et al., 2006Go), although some authors (e.g., Pichi Sermolli, 1977Go) have further dissected this into smaller families, Danaeaceae C. Agardh, Angiopteridaceae Fée ex J. Bommer, and Christenseniaceae Ching (= Kaulfussiaceae Campb.) with varying circumscriptions (Murdock et al., 2006bGo). Marattiaceae s.l. is typically treated as having four genera: Danaea Sm., restricted to the neotropics; Christensenia Maxon, restricted to Southeast Asia; Marattia Sw., with a pantropical distribution; and Angiopteris Hoffm., native to Madagascar, Asia, and Oceania (introduced in Hawaii and the American tropics). Some authors have recognized up to four poorly defined genera segregated from Angiopteris s.l.: Archangiopteris Christ and Giesenh., Macroglossum Copel., Protangiopteris Hayata, and Protomarattia Hayata, all with highly restricted ranges in Asia. Presl (1845)Go subdivided Marattia into five genera (Marattia,Discostegia C. Presl, Gymnotheca C. Presl, Stibasia C. Presl, and Eupodium J. Sm.), a taxonomy also followed by De Vriese and Harting (1853)Go in their monograph of the family; however, more recent treatments (e.g., Copeland, 1947Go; Rolleri et al., 2003Go) have treated these genera within a more broadly circumscribed Marattia (referred to in this paper as Marattia s.l.).

While marattioids are united as a group by multiple distinctive characters, they are at the same time morphologically diverse. Danaea is the only reproductively dimorphic genus in the family and is also distinguished from other genera by long synangia that are sunken into a thickened lamina of the fertile fronds, sporangia that dehisce via apical pores, and large peltate scales. Christensenia is unique among marattioids in having palmately lobed to pedately compound blades, reticulate venation, circular synangia that appear more or less scattered on the abaxial side of the frond. Marattia is the most widespread and morphologically diverse genus, the paleotropical species having sessile synangia, most neotropical species having sessile or slightly raised synangia with lower pinnae more divided, and the neotropical M. laevis Sm. complex having distinctive stalked synangia. Angiopteris is distinguished by having sporangia that are fused only at the base and by having monolete spores; all other genera have trilete spores.

While there are some taxonomic difficulties with nearly all lineages, marattioids present an especially challenging combination of problems for a taxonomist, analogous in many ways to tree ferns, and these problems have contributed to a chronically unsettled species-level taxonomy. Many marattioids grow to be extremely large with some species having fronds over 6 m long. The large size of marattioids presents problems for taxonomists due to the necessarily fragmentary nature of most collections, some of which, including type specimens, consist of only a few small pinnules. Apart from size, separating plasticity from heritable differences is perhaps the greatest challenge for systematists working on marattioid ferns. Within populations and even within individuals one can find striking differences in morphology. Light, temperature, and frond age have a strong effect on all aspects of frond morphology, and irregular frond division is common in the family (Backer 1929Go; Asama, 1960Go; Holttum, 1978Go; A. Murdock, personal observation); herbarium specimens commonly present only a small picture of the overall variability of any given plant. Some investigators have misinterpreted hypervariability as taxonomically informative and have named numerous poorly differentiated species, making this group nomenclaturally burdensome. Additionally, the size, slow growth rate, and the difficulty of obtaining and propagating live plants make common garden experiments and even moderately large living collections impractical.

Because of the taxonomic difficulties in marattioid ferns, species diversity estimates have varied widely among investigators. In Angiopteris, Ching (1959)Go recognized 62 species and an additional nine species of Archangiopteris just for China, while Rolleri (2003)Go recognized only 10 species worldwide; others have only recognized one polymorphic species in Angiopteris (e.g., Hooker and Baker, 1868Go; Backer, 1929Go). Nearly all species of Angiopteris have at one time been synonymized with A. evecta (type from Tahiti), and nearly all Old World Marattia species have at one time been synonymized with M. fraxinea (type from Réunion). One result of this taxonomic uncertainty is that assessment of conservation status of marattioid taxa is problematic. Many species have been construed as highly restricted endemics (Fu and Jin, 1992Go; Hsu et al., 2000Go; Kingston et al., 2004Go), and some are threatened by deforestation, plant invasion, and predation by feral pigs (Warshauer, 1979Go; Pickard, 1983Go; Brownsey and Smith-Dodsworth, 2000Go). These factors make a strong argument in favor of using molecular data and the resulting phylogenetic inferences to make informed taxonomic and conservation decisions for marattioid ferns.

Previous work
A morphological analysis by Hill and Camus (1986Go) was the first study to examine the evolutionary relationships of marattioid ferns in a cladistic framework and includes one of the most detailed discussions of morphological characters in the group to date. In addition to extant taxa, Hill and Camus (1986)Go included several well-characterized fossil taxa; Liu et al. (2000)Go used a similar approach, but they were more concerned with interrelationships of fossil taxa. Hill and Camus (1986)Go, using Ophioglossales as an outgroup, resolved Christensenia as sister to a clade containing all other marattioids and resolved Angiopteris as monophyletic only if defined in a broad sense to include Archangiopteris,Macroglossum, and Protomarattia. Both Danaea and Marattia were resolved as monophyletic; however, relationships among Danaea,Marattia, and Angiopteris s.l. remained unresolved. However, Hill and Camus (1986Go, p. 276) stated that Marattia,Macroglossum,Angiopteris, and Archangiopteris were all found to be paraphyletic in their analysis, despite the fact that their figures show only Angiopteris to be paraphyletic, and then only if narrowly circumscribed. Additionally, they implied that Marattia laevis should be segregated as a more narrowly defined genus based on cladistic principles and then discussed multiple alternatives to recognizing a paraphyletic Marattia; however, these conclusions and discussions are not consistent with their trees, where Marattia s.l. is shown as monophyletic.

Apart from recent phylogenetic work on Danaea (Christenhusz et al., 2008Go) and the complete plastid genome sequence of Angiopteris evecta (Roper et al., 2007), relatively few marattioid ferns have been sequenced. Species that have been sequenced have mostly been part of large-scale phylogenetic studies (e.g., Hasebe et al., 1993Go; Manhart, 1994Go, Manhart, 1995Go; Pryer et al., 2001aGo, 2004Go; Wikström and Pryer, 2005Go; Qiu et al., 2006Go) with minimal sampling of marattioids. There have also been several recent papers focused specifically on Chinese and Taiwanese species of Angiopteris or Archangiopteris (Hsu et al., 2000Go; Chiang et al., 2002Go; Li and Lu, 2007Go; S. Lin, Shenzhen Fairy Lake Botanical Gardens, unpublished manuscript), which have uncovered very little genetic differentiation between the many named species in the region. No molecular analysis has included the genus Christensenia, apart from Christenhusz et al. (2008Go), where it was used as an outgroup to Danaea. Inclusion of Christensenia in nucleotide sequence data sets could potentially improve the resolution of the placement of marattioid ferns in relation to other early-diverging fern lineages.

For this study, a phylogeny of marattioid ferns was inferred using nucleotide sequence data and morphology. This phylogeny helps to clarify the poorly understood taxonomy of marattioids and to interpret the morphological evolution and biogeographic history of the group. The understanding of polarity of morphological evolution in marattioid ferns is crucial for understanding relationships of extant taxa to putative fossil relatives, dating divergence times, and assessing morphological homologies with related lineages. Because confidence in rooting is essential for determining polarity and rooting can be difficult in groups that are phylogenetically isolated from other lineages, rooting hypotheses and outgroup relationships were explored using a multifaceted approach designed to isolate potential sources of bias and produce a consensus root position.

MATERIALS AND METHODS

Sampling
Sampling for this study represented the broad range of morphology and geographic distribution in marattioids, including multiple individuals from all currently recognized genera (Christensenia, Danaea, Angiopteris, Archangiopteris, Macroglossum, and Marattia) as well as individuals from all previously recognized genera and subgenera within Marattia s.l. (i.e., Discostegia C. Presl, Eupodium J. Sm., Gymnotheca C. Presl, Stibasia C. Presl, and subgenus Mesosorus Rosenst.). Generic names are used in the sense of Copeland (1947)Go. Species names used are based on regional floristic treatments, taxonomic revisions, and the taxonomic opinions of the author based on examination of types and numerous other collections. Names shown in quotes are those that are used in regional floras but are likely to change in upcoming revisions (e.g., M. Christenhusz, University of Turku, unpublished manuscript, for Danaea).

Morphology
The previously published morphological analysis by Hill and Camus (1986)Go provided the starting point for building a new morphological matrix for marattioid ferns. Many characters used by Hill and Camus were at least partially nonindependent, and character states were questionably discrete in some cases, so the data set was critically reviewed and extensively recoded. Many characters useful for taxonomic purposes were excluded due to a lack of independence or a nondiscrete nature. While there are several cladistic analyses of fossils that included some marattioid ferns (e.g., Rothwell, 1999Go; Rothwell and Nixon, 2006Go), these contained few additional characters that were useful for resolution within marattioid ferns. Additional characters have been added from a variety of literature sources and from observation of specimens and living material.

Phylogenetic markers
Multiple regions of the chloroplast genome, including rps4 plus the rps4–trnSGGA intergenic spacer (rps4–trnS) (~900 bp), and the region from trnSGCU to trnGUUC (including the trnG intron) (trnSGG) (~1600–2100 bp) were sequenced to provide resolution within the ingroup. Eighty-eight accessions were sequenced for rps4–trnS and trnSGG for maximal resolution within marattioid ferns. The chloroplast region rps4–trnS was first used in fern phylogenetics by Smith and Cranfill (2002)Go and has been widely used in other groups of green plants. Apart from Small et al. (2005)Go, a study that assessed the potential utility of numerous chloroplast noncoding regions for fern phylogenetics, trnSGG has not been widely used in fern phylogenetic studies. It appears, however, to be a highly promising phylogenetic marker for lower-level phylogeny reconstruction in many groups of plants (Hamilton, 1999Go; Shaw et al., 2005Go, 2007Go; Szövényi et al., 2006Go). In most nonflowering plants, trnSGG includes multiple short coding regions (trnSGCU, psaM, ycf12, and trnGUUC), which serve as positions for primer design and assist with sequence alignment. Three other chloroplast regions with relatively slow sequence evolution, atpB (~1270 bp), rbcL (~1370 bp), and rps4 (~600 bp) were sequenced to provide resolution at deeper levels and assist with outgroup rooting. Because rbcL, atpB, and rps4 had limited variability in the ingroup, a smaller data set of 36 taxa (with 22 ingroup taxa selected to represent all major clades, and 14 outgroup taxa) was assembled for outgroup rooting.

Primers
Primers used in this study are shown in Table 1, with corresponding primer locations for trnSGG shown in Fig. 1. Novel primers were designed when necessary after initial sequencing efforts using published primer sequences and comparison with sequences in GenBank. The forward PCR primer for trnSGG from Shaw et al. (2005Go) is not optimal for pteridophytes because of base mismatches, but the reverse primer is generally applicable. Nor are the updated primers in Shaw et al. (2007)Go, designed for angiosperms, not for pteridophytes because of base mismatches in both the forward and reverse primers and the location of the forward primer over an insertion found in the trnS gene in Equisetum. For trnSGG, forward and reverse PCR primers trnS-GCU-F1 and 3'-trnG-UUC, and internal sequencing primers 5'-trnG2G and 5'-trnG2S can be used for most nonflowering plants; additional internal primers were designed specifically for use in Marattiaceae and may work in related fern families; with minor modifications, these could be used in more distantly related ferns and other nonflowering plant lineages. Locations of primer sites in trnSGG are shown in Fig. 1. Internal primers for rbcL were designed specifically for this study and will likely have limited applicability in other lineages; other primers used for rbcL and atpB are variations on those used in other fern studies (e.g., Hasebe et al., 1993Go; Wolf, 1997Go; Pryer et al., 2001a, bGo) but have been redesigned for use in Marattiaceae and other early-diverging fern lineages.


View this table:
[in this window]
[in a new window]

 
Table 1. Primers used in this study. PCR amplification primers are in boldface; direction indicated by F (forward) or R (reverse). Main primers used for sequencing trnSGG indicated by an asterisk (*). Primer location for trnSGG, as shown in Fig. 1, is in boldface followed in brackets by the coordinate within trnSGG of the 5'primer end using the Angiopteris cp genome (Roper et al., 2007) as the standard.

 

Figure 1
View larger version (4K):
[in this window]
[in a new window]

 
Fig. 1. Map of the chloroplast trnSGG region and location of primers used in this study. Primer locations and direction correspond with Table 1; black circles indicate PCR primers; white circles indicate internal sequencing primers. Location of a gene above or below the line indicates relative direction of transcription.

 
Sequencing
Extraction of genomic DNA was performed using the Qiagen DNeasy Plant Mini Kit (Valencia, California, USA) from fresh, silica-dried, or recent herbarium specimen leaf material. In general, fresh or silica-dried material was needed; amplification of PCR products from herbarium specimens was often problematic. PCR was performed on 1:1 to 1:50 dilutions of total genomic extraction, using Bioneer PCR Premix Tubes (Bioneer, Alameda, California, USA), Amersham PCR Beads, or AmpliTaq Gold (Applied Biosystems, Foster City, California, USA) (specific reaction mixes available on request from the author). PCR was performed on an MJ Research (Waltham, Massachusetts, USA) PTC-200 thermal cycler with protocols following standard PCR protocols for atpB and rbcL, the protocol of Smith and Cranfill (2002)Go for rps4–trnS, and the protocol of Shaw et al. (2005)Go for trnSGG. Touchdown PCR (Don et al., 1991Go) was used to amplify several samples that were problematic under standard protocols. PCR product was visualized and quantified on agarose gels stained with ethidium bromide and subsequently purified using USB ExoSAP-IT or USB PCR Pre-Sequencing Kit (USB Corp., Cleveland, Ohio, USA) prior to sequencing. Nucleotide sequencing was done on Applied Biosystems 96 capillary 3730xl DNA Analyzers or an Applied Biosystems 377 sequencer.

Sequence lignment and phylogenetic analysis
Contigs of sequences were constructed and proofread using the program Sequencher (Gene Codes, Ann Arbor, Michigan, USA). The rooting data sets were supplemented with sequences from GenBank (see Appendix 1); outgroup taxa were mostly single species, but in some cases were a composite of sequences from two related species when necessary due to availability of sequences and plant material. Coding regions were confidently alignable by eye. Where indels occurred in coding regions, genes were translated into amino acids using the program MacClade version 4.08 (Maddison and Maddison, 2005Go) to accurately identify the triplets involved in the indels. Noncoding regions were aligned first by the program ClustalW (http://www.ebi.ac.uk/clustalw/) and then fixed by eye where obvious alignment errors occurred. Areas of uncertain alignment and complex nested indels were excluded in final analyses. Microstructural characters were coded as present/absent where discrete, following Graham et al. (2000Go). Identical taxa (i.e., with a pairwise difference of 0 in informative characters), of which there were several, particularly in Angiopteris, were merged prior to analysis when necessary for longer analyses.

Data from each gene region were analyzed separately and combined into three data sets: {1} rps4–trnS+trnSGG; {2} rps4–trnS+trnSGG with morphology and coded indels; and {3} atpB+rbcL+rps4 (for rooting analysis)—these data sets will be referred to later by their respective numbers inside the braces.

Maximum parsimony (MP), maximum likelihood (ML), and Bayesian inference using a Markov chain Monte Carlo approach (BI) were used to analyze phylogenetic data in this study using PAUP* version 4.0b10 (Swofford, 2002Go) and MrBayes version 3.1 (Huelsenbeck and Ronquist, 2001Go; Ronquist and Huelsenbeck, 2003Go). The programs ModelTest version 3.7 (Posada and Crandall, 1998Go) and MrModelTest version 2.2 (Nylander, 2004Go) were used to select models for ML and BI analyses, respectively, with model selection based on the Akaike information criterion (AIC). For MP and ML analyses, all data sets were analyzed using a heuristic search with 10 random sequence addition replicates and tree-bisection-reconnection (TBR) branch swapping. Starting trees for full ingroup MP analyses were obtained using 100000 random addition sequence replicates with no swapping and MulTrees unselected; the best trees found were then used as starting trees for complete TBR branch swapping with MulTrees selected. MP bootstrap analyses were run using heuristic searches with 10 random addition sequence replicates or were limited to 10 million rearrangements per replicate when necessary to complete an analysis. Decay analyses used commands generated in MacClade, which were then executed in PAUP*. BI analyses were continued until apparent stationarity, diagnosed using the program Tracer version 1.3 (Rambaut and Drummond, 2003Go). The BI analysis of all data sets was based on 3000000 generations, four parallel chains, sampled every 1000 generations, with a burn-in of 500000 generations excluded from the final analysis. Tests for the presence of a molecular clock in the ingroup were performed using a likelihood ratio test (Felsenstein, 1981Go) for all genes separately and combined in PAUP*.

Rooting
One hundred random outgroup sequences used in rooting trials were generated in MacClade using the average base composition of the ingroup (which did not significantly differ from the outgroups). Each of the 100 random sequences was used individually as an outgroup to the rooting data set {3} and analyzed under MP, and the first 20 sequences were analyzed under ML, to assess tendency of random outgroups to attach to particular ingroup branches. In both MP and ML random rooting trials, ingroup topology was constrained to match that of Fig. 2 using a backbone constraint tree, and random sequence addition was used. Shimodaira–Hasegawa tests (Shimodaira and Hasegawa, 1999Go; Goldman et al., 2000Go) comparing the likelihood scores, and Wilcoxon signed-rank tests (Templeton, 1983Go) comparing the parsimony tree length scores of 41 constraint trees corresponding to all possible root attachment points within the atpB+rbcL+rps4 data set were performed with three possible outgroups (Equisetaceae, Ophioglossaceae, Osmundaceae) following the general procedures of Graham et al. (2002Go) and Buckley et al. (2001)Go in PAUP* using RELL estimation. Constraint trees for all root positions analyzed were generated using MacClade. Lundberg root vectors were created by reconstructing the nodal states for a polytomy containing all outgroup taxa in MacClade with equivocal reconstructions coded as unknown. Midpoint and Lundberg rootings were calculated using PAUP*.


Figure 2
View larger version (28K):
[in this window]
[in a new window]

 
Fig. 2. Phylogeny of Marattiaceae (ingroup only) based on maxiumum likelihood (ML) analysis of combined trnSGG and rps4–trnS{1}. Support values: Bayesian posterior probabilities / maxiumum parsimony (MP) bootstrap / decay index. "+" = 100; "-" = <50 (0 for decay index values). The branch indicated by *, supported only by MP analysis, collapses in the ML tree but is preserved here for explanatory purposes. Genera recognized in a taxonomic revision by Murdock (2008) indicated by branches with support values in ovals. Hypothesized root placement indicated by a dashed line (see Discussion, Rooting).

 
RESULTS

Phylogenetic relationships
Sequences of trnSGG proved to be the most informative within the 88 taxon ingroup data set, with 483 variable characters (367 parsimony informative) of an aligned length of 1992 characters; rps4–trnS had 246 variable characters (213 parsimony informative) of an aligned length of 992 characters. Both regions contained microstructural characters (indels and repeat motifs): 56 in trnSGG and 14 in rsp4–trnS (with only one of these in the coding region of rps4). In the 36-taxon rooting data set {3} (22 ingroup taxa and 14 outgroup taxa), rbcL had 566 variable characters (432 parsimony informative) of an aligned length of 1371 characters; atpB had 585 variable characters (426 parsimony informative) of an aligned length of 1269 characters; rps4 had 331 variable characters (266 parsimony informative) of an aligned length of 523 characters.

Support for the monophyly of Marattiaceae was consistently high in all analyses that included outgroups (100% Bayesian posterior probability and MP bootstrap), regardless of root placement within the family (described later). Support for internal branches was variable: in general, clades of closely related (and morphologically similar taxa) were highly supported on long branches with 100% Bayesian posterior probability, 100% MP bootstrap, and very high decay index values ≥43; however, support for relationships between and within these clades was often much lower (Fig. 2). Major clades that received high support included Danaea,Christensenia, and Angiopteris s.l. (i.e., including Archangiopteris and Macroglossum). Marattia s.l. was resolved as paraphyletic in all analyses, regardless of optimality criterion or inclusion of morphology. Marattia s.s., i.e., the clade that includes the type Marattia alata Sw., was resolved as sister to Angiopteris s.l. (91% Bayesian posterior probability, 81% MP bootstrap, decay index = 6). Two groups of taxa generally treated in Marattia s.l., the neotropical Marattia laevis complex, which corresponds to Eupodium J. Sm., and a clade of the paleotropical taxa, were both strongly supported as monophyletic; these two clades (referred to as "Marattia A" and "Marattia B," respectively) were resolved as sisters in all molecular analyses; however, support for this relationship was lower (97% Bayesian posterior probability, 68% MP bootstrap, decay index = 6).

While MP analyses of data set {1} (with or without inclusion of outgroup taxa) and trnSGG alone were congruent with the topology produced by all ML and BI analyses, other MP analyses produced an alternate topology. In MP analyses of data set {2}, and rps4–trnS alone, Christensenia was resolved as sister to a clade of Marattia s.l. + Angiopteris s.l. Relaxation of parsimony by only a single step is all that is required to produce the ML topology, indicating the low level of support for either placement of Christensenia under MP. Apart from the change in position of Christensenia, all other relationships were congruent with the topology shown in Fig. 2, with bootstrap support values varying slightly but without showing any trend for increased or decreased support.

Relative to related lineages, even the longest branches within Marattiaceae are short (Fig. 3), suggesting either recent radiation or reduced rates of molecular evolution in the family relative to other lineages, as proposed by Soltis et al. (2002)Go. While trnSGG proved to be the most variable region used in this study, this region is far more conserved in Marattiaceae than in other groups, including the filmy ferns (Hymenophyllaceae) (J. Nitta, University of California, Berkeley, unpublished manuscript); all regions used in this study proved to be much less informative than would be expected given the variation found in numerous studies of much younger lineages using the same markers. Tests for the presence of a molecular clock in the ingroup were rejected by a likelihood ratio test for all genes separately and combined. The relatively short length of branches, in combination with the distance from the nearest outgroups, likely contribute to the difficulty in rooting the phylogeny (see section Rooting).


Figure 3
View larger version (17K):
[in this window]
[in a new window]

 
Fig. 3. Maxiumum likelihood (ML) phylogram of Marattiaceae and related lineages based on combined analysis of atpB+rbcL+rps4{3}, showing short branch lengths within Marattiaceae. Root position shown is recovered when Equisetum is included in the outgroup and when the GTR+I+G model is implemented in ML; root position within Marattiaceae varies with outgroup and model choice.

 
Morphology
The MP analysis of the morphological data set produced a topology that is incongruent with molecular trees. While many elements agree with the molecular results (e.g., Marattia paraphyletic; Danaea, Christensenia,Eupodium,Ptisana, and Angiopteris monophyletic), the interrelationships of these groups differed, albeit with very few characters supporting any branch. Character descriptions, morphological matrix, and cladogram are included in Appendix S1 (see Supplemental Data with online version of this article).

Rooting
Rooting of the phylogeny was inconsistent between data sets, model choice, and optimality criteria. Additionally, both rooting and ingroup topology were affected by outgroup choice and analytical method. Figure 4 summarizes the rooting positions resolved by various analyses on the combined atpB+rbcL+rps4 data set. As a convention, root positions in this paper are referred to by the clade subtended by the branch to which the root attaches, e.g., "the branch subtending Danaea" and will be followed by the root position number in Fig. 5 in brackets, e.g., [1].


Figure 4
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 4. Unrooted network, with topology and branch lengths from maxiumum likelihood (ML) analysis of ingroup taxa. Root positions reconstructed by various analyses of atpB+rbcL+rps4 {3} with three possible outgroups (Equisetaceae [Equ], Ophioglossaceae [Oph], or Osmundaceae [Osm]) are indicated with arrows. Numbers correspond to root positions in Fig. 5. Rooting with Equisetaceae or any larger group containing Equisetum produces the same rooting results.

 

Figure 5
View larger version (26K):
[in this window]
[in a new window]

 
Fig. 5. The 41 possible root attachment points compared in Fig. 6. Root position 1, shown in a black circle, is the consensus root position.

 
Parsimony rooting
With the exception of atpB analyzed separately, MP consistently resolved the root on the branch subtending Danaea [1], the longest internal branch, in analyses of data set {1}, data set {3}, rbcl alone, and rps4 alone. To test for the possibility of long-branch attraction (LBA), the ingroup portion of data set {3} was rooted with 100 random outgroups with base composition equal to the average composition of the ingroup taxa. In all 100 trials, MP attached the root to the branch subtending Danaea [1]. In 76% of the random outgroup rooting trials, Christensenia was resolved on the branch between Danaea and the rest of the taxa (i.e., sister to all other taxa), but in 24% of the trials Christensenia was separated from Danaea and placed in the same position recovered by ingroup analyses (i.e., in a clade with Marattia s.s. and Angiopteris). When Danaea was excluded from the analysis, MP resolved the root on the branch subtending Christensenia [24] in 97% of the trials. Lundberg rooting (Lundberg, 1972Go) agrees with all other MP rooting results, attaching the root to the branch subtending Danaea [1], as do results of analyses in which third position nucleotides were downweighted or excluded.

Wilcoxon signed-rank tests (Templeton, 1983Go), comparing the parsimony tree length scores of 41 constraint trees corresponding to all possible root attachment points within data set {3}, were performed with three possible outgroups (Equisetaceae, Ophioglossaceae, Osmundaceae) (Figs. 5, 6). In the three analyses, only two branches were found to be not significantly less parsimonious (p≤ 0.05) than the MP root placement (the branch subtending Christensenia [24] and the branch subtending "Marattia A" [7]).


Figure 6
View larger version (25K):
[in this window]
[in a new window]

 
Fig. 6. The 41 possible root positions for the ingroup topology based on atpB+rbcL+rps4 {3} tested with three possible outgroups (A–C). Values on the x-axis correspond to the 41 root positions in Fig. 5 ranked by likelihood. Root positions that were significantly less optimal than the ML root position using a Shimodaira–Hasegawa test are marked with an asterisk (*). Root positions not significantly less parsimonious than the most parsimonious position using a Wilcoxon signed-rank test are marked with a plus-sign (+). "ML" and "MP" indicate the maximum likelihood and most parsimonious root positions.

 
Likelihood and Bayesian rooting
The unconstrained ML analysis of data set {3} under a GTR+G+I model of evolution attached the root to the branch subtending Christensenia [24] when Equisetum is included in the outgroup. However, when Equisetum is excluded from the analysis, rooting changes. In global analyses, ML and BI analyses of data set {3} resolved Marattiaceae and Equistaceae as a clade (99% posterior probability) sister to leptosporangiate ferns (98% posterior probability), in agreement with other previous studies (Pryer et al., 2001, 2004; Wikström and Pryer, 2005Go). Despite the fact that Equisetaceae and Marattiaceae share a most recent common ancestor in these analyses, the Equisetum crown group is on such a long branch that Osmundaceae, Ophioglossaceae, and Psilotaceae are closer to Marattiaceae in terms of branch length. When either Osmundaceae or Ophioglossaceae are used as outgroups, unconstrained ML analyses (under a GTR+G+I model of evolution) attach the root to the branch subtending Marattia s.s. [26]. Under simpler models of evolution, however, rooting became consistent, with the root attached to the branch below Danaea (in agreement with MP) in all outgroup combinations with both gene data sets. Midpoint rooting and enforcement of a molecular clock assumption on both gene data sets also resolved the root in agreement with MP. The BI analyses of both gene data sets resolved a basal polytomy (not shown).

As a comparison to the MP random rooting analysis, the ingroup portion of data set {3} was rooted with 20 random outgroups with base composition equal to the average composition of the ingroup taxa and analyzed under ML. In all trials, multiple equally likely trees were recovered, resulting in a total of 313 trees for the 20 trial runs. Each of the 41 possible root positions (Fig. 5) was resolved at least once; frequency of root attachment resolution was nonrandom and scaled approximately to branch length.

Shimodaira–Hasegawa (SH) tests, comparing the likelihood scores of 41 constraint trees corresponding to all possible root attachment points within data set {3}, were performed with three possible outgroups (Equisetaceae, Ophioglossaceae, Osmundaceae); the results are shown in Figs. 5 and6 In the SH tests the majority of root positions were not significantly less optimal than the ML root position. Analyses rooted with the nearest outgroup, Osmundaceae, could exclude 14 of 41 possible root positions (p ≤ 0.05); Ophioglossaceae and Equisetaceae excluded six and four root positions, respectively, all of which were redundant with those excluded by rooting with Osmundaceae. In all tests, the root position resolved by MP was not significantly less optimal than the ML root position.

DISCUSSION

Phylogenetic relationships
While it was hoped that sequence data from Christensenia might break up the long branch leading to marattioid ferns and help resolve relationships among the early-diverging fern lineages, this was not the case. Christensenia could be sister to all other marattioid ferns, depending on root placement (see Discussion, section Rooting) however, nucleotide sequences from Christensenia are highly similar to those from other marattioids, and Christensenia is in no way a genetic "missing link" between marattioids and other early-diverging fern lineages.

Long-branch attraction, but not long-branch repulsion (LBR), appears to be affecting the analyses: the ML ingroup topology was consistent with or without outgroups, and the MP topology was either congruent with the ML topology or produced an alternate (but equally poorly supported) topology when outgroups were included. If LBR were affecting the ML analysis, MP analysis would not be expected to produce the same topology because parsimony has the opposite bias. Additionally, removing Danaea and/or Christensenia from the analysis had no effect on the ML topology. Under MP, real outgroups performed no differently than random outgroups. While this result does not indicate that the root chosen by MP is incorrect, it does provide evidence that LBA may be occluding the true signal. In contrast, when random outgroups were used to root the ingroup under ML, each of the 41 possible root positions was selected at least once, with a weak tendency to attach the root to longer internal branches. While LBA may be biasing MP analyses, the root position selected by MP corresponds to that chosen by midpoint rooting, enforcement of molecular clock assumptions under ML, MP Lundberg rooting, and morphology (Rothwell and Nixon, 2006Go; Murdock et al., 2006aGo); it is this consensus rooting that is shown in Fig. 2.

The phylogeny of marattioid ferns presented here has several important impacts on the understanding of the taxonomy of marattioids. Marattiaceae s.l. is resolved as monophyletic with strong support, in agreement with earlier phylogenies (e.g., Pryer et al., 2001; Des Marais et al., 2003Go; Pryer et al., 2004; Wikström and Pryer, 2005Go). Marattia s.l. is paraphyletic based on the evidence presented in this study, and is resolved in all analyses (excluding morphology alone) as three well-supported clades: Marattia s.s. (neotropical taxa and M. douglasii of Hawaii), "Marattia A," which corresponds to Eupodium (the neotropical Marattia laevis complex with stalked synangia), and "Marattia B," a clade comprising the paleotropical taxa. A taxonomic revision of Marattiaceae based on the phylogenetic results presented here, including new combinations and a more detailed discussion of nomenclature can be found in a forthcoming paper (Murdock, 2008).

The resolution of Danaea as sister to the other marattioid ferns (i.e., the consensus root position shown in Fig. 2) agrees with all previous molecular phylogenies (Pryer et al., 2001, 2004; Des Marais et al., 2003Go; Murdock, 2005Go; Wikström and Pryer, 2005Go; Rothwell and Nixon, 2006Go; Qiu et al., 2006Go; Schuettpelz et al., 2006Go), but disagrees with the morphological analysis of Hill and Camus (1986Go). The sampling in this study is insufficient to make any clear conclusions about relationships of species groups within Danaea; Christenhusz et al. (2008)Go, contains more comprehensive sampling within Danaea than is presented here.

Christensenia has been previously described as the most primitive of the marattioid ferns (e.g., Campbell, 1911Go; Hill and Camus, 1986Go); the results of this study do not support this hypothesis. Christensenia is certainly distinctive, and the unique morphology of this genus is mirrored in the nucleotide sequences: while similar to Angiopteris and Marattia s.l., the sequences of Christensenia contain multiple autapomorphies, mostly in the form of unique indels that provide no grouping information within the family. The origin of Christensenia from within a clade of pinnately compound, free-veined taxa may seem surprising, but is the inevitable conclusion unless Christensenia was resolved as sister to all other marattioid taxa. Asama (1960)Go and Hill and Camus (1986)Go hypothesized that the leaf morphology and venation of Christensenia resulted from the compression of a rachis of a complex pinnate frond; the position in the marattioid phylogeny and the disposition of sori in irregular rows on both sides of main vein branches support this hypothesis. Christensenia has likely undergone significant morphological change as adaptation to extreme low-light environments: low, creeping habit; large, raised stomata; and high levels of chlorophyll concentrated in cells on the adaxial leaf surface (with abaxial surface cells possibly providing internal light reflection) (Nasrulhaq-Boyce and Mohamed, 1987Go; Rolleri, 1993Go).

Because of the more or less free sporangia and trilete spores, it might be presumed that Angiopteris represents a plesiomorphic form among extant marattioids. Bierhorst (1971, p. 368)Go hypothesized the opposite case, arguing that the semifree state seen in Angiopteris sporangia is secondarily derived from a fully synangiate state: "The sporangia of Angiopteris, although discrete at maturity, have no morphological integrity when initiated.... In the Marattiales, the difference between free sporangia and synangia is directly attributable to differential growth of the intersporangial areas." The phylogeny presented here supports Bierhorst’s hypothesis: partially free sporangia and trilete spores are a derived state in extant marattioids. Fossil marattioids had both monolete and trilete spores (some fossils showing a mixture of the two in the same sporangium), and fused to free sporangia. Because monolete, ellipsoid spores are generally held to be a derived state in fern lineages, with the tetrahedral spores of the lycophyte and bryophyte lineages (and the trilete marks found in the majority of these lineages) assumed to be the plesiomorphic state, the trilete, tetrahedral spores in Angiopteris may represent a reversal to the ancestral state (Tryon and Lugardon, 1991Go; Kenrick and Crane, 1997Go). However, given that basal fern lineages have both monolete and trilete spores, and fossil marattioids also produce both types of spores, the polarity of the character state change in spore morphology is difficult to interpret. The finding that Angiopteris is monophyletic only if it is circumscribed to include Macroglossum and Archangiopteris (which also encompasses the rarely recognized Protomarattia and Protangiopteris) agrees with the results of Hill and Camus (1986)Go, and the morphological arguments presented in Camus (1989)Go and Chang (1975)Go.

Marattia douglasii
Hooker and Arnott (1832)Go originally identified Marattia douglasii as M. alata, known only from the American tropics. Later authors either considered this species to represent an entirely new genus, Stibasia (e.g., Presl, 1845Go; de Vriese and Harting, 1853Go), or to be more closely related to the M. melanesica complex in Southeast Asia (Kuhn, 1889Go; Copeland, 1947Go). Nucleotide sequence data support Hooker and Arnott’s original interpretation that M. douglasii is related to M. alata. The most likely biogeographic scenario that led to this distribution pattern is that the progenitor of M. douglasii dispersed to Hawaii from the American tropics, possibly carried by a tropical storm track (Ramp Neale et al., 2007Go). However, given the current sampling and topology, it remains a remote possibility that Marattia s.s. originated in Hawaii and dispersed to the neotropics. Further sampling within Marattia s.s. and clarification of fossil affinities may help answer this question.

Marattia purpurascens
One of the rarest species in the family, M. purpurascens, is endemic to Ascension Island and is known from only a few remaining populations (Cronk, 1980Go; Gray et al., 2005Go). While morphologically similar to the African M. fraxinea, the sequences from this species are clearly distinct from the South African M. fraxinea sequenced in this study. It is likely, given the wide range of morphological variation in M. fraxinea, that several species are subsumed under this name; while this has been previously proposed, species delimitation has proven difficult due to intermediate and overlapping morphologies (Pichi Sermolli, 1969Go). Because Ascension Island is quite young, ~1 million years old, and the habitat for M. purpurascens on Ascension Island is probably much younger than that, given the relatively recent volcanic activity (Nielson and Sibbett, 1996Go; Ashmole and Ashmole, 2000Go), it would be remarkable (albeit not impossible) for the observed sequence divergence to have arisen in that time. Additional sampling from within the M. fraxinea complex, particularly from West African populations, is needed to clarify the origin of M. purpurascens and other lineages in the M. fraxinea complex.

Marattia howeana
Marattia howeana, another rare, endangered species, is endemic to Lord Howe Island. Threatened due to habitat destruction and grazing by introduced pigs, the few remaining known populations of M.howeana are being protected and monitored. While treated in the past as a variety of M. fraxinea, M. howeana is both morphologically and genetically distinct from M. fraxinea as well as nearby Australasian species.

Angiopteris evecta
This species, the type of Angiopteris, is extremely widespread and encompasses a wide range of morphological variation. The morphological plasticity of this genus is such that nearly all geographic and taxonomic distinctions break down upon close scrutiny. While there are indications of phylogenetic and biogeographic structure within the A. evecta complex, the markers used for this study are insufficiently informative to clearly elucidate relationships in this clade. Angiopteris longifolia, treated in most floras as a synonym of A. evecta, is distinguished mainly by having sori extending to the apex of the pinnule (as opposed to stopping well below the apex in A. evecta). In the Society Islands, these two forms can be seen in mixed populations with no other apparent differences (A. Murdock, personal observation). The form introduced to the Hawaiian Islands, purportedly from Samoa, is comparable to A. longifolia. While there is little basis for maintaining these as separate species, particularly given the range of morphology that A. evecta encompasses (as currently understood), it is interesting to note that MP analysis groups A. "longifolia" from Tahaa (Society Islands) with the three accessions from the Hawaiian islands, albeit with poor support (Fig. 2).

Morphology
Morphology suffers from the same general problem as the molecular data: the major clades are highly supported by multiple characters; however, there are very few characters that can inform interrelationships of these clades, and there is often homoplasy in these few characters. While inclusion of morphology in phylogenetic analyses has increased phylogenetic resolution in other plant groups (e.g., Doyle, et al., 1994Go; Mishler et al., 1994Go), in Marattiaceae the clades that received strong support from the morphology and microstructural data are clades that were already strongly supported by sequence data.

Rooting
Rooting can be problematic in cases where the clade in question is distantly related to the closest available outgroup and when internal branches of the clade are relatively short (Graham et al., 2002Go; Schuettpelz and Hoot, 2006Go); Marattiaceae is an excellent example of just such a clade. Because accurate root placement is essential for making conclusions about character evolution, estimating divergence dates, and for testing most other tree-based hypotheses, it is prudent to check for possible rooting biases even when relatively close outgroups are available.

Multiple methods for dealing with problematic rooting and long-branch attraction (LBA) have been proposed by previous authors (see e.g., Wheeler, 1990Go; Nixon and Carpenter, 1993Go1996Go; Mishler, 1994Go; Sanderson et al., 2000Go; Sanderson and Shaffer, 2002Go; Graham et al., 2002Go; Huelsenbeck et al., 2002Go). Lundberg (1972)Go suggested a method in which the ingroup is rooted using a vector statement representing a hypothetical ancestor. Lundberg rooting is of particular utility when outgroups are readily identifiable but are so distant that they are affecting the topology of the ingroup taxa; however, this is not really the case in the current study, as the ingroup topology is consistent with different outgroups, only the root position changes. Moreover, the composition of the Lundberg vector has a strong effect on rooting, and the best method of determining this vector is not always clear when there are multiple outgroups to choose from and their character composition is highly heterogeneous.

If a data set is clock-like, midpoint rooting or enforcing a molecular clock in a ML analysis are robust methods for selecting a root position. Huelsenbeck et al. (2002)Go found that even if data violates clock assumptions, enforcing a molecular clock may still identify the correct root. Midpoint rooting may choose the correct root under even wider ranges of rate heterogeneity (Swofford et al., 1996Go). However, with both midpoint rooting and molecular clock enforcement, if clock assumptions are violated by the data, confidence in these rootings is impossible to assess, and alternative criteria are always desirable. While nonreversible character change models and asymmetric step-matrices have been used for rooting trees, these methods do not appear to provide consistently accurate results and often require unjustifiable assumptions (Yang, 1994Go; Huelsenbeck et al., 2002Go).

Neither MP nor ML is immune from the effects of LBA; MP tends to be more strongly affected, although ML can suffer similarly if models are poorly suited to the data due to violated assumptions or mismatch in the selection process (Sullivan and Swofford, 2001Go; Swofford et al., 2001Go; Steel, 2005Go). Additionally, ML might in some cases produce positively misleading results in the form of long-branch repulsion (LBR) (Siddall, 1998Go; Siddall and Whiting, 1999Go). While simulation studies have disagreed on the question of whether MP or ML are better suited to dealing with the challenges of LBA, the general consensus is that MP will outperform ML when long-branch taxa are sister in the true phylogeny (Kolaczkowski and Thornton, 2004Go; Bergsten, 2005Go); as Pol and Siddall (2001)Go stated, it is not clear "whether parsimony outperforms likelihood for reasons other than succumbing to its own vices at the right time." The most effective measures for preventing LBA and LBR are to (1) use a model-based method with a good-fitting model, (2) increase taxon sampling as much as possible to break up the long branches, and (3) choose appropriately varying character sets in which the branches are not particularly long (Mishler, 2005Go). While model-based methods (e.g., ML and BI) are readily available and widely applied, we currently have very few tools to accurately assess a priori the fit of models to data and to determine how model mismatch might affect the results (but see Goremykin and Hellwig, 2006Go). Additionally, increasingly parameter-rich and better-fitting models are not necessarily better for phylogenetic inference (Tarrio et al, 2000Go; Steel, 2005Go; Kelchner and Thomas, 2006Go). Increased sampling will be maximally effective for rooting if additional taxa fall along long branches; additional sampling within crown groups will likely have little effect. Unfortunately, in many cases, including Marattiaceae, extinction has removed the possibility of sampling below the crown group for molecular analyses.

The rooting problem encountered in Marattiaceae is similar to that found in several previous studies (e.g., Tarrio et al., 2000Go; Wiens and Hollingsworth, 2000Go; Graham et al., 2002Go; Schuettpelz and Hoot, 2006Go); rooting techniques used in this paper may prove useful in similar situations regardless of the organisms in question. In Marattiaceae, the rate heterogeneity and the absence of any near outgroup are the most severe issues; nucleotide compositional differences between the ingroup and the outgroup do not appear to be a factor as they were in Drosophila (Tarrio et al., 2000Go). Because model selection and outgroup composition have a strong effect on root placement, methods such as midpoint rooting, enforcement of a molecular clock, or Lundberg rooting may be the best options for rooting the family. All three of these methods resolve the same root position, the branch subtending Danaea [1], which is the same position found by MP and ML under simpler models. In addition, when the more parameter-rich models used in ML select alternative root positions, these positions are never significantly more optimal than root position [1] on the branch subtending Danaea.

Other possible data sources for determining the root of any clade include fossil morphology and external evidence such as stratigraphy and historical biogeography. While not always available, these are valuable sources of information on rooting and character polarity that should not be overlooked. In the case of Marattiaceae, fossil morphology, as currently understood, does not provide clear rooting information because of large amounts of missing data, seemingly random distribution of character states in fossil taxa, and the absence of Tertiary fossils that could be more directly compared with extant taxa. If the relationships of fossil marattioids and polarity of character change were more clearly understood, it is likely that more information on the rooting of the extant marattioid ferns could be obtained. Additionally, genomic structural data or gene duplications may provide clearer rooting information for the family in the future.

Appendix 1. Voucher information and GenBank accession numbers for taxa used in this study. A dash (—) indicates the region was not sampled. Voucher specimens are deposited in the following herbaria: UC = University of California Berkeley; TUR = University of Turku, Finland; NY = New York Botanic Garden; PTBG = National Tropical Botanical Garden, Lauai, Kauai; K = Royal Botanic Gardens, Kew, E = Royal Botanic Gardens, Edinburgh. Place of origin indicated for all taxa, including plants collected from cultivation for this study. Plants with AGM collection numbers were collected by the author in the field or are from the living collections of the University of California Botanical Garden, the private collections of R. Whitehead, or Humboldt State University (HSC).


View this table:
[in this window]
[in a new window]

 
 

FOOTNOTES

1 The author thanks B. Mishler, A. Smith, K. Will, T. Carlson, and others for manuscript comments and plant material: M. Frantz, J. Strother, S. Lin, M. Lehnert, J. Game, the Mishler and Baldwin laboratories, D. Kelch, K. Pryer, E. Schuettpelz, H. Schneider, J. Metzgar, M. Windham, N. Nagalingum, P. Korall, A. Grusz, G. Theseira, Forest Research Institute Malaysia, D. Palmer, P. Bily and the Nature Conservancy (Maui), M. Christenhusz, D. Lorence and N.T.B.G., S. Stroud and the Ascension Island Conservation Centre, M. A. H. Mohamed and University of Malaya, R.B.G. Kew, R.B.G. Edinburgh, D. Walker, B. Weigle, R. Whitehead, K. Roux and S.A.N.B.I., Xishuangbana Botanic Garden, University of California Botanical Garden, T. Ranker, S. Graham, H. Rai, T. Motley, and D. Barrington. This research was a portion of the author’s doctoral dissertation research, which was supported by NSF Doctoral Dissertation Improvement Grant (DEB-0608497), NSF ATOL Grant (DEB-0228729), University of California Pacific Rim Foundation, Polynesia Education and Research Laboratories Research Fellowship, and the University of California, Berkeley Department of Integrative Biology Summer Research and Hansen Travel Grants. Back

2 E-mail: murdock{at}berkeley.edu Back

LITERATURE CITED

Asama, K. 1960. Evolution of leaf forms through the ages explained by the successive retardation and neoteny. Science Reports, Tohoku University, Sendai, Japan, 2nd series (Geology), vol. 4: 252–280.

Ashmole, P., AND M. Ashmole. 2000. St. Helena and Ascension Island: A natural history. Oswestry, Shropshire, UK.

Backer, C. A. 1929. The problem of Krakatao as seen by a botanist. Visser and Co., Weltevreden, Java; Martinus Nijhoff, The Hague, Netherlands.

Bergsten, J. 2005. A review of long-branch attraction. Cladistics 21: 163–193.[CrossRef][Web of Science]

Bierhorst, D. W. 1971. Morphology of vascular plants. Macmillan, New York, New York, USA.

Brownsey, P. J., AND J. C. Smith-Dodsworth. 2000. New Zealand Ferns and Allied Plants, 2nd ed. David Bateman, Auckland, New Zealand.

Buckley, T. R., C. Simon, H. Shimodaira, AND G. K. Chambers. 2001. Evaluating hypotheses on the origin and evolution of the New Zealand alpine cicadas (Maoricicada) using multiple-comparison tests of tree topology. Molecular Biology and Evolution 18: 223–234.[Abstract/Free Full Text]

Burleigh, J. G., AND S. Mathews. 2004. Phylogenetic signal in nucleotide data from seed plants: Implications for resolving the seed plant tree of life. American Journal of Botany 91: 1599–1613.[Abstract/Free Full Text]

Campbell, D. H. 1911. The Eusporangiatae: The comparative morphology of the Ophioglossaceae and Marattiaceae. Carnegie Institution of Washington, Publication 140, New York, New York, USA.

Camus, J. M. 1989. The limits and affinities of Marattialean fern genera in China and the west Pacific. In K. H. Shing, and K. U. Kramer [eds.], Proceedings of the International Symposium on Systematic Pteridology, Beijing, China, 1988. 31–37. China Science and Technology Press, Beijing, China.

Chang, C.-Y. 1975. Morphology of Archangiopteris Christ. &Gies. and its relationship with Angiopteris Hoffm. Acta Botanica Sinica 15: 235–247 [English translation of Acta Botanica Sinica 15: 261–276. 1973].

Chiang, T. Y., Y. C. Chiang, C. H. Chou, Y. P. Cheng, AND W. L. Chiou. 2002. Phylogeography and conservation of Archangiopteris somai Hayata and A. itoi Shieh based on nucleotide variation of the atp B-rbc L intergenic spacer of chloroplast DNA. Fern Gazette 16: 335 (Abstract).

Ching, R. C. 1959. Pteridophyta: Ophioglossaceae–Oleandraceae. In S. Chien, and W. Chun [eds.], Flora Reipublicae Popularis Sinicae, vol. 2. Academia Sinica, Peking [Beijing], China.

Christenhusz, M. J. M., H. Tuomisto, J. S. Metzgar, AND K. M. Pryer. 2008. Evolutionary relationships within the neotropical, eusporangiate fern genus Danaea (Marattiaceae). Molecular Phylogenetics and Evolution 46: 34–48.[CrossRef][Web of Science][Medline]

Collinson, M. E. 2001. Cainozoic ferns and their distribution. Brittonia 53: 173–235.[Web of Science]

Copeland, E. B. 1947. Genera filicum: The genera of ferns. [Annales Cryptogamici et Phytopathologici 5]. Chronica Botanica, Waltham, Massachusetts, USA.

Cronk, Q. C. B. 1980. Extinction and survival in the endemic vascular flora of Ascension Island. Biological Conservation 17: 207–219.[CrossRef][Web of Science]

De Vriese, W. H., AND P. Harting. 1853. Monographie des Marattiacées. Noothoven van Goor, Arnz Co., Leiden, Netherlands.

Des Marais, D. L., A. R. Smith, D. M. Britton, AND K. M. Pryer. 2003. Phylogenetic relationships of extant horsetails, Equisetum, based on chloroplast DNA sequence data (rbcL and trnL-F). International Journal of Plant Sciences 164: 737–751.[CrossRef][Web of Science]

Don, R. H., P. T. Cox, B. J. Wainwright, K. Baker, AND J. S. Mattick. 1991. ‘Touchdown’PCR to circumvent spurious priming during gene amplification. Nucleic Acids Research 19: 4008.[Free Full Text]

Doyle, J. A., M. J. Donoghue, AND E. A. Zimmer. 1994. Integration of morphological and ribosomal RNA data on the origin of angiosperms. Annals of the Missouri Botanical Garden 81: 419–450.[CrossRef][Web of Science]

Embley, T. M., AND W. Martin. 2006. Eukaryotic evolution, changes and challenges. Nature 440: 623–630.[CrossRef][Medline]

Felsenstein, J. 1981. Evolutionary trees from DNA sequences: A maximum likelihood approach. Journal of Molecular Evolution 17: 368–376.[CrossRef][Web of Science][Medline]

Fu, L. K., AND J. M. Jin [eds.] 1992. China plant red data book: Rare and endangered plants, vol. 1. Science Press, Beijing, China.

Giribet, G. 2002. Current advances in the phylogenetic reconstruction of metazoan evolution. A new paradigm for the Cambrian explosion? Molecular Phylogenetics and Evolution 24: 345–357.[CrossRef][Web of Science][Medline]

Goldman, N., J. P. Anderson, AND A. G. Rodrigo. 2000. Likelihood-based tests of topologies in phylogenetics. Systematic Biology 49: 652–670.[Abstract/Free Full Text]

Goremykin, V. V., AND F. H. Hellwig. 2006. A new test of phylogenetic model fitness addresses the issue of the basal angiosperm phylogeny. Gene 381: 81–91.[CrossRef][Medline]

Graham, S. W., R. G. Olmstead, AND S. C. H. Barrett. 2002. Rooting phylogenetic trees with distant outgroups: a case study from the commelinoid monocots. Molecular Biology and Evolution 19: 1769–1781.[Abstract/Free Full Text]

Graham, S. W., P. A. Reeves, A. C. E. Burns, AND R. G. Olmstead. 2000. Microstructural changes in noncoding chloroplast DNA: Interpretation, evolution, and utility of indels and inversions in basal angiosperm phylogenetic inference. International Journal of Plant Sciences 161 supplement: S83–S96.[CrossRef][Web of Science]

Gray, A., T. Pelembe, AND S. Stroud. 2005. The conservation of the endemic vascular flora of Ascension Island and threats from alien species. Oryx 39: 1–5.[CrossRef]

Hamilton, M. B. 1999. Tropical tree gene flow and seed dispersal: Deforestation affects the genetic structure of the surviving forest fragments. Nature 401: 129–130.[CrossRef][Web of Science]

Hasebe, M., M. Ito, R. Kofuji, K. Ueda, AND K. Iwatsuki. 1993. Phylogenetic relationships of ferns deduced from rbcL gene sequence. Journal of Molecular Evolution 37: 476–482.[CrossRef][Web of Science][Medline]

Hill, C. R., AND J. M. Camus. 1986. Evolutionary cladistics of marattialean ferns. Bulletin of the British Museum (Natural History). Botany. 14: 219–300.

Holttum, R. E. 1978. The morphology and taxonomy of Angiopteris (Marattiaceae) with description of a new species. Kew Bulletin 32: 587–594.[CrossRef]

Hooker, W. J., AND G. A. W. Arnott. 1832. The botany of Captain Beechey’s voyage; comprising an account of the plants collected by Messrs. Lay and Collie, and other officers of the expedition, during the voyage to the Pacific and Behring’s Strait, performed in His Majesty’s ship Blossom, under the command of Captain F. W. Beechey ...in the years 1825, 26, 27, and 28, part 3. H. G. Bohn, London, UK.

Hooker, W. J., AND J. G. Baker. 1868. Synopsis filicum; or a synopsis of all known ferns, including the Osmundaceae, Schizaeaceae, Marattiaceae, and Ophioglossaceae (chiefly derived from the Kew Herbarium). Robert Hardwicke, London, UK.

Hsu, T. W., S. J. Moore, AND T. Y. Chiang. 2000. Low RAPD polymorphism in Archangiopteris itoi, a rare and endemic fern in Taiwan. Botanical Bulletin of Academia Sinica 41: 15–18.[Web of Science]

Huelsenbeck, J. P., AND F. Ronquist. 2001. MRBAYES, Bayesian inference of phylogenetic trees. Bioinformatics (Oxford, England) 17: 754–755.[CrossRef]

Huelsenbeck, J. P., J. P. Bollback, AND A. M. Levine. 2002. Inferring the root of a phylogenetic tree. Systematic Biology 51: 32–43.[Abstract/Free Full Text]

Kelch, D. G., A. Driskell, AND B. D. Mishler. 2004. Inferring phylogeny using genomic characters: A case study using land plant plastomes. In B. Goffinet, V. Hollowell, and R. Magill [eds.], Molecular systematics of bryophytes. Monographs in systematic botany from the Missouri Botanical Garden, vol. 98, 3–11. Missouri Botanical Garden Press, St. Louis, Missouri, USA.

Kelchner, S. A., AND M. A. Thomas. 2006. Model use in phylogenetics: Nine key questions. Trends in Ecology &Evolution 22: 87–94.[CrossRef][Web of Science][Medline]

Kenrick, P., AND P. R. Crane. 1997. The origin and early diversification of land plants: A cladistic study. Smithsonian Institution, Washington, D.

Kingston, N., S. Waldren, AND N. Smyth. 2004. Conservation genetics and ecology of Angiopteris chauliodonta Copel. (Marattiaceae), a critically endangered fern from Pitcairn Island, South Central Pacific Ocean. Biological Conservation 117: 309–319.[CrossRef][Web of Science]

Kolaczkowski, B., AND J. W. Thornton. 2004. Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature 431: 980–984.[CrossRef][Medline]

Kuhn, M. 1889. Farne. In H. G. A. Engler [ed.], Die Forschungsreise S.M.S. "Gazelle"in den Jahren 1874 bis 1876 unter Kommando des Kapitän zur See Freiherrn von Schleinitz herausgegeben von dem Hydrographischen Amt des Reichs-Marine-Amts. IV. Theil. Botanik. Ernst Siegfried Mittler und Sohn, Berlin, Germany.

Li, C.-X., AND S.-G. Lu. 2007. Phylogeny and divergence of Chinese Angiopteridaceae based on chloroplast DNA sequence data (rbcL and trnL-F). Chinese Science Bulletin 52: 91–97.[CrossRef][Web of Science]

Liu, Z.-H., J. Hilton, AND C.-S. Li. 2000. Review on the origin, evolution and phylogeny of Marattiales. Chinese Bulletin of Botany 17: 39–52.

Lundberg, J. G. 1972. Wagner networks and ancestors. Systematic Zoology 21: 398–413.[Abstract]

Maddison, D. R., AND W. P. Maddison. 2005. MacClade, version 4. 08. Sinauer, Sunderland, Massachusetts, USA.

Manhart, J. R. 1994. Phylogenetic analysis of green plant rbcL sequences. Molecular Phylogenetics and Evolution 3: 114–127.[CrossRef][Medline]

Manhart, J. R. 1995. Chloroplast 16S rDNA sequences and phylogenetic relationships of fern allies and ferns. American Fern Journal 85: 182–192.[CrossRef][Web of Science]

Millay, M. A. 1979. Studies of Paleozoic Marattialeans: A monograph of the American species of Scolecopteris. Palaeontographica. Abteilung B 169: 1–69.

Mishler, B. D. 1994. Cladistic analysis of molecular and morphological data. American Journal of Physical Anthropology 94: 143–156.[CrossRef][Web of Science][Medline]

Mishler, B. D. 2000. Deep phylogenetic relationships among "plants"and their implications for classification. Taxon 49: 661–683.[CrossRef][Web of Science]

Mishler, B. D. 2005. The logic of the data matrix in phylogenetic analysis. In V. A. Albert [ed.], Parsimony, phylogeny, and genomics, 57–70. Oxford University Press, Oxford, UK.

Mishler, B. D., L. A. Lewis, M. A. Buchheim, K. S. Renzaglia, D. J. Garbary, C. F. Delwiche, F. W. Zechman, T. S. Kantz, AND R. L. Chapman. 1994. Phylogenetics of the "green algae"and "bryophytes.". Annals of the Missouri Botanical Garden 81: 451–483.[CrossRef][Web of Science]

Murdock, A. G. 2005. Molecular evolution and phylogeny of marattioid ferns, an ancient lineage of land plants. Botany 2005: Annual Meeting of the Botanical Society of America, Austin, Texas, USA [online abstract]. Website http://www.2005.botanyconference.org/engine/search/index.php?func=detail&aid=41.

Murdock, A. G. In press. A taxonomic revision of the eusporangiate fern family Marattiaceae, with description of a new genus Ptisana. Taxon.

Murdock, A. G., A. R. Smith, B. D. Mishler, AND K. S. Renzaglia. 2006a. Evolutionary origin of the leptosporangiate ferns: Sorting through conflicting evidence. Botany 2006: Annual Meeting of the Botanical Society of America, Chico, California, USA [online abstract]. Website http://www.2006.botanyconference.org/engine/search/index.php?func=detail&aid=769.

Murdock, A. G., J. L. Reveal, AND A. Doweld. 2006b. (1746) Proposal to conserve the name Marattiaceae against Danaeaceae (Pteridophyta). Taxon 55: 1040–1042.

Nadot, S., R. Bajon, AND B. Lejeune. 1994. The chloroplast gene rps4 as a tool for the study of Poaceae phylogeny. Plant Systematics and Evolution 191: 27–38.[CrossRef][Web of Science]

Nasrulhaq-Boyce, A., AND M. A. Mohamed. 1987. Photosynthetic and respiratory characteristics of Malayan sun and shade ferns. New Phytologist 105: 81–88.[CrossRef][Web of Science]

Nickrent, D. L., C. L. Parkinson, J. D. Palmer, AND R. J. Duff. 2000. Multigene phylogeny of land plants with special reference to bryophytes and the earliest land plants. Molecular Biology and Evolution 17: 1885–1895.[Abstract/Free Full Text]

Nielson, D. L., AND B. S. Sibbett. 1996. Geology of Ascension Island, South Atlantic Ocean. Geothermics 25: 427–448.[CrossRef][Web of Science]

Nixon, K. C., AND J. M. Carpenter. 1993. On outgroups. Cladistics 9: 413–426.[CrossRef][Web of Science]

Nixon, K. C., AND J. M. Carpenter. 1996. On simultaneous analysis. Cladistics 12: 221–241.[CrossRef][Web of Science]

Nylander, J. A. A. 2004. MrModeltest v2. Program distributed by the author. Evolutionary Biology Centre, Uppsala University, Sweden.

Pichi Sermolli, R. E. G. 1969. Adumbratio Florae Aethiopicae. 16. Marattiaceae. Webbia 23: 329–351.[Medline]

Pichi Sermolli, R. E. G. 1977. Attempt at the classification of pteridophytes based on phylogeny. Webbia 31: 313–512.

Pickard, J. 1983. Rare or threatened vascular plants of Lord Howe Island. Biological Conservation 27: 125–139.[CrossRef][Web of Science]

Pol, D., AND M. E. Siddall. 2001. Biases in maximum likelihood and parsimony: A simulation approach to a 10-taxon case. Cladistics 17: 266–281.[CrossRef][Web of Science]

Posada, D., AND K. A. Crandall. 1998. Modeltest: Testing the model of DNA substitution. Bioinformatics (Oxford, England) 14: 817–818.[CrossRef]

Presl, C. B. 1845. Supplementum tentaminis pteridographiae, continens genera et species ordinum dictorum Marattiaceae, Ophioglossaceae, Osmundaceae, Schizaeaceae et Lygodiaceae. Haase, Prague, Czech Republic.

Pryer, K. M., H. Schneider, A. R. Smith, R. Cranfill, P. G. Wolf, J. Hunt, AND S. D. Sipes. 2001. Horsetails and ferns are a monophyletic group and the closest living relatives to seed plants. Nature 409: 618–622.[CrossRef][Medline]

Pryer, K. M., E. Schuettpelz, P. G. Wolf, H. Schneider, A. R. Smith, AND R. Cranfill. 2004. Phylogeny and evolution of ferns (monilophytes) with a focus on the early leptosporangiate divergences. American Journal of Botany 91: 1582–1598.[Abstract/Free Full Text]

Qiu, Y. L., L. Li, B. Wang, Z. Chen, V. Knoop, M. Groth-Malonek, O. Dombrovska, J. Lee, L. Kent, J. Rest, G. F. Estabrook, T. A. Hendry, D. W. Taylor, C. M. Testa, M. Ambros, B. Crandall-Stotler, R. J. Duff, M. Stech, W. Frey, D. Quandt, AND C. C. Davis. 2006. The deepest divergences in land plants inferred from phylogenomic evidence. Proceedings of the National Academy of Sciences, USA 103: 15511–15516.[Abstract/Free Full Text]

Rambaut, A., AND A. Drummond. 2003. Tracer v. 1.3: MCMC trace file analyser. Oxford University, Oxford, UK.

Ramp Neale, J. M., R. B. Neale, AND T. A. Ranker. 2007. Identifying the dominant source regions for the colonization of Hawaiian fern species using a trajectory and dispersal model. Botany and Plant Biology 2007 Joint Congress, Chicago, Illinois, USA [online abstract]. Website http://www.2007.botanyconference.org/engine/search/index.php?func=detail&aid=1722.

Rolleri, C. H. 1993. Revision of the genus Christensenia. American Fern Journal 83: 3–19.[CrossRef][Web of Science]

Rolleri, C. H. 2003. Caracteres diagnósticos y taxonomía en el género Angiopteris Hoffm. (Marattiaceae Bercht. &J. S. Presl) II, Sinopsis de las especies. Revista del Museo de La Plata. Sección Botánica 16: 1–23.

Rolleri, C. H., M. C. Lavalle, A. Mengascini, AND M. Rodríguez. 2003. Sistemática de los helechos maratiaceae [sic] (Marattiales-Marattiaceae). Revista del Museo de La Plata. Sección Botánica 16: 1–21.

Ronquist, F., AND J. P. Huelsenbeck. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics (Oxford, England) 19: 1572–1574.[CrossRef]

Roper, J. M., S. K. Hansen, P. G. Wolf, K. G. Karol, D. F. Mandoli, K. D. E. Everett, J. Kuehl, AND J. L. Boore. 2007. The complete plastid genome sequence of Angiopteris evecta (G. Forst.) Hoffm. (Marattiaceae). American Fern Journal 97: 95–106.[CrossRef][Web of Science]

Rothwell, G. W. 1999. Fossils and ferns in the resolution of land plant phylogeny. Botanical Review 65: 188–218.[CrossRef]

Rothwell, G. W., AND K. C. Nixon. 2006. How does the inclusion of fossil data change our conclusions about the phylogenetic history of euphyllophytes? International Journal of Plant Sciences 167: 737–749.[CrossRef][Web of Science]

Sanderson, M. J., AND H. B. Shaffer. 2002. Troubleshooting molecular phylogenetic analyses. Annual Review of Ecology and Systematics 33: 49–72.[CrossRef][Web of Science]

Sanderson, M. J., M. F. Wojciechowski, J.-M. Hu, T. S. Khan, AND S. G. Brady. 2000. Error, bias, and long-branch attraction in data for two chloroplast photosystem genes in seed plants. Molecular Biology and Evolution 17: 782–792.[Abstract/Free Full Text]

Schuettpelz, E., AND S. B. Hoot. 2006. Inferring the root of Isoëtes: Exploring alternatives in the absence of an acceptable outgroup. Systematic Botany 31: 258–270.[CrossRef][Web of Science]

Schuettpelz, E., P. Korall, AND K. M. Pryer. 2006. Plastid atpA data provide improved support for deep relationships among ferns. Taxon 55: 897–906.[Web of Science]

Schuettpelz, E., AND K. M. Pryer. 2007. Fern phylogeny inferred from 400 leptosporangiate species and three plastid genes. Taxon 56: 1037–1050.[Web of Science]

Schneider, H. 2007. Plant morphology as the cornerstone to the integration of fossils and extant taxa in phylogenetic analyses. Species. Phylogeny and Evolution 1: 65–71.

Shaw, J., E. B. Lickey, J. T. Beck, S. B. Farmer, W. Liu, J. Miller, K. C. Siripun, C. T. Winder, E. E. Schilling, AND R. L. Small. 2005. The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. American Journal of Botany 92: 142–166.[Abstract/Free Full Text]

Shaw, J., E. B. Lickey, E. E. Schilling, AND R. L. Small. 2007. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. American Journal of Botany 94: 275–288.[Abstract/Free Full Text]

Shimodaira, H., AND M. Hasegawa. 1999. Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Molecular Biology and Evolution 16: 1114–1116.[Web of Science]

Siddall, M. E. 1998. Success of parsimony in the four-taxon case: Long-branch repulsion by likelihood in the Farris Zone. Cladistics 14: 209–220.[CrossRef][Web of Science]

Siddall, M. E., AND M. F. Whiting. 1999. Long-branch abstractions. Cladistics 15: 9–24.[CrossRef][Web of Science]

Small, R. L., E. B. Lickey, J. Shaw, AND W. J. Hauk. 2005. Amplification of noncoding chloroplast DNA for phylogenetic studies in lycophytes and monilophytes with a comparative example of relative phylogenetic utility from Ophioglossaceae. Molecular Phylogenetics and Evolution 36: 509–522.[CrossRef][Web of Science][Medline]

Smith, A. R. 1995. Non-molecular phylogenetic hypotheses for ferns. American Fern Journal 85: 104–122.[CrossRef][Web of Science]

Smith, A. R., AND R. B. Cranfill. 2002. Intrafamilial relationships of the thelypteroid ferns (Thelypteridaceae). American Fern Journal 92: 131–149.[CrossRef][Web of Science]

Smith, A. R., K. M. Pryer, E. Schuettpelz, P. Korall, H. Schneider, AND P. G. Wolf. 2006. A classification for extant ferns. Taxon 55: 705–731.[Web of Science]

Soltis, P. S., D. E. Soltis, V. Savolainen, P. R. Crane, AND T. G. Barraclough. 2002. Rate heterogeneity among lineages of tracheophytes: Integration of molecular and fossil data and evidence for molecular living fossils. Proceedings of the National Academy of Sciences, USA 99: 4430–4435.[Abstract/Free Full Text]

Steel, M. 2005. Should phylogenetic models be trying to ‘fit an elephant’? Trends in Genetics 21: 307–309.[CrossRef][Web of Science][Medline]

Sullivan, J., AND D. L. Swofford. 2001. Should we use model-based methods for phylogenetic inference when we know that assumptions about among-site rate variation and nucleotide substitution pattern are violated? Systematic Biology 50: 723–729.[Free Full Text]

Swofford, D. L. 2002. PAUP*: Phylogenetic analysis using parsimony (*and other methods). Sinauer, Sunderland, Massachusetts, USA.

Swofford, D., G. Olsen, P. Waddell, AND D. M. Hillis. 1996. Phylogenetic inference. In D. Hillis, C. Moritz, and B. Mable [eds.], Molecular systematics, 2nd ed., 407–511. Sinauer, Sunderland, Massachusetts, USA.

Swofford, D. L., P. J. Waddell, J. P. Huelsenbeck, P. G. Foster, P. O. Lewis, AND J. S. Rogers. 2001. Bias in phylogenetic estimation and its relevance to the choice between parsimony and likelihood methods. Systematic Biology 50: 525–539.[Free Full Text]

Szövényi, P., Z. Hock, E. Urmi, AND J. J. Schneller. 2006. Contrasting phylogeographic patterns in Sphagnum fimbriatum and Sphagnum squarrosum (Bryophyta, Sphagnopsida) in Europe. New Phytologist 172: 784–794.[CrossRef][Web of Science][Medline]

Tarrio, R., F. Rodríguez-Trelles, AND F. J. Ayala. 2000. Tree rooting with outgroups when they differ in their nucleotide composition from the ingroup: The Drosophila saltans and willistoni groups, a case study. Molecular Phylogenetics and Evolution 16: 344–349.[CrossRef][Web of Science][Medline]

Taylor, T. N., AND E. L. Taylor. 1993. The biology and evolution of fossil plants. Prentice Hall, Englewood Cliffs, New Jersey, USA.

Templeton, A. R. 1983. Phylogenetic inferences from restriction endonucleases cleavage site maps with particular reference to the evolution of humans and apes. Evolution 37: 221–244.[CrossRef][Web of Science]

Tryon, A. F., AND B. Lugardon. 1991. Spores of the Pteridophyta: Surface, wall structure, and diversity based on electron microscope studies. Springer-Verlag, New York, New York, USA.

Warshauer, F. R. 1979. An overview of the feral pig problem in the Hawaii Volcanoes National Park. In R. M. Linn [ed.], Proceedings of the 2nd Conference of Scientific Research in National Parks, 1–5. National Parks and Wildlife Service, San Francisco, California, USA.

Wheeler, W. C. 1990. Nucleic acid sequence phylogeny and random outgroups. Cladistics 6: 363–367.[CrossRef][Web of Science]

Wiens, J. J., AND B. D. Hollingsworth. 2000. War of the iguanas: Conflicting molecular and morphological phylogenies and long-branch attraction in iguanid lizards. Systematic Biology 49: 143–159.[Abstract/Free Full Text]

Wikström, N., AND K. Pryer. 2005. Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails. Molecular Phylogenetics and Evolution 36: 484–493.[CrossRef][Web of Science][Medline]

Wolf, P. G. 1997. Evaluation of atpB nucleotide sequences for phylogenetic studies of ferns and other pteridophytes. American Journal of Botany 84: 1429–1440.[Abstract]

Wolf, P. G., K. G. Karol, D. F. Mandoli, J. Kuehl, K. Arumuganathan, M. W. Ellis, B. D. Mishler, D. G. Kelch, R. G. Olmstead, AND J. L. Boore. 2005. The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae). Gene 350: 117–128.[CrossRef][Web of Science][Medline]

Yang, Z. 1994. Estimating the pattern of nucleotide substitution. Journal of Molecular Evolution 39: 105–111.[Web of Science][Medline]


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Facebook Facebook   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?


This article has been cited by other articles:


Home page
Am. J. Bot.Home page
S. W. Graham and W. J. D. Iles
Different gymnosperm outgroups have (mostly) congruent signal regarding the root of flowering plant phylogeny
Am. J. Botany, January 1, 2009; 96(1): 216 - 227.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Murdock, A. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Murdock, A. G.
Agricola
Right arrow Articles by Murdock, A. G.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS