Cospeciation between bacterial endosymbionts (Buchnera) and a recent radiation of aphids (Uroleucon) and pitfalls of testing for phylogenetic congruence

Abstract.— Previous studies of phylogenetic congruence between aphids and their symbiotic bacteria ( Buchnera ) supported long-term vertical transmission of symbionts. However, those studies were based on distantly related aphids and would not have revealed horizontal transfer of symbionts among closely related hosts. Aphid species of the genus Uroleucon are closely related phylogenetically and overlap in geographic ranges, habitats, and parasitoids. To examine support for congruence of phylogenies of Buchnera and Uroleucon , sequences from four mitochondrial, one nuclear, and one endosymbiont gene ( trpB ) were obtained. Congruence of phylogenies based on pooled aphid genes with phylogenies based on trpB was highly significant: Most nodes resolved by trpB corresponded to nodes resolved by the pooled aphid genes. Furthermore, no nodes were both inconsistent between the trees and strongly supported in both trees. Two kinds of analyses testing the null hypothesis of perfect congruence between pairwise combinations of datasets and tree topologies were performed: the Kishino-Hasegawa test and the likelihood-ratio test. Both tests indicated significant disagreement among most pairwise combinations of mitochondrial, nuclear, and symbiont datasets. Because rampant recombination among mitochondrial genomes of different aphid species is unlikely, inaccurate assumptions in the evolutionary models underlying these tests appear to be causing the hypothesis of a shared history to be incorrectly rejected. Moreover, trpB was more consistent with the aphid genes as a set than any single aphid gene was with the others, suggesting that the symbionts show the same phylogeny as the aphids. Overall, analyses support the interpretation that symbionts and aphids have undergone strict cospeciation, with no horizontal transmission of symbionts even among closely related, ecologically similar aphid hosts.     

Using phylogeographic analyses of gene trees to test species status and processes

A gene tree is an evolutionary reconstruction of the genealogical history of the genetic variation found in a sample of homologous genes or DNA regions that have experienced little or no recombination. Gene trees have the potential of straddling the interface between intra- and interspecific evolution. It is precisely at this interface that the process of speciation occurs, and gene trees can therefore be used as a powerful tool to probe this interface. One application is to infer species status. The cohesion species is defined as an evolutionary lineage or set of lineages with genetic exchangeability and/or ecological interchangeability. This species concept can be phrased in terms of null hypotheses that can be tested rigorously and objectively by using gene trees. First, an overlay of geography upon the gene tree is used to test the null hypothesis that the sample is from a single evolutionary lineage. This phase of testing can indicate that the sampled organisms are indeed from a single lineage and therefore a single cohesion species. In other cases, this null hypothesis is not rejected due to a lack of power or inadequate sampling. Alternatively, this null hypothesis can be rejected because two or more lineages are in the sample. The test can identify lineages even when hybridization and lineage sorting occur. Only when this null hypothesis is rejected is there the potential for more than one cohesion species. Although all cohesion species are evolutionary lineages, not all evolutionary lineages are cohesion species. Therefore, if the first null hypothesis is rejected, a second null hypothesis is tested that all lineages are genetically exchangeable and/or ecologically interchangeable. This second test is accomplished by direct contrasts of previously identified lineages or by overlaying reproductive and/or ecological data upon the gene tree and testing for significant transitions that are concordant with the previously identified lineages. Only when this second null hypothesis is rejected is a lineage elevated to the status of cohesion species. By using gene trees in this manner, species can be identified with objective, a priori criteria with an inference procedure that automatically yields much insight into the process of speciation. When one or more of the null hypotheses cannot be rejected, this procedure also provides specific guidance for future work that will be needed to judge species status.     

When does the incongruence length difference test fail?

This paper examines the efficiency of the incongruence length difference test (ILD) proposed by Farris et al. (1994) for assessing the incongruence between sets of characters. DNA sequences were simulated under various evolutionary conditions: (1) following symmetric or asymmetric trees, (2) with various mutation rates, (3) with constant or variable evolutionary rates along the branches, and (4) with different among-site substitution rates. We first compared two sets of sequences generated along the same tree and under the same evolutionary conditions. The probability of a Type-I error (wrongly rejecting the true hypothesis of congruence) was substantially below the standard 5% level of significance given by the ILD test; this finding indicates that the choice of the 5% level is rather conservative in this case. We then compared two data sets, still generated along the same tree, but under different evolutionary conditions (constant vs. variable evolutionary rate, homogeneity vs. heterogeneity rate of substitution). Under these conditions, the probability of rejecting the true hypothesis of congruence was greater than the 5% given by the ILD test and increased with the number of sites and the degree to which the tree was asymmetric. Finally, the comparison of the two data sets, simulated under contrasting tree structures (symmetric vs. asymmetric) but under the same evolutionary conditions, led us to reject the hypothesis of congruence, albeit weakly, particularly when the number of informative sites was low and among-site substitution rate heterogeneous. We conclude that the ILD test has only limited power to detect incongruence caused by differences in the evolutionary conditions or in the tree topology, except when numerous characters are present and the substitution rate is homogeneous from site to site.     

Cophylogeny of the Ficus microcosm

The various mutualistic and antagonistic symbioses between fig trees (Ficus: Moraceae) and chalcid wasps comprise a community in microcosm. Phylogenetic estimates of figs and fig wasps show general topological correspondence, making the microcosm a model system for cophylogeny. Incongruence between phylogenies from associated organisms can be reconciled through a combination of evolutionary events. Cophylogeny mapping reconciles phylogenies by embedding an associate tree into a host tree, finding the optimal combinations of events capable of explaining incongruence and evaluating the level of codivergence. This review addresses the results of cophylogeny analysis concerning Ficus and discusses the plausibility of different evolutionary events. Five different associations encompassing fig-pollinator, fig-parasite and pollinator-parasitoid interactions are reconciled. The method improves on previous comparisons by employing 'jungles' to provide an exhaustive and quantitative analysis of cophylogeny. A jungle is a mechanism for inferring host switches and obtaining all potentially optimal solutions to the reconciliation problem. The results support the consensus that figs codiverge significantly with pollinators but not non-pollinators. However, pollinators still appear to have switched between hosts in contradiction to the traditional model of faithful codivergence. This emphasises the growing realisation that evolutionary transitions in the microcosm are more flexible than previously thought and host specificity is necessary but not sufficient for codivergence. The importance of sampling strategy is emphasised by the influence of taxon set on the fig-pollinator and fig-parasite jungles. Spurious significant results for fig-parasite and fig-parasitoid jungles indicate that the choice of congruence measure influences significance; the total number of events required to reconcile two trees ('total cost') is not a good measure of congruence when switches cannot be realistically weighted.     

A cophylogenetic perspective of RNA-virus evolution

The extent to which viruses and their hosts codiverge remains an open question, given that numerous cases of both "cospeciation" and horizontal switching have recently been documented. DNA viruses that form persistent infections are thought to be the most likely candidates for phylogenetic congruence. Phylogenetic reconciliation analysis was used to compare established phylogenies for four RNA viruses and their hosts. The analysis employs a cophylogeny mapping technique, implemented in TreeMap v2.0, to find the most parsimonious combinations of evolutionary events able to reconcile any incongruence. This technique is guaranteed to recover all potentially optimal solutions to the reconciled tree and specifically tests the null hypothesis that an associate phylogeny is no more congruent with a host phylogeny than would be a random tree with the same taxon set. Phylogenies for Hantavirus, Spumavirus, and avian sarcoma leukosis virus were found to be significantly similar to their host trees, whereas Lyssavirus and Arenavirus displayed no significant congruence. These results demonstrate that RNA viruses are able to form stable associations with their hosts over evolutionary time scales and that the details of such associations are consistent with persistent infection being a necessary but not sufficient precondition.     

NULL MODELS FOR THE NUMBER OF EVOLUTIONARY STEPS IN A CHARACTER ON A PHYLOGENETIC TREE

Random trees and random characters can be used in null models for testing phylogenetic hypothesis. We consider three interpretations of random trees: first, that trees are selected from the set of all possible trees with equal probability; second, that trees are formed by random speciation or coalescence (equivalent); and third, that trees are formed by a series of random partitions of the taxa. We consider two interpretations of random characters: first, that the number of taxa with each state is held constant, but the states are randomly reshuffled among the taxa; and second, that the probability each taxon is assigned a particular state is constant from one taxon to the next. Under null models representing various combinations of randomizations of trees and characters, exact recursion equations are given to calculate the probability distribution of the number of character state changes required by a phylogenetic tree. Possible applications of these probability distributions are discussed. They can be used, for example, to test for a panmictic population structure within a species or to test phylogenetic inertia in a character's evolution. Whether and how a null model incorporates tree randomness makes little difference to the probability distribution in many but not all circumstances. The null model's sense of character randomness appears more critical. The difficult issue of choosing a null model is discussed.     

Predicting protein-protein interaction by searching evolutionary tree automorphism space

MOTIVATION: Uncovering the protein-protein interaction network is a fundamental step in the quest to understand the molecular machinery of a cell. This motivates the search for efficient computational methods for predicting such interactions. Among the available predictors are those that are based on the co-evolution hypothesis "evolutionary trees of protein families (that are known to interact) are expected to have similar topologies". Many of these methods are limited by the fact that they can handle only a small number of protein sequences. Also, details on evolutionary tree topology are missing as they use similarity matrices in lieu of the trees. RESULTS: We introduce MORPH, a new algorithm for predicting protein interaction partners between members of two protein families that are known to interact. Our approach can also be seen as a new method for searching the best superposition of the corresponding evolutionary trees based on tree automorphism group. We discuss relevant facts related to the predictability of protein-protein interaction based on their co-evolution. When compared with related computational approaches, our method reduces the search space by approximately 3 x 10(5)-fold and at the same time increases the accuracy of predicting correct binding partners.     

A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between Type 2 Diabetes and Prostate Cancer

There is increasing evidence that pleiotropy, the association of multiple traits with the same genetic variants/loci, is a very common phenomenon. Cross-phenotype association tests are often used to jointly analyze multiple traits from a genome-wide association study (GWAS). The underlying methods, however, are often designed to test the global null hypothesis that there is no association of a genetic variant with any of the traits, the rejection of which does not implicate pleiotropy. In this article, we propose a new statistical approach, PLACO, for specifically detecting pleiotropic loci between two traits by considering an underlying composite null hypothesis that a variant is associated with none or only one of the traits. We propose testing the null hypothesis based on the product of the Z-statistics of the genetic variants across two studies and derive a null distribution of the test statistic in the form of a mixture distribution that allows for fractions of variants to be associated with none or only one of the traits. We borrow approaches from the statistical literature on mediation analysis that allow asymptotic approximation of the null distribution avoiding estimation of nuisance parameters related to mixture proportions and variance components. Simulation studies demonstrate that the proposed method can maintain type I error and can achieve major power gain over alternative simpler methods that are typically used for testing pleiotropy. PLACO allows correlation in summary statistics between studies that may arise due to sharing of controls between disease traits. Application of PLACO to publicly available summary data from two large case-control GWAS of Type 2 Diabetes and of Prostate Cancer implicated a number of novel shared genetic regions: 3q23 (ZBTB38), 6q25.3 (RGS17), 9p22.1 (HAUS6), 9p13.3 (UBAP2), 11p11.2 (RAPSN), 14q12 (AKAP6), 15q15 (KNL1) and 18q23 (ZNF236).     

Differences among Yanomama Indian villages: Do the patterns of allele frequencies,anthropometrics and map locations correspond?

In order to determine the degree of correspondence between sets of multivariate observations based on different kinds of traits, two new methods, derived from fundamentally different notions of “correspondence,” are adopted here and compared. Using networks or trees to represent contemporary relationships, the first method tests the similarity of the cluster or hierarchic structures implicit in two sets of data. The second approach tests the departure from perfect geometric congruence or superimposability. Computer simulation was used to generate the distributions needed for significance tests under the null hypothesis. By the first technique, we find significant correspondence among the cluster structures for geographic, allele frequency, and anthropometric data on 19 Yanomama Indian villages. The results are similar and more precise for a subset consisting of seven villages. Some of these results differ from the conclusions which would be reached with the conventional correlations based upon entries in distance tables. The direct test of congruence, used only for the data on the subset of seven villages, gives results which differ substantially from those based on cluster-structure. There are, however, similarities between the measure of congruence and the simple correlations based on entries in the distance tables. The significant correspondences observed call for some explanation. Cultural and demographic features determine the particular non-random allocation of individuals to village fragments when a village splits. These social phenomena are invoked in tentative explanation of the agreement among historical, biological, and geographic relationships of villages.     

Genomic rearrangements in trypanosomatids: an alternative to the "one gene" evolutionary hypotheses?

Most molecular trees of trypanosomatids are based on point mutations within DNA sequences. In contrast, there are very few evolutionary studies considering DNA (re) arrangement as genetic characters. Waiting for the completion of the various parasite genome projects, first information may already be obtained from chromosome size-polymorphism, using the appropriate algorithms for data processing. Three illustrative models are presented here. First, the case of Leishmania (Viannia) braziliensis/L. (V.) peruviana is described. Thanks to a fast evolution rate (due essentially to amplification/deletion of tandemly repeated genes), molecular karyotyping seems particularly appropriate for studying recent evolutionary divergence, including eco-geographical diversification. Secondly, karyotype evolution is considered at the level of whole genus Leishmania. Despite the fast chromosome evolution rate, there is qualitative congruence with MLEE- and RAPD-based evolutionary hypotheses. Significant differences may be observed between major lineages, likely corresponding to major and less frequent rearrangements (fusion/fission, translocation). Thirdly, comparison is made with Trypanosoma cruzi. Again congruence is observed with other hypotheses and major lineages are delineated by significant chromosome rearrangements. The level of karyotype polymorphism within that "species" is similar to the one observed in "genus" Leishmania. The relativity of the species concept among these two groups of parasites is discussed.     

TOPD/FMTS: a new software to compare phylogenetic trees

SUMMARY: TOPD/FMTS has been developed to evaluate similarities and differences between phylogenetic trees. The software implements several new algorithms (including the Disagree method that returns the taxa, that disagree between two trees and the Nodal method that compares two trees using nodal information) and several previously described methods (such as the Partition method, Triplets or Quartets) to compare phylogenetic trees. One of the novelties of this software is that the FMTS (From Multiple to Single) program allows the comparison of trees that contain both orthologs and paralogs. Each option is also complemented with a randomization analysis to test the null hypothesis that the similarity between two trees is not better than chance expectation. AVAILABILITY: The Perl source code of TOPD/FMTS is available at http://genomes.urv.es/topd.     

Testing the generalized neutrality hypothesis

Various tests of the hypothesis of selective neutrality based on gene frequency are now available. These tests take as null hypothesis the concept of “strict neutrality”: all new mutants are required to be selectively identical to each other. For evolutionary questions, however, (as opposed to those of genetic polymorphism), a wider null hypothesis might be of interest. Since deleterious alleles have essentially no evolutionary importance, one might wish to test the null hypothesis that only neutral or deleterious mutations occur. The principal alternative to this hypothesis is that there exists heterotic selection of some form for some alleles tending to maintain a level of genetic polymorphism higher than that under neutrality. In this paper an assessment is made of the usefulness of a test of strict neutrality first proposed by this author (Ewens, 1972) as a test of null hypothesis of “generalized neutrality,” i.e. that only neutral or deleterious alleles occur. At the same time some remarks will be made about estimation of the fundamental parameter θ defining these processes.     

Macroecology of live-bearing in fishes: latitudinal and depth range comparisons with egg-laying relatives

The habitat heterogeneity hypothesis states that an increase in habitat heterogeneity leads to an increase in species diversity. We tested this hypothesis for a community of small mammals in the semiarid, sand-shinnery-oak ecosystem of the southwestern United States. We used indices of differentiation diversity to quantify differences between two habitat types (blowouts in a sand-shinnery-oak matrix) in terms of species diversity. The Wilson-Shmida index (β_T) considers species composition only, whereas the Morisita-Horn index (C_mH) also takes species abundances into account. We constructed null models to test the hypothesis that differentiation diversity between habitat types is greater than that produced by stochastic processes. Two models were constructed, one based on the random placement of species and one based on the random placement of individuals. No evidence supported the hypothesis that habitat heterogeneity enhances diversity of a landscape by increasing the number of species in an area. Indeed, paired habitats were more similar than chance alone would dictate in terms of species identities. In contrast, habitat heterogeneity affects diversity by significantly altering the relative proportions of species in contrasting habitat types. Because seeds differentially accumulate at the interface between blowouts and matrix, the high productivity of the edge may actually homogenize habitat types in terms of species richness. Nonetheless, blowouts might best be considered to be microhabitats which enhance or complement the value of the matrix even though the species which use either habitat type are identical.     

Assessing progress in systematics with continuous jackknife function analysis

Systematists expect their hypotheses to be asymptotically precise. As the number of phylogenetically informative characters for a set of taxa increases, the relationships implied should stabilize on some topology. If true, this increasing stability should clearly manifest itself if an index of congruence is plotted against the accumulating number of characters. Continuous jackknife function (CJF) analysis is a new graphical method that portrays the extent to which available data converge on a specified phylogenetic hypothesis, the reference tree. The method removes characters with increasing probability, analyzes the rarefied data matrices phylogenetically, and scores the clades shared between each of the resulting trees and the reference tree. As more characters are removed, the number of shared clades must decrease, but the rate of decrease will depend on how decisively the data support the reference tree. Curves for stable phylogenies are clearly asymptotic with nearly 100% congruence for a substantial part of the curve. Less stable phylogenies lose congruent nodes quickly as characters are excluded, resulting in a more linear or even a sigmoidal relationship. Curves can be interpreted as predictors of whether the addition of new data of the same type is likely to alter the hypothesis under test. Continuous jackknife function analysis makes statistical assumptions about the collection of character data. To the extent that CJF curves are sensitive to violations of unbiased character collection, they will be misleading as predictors. Convergence of data on a reference tree does not guarantee historical accuracy, but it does predict that the accumulation of further data under the sampling model will not lead to rapid changes in the hypothesis.     

Geographical and ecological patterns of range size in North American trees

We analysed the degree to which spatial, ecological and phylogenetic factors influenced geographic gradients in range size for all trees native to North America. We compared observed range sizes, measured in either one (latitudinal or longitudinal extent) or two dimensions (range area), with three null models that respected constraints imposed by the geometry of the continent; we tested whether climatic gradients could account for increasing range size with latitude – a Rapoport effect – in North American trees; and whether variation in range size was associated with phylogeny or ecological characteristics of species. We documented an excess of species with small ranges and a lack of widely distributed species compared with null expectations both with and without environmental constraints. We found evidence for a Rapoport effect in North American trees, at different taxonomic levels and for different groups defined by ecological characteristics. This pattern is not an outcome constrained by continental shape and does not depend on the metric for range size, but it is consistent with the climatic variability hypothesis proposed to explain the Rapoport pattern. Finally we showed that trees east of the Rocky Mountains have larger ranges on average than trees to the west or in Mexico and that the proportion of deciduous compared to evergreen species increases with range size. Variation in range size among North American trees is not spatially random, and has a phylogenetic and ecological signal. Consistent with the climatic variability hypothesis, range size of North American trees increases with increasing variability in annual temperature range at higher latitudes.     

A phylogenetic test of multiple proposals for the origins of the East Indies coral reef biota

Aim  We analysed data on Indo-Pacific coral reef taxa to test four mechanisms proposed for the origins of the biodiversity centre in the East Indies Triangle: (1) the centre of origin hypothesis, (2) the centre of accumulation hypothesis, (3) the centre of overlap hypothesis, and (4) the bioaccumulation hypothesis.
Location  The Indian and western Pacific oceans.
Methods  The data set consisted of eight clades of fishes, four of corals and one of molluscs, consisting of a total of 95 species, with 72 informative nodes, in 29 areas. We used phylogenetic analysis for comparing trees (PACT) to reconstruct reticulate area relationships.
Results  An analysis using PACT produced very little congruence between area relationships among the histories of the taxa examined. Only two sets of two taxa each showed substantial congruence; several other taxa showed partial congruence.
Main conclusions  Despite numerous 'explanatory' claims for the incredible biodiversity of the marine macrofauna of the East Indies Triangle, the patterns obtained here did not accord with any of the hypotheses proposed. It is possible that, with the addition of more taxa, additional patterns would emerge. Much more systematic work within the East Indies is required to resolve this problem.     

How Bayes tests of molecular phylogenies compare with frequentist approaches

MOTIVATION: The desire to compare molecular phylogenies has stimulated the design of numerous tests. Most of these tests are formulated in a frequentist framework, and it is not known how they compare with Bayes procedures. I propose here two new Bayes tests that either compare pairs of trees (Bayes hypothesis test, BHT), or test each tree against an average of the trees included in the analysis (Bayes significance test, BST). RESULTS: The algorithm, based on a standard Metropolis-Hastings sampler, integrates nuisance parameters out and estimates the probability of the data under each topology. These quantities are used to estimate Bayes factors for composite vs. composite hypotheses. Based on two data sets, the BHT and BST are shown to construct similar confidence sets to the bootstrap and the Shimodaira Hasegawa test, respectively. This suggests that the known difference among previous tests is mainly due to the null hypothesis considered.     

Congruence of Morphological and Molecular Phylogenies

When phylogenetic trees constructed from morphological and molecular evidence disagree (i.e. are incongruent) it has been suggested that the differences are spurious or that the molecular results should be preferred a priori. Comparing trees can increase confidence (congruence), or demonstrate that at least one tree is incorrect (incongruence). Statistical analyses of 181 molecular and 49 morphological trees shows that incongruence is greater between than within the morphological and molecular partitions, and this difference is significant for the molecular partition. Because the level of incongruence between a pair of trees gives a minimum bound on how much error is present in the two trees, our results indicate that the level of error may be underestimated by congruence within partitions. Thus comparisons between morphological and molecular trees are particularly useful for detecting this incongruence (spurious or otherwise). Molecular trees have higher average congruence than morphological trees, but the difference is not significant, and both within- and between-partition incongruence is much lower than expected by chance alone. Our results suggest that both molecular and morphological trees are, in general, useful approximations of a common underlying phylogeny and thus, when molecules and morphology clash, molecular phylogenies should not be considered more reliable a priori.     

Dissimilar rates in molecular evolution

In this work we present an evolutionary tree based on the differences in the physico-chemical properties involved in amino acid substitutions, instead of considering, for its construction, only the number of changes between species. Phylogenetic trees were constructed from the differences in bulkiness, refractivity index, hydrophobicity, polarity and optical rotation of 9 vertebrate calcitonins. A correlation of the form y = a xb was found between the number of changes (x) and the differences in any given physico-chemical property (y). This correlation implies that the evolutionary time can not be evaluated directly from the number of changes between species.     

Dissimilar rates in molecular evolution

In this work we present an evolutionary tree based on the differences in the physico-chemical properties involved in amino acid substitutions, instead of considering, for its construction, only the number of changes between species. Phylogenetic trees were constructed from the differences in bulkiness, refractivity index, hydrophobicity, polarity and optical rotation of 9 vertebrate calcitonins. A correlation of the form y=a x^b was found between the number of changes (x) and the differences in any given physico-chemical property (y). This correlation implies that the evolutionary time can not be evaluated directly from the number of changes between species.