Exploring the relationship between sequence similarity and accurate phylogenetic trees

We have characterized the relationship between accurate phylogenetic reconstruction and sequence similarity, testing whether high levels of sequence similarity can consistently produce accurate evolutionary trees. We generated protein families with known phylogenies using a modified version of the PAML/EVOLVER program that produces insertions and deletions as well as substitutions. Protein families were evolved over a range of 100-400 point accepted mutations; at these distances 63% of the families shared significant sequence similarity. Protein families were evolved using balanced and unbalanced trees, with ancient or recent radiations. In families sharing statistically significant similarity, about 60% of multiple sequence alignments were 95% identical to true alignments. To compare recovered topologies with true topologies, we used a score that reflects the fraction of clades that were correctly clustered. As expected, the accuracy of the phylogenies was greatest in the least divergent families. About 88% of phylogenies clustered over 80% of clades in families that shared significant sequence similarity, using Bayesian, parsimony, distance, and maximum likelihood methods. However, for protein families with short ancient branches (ancient radiation), only 30% of the most divergent (but statistically significant) families produced accurate phylogenies, and only about 70% of the second most highly conserved families, with median expectation values better than 10(-60), produced accurate trees. These values represent upper bounds on expected tree accuracy for sequences with a simple divergence history; proteins from 700 Giardia families, with a similar range of sequence similarities but considerably more gaps, produced much less accurate trees. For our simulated insertions and deletions, correct multiple sequence alignments did not perform much better than those produced by T-COFFEE, and including sequences with expressed sequence tag-like sequencing errors did not significantly decrease phylogenetic accuracy. In general, although less-divergent sequence families produce more accurate trees, the likelihood of estimating an accurate tree is most dependent on whether radiation in the family was ancient or recent. Accuracy can be improved by combining genes from the same organism when creating species trees or by selecting protein families with the best bootstrap values in comprehensive studies.     

COSPECIATION ANALYSIS OF AN OBLIGATE POLLINATION MUTUALISM: HAVEGLOCHIDION TREES (EUPHORBIACEAE) AND POLLINATING EPICEPHALA MOTHS(GRACILLARIIDAE) DIVERSIFIED IN PARALLEL?

Species-specific obligate pollination mutualism between Glochidion trees (Euphorbiaceae) and Epicephala moths (Gracillariidae) involves a large number of interacting species and resembles the classically known fig-fig wasp and yucca-yucca moth associations. To assess the extent of parallel cladogenesis in Glochidion-Epicephala association, we reconstruct phylogenetic relationships of 18 species of Glochidion using nuclear ribosomal DNA sequences (internal and external transcribed spacers) and those of the corresponding 18 Epicephala species using mitochondrial (the cytochrome oxidase subunit I gene) and nuclear DNA sequences (the arginine kinase and elongation factor-1alpha genes). Based on the obtained phylogenies, we determine whether Glochidion and Epicephala have undergone parallel diversification using several different methods for investigating the level of cospeciation between phylogenies. These tests indicate that there is generally a greater degree of correlation between Glochidion and Epicephala phylogenies than expected in a random association, but the results are sensitive to selection of different phylogenetic hypotheses and analytical methods for evaluating cospeciation. Perfect congruence between phylogenies is not found in this association, which likely resulted from host shift by the moths. The observed significant discrepancy between Glochidion and Epicephala phylogenies implies that the one-to-one specificity between the plants and moths has been maintained through a complex speciation process or that there is an underestimated diversity of association between Glochidion trees and Epicephala moths.     

COMPARISONS OF OBSERVED PHYLOGENETIC TOPOLOGIES WITH NULL EXPECTATIONS AMONG THREE MONOPHYLETIC LINEAGES

Three null models have been proposed to predict the relative frequencies of topologies of phylogenetic trees. One null model assumes each distinguishable n-member tree is equally likely (proportional-to-distinguishable-arrangements model). A second model assumes that each topological type is equally likely (equiprobable model). A third model assumes that the probability of each topological type is determined by random speciation (Markov model). We sampled published phylogenetic trees from three major groups of organisms: division Angiospermae, class Insecta, and superclass Tetrapoda. Our sampling was more restricted than previous studies and was designed to test whether observed topological frequencies were distinguishable from those predicted by the three null models. The pattern of evolution reflected in five-member phylogenetic trees is different from predictions of the equiprobable and Markov model but is indistinguishable from the proportional-to-distinguishable-arrangements model. This indicates that 1) speciation (and/or extinction) is not equally likely among all taxa, even for small phylogenies; or 2) systematists' attempts at reconstructing small phylogenies are, on average, indistinguishable from those expected if they had merely selected a tree at random from the pool of all possible trees. The topology frequencies were not different among the three groups of organisms, suggesting that factors shaping patterns of speciation and extinction are consistent among major taxonomic groups.     

Phylogeny imbalance: taxonomic level matters

Two lines of evidence indicate that the degree of symmetry in phylogenetic topologies differs at different hierarchical levels. First, in a set of 61 phylogenies with superspecific taxa as their terminals, trees were on average more unbalanced (asymmetric) when the species richness of terminals was considered than when it was not. Second, nodes with a given number of higher taxa descended from them were on average more unbalanced than were nodes with the same number of species as descendants. There are several possible reasons--some biological, some artifactual--for the differences. Whatever the reason, these results caution against treating species-level and higher level phylogenies as equivalent when considering tree shape. The imbalance measure adopted here permits the use of trees that contain polytomies, facilitating a larger sample than has been achieved previously.     

Random sampling of constrained phylogenies: conducting phylogenetic analyses when the phylogeny is partially known

Statistical randomization tests in evolutionary biology often require a set of random, computer-generated trees. For example, earlier studies have shown how large numbers of computer-generated trees can be used to conduct phylogenetic comparative analyses even when the phylogeny is uncertain or unknown. These methods were limited, however, in that (in the absence of molecular sequence or other data) they allowed users to assume that no phylogenetic information was available or that all possible trees were known. Intermediate situations where only a taxonomy or other limited phylogenetic information (e.g., polytomies) are available are technically more difficult. The current study describes a procedure for generating random samples of phylogenies while incorporating limited phylogenetic information (e.g., four taxa belong together in a subclade). The procedure can be used to conduct comparative analyses when the phylogeny is only partially resolved or can be used in other randomization tests in which large numbers of possible phylogenies are needed.     

Molecular phylogenetics of the Espeletia complex (Asteraceae): evidence from nrDNA ITS sequences on the closest relatives of an Andean adaptive radiation

The subtribe Espeletiinae (Asteraceae, Heliantheae) comprises morphologically and ecologically diverse plants endemic to the tropical montane paramos of the Andes of Venezuela, Colombia, and Ecuador. Though the ecophysiology and ecology of this adaptive radiation have been well studied, relationships among taxa in the subtribe and between the subtribe and other taxa in the Heliantheae are poorly known. In this study, sequences from the internal transcribed spacer (ITS) region of nuclear ribosomal DNA are used to test previous hypotheses about the phylogenetic position of the Espeletiinae within the Heliantheae and to determine which taxa are the subtribe's closest relatives. Gene phylogenies based on maximum parsimony analyses reveal that the Espeletiinae clade is nested well within the subtribe Melampodiinae and thus should be considered a monophyletic complex of species, not a separate subtribe. The most parsimonious gene trees suggest that the genus Ichthyothere may be the sister taxon to the Espeletia complex and that the genus Smallanthus and a species of Rumfordia are likely among the complex's other closest living relatives. These data offer preliminary insights into the origins of this adaptive radiation and the broader phylogenetic context in which it occurred.     

Euclidean nature of phylogenetic distance matrices

Phylogenies are fundamental to comparative biology as they help to identify independent events on which statistical tests rely. Two groups of phylogenetic comparative methods (PCMs) can be distinguished: those that take phylogenies into account by introducing explicit models of evolution and those that only consider phylogenies as a statistical constraint and aim at partitioning trait values into a phylogenetic component (phylogenetic inertia) and one or multiple specific components related to adaptive evolution. The way phylogenetic information is incorporated into the PCMs depends on the method used. For the first group of methods, phylogenies are converted into variance-covariance matrices of traits following a given model of evolution such as Brownian motion (BM). For the second group of methods, phylogenies are converted into distance matrices that are subsequently transformed into Euclidean distances to perform principal coordinate analyses. Here, we show that simply taking the elementwise square root of a distance matrix extracted from a phylogenetic tree ensures having a Euclidean distance matrix. This is true for any type of distances between species (patristic or nodal) and also for trees harboring multifurcating nodes. Moreover, we illustrate that this simple transformation using the square root imposes less geometric distortion than more complex transformations classically used in the literature such as the Cailliez method. Given the Euclidean nature of the elementwise square root of phylogenetic distance matrices, the positive semidefinitiveness of the phylogenetic variance-covariance matrix of a trait following a BM model, or related models of trait evolution, can be established. In that way, we build a bridge between the two groups of statistical methods widely used in comparative analysis. These results should be of great interest for ecologists and evolutionary biologists performing statistical analyses incorporating phylogenies.     

ACCURACY OF ESTIMATED PHYLOGENIES: EFFECTS OF TREE TOPOLOGY AND EVOLUTIONARY MODEL

A simulation study was carried out to investigate the relative importance of tree topology (both balance and stemminess), evolutionary rates (constant, varying among characters, and varying among lineages), and evolutionary models in determining the accuracy with which phylogenetic trees can be estimated. The three evolutionary context models were phyletic (characters can change at each simulated time step), speciational (changes are possible only at the time of speciation into two daughter lineages), and punctuational (changes occur at the time of speciation but only in one of the daughter lineages). UPGMA clustering and maximum parsimony (“Wagner trees”) methods for estimating phylogenies were compared. All trees were based on eight recent OTUs. The three evolutionary context models were found to have the largest influence on accuracy of estimates by both methods. The next most important effect was that of the stemminess × context interaction. Maximum parsimony and UPGMA performed worst under the punctuational models. Under the phyletic model, trees with high stemminess values could be estimated more accurately and balanced trees were slightly easier to estimate than unbalanced ones. Overall, maximum parsimony yielded more accurate trees than UPGMA—but that was expected for these simulations since many more characters than OTUs were used. Our results suggest that the great majority of estimated phylogenetic trees are likely to be quite inaccurate; they underscore the inappropriateness of characterizing current phylogenetic methods as being for reconstruction rather than for estimation.     

Assessing the performance of single-copy genes for recovering robust phylogenies

Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.     

No substitute for real data: A cautionary note on the use of phylogenies from birth–death polytomy resolvers for downstream comparative analyses

The statistical estimation of phylogenies is always associated with uncertainty, and accommodating this uncertainty is an important component of modern phylogenetic comparative analysis. The birth–death polytomy resolver is a method of accounting for phylogenetic uncertainty that places missing (unsampled) taxa onto phylogenetic trees, using taxonomic information alone. Recent studies of birds and mammals have used this approach to generate pseudoposterior distributions of phylogenetic trees that are complete at the species level, even in the absence of genetic data for many species. Many researchers have used these distributions of phylogenies for downstream evolutionary analyses that involve inferences on phenotypic evolution, geography, and community assembly. I demonstrate that the use of phylogenies constructed in this fashion is inappropriate for many questions involving traits. Because species are placed on trees at random with respect to trait values, the birth–death polytomy resolver breaks down natural patterns of trait phylogenetic structure. Inferences based on these trees are predictably and often drastically biased in a direction that depends on the underlying (true) pattern of phylogenetic structure in traits. I illustrate the severity of the phenomenon for both continuous and discrete traits using examples from a global bird phylogeny.     

Phylotranscriptomics of Swertiinae (Gentianaceae) reveals that key floral traits are not phylogenetically correlated

Establishing how lineages with similar traits are phylogenetically related remains critical for understanding the origin of biodiversity on Earth. Floral traits in plants are widely used to explore phylogenetic relationships and to delineate taxonomic groups. The subtribe Swertiinae (Gentianaceae) comprises more than 350 species with high floral diversity ranging from rotate to tubular corollas and possessing diverse nectaries. Here we performed phylogenetic analysis of 60 species from all 15 genera of the subtribe Swertiinae sensu Ho and Liu, representing the range of floral diversity, using data from the nuclear and plastid genomes. Extensive topological conflicts were present between the nuclear and plastome trees. Three of the 15 genera represented by multiple species are polyphyletic in both trees. Key floral traits including corolla type, absence or presence of lobe scales, nectary type, nectary position, and stigma type are randomly distributed in the nuclear and plastome trees without phylogenetic correlation. We also revealed the likely ancient hybrid origin of one large clade comprising 10 genera with diverse floral traits. These results highlight the complex evolutionary history of this subtribe. The phylogenies constructed here provide a basic framework for further exploring the ecological and genetic mechanisms underlying both species diversification and floral diversity.     

Diversification and biogeographical history of Neotropical plethodontid salamanders

The Neotropical bolitoglossine salamanders represent an impressive adaptive radiation, comprising roughly 40% of global salamander species diversity. Despite decades of morphological studies and molecular work, a robust multilocus phylogenetic hypothesis based on DNA sequence data is lacking for the group. We estimated species trees based on multilocus nuclear and mitochondrial data for all major lineages within the bolitoglossines, and used our new phylogenetic hypothesis to test traditional biogeographical scenarios and hypotheses of morphological evolution in the group. In contrast to previous phylogenies, our results place all Central American endemic genera in a single clade and suggest that Central America played a critical role in the early biogeographical history of the group. The large, predominantly Mexican genus Pseudoeurycea is paraphyletic, and analyses of the nuclear data place two lineages of Pseudoeurycea as the sister group of Bolitoglossa. Our phylogeny reveals extensive homoplasy in morphological characters, which may be the result of truncation or alteration of a shared developmental trajectory. We used our phylogenetic results to revise the taxonomy of the genus Pseudoeurycea. © 2015 The Linnean Society of London     

The dynamics of preferential host switching: Host phylogeny as a key predictor of parasite distribution*

Parasites often jump to and become established in a new host species. There is much evidence that the probability of such host shifts decreases with increasing phylogenetic distance between donor and recipient hosts, but the consequences of such preferential host switching remain little explored. We develop a computational model to investigate the dynamics of parasite host shifts in the presence of this phylogenetic distance effect. In this model, a clade of parasites evolves on an evolving clade of host species where parasites can cospeciate with their hosts, switch to new hosts, speciate within hosts or become extinct. Our model predicts that host phylogenies are major determinants of parasite distributions across trees. In particular, we predict that trees consisting of few large clades of host species and those with fast species turnover should harbor more parasites than trees with many small clades and those that diversify more slowly. Within trees, large clades are predicted to exhibit a higher fraction of infected species than small clades. We discuss our results in the light of recent cophylogenetic studies in a wide range of host–parasite systems.     

When can host shifts produce congruent host and parasite phylogenies? A simulation approach

Congruence between host and parasite phylogenies is often taken as evidence for cospeciation. However, 'pseudocospeciation', resulting from host-switches followed by parasite speciation, may also generate congruent trees. To investigate this process and the conditions favouring its appearance, we here simulated the adaptive radiation of a parasite onto a new range of hosts. A very high congruence between the host tree and the resulting parasite trees was obtained when parasites switched between closely related hosts. Setting a shorter time lag for speciation after switches between distantly related hosts further increased the degree of congruence. The shape of the host tree, however, had a strong impact, as no congruence could be obtained when starting with highly unbalanced host trees. The strong congruences obtained were erroneously interpreted as the result of cospeciations by commonly used phylogenetic software packages despite the fact that all speciations resulted from host-switches in our model. These results highlight the importance of estimating the age of nodes in host and parasite phylogenies when testing for cospeciation and also demonstrate that the results obtained with software packages simulating evolutionary events must be interpreted with caution.     

TESTING ADAPTIVE RADIATION AND KEY INNOVATION HYPOTHESES IN SPIDERS

We combine statistical and phylogenetic approaches to test the hypothesis that adaptive radiation and key innovation have contributed to the diversity of the order Araneae. The number of unbalanced araneid clades (those whose species numbers differ by 90% or more) exceeds the number predicted by a null Markovian model. The current phylogeny of spider families contains 74 bifurcating nodes, of which 31 are unbalanced. As this is significantly more than the 14.8 expected unbalanced nodes, some of the diversity within the Araneae can be attributed to some deterministic cause (e.g., adaptive radiation). One of the more highly unbalanced (97%) bifurcations divides the orb-weaving spiders into the Deinopoidea and the larger Araneoidea. A simple statistical model shows that the inequality in diversity between the Deinopoidea and the Araneoidea is significant, and that it is associated with the replacement of primitive cribellar capture thread by viscous adhesive thread and a change from a horizontal to a vertical orb-web orientation. These changes improve an orb-web's ability to intercept and retain prey and expand the adaptive zone that orb-weaving spiders can occupy and are, therefore, considered to be “key innovations.”     

Whole-tree methods for detecting differential diversification rates

Prolific cladogenesis, adaptive radiation, species selection, key innovations, and mass extinctions are a few examples of biological phenomena that lead to differential diversification among lineages. Central to the study of differential diversification rates is the ability to distinguish chance variation from that which requires deterministic explanation. To detect diversification rate variation among lineages, we propose a number of methods that incorporate information on the topological distribution of species diversity from all internal nodes of a phylogenetic tree. These whole-tree methods (M(Pi), M(Sigma), and M(R)) are explicitly connected to a null model of random diversification--the equal-rates Markov (ERM) random branching model--and an alternative model of differential diversification: M(Pi) is based on the product of individual nodal ERM probabilities; M(Sigma) is based on the sum of individual nodal ERM probabilities, and M(R) is based on a transformation of ERM probabilities that corresponds to a formalized system that orders trees by their relative symmetry. These methods have been implemented in a freely available computer program, SYMMETREE, to detect clades with variable diversification rates, thereby allowing the study of biological processes correlated with and possibly causal to shifts in diversification rate. Application of these methods to several published phylogenies demonstrates their ability to contend with relatively large, incompletely resolved trees. These topology-based methods do not require estimates of relative branch lengths, which should facilitate the analysis of phylogenies, such as supertrees, for which such data are unreliable or unavailable.     

PHYLOGENETIC HYPOTHESES UNDER THE ASSUMPTION OF NEUTRAL QUANTITATIVE-GENETIC VARIATION

There are many situations in which the only available characters for reconstructing phylogenies are morphological. Those traits that are subject only to the forces of mutation and random genetic drift can be used to obtain unbiased estimates of phylogenetic relationships. However, the accurate recovery of a phylogeny from information on neutral characters requires the procurement of data for a large number of independent traits, individuals, and populations. Phylogenetic trees fit to data from more than five species will almost always contain topological errors, even with very large data sets. The population-genetic consequences of the neutral model are reviewed, and some statistical methods for testing whether the diversification of a phylogeny is compatible with such a model are outlined. The theory is then applied to a very large data set on cranial morphology in modern man. The results are consistent with the hypothesis that interracial differences in human skull dimensions are a simple consequence of random drift and mutation.     

The Use (and Misuse) of Phylogenetic Trees in Comparative Behavioral Analyses

Phylogenetic comparative methods play a critical role in our understanding of the adaptive origin of primate behaviors. To incorporate evolutionary history directly into comparative behavioral research, behavioral ecologists rely on strong, well-resolved phylogenetic trees. Phylogenies provide the framework on which behaviors can be compared and homologies can be distinguished from similarities due to convergent or parallel evolution. Phylogenetic reconstructions are also of critical importance when inferring the ancestral state of behavioral patterns and when suggesting the evolutionary changes that behavior has undergone. Improvements in genome sequencing technologies have increased the amount of data available to researchers. Recently, several primate phylogenetic studies have used multiple loci to produce robust phylogenetic trees that include hundreds of primate species. These trees are now commonly used in comparative analyses and there is a perception that we have a complete picture of the primate tree. But how confident can we be in those phylogenies? And how reliable are comparative analyses based on such trees? Herein, we argue that even recent molecular phylogenies should be treated cautiously because they rely on many assumptions and have many shortcomings. Most phylogenetic studies do not model gene tree diversity and can produce misleading results, such as strong support for an incorrect species tree, especially in the case of rapid and recent radiations. We discuss implications that incorrect phylogenies can have for reconstructing the evolution of primate behaviors and we urge primatologists to be aware of the current limitations of phylogenetic reconstructions when applying phylogenetic comparative methods.     

Phylogenetic Resolution and Quantifying the Phylogenetic Diversity and Dispersion of Communities

Conservation biologists and community ecologists have increasingly begun to quantify the phylogenetic diversity and phylogenetic dispersion in species assemblages. In some instances, the phylogenetic trees used for such analyses are fully bifurcating, but in many cases the phylogenies being used contain unresolved nodes (i.e. polytomies). The lack of phylogenetic resolution in such studies, while certainly not preferred, is likely to continue particularly for those analyzing diverse communities and datasets with hundreds to thousands of taxa. Thus it is imperative that we quantify potential biases and losses of statistical power in studies that use phylogenetic trees that are not completely resolved. The present study is designed to meet both of these goals by quantifying the phylogenetic diversity and dispersion of simulated communities using resolved and gradually ‘unresolved’ phylogenies. The results show that: (i) measures of community phylogenetic diversity and dispersion are generally more sensitive to loss of resolution basally in the phylogeny and less sensitive to loss of resolution terminally; and (ii) the loss of phylogenetic resolution generally causes false negative results rather than false positives.     

Inferring epidemic contact structure from phylogenetic trees

Contact structure is believed to have a large impact on epidemic spreading and consequently using networks to model such contact structure continues to gain interest in epidemiology. However, detailed knowledge of the exact contact structure underlying real epidemics is limited. Here we address the question whether the structure of the contact network leaves a detectable genetic fingerprint in the pathogen population. To this end we compare phylogenies generated by disease outbreaks in simulated populations with different types of contact networks. We find that the shape of these phylogenies strongly depends on contact structure. In particular, measures of tree imbalance allow us to quantify to what extent the contact structure underlying an epidemic deviates from a null model contact network and illustrate this in the case of random mixing. Using a phylogeny from the Swiss HIV epidemic, we show that this epidemic has a significantly more unbalanced tree than would be expected from random mixing.