Phylogenetics of flowering plants based on combined analysis of plastid atpB and rbcL gene sequences

Inference of higher-order relationships in the cycads from a large chloroplast data set

We investigated higher-order relationships in the cycads, an ancient group of seed-bearing plants, by examining a large portion of the chloroplast genome from seven species chosen to exemplify our current understanding of taxonomic diversity in the order. The regions considered span approximately 13.5 kb of unaligned data per taxon, and comprise a diverse range of coding sequences, introns and intergenic spacers dispersed throughout the plastid genome. Our results provide substantial support for most of the inferred backbone of cycad phylogeny, and weak evidence that the sister-group of the cycads among living seed plants is Ginkgo biloba. Cycas (representing Cycadaceae) is the sister-group of the remaining cycads; Dioon is part of the next most basal split. Two of the three commonly recognized families of cycads (Zamiaceae and Stangeriaceae) are not monophyletic; Stangeria is embedded within Zamiaceae, close to Zamia and Ceratozamia, and not closely allied to the other genus of Stangeriaceae, Bowenia. In contrast to the other seed plants, cycad chloroplast genomes share two features with Ginkgo: a reduced rate of evolution and an elevated transition:transversion ratio. We demonstrate that the latter aspect of their molecular evolution is unlikely to have affected inference of cycad relationships in the context of seed-plant wide analyses.

Evolution of carnivory in Lentibulariaceae and the Lamiales

As a basis for analysing the evolution of the carnivorous syndrome in Lentibulariaceae (Lamiales), phylogenetic reconstructions were conducted based on coding and non-coding chloroplast DNA (matK gene and flanking trnK intron sequences, totalling about 2.4 kb). A dense taxon sampling including all other major lineages of Lamiales was needed since the closest relatives of Lentibulariaceae and the position of "proto-carnivores" were unknown. Tree inference using maximum parsimony, maximum likelihood, and Bayesian approaches resulted in fully congruent topologies within Lentibulariaceae, whereas relationships among the different lineages of Lamiales were only congruent between likelihood and Bayesian optimizations. Lentibulariaceae and their three genera (Pinguicula, Genlisea, and Utricularia) are monophyletic, with Pinguicula being sister to a Genlisea-Utricularia clade. Likelihood and Bayesian trees converge on Bignoniaceae as sister to Lentibulariaceae, albeit lacking good support. The "proto-carnivores" (Byblidaceae, Martyniaceae) are found in different positions among other Lamiales but not as sister to the carnivorous Lentibulariaceae, which is also supported by Khishino-Hasegawa tests. This implies that carnivory and its preliminary stages ("proto-carnivores") independently evolved more than once among Lamiales. Ancestral states of structural characters connected to the carnivorous syndrome are reconstructed using the molecular tree, and a hypothesis on the evolutionary pathway of the carnivorous syndrome in Lentibulariaceae is presented. Extreme DNA mutational rates found in Utricularia and Genlisea are shown to correspond to their unusual nutritional specialization, thereby hinting at a marked degree of carnivory in these two genera.

The phylogenetic history of Selaginellaceae based on DNA sequences from the plastid and nucleus: extreme substitution rates and rate heterogeneity

Molecular phylogenetic research on Selaginellaceae has focused on the plastid gene rbcL, which in this family has unusually high substitution rates. Here we develop a molecular data set from the nuclear 26S ribosomal DNA gene with the aim of evaluating and extending the results of previous phylogenetic research. The 26S rDNA and the rbcL regions were sequenced for a sample of 23 species, which represent the main elements of species diversity in the family. The data were analysed independently and in combination using both maximum parsimony and Bayesian inference. Although several between genome differences were found, the general pattern of relationships uncovered by all analyses was very similar. Results corroborate the previous study supporting new groupings not previously recognised on morphological grounds. Substitution rates in the 26S rDNA were also found to be high (26% informative) for the region analysed, but lower than for rbcL (37% informative). These data indicate that high substitution rates might be widespread in all three genomes (i.e., plastid, mitochondrion, and nucleus).

Amborella not a "basal angiosperm"? Not so fast

The sequence of the plastid genome of Amborella trichopoda, the putative sister to all other extant angiosperms, was recently reported (Molecular Biology and Evolution 20: 1499-1505). Goremykin et al. used sequence data for 61 plastid genes from Amborella and 12 other embryophytes in phylogenetic analyses and concluded that Amborella is not the sister to the remaining flowering plants; the monocots instead occupy this position. The authors attributed their results, which differ substantially from all recent phylogenetic analyses of angiosperms, to the increased character sampling (30 017 nucleotides in their aligned matrix) in their analysis relative to published studies that included fewer genes but more taxa. We hypothesized that the difference in topology is not due to limited character sampling in previous studies but to limited taxon sampling in the analysis by Goremykin et al. To test this, we conducted a series of phylogenetic analyses using a three-gene, 12 (or more)-taxon data set to evaluate the topological effects of (i) including three vs. 61 genes for (nearly) the same set of taxa, (ii) analyzing different codon positions, (iii) substituting representatives of other basal lineages for Amborella, (iv) replacing the grasses used to represent the monocots with other monocots, selected either for their phylogenetic position or randomly, and (v) adding other basal taxa-Nymphaea, Austrobaileya, magnoliids, and monocots-to the 12-taxon data set. Our results demonstrate that the "monocots basal" topology obtained by Goremykin et al. is not due to increased character sampling of the plastid genome; their topology was obtained using only two plastid genes or two plastid genes and one nuclear gene. This topology was also retained when either Nymphaea or Austrobaileya was substituted for Amborella, demonstrating that any of the three basal lineages will attach to Calycanthus for lack of any other close branch. Furthermore, the "monocots basal" topology is not robust to changes in sampling of monocots. Simply adding Oncidium, for example, places Amborella sister to the other angiosperms. Thus, limited taxon sampling, focusing on organisms with complete genome sequences, can lead to artifactual results.

The sequence of the largest subunit of RNA polymerase II is a useful marker for inferring seed plant phylogeny

We used RT-PCR to sequence approximately 3 kb of the gene coding for the largest subunit of RNA polymerase II (rpb1) from nine land plants. Our results show that plant rpb1 genes all have a similar GC-content and that their amino acid sequences evolve at a similar rate in most species we examined, except for the Arabidopsis thaliana and rice sequences which evolve faster. This gene also exists as a single copy in most species and contains enough phylogenetically informative sites to resolve the evolutionary relationships among seed plants. Protein maximum parsimony, as well as neighbor-joining and maximum likelihood analyses of DNA and protein sequences, all generated identical tree topologies with similar strong support values at each node. The angiosperms are a clade comprising Amborella as a sister group to all other angiosperms, followed by Nymphaea, Magnolia, Arabidopsis, and a monocot clade containing maize and rice. The gymnosperms also form a monophyletic clade with Welwitschia and pine grouped together and sister to a Cycas and Zamia clade. These findings concur with recent studies that refute the Anthophyte Hypothesis and place Amborella at the base of the angiosperm tree. These rpb1 sequences also give a more consistent picture of seed plant relationships than similar analyses performed on data sets made of 18S rDNA, atpB, and rbcL sequences from the same species. These sequences therefore show great promise to help further resolve the phylogenetic relationships of seed plants.

Darwin's abominable mystery: Insights from a supertree of the angiosperms

Angiosperms are among the major terrestrial radiations of life and a model group for studying patterns and processes of diversification. As a tool for future comparative studies, we compiled a supertree of angiosperm families from published phylogenetic studies. Sequence data from the plastid rbcL gene were used to estimate relative timing of branching events, calibrated by using robust fossil dates. The frequency of shifts in diversification rate is largely constant among time windows but with an apparent increase in diversification rates within the more recent time frames. Analyses of species numbers among families revealed that diversification rate is a labile attribute of lineages at all levels of the tree. An examination of the top 10 major shifts in diversification rates indicates they cannot easily be attributed to the action of a few key innovations but instead are consistent with a more complex process of diversification, reflecting the interactive effects of biological traits and the environment.

Elatinaceae are sister to Malpighiaceae; Peridiscaceae belong to Saxifragales

Phylogenetic data from plastid (ndhF and rbcL) and nuclear (PHYC) genes indicate that, within the order Malpighiales, Elatinaceae are strongly supported as sister to Malpighiaceae. There are several putative morphological synapomorphies for this clade; most notably, they both have a base chromosome number of X = 6 (or some multiple of three or six), opposite or whorled leaves with stipules, unicellular hairs (also uniseriate in some Elatinaceae), multicellular glands on the leaves, and resin (Elatinacae) or latex (Malpighiaceae). Further study is needed to determine if these features are synapomorphic within the order. Malpighiaceae have previously been inferred as sister to Peridiscaceae based on rbcL sequence data, but the rbcL sequence of Whittonia is a chimera of two sequences, neither of which appears to be Whittonia. Our data from plastid (atpB, rbcL) and nuclear (18S rDNA) genes instead place Peridiscaeace as a member of the Saxifragales.

Mitochondrial DNA sequences reveal the photosynthetic relatives of Rafflesia, the world's largest flower

All parasites are thought to have evolved from free-living ancestors. However, the ancestral conditions facilitating the shift to parasitism are unclear, particularly in plants because the phylogenetic position of many parasites is unknown. This is especially true for Rafflesia, an endophytic holoparasite that produces the largest flowers in the world and has defied confident phylogenetic placement since its discovery >180 years ago. Here we present results of a phylogenetic analysis of 95 species of seed plants designed to infer the position of Rafflesia in an evolutionary context using the mitochondrial gene matR (1,806 aligned base pairs). Overall, the estimated phylogenetic tree is highly congruent with independent analyses and provides a strongly supported placement of Rafflesia with the order Malpighiales, which includes poinsettias, violets, and passionflowers. Furthermore, the phylogenetic placement of Mitrastema, another enigmatic, holoparasitic angiosperm with the order Ericales (which includes blueberries and persimmons), was obtained with these data. Although traditionally classified together, Rafflesia and Mitrastema are only distantly related, implying that their endoparasitic habits result from convergent evolution. Our results indicate that the previous significant difficulties associated with phylogenetic placement of holoparasitic plants may be overcome by using mitochondrial DNA so that a broader understanding of the origins and evolution of parasitism may emerge.

18S gene trees are positively misleading for monocot/dicot phylogenetics

Monocots are consistently paraphyletic in 18S gene trees in all studies to date. This anomaly is generally expressed in the phylogenetic associations of two lineages, that of Acoraceae, which is excluded from the monocots, and Ceratophyllaceae, which sometimes clusters within the monocots. Six explanations for these unexpected results are proposed: (1) erroneous published sequences, (2) actual paraphyly of monocots with dicots, (3) insufficient taxon density among relevant taxa; (4) long-branch attraction between selected taxa, (5) an anomalous 18S paralogue in the basal species of monocots, and (6) differential lineage sorting or other molecular evolutionary events. Tests of each of these hypotheses were conducted. For the first five hypotheses, test results refuted the hypothesis. Support for the sixth hypothesis, differential lineage sorting, is the strongest. Since lineage sorting seems to be restricted to a few species, identification and removal of these species prior to performing separate or combined phylogenetic analysis of DNA data incorporating sequences from 18S may be advisable.

Angiosperm phylogeny based on matK sequence information

Plastid matK gene sequences for 374 genera representing all angiosperm orders and 12 genera of gymnosperms were analyzed using parsimony (MP) and Bayesian inference (BI) approaches. Traditionally, slowly evolving genomic regions have been preferred for deep-level phylogenetic inference in angiosperms. The matK gene evolves approximately three times faster than the widely used plastid genes rbcL and atpB. The MP and BI trees are highly congruent. The robustness of the strict consensus tree supercedes all individual gene analyses and is comparable only to multigene-based phylogenies. Of the 385 nodes resolved, 79% are supported by high jackknife values, averaging 88%. Amborella is sister to the remaining angiosperms, followed by a grade of Nymphaeaceae and Austrobaileyales. Bayesian inference resolves Amborella + Nymphaeaceae as sister to the rest, but with weak (0.42) posterior probability. The MP analysis shows a trichotomy sister to the Austrobaileyales representing eumagnoliids, monocots + Chloranthales, and Ceratophyllum + eudicots. The matK gene produces the highest internal support yet for basal eudicots and, within core eudicots, resolves a crown group comprising Berberidopsidaceae/Aextoxicaceae, Santalales, and Caryophyllales + asterids. Moreover, matK sequences provide good resolution within many angiosperm orders. Combined analyses of matK and other rapidly evolving DNA regions with available multigene data sets have strong potential to enhance resolution and internal support in deep level angiosperm phylogenetics and provide additional insights into angiosperm evolution.

A decade of progress in plant molecular phylogenetics

Over the past decade, botanists have produced several thousand phylogenetic analyses based on molecular data, with particular emphasis on sequencing rbcL, the plastid gene encoding the large subunit of Rubisco (ribulose bisphosphate carboxylase). Because phylogenetic trees retrieved from the three plant genomes (plastid, nuclear and mitochondrial) have been highly congruent, the "Angiosperm Phylogeny Group" has used these DNA-based phylogenetic trees to reclassify all families of flowering plants. However, in addition to taxonomy, these major phylogenetic efforts have also helped to define strategies to reconstruct the "tree of life", and have revealed the size of the ancestral plant genome, uncovered potential candidates for the ancestral flower, identified molecular living fossils, and linked the rate of neutral substitutions with species diversity. With an increased interest in DNA sequencing programmes in non-model organisms, the next decade will hopefully see these phylogenetic findings integrated into new genetic syntheses, from genomes to taxa.

Floral development and morphology of Vochysiaceae. I. The structure of the gynoecium

Vochysiaceae are divided into two tribes on the basis of ovary structure (superior trilocular or inferior unilocular). The superior trilocular ovary has been considered basal in the family, and the term "pseudomonomerous" was used to indicate the presumed evolutionary derivation of the unilocular condition from the trilocular. However, recent evidence that Vochysiaceae are Myrtalean suggests that the superior ovary may be secondarily derived. In addition, published figures cast doubt on the interpretation of the putatively unilocular ovaries. To understand these features, floral ontogeny and anatomy were examined using scanning electron microscopy and serial sectioning. In all taxa examined, the ovary develops in an epigynous fashion, on a concave floral apex, supporting the hypothesis that the superior ovary is secondarily derived. Subsequent to initiation of the ovary, differential growth results in ovaries that are superior, inferior, or partly inferior in different genera. Sections of floral buds of the two unilocular genera, Erisma and Erismadelphus, show aborted locules in the latter but not in the former. The application of the term "pseudomonomerous" to both genera obscures this significant difference. The position of the placenta in the truly unilocular genus varies among species, suggesting a character transformation series from multilocular through intermediates to truly unilocular.

Floral development and morphology of Vochysiaceae. II. The position of the single fertile stamen

A conspicuous feature of Vochysiaceae flowers is their single fertile stamen. In some genera, the stamen is in front of a petal, whereas in others it is in front of a sepal. This difference has been attributed to two independent reductions, which implies the stamen is not homologous across the family. The observation that genera with an antesepalous stamen have only one petal, whereas those with an antepetalous stamen have three or five petals, led us to the hypothesis that in all genera the stamen arises in an antepetalous position, but that it is displaced during development in single-petaled taxa. We examined developing buds of five genera using scanning electron microscopy and serial sectioning and conclude that the stamen in all genera is fundamentally antepetalous. The stamen is not displaced. The petal, however, appears to be displaced in some genera. Further, the position of the fertile stamen in Erisma has been misinterpreted. We discuss the evolution of the androecium in this family in a phylogenetic context and consider the significance of symmetry and of loss vs. suppression in the development of Vochysiaceae flowers.

Duplication and Diversification in the APETALA1/FRUITFULL Floral Homeotic Gene Lineage: Implications for the Evolution of Floral Development

Phylogenetic analyses of angiosperm MADS-box genes suggest that this gene family has undergone multiple duplication events followed by sequence divergence. To determine when such events have taken place and to understand the relationships of particular MADS-box gene lineages, we have identified APETALA1/FRUITFULL-like MADS-box genes from a variety of angiosperm species. Our phylogenetic analyses show two gene clades within the core eudicots, euAP1 (including Arabidopsis APETALA1 and Antirrhinum SQUAMOSA) and euFUL (including Arabidopsis FRUITFULL). Non-core eudicot species have only sequences similar to euFUL genes (FUL-like). The predicted protein products of euFUL and FUL-like genes share a conserved C-terminal motif. In contrast, predicted products of members of the euAP1 gene clade contain a different C terminus that includes an acidic transcription activation domain and a farnesylation signal. Sequence analyses indicate that the euAP1 amino acid motifs may have arisen via a translational frameshift from the euFUL/FUL-like motif. The euAP1 gene clade includes key regulators of floral development that have been implicated in the specification of perianth identity. However, the presence of euAP1 genes only in core eudicots suggests that there may have been changes in mechanisms of floral development that are correlated with the fixation of floral structure seen in this clade.

Phylogenetic analyses of Cornales based on 26S rRNA and combined 26S rDNA-MATK-RBCL sequence data

Nuclear 26S rDNA sequences were used to corroborate and test previously published matK-rbcL-based hypotheses of phylogenetic relationships in Cornales. Sequences were generated for 53 taxa including Alangium, Camptotheca, Cornus, Curtisia, Davidia, Diplopanax, Mastixia, Nyssa, and four families: Grubbiaceae, Hydrangeaceae, Hydrostachyaceae, and Loasaceae. Fifteen taxa from asterids were used as outgroups. The 26S rDNA sequences were initially analyzed separately and then combined with matK-rbcL sequences, using both parsimony and maximum likelihood methods. Eight strongly supported major clades were identified within Cornales by all analyses: Cornus, Alangium, nyssoids (Nyssa, Davidia, and Camptotheca), mastixioids (Mastixia and Diplopanax), Hydrangeaceae, Loasaceae, Grubbia-Curtisia, and Hydrostachys. However, relationships among the major lineages are not strongly supported in either 26S rDNA or combined 26S rDNA-matK-rbcL topologies, except for the sister relationships between Cornus and Alangium and between nyssoids and mastixioids in the tree from combined data. Discrepancies in relationships among major lineages, especially the placement of the long-branched Hydrostachys, were found between parsimony and maximum likelihood trees in all analyses. Incongruence between the 26S rDNA and matK-rbcL data sets was suggested, where Hydrangeaceae was found to be largely responsible for the incongruence. The long branch of Hydrostachys revealed in previous analyses was reduced significantly with more sampling. Maximum likelihood analysis of combined 26S rDNA-matK-rbcL sequences suggested that Hydrostachys might be sister to the remainder of Cornales, that Cornus-Alangium are sisters, that nyssoids-mastixioids are sisters, and that Hydrangeaceae-Loasaceae are sisters, consistent with previous analyses of matK-rbcL sequence data.

The challenge of constructing large phylogenetic trees

The amount of sequence data available to reconstruct the evolutionary history of genes and species has increased 20-fold in the past decade. Consequently the size of phylogenetic analyses has grown as well, and phylogenetic methods, algorithms and their implementations have struggled to keep pace. Computational and other challenges raised by this burgeoning database emerge at several stages of analysis, from the optimal assembly of large data matrices from sequence databases, to the efficient construction of trees from these large matrices and the piece-wise assembly of 'supertrees' from those trees in turn. A final challenge is posed by the difficulty of visualizing and making inferences from trees that might soon routinely contain thousands of species.

Gunnerales are sister to other core eudicots: implications for the evolution of pentamery

Phylogenetic relationships among many lineages of angiosperms have been clarified via the analysis of large molecular data sets. However, with a data set of three genes (18S rDNA, rbcL, and atpB), relationships among lineages of core eudicots (Berberidopsidales, Caryophyllales, Gunnerales, Santalales, Saxifragales, asterids, rosids) remain essentially unresolved. We added 26S rDNA sequences to a three-gene matrix for 201 eudicots (8430 base pair aligned nucleotides per taxon). Parsimony analyses provided moderate (84%) jackknife support for Gunnerales, which comprise the two enigmatic families Gunneraceae and Myrothamnaceae, as sister to all other core eudicots. This position of Gunnerales has important implications for floral evolution. A dimerous or trimerous perianth is frequently encountered in early-diverging eudicots (e.g., Buxaceae, Proteales, Ranunculales, Trochodendraceae), whereas in core eudicots, pentamery predominates. Significantly, dimery is found in Gunneraceae and perhaps Myrothamnaceae (the merosity of the latter has also been interpreted as labile). Parsimony reconstructions of perianth merosity demonstrate lability among early-diverging eudicots and further indicate that a dimerous perianth could be the immediate precursor to the pentamerous condition characteristic of core eudicots. Thus, the developmental canalization that yielded the pentamerous condition of core eudicots occurred after the node leading to Gunnerales.

Relationships among seed plants inferred from highly conserved genes: sorting conflicting phylogenetic signals among ancient lineages

Phylogenetic studies based on different types and treatment of data provide substantially conflicting hypotheses of relationships among seed plants. We conducted phylogenetic analyses of sequences of two highly conserved chloroplast genes, psaA and psbB, for a comprehensive taxonomic sample of seed plants and land plants. Parsimony analyses of two different codon position partitions resulted in well-supported, but significantly conflicting, phylogenetic trees. First and second codon positions place angiosperms and gymnosperms as sister clades and Gnetales as sister to Pinaceae. Third positions place Gnetales as sister to all other seed plants. Maximum likelihood trees for the two partitions are also in conflict. Relationships among the main seed plant clades according to first and second positions are similar to those found in parsimony analysis for the same data, but the third position maximum likelihood tree is substantially different from the corresponding parsimony tree, although it agrees partially with the first and second position trees in placing Gnetales as the sister group of Pinaceae. Our results document high rate heterogeneity among lineages, which, together with the greater average rate of substitution for third positions, may reduce phylogenetic signal due to long-branch attraction in parsimony reconstructions. Whereas resolution of relationships among major seed plant clades remains pending, this study provides increased support for relationships within major seed plant clades.

Phylogeny of seed plants based on evidence from eight genes

Relationships among the five groups of extant seed plants (cycads, Ginkgo, conifers, Gnetales, and angiosperms) remain uncertain. To explore relationships among groups of extant seed plants further and to attempt to explain the conflict among molecular data sets, we assembled a data set of four plastid (cpDNA) genes (rbcL, atpB, psaA, and psbB), three mitochondrial (mtDNA) genes (mtSSU, coxI, and atpA), and one nuclear gene (18S rDNA) for 19 exemplars representing the five groups of living seed plants. Analyses of the combined eight-gene data set (15 772 base pairs/taxon) with maximum parsimony (MP), maximum likelihood (ML), and Bayesian approaches reveal a gymnosperm clade that is sister to angiosperms. Within the gymnosperms, a conifer clade includes Gnetales as sister to Pinaceae. Cycads and Ginkgo are either successive sisters to this conifer clade (including Gnetales) or a clade that is sister to conifers and Gnetales. All analyses of the mtDNA partition and ML analyses of the nuclear partition yield very similar topologies. However, MP analyses of the combined cpDNA genes place Gnetales as sister to all other seed plants with strong bootstrap support, whereas ML and Bayesian analyses of the cpDNA data set place Gnetales as sister to Pinaceae. Maximum parsimony and ML analyses of first and second codon positions of the cpDNA partiation also place Gnetales as sister to Pinaceae. In contrast, MP analyses of third codon positions place Gnetales as sister to other seed plants, although ML analyses of third codon positions place Gnetales with Pinaceae. Thus, most of the discrepancies in seed plant topologies involve third codon positions of cpDNA genes. The likelihood ratio (LR) and Shimodaira-Hasegasa (SH) tests were applied to the cpDNA data. The preferred topology based on the LR test is that Gnetales are sister to Pseudotsuga. The SH test based on first and second codon and all three codon positions indicated that there is no significant difference between the best topology (Gnetales sister to Pseudotsuga) and Gnetales sister to a conifer clade. However, there is a significant difference between the best topology and topologies in which Gnetales are sister to the rest of the seed plants or Gnetales sister to angiosperms.

The origin of the apple subfamily (Maloideae; Rosaceae) is clarified by DNA sequence data from duplicated GBSSI genes

For 70 yr the leading hypothesis for the origin of the Maloideae has involved wide hybridization between ancestors of two other subfamilies. The basis of this hypothesis is that Maloideae have a base chromosome number of 17, whereas other Rosaceae are mostly x = 7, 8, or 9. To investigate this hypothesis we cloned and sequenced approximately 1.8 kilobases from the 5' portion of granule-bound starch synthase (GBSSI, or waxy) genes for 89 clones from 32 Rosaceae genera. Previous studies demonstrate the presence of two copies in all Rosaceae (GBSSI-1 and GBSSI-2) and four in Maloideae (GBSSI-1A, GBSSI-1B, GBSSI-2A, and GBSSI-2B). Parsimony and maximum likelihood analyses nest Gillenia, a genus of the southeastern United States with a base chromosome number of 9, within either Maloideae GBSSI-1 or GBSSI-2. Monophyly of Maloideae plus Gillenia is well supported by bootstrap values, loss of the sixth intron in all GBSSI-1 sequences, intron alignability between genera, and numerous nonmolecular characters. Our results falsify the wide-hybridization hypothesis and are consistent with a polyploid origin involving only members of a lineage that contained the ancestors of Gillenia. Under this hypothesis, the subfamily originated in North America, and the high Maloideae chromosome number arose via aneuploidy from x = 18.

Snake phylogeny: evidence from nuclear and mitochondrial genes

We constructed phylogenies of snakes from the c-mos and cytochrome b genes using conventional phylogenetic methods as well as the relatively new method of Bayesian inference. For all methods, there was excellent congruence between the c-mos and cytochrome b genes, implying a high level of support for the shared clades. Our results agree with previous studies in two important respects: first, that the scolecophidians and alethinophidians are monophyletic sister clades; and second, that the Colubroidea is a monophyletic group with the Acrochordidae as its sister clade. However, our results differ from previous studies in the finding that Loxocemus and Xenopeltis cluster with pythons. An additional noteworthy result from our data is that the genera Exiliboa and Ungaliophis, often placed with Tropidophis (and Trachyboa, not included in the present study) in the Tropidophiidae, are in reality boids.

Phylogenetics of asterids based on 3 coding and 3 non-coding chloroplast DNA markers and the utility of non-coding DNA at higher taxonomic levels

Asterids comprise 1/4-1/3 of all flowering plants and are classified in 10 orders and >100 families. The phylogeny of asterids is here explored with jackknife parsimony analysis of chloroplast DNA from 132 genera representing 103 families and all higher groups of asterids. Six different markers were used, three of the markers represent protein coding genes, rbcL, ndhF, and matK, and three other represent non-coding DNA; a region including trnL exons and the intron and intergenic spacers between trnT (UGU) to trnF (GAA); another region including trnV exons and intron, trnM and intergenic spacers between trnV (UAC) and atpE, and the rps16 intron. The three non-coding markers proved almost equally useful as the three coding genes in phylogenetic reconstruction at the high level of orders and families in asterids, and in relation to the number of aligned positions the non-coding markers were even more effective. Basal interrelationships among Cornales, Ericales, lamiids (new name replacing euasterids I), and campanulids (new name replacing euasterids II) are resolved with strong support. Family interrelationships are fully or almost fully resolved with medium to strong support in Cornales, Garryales, Gentianales, Solanales, Aquifoliales, Apiales, and Dipsacales. Within the three large orders Ericales, Lamiales, and Asterales, family interrelationships remain partly unclear. The analysis has contributed to reclassification of several families, e.g., Tetrameristaceae, Ebenaceae, Styracaceae, Montiniaceae, Orobanchaceae, and Scrophulariaceae (by inclusion of Pellicieraceae, Lissocarpaceae, Halesiaceae, Kaliphoraceae, Cyclocheilaceae, and Myoporaceae+Buddlejaceae, respectively), and to the placement of families that were unplaced in the APG-system, e.g., Sladeniaceae, Pentaphylacaceae, Plocospermataceae, Cardiopteridaceae, and Adoxaceae (in Ericales, Ericales, Lamiales, Aquifoliales, and Dipsacales, respectively), and Paracryphiaceae among campanulids. Several families of euasterids remain unclassified to order.

Relationships within Cornales and circumscription of Cornaceae-matK and rbcL sequence data and effects of outgroups and long branches

Phylogenetic relationships in Cornales were assessed using sequences rbcL and matK. Various combinations of outgroups were assessed for their suitability and the effects of long branches and outgroups on tree topology were examined using RASA 2.4 prior to conducting phylogenetic analyses. RASA identified several potentially problematic taxa having long branches in individual data sets that may have obscured phylogenetic signal, but when data sets were combined RASA no longer detected long branch problems. t(RASA) provides a more conservative measurement for phylogenetic signal than the PTP and skewness tests. The separate matK and rbcL sequence data sets were measured as not containing phylogenetic signal by RASA, but PTP and skewness tests suggested the reverse [corrected]. Nonetheless, the matK and rbcL sequence data sets suggested relationships within Cornales largely congruent with those suggested by the combined matK-rbcL sequence data set that contains significant phylogenetic signal as measured by t(RASA), PTP, and skewness tests. Our analyses also showed that a taxon having a long branch on the tree may not be identified as a "long-branched" taxon by RASA. The long branches identified by RASA had little effect on the arrangement of other taxa in the tree, but the placements of the long-branched taxa themselves were often problematic. Removing the long-branched taxa from analyses generally increased bootstrap support, often substantially. Use of non-optimal outgroups (as identified by RASA) decreased phylogenetic resolution in parsimony analyses and suggested different relationships in maximum likelihood analyses, although usually weakly supported clades (less than 50% support) were impacted. Our results do not recommend using t(RASA) as a sole criterion to discard data or taxa in phylogenetic analyses, but t(RASA) and the taxon variance ratio obtained from RASA may be useful as a guide for improved phylogenetic analyses. Results of parsimony and ML analyses of the sequence data using optimal outgroups suggested by RASA revealed four major clades within Cornales: (1) Curtisia-Grubbia, (2) Cornus-Alangium, (3) Nyssa-Camptotheca-Davidia-Mastixia-Diplopanax, and (4) Hydrangeaceae-Loasaceae, with clades (2) and (3) forming a monophyletic group sister to clade (4) and clade (1) sister to the remainder of Cornales. However, there was not strong bootstrap support for relationships among the major clades. The placement of Hydrostachys could not be reliably determined, although most analyses place the genus within Hydrangeaceae; ML analyses, for example, placed the genus as the sister of Hydrangeeae. Our results supported a Cornales including the systematically problematic Hydrostachys, a Cornaceae consisting of Cornus and Alangium, a Nyssaceae consisting of Nyssa and Camptotheca, a monogeneric Davidiaceae, a Mastixiaceae consisting of Mastixia and Diplopanax, and an expanded Grubbiaceae consisting of Grubbia and Curtisia, and two larger families, Hydrangeaceae and Loasaceae.

One hundred and seventeen clades of euagarics

This study provides a first broad systematic treatment of the euagarics as they have recently emerged in phylogenetic systematics. The sample consists of 877 homobasidiomycete taxa and includes approximately one tenth (ca. 700 species) of the known number of species of gilled mushrooms that were traditionally classified in the order Agaricales. About 1000 nucleotide sequences at the 5(') end of the nuclear large ribosomal subunit gene (nLSU) were produced for each taxon. Phylogenetic analyses of nucleotide sequence data employed unequally weighted parsimony and bootstrap methods. Clades revealed by the analyses support the recognition of eight major groups of homobasidiomycetes that cut across traditional lines of classification, in agreement with other recent phylogenetic studies. Gilled fungi comprise the majority of species in the euagarics clade. However, the recognition of a monophyletic euagarics results in the exclusion from the clade of several groups of gilled fungi that have been traditionally classified in the Agaricales and necessitates the inclusion of several clavaroid, poroid, secotioid, gasteroid, and reduced forms that were traditionally classified in other basidiomycete orders. A total of 117 monophyletic groups (clades) of euagarics can be recognized on the basis on nLSU phylogeny. Though many clades correspond to traditional taxonomic groups, many do not. Newly discovered phylogenetic affinities include for instance relationships of the true puffballs (Lycoperdales) with Agaricaceae, of Panellus and the poroid fungi Dictyopanus and Favolaschia with Mycena, and of the reduced fungus Caripia with Gymnopus. Several clades are best supported by ecological, biochemical, or trophic habits rather than by morphological similarities.

Bayesian inference of phylogeny and its impact on evolutionary biology

As a discipline, phylogenetics is becoming transformed by a flood of molecular data. These data allow broad questions to be asked about the history of life, but also present difficult statistical and computational problems. Bayesian inference of phylogeny brings a new perspective to a number of outstanding issues in evolutionary biology, including the analysis of large phylogenetic trees and complex evolutionary models and the detection of the footprint of natural selection in DNA sequences.

Phylogenetic relationships of Theaceae inferred from chloroplast DNA sequence data

Tribal and generic relationships within Theaceae were investigated using cladistic analyses of chloroplast-encoded rbcL and matK + flanking intergenic spacer region data. Molecular data were employed because recent morphological and anatomical studies of tea (Camellia sinensis) and related plant species provide conflicting support for tribal and generic relationships within the family. Parsimony analyses of separate and combined data consistently identify three strongly supported lineages: Theeae, Stewartieae, and Gordonieae. These data support the broad generic circumscription of Camellia and Stewartia but do not support the recognition of Gordonia sensu lato. Gordonia lasianthus and Gordonia brandegeei are the basal clade in Gordonieae, a position far removed from all other representatives of Gordonia sensu lato (Polyspora and Laplacea) included in this study. This phylogeny most closely mirrors Airy-Shaw's tribe Camellieae [= Theeae] and his two subtribes Stewartiinae and Gordoniinae, first published in 1936. We recognize all three major lineages at the tribal level, although there is weak statistical support for a sister relationship between Gordonieae and Theeae. We also find statistical support for the recognition of the two former subfamilies Theoideae and Ternstroemioideae as two separate families, Theaceae and Ternstroemiaceae.

Evolution of the angiosperms: calibrating the family tree

Growing evidence of morphological diversity in angiosperm flowers, seeds and pollen from the mid Cretaceous and the presence of derived lineages from increasingly older geological deposits both imply that the timing of early angiosperm cladogenesis is older than fossil-based estimates have indicated. An alternative to fossils for calibrating the phylogeny comes from divergence in DNA sequence data. Here, angiosperm divergence times are estimated using non-parametric rate smoothing and a three-gene dataset covering ca. 75% of all angiosperm families recognized in recent classifications. The results provide an initial hypothesis of angiosperm diversification times. Using an internal calibration point, an independent evaluation of angiosperm and eudicot origins is performed. The origin of the crown group of extant angiosperms is indicated to be Early to Middle Jurassic (179-158 Myr), and the origin of eudicots is resolved as Late Jurassic to mid Cretaceous (147-131 Myr). Both estimates, despite a conservative calibration point, are older than current fossil-based estimates.

Molecular systematics of the western rattlesnake, Crotalus viridis (Viperidae), with comments on the utility of the D-loop in phylogenetic studies of snakes

Crotalus viridis, the western rattlesnake, ranges throughout western North America and has been divided into at least eight subspecies. However, the validity of and relationships among these subspecies and the monophyly of C. viridis as a whole are questionable. We used mitochondrial DNA sequence data from the D-loop region and ND2 gene to examine the relationships among 26 populations of C. viridis and to test the monophyly of this species. These data were analyzed separately and combined using maximum-likelihood and maximum-parsimony. The C. viridis group was monophyletic in all combined analyses, consisting of two strongly divergent clades. We recommend that these clades be recognized as two distinct evolutionary species: C. viridis and C. oreganus. Crotalus viridis should be restricted to the subspecies viridis and nuntius and the remaining subspecies be assigned to the species C. oreganus. Our data do not allow strong evaluation of the validity of the subspecies. We found that the ND2 gene had greater sequence divergences among closely related individuals than the D-loop region, but this relationship reversed at higher levels of divergence. This pattern is apparently due to: (1) ND2 third positions evolving faster than the D-loop but becoming saturated at higher levels of divergence, and (2) the D-loop evolving faster than ND2 second (and possibly first) positions. Our results suggest that the ND2 gene is preferable for examining intraspecific relationships and the D-loop may better resolve relationships between species of snakes. The latter result is contrary to the common perception of the phylogenetic utility of the D-loop. Another unusual result is that the 145 bp spacer region, adjacent to the 5' end of the light strand of the D-loop, provides greater phylogenetic resolution than the 1030 bp D-loop.

Molecular systematics of Iridaceae: evidence from four plastid DNA regions

Iridaceae are one of the largest families of Lilianae and probably also among the best studied of monocotyledons. To further evaluate generic, tribal, and subfamilial relationships we have produced four plastid DNA data sets for 57 genera of Iridaceae plus outgroups: rps4, rbcL (both protein-coding genes), the trnL intron, and the trnL-F intergenic spacer. All four matrices produce similar although not identical trees, and we thus analyzed them in a combined analysis, which produced a highly resolved and well-supported topology, in spite of the fact that the partition homogeneity test indicated strong incongruence. In each of the individual trees, some genera or groups of genera are misplaced relative to morphological cladistic studies, but the combined analysis produced a pattern much more similar to these previous ideas of relationships. In the combined tree, all subfamilies were resolved as monophyletic, except Nivenioideae that formed a grade in which Ixioideae were embedded. Achlorophyllous Geosiris (sometimes referred to Geosiridaceae or Burmanniaceae) fell within the nivenioid grade. Most of the tribes were monophyletic, and Isophysis (Tasmanian) was sister to the rest of the family; Diplarrhena (Australian) fell in a well-supported position as sister to Irideae/Sisyrinchieae/Tigridieae/Mariceae (i.e., Iridoideae); Bobartia of Sisyrinchieae is supported as a member of Irideae. The paraphyly of Nivenioideae is suspicious due to extremely high levels of sequence divergence, and when they were constrained to be monophyletic the resulting trees were only slightly less parsimonious (<1.0%). However, this subfamily also lacks clear morphological synapomorphies and is highly heterogeneous, so it is difficult to develop a strong case on nonmolecular grounds for their monophyly.

Molecular systematics of Malpighiaceae: evidence from plastid rbcL and matK sequences

Phylogenetic analyses of DNA nucleotide sequences from the plastid genes rbcL and matK were employed to investigate intergeneric relationships within Malpighiaceae. Cladistic relationships generated from the independent data matrices for the family are generally in agreement with those from the combined matrix. At the base of Malpighiaceae are several clades mostly representing genera from a paraphyletic subfamily Byrsonimoideae. Intergeneric relationships among these byrsonimoid malpighs are well supported by the bootstrap, and the tribe Galphimeae is monophyletic. There is also a well-supported clade of genera corresponding to tribes Banisterieae plus Gaudichaudieae present in all trees, and many of the relationships among these banisterioid malpighs are well supported by the bootstrap. However, tribes Hiraeae and Tricomarieae (the hiraeoid malpighs) are paraphyletic and largely unresolved. Species of Mascagnia are distributed throughout these hiraeoid clades, confirming the suspected polyphyly of this large genus. Optimization of selected morphological characters on these trees demonstrates clear phylogenetic trends such as the evolution of globally symmetrical from radially symmetrical pollen, increased modification and sterilization of stamens, and switch from base chromosome number n = 6 to n = 10.

Was the ANITA rooting of the angiosperm phylogeny affected by long-branch attraction? Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya

Five groups of basal angiosperms, Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya (ANITA), were identified in several recent studies as representing a series of the earliest-diverging lineages of the angiosperm phylogeny. All of these studies except one employed a multigene analysis approach and used gymnosperms as the outgroup to determine the ingroup topology. The high level of divergence between gymnosperms and angiosperms, however, has long been implicated in the difficulty of reconstructing relationships at the base of angiosperm phylogeny using DNA sequences, for fear of long-branch attraction (LBA). In this study, we replaced the gymnosperm sequences from the five-gene matrix (mitochondrial atp1 and matR, plastid atpB and rbcL, and nuclear 18S rDNA) used in our earlier study with four categories of divergent sequences--random sequences with equal base frequencies or equally AT- and GC-rich contents, homopolymers and heteropolymers, misaligned gymnosperm sequences, and aligned lycopod and bryophyte sequences--to evaluate whether the gymnosperms were an appropriate outgroup to angiosperms in our earlier study that identified the ANITA rooting. All 24 analyses performed rooted the angiosperm phylogeny at either Acorus or Alisma (or Alisma-Triglochin-Potamogeton in one case due to use of a slightly different alignment) and placed the monocots as a basal grade, producing genuine LBA results. These analyses demonstrate that the identification of ANITA as the basalmost extant angiosperms was based on historical signals preserved in the gymnosperm sequences and that the gymnosperms were an appropriate outgroup with which to root the angiosperm phylogeny in the multigene sequence analysis. This strategy of evaluating the appropriateness of an outgroup using artificial sequences and a series of outgroups with increments of divergence levels can be applied to investigations of phylogenetic patterns at the bases of other major clades, such as land plants, animals, and eukaryotes.

Phylogeny of the Celastraceae inferred from 26S nuclear ribosomal DNA, phytochrome B, rbcL, atpB, and morphology

Phylogenetic relationships within Celastraceae (spindle-tree family) were inferred from nucleotide sequence characters from the 5' end of 26S nuclear ribosomal DNA (including expansion segments D1-D3; 84 species sampled), phytochrome B (58 species), rbcL (31 species), atpB (23 species), and morphology (94 species). Among taxa of questionable affinity, Forsellesia is a member of Crossosomataceae, and Goupia is excluded from Celastraceae. However, Brexia, Canotia, Lepuropetalon, Parnassia, Siphonodon, and Stackhousiaceae are supported as members of Celastraceae. Gymnosporia and Tricerma are distinct from Maytenus, Cassine is supported as distinct from Elaeodendron, and Dicarpellum is distinct from Salacia. Catha, Maytenus, and Pristimera are not resolved as natural genera. Hippocrateaceae (including Plagiopteron and Lophopetalum) are a clade nested within a paraphyletic Celastraceae. These data also suggest that the Loesener's classification of Celastraceae sensu stricto and Hallé's classification of Hippocrateaceae are artificial. The diversification of the fruit and aril within Celastraceae appears to be complex, with multiple origins of most fruit and aril forms.

Feeding specialization and host-derived chemical defense in Chrysomeline leaf beetles did not lead to an evolutionary dead end

Combination of molecular phylogenetic analyses of Chrysomelina beetles and chemical data of their defensive secretions indicate that two lineages independently developed, from an ancestral autogenous metabolism, an energetically efficient strategy that made the insect tightly dependent on the chemistry of the host plant. However, a lineage (the interrupta group) escaped this subordination through the development of a yet more derived mixed metabolism potentially compatible with a large number of new host-plant associations. Hence, these analyses on leaf beetles document a mechanism that can explain why high levels of specialization do not necessarily lead to "evolutionary dead ends."

Ancient flowering plants: DNA sequences and angiosperm classification

Phylogenetic analyses of gene sequences provide a clear pattern of which extant flowering plant genera diversified earliest. Combined with complete genomic sequences, these data will vastly improve understanding of the genetic basis of plant diversity.

The granule-bound starch synthase (GBSSI) gene in the Rosaceae: multiple loci and phylogenetic utility

We sampled the 5' end of the granule-bound starch synthase gene (GBSSI or waxy) in Rosaceae, sequencing 108 clones from 18 species in 14 genera representing all four subfamilies (Amygdaloideae, Maloideae, Rosoideae, and Spiraeoideae), as well as four clones from Rhamnus catharticus (Rhamnaceae). This is the first phylogenetic study to use the 5' portion of this nuclear gene. Parsimony and maximum-likelihood analyses of 941 bases from seven complete and two partial exons demonstrate the presence of two loci (GBSSI-1 and GBSSI-2) in the Rosaceae. Southern hybridization analyses with locus-specific probes confirm that all four Rosaceae subfamilies have at least two GBSSI loci, even though only one locus has been reported in all previously studied diploid flowering plants. Phylogenetic analyses also identify four clades representing four loci in the Maloideae. Phylogenetic relationships inferred from GBSSI sequences are largely compatible with those from chloroplast (cpDNA: ndhF, rbcL) and nuclear ribosomal internal transcribed spacer (nrITS) DNA. Large clades are marked by significant intron variation: a long first intron plus no sixth intron in Maloideae GBSSI-1, a long fourth intron in Rosoideae GBSSI-1, and a GT to GC mutation in the 5' splice site of the fourth intron in all GBSSI-2 sequences. Our data do not support the long-held hypothesis that Maloideae originated from an ancient hybridization between amygdaloid and spiraeoid ancestors. Instead, Spiraeoideae genera (Kageneckia and Vauquelinia) are their closest relatives in all four GBSSI clades.

Subtribal and generic relationships of Maxillarieae (Orchidaceae) with emphasis on Stanhopeinae: combined molecular evidence

The monophyly of and phylogenetic relationships within the orchid tribe Maxillarieae Pfitzer were evaluated using parsimony analyses of combined nuclear ribosomal and plastid DNA sequence data of ITS 1 and 2, matK, and the trnL intron and the trnL-F intergene spacer. Each of the separate analyses produced highly congruent but weakly supported patterns (by the bootstrap), so these were combined in a single analysis. Analysis of 90 ingroup taxa (representing ∼35% of currently recognized genera) and four outgroup taxa produced resolved and highly supported cladograms. Based on the cladograms, we recognize six subtribes: Eriopsidinae, Oncidiinae (including Pachyphyllinae, Ornithocephalinae, and Telipogoninae), Stanhopeinae, Coeliopsidinae, Maxillariinae (including Lycastinae and Bifrenariinae), and Zygopetalinae (including Cryptarrheninae, Dichaeinae, Huntleyinae, and Warreinae). Stanhopeinae were sampled most intensively; their generic relationships were highly resolved in the analysis and largely agree with currently accepted generic concepts based on morphology. Coeliopsidinae (Coeliopsis, Lycomormium, Peristeria) are sister to Stanhopeinae. Correlations are drawn among phylogeny, pollination mechanisms, and life history traits.

Utility of 17 chloroplast genes for inferring the phylogeny of the basal angiosperms

Sequences from 14 slowly evolving chloroplast genes (including three highly conserved introns) were obtained for representative basal angiosperm and seed-plant taxa, using novel primers described here. These data were combined with published sequences from atpB, rbcL, and newly obtained sequences from ndhF. Combined data from these 17 genes permit sturdy, well-resolved inference of major aspects of basal angiosperm relationships, demonstrating that the new primers are valuable tools for sorting out the deepest events in flowering plant phylogeny. Sequences from the inverted repeat (IR) proved to be particularly reliable (low homoplasy, high retention index). Representatives of Cabomba and Illicium were the first two successive branches of the angiosperms in an initial sampling of 19 exemplar taxa. This result was strongly supported by bootstrap analysis and by two small insertion/deletion events in the slowly evolving introns. Several paleoherb groups (representatives of Piperales) formed a strongly supported clade with taxa representing core woody magnoliids (Laurales, Magnoliales, and Winteraceae). The monophyly of the sampled eudicots and monocots was also well supported. Analyses of three major partitions of the data showed many of the same clades and supported the rooting seen with all the data combined. While Amborella trichopoda was supported as the sister group of the remaining angiosperms when we added Amborella and Nymphaea odorata to the analysis, a strongly conflicting rooting was observed when Amborella alone was added.

Recent progress in reconstructing angiosperm phylogeny

In the past year, the study of angiosperm phylogeny has moved from tentative inferences based on relatively small data matrices into an era of sophisticated, multigene analyses and significantly greater confidence. Recent studies provide both strong statistical support and mutual corroboration for crucial aspects of angiosperm phylogeny. These include identifying the earliest extant lineages of angiosperms, confirming Amborella as the sister of all other angiosperms, confirming some previously proposed lineages and redefining other groups consistent with their phylogeny. This phylogenetic framework enables the exploration of both genotypic and phenotypic diversification among angiosperms.

Phylogenetic relationships within the Gentianales based on NDHF and RBCL sequences, with particular reference to the Loganiaceae

Phylogenetic relationships in the Gentianales with focus on Loganiaceae sensu lato are evaluated using parsimony analyses of nucleotide sequence data from the plastid genes rbcL and ndhF. Inter- and intrafamilial relationships in the Gentianales, which consist of the families Apocynaceae (including Asclepiadaceae), Gelsemiaceae, Gentianaceae, Loganiaceae, and Rubiaceae, are studied and receive increased support from the combination of rbcL and ndhF data, which indicate that the family Rubiaceae forms the sister group to the successively nested Gentianaceae, Apocynaceae, and Loganiaceae, all of which are well supported. The family Gelsemiaceae forms a distinct, supported group sister to Apocynaceae. The Loganiaceae sensu stricto form a strongly supported group consisting of 13 genera: Antonia, Bonyunia, Gardneria, Geniostoma, Labordia, Logania, Mitrasacme, Mitreola, Neuburgia, Norrisia, Spigelia, Strychnos, and Usteria. These genera form two well-supported lineages. Several members of Loganiaceae sensu Leeuwenberg and Leenhouts, i.e., Androya, Peltanthera, Plocosperma, Polypremum, and Sanango are clearly not members of the Gentianales. The earlier exclusion of Buddlejaceae (including Buddleja, Emorya, Gomphostigma, and Nicodemia) as well as the reclassification of the genera Nuxia and Retzia to Stilbaceae of the Lamiales are all well supported.