Summarizing a posterior distribution of trees using agreement subtrees

Bayesian inference of phylogeny is unique among phylogenetic reconstruction methods in that it produces a posterior distribution of trees rather than a point estimate of the best tree. The most common way to summarize this distribution is to report the majority-rule consensus tree annotated with the marginal posterior probabilities of each partition. Reporting a single tree discards information contained in the full underlying distribution and reduces the Bayesian analysis to simply another method for finding a point estimate of the tree. Even when a point estimate of the phylogeny is desired, the majority-rule consensus tree is only one possible method, and there may be others that are more appropriate for the given data set and application. We present a method for summarizing the distribution of trees that is based on identifying agreement subtrees that are frequently present in the posterior distribution. This method provides fully resolved binary trees for subsets of taxa with high marginal posterior probability on the entire tree and includes additional information about the spread of the distribution.     

Trees of trees: an approach to comparing multiple alternative phylogenies

Phylogenetic analysis very commonly produces several alternative trees for a given fixed set of taxa. For example, different sets of orthologous genes may be analyzed, or the analysis may sample from a distribution of probable trees. This article describes an approach to comparing and visualizing multiple alternative phylogenies via the idea of a "tree of trees" or "meta-tree." A meta-tree clusters phylogenies with similar topologies together in the same way that a phylogeny clusters species with similar DNA sequences. Leaf nodes on a meta-tree correspond to the original set of phylogenies given by some analysis, whereas interior nodes correspond to certain consensus topologies. The construction of meta-trees is motivated by analogy with construction of a most parsimonious tree for DNA data, but instead of using DNA letters, in a meta-tree the characters are partitions or splits of the set of taxa. An efficient algorithm for meta-tree construction is described that makes use of a known relationship between the majority consensus and parsimony in terms of gain and loss of splits. To illustrate these ideas meta-trees are constructed for two datasets: a set of gene trees for species of yeast and trees from a bootstrap analysis of a set of gene trees in ray-finned fish. A software tool for constructing meta-trees and comparing alternative phylogenies is available online, and the source code can be obtained from the author.     

Imputing supertrees and supernetworks from quartets

Inferring species phylogenies is an important part of understanding molecular evolution. Even so, it is well known that an accurate phylogenetic tree reconstruction for a single gene does not always necessarily correspond to the species phylogeny. One commonly accepted strategy to cope with this problem is to sequence many genes; the way in which to analyze the resulting collection of genes is somewhat more contentious. Supermatrix and supertree methods can be used, although these can suppress conflicts arising from true differences in the gene trees caused by processes such as lineage sorting, horizontal gene transfer, or gene duplication and loss. In 2004, Huson et al. (IEEE/ACM Trans. Comput. Biol. Bioinformatics 1:151-158) presented the Z-closure method that can circumvent this problem by generating a supernetwork as opposed to a supertree. Here we present an alternative way for generating supernetworks called Q-imputation. In particular, we describe a method that uses quartet information to add missing taxa into gene trees. The resulting trees are subsequently used to generate consensus networks, networks that generalize strict and majority-rule consensus trees. Through simulations and application to real data sets, we compare Q-imputation to the matrix representation with parsimony (MRP) supertree method and Z-closure, and demonstrate that it provides a useful complementary tool.     

A molecular phylogeny of enteric bacteria and implications for a bacterial species concept

A molecular phylogeny for seven taxa of enteric bacteria (Citrobacter freundii, Enterobacter cloacae, Escherichia coli, Hafnia alvei, Klebsiella oxytoca, Klebsiella pneumoniae, and Serratia plymuthica) was made from multiple isolates per taxa taken from a collection of environmental enteric bacteria. Sequences from five housekeeping genes (gapA, groEL, gyrA, ompA, and pgi) and the 16S rRNA gene were used to infer individual gene trees and were concatenated to infer a composite molecular phylogeny for the species. The isolates from each taxa formed tight species clusters in the individual gene trees, suggesting the existence of 'genotypic' clusters that correspond to traditional species designations. These sequence data and the resulting gene trees and consensus tree provide the first data set with which to assess the utility of the recently proposed core genome hypothesis (CGH). The CGH provides a genetically based approach to applying the biological species concept to bacteria.     

The systematic utility of floral and vegetative fragrance in two genera of nyctaginaceae

We examined relationships between fragrance and phylogeny using a number of approaches to coding fragrance data and comparing the hierarchical information in fragrance data with the phylogenetic signal in a DNA sequence data set. We first used distance analyses to determine which coding method(s) best distinguishes species while grouping conspecifics. Results suggest that interspecific differences in fragrance composition were maximized by coding as presence/absence of fragrance compounds and biosynthetic pathways rather than when quantitative information was also included. Useful systematic information came from both compounds and pathways and from fragrance emitted by both floral and vegetative tissues. The coding methods that emerged from the distance analyses as best distinguishing species were then adapted for use in phylogenetic analysis. Although hierarchical signal among fragrance data sets was congruent, this signal was highly incongruent with the phylogenetic signal in the DNA sequence data. Notably, topologies inferred from fragrance data sets were congruent with the DNA topology only in the most distal portions (e.g., sister group pairs or closely related species that had similar fragrance profiles were often recovered by analyses of fragrance). Examination of consistency and retention indices for individual fragrance compounds and pathways as optimized onto one of the most-parsimonious trees inferred from DNA data revealed that although most compounds were homoplastic, some compounds were perfectly congruent with the DNA phylogeny. In particular, compounds and pathways found in a few taxa were less homoplastic than those found in many taxa. Pathways that synthesize few volatiles also seem to have lower homoplasy than those that produce many. Although fragrance data as a whole may not be useful in phylogeny reconstruction, these data can provide additional support for clades reconstructed with other types of characters. Factors other than phylogeny, including pollinator interactions, also likely influence fragrance composition.     

Phylogenetic relationships of spatangoid sea urchins (Echinoidea): taxon sampling density and congruence between morphological and molecular estimates

A phylogeny for 21 species of spatangoid sea urchins is constructed using data from three genes and results compared with morphology-based phylogenies derived for the same taxa and for a much larger sample of 88 Recent and fossil taxa. Different data sets and methods of analysis generate different phylogenetic hypotheses, although congruence tests show that all molecular approaches produce trees that are congruent with each other. By contrast, the trees generated from morphological data differ significantly according to taxon sampling density and only those with dense sampling (after a posteriori weighting) are congruent with molecular estimates. With limited taxon sampling, secondary reversals in deep-water taxa are interpreted as plesiomorphies, pulling them to a basal position. The addition of fossil taxa with their unique character combinations reveals hidden homoplasy and generates a phylogeny that is compatible with molecular estimates. As homoplasy levels were found to be broadly similar across different anatomical structures in the echinoid test, no one suite of morphological characters can be considered to provide more reliable phylogenetic information. Some traditional groupings are supported, including the grouping of Loveniidae, Brissidae and Spatangidae within the Micrasterina, but the Asterostomatidae is shown to be polyphyletic with members scattered amongst at least five different clades. As these are mostly deep-sea taxa, this finding implies multiple independent invasions into the deep sea.     

Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate

Phylogenies based on mitochondrial DNA (mtDNA) may represent gene trees that may not be congruent with the equivalent species tree. One solution to this problem is to include additional, independent loci from the nuclear genome. Sequence data from the seventh intron of the beta-fibrinogen gene were generated for 25 specimens of vipers, including 8 nominal species of the Trimeresurus complex of Asian pit vipers. Phylogenetic trees were generated using maximum-parsimony and maximum-likelihood methods. The taxonomic level at which the intron provided significant phylogenetic information was examined and the trees were compared to those produced from previously obtained mtDNA cytochrome b sequences. A variety of different approaches (separate analyses, conditional data combination, and consensus) were used in an attempt to provide a sound organismal phylogeny based on both nuclear and mtDNA data sets. We discuss the implications for the gene tree-species tree debate and its particular relevance to medically important organisms.     

Total evidence, consensus, and bat phylogeny: A distance-based approach

Resolution of the total evidence (i.e., character congruence) versus consensus (i.e., taxonomic congruence) debate has been impeded by (1) a failure to employ validation methods consistently across both tree-building and consensus analyses, (2) the incomparability of methods for constructing as opposed to those for combining trees, and (3) indifference to aspects of trees other than their topologies. We demonstrate a uniform, distance-based approach which allows for comparability among the results of character- and taxonomic-congruence studies, whether or not an identical suite of taxa has been included in all contributing data sets. Our results indicate that total-evidence and consensus trees differ little in topology if branch lengths are taken into account when combining two or more trees. In addition, when character-state data are converted to distances, our method permits their combination with information produced by techniques which generate distances directly. Moreover, treating all data sets or trees as distance matrices avoids the problem that different numbers of characters in contributing studies may confound the conclusions of a total-evidence or consensus analysis. Our protocol is illustrated with an example involving bats, in which the three component studies based on serology, DNA hybridization, and anatomy imply distinct phylogenies. However, the total-evidence and consensus trees support a fourth, somewhat different, topology resolved at all but one node and which conforms closely to the currently accepted higher category classification of Chiroptera.     

Combined molecular phylogenetic analysis of the Orthoptera (Arthropoda,Insecta) and implications for their higher systematics

A phylogenetic analysis of mitochondrial and nuclear rDNA sequences from species of all the superfamilies of the insect order Orthoptera (grasshoppers, crickets, and relatives) confirmed that although mitochondrial sequences provided good resolution of the youngest superfamilies, nuclear rDNA sequences were necessary to separate the basal groups. To try to reconcile these data sets into a single, fully resolved orthopteran phylogeny, we adopted consensus and combined data strategies. The consensus analysis produced a partially resolved tree that lacked several well-supported features of the individual analyses. However, this lack of resolution was explained by an examination of resampled data sets, which identified the likely source of error as the relatively short length of the individual mitochondrial data partitions. In a subsequent comparison in which the mitochondrial sequences were initially combined, we observed less conflict. We then used two approaches to examine the validity of combining all of the data in a single analysis: comparative analysis of trees recovered from resampled data sets, and the application of a randomization test. Because the results did not point to significant levels of heterogeneity in phylogenetic signal between the mitochondrial and nuclear data sets, we therefore proceeded with a combined analysis. Reconstructing phylogenies under the minimum evolution and maximum likelihood optimality criteria, we examined monophyly of the major orthopteran groups, using nonparametric and parametric bootstrap analysis and Kishino-Hasegawa tests. Our analysis suggests that phylogeny reconstruction under the maximum likelihood criteria is the most discriminating approach for the combined sequences. The results indicate, moreover, that the caeliferan Pneumoroidea and Pamphagoidea, as previously suggested, are polyphyletic. The Acridoidea is redefined to include all pamphagoid families other than the Pyrgomorphidae, which we propose should be accorded superfamily status.     

Phylogenomic analysis suggests Coreidae and Alydidae (Hemiptera: Heteroptera) are not monophyletic

Next‐generation sequencing technologies (NGS) allow systematists to amass a wealth of genomic data from non‐model species for phylogenetic resolution at various temporal scales. However, phylogenetic inference for many lineages dominated by non‐model species has not yet benefited from NGS, which can complement Sanger sequencing studies. One such lineage, whose phylogenetic relationships remain uncertain, is the diverse, agriculturally important and charismatic Coreoidea (Hemiptera: Heteroptera). Given the lack of consensus on higher‐level relationships and the importance of a robust phylogeny for evolutionary hypothesis testing, we use a large data set comprised of hundreds of ultraconserved element (UCE) loci to infer the phylogeny of Coreoidea (excluding Stenocephalidae and Hyocephalidae), with emphasis on the families Coreidae and Alydidae. We generated three data sets by including alignments that contained loci sampled for at least 50%, 60%, or 70% of the total taxa, and inferred phylogeny using maximum likelihood and summary coalescent methods. Twenty‐six external morphological features used in relatively comprehensive phylogenetic analyses of coreoids were also re‐evaluated within our molecular phylogenetic framework. We recovered 439–970 loci per species (16%–36% of loci targeted) and combined this with previously generated UCE data for 12 taxa. All data sets, regardless of analytical approach, yielded topologically similar and strongly supported trees, with the exception of outgroup relationships and the position of Hydarinae. We recovered a monophyletic Coreoidea, with Rhopalidae highly supported as the sister group to Alydidae + Coreidae. Neither Alydidae nor Coreidae were monophyletic; the coreid subfamilies Hydarinae and Pseudophloeinae were recovered as more closely related to Alydidae than to other coreid subfamilies. Coreinae were paraphyletic with respect to Meropachyinae. Most morphological traits were homoplastic with several clades defined by few, if any, synapomorphies. Our results demonstrate the utility of phylogenomic approaches in generating robust hypotheses for taxa with long‐standing phylogenetic problems and highlight that novel insights may come from such approaches.     

Reconstructing the duplication history of tandemly repeated genes

We present a novel approach to deal with the problem of reconstructing the duplication history of tandemly repeated genes that are supposed to have arisen from unequal recombination. We first describe the mathematical model of evolution by tandem duplication and introduce duplication histories and duplication trees. We then provide a simple recursive algorithm which determines whether or not a given rooted phylogeny can be a duplication history and another algorithm that simulates the unequal recombination process and searches for the best duplication trees according to the maximum parsimony criterion. We use real data sets of human immunoglobulins and T-cell receptors to validate our methods and algorithms. Identity between most parsimonious duplication trees and most parsimonious phylogenies for the same data, combined with the agreement with additional knowledge about the sequences, such as the presence of polymorphisms, shows strong evidence that our reconstruction procedure provides good insights into the duplication histories of these loci.     

Untangling long branches: identifying conflicting phylogenetic signals using spectral analysis, neighbor-net, and consensus networks

Long-branch attraction is a well-known source of systematic error that can mislead phylogenetic methods; it is frequently invoked post hoc, upon recovering a different tree from the one expected based on prior evidence. We demonstrate that methods that do not force the data onto a single tree, such as spectral analysis, Neighbor-Net, and consensus networks, can be used to detect conflicting signals within the data, including those caused by long-branch attraction. We illustrate this approach using a set of taxa from three unambiguously monophyletic families within the Pelecaniformes: the darters, the cormorants and shags, and the gannets and boobies. These three families are universally acknowledged as forming a monophyletic group, but the relationship between the families remains contentious. Using sequence data from three mitochondrial genes (12S, ATPase 6, and ATPase 8) we demonstrate that the relationship between these three families is difficult to resolve because they are separated by a short internal branch and there are conflicting signals due to long-branch attraction, which are confounded with nonhomogeneous sequence evolution across the different genes. Spectral analysis, Neighbor-Net, and consensus networks reveal conflicting signals regarding the placement of one of the darters, with support found for darter monophyly, but also support for a conflicting grouping with the outgroup, pelicans. Furthermore, parsimony and maximum-likelihood analyses produced different trees, with one of the two most parsimonious trees not supporting the monophyly of the darters. Monte Carlo simulations, however, were not sensitive enough to reveal long-branch attraction unless the branches are longer than those actually observed. These results indicate that spectral analysis, Neighbor-Net, and consensus networks offer a powerful approach to detecting and understanding the source of conflicting signals within phylogenetic data.     

Genetic algorithm for large-scale maximum parsimony phylogenetic analysis of proteins

Inferring phylogeny is a difficult computational problem. For example, for only 13 taxa, there are more then 13 billion possible unrooted phylogenetic trees. Heuristics are necessary to minimize the time spent evaluating non-optimal trees. We describe here an approach for heuristic searching, using a genetic algorithm, that can reduce the time required for weighted maximum parsimony phylogenetic inference, especially for data sets involving a large number of taxa. It is the first implementation of a weighted maximum parsimony criterion using amino acid sequences. To validate the weighted criterion, we used an artificial data set and compared it to a number of other phylogenetic methods. Genetic algorithms mimic the natural selection's ability to solve complex problems. We have identified several parameters affecting the genetic algorithm. Methods were developed to validate these parameters, ensuring optimal performance. This approach allows the construction of phylogenetic trees with over 200 taxa in practical time on a regular PC.     

Combining data sets with different phylogenetic histories

The possibility that two data sets may have different underlying phylogenetic histories (such as gene trees that deviate from species trees) has become an important argument against combining data in phylogenetic analysis. However, two data sets sampled for a large number of taxa may differ in only part of their histories. This is a realistic scenario and one in which the relative advantages of combined, separate, and consensus analysis become much less clear. I propose a simple methodology for dealing with this situation that involves (1) partitioning the available data to maximize detection of different histories, (2) performing separate analyses of the data sets, and (3) combining the data but considering questionable or unresolved those parts of the combined tree that are strongly contested in the separate analyses (and which therefore may have different histories) until a majority of unlinked data sets support one resolution over another. In support of this methodology, computer simulations suggest that (1) the accuracy of combined analysis for recovering the true species phylogeny may exceed that of either of two separately analyzed data sets under some conditions, particularly when the mismatch between phylogenetic histories is small and the estimates of the underlying histories are imperfect (few characters, high homoplasy, or both) and (2) combined analysis provides a poor estimate of the species tree in areas of the phylogenies with different histories but gives an improved estimate in regions that share the same history. Thus, when there is a localized mismatch between the histories of two data sets, the separate, consensus, and combined analyses may all give unsatisfactory results in certain parts of the phylogeny. Similarly, approaches that allow data combination only after a global test of heterogeneity will suffer from the potential failings of either separate or combined analysis, depending on the outcome of the test. Excision of conflicting taxa is also problematic, in that doing so may obfuscate the position of conflicting taxa within a larger tree, even when their placement is congruent between data sets. Application of the proposed methodology to molecular and morphological data sets for Sceloporus lizards is discussed.     

Probabilistic Analysis Indicates Discordant Gene Trees in Chloroplast Evolution

Abstract Analyses of whole-genome data often reveal that some genes have evolutionary histories that diverge from the majority phylogeny estimated for the entire genome. We present a probabilistic model that deals with heterogeneity among gene trees, implement it via the Gibbs sampler, and apply it to the plastid genome. Plastids and their genomes are transmitted as a single block without recombination, hence homogeneity among gene trees within this genome is expected. Nevertheless, previous work has revealed clear heterogeneity among plastid genes (e.g., Delwiche and Palmer 1996). Other studies, using whole plastid genomes of various algae and land plants, found little additional heterogeneity (Martin et al. 1998; Adachi et al. 2000). We augment the earlier studies by using a data set of 14 taxa: 6 land plants, 2 green algae, a diatom, 2 red algae and a cryptophyte, the cyanelle of the glaucocystophyte Cyanophora, and the blue–green alga Synechocystis as an outgroup. Contrary to the earlier analyses, we cannot find even a single, dominant consensus tree. Therefore, we formulate a probabilistic model that divides the genes into two sets: those that follow the consensus tree and those that have independent gene trees. No particular tree is supported by more than three-fourths of the genes. But the set of genes that follows a certain tree is fairly independent of data processing and the method of analysis. With one possible exception, we find no evidence for collinear or functionally related genes to follow similar trees. The phylogenetic pattern also seems independent of bias in amino acid composition. Among possible explanations for the observed phenomenon, the hypothesis that different genes have different covarion structures is difficult to assess. But gene duplication may be possible through the inverted or direct repeat regions, while horizontal gene transfer seems less likely. In contrast to green algae and land plants, inverted repeat regions in red algae and in Cyanophora show abundant differences among the copies. Thus, genes may get duplicated when they are recruited into the inverted repeat region and one of the two copies may be lost after leaving the inverted repeat region.     

Alignments of mitochondrial genome arrangements: applications to metazoan phylogeny

Mitochondrial genomes provide a valuable dataset for phylogenetic studies, in particular of metazoan phylogeny because of the extensive taxon sample that is available. Beyond the traditional sequence-based analysis it is possible to extract phylogenetic information from the gene order. Here we present a novel approach utilizing these data based on cyclic list alignments of the gene orders. A progressive alignment approach is used to combine pairwise list alignments into a multiple alignment of gene orders. Parsimony methods are used to reconstruct phylogenetic trees, ancestral gene orders, and consensus patterns in a straightforward approach. We apply this method to study the phylogeny of protostomes based exclusively on mitochondrial genome arrangements. We, furthermore, demonstrate that our approach is also applicable to the much larger genomes of chloroplasts.     

Uncovering hidden phylogenetic consensus in large data sets

Many of the steps in phylogenetic reconstruction can be confounded by “rogue” taxa—taxa that cannot be placed with assurance anywhere within the tree, indeed, whose location within the tree varies with almost any choice of algorithm or parameters. Phylogenetic consensus methods, in particular, are known to suffer from this problem. In this paper, we provide a novel framework to define and identify rogue taxa. In this framework, we formulate a bicriterion optimization problem, the relative information criterion, that models the net increase in useful information present in the consensus tree when certain taxa are removed from the input data. We also provide an effective greedy heuristic to identify a subset of rogue taxa and use this heuristic in a series of experiments, with both pathological examples from the literature and a collection of large biological data sets. As the presence of rogue taxa in a set of bootstrap replicates can lead to deceivingly poor support values, we propose a procedure to recompute support values in light of the rogue taxa identified by our algorithm; applying this procedure to our biological data sets caused a large number of edges to move from “unsupported” to “supported” status, indicating that many existing phylogenies should be recomputed and reevaluated to reduce any inaccuracies introduced by rogue taxa. We also discuss the implementation issues encountered while integrating our algorithm into RAxML v7.2.7, particularly those dealing with scaling up the analyses. This integration enables practitioners to benefit from our algorithm in the analysis of very large data sets (up to 2,500 taxa and 10,000 trees, although we present the results of even larger analyses).     

Phylogenetic Relationships of Extinct Cetartiodactyls: Results of Simultaneous Analyses of Molecular, Morphological, and Stratigraphic Data

Although some recent morphological and molecular studies agree that Cetacea is closely related to Hippopotamidae, there is little consensus on the phylogeny within Cetartiodactyla. We addressed this problem by conducting two analyses: (1) a simultaneous cladistic analysis of intrinsic data (morphology and molecules) and (2) a stratocladistic analysis, which included morphological, molecular, and stratigraphic data. Unlike previous simultaneous analyses, we had the opportunity to include data from the recently described hindlimbs of protocetid and pakicetid cetaceans. Our intrinsic dataset includes 73 taxa scored for 8,229 informative characters, of which 208 are morphological and 8,021 molecular. Both analyses supported the exclusion of Mesonychia from Cetartiodactyla and a close phylogenetic relationship between Hippopotamidae and Cetacea. Many polytomies in the strict consensus of the most parsimonious trees for the intrinsic dataset can be attributed to differing positions for Raoellidae, which in some trees is the sister-group to Cetacea. Pruning Raoellidae and 18 other taxa from all most parsimonious produced a fully resolved agreement subtree, which indicates that the Old World taxa Cebochoerus and Mixtotherium are successive stem taxa to Whippomorpha (i.e., Cetacea + Hippopotamidae). The main result of adding stratigraphic information to the intrinsic dataset was that we found fewer most parsimonious trees, which in most respects were congruent with a subset of the shortest trees for the intrinsic dataset. Our stratocladistic analysis supports species of Diacodexis as the most basal cetartiodactyls, a clade of suiform cetartiodactyls, a monophyletic Tylopoda that includes Protoceratidae, and a monophyletic Carnivora. We were unable to identify any pre-Miocene stem taxa to Hippopotamidae, thus its ghost lineage is still 39 million years long. The relatively low Bremer support for many nodes in our trees indicates that our phylogenetic hypotheses should be subjected to further testing.     

Missing data,incomplete taxa,and phylogenetic accuracy

The problem of missing data is often considered to be the most important obstacle in reconstructing the phylogeny of fossil taxa and in combining data from diverse characters and taxa for phylogenetic analysis. Empirical and theoretical studies show that including highly incomplete taxa can lead to multiple equally parsimonious trees, poorly resolved consensus trees, and decreased phylogenetic accuracy. However, the mechanisms that cause incomplete taxa to be problematic have remained unclear. It has been widely assumed that incomplete taxa are problematic because of the proportion or amount of missing data that they bear. In this study, I use simulations to show that the reduced accuracy associated with including incomplete taxa is caused by these taxa bearing too few complete characters rather than too many missing data cells. This seemingly subtle distinction has a number of important implications. First, the so-called missing data problem for incomplete taxa is, paradoxically, not directly related to their amount or proportion of missing data. Thus, the level of completeness alone should not guide the exclusion of taxa (contrary to common practice), and these results may explain why empirical studies have sometimes found little relationship between the completeness of a taxon and its impact on an analysis. These results also (1) suggest a more effective strategy for dealing with incomplete taxa, (2) call into question a justification of the controversial phylogenetic supertree approach, and (3) show the potential for the accurate phylogenetic placement of highly incomplete taxa, both when combining diverse data sets and when analyzing relationships of fossil taxa.     

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.