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

Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods

BACKGROUND AND AIMS

Despite a recent new classification, a stable phylogeny for the cycads has been elusive, particularly regarding resolution of Bowenia, Stangeria and Dioon. In this study, five single-copy nuclear genes (SCNGs) are applied to the phylogeny of the order Cycadales. The specific aim is to evaluate several gene tree-species tree reconciliation approaches for developing an accurate phylogeny of the order, to contrast them with concatenated parsimony analysis and to resolve the erstwhile problematic phylogenetic position of these three genera.

METHODS

DNA sequences of five SCNGs were obtained for 20 cycad species representing all ten genera of Cycadales. These were analysed with parsimony, maximum likelihood (ML) and three Bayesian methods of gene tree-species tree reconciliation, using Cycas as the outgroup. A calibrated date estimation was developed with Bayesian methods, and biogeographic analysis was also conducted.

KEY RESULTS

Concatenated parsimony, ML and three species tree inference methods resolve exactly the same tree topology with high support at most nodes. Dioon and Bowenia are the first and second branches of Cycadales after Cycas, respectively, followed by an encephalartoid clade (Macrozamia-Lepidozamia-Encephalartos), which is sister to a zamioid clade, of which Ceratozamia is the first branch, and in which Stangeria is sister to Microcycas and Zamia.

CONCLUSIONS

A single, well-supported phylogenetic hypothesis of the generic relationships of the Cycadales is presented. However, massive extinction events inferred from the fossil record that eliminated broader ancestral distributions within Zamiaceae compromise accurate optimization of ancestral biogeographical areas for that hypothesis. While major lineages of Cycadales are ancient, crown ages of all modern genera are no older than 12 million years, supporting a recent hypothesis of mostly Miocene radiations. This phylogeny can contribute to an accurate infrafamilial classification of Zamiaceae.

Different gymnosperm outgroups have (mostly) congruent signal regarding the root of flowering plant phylogeny

We examined multiple plastid genes from a diversity of gymnosperm lineages to explore the consistency of signal among different outgroups for rooting flowering plant phylogeny. For maximum parsimony (MP), most outgroups attach on a branch of the underlying ingroup tree that leads to Amborella. Maximum likelihood (ML) analyses either root angiosperms on a nearby branch or find split support for these neighboring root placements, depending on the outgroup. The inclusion of two species of Hydatellaceae, recently recognized as an ancient line of angiosperms, does not aid in inference of the root. Cost profiles for placing the root in suboptimal locations are highly correlated across most outgroup comparisons, even comparing MP and ML profiles. Those for Gnetales are the most deviant of all those considered. This divergent outgroup either attaches on a long eudicot branch with moderate bootstrap support in MP analyses or supports no particular root location in ML analysis. Removing the most rapidly evolving sites in rate classifications based on two divergent angiosperm root placements with Gnetales yields strongly conflicting root placements in MP analysis, despite substantial overlap in the estimated sets of conservative sites. However, the generally high consistency in rooting signal among distantly related gymnosperm clades suggests that the long branch connecting angiosperms to their extant relatives may not interfere substantially with inference of the angiosperm root.

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.

Experimental Design Criteria in Phylogenetics: Where to Add Taxa

Accurate phylogenetic inference is a topic of intensive research and debate and has been studied in response to many different factors: for example, differences in the method of reconstruction, the shape of the underlying tree, the substitution model, and varying quantities and types of data. Investigating whether the conditions used might lead to inaccurate inference has been attempted through elaborate data exploration but less attention has been given to creating a unified methodology to enable experimental designs in phylogenetic analysis to be improved and so avoid suboptimal conditions. Experimental design has been part of the field of statistics since the seminal work of Fisher in the early 20th century and a large body of literature exists on how to design optimum experiments. Here we investigate the use of the Fisher information matrix to decide between candidate positions for adding a taxon to a fixed topology, and introduce a parameter transformation that permits comparison of these different designs. This extension to Goldman (1998. Proc. R. Soc. Lond. B. 265: 1779-1786) thus allows investigation of "where to add taxa" in a phylogeny. We compare three different measures of the total information for selecting the position to add a taxon to a tree. Our methods are illustrated by investigating the behavior of the three criteria when adding a branch to model trees, and by applying the different criteria to two biological examples: a simplified taxon-sampling problem in the balsaminoid Ericales and the phylogeny of seed plants.

Duplication and paralog sorting of RPB2 and RPB1 genes in core eudicots

RPB1 and RPB2, which encode the largest and second largest subunits of RNA polymerase II, respectively, are essential single copy genes in fungi, animals and most plants. Two paralogs of the RPB2 gene have been found in some groups of angioperms [Oxelman, B., Yoshikawa, N., McConaughy, B.L., Luo, J., Denton, A.L., Hall, B.D., 2004. RPB2 gene phylogeny in flowering plants, with particular emphasis on asterids. Mol. Phylogenet. Evol. 32, 462-479]. Here, we report the results of experiments designed to identify the evolutionary origin of the RPB2 duplicate copies. Through careful sampling and phylogenetic analysis, we were able to construct the RPB2 gene tree in angiosperms and infer the phylogenetic positions of the gene duplication and gene loss events that occurred. Our study shows that an RPB2 gene duplication occurred early in core eudicot evolution, at or near the time of the Buxaceae/Trochodendraceae divergence. Subsequently, multiple gene duplication and paralog sorting events happened independently in different core eudicot taxa. Differential expression of the two RPB2 gene paralogs may explain the preservation of both paralogs in the asterids. One gene (RPB2-i) accounts for most of the RPB2 mRNA made in the flower organs while the other gene (RPB2-d) is predominantly used in the vegetative tissues. We also found two paralogs of the RPB1 gene in some core eudicot species. The RPB1 gene duplication occurred before core eudicot divergence, around the time of RPB2 gene duplication. Several independent RPB1 paralog sorting events happened in different core eudicot taxa; their occurrence was independent of the RPB2 paralog sorting events. Our results suggest that a polyploidization event happened at or near the time of the Buxaceae/Trochodendraceae divergence. We propose that this polyploidization and the partial diploidization processes thereafter may have been the driving force of core eudicot radiation.

Phylogenetic utility of protein (RPB2, β-tubulin) and ribosomal (LSU, SSU) gene sequences in the systematics of Sordariomycetes (Ascomycota, Fungi)

The Sordariomycetes is an important group of fungi whose taxonomic relationships and classification is obscure. There is presently no multi-gene molecular phylogeny that addresses evolutionary relationships among different classes and orders. In this study, phylogenetic analyses with a broad taxon sampling of the Sordariomycetes were conducted to evaluate the utility of four gene regions (LSU rDNA, SSU rDNA, beta-tubulin and RPB2) for inferring evolutionary relationships at different taxonomic ranks. Single and multi-gene genealogies inferred from Bayesian and Maximum Parsimony analyses were compared in individual and combined datasets. At the subclass level, SSU rDNA phylogenies demonstrate their utility as a marker to infer phylogenetic relationships at higher levels. All analyses with SSU rDNA alone, combined LSU rDNA and SSU rDNA, and the combined 28 S rDNA, SSU rDNA and RPB2 datasets resulted in three subclasses: Hypocreomycetidae, Sordariomycetidae and Xylariomycetidae, which correspond well to established morphological classification schemes. At the ordinal level, the best resolved phylogeny was obtained from the combined LSU rDNA and SSU rDNA datasets. Individually, the RPB2 gene dataset resulted in significantly higher number of parsimony informative characters. Our results supported the recent separation of Boliniaceae, Chaetosphaeriaceae and Coniochaetaceae from Sordariales and placement of Coronophorales in Hypocreomycetidae. Microascales was found to be paraphyletic and Ceratocystis is phylogenetically associated to Faurelina, while Microascus and Petriella formed another clade and basal to other members of Halosphaeriales. In addition, the order Lulworthiales does not appear to fit in any of the three subclasses. Congruence between morphological and molecular classification schemes is discussed.

Seed plant phylogeny: gnetophytes are derived conifers and a sister group to Pinaceae

The phylogenetic position of gnetophytes has long been controversial. We sequenced parts of the genes coding for the largest subunit of nuclear RNA polymerase I, II, and III and combined these sequences with those of four chloroplast genes, two mitochondrial genes, and 18S rRNA genes to address this issue. Both maximum likelihood and maximum parsimony analyses of the sites not affected by high substitution levels strongly support a phylogeny where gymnosperms and angiosperms are monophyletic, where cycads are at the base of gymnosperm tree and are followed by ginkgos, and where gnetophytes are grouped within conifers as the sister group of pines. The evolution of several morphological and molecular characters of gnetophytes and conifers will therefore need to be reinterpreted.

Expression of floral MADS-box genes in basal angiosperms: implications for the evolution of floral regulators

The ABC model of floral organ identity is based on studies of Arabidopsis and Antirrhinum, both of which are highly derived eudicots. Most of the genes required for the ABC functions in Arabidopsis and Antirrhinum are members of the MADS-box gene family, and their orthologs are present in all major angiosperm lineages. Although the eudicots comprise 75% of all angiosperms, most of the diversity in arrangement and number of floral parts is actually found among basal angiosperm lineages, for which little is known about the genes that control floral development. To investigate the conservation and divergence of expression patterns of floral MADS-box genes in basal angiosperms relative to eudicot model systems, we isolated several floral MADS-box genes and examined their expression patterns in representative species, including Amborella (Amborellaceae), Nuphar (Nymphaeaceae) and Illicium (Austrobaileyales), the successive sister groups to all other extant angiosperms, plus Magnolia and Asimina, members of the large magnoliid clade. Our results from multiple methods (relative-quantitative RT-PCR, real-time PCR and RNA in situ hybridization) revealed that expression patterns of floral MADS-box genes in basal angiosperms are broader than those of their counterparts in eudicots and monocots. In particular, (i) AP1 homologs are generally expressed in all floral organs and leaves, (ii) AP3/PI homologs are generally expressed in all floral organs and (iii) AG homologs are expressed in stamens and carpels of most basal angiosperms, in agreement with the expectations of the ABC model; however, an AG homolog is also expressed in the tepals of Illicium. The broader range of strong expression of AP3/PI homologs is inferred to be the ancestral pattern for all angiosperms and is also consistent with the gradual morphological intergradations often observed between adjacent floral organs in basal angiosperms.

Improving phylogenetic inference of mushrooms with RPB1 and RPB2 nucleotide sequences (Inocybe; Agaricales)

Approximately 3000 bp across 84 taxa have been analyzed for variable regions of RPB1, RPB2, and nLSU-rDNA to infer phylogenetic relationships in the large ectomycorrhizal mushroom genus Inocybe (Agaricales; Basidiomycota). This study represents the first effort to combine variable regions of RPB1 and RPB2 with nLSU-rDNA for low-level phylogenetic studies in mushroom-forming fungi. Combination of the three loci increases non-parametric bootstrap support, Bayesian posterior probabilities, and resolution for numerous clades compared to separate gene analyses. These data suggest the evolution of at least five major lineages in Inocybe-the Inocybe clade, the Mallocybe clade, the Auritella clade, the Inosperma clade, and the Pseudosperma clade. Additionally, many clades nested within each major lineage are strongly supported. These results also suggest the family Crepiodataceae sensu stricto is sister to Inocybe. Recognition of Inocybe at the family level, the Inocybaceae, is recommended.

The origin and diversification of angiosperms

The angiosperms, one of five groups of extant seed plants, are the largest group of land plants. Despite their relatively recent origin, this clade is extremely diverse morphologically and ecologically. However, angiosperms are clearly united by several synapomorphies. During the past 10 years, higher-level relationships of the angiosperms have been resolved. For example, most analyses are consistent in identifying Amborella, Nymphaeaceae, and Austrobaileyales as the basalmost branches of the angiosperm tree. Other basal lineages include Chloranthaceae, magnoliids, and monocots. Approximately three quarters of all angiosperm species belong to the eudicot clade, which is strongly supported by molecular data but united morphologically by a single synapomorphy-triaperturate pollen. Major clades of eudicots include Ranunculales, which are sister to all other eudicots, and a clade of core eudicots, the largest members of which are Saxifragales, Caryophyllales, rosids, and asterids. Despite rapid progress in resolving angiosperm relationships, several significant problems remain: (1) relationships among the monocots, Chloranthaceae, magnoliids, and eudicots, (2) branching order among basal eudicots, (3) relationships among the major clades of core eudicots, (4) relationships within rosids, (5) relationships of the many lineages of parasitic plants, and (6) integration of fossils with extant taxa into a comprehensive tree of angiosperm phylogeny.

Phylogenetic signal in nucleotide data from seed plants: implications for resolving the seed plant tree of life

Effects of taxonomic sampling and conflicting signal on the inference of seed plant trees supported in previous molecular analyses were explored using 13 single-locus data sets. Changing the number of taxa in single-locus analyses had limited effects on log likelihood differences between the gnepine (Gnetales plus Pinaceae) and gnetifer (Gnetales plus conifers) trees. Distinguishing among these trees also was little affected by the use of different substitution parameters. The 13-locus combined data set was partitioned into nine classes based on substitution rates. Sites evolving at intermediate rates had the best likelihood and parsimony scores on gnepine trees, and those evolving at the fastest rates had the best parsimony scores on Gnetales-sister trees (Gnetales plus other seed plants). When the fastest evolving sites were excluded from parsimony analyses, well-supported gnepine trees were inferred from the combined data and from each genomic partition. When all sites were included, Gnetales-sister trees were inferred from the combined data, whereas a different tree was inferred from each genomic partition. Maximum likelihood trees from the combined data and from each genomic partition were well-supported gnepine trees. A preliminary stratigraphic test highlights the poor fit of Gnetales-sister trees to the fossil data.

Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?

Background

Numerous studies, using in aggregate some 28 genes, have achieved a consensus in recognizing three groups of plants, including Amborella, as comprising the basal-most grade of all other angiosperms. A major exception is the recent study by Goremykin et al. (2003; Mol. Biol. Evol. 20:1499–1505), whose analyses of 61 genes from 13 sequenced chloroplast genomes of land plants nearly always found 100% support for monocots as the deepest angiosperms relative to Amborella, Calycanthus, and eudicots. We hypothesized that this conflict reflects a misrooting of angiosperms resulting from inadequate taxon sampling, inappropriate phylogenetic methodology, and rapid evolution in the grass lineage used to represent monocots.

Results

We used two main approaches to test this hypothesis. First, we sequenced a large number of chloroplast genes from the monocot Acorus and added these plus previously sequenced Acorus genes to the Goremykin et al. (2003) dataset in order to explore the effects of altered monocot sampling under the same analytical conditions used in their study. With Acorus alone representing monocots, strongly supported Amborella-sister trees were obtained in all maximum likelihood and parsimony analyses, and in some distance-based analyses. Trees with both Acorus and grasses gave either a well-supported Amborella-sister topology or else a highly unlikely topology with 100% support for grasses-sister and paraphyly of monocots (i.e., Acorus sister to "dicots" rather than to grasses). Second, we reanalyzed the Goremykin et al. (2003) dataset focusing on methods designed to account for rate heterogeneity. These analyses supported an Amborella-sister hypothesis, with bootstrap support values often conflicting strongly with cognate analyses performed without allowing for rate heterogeneity. In addition, we carried out a limited set of analyses that included the chloroplast genome of Nymphaea, whose position as a basal angiosperm was also, and very recently, challenged.

Conclusions

These analyses show that Amborella (or Amborella plus Nymphaea), but not monocots, is the sister group of all other angiosperms among this limited set of taxa and that the grasses-sister topology is a long-branch-attraction artifact leading to incorrect rooting of angiosperms. These results highlight the danger of having lots of characters but too few and, especially, molecularly divergent taxa, a situation long recognized as potentially producing strongly misleading molecular trees. They also emphasize the importance in phylogenetic analysis of using appropriate evolutionary models.

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.