Bootstrap confidence levels for phylogenetic trees

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783–791] suggested the use of the bootstrap to answer this question. Felsenstein’s method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein’s method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

Matchings and phylogenetic trees

This paper presents a natural coordinate system for phylogenetic trees using a correspondence with the set of perfect matchings in the complete graph. This correspondence produces a distance between phylogenetic trees, and a way of enumerating all trees in a minimal step order. It is useful in randomized algorithms because it enables moves on the space of trees that make random optimization strategies “mix” quickly. It also promises a generalization to intermediary trees when data are not decisive as to their choice of tree, and a new way of constructing Bayesian priors on tree space.

Bootstrap confidence levels for phylogenetic trees.

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783-791] suggested the use of the bootstrap to answer this question. Felsenstein's method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein's method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

Phylogenetic relationships between the Acantharea and the Polycystinea: A molecular perspective on Haeckel’s Radiolaria

Polycystine radiolaria are among few protistan groups that possess a comprehensive fossil record available for study by micropaleontologists. The Polycystinea and the Acantharea, whose skeletons do not become fossilized, were once members of the class “Radiolaria” (“Radiolaria” sensu lato: Polycystinea, Phaeodarea, and Acantharea) originally proposed by Haeckel but are now included in the superclass Actinopoda. Phylogenetic relationships within this superclass remain largely enigmatic. We investigated the evolutionary relationship of the Acantharea and the Polycystinea to other protists using phylogenetic analyses of 16S-like ribosomal RNA (rRNA) coding regions. We circumvented the need to culture these organisms by collecting and maintaining reproductive stages that contain many copies of their genomic DNA. This strategy facilitated extraction of genomic DNA and its purification from symbiont and prey DNA. Phylogenetic trees inferred from comparisons of 16S-like coding regions do not support a shared history between the Acantharea and the Polycystinea. However, the monophyly of the Acantharea and the separate monophyly of the Polycystinea (Spumellarida) are well supported by our molecular-based trees. The acantharian lineage branches among crown organisms whereas the polycystine lineage diverges before the radiation of the crown groups. We conclude that the Actinopoda does not represent a monophyletic evolutionary assemblage and recommend that this taxonomic designation be discarded.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Phylogenetic trees for groups of closely related species often have different topologies, depending on the genes used. One explanation for the discordant topologies is the persistence of polymorphisms through the speciation phase, followed by differential fixation of alleles in the resulting species. The existence of transspecies polymorphisms has been documented for alleles maintained by balancing selection but not for neutral alleles. In the present study, transspecific persistence of neutral polymorphisms was tested in the endemic haplochromine species flock of Lake Victoria cichlid fish. Putative noncoding region polymorphisms were identified at four randomly selected nuclear loci and tested on a collection of 12 Lake Victoria species and their putative riverine ancestors. At all loci, the same polymorphism was found to be present in nearly all the tested species, both lacustrine and riverine. Different polymorphisms at these loci were found in cichlids of other East African lakes (Malawi and Tanganyika). The Lake Victoria polymorphisms must have therefore arisen after the flocks now inhabiting the three great lakes diverged from one another, but before the riverine ancestors of the Lake Victoria flock colonized the Lake. Calculations based on the mtDNA clock suggest that the polymorphisms have persisted for about 1.4 million years. To maintain neutral polymorphisms for such a long time, the population size must have remained large throughout the entire period.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Plasmon analyses of Triticum (wheat) and Aegilops: PCR–single-strand conformational polymorphism (PCR-SSCP) analyses of organellar DNAs

To investigate phylogenetic relationships among plasmons in Triticum and Aegilops, PCR–single-strand conformational polymorphism (PCR-SSCP) analyses were made of 14.0-kb chloroplast (ct) and 13.7-kb mitochondrial (mt)DNA regions that were isolated from 46 alloplasmic wheat lines and one euplasmic line. These plasmons represent 31 species of the two genera. The ct and mtDNA regions included 10 and 9 structural genes, respectively. A total of 177 bands were detected, of which 40.6% were variable. The proportion of variable bands in ctDNA (51.1%) was higher than that of mtDNA (28.9%). The phylogenetic trees of plasmons, derived by two different models, indicate a common picture of plasmon divergence in the two genera and suggest three major groups of plasmons (Einkorn, Triticum, and Aegilops). Because of uniparental plasmon transmission, the maternal parents of all but one polyploid species were identified. Only one Aegilops species, Ae. speltoides, was included in the Triticum group, suggesting that this species is the plasmon and B and G genome donor of all polyploid wheats. ctDNA variations were more intimately correlated with vegetative characters, whereas mtDNA variations were more closely correlated with reproductive characters. Plasmon divergence among the diploids of the two genera largely paralleled genome divergence. The relative times of origin of the polyploid species were inferred from genetic distances from their putative maternal parents.

Gene content phylogeny of herpesviruses

Clusters of orthologous groups [COGs; Tatusov, R. L., Koonin, E. V. & Lipman, D. J. (1997) Science 278, 631–637] were identified for a set of 13 completely sequenced herpesviruses. Each COG represented a family of gene products conserved across several herpes genomes. These families were defined without using an arbitrary threshold criterion based on sequence similarity. The COG technique was modified so that variable stringency in COG construction was possible. High stringencies identify a core set of highly conserved genes. Varying COG stringency reveals differences in the degree of conservation between functional classes of genes. The COG data were used to construct whole-genome phylogenetic trees based on gene content. These trees agree well with trees based on other methods and are robust when tested by bootstrap analysis. The COG data also were used to construct a reciprocal tree that clustered genes with similar phylogenetic profiles. This clustering may give clues to genes with related functions or with related histories of acquisition and loss during herpesvirus evolution.

A phylogenomic approach to microbial evolution

To study the origin and evolution of biochemical pathways in microorganisms, we have developed methods and software for automatic, large-scale reconstructions of phylogenetic relationships. We define the complete set of phylogenetic trees derived from the proteome of an organism as the phylome and introduce the term phylogenetic connection as a concept that describes the relative relationships between taxa in a tree. A query system has been incorporated into the system so as to allow searches for defined categories of trees within the phylome. As a complement, we have developed the pyphy system for visualising the results of complex queries on phylogenetic connections, genomic locations and functional assignments in a graphical format. Our phylogenomics approach, which links phylogenetic information to the flow of biochemical pathways within and among microbial species, has been used to examine more than 8000 phylogenetic trees from seven microbial genomes. The results have revealed a rich web of phylogenetic connections. However, the separation of Bacteria and Archaea into two separate domains remains robust.

The RDP-II (Ribosomal Database Project)

The Ribosomal Database Project (RDP-II), previously described by Maidak et al. [Nucleic Acids Res. (2000), 28, 173–174], continued during the past year to add new rRNA sequences to the aligned data and to improve the analysis commands. Release 8.0 (June 1, 2000) consisted of 16 277 aligned prokaryotic small subunit (SSU) rRNA sequences while the number of eukaryotic and mitochondrial SSU rRNA sequences in aligned form remained at 2055 and 1503, respectively. The number of prokaryotic SSU rRNA sequences more than doubled from the previous release 14 months earlier, and ~75% are longer than 899 bp. An RDP-II mirror site in Japan is now available (http://wdcm.nig.ac.jp/RDP/html/index.html). RDP-II provides aligned and annotated rRNA sequences, derived phylogenetic trees and taxonomic hierarchies, and analysis services through its WWW server (http://rdp.cme.msu.edu/). Analysis services include rRNA probe checking, approximate phylogenetic placement of user sequences, screening user sequences for possible chimeric rRNA sequences, automated alignment, production of similarity matrices and services to plan and analyze terminal restriction fragment polymorphism experiments. The RDP-II email address for questions and comments has been changed from curator@cme.msu.edu to rdpstaff@msu.edu.

PALI—a database of Phylogeny and ALIgnment of homologous protein structures

PALI (release 1.2) contains three-dimensional (3-D) structure-dependent sequence alignments as well as structure-based phylogenetic trees of homologous protein domains in various families. The data set of homologous protein structures has been derived by consulting the SCOP database (release 1.50) and the data set comprises 604 families of homologous proteins involving 2739 protein domain structures with each family made up of at least two members. Each member in a family has been structurally aligned with every other member in the same family (pairwise alignment) and all the members in the family are also aligned using simultaneous super­position (multiple alignment). The structural alignments are performed largely automatically, with manual interventions especially in the cases of distantly related proteins, using the program STAMP (version 4.2). Every family is also associated with two dendrograms, calculated using PHYLIP (version 3.5), one based on a structural dissimilarity metric defined for every pairwise alignment and the other based on similarity of topologically equivalent residues. These dendrograms enable easy comparison of sequence and structure-based relationships among the members in a family. Structure-based alignments with the details of structural and sequence similarities, superposed coordinate sets and dendrograms can be accessed conveniently using a web interface. The database can be queried for protein pairs with sequence or structural similarities falling within a specified range. Thus PALI forms a useful resource to help in analysing the relationship between sequence and structure variation at a given level of sequence similarity. PALI also contains over 653 ‘orphans’ (single member families). Using the web interface involving PSI_BLAST and PHYLIP it is possible to associate the sequence of a new protein with one of the families in PALI and generate a phylogenetic tree combining the query sequence and proteins of known 3-D structure. The database with the web interfaced search and dendrogram generation tools can be accessed at http://pa uling.mbu.iisc.ernet.in/~pali.

PlantsP: a functional genomics database for plant phosphorylation

The PlantsP database is a curated database that combines information derived from sequences with experimental functional genomics information. PlantsP focuses on plant protein kinases and protein phosphatases. The database will specifically provide a resource for information on a collection of T-DNA insertion mutants (knockouts) in each protein kinase and phosphatase in Arabidopsis thaliana. PlantsP also provides a curated view of each protein that includes a comprehensive annotation of functionally related sequence motifs, sequence family definitions, alignments and phylogenetic trees, and descriptive information drawn directly from the literature. PlantsP is available at http://PlantsP.sdsc.edu.

Phylogenetics and the origin of species

A recent criticism that the biological species concept (BSC) unduly neglects phylogeny is examined under a novel modification of coalescent theory that considers multiple, sex-defined genealogical pathways through sexual organismal pedigrees. A competing phylogenetic species concept (PSC) also is evaluated from this vantage. Two analytical approaches are employed to capture the composite phylogenetic information contained within the braided assemblages of hereditary pathways of a pedigree: (i) consensus phylogenetic trees across allelic transmission routes and (ii) composite phenograms from quantitative values of organismal coancestry. Outcomes from both approaches demonstrate that the supposed sharp distinction between biological and phylogenetic species concepts is illusory. Historical descent and reproductive ties are related aspects of phylogeny and jointly illuminate biotic discontinuity.

Perspectives on archaeal diversity, thermophily and monophyly from environmental rRNA sequences.

Phylogenetic analysis of ribosomal RNA sequences obtained from uncultivated organisms of a hot spring in Yellowstone National Park reveals several novel groups of Archaea, many of which diverged from the crenarchaeal line of descent prior to previously characterized members of that kingdom. Universal phylogenetic trees constructed with the addition of these sequences indicate monophyly of Archaea, with modest bootstrap support. The data also show a specific relationship between low-temperature marine Archaea and some hot spring Archaea. Two of the environmental sequences are enigmatic: depending upon the data set and analytical method used, these sequences branch deeply within the Crenarchaeota, below the bifurcation between Crenarchaeota and Euryarchaeota, or even as the sister group to Eukaryotes. If additional data confirm either of the latter two placements, then the organisms represented by these ribosomal RNA sequences would merit recognition as a new kingdom, provisionally named "Korarchaeota."