Fast statistical tests for detecting heterotachy in protein evolution

The w statistic introduced by Lockhart et al. (1998. A covariotide model explains apparent phylogenetic structure of oxygenic photosynthetic lineages. Mol Biol Evol. 15:1183-1188) is a simple and easily calculated statistic intended to detect heterotachy by comparing amino acid substitution patterns between two monophyletic groups of protein sequences. It is defined as the difference between the fraction of varied sites in both groups and the fraction of varied sites in each group. The w test has been used to distinguish a covarion process from equal rates and rates variation across sites processes. Using simulation we show that the w test is effective for small data sets and for data sets that have low substitution rates in the groups but can have difficulties when these conditions are not met. Using site entropy as a measure of variability of a sequence site, we modify the w statistic to a w' statistic by assigning as varied in one group those sites that are actually varied in both groups but have a large entropy difference. We show that the w' test has more power to detect two kinds of heterotachy processes (covarion and bivariate rate shifts) in large and variable data. We also show that a test of Pearson's correlation of the site entropies between two monophyletic groups can be used to detect heterotachy and has more power than the w' test. Furthermore, we demonstrate that there are settings where the correlation test as well as w and w' tests do not detect heterotachy signals in data simulated under a branch length mixture model. In such cases, it is sometimes possible to detect heterotachy through subselection of appropriate taxa. Finally, we discuss the abilities of the three statistical tests to detect a fourth mode of heterotachy: lineage-specific changes in proportion of variable sites.     

Construction of a Phylogenetic Tree of Photosynthetic Prokaryotes Based on Average Similarities of Whole Genome Sequences

Phylogenetic trees have been constructed for a wide range of organisms using gene sequence information, especially through the identification of orthologous genes that have been vertically inherited. The number of available complete genome sequences is rapidly increasing, and many tools for construction of genome trees based on whole genome sequences have been proposed. However, development of a reasonable method of using complete genome sequences for construction of phylogenetic trees has not been established. We have developed a method for construction of phylogenetic trees based on the average sequence similarities of whole genome sequences. We used this method to examine the phylogeny of 115 photosynthetic prokaryotes, i.e., cyanobacteria, Chlorobi, proteobacteria, Chloroflexi, Firmicutes and nonphotosynthetic organisms including Archaea. Although the bootstrap values for the branching order of phyla were low, probably due to lateral gene transfer and saturated mutation, the obtained tree was largely consistent with the previously reported phylogenetic trees, indicating that this method is a robust alternative to traditional phylogenetic methods.     

Heterotachy processes in rhodophyte-derived secondhand plastid genes: Implications for addressing the origin and evolution of dinoflagellate plastids

Serial transfer of plastids from one eukaryotic host to another is the key process involved in evolution of secondhand plastids. Such transfers drastically change the environment of the plastids and hence the selection regimes, presumably leading to changes over time in the characteristics of plastid gene evolution and to misleading phylogenetic inferences. About half of the dinoflagellate protists species are photosynthetic and unique in harboring a diversity of plastids acquired from a wide range of eukaryotic algae. They are therefore ideal for studying evolutionary processes of plastids gained through secondary and tertiary endosymbioses. In the light of these processes, we have evaluated the origin of 2 types of dinoflagellate plastids, containing the peridinin or 19'-hexanoyloxyfucoxanthin (19'-HNOF) pigments, by inferring the phylogeny using "covarion" evolutionary models allowing the pattern of among-site rate variation to change over time. Our investigations of genes from secondary and tertiary plastids derived from the rhodophyte plastid lineage clearly reveal "heterotachy" processes characterized as stationary covarion substitution patterns and changes in proportion of variable sites across sequences. Failure to accommodate covarion-like substitution patterns can have strong effects on the plastid tree topology. Importantly, multigene analyses performed with probabilistic methods using among-site rate and covarion models of evolution conflict with proposed single origin of the peridinin- and 19'-HNOF-containing plastids, suggesting that analysis of secondhand plastids can be hampered by convergence in the evolutionary signature of the plastid DNA sequences. Another type of sequence convergence was detected at protein level involving the psaA gene. Excluding the psaA sequence from a concatenated protein alignment grouped the peridinin plastid with haptophytes, congruent with all DNA trees. Altogether, taking account of complex processes involved in the evolution of dinoflagellate plastid sequences (both at the DNA and amino acid level), we demonstrate the difficulty of excluding independent, tertiary origin for both the peridinin and 19'-HNOF plastids involving engulfment of haptophyte-like algae. In addition, the refined topologies suggest the red algal order, Porphyridales, as the endosymbiont ancestor of the secondary plastids in cryptophytes, haptophytes, and heterokonts.     

Identification of microorganisms by PCR amplification and sequencing of a universal amplified ribosomal region present in both prokaryotes and eukaryotes

The small ribosomal subunit contains 16S rRNA in prokaryotes and 18S rRNA in eukaryotes. Even though it has been known that some small ribosomal sequences are conserved in 16S rRNA and 18S rRNA molecules, they have been used separately for taxonomic and phylogenetic studies. Here, we report the existence of two highly conserved ribosomal sequences in all organisms that allow the amplification of a zone containing approximately 495 bp in prokaryotes and 508 bp in eukaryotes which we have named the "Universal Amplified Ribosomal Region" (UARR). Amplification and sequencing of this zone is possible using the same two universal primers (U1F and U1R) designed on the basis of two highly conserved ribosomal sequences. The UARR encompasses the V6, V7 and V8 domains from SSU rRNA in both prokaryotes and eukaryotes. The internal sequence of this zone in prokaryotes and eukaryotes is variable and the differences become less marked on descent from phyla to species. Nevertheless, UARR sequence allows species from the same genus to be differentiated. Thus, by UARR sequence analysis the construction of universal phylogenetic trees is possible, including species of prokaryotic and eukaryotic microorganisms together. Single isolates of prokaryotic and eukaryotic microorganisms from different sources can be identified by amplification and sequencing of UARR using the same pair of primers and the same PCR and sequencing conditions.     

Testing a covariotide model of DNA substitution

The concomitantly variable codons hypothesis of DNA substitution argues that at any time only a fraction of the codons in a gene are capable of accepting a mutation. However, as mutations are fixed at some positions in a gene, the sites that are potentially variable also change because of changed functional constraints. This hypothesis has been termed the "covarion" hypothesis or when the model is applied to nucleotides, the "covariotide" hypothesis. The covarion-covariotide model has proven to be remarkably difficult to test. Here I examine a covariotide hypothesis for 11 genes using a likelihood ratio test. I show that in nine of the genes examined a covariotide model provides a better explanation of the data than a model that does not allow constraints to change over time.     

Species boundaries and phylogenetic relationships within the green algal genus Codium (Bryopsidales) based on plastid DNA sequences

Despite the potential model role of the green algal genus Codium for studies of marine speciation and evolution, there have been difficulties with species delimitation and a molecular phylogenetic framework was lacking. In the present study, 74 evolutionarily significant units (ESUs) are delimited using 227 rbcL exon 1 sequences obtained from specimens collected throughout the genus' range. Several morpho-species were shown to be poorly defined, with some clearly in need of lumping and others containing pseudo-cryptic diversity. A phylogenetic hypothesis of 72 Codium ESUs is inferred from rbcL exon 1 and rps3-rpl16 sequence data using a conventional nucleotide substitution model (GTR+Gamma+I), a codon position model and a covariotide (covarion) model, and the fit of a multitude of substitution models and alignment partitioning strategies to the sequence data is reported. Molecular clock tree rooting was carried out because outgroup rooting was probably affected by phylogenetic bias. Several aspects of the evolution of morphological features of Codium are discussed and the inferred phylogenetic hypothesis is used as a framework to study the biogeography of the genus, both at a global scale and within the Indian Ocean.     

Evolution of glutamine synthetase in heterokonts: evidence for endosymbiotic gene transfer and the early evolution of photosynthesis

Although the endosymbiotic evolution of chloroplasts through primary and secondary associations is well established, the evolutionary timing and stability of the secondary endosymbiotic events is less well resolved. Heterokonts include both photosynthetic and nonphotosynthetic members and the nonphotosynthetic lineages branch basally in phylogenetic reconstructions. Molecular and morphological data indicate that heterokont chloroplasts evolved via a secondary endosymbiosis, involving a heterotrophic host cell and a photosynthetic ancestor of the red algae and this endosymbiotic event may have preceded the divergence of heterokonts and alveolates. If photosynthesis evolved early in this lineage, nuclear genomes of the nonphotosynthetic groups may contain genes that are not essential to photosynthesis but were derived from the endosymbiont genome through gene transfer. These genes offer the potential to trace the evolutionary history of chloroplast gains and losses within these lineages. Glutamine synthetase (GS) is essential for ammonium assimilation and glutamine biosynthesis in all organisms. Three paralogous gene families (GSI, GSII, and GSIII) have been identified and are broadly distributed among prokaryotic and eukaryotic lineages. In diatoms (Heterokonta), the nuclear-encoded chloroplast and cytosolic-localized GS isoforms are encoded by members of the GSII and GSIII family, respectively. Here, we explore the evolutionary history of GSII in both photosynthetic and nonphotosynthetic heterokonts, red algae, and other eukaryotes. GSII cDNA sequences were obtained from two species of oomycetes by polymerase chain reaction amplification. Additional GSII sequences from eukaryotes and bacteria were obtained from publicly available databases and genome projects. Bayesian inference and maximum likelihood phylogenetic analyses of GSII provided strong support for the monophyly of heterokonts, rhodophytes, chlorophytes, and plants and strong to moderate support for the Opisthokonts. Although the phylogeny is reflective of the unikont/bikont division of eukaryotes, we propose based on the robustness of the phylogenetic analyses that the heterokont GSII gene evolved via endosymbiotic gene transfer from the nucleus of the red-algal endosymbiont to the nucleus of the host. The lack of GSIII sequences in the oomycetes examined here further suggests that the GSIII gene that functions in the cytosol of photosynthetic heterokonts was replaced by the endosymbiont-derived GSII gene.     

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and,with the exception of endosymbiotic gene transfers,few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence–absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.     

Phylogenetic positions of several amitochondriate protozoa-Evidence from phylogenetic analysis of DNA topoisomerase II

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their 'lacking mitochondria' and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of 'long-branch attraction' appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Topological Estimation Biases with Covarion Evolution

Covarion processes allow changes in evolutionary rates at sites along the branches of a phylogenetic tree. Covarion-like evolution is increasingly recognized as an important mode of protein evolution. Several recent reports suggest that maximum likelihood estimation employing covarion models may support different optimal topologies than estimation using standard rates-across-sites (RAS) models. However, it remains to be demonstrated that ignoring covarion evolution will generally result in topological misestimation. In this study we performed analytical and theoretical studies of limiting distances under the covarion model and four-taxon tree simulations to investigate the extent to which the covarion process impacts on phylogenetic estimation. In particular, we assessed the limits of an RAS model-based maximum likelihood method to recover the phylogenies when the sequence data were simulated under the covarion processes. We find that, when ignored, covarion processes can induce systematic errors in phylogeny reconstruction. Surprisingly, when sequences are evolved under a covarion process but an RAS model is used for estimation, we find that a long branch repel bias occurs.     

Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplo-monads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be poly-phyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

A mixed branch length model of heterotachy improves phylogenetic accuracy

Evolutionary relationships are typically inferred from molecular sequence data using a statistical model of the evolutionary process. When the model accurately reflects the underlying process, probabilistic phylogenetic methods recover the correct relationships with high accuracy. There is ample evidence, however, that models commonly used today do not adequately reflect real-world evolutionary dynamics. Virtually all contemporary models assume that relatively fast-evolving sites are fast across the entire tree, whereas slower sites always evolve at relatively slower rates. Many molecular sequences, however, exhibit site-specific changes in evolutionary rates, called "heterotachy." Here we examine the accuracy of 2 phylogenetic methods for incorporating heterotachy, the mixed branch length model--which incorporates site-specific rate changes by summing likelihoods over multiple sets of branch lengths on the same tree--and the covarion model, which uses a hidden Markov process to allow sites to switch between variable and invariable as they evolve. Under a variety of simple heterogeneous simulation conditions, the mixed model was dramatically more accurate than homotachous models, which were subject to topological biases as well as biases in branch length estimates. When data were simulated with strong versions of the types of heterotachy observed in real molecular sequences, the mixed branch length model was more accurate than homotachous techniques. Analyses of empirical data sets confirmed that the mixed branch length model can improve phylogenetic accuracy under conditions that cause homotachous models to fail. In contrast, the covarion model did not improve phylogenetic accuracy compared with homotachous models and was sometimes substantially less accurate. We conclude that a mixed branch length approach, although not the solution to all phylogenetic errors, is a valuable strategy for improving the accuracy of inferred trees.     

Testing the covarion hypothesis of molecular evolution

The covarion hypothesis of molecular evolution states that the fixation of mutations may alter the probability that any given position will fix the next change. Tests of this hypothesis using the divergence of real sequences are compromised because models of rate variation among sites (e.g., the gamma version of the one-parameter equation) predict sequence divergence values similar to those for the covarion process. This study therefore focuses on the extent to which the varied and unvaried codons of two well-diverged taxa are the same, because fewer are expected by the covarion hypothesis than by the gamma model. The data for these tests are the protein sequences of Cu, Zn superoxide dismutase (SOD) for mammals and plants. Simulation analyses show that the covarion hypothesis makes better predictions about the frequencies of varied and unhit positions in common between these two taxa than does the gamma version of the one-parameter model. Furthermore, the analysis of SOD tertiary structure demonstrates that mammal and plant variabilities are distributed differently on the protein. These results support the conclusions that the variable and invariable codons of mammal and plant SODs are different and that the covarion model explains the evolution of this protein better than the gamma version of the one-parameter process. Unlike other models, the covarion hypothesis accounts for rate fluctuations among positions over time, which is an important parameter of molecular evolution.      

Testing for covarion-like evolution in protein sequences

The covarion hypothesis of molecular evolution proposes that selective pressures on an amino acid or nucleotide site change through time, thus causing changes of evolutionary rate along the edges of a phylogenetic tree. Several kinds of Markov models for the covarion process have been proposed. One model, proposed by Huelsenbeck (2002), has 2 substitution rate classes: the substitution process at a site can switch between a single variable rate, drawn from a discrete gamma distribution, and a zero invariable rate. A second model, suggested by Galtier (2001), assumes rate switches among an arbitrary number of rate classes but switching to and from the invariable rate class is not allowed. The latter model allows for some sites that do not participate in the rate-switching process. Here we propose a general covarion model that combines features of both models, allowing evolutionary rates not only to switch between variable and invariable classes but also to switch among different rates when they are in a variable state. We have implemented all 3 covarion models in a maximum likelihood framework for amino acid sequences and tested them on 23 protein data sets. We found significant likelihood increases for all data sets for the 3 models, compared with a model that does not allow site-specific rate switches along the tree. Furthermore, we found that the general model fit the data better than the simpler covarion models in the majority of the cases, highlighting the complexity in modeling the covarion process. The general covarion model can be used for comparing tree topologies, molecular dating studies, and the investigation of protein adaptation.     

Comparison of dissimilarity patterns of E coli, yeast and mammalian tRNAs.

A number of experimental approaches have been developed for identification of recognition (identity) sites in tRNAs. Along with them a theoretical methodology has been proposed by McClain et al that is based on concomitant analysis of all tRNA sequences from a given species. This approach allows an evaluation of nucleotide combinations present in isoacceptor tRNAs specific for the given amino acid, and not present in equivalent positions in cloverleaf structure in other tRNAs of the same organism. These elements predicted from computer analysis of the databank could be tested experimentally for their participation in forming recognition sites. The correlation between theoretical predictions and experimental data appeared promising. The aim of the present work consisted of introducing further improvements into McClain's procedure by: i), introducing into analysis a variable region in tRNAs which had not been previously considered; to accomplish this, 'normalization' of variable nucleotides was suggested, based on primary and tertiary structures of tRNAs; ii), developing a new procedure for comparison of patterns for synonymous and non-synonymous tRNAs from different organisms; iii), analysis of 3- and 4-positional contacts between tRNAs and enzymes in addition to a formerly used 2-positional model. A systematic application of McClain's procedure to mammalian, yeast and E coli tRNAs led to the following results: i), imitancy patterns for non-synonymous tRNAs of any amino acid specificity and from any organisms analysed so far overlap by no more than 30%, providing a structural basis for discrimination with high fidelity between cognate and non-cognate tRNAs; ii), the predicted identity sites are non-randomly distributed within tRNA molecules; the dominant role is ascribed to only two regions--anticodon and amino acid stem which are located far apart from one another at extremes of all tRNA molecules; iii), the imitancy patterns for synonymous tRNAs in lower (yeast) and higher (mammalian) eukaryotes are similar but not identical; iv), distribution of predicted identity sites in the cloverleaf structure in prokaryotes and eukaryotes is essentially different: in eubacterial tRNAs the major role in recognition plays anticodon and/or amino acid acceptor stem, whereas in eukaryotic (both unicellular and multicellular) tRNAs the remaining part of the molecules is also involved in recognition; v), the imitancy patterns of synonymous tRNAs from prokaryotes and eukaryotes are dissimilar, this observation leads to the prediction that the tRNA identity sites for the same amino acid in prokaryotes and eukaryotes may differ.     

Species-specific patterns of DNA bending and sequence.

Nucleotide sequences in the GenEMBL database were analyzed using strategies designed to reveal species-specific patterns of DNA bending and DNA sequence. The results uncovered striking species-dependent patterns of bending with more variations among individual organisms than between prokaryotes and eukaryotes. The frequency of bent sites in sequences from different bacteria was related to genomic A + T content and this relationship was confirmed by electrophoretic analysis of genomic DNA. However, base composition was not an accurate predictor for DNA bending in eukaryotes. Sequences from C. elegans exhibited the highest frequency of bent sites in the database and the RNA polymerase II locus from the nematode was the most bent gene in GenEMBL. Bent DNA extended throughout most introns and gene flanking segments from C.elegans while exon regions lacked A-tract bending characteristics. Independent evidence for the strong bending character of this genome was provided by electrophoretic studies which revealed that a large number of the fragments from C.elegans DNA exhibited anomalous gel mobilities when compared to genomic fragments from over 20 other organisms. The prevalence of bent sites in this genome enabled us to detect selectively C.elegans sequences in a computer search of the database using as probes C.elegans introns, bending elements, and a 20 nucleotide consensus sequence for bent DNA. This approach was also used to provide additional examples of species-specific sequence patterns in eukaryotes where it was shown that (A) greater than or equal to 10 and (A.T) greater than or equal to 5 tracts are prevalent throughout the untranslated DNA of D.discodium and P.falciparum, respectively. These results provide new insight into the organization of eukaryotic DNA because they show that species-specific patterns of simple sequences are found in introns and in other untranslated regions of the genome.     

Estimation of phylogenetic inconsistencies in the three domains of life

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.     

Nonuniformity of nucleotide substitution rates in molecular evolution: Computer simulation and analysis of 5S ribosomal RNA sequences

Summary The effects of temporal (among different branches of a phylogeny) and spatial (among different nucleotide sites within a gene) nonuniformities of nucleotide substitution rates on the construction of phylogenetic trees from nucleotide sequences are addressed. Spatial nonuniformity may be estimated by using Shannon's (1948) entropy formula to measure the Relative Nucleotide Variability (RNV) at each nucleotide site in an aligned set of sequences; this is demonstrated by a comparative analysis of 5S rRNAs. New methods of constructing phylogenetic trees are proposed that augment the Unweighted Pair-Group Using Arithmetic Averages (UPGMA) algorithm by estimating and compensating for both spatial and temporal nonuniformity in substitution rates. These methods are evaluated by computer simulations of 5S rRNA evolution that include both kinds of nonuniformities. It was found that the proposed Reference Ratio Method improved both the ability to reconstruct the correct topology of a tree and also the estimation of branch lengths as compared to UPGMA. A previous method (Farris et al. 1970; Klotz et al. 1979; Li 1981) was found to be less successful in reconstructing topologies when there is high probability of multiple mutations at some sites. Phylogenetic analyses of 5S rRNA sequences support the endosymbiotic origins of both chloroplasts and mitochondria, even though the latter exhibit an accelerated rate of nucleotide substitution. Phylogenetic trees also reveal an adaptive radiation within the eubacteria and another within the eukaryotes for the origins of most major phyla within each group during the Precambrian era.     

On the origin of prokaryotes

It has been a tacit assumption of evolutionary theory that the closest surviving relatives of the first cellular organisms are to be found among prokaryotes. This paper draws attention to the fact that many stages of evolution appear to have been accompanied by physical loss of superfluous DNA. It is postulated that the genomes of prokaryotes—where almost every gene is represented by one copy only—represent the results of this process carried to its extreme. On this basis certain features of very early evolution which have been eliminated from prokaryotes may survive in eukaryotes. If correct, the hypothesis would require a careful re-evaluation of the assumptions underlying use of some sequence data to construct phylogenetic trees.     

Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, En-tamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group—Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.