Multigene Phylogeny of Land Plants with Special Reference to Bryophytes and the Earliest Land Plants

Daniel L. Nickrent, Christopher L. Parkinson1,, Jeffrey D. Palmer and R. Joel Duff2,

*Department of Plant Biology and Center for Systematic Biology, Southern Illinois University at Carbondale; and
{dagger}Department of Biology, Indiana University at Bloomington


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
A widely held view of land plant relationships places liverworts as the first branch of the land plant tree, whereas some molecular analyses and a cladistic study of morphological characters indicate that hornworts are the earliest land plants. To help resolve this conflict, we used parsimony and likelihood methods to analyze a 6,095-character data set composed of four genes (chloroplast rbcL and small-subunit rDNA from all three plant genomes) from all major land plant lineages. In all analyses, significant support was obtained for the monophyly of vascular plants, lycophytes, ferns (including Psilotum and Equisetum), seed plants, and angiosperms. Relationships among the three bryophyte lineages were unresolved in parsimony analyses in which all positions were included and weighted equally. However, in parsimony and likelihood analyses in which rbcL third-codon-position transitions were either excluded or downweighted (due to apparent saturation), hornworts were placed as sister to all other land plants, with mosses and liverworts jointly forming the second deepest lineage. Decay analyses and Kishino-Hasegawa tests of the third-position-excluded data set showed significant support for the hornwort-basal topology over several alternative topologies, including the commonly cited liverwort-basal topology. Among the four genes used, mitochondrial small-subunit rDNA showed the lowest homoplasy and alone recovered essentially the same topology as the multigene tree. This molecular phylogeny presents new opportunities to assess paleontological evidence and morphological innovations that occurred during the early evolution of terrestrial plants.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
Understanding the evolutionary history of bryophytes (mosses, liverworts, and hornworts) is compelling because they represent the first green plants to colonize the terrestrial environment, in the course of which they evolved numerous important adaptations, including alternation of gametophytic (haploid) and sporophytic (diploid) generations, elaborate gametophytes, specialized gametangia, and desiccation-resistant spore walls (Renzaglia et al. 2000 Citation ). Spore microfossils of terrestrial (nonalgal) plants are known from the middle Ordovician (470 MYA) (Gray 1993Citation ); however, the earliest spore records identified as liverworts are from the lower Silurian (430 MYA) (Taylor 1995Citation ) and the earliest spore records for hornworts are from the upper Silurian (420 MYA) (Richardson 1985Citation ). Furthermore, unequivocal macrofossils of liverworts (Pallavicinites) and hornworts (Notothylacites) date to the upper Devonian (370 MYA) (Hueber 1961Citation ) and upper Cretaceous (70–90 MYA) (Nemejc and Pacltova 1974Citation ), respectively, i.e., long after their initial appearance. Determining which bryophyte group is the most ancient among extant land plants has proven difficult, and a myriad of phylogenetic hypotheses have been advanced with fundamentally different tree topologies. The long-standing concept of bryophyte monophyly was first challenged by morphological cladistic analyses which placed liverworts as the sister to all other land plants (Mishler and Churchill 1984Citation ). Analyses combining morphological characters with rRNA sequences (Mishler et al. 1994Citation ) and analyses of the chloroplast gene rbcL (Lewis, Mishler, and Vilgalys 1997Citation ) have generally supported a "liverworts-basal" topology (LBT) and placed either mosses or hornworts as sister to vascular plants (tracheophytes). More recently, the LBT was deemed the most parsimonious explanation to account for the presence and absence of three mitochondrial introns (Qiu et al. 1998Citation ).

Despite its wide acceptance, the LBT has been challenged by data derived from spermatozoid ultrastructure and spermatogenesis (Renzaglia and Duckett 1991Citation ) and, more recently (Garbary and Renzaglia 1998Citation ), by an extensive cladistic analysis of morphological and ultrastructural data that yielded a monophyletic moss-liverwort clade, with hornworts sister to all land plants. This "hornworts-basal" topology (HBT) was also obtained in analyses of nuclear small-subunit (nuSSU) rDNA (Hedderson, Chapman, and Rootes 1996Citation ), combined analysis of five chloroplast genes (Nishiyama and Kato 1999Citation ), and analyses of mitochondrial cox3 (Malek et al. 1996Citation ), nad5 (Beckert et al. 1999Citation ), and mitochondrial small-subunit (mtSSU) rDNA (Duff and Nickrent 1999 Citation ). In an effort to distinguish between the LBT and the HBT and to provide an improved understanding of land plant phylogeny, we generated additional sequence data to allow the construction of a multigene data set representing all three plant genomes and the same suite of exemplar taxa.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
Ingroup and Outgroup Selection
Sequences were obtained from chloroplast SSU (cpSSU) rDNA, chloroplast rbcL, nuSSU rDNA, and mtSSU rDNA for representatives of all major extant land plant lineages. Sampling included two hornworts (Phaeoceros and Megaceros), two liverworts (Marchantia and Calypogeia), three mosses (Takakia, Polytrichum, and Sphagnum), four lycophytes (Selaginella, Isoëtes, Huperzia, and Diphasiastrum), six ferns (Psilotum, Ophioglossum, Botrychium, Angiopteris, Equisetum, and Adiantum), six gymnosperms (Cycas, Ginkgo, Juniperus, Pinus, Gnetum, and Ephedra), and five angiosperms (Nymphaea, Magnolia, Zea, Nicotiana, and Glycine). Chara and Spirogyra (green algae) were used as outgroups. Taxa were initially selected based on the presence of one or more of the four gene sequences in GenBank, as well as a desire to represent a wide range of taxonomic diversity. From this list, sequences from 30 genera were obtained from GenBank or generated for this study. For 10 of the 28 ingroup taxa, the same species was used for all four gene sequences (table 1 ). Placeholder species were used for the remaining 18 taxa: 8 cases with two species, 11 cases with three species, and for Ephedra, four different species were used. Of the 119 DNA sequences assembled, 27 were newly generated. As previously reported (Duff and Nickrent 1999Citation ), PCR amplifications of mtSSU rDNA from Coleochaete and Selaginella failed; hence, Coleochaete was not used as an outgroup, and Selaginella represented the only case of a missing ingroup sequence in this study. Complete voucher information for newly generated sequences, full scientific names, and GenBank numbers for sequences of all taxa are given in table 1 .


View this table:
[in this window]
[in a new window]
 
Table 1 GenBank Accession Numbers for Taxa Used in Multigene Analysis of Land Plant Phylogeny

 
Molecular Methods
The protocols used here for extracting genomic DNA, PCR, cloning, and sequencing have previously been reported (Nickrent 1994Citation ; Qiu et al. 1998Citation ; Duff and Nickrent 1999Citation ; Chaw et al. 2000Citation ). The primers used to amplify and sequence the newly generated sequences were cited in the following: rbcL (Soltis and Soltis 1998Citation ), cpSSU rDNA (Nickrent, Duff, and Konings 1997Citation ), nuSSU rDNA (Nickrent and Starr 1994Citation ), and mtSSU rDNA (Duff and Nickrent 1999Citation ).

Alignments
Separate gene and multigene alignments were conducted manually using SeqApp (Gilbert 1993Citation ). Alignment of rbcL and cpSSU rDNA sequences was trivial, given essentially no length variation, whereas for nuclear and mitochondrial SSU rDNA sequences, alignments were guided by reference to published higher-order structures (Nickrent and Soltis 1995Citation ; Duff and Nickrent 1999Citation ). Terminal priming regions were excluded from all sequences, as were regions where alignment was ambiguous or where extensive length variation occurred (e.g., introns in mtSSU rDNA; see Duff and Nickrent 1999Citation ). There were (with reference to Glycine sequences) 51 such sites for rbcL, 63 for cpSSU rDNA, 133 for nuSSU rDNA, and 633 for mtSSU rDNA, for a total of 909 sites excluded from the multigene data set. The key to excluded regions and the final alignment (6,095 positions: rbcL = 1,351, cpSSU rDNA = 1,480, nuSSU rDNA = 1,720, and mtSSU rDNA = 1,544) can be found at http://www.science.siu.edu/landplants/Alignments/Alignments.html.

Parsimony and Likelihood Analyses
There currently exists conflicting information in the literature as to whether rbcL third codon positions are saturated, or, more specifically, whether the process of substitution has become so randomized as to preclude the use of this position for ascertaining deep phylogenetic relationships. Evidence that third-codon-position transitions are effectively saturated has been presented by Goremykin et al. (1996)Citation and Chaw et al. (2000)Citation . Conversely, Lewis, Mishler, and Vilgalys (1997)Citation showed that trees generated from only third positions are nearly identical to trees produced from all sites and that these trees agree with independent evidence. Hence, they concluded that, despite being highly variable, third positions retain phylogenetic signal across all land plants. A similar conclusion was reached by Källersjö et al. (1998)Citation . To test for saturation, we plotted uncorrected p-distances from first-, second-, and third-position transitions and transversions against corrected F84 distances (Felsenstein 1984Citation ).

Four versions of the multigene data set were constructed and analyzed with either maximum parsimony (MP) or maximum likelihood (ML). In the first version, all positions, including rbcL third codon positions, were included and were weighted equally. Analyses of this version were called MP+3Ti and ML-rates. The second version also included all positions, but differential weighting (ranging from 0.5 to 1.0) was applied to third codon positions. Analyses of this version were dubbed MP+weights. The third version excluded rbcL third-codon-position transitions, which was accomplished by recoding all third-position nucleotides as either R (A or G) or Y (C or T). Analyses of this version were called MP-3Ti and ML-3Ti. In the fourth version, all positions were included, and site-to-site rate variation was accounted for using the program DNArates (S. Pract, R. Overbeek, and G. Olsen, personal communication; available at the Ribosomal Database Project website). Analysis of this version was called ML+rates.

All MP analyses were conducted using PAUP* (Swofford 1998Citation ). For these analyses, the separate and combined gene matrices were analyzed either with all characters equally weighted or with unequal weighting (see above). Heuristic search strategies with random taxon addition (100 replicates), MULPARS on, and tree bisection-reconnection (TBR) branch swapping was employed. MP bootstrap (BS) support was assessed using PAUP*, employing a heuristic search strategy with 100 replications (1,000 for the multigene data set).

ML analyses employed fastDNAml, version 1.06, rewritten to run in parallel, on 8–16 nodes of the STARRS/SP cluster (IBM RS/6000 computer platform, located at Indiana University). The F84 model (Felsenstein 1984Citation ) as implemented in PHYLIP (Felsenstein 1993Citation ) was used, with the initial transition/transversion (ti/tv) ratio estimated using PUZZLE, version 4.02, under the Tamura-Nei model of evolution with parameter estimation set to "approximate" (Strimmer and von Haeseler 1996Citation ). Ten initial ML trees were inferred for each individual and the combined data set by randomizing "input" order with jumble and using "global" swapping across all nodes (equivalent to subtree-pruning-regrafting of PAUP*). The optimal tree (best log likelihood score) was then input into PAUP* to reoptimize the ti/tv ratio using a model that incorporated variability in rates of change. We used the F84 evolutionary model assuming a discrete gamma distribution with four categories of site-to-site rate variability. The resulting ti/tv ratio was used to infer a new tree as above, further optimizing branch lengths. This tree and the optimized ti/tv ratio were then used to estimate evolutionary rates of change for each sequence position by partitioning the sites into 35 "rate" categories using DNArates. A new ML tree, incorporating the rate categories and the reoptimized ti/tv ratio, was then inferred. This new optimal tree was then used for a second round of rate estimation and tree inference. This process was iterated until a stable topology was achieved.

ML bootstrapping with and without rates used SEQBOOT (PHYLIP) to generate 100 pseudoreplicate data sets that were analyzed using fastDNAml, and the resulting bootstrap numbers were generated using CONSENSE (PHYLIP). The individual trees from the 100 pseudoreplicates were analyzed as above using fastDNAml, randomizing the input order, and allowing swapping across all nodes. Each initial "no-rates" topology was then used to infer rate categories incorporating the optimized ti/tv ratio. A second round of rate category generation and tree inference was performed, and the final 100 trees were imported into CONSENSE to generate bootstrap node values.

Alternative topologies, evaluated under ML and MP, were determined using the parametric K-H test (Kishino and Hasegawa 1989) and the nonparametric Templeton (1983)Citation ranked-sums test as implemented in PAUP*. Decay analyses (Bremer 1994) were performed with MP using PAUP* and AUTODECAY (Eriksson 1998).


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
Composition Bias and Saturation
Nucleotide compositions for the entire multigene data set varied little among taxa (A = 0.253, C = 0.228, G = 0.292, T = 0.227), with an average A+T composition of 48%. A chi-square test of homogeneity of base frequencies across taxa gave 45.607 (df = 87, P = 0.9999); hence, nucleotide composition was homogeneous. Homogeneity was also seen for all of the separate gene partitions with the exception of rbcL. Chi-square tests of homogeneity showed that rbcL first and second codon positions were homogeneous, whereas third positions (with A+T averaging 72% and ranging from 46% in Selaginella to 85% in Marchantia) were significantly heterogeneous across taxa. These results agree with those reported for rbcL by Lewis, Mishler, and Vilgalys (1997)Citation and Nishiyama and Kato (1999)Citation and therefore require the use of sequence evolution models that do not assume compositional equality. In tests for saturation using distance plots (fig. 1 ), first- and second-position transitions and transversions and third-position transversions were significantly linear, whereas third-position transitions were significantly nonlinear. This result parallels those of Goremykin et al. (1996)Citation and Chaw et al. (2000)Citation and, we believe, indicates third-position transitional saturation for rbcL.



View larger version (30K):
[in this window]
[in a new window]
 
Fig. 1.—Plots of uncorrected pairwise sequence divergence (p-distances) versus corrected distances (F84-distance) for first, second, and third codon positions in rbcL. Transitional (ti) and transversional (tv) substitutions are combined for first and second positions, whereas they are plotted separately for third positions. Deviation from the diagonal trend line, as is evident for third-codon-position transitions, indicates saturation. A total of 427 pairwise comparisons are represented in each panel. Eight ingroup-to-algal-outgroup comparisons for third-position transitions were excluded because the observed dissimilarity exceeded 0.75

 
Phylogenetic Relationships Recovered Using All Analytical Methods
Of the 6,095 characters in the alignment, 1,037 (17%) were variable but parsimony-uninformative, and 1,219 (20%) were parsimony-informative. For the multigene data set, MP+3Ti analysis resulted in one tree of length 7,074 (fig. 2 ). MP-3Ti analysis yielded a single tree of length 5,252 (fig. 3 ). These two trees are highly congruent for relationships among tracheophytes but differed with regard to bryophyte relationships (see below).



View larger version (33K):
[in this window]
[in a new window]
 
Fig. 2.—Land plant relationships derived from an MP+3Ti analysis of the four-gene data set. A single tree (length = 7,074) was obtained with a consistency index (excluding uninformative sites) of 0.3900, a retention index of 0.4756, and a rescaled consistency index of 0.2332. Values above the branches are bootstrap percentages derived from 100 replications

 


View larger version (39K):
[in this window]
[in a new window]
 
Fig. 3.—Land plant relationships derived from an MP-3Ti analysis of the four-gene data set. A single tree (length = 5,251) was obtained with a consistency index (excluding uninformative sites) of 0.4312, a retention index of 0.5275, and a rescaled consistency index of 0.2941. Branch lengths shown are derived from the ML+rates analysis (log likelihood = -37,923.99574, ti/tv = 2.25). Numbers at nodes (left to right): ML-rates BS values, ML+rates BS values, MP-3Ti BS values, and decay index. Nodes not present are indicated by dashes

 
As shown in figures 2 and 3 and table 2 , MP+3Ti, MP-3Ti, ML-rates, and ML+rates analyses all yielded high BS values for the majority of land plant clades, including angiosperms, gymnosperms, seed plants, ferns (including Psilotum and Equisetum), lycophytes, tracheophytes, liverworts, and hornworts. When BS support >70% is considered significant (Hillis and Bull 1993Citation ), the following percentages of total nodes were supported as significant among the various analyses: MP+3Ti (70%), MP-3Ti (78%), ML-rates (74%), and ML+rates (81%). The MP and ML analyses gave strong support for a clade composed of all vascular plants (tracheophytes), although this was not the case among the nonmitochondrial genes analyzed individually. Moderate (MP) to high (ML) BS values were obtained for lycophyte monophyly. All methods gave high BS support for lycophytes as the sister to all other vascular plants, consistent with other molecular (Raubeson and Jansen 1992Citation ) and morphological (Kenrick and Crane 1997Citation ) analyses. The two heterosporous genera (Selaginella and Isoëtes) are sisters, as are the two homosporous genera (Huperzia and Diphasiastrum); these relationships generally follow current classifications (Kenrick and Crane 1997Citation ). The monophyly of ferns, here including the horsetail Equisetum and the whiskfern Psilotum, has high BS and decay analysis support. Within this clade, Psilotum is strongly supported as sister to the eusporangiate ferns Botrychium and Ophioglossum, a relationship now well documented (Manhart 1994Citation ; Pryer, Smith, and Skog 1995Citation ; Wolf 1997Citation ). Weak support exists for the placement of Angiopteris as sister to Equisetum and Adiantum, although this topology was also obtained in a combined analysis of rbcL and morphological characters (Pryer, Smith, and Skog 1995Citation ). Seed plants are clearly resolved as monophyletic, with moderate (MP) to strong (ML) support for a clade of extant gymnosperms. A clade of gnetophytes (Gnetum and Ephedra) and conifers (Pinus and Juniperus) is very strongly supported (98%–100% BS, decay value of 24 steps), in agreement with two recent multigene studies of seed plant phylogeny (Bowe, Coat, and dePamphilis 2000Citation ; Chaw et al. 2000Citation ). Angiosperms are strongly supported as monophyletic in all analyses and with all separate gene partitions (table 2 ).


View this table:
[in this window]
[in a new window]
 
Table 2 Bootstrapa Percentages for Selected Land Plant Clades

 
Phylogenetic Relationships that Differed Among the Various Analyses
Several differences can be seen when the MP+3Ti tree (fig. 2 ) and the MP-3Ti tree (fig. 3 ) are compared. Euphyllophytes (vascular plants exclusive of lycophytes; Kenrick and Crane 1997Citation ) are strongly supported using ML+rates and MP-3Ti. Different relationships among the five angiosperm taxa were obtained following analyses of the multigene data using MP and ML. For MP+3Ti analyses, Zea was sister to the rest, albeit with low BS support (fig. 2 ). The MP-3Ti, MP+weights, ML-rates, and ML+rates analyses all resolved Nymphaea as basalmost with higher BS support, a result consistent with recent global angiosperm analyses (Parkinson, Adams, and Palmer 1999Citation ; Qiu et al. 1999Citation ; P. S. Soltis, D. E. Soltis, and Chase 1999Citation ; Bowe, Coat, and dePamphilis 2000Citation ; Chaw et al. 2000Citation ; D. E. Soltis, P. S. Soltis, and Chase 2000Citation ). Both the MP+3Ti and the MP-3Ti analyses placed Juniperus as sister to the gnetophytes, although in both cases BS support for this relationship was low. In contrast, ML both with and without rates placed Pinus as sister to the gnetophytes, with Juniperus sister to that entire clade.

The most striking differences between the MP+3Ti and MP-3Ti analyses involved bryophyte relationships. The single shortest tree recovered from the multigene data set using MP+3Ti (fig. 2 ) showed the liverworts Marchantia and Calypogeia as sister to all other land plants (the LBT), with paraphyletic mosses and the hornworts Megaceros and Phaeoceros sister to the tracheophytes. None of these bryophyte relationships received BS support >50%. Mosses were not supported as monophyletic using any of the nonmitochondrial partitions separately (table 2 ), but only achieved monophyly in the multigene data sets in which rbcL third position transitions were excluded or downweighted (MP+weights). A similar result was obtained for the sister group relationship between mosses and liverworts. One topological difference occurred within the moss clade when using ML-rates as opposed to ML+rates and MP-3Ti. The former method placed Polytrichum as sister to a Takakia/Sphagnum clade (tree not shown), whereas the latter methods placed Takakia as sister to Sphagnum/Polytrichum (fig. 3 ). For the MP+weights analysis, the HBT was obtained using all weights up to 0.93, whereas the LBT was retrieved with weights above this threshold.

In general, BS support for most clades was equivalent or higher for ML+rates than for ML-rates; this was true for the multigene data set and, to a lesser extent, for analyses of individual genes (table 2 ). Such an increase in BS support can be seen for the HBT, where the ML+rates and ML-rates BS values were 62% and 78%, respectively. Three clades in the multigene analysis showed marked differences in BS values between ML+rates and MP: gymnosperms, lycophytes, and HBT (table 2 ). For the HBT, the lower BS values obtained when using ML-rates (62%) or ML+rates (78%) as compared with MP (95%) may be attributed to the influence of the rbcL partition. Of the four partitions, only mtSSU and cpSSU rDNA analyzed separately retrieved the HBT (but with low to moderate BS support). Using ML and MP methods, combining the three SSU rDNA partitions gave BS values of 96%–97% for the HBT. That combined analyses can result in increased resolution and BS support for particular clades that were not well supported following individual analyses has been documented (Kluge 1989Citation ; Mishler et al. 1994Citation ; Soltis et al. 1998Citation ; Soltis, Soltis, and Chase 1999Citation ). Moreover, clades that were not present following analyses of individual genes can appear and gain significant BS support upon combined analyses. Examples here include the euphyllophyte, fern, and gymnosperm clades using rbcL, cpSSU, and nuSSU rDNA separately and in combination.

Decay Analyses and Hypothesis Testing
Decay analysis showed that 12 additional steps were required to collapse the hornworts-basal node (fig. 3 ). Five alternative topologies were tested with the multigene data set using the K-H test with ML and MP (fig. 4AF ). Two frequently cited LBTs depict either hornworts (fig. 4C ) or mosses (fig. 4D ) as sister to tracheophytes, but both were rejected as statistically worse (P < 0.05 level) than the ML tree. Similarly, the paraphyletic moss topology (fig. 4F ) was significantly worse than the MP and ML trees. An alternative HBT (fig. 4B ), in which liverworts and mosses were not sister taxa, added six steps to the most-parsimonious MP tree and could not be rejected relative to the ML tree. Finally, a monophyletic-bryophytes topology, although 14 steps longer than the MP tree, also could not be rejected relative to the ML tree. The K-H and Templeton tests conducted using MP as the criterion showed that all alternative trees (BF) were significantly worse than tree A (P < 0.05).



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 4.—Alternative topologies imposed on the multigene data set: A (Hedderson, Chapman, and Rootes 1996Citation ; Garbary and Renzaglia 1998Citation ; Duff and Nickrent 1999; Nishiyama and Kato 1999Citation ; present work), B (Mishler et al. 1992Citation ), C (Lewis, Mishler, and Vilgalys 1997Citation ), D (Mishler and Churchill 1984Citation ), E (Garbary, Renzaglia, and Duckett 1993Citation ), and F (Malek et al. 1996Citation ). Numbers below trees represent additional steps compared with the shortest MP tree (length = 5,251). Results from K-H ML tests: * = statistically significantly different from best tree (A) at P < 0.05; — = not statistically significantly different. K-H MP and Templeton tests resulted in all alternate topologies being significantly different from tree A at P < 0.05. O = outgroup; H = hornworts; L = liverworts; M = mosses; Sp = Sphagnum; T = tracheophytes

 

    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
Outgroup Selection and Comparison Among the SSU Gene Partitions
This study is the first to sample among all major land plant lineages and across four genes representing the three plant subcellular genomes. Although sampling among tracheophytes was adequate, it is clear that additional bryophyte sequences are needed to place in context the results reported here. With the mtSSU rDNA data partition, the choice of outgroup is critical. In a previous analysis that utilized mtSSU rDNA (Duff and Nickrent 1999Citation ), basal land plant relationships were not unambiguously resolved with MP (although the HBT was retrieved with ML). This result was likely a consequence of using an algal outgroup (Prototheca) whose mtSSU sequence was much more distantly related to land plant sequences than the charophycean algal sequences used here. Pairwise comparisons of Kimura two-parameter distances using mtSSU rDNA gave the following results: Chara to Phaeoceros = 0.114, Chara to Glycine = 0.134, Glycine to Phaeoceros = 0.110, Prototheca to Phaeoceros = 0.342, Prototheca to Glycine = 0.364, and Prototheca to Chara = 0.347. These distances show that Chara is patristically closer to land plants than Prototheca (and hence probably a more appropriate outgroup), and branch lengths from Chara to ingroup taxa were not disproportionately long so as to perturb phylogenetic relationships (Felsenstein 1978Citation ).

Lower homoplasy in mitochondrial genes may derive from generally lower mutation rates as compared with nuclear or chloroplast genes (Wolfe, Li, and Sharp 1987Citation ; Palmer 1990Citation ; Wolfe 1996Citation ) and further highlights their utility in addressing ancient land plant relationships. Chloroplast SSU rDNA had the second lowest amount of homoplasy. Although this gene alone contains insufficient phylogenetic signal to resolve many clades, it did recover the liverwort clade and the HBT, two relationships that were not obtained using the longer and more variable nuSSU rDNA alone (table 2 ). The poor overall performance of nuSSU rDNA in resolving land plant relationships cannot be attributed to insufficient taxon sampling in this study, as similarly poor results were obtained in a recent study that utilized 93 ingroup sequences (Soltis et al. 1999Citation ).

Influence of the rbcL Partition on Resulting Trees
Chloroplast rbcL alone recovered a number of land plant clades but contained the highest level of homoplasy among the four gene partitions (table 2 ). When rbcL was analyzed alone using MP with all positions equally weighted, three trees were obtained that supported the LBT (62% BS support; results not shown). Removal of all third positions resulted in extensive loss of resolution across the entire tree, in agreement with Lewis, Mishler, and Vilgalys (1997)Citation and Källersjö et al. (1998)Citation . Removal of third-position transitions resulted in a tree with no BS support for any of the major land plant clades. These exercises demonstrated that most of the phylogenetic signal in rbcL resides in third codon positions, specifically, third-position transitions. When the rbcL partition was included in the multigene data set and third codon positions were weighted 0.94 or less, the HBT was obtained, thus demonstrating how these competing topologies are sensitive to the degree of influence imparted by third-position differences.

Evolutionary Rates
Qiu and Palmer (1999)Citation list several factors that might potentially confound the problem of determining the basalmost land plant lineage: (1) a wide "evolutionary gap" (both with morphological and molecular characters) between charophyte outgroups and land plants that eliminates synapomorphies; (2) long-branch attraction problems between ingroup and outgroup; (3) ancient and rapid radiation of bryophyte lineages at the base of the land plant tree, leaving little phylogenetic signal; and (4) extinction. As discussed above, the sequence divergence between Chara and a hornwort is equivalent to that between a hornwort and a flowering plant, thus demonstrating that the "evolutionary gap" between ingroup and outgroup is not disproportionate, at least for slowly evolving genes such as mtSSU rDNA. As shown in figures 2 and 3 , both the gnetophyte and the Selaginella clades have increased substitution rates, as indicated by their branch lengths, but in neither case can this be attributed to low taxon sampling. The fast evolutionary rate of gnetophytes has been well documented; however, we do not feel the position of the gnetophyte clade as sister to conifers is artifactual given that this result is in agreement with several recent molecular phylogenetic studies (Bowe, Coat, and dePamphilis 2000Citation ; Chaw et al. 2000Citation ; Sanderson et al. 2000Citation ). Similarly, despite the long branch leading to Selaginella, its position as sister to Isoëtes when ML+rates was employed is fully in line with molecular and morphological analyses. Long-branch effects can be demonstrated with the nuSSU rDNA partition when analyzed alone with MP. Here, Selaginella migrates to the base of the tree, near liverworts. Among bryophytes, there is some evidence of lineage-specific and gene-specific substitution rate heterogeneity. For example, transition rates in complex thalloid liverwort rbcL are about half those seen in leafy and simple thalloid liverworts (Lewis, Mishler, and Vilgalys 1997Citation ), whereas the opposite is true for nuSSU rDNA (Capesius and Bopp 1997Citation ). There is no indication in the present analysis, however, that relationships within and among hornworts, liverworts, and mosses are being influenced by long-branch attraction artifacts.

Mitochondrial Introns
The presence/absence of three mitochondrial introns (Qiu et al. 1998Citation ) has been interpreted as supporting the LBT, as all three introns "are present, with occasional losses, in mosses, hornworts and all major lineages of vascular plants, but are entirely absent from liverworts, green algae and all other eukaryotes." How can this result be reconciled with the findings reported here? If the HBT is correct, these introns must have been acquired in the common ancestor of all land plants and then lost in the common ancestor of liverworts. This scenario is more likely if liverworts are monophyletic, as supported by our analyses (fig. 1 ), than if they are paraphyletic, as suggested by some previous studies that sampled more broadly within liverworts (Capesius and Bopp 1997Citation ; Lewis, Mishler, and Vilgalys 1997Citation ). Postulating a single additional loss of each intron is not unreasonable given that all three introns are already known to have been lost in two different lineages (Isoëtes and Ephedra) and that two of the three introns have been lost many times in land plant evolution (Qiu et al. 1998Citation ; Qiu and Palmer, personal communication). Addition of the three intron characters to the multigene matrix does not change the MP-3Ti tree topology, and BS support for the HBT remains high (92%; data not shown).

Reinterpretation of Morphological Features
Given the strong support for the HBT, reinterpretation of key morphological and anatomical innovations is now required. Several ultrastructural features shared between hornworts and charophycean algae (such as Coleochaete) support the concept that hornworts are the earliest land plant lineage. Unique among land plants, hornwort chloroplasts contain pyrenoids (localized groupings of RUBISCO) that are essentially identical to those in green algae (Vaughn et al. 1992Citation ) and also share with algae channel thylakoids. The HBT thus requires that these chloroplast structures have been lost in other land plants. Conversely, hornworts lack chloroplast granal end membranes and staggered spermatozoid flagella, which presumably evolved subsequent to the divergence of hornworts and other land plants (Garbary and Renzaglia 1998Citation ). If the HBT and the sister group status of liverworts and mosses are correct, the presence of stomates in some hornworts and mosses and all vascular plants may indicate that these structures are not homologous (Renzaglia et al. 2000Citation ). If moss and vascular plant stomates are homologous, then stomates were lost in the ancestor of liverworts. Similarly, a loss in liverworts is invoked for the perine layer in spores, a feature shared by mosses and vascular plants (Garbary and Renzaglia 1998Citation ).

As in most analyses to date, the paraphyly of bryophytes is demonstrated in this multigene analysis. The sister group relationship between mosses and liverworts has received support from data derived from spermatogenesis (Garbary, Renzaglia, and Duckett 1993Citation ), sperm ultrastructure (Maden et al. 1997Citation ), and overall morphology (Garbary and Renzaglia 1998Citation ), as well as molecular characters such as nuclear SSU rDNA (Hedderson, Chapman, and Cox 1998Citation ) and mitochondrial SSU rDNA (Duff and Nickrent 1999 Citation ) and an analysis of five chloroplast genes (Nishiyama and Kato 1999Citation ). The morphological and ultrastructural synapomorphies that support a moss/liverwort clade are (1) the presence of gametophytic slime papillae, (2) superficial archegonia, (3) dimorphic basal bodies on spermatozoa, and (4) apertures associated with the spline on spermatozoa (Garbary and Renzaglia 1998Citation ). Takakia was first described as a liverwort; however, substantial morphological (Renzaglia, McFarland, and Smith 1997Citation ) and molecular (Hedderson, Chapman, and Rootes 1996Citation ; Duff and Nickrent 1999Citation ) evidence now places the genus within the mosses, a result supported by this multigene analysis. Additional sequencing to allow inclusion of putatively basal liverworts such as Haplomitrium and Monoclea in future multigene analyses is required to test the robustness of the moss-liverwort clade.


    Conclusions and Future Work
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
Although both a thalloid and a leafy liverwort were included in this study, increased sampling of this group and other bryophytes is required to represent more adequately the overall taxonomic diversity of bryophytes; work is currently underway to remedy this situation (unpublished data). In addition to ingroup taxon density, inclusion of other green algal outgroup taxa is critical for future land plant phylogenetic studies. As mentioned in Materials and Methods, a mtSSU rDNA sequence for Coleochaete has yet to be obtained. When sequences from three genes from Coleochaete are included in the four-gene data set (with its mtSSU rDNA scored as missing), an MP+3Ti analysis yields the HBT (73% BS support) as opposed to the LBT, which is obtained without these Coleochaete sequences. As the taxonomic breadth of sequences increases, other genes, such as mitochondrial nad5 (Beckert et al. 1999Citation ), cox3 (Malek et al. 1996Citation ), and chloroplast photosystem genes (Nishiyama and Kato 1999Citation ), can be added to the present four-gene data set. A preliminary analysis in which nad5 was added to the four-gene matrix reported here (18 ingroup taxa, half of which were placeholders) gave a tree fully congruent with the one shown in figure 3 , with the HBT receiving 90% BS support (unpublished data). The present study, and several recent ones on angiosperm and gymnosperm evolution, demonstrate the tremendous utility of mitochondrial genes. That these genes often give strong support for many noncontroversial relationships provides a compelling reason to consider seriously those relationships inferred from these genes that are not recovered using other data sets.

As stated in the Introduction, a goal of this study was to determine whether molecular and morphological data better support the LBT or the HBT. It appears that most support for the former hypothesis derives from the chloroplast gene rbcL, specifically, third-codon-position transitions of this gene. In contrast, the HBT is supported by the present study, as well as combined analyses of rbcL and other chloroplast genes (Nishiyama and Kato 1999Citation ), mitochondrial genes (Malek et al. 1996Citation ; Beckert et al. 1999Citation ; Duff and Nickrent 1999Citation ), nuclear SSU rDNA (Hedderson, Chapman, and Cox 1998Citation ), and morphology (Garbary and Renzaglia 1998Citation ). Decay analyses indicate that all the alternative topologies tested using MP result in longer trees, and K-H tests show that these alternatives, including the LBT, are statistically different from the HBT. Thus, we conclude that the HBT hypothesis is most consistent with all of the available evidence (morphological and molecular), is not in conflict with presence/absence data on mitochondrial introns, and therefore provides the best explanation of early land plant evolution. Further refinement and resolution of this tree will provide a firm phylogenetic foundation to help interpret the evolution of plant morphology, ultrastructure, development, and biochemistry.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 
We acknowledge Y.-L. Qiu for several plant DNAs and Karen Renzaglia for insightful discussions and encouragement. This research was supported by grants from the NSF (MCB-9808752 to D.L.N.) and the NIH (GM-35087 to J.D.P.) and by an NIH fellowship (GM-19225 to C.L.P).


    Footnotes
 
Pamela Soltis, Reviewing Editor

1 Present address: Department of Biology, University of Central Florida. Back

2 Present address: Department of Biology, University of Akron. Back

3 Abbreviations: HBT, hornwort basal topology; LBT, liverwort basal topology. Back

4 Keywords: land plant phylogeny embryophyte hornwort liverwort moss rbcL, small-subunit ribosomal DNA Back

5 Address for correspondence and reprints: Daniel L. Nickrent, Department of Plant Biology and Center for Systematic Biology, Southern Illinois University, Carbondale, Illinois 62901-6509. E-mail: nickrent{at}plant.siu.edu Back


    literature cited
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions and Future Work
 Acknowledgements
 literature cited
 

    Beckert, S., S. Steinhauser, H. Muhle, and V. Knoop. 1999. A molecular phylogeny of bryophytes based on nucleotide sequences of the mitochondrial nad5 gene. Plant Syst. Evol. 218:179–192[ISI]

    Bowe, L. M., G. Coat, and C. W. dePamphilis. 2000. Phylogeny of seed plants based on all three plant genomic compartments: extant gymnosperms are monophyletic and Gnetales derived conifers. Proc. Natl. Acad. Sci. USA 97:4092–4097

    Bremer, K. 1994. Branch support and tree stability. Cladistics 10:295–304

    Capesius, I., and M. Bopp. 1997. New classification of liverworts based on molecular and morphological data. Plant Syst. Evol. 207:87–97[ISI]

    Chaw, S.-M., C. L. Parkinson, Y. Cheng, T. M. Vincent, and J. D. Palmer. 2000. Seed plant phylogeny inferred from all three plant genomes: monophyly of extant gymnosperms and origin of Gnetales from conifers. Proc. Natl. Acad. Sci. USA 97:4086–4091

    Duff, R. J., and D. L. Nickrent. 1999. Phylogenetic relationships of land plants using mitochondrial small-subunit rDNA sequences. Am. J. Bot. 86:372–386[Abstract/Free Full Text]

    Eriksson, T. 1998. AutoDecay, version 4.0. Bergius Foundation, Royal Swedish Academy of Sciences, Stockholm

    Felsenstein, J. 1978. Cases in which parsimony or compatibility will be positively misleading. Syst. Zool. 27:401–410[ISI]

    ———. 1984. Distance methods for inferring phylogenies: a justification. Evolution 38:16–24

    ———. 1993. PHYLIP (phylogeny inference package). Version 3.5. Distributed by the author, Department of Genetics, University of Washington, Seattle

    Garbary, D. J., and K. S. Renzaglia. 1998. Bryophyte phylogeny and the evolution of land plants: evidence from development and ultrastructure. Pp. 45–63 in J. W. Bates, N. W. Ashton, and J. G. Duckett, eds. Bryology for the twenty-first century. Maney and the British Bryological Society, Leeds, England

    Garbary, D. J., K. S. Renzaglia, and J. G. Duckett. 1993. The phylogeny of land plants: a cladistic analysis based on male gametogenesis. Plant Syst. Evol. 188:237–269[ISI]

    Gilbert, D. G. 1993. SeqApp. Version 1.9a157. Biocomputing Office, Biology Department, Indiana University, Bloomington

    Goremykin, V., V. Bobrova, J. Pahnke, A. Troitsky, A. Antonov, and W. Martin. 1996. Noncoding sequences from the slowly evolving chloroplast inverted repeat in addition to rbcL data do not support Gnetalean affinities of angiosperms. Mol. Biol. Evol. 13:383–396[Abstract]

    Gray, J. 1993. Major Paleozoic land plant evolutionary bio-events. Palaeogeogr. Palaeoclimatol. Palaeoecol. 104:153–169[ISI]

    Hedderson, T. A., R. Chapman, and C. J. Cox. 1998. Bryophytes and the origins and diversification of land plants: new evidence from molecules. Pp. 65–77 in J. W. Bates, N. W. Ashton, and J. G. Duckett, eds, Bryology for the twenty-first century. Maney and the British Bryological Society, Leeds, England

    Hedderson, T. A., R. L. Chapman, and W. L. Rootes. 1996. Phylogenetic relationships of bryophytes inferred from nuclear-encoded rRNA gene sequences. Plant Syst. Evol. 200:213–224[ISI]

    Hillis, D. M., and J. J. Bull. 1993. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst. Biol. 42:182–191[ISI]

    Hueber, F. M. 1961. Hepaticites devonicus, a new fossil liverwort from the Devonian of New York. Ann. Mo. Bot. Gard. 48:125–132

    Källersjö, M., J. S. Farris, M. W. Chase, B. Bremer, M. F. Fay, C. J. Humphries, G. Petersen, O. Seberg, and K. Bremer. 1998. Simultaneous parsimony jackknife analysis of 2538 rbcL DNA sequences reveals support for major clades of green plants, seed plants and flowering plants. Plant Syst. Evol. 213:259–287[ISI]

    Kenrick, P., and P. R. Crane. 1997. The origin and early diversification of land plants: a cladistic study. Smithsonian Institute Press, Washington, D.C

    Kishino, H., and M. Hasegawa. 1989. Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order of Hominoidea. J. Mol. Evol. 29:170–179[ISI][Medline]

    Kluge, A. G. 1989. A concern for evidence and a phylogenetic hypothesis of relationships among Epicrates (Boidae, Serpentes). Syst. Zool. 38:7–25[ISI]

    Lewis, L. A., B. D. Mishler, and R. Vilgalys. 1997. Phylogenetic relationships of the liverworts (Hepaticae), a basal embryophyte lineage, inferred from nucleotide sequence data of the chloroplast gene rbcL. Mol. Phylogenet. Evol. 7:377–393

    Maden, A. R., D. P. Whittier, D. J. Garbary, and K. S. Renzaglia. 1997. Ultrastructure of the spermatozoid of Lycopodiella lateralis (Lycopodiaceae). Can. J. Bot. 75:1728–1738[ISI]

    Malek, O., K. Lättig, R. Hiesel, A. Brennicke, and V. Knoop. 1996. RNA editing in bryophytes and a molecular phylogeny of land plants. EMBO J. 15:1403–1411[Abstract]

    Manhart, J. R. 1994. Phylogenetic analysis of green plant rbcL sequences. Mol. Phylogenet. Evol. 3:114–127[Medline]

    Mishler, B. D., and S. P. Churchill. 1984. A cladistic approach to the phylogeny of the bryophytes. Brittonia 36:406–424

    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. Phylogenetic relationships of the "green algae" and "bryophytes." Ann. Mo. Bot. Gard. 81:451–483

    Mishler, B. D., P. H. Thrall, J. S. J. Hopple, E. De Luna, and R. Vilgalys. 1992. A molecular approach to the phylogeny of bryophytes: cladistic analysis of chloroplast-encoded 16S and 23S ribosomal RNA genes. Bryologist 95:172–180

    Nemejc, F., and B. Pacltova. 1974. Hepaticae in the Senonian of South Bohemia. Paleobotanist 21:23–26

    Nickrent, D. L. 1994. From field to film: rapid sequencing methods for field collected plant species. BioTechniques 16:470–475

    Nickrent, D. L., R. J. Duff, and D. A. M. Konings. 1997. Structural analyses of plastid-derived 16S rRNAs in holoparasitic angiosperms. Plant Mol. Biol. 34:731–743[ISI][Medline]

    Nickrent, D. L., and D. E. Soltis. 1995. A comparison of angiosperm phylogenies based upon complete 18S rDNA and rbcL sequences. Ann. Mo. Bot. Gard. 82:208–234

    Nickrent, D. L., and E. M. Starr. 1994. High rates of nucleotide substitution in nuclear small-subunit (18S) rDNA from holoparasitic flowering plants. J. Mol. Evol. 39:62–70[ISI][Medline]

    Nishiyama, T., and M. Kato. 1999. Molecular phylogenetic analysis among bryophytes and tracheophytes based on combined data of plastid coded genes and the 18S rRNA gene. Mol. Biol. Evol. 16:1027–1036[Abstract]

    Palmer, J. D. 1990. Contrasting modes and tempos of genome evolution in land plant organelles. Trends Genet. 6:115–120[ISI][Medline]

    Parkinson, C. L., K. L. Adams, and J. D. Palmer. 1999. Multigene analyses identify the three earliest lineages of extant flowering plants. Curr. Biol. 9:1485–1488[ISI][Medline]

    Pryer, K. M., A. R. Smith, and J. E. Skog. 1995. Phylogenetic relationships of extant ferns based on evidence from morphology and rbcL sequences. Am. Fern J. 85:205–282[ISI]

    Qiu, Y., Y. Lee, F. Bernasconi-Quadroni, D. E. Soltis, P. S. Soltis, M. Zanis, E. A. Zimmer, Z. Chen, V. Savolainen, and M. W. Chase. 1999. The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes. Nature 402:404–409

    Qiu, Y.-L., Y. Cho, J. C. Cox, and J. D. Palmer. 1998. The gain of three mitochondrial introns identifies liverworts as the earliest land plants. Nature 394:671–674

    Qiu, Y.-L., and J. D. Palmer. 1999. Phylogeny of early land plants: insights from genes and genomes. Trends Plant Sci. 4:26–30[ISI][Medline]

    Raubeson, L. A., and R. K. Jansen. 1992. Chloroplast DNA evidence on the ancient evolutionary split in vascular land plants. Science 255:1697–1699

    Renzaglia, K. S., and J. G. Duckett. 1991. Towards an understanding of the differences between the blepharoplast of mosses and liverworts, comparisons with hornworts, biflagellated lycopods and charophytes: a numerical analysis. New Phytol. 117:187–208[ISI]

    Renzaglia, K. S., R. J. Duff, D. L. Nickrent, and D. J. Garbary. 2000. Vegetative and reproductive innovations of early land plants: implications for a unified phylogeny. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355:769–793[ISI][Medline]

    Renzaglia, K. S., K. D. McFarland, and D. K. Smith. 1997. Anatomy and ultrastructure of the sporophyte of Takakia ceratophylla (Bryophyta). Am. J. Bot. 84:1337–1350[Abstract]

    Richardson, J. B. 1985. Lower Palaeozoic sporomorphs: their stratigraphical distribution and possible affinities. Philos. Trans. R. Soc. Lond. B Biol. Sci. 309:201–205[ISI]

    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. Mol. Biol. Evol. 17:782–797[Abstract/Free Full Text]

    Soltis, D. E., and P. S. Soltis. 1998. Choosing an approach and an appropriate gene for phylogenetic analysis. Pp. 1–42 in D. E. Soltis, P. S. Soltis, and J. J. Doyle, eds. Molecular systematics of plants II. DNA sequencing. Kluwer Academic Publishers, Boston

    Soltis, D. E., P. S. Soltis, M. W. Chase et al. (13 co-authors). 2000. Angiosperm phylogeny inferred from 18S rDNA, rbcL, and atpB sequences. Bot. J. Linn. Soc. 133:381–461[ISI]

    Soltis, D. E., P. S. Soltis, M. E. Mort, M. W. Chase, V. Savolainen, S. B. Hoot, and C. M. Morton. 1998. Inferring complex phylogenies using parsimony: an empirical approach using three large DNA data sets for angiosperms. Syst. Biol. 47:32–42[ISI][Medline]

    Soltis, P. S., D. E. Soltis, and M. W. Chase. 1999. Angiosperm phylogeny inferred from multiple genes as a tool for comparative biology. Nature 402:402–407

    Soltis, P. S., D. E. Soltis, P. G. Wolf, D. L. Nickrent, S.-M. Chaw, and R. L. Chapman. 1999. The phylogeny of land plants inferred from 18S rDNA sequences: pushing the limits of rDNA signal? Mol. Biol. Evol. 16:1774–1784[Abstract/Free Full Text]

    Strimmer, K., and A. von Haeseler. 1996. Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies. Mol. Biol. Evol. 13:964–969[Free Full Text]

    Swofford, D. L. 1998. PAUP*: phylogenetic analysis using parsimony. Version 4.0 b2a. Sinauer, Sunderland, Mass

    Taylor, W. A. 1995. Spores in earliest land plants. Nature 373:393–392

    Templeton, A. R. 1983. Phylogenetic inference from restriction endonuclease cleavage site maps with particular reference to the evolution of humans and the apes. Evolution 37:221–244

    Vaughn, K. C., R. Ligrone, H. A. Owens, J. Hasegawa, E. O. Campbell, K. S. Renzaglia, and J. Monge-Najera. 1992. The Anthocerote chloroplast: a review. New Phytol. 120:169–190[ISI]

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

    Wolfe, K., W.-H. Li, and P. Sharp. 1987. Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs. Proc. Natl. Acad. Sci. USA 84:9054–9058

    Wolfe, K. H. 1996. Molecular evolution of plants: more genomes, less generalities. Pp. 45–57 in B. Andersson, A. H. Salter, and J. Barber, eds. Molecular genetics of photosynthesis. IRL Press, Oxford, England

Accepted for publication August 23, 2000.