Search

Background: Cystic or alveolar echinococcosis caused by the larval stages of Echinococcus spp. is a very severe zoonotic helminth infection. Echinococcus shiquicus is a newly discovered species that has only been reported in the Qinghai and Sichuan provinces of the Qinghai-Tibet plateau, China where, to date, it has only been confirmed in Tibetan foxes and wild small mammal populations of the Tibetan plateau. Information on its genetic and evolutionary diversity is scanty. The aim of this study was to investigate the prevalence of E. shiquicus in plateau pikas (Ochotona curzoniae), a known intermediate host, and to determine the genetic variation and phylogenetic relationship of the E. shiquicus population in the Tibet region of China based on mitochondrial DNA., Methods: Echinococcus shiquicus samples were collected from Damxung and Nyêmo counties (located in Tibet Autonomous Region, China). The mitochondrial cox1 and nad1 gene sequences were analyzed, and the genetic diversity and epidemiology of E. shiquicus in the region were discussed based on the results., Results: The prevalence of E. shiquicus in pikas in Damxung and Nyêmo counties was 3.95% (6/152) and 6.98% (9/129), respectively. In combination with previous public sequence data, the haplotype analysis revealed 12 haplotypes (H) characterized by two distinct clusters (I and II), and a sequence distance of 99.1-99.9% from the reference haplotype (H1). The diversity and neutrality indices for the entire E. shiquicus populations were: haplotype diversity (Hd) ± standard deviation (SD) 0.862 ± 0.035; nucleotide diversity (Hd ± SD) 0.0056 ± 0.0003; Tajima's D 0.876 (P > 0.05); and Fu's F 6.000 (P > 0.05)., Conclusions: This was the first analysis of the newly discovered E. shiquicus in plateau pikas in the Tibet Autonomous Region of China. The neutrality indices suggest a deficiency of alleles, indicative of a recent population bottleneck.

A novel nocardioform strain, CICC 11023 T , was isolated from a tissue biopsy of neck lesions of a patient with primary cutaneous nocardiosis and characterized to establish its taxonomic position. The morphological, biochemical, physiological and chemotaxonomic properties of strain CICC 11023 T were consistent with classification in the genus Nocardia . Whole-cell hydrolysates were rich in meso-diaminopimelic acid, galactose, arabinose and fructose. Mycolic acids were present. The major polar lipids were diphosphatidylglycerol, phosphatidylethanolamine, one unidentified phospholipid and two unidentified lipids, and the predominant menaquinone was cyclo MK-8 (H4, ω-cyclo). The main fatty acids (>5 %) were C 18 : 0 10-methyl (TBSA), C 16 : 0 , summed feature 4 (C 16 : 1 trans 9/C 15 : 0 iso 2OH), C 15 : 0 and C 17 : 0 10-methyl. Phylogenetic analyses based on 16S rRNA gene sequences revealed that the isolate is most closely related (>98 % similarity) to the type strains Nocardia ninae OFN 02.72 T , Nocardia iowensis UI 122540 T and Nocardia alba YIM 30243 T , and phylogenetic analysis of gyrB gene sequences showed similarity (89.1-92.2 %) to Nocardia vulneris NBRC 108936 T , Nocardia brasiliensis IFM 0236 T and Nocardia exalbida IFM 0803 T . DNA-DNA hybridization results for strain CICC 11023 T compared to Nocardia type strains ranged from 20.4 to 35.4 %. The genome of strain CICC 11023 T was 8.78 Mbp with a G+C content of 67.4 mol% overall. The average nucleotide identity (ANI) values between strain CICC 11023 T and N. alba YIM 30243 T were low (OrthoANIu=77.47 %), and the ANI values between strain CICC 11023 T and N. vulneris NBRC 108936 T were low (OrthoANIu=83.75 %). Consequently, strain CICC 11023 T represents a novel Nocardia species on the basis of this polyphasic study, for which the name Nocardia colli sp. nov. is proposed. The type strain is CICC 11023 T (=KCTC 39837 T ).

This study determined the first complete mitochondrial genome (mitogenome) of a stonefly, Scopura longa (Plecoptera: Scopuridae), and reconstructed a phylogeny based on two datasets of mitogenomes in eighteen available stoneflies to examine the relationships among Plecoptera. The complete mitogenome of S. longa is a circular molecule of 15,798bp in size. It contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). Most PCGs used standard ATN start codons and TAN stop codons. All tRNAs could be folded as typical cloverleaf secondary structures except tRNA Ser(AGN) , with the dihydrouridine (DHU) arm missing. Tandem repeat regions, two potential stem-loop (SL) structures, Poly T structures and G+C-rich regions are detected in the control region. Finally, the phylogenetic relationships among the families within the Arctoperlaria were reconstructed. The topological structures of the two trees were almost identical. The present phylogenetic analysis shows that S. longa belongs to the infraorder Euholognatha. The monophyly of each family is generally well supported based on nucleotide sequences. The Pteronarcyidae is sister to ((Peltoperlidae+Styloperlidae)+(Perlidae+(Perlodidae+Chloroperlidae)))., (Copyright © 2018 Elsevier B.V. All rights reserved.)

The core of the Vibrio Harveyi clade contains V. harveyi, V. campbellii, V. owensii, V. jasicida, and V. rotiferianus. They are well recognized aquatic animal pathogens, but misclassification has been common due to similarities in their rDNA sequences and phenotypes. To better understand their evolutionary relationships and functional features, we sequenced a shrimp pathogen strain V. harveyi 1114GL, reclassified it as V. campbellii and compared this and 47 other sequenced Vibrio genomes in the Harveryi clade. A phylogeny based on 1,775 genes revealed that both V. owensii and V. jasicida were closer to V. campbellii than to V. harveyi and that V. campbellii strains can be divided into two distinct groups. Species-specific genes such as intimin and iron acquisition genes were identified in V. campbellii. In particular, the 1114GL strain contains two bacterial immunoglobulin-like genes for cell adhesion with 22 Big_2 domains that have been extensively reshuffled and are by far the most expanded among all species surveyed in this study. The 1114GL strain differed from ATCC BAA-1116 by ~9% at the synonymous sites, indicating high diversity within V. campbellii. Our study revealed the characteristics of V. campbellii in the Harveyi clade and the genetic basis for their wide-spread pathogenicity., Competing Interests: The authors declare no competing financial interests.

The enzymes catalyzing the first two reactions in the sulfate activation pathway, ATP-sulfurylase (S) and APS-kinase (K), are fused as 'KS' in animals but are fused as 'SK' in select bacteria and fungi. We have discovered a novel triple fusion protein of K, S, and pyrophosphatase (P) in several protozoan genomes within the Stramenopile lineage. These triple domain fusion proteins led us to hypothesize that pyrophosphatase enzymes and sulfate activation enzymes physically interact to impact the thermodynamics of the sulfate activation pathway. In support of this hypothesis, we demonstrate through biochemical assays that separately encoded KS and P proteins physically interact and that KS/P complexes activate more sulfate than KS alone. We also conclude on the basis of phylogenetic analyses that all known KS fusions originate from a single fusion event early in the eukaryotic lineage. Strikingly, our analyses support the same conclusion for all known SK fusions. These observations indicate that the patchwork of fused and nonfused S and K genes observed in modern-day eukaryotes and prokaryotes are the result of the two ancestral fusion genes evolving by an assortment of gene fissions, duplications, deletions, and horizontal transfers in different lineages. Our integrative use of genomics, phylogenetics, and biochemistry to characterize pyrophosphatase as a new member of the sulfate activation pathway should be effective at identifying new protein members and connections in other molecular pathways.

Many Saccharomyces cerevisiae duplicate genes that were derived from an ancient whole-genome duplication (WGD) unexpectedly show a small synonymous divergence (K(S)), a higher sequence similarity to each other than to orthologues in Saccharomyces bayanus, or slow evolution compared with the orthologue in Kluyveromyces waltii, a non-WGD species. This decelerated evolution was attributed to gene conversion between duplicates. Using approximately 300 WGD gene pairs in four species and their orthologues in non-WGD species, we show that codon-usage bias and protein-sequence conservation are two important causes for decelerated evolution of duplicate genes, whereas gene conversion is effective only in the presence of strong codon-usage bias or protein-sequence conservation. Furthermore, we find that change in mutation pattern or in tDNA copy number changed codon-usage bias and increased the K(S) distance between K. waltii and S. cerevisiae. Intriguingly, some proteins showed fast evolution before the radiation of WGD species but little or no sequence divergence between orthologues and paralogues thereafter, indicating that functional conservation after the radiation may also be responsible for decelerated evolution in duplicates.

It is firmly believed that ancestral primates were nocturnal, with nocturnality having been maintained in most prosimian lineages. Under this traditional view, the opsin genes in all nocturnal prosimians should have undergone similar degrees of functional relaxation and accumulated similar extents of deleterious mutations. This expectation is rejected by the short-wavelength (S) opsin gene sequences from 14 representative prosimians. We found severe defects of the S opsin gene only in lorisiforms, but no defect in five nocturnal and two diurnal lemur species and only minor defects in two tarsiers and two nocturnal lemurs. Further, the nonsynonymous-to-synonymous rate ratio of the S opsin gene is highest in the lorisiforms and varies among the other prosimian branches, indicating different time periods of functional relaxation among lineages. These observations suggest that the ancestral primates were diurnal or cathemeral and that nocturnality has evolved several times in the prosimians, first in the lorisiforms but much later in other lineages. This view is further supported by the distribution pattern of the middle-wavelength (M) and long-wavelength (L) opsin genes among prosimians.

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage., (2002 Elsevier Science (USA).)

In the reconstruction of a large phylogenetic tree, the most difficult part is usually the problem of how to explore the topology space to find the optimal topology. We have developed a "divide-and-conquer" heuristic algorithm in which an initial neighbor-joining (NJ) tree is divided into subtrees at internal branches having bootstrap values higher than a threshold. The topology search is then conducted by using the maximum-likelihood method to reevaluate all branches with a bootstrap value lower than the threshold while keeping the other branches intact. Extensive simulation showed that our simple method, the neighbor-joining maximum-likelihood (NJML) method, is highly efficient in improving NJ trees. Furthermore, the performance of the NJML method is nearly equal to or better than existing time-consuming heuristic maximum-likelihood methods. Our method is suitable for reconstructing relatively large molecular phylogenetic trees (number of taxa >/= 16).

The minimum sum of branch lengths (S), or the minimum evolution (ME) principle, has been shown to be a good optimization criterion in phylogenetic inference. Unfortunately, the number of topologies to be analyzed is computationally prohibitive when a large number of taxa are involved. Therefore, simplified, heuristic methods, such as the neighbor-joining (NJ) method, are usually employed instead. The NJ method analyzes only a small number of trees (compared with the size of the entire search space); so, the tree obtained may not be the ME tree (for which the S value is minimum over the entire search space). Different compromises between very restrictive and exhaustive search spaces have been proposed recently. In particular, the "stepwise algorithm" (SA) utilizes what is known in computer science as the "beam search," whereas the NJ method employs a "greedy search." SA is virtually guaranteed to find the ME trees while being much faster than exhaustive search algorithms. In this study we propose an even faster method for finding the ME tree. The new algorithm adjusts its search exhaustiveness (from greedy to complete) according to the statistical reliability of the tree node being reconstructed. It is also virtually guaranteed to find the ME tree. The performances and computational efficiencies of ME, SA, NJ, and our new method were compared in extensive simulation studies. The new algorithm was found to perform practically as well as the SA (and, therefore, ME) methods and slightly better than the NJ method. For searching for the globally optimal ME tree, the new algorithm is significantly faster than existing ones, thus making it relatively practical for obtaining all trees with an S value equal to or smaller than that of the NJ tree, even when a large number of taxa is involved., (Copyright 2000 Academic Press.)

The phylogenetic relationships among worldwide species of genus Ochotona were investigated by sequencing mitochondrial cytochrome b and ND4 genes. Parsimony and neighbor-joining analyses of the sequence data yielded congruent results that strongly indicated three major clusters: the shrub-steppe group, the northern group, and the mountain group. The subgeneric classification of Ochotona species needs to be revised because each of the two subgenera in the present classification contains species from the mountain group. To solve this taxonomic problem so that each taxon is monophyletic, i.e. , represents a natural clade, Ochotona could be divided into three subgenera, one for the shrub-steppe species, a second for the northern species, and a third for the mountain species. The inferred tree suggests that the differentiation of this genus in the Palearctic Region was closely related to the gradual uplifting of the Tibet (Qinghai-Xizang) Plateau, as hypothesized previously, and that vicariance might have played a major role in the differentiation of this genus on the Plateau. On the other hand, the North American species, O. princeps, is most likely a dispersal event, which might have happened during the Pliocene through the opening of the Bering Strait. The phylogenetic relationships within the shrub-steppe group are worth noting in that instead of a monophyletic shrub-dwelling group, shrub dwellers and steppe dwellers are intermingled with each other. Moreover, the sequence divergence within the sister taxa of one steppe dweller and one shrub dweller is very low. These findings support the hypothesis that pikas have entered the steppe environment several times and that morphological similarities within steppe dwellers were due to convergent evolution., (Copyright 2000 Academic Press.)

The phylogenetic relationships of nine genera in four tribes of the family Brassicaceae were estimated from the sequences of the internal transcribed spacer region (ITS) of the 18S-25S nuclear ribosomal DNA. The entire ITS region of 16 accessions belonging to 10 species of seven genera was sequenced. Eight published sequences of Brassicaceae were also used. A total of 27 sequences were included in this study; four of them were found to be pseudogenes. Both the neighbor-joining and the parsimony trees suggest that the nine genera can be divided into three groups: (1) Arabidopsis, Cardaminopsis, Capsella, and Lepidium; (2) Rorippa and Cardamine; and (3) Brassica, Sinapis, and Raphanus. In contradiction to the proposal that Cardaminopsis and Arabidopsis be put into an expanded tribe Arabideae, our data show that these two genera are more closely related to Capsella and Lepidium (tribe Lepidieae) than to Rorippa and Cardamine (tribe Arabideae). Further, our data show that within the tribe Brassiceae, Raphanus is more closely related to B. nigra than to the B. oleracea/B. rapa clade. This result is in agreement with the nuclear data obtained in several studies, but is in conflict with the RFLP data of mitochondrial and chloroplast DNA. As pointed out by previous authors, it is possible that Raphanus is a hybrid between the B. nigra and B. oleracea/B. rapa lineages with the latter as the maternal parent., (Copyright 1999 Academic Press.)

Coalescent theory represents the most significant progress in theoretical population genetics in the past two decades of this century. It is now widely recognized as a cornerstone for rigorous statistical analyses of molecular data from populations. In the future, challenges from the rapidly expanding body of molecular data will continue to inject fresh blood into the development of coalescent theory. As the 21st century looms, we review recent progress in coalescent theory and discuss some prospects for future developments, particularly with regard to human population genetics., (Copyright 1999 Academic Press.)

The relationship between change in genomic GC content and protein evolution in bacteria was studied by simple correlational analysis (at the genus level) and by Felsenstein's (1985) independent contrast test. We first used the dnaA gene in bacteria as an example to show (1) that the amino acid composition of a protein can be dramatically affected by mutational pressure (the genomic GC content), (2) that surprisingly, deleting relatively closely-related genera may increase rather than decrease the correlation between genomic GC content and amino acid composition, and (3) that most unexpectedly, as the genomic GC content increases, both strongly hydrophobic and strongly hydrophilic amino acids tend to change to ambivalent amino acids, suggesting that the majority of these amino acid substitutions are not caused by positive Darwinian selection. These patterns were then also shown to hold for the 14 other genes studied, indicating their generality for the evolution of bacterial proteins. As directional mutation pressure can affect the amino acid composition of proteins, it may mislead phylogenetic inference, even if protein instead of DNA sequences are used.

The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, we developed formulas for correcting the estimation biases of the paralinear and LogDet distances, i.e., the bias-corrected distance is estimated by dc = d - 2var(d), where d and var(d) are the estimated distance and sampling variance, respectively. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, we developed a method for estimating the variance-covariance matrix of paralinear distances, so that statistical tests of DNA phylogenies can be conducted in the nonstationary case. Third, a new LogDet-based method for testing the molecular clock hypothesis was developed under nonstationary nucleotide frequencies.

To determine the evolutionary positions of the conifer genera Amentotaxus, Phyllocladus, and Nageia, we obtained 18S rRNA sequences from 11 new taxa representing the major living orders and families of gymnosperms. With the published Chlamydomonas as an outgroup, phylogenetic analyses of our new data and available sequences indicate that (1) the Gnetales form a monophyletic group, which is an outgroup to the conifers, (2) the conifers are monophyletic, (3) Taxaceae, Cephalotaxaceae, Cupressaceae, and Taxodiaceae form a monophyletic group, (4) Amentotaxus is closer to Torreya than to Cephalotaxus, suggesting that Amentotaxus is better to be classified as a member of Taxaceae, (5) Phyllocladus, Dacrycarpus, Podocarpus, and Nageia form a monophyletic group, and (6) Pinaceae is an out-group to the other families of conifers. Our finding that Phyllocladus is a sister group of the Podocarpaceae disagrees with the suggestion that the phylloclade of the genus is an ancient structure and that the genus is a terminal taxon within the Podocarpaceae. The genus Nageia is more closely related to Podocarpus than to Dacrycarpus and was derived from within the Podocarpaceae. In conclusion, our data indicate that in conifers, the uniovulate cone occurred independently in Taxacaeae and Cephalotaxaceae, and in Podocarpaceae after the three families separated from Pinaceae, and support the hypothesis that the uniovulate cone is derived from reduction of a multiovulate cone.

The bootstrap is a statistical technique that is widely used to assess confidence limits on phylogenies. We show that the power of the bootstrap test is lower than those of the C and S tests suggested by Felsenstein, unless the critical value employed in the bootstrap test is correctly selected. If the 95% critical value is used, the bootstrap proportions are underestimates of the confidence level when the number of possible alternative topologies is three or more; the degree of underestimation increases with the number of competing alternative topologies. To overcome this problem, we propose the complete-and-partial bootstrap technique as a method for obtaining an unbiased estimate of the confidence level. The method is based on a multinomial model of many alternatives among which the choice is to be made. The complete-and-partial bootstrap technique can be used to estimate the effective number of competing alternative topologies and the confidence level of the monophyly of a particular group of taxa or of an inferred tree topology. This approach can be used with the maximum parsimony or neighbor-joining tree reconstruction method.

The evolutionary position of the yew family, Taxaceae, has been very controversial. Some plant taxonomists strongly advocate excluding Taxaceae from the conifer order and raising its taxonomic status to a new order or even class because of its absence of seed cones, contrary to the case in the majority of conifers. However, other authors believe that the Taxaceae are not fundamentally different from the rest of the conifers except in that they possess the most reduced solitary-ovule cones. To resolve the controversy, we have sequenced the 18S rRNA genes from representative gymnosperms: Taxus mairei (Taxaceae), Podocarpus nakaii (Podocarpaceae), Pinus luchuensis (Pinaceae), and Ginkgo biloba (Ginkgoales). Our phylogenetic analysis of the new sequence data with the published 18S rRNA sequence of Zamia pumila (a.cycad) as an outgroup strongly indicates that Taxus, Pinus, and Podocarpus form a monophyletic group with the exclusion of Ginkgo and that Taxus is more closely related to Pinus than to Podocarpus. Therefore, Taxaceae should be classified as a family of Coniferales. Our finding that Taxaceae, Pinaceae, and Podocarpaceae form a clade contradicts both the view that the uniovulate seed of Taxaceae is a primitive character and the view that the Taxaceae are descendants of the Podocarpaceae. Rather, the uniovulate seed of Taxaceae and that of some species of Podocarpus appear to have different origins, probably all reduced from multiovulate cones.

We have sequenced the mitochondrial cytochrome b gene from the guinea pig, the African porcupine, and a South American opossum. A phylogenetic analysis, which includes 22 eutherian and four other vertebrate cytochrome b sequences, indicates that the guinea pig and the porcupine constitute a natural clade (Hystricomorpha) that is not a sister group to the clade of mice and rats (Myomorpha). Therefore, the hypothesis that the Rodentia is paraphyletic receives additional support. The artiodactyls, the perissodactyls, and the cetaceans form a group that is separated from the primates and the rodents. The 26 sequences are used to study the structure/function relationships in cytochrome b, whose function is electron transport. Most of the amino acid residues involved in the two reaction centers are well conserved in evolution. The four histidines that are believed to ligate the two hemes are invariant among the 26 sequences, but their nearby residues are not well conserved in evolution. The eight transmembrane domains represent some of the most divergent regions in the cytochrome b sequence. The rate of nonsynonymous substitution is considerably faster in the human and elephant lineages than in other eutherian lineages; the faster rate might be due to coevolution between cytochrome b and cytochrome c.

The statistical properties of sample estimation and bootstrap estimation of phylogenetic variability from a sample of nucleotide sequences are studied by using model trees of three taxa with an outgroup and by assuming a constant rate of nucleotide substitution. The maximum-parsimony method of tree reconstruction is used. An analytic formula is derived for estimating the sequence length that is required if P, the probability of obtaining the true tree from the sampled sequences, is to be equal to or higher than a given value. Bootstrap estimation is formulated as a two-step sampling procedure: (1) sampling of sequences from the evolutionary process and (2) resampling of the original sequence sample. The probability that a bootstrap resampling of an original sequence sample will support the true tree is found to depend on the model tree, the sequence length, and the probability that a randomly chosen nucleotide site is an informative site. When a trifurcating tree is used as the model tree, the probability that one of the three bifurcating trees will appear in > or = 95% of the bootstrap replicates is < 5%, even if the number of bootstrap replicates is only 50; therefore, the probability of accepting an erroneous tree as the true tree is < 5% if that tree appears in > or = 95% of the bootstrap replicates and if more than 50 bootstrap replications are conducted. However, if a particular bifurcating tree is observed in, say, < 75% of the bootstrap replicates, then it cannot be claimed to be better than the trifurcating tree even if > or = 1,000 bootstrap replications are conducted. When a bifurcating tree is used as the model tree, the bootstrap approach tends to overestimate P when the sequences are very short, but it tends to underestimate that probability when the sequences are long. Moreover, simulation results show that, if a tree is accepted as the true tree only if it has appeared in > or = 95% of the bootstrap replicates, then the probability of failing to accept any bifurcating tree can be as large as 58% even when P = 95%, i.e., even when 95% of the samples from the evolutionary process will support the true tree. Thus, if the rate-constancy assumption holds, bootstrapping is a conservative approach for estimating the reliability of an inferred phylogeny for four taxa.

The statistical properties of sample estimation and bootstrap estimation of phylogenetic variability from a sample of nucleotide sequences were studied by considering model trees of three taxa with an outgroup. The cases of constant and varying rates of nucleotide substitution were compared. From sequences obtained by simulation, phylogenetic trees were constructed by using the maximum parsimony (MP) and neighbor-joining (NJ) methods. The effectiveness and consistency of the MP method were studied in terms of proportions of informative sites. The results of simulation showed that bootstrap estimation of the confidence level for an inferred phylogeny can be used even under unequal rates of evolution if the rate differences are not large so that the MP method is not misleading. The condition under which the MP method becomes misleading (inconsistent) is more stringent for slowly evolving sequences than for rapidly evolving ones, and it also depends on the length of the internal branch. If the rate differences are large so that the MP method becomes consistently misleading, then bootstrap estimation will reinforce an erroneous conclusion on topology. Similar conclusions apply to the NJ method with uncorrected distances. The NJ method with corrected distances performs poorly when the sequence length is short but can avoid the inconsistency problem if the sequence length is long and if the distances can be estimated accurately.

An analytical method is presented for constructing linear invariants. All linear invariants of a k-species tree can be derived from those of (k-1)-species trees using this method. The new method is simpler than that of Cavender, which relies on numerical computations. Moreover, the new method provides a convenient tool to study the relationships between linear invariants of the same tree or of different trees. All linear invariants of trees of up to five species are derived in this study. For four species, there are 16 independent linear invariants for each of the three possible unrooted trees, 14 of which are shared by two unrooted trees and 12 of these are shared by all three unrooted trees; the last types of linear invariants can be used to construct tests on the assumptions about nucleotide substitutions. The number of linear invariants for a tree is found to increase rapidly with the number of species.

1. Molecular data indicate that caviomorphs (guinea-pig-like rodents) and myomorphs (rat-like rodents) are not monophyletic. 2. Rather, the evolutionary lineage leading to the guinea-pig may have branched off prior to the divergence among myomorphs, lagomorphs, primates, chiropterans, artiodactyls, and carnivores. 3. Thus, the guinea-pig lineage probably represents an ancient eutherian lineage, and should be conferred an independent ordinal status. 4. The gundis (Ctenodactylidae) also seem to have branched off before the divergence among myomorphs, primates, and artiodactyls, but after the divergence of the guinea-pig. 5. Therefore, the order Rodentia as defined at the present time is in all probability a paraphyletic group devoid of taxonomic validity. 6. Previous claims pertaining to large differences in the rate of molecular evolution between guinea-pigs and myomorphs may have been exaggerated in many cases as a result of the erroneous phylogenetic position attributed to the guinea-pig. 7. The average rate of amino acid replacement in the guinea-pig is comparable to that in the rat and the mouse. 8. Protein-coding genes of myomorphs and caviomorphs evole, on average, about two times faster than their counterparts in gundis and humans.

The necessary and sufficient conditions for the existence of linear invariants under semigroups of probability transition matrices are derived. It is found that a biologically meaningful nucleotide substitution model has linear invariants if and only if it is a submodel of one of the three most general models, which include the so-called balanced and unbalanced transversion models. Each of these three general models is a nucleotide substitution model with six parameters.

Invariants are functions of the probabilities of state configurations among lineages, with expected values equal to zero under certain phylogenies. For two-state sequences, the existence of certain quadratic invariants requires a symmetric substitution model. For sequences with more than two states, the necessary condition for the existence of certain quadratic invariants in terms of independent events is much stronger than symmetry. For DNA sequences, only three parameters are allowed in the substitution model, which includes Kimura's two-parameter model as a special case.

The guinea-pig (Cavia porcellus), traditionally classified as a New World hystricomorph rodent, often shows anomalous morphological and molecular features in comparison with other eutherian mammals. For example, its insulin differs from that of other mammals in anabolic and growth-promoting activities and in its capability to form hexamers. Indeed, the literature about the molecular evolution of guinea-pigs abounds in references to 'convergent evolution', 'extremely rapid rates of substitution', and 'unique evolutionary mechanisms'. These claims are based on the assumption that the guinea-pig is a rodent. Our phylogenetic analyses of amino-acid sequence data, however, imply that the guinea-pig diverged before the separation of the primates and the artiodactyls from the myomorph rodents (rats and mice). If true, then the myomorphs and the caviomorphs do not constitute a natural clade, and the Caviomorpha (or the Histricomorpha) should be elevated in taxonomical rank and regarded as a separate mammalian order distinct from the Rodentia. If, as suggested by recent data, the myomorphs branched off before the divergence among the carnivores, lagomorphs, artiodactyls and primates, then the new order would represent an early divergence in eutherian radiation.

Phylogenetic analysis of DNA sequences from primates, rodents, lagomorphs, artiodactyls, carnivores, and birds strongly suggests that the order Rodentia is an outgroup to the other four mammalian orders and that Artiodactyla and Carnivora belong to a superordinal clade. Further, there is strong evidence against the Glires concept, which unites Lagomorpha and Rodentia. The radiation among Lagomorpha, Primates, and Artiodactyla--Carnivora is very bush-like, but there is some evidence that Lagomorpha has branched off first. Thus, the branching sequence for these five orders of mammals seems to be Rodentia, Lagomorpha, Primates, Artiodactyla, and Carnivora. The branching date for Rodentia could be as early as 100 million years ago. The rate of nucleotide substitution in the rodent lineage is shown to be at least 1.5 times higher than those in the other four mammalian lineages.

How many primary lineages of life exist and what are their evolutionary relationships? These are fundamental but highly controversial issues. Woese and co-workers propose that archaebacteria, eubacteria and eukaryotes are the three primary lines of descent and their relationships can be represented by Fig. 1a (the 'archaebacterial tree') if one neglects the root of the tree. In contrast, Lake claims that archaebacteria are paraphyletic, and he groups eocytes (extremely thermophilic, sulphur-dependent bacteria) with eukaryotes, and halobacteria with eubacteria (the 'eocyte tree', Fig. 1b). Lake's view has gained considerable support as a result of an analysis of small subunit ribosomal RNA sequence data by a new approach, the evolutionary parsimony method. Here we report that analysis of small subunit data by the neighbour-joining and maximum parasimony methods favours the archaebacterial tree and that computer simulations using either the archaebacterial or the eocyte tree as a model tree show that the probability of recovering the model tree is very high (greater than 90 per cent) for both the neighbour-joining and maximum parsimony methods but is relatively low for the evolutionary parsimony method. Moreover, analysis of large subunit rRNA sequences by all three methods strongly favours the archaebacterial tree.

A simple approach to testing the significance of the branching order, estimated from protein or DNA sequence data, of three taxa is proposed. The branching order is inferred by the transformed-distance method, under the assumption that one or two outgroups are available, and the branch lengths are estimated by the least-squares method. The inferred branching order is considered significant if the estimated internodal distance is significantly greater than zero. To test this, a formula for the variance of the internodal distance has been developed. The statistical test proposed has been checked by computer simulation. The same test also applies to the case of four taxa with no outgroup, if one considers an unrooted tree. Formulas for the variances of internodal distances have also been developed for the case of five taxa. Conditions are given under which it is more efficient to add the sequence of a fifth taxon than to do 25% more nucleotide sequencing in each of the original four. A method is presented for combining analyses of disparate data to get a single P value. Finally, the test, applied to the human-chimpanzee-gorilla problem, shows that the issue is not yet resolved.

Analysis of the rate of nucleotide substitution at silent sites in Drosophila genes reveals three main points. First, the silent rate varies (by a factor of two) among nuclear genes; it is inversely related to the degree of codon usage bias, and so selection among synonymous codons appears to constrain the rate of silent substitution in some genes. Second, mitochondrial genes may have evolved only as fast as nuclear genes with weak codon usage bias (and two times faster than nuclear genes with high codon usage bias); this is quite different from the situation in mammals where mitochondrial genes evolve approximately 5-10 times faster than nuclear genes. Third, the absolute rate of substitution at silent sites in nuclear genes in Drosophila is about three times higher than the average silent rate in mammals.

A simple method is proposed for constructing phylogenetic trees from distance matrices. The procedure for constructing tree topologies is similar to that of the unweighted pair-group method (UPG method) but makes corrections for unequal rates of evolution among lineages. The procedure for estimating branch lengths is the same as that of the Fitch and Margoliash method (F-M method) except that it allows no negative branch lengths. The performance of the present procedure for the construction of tree topologies is compared with that of the UPG method, the F-M method, Farris' method, and the modified Farris method by using Tateno's simulation outputs for nucleotide sequence divergence and his results for the performances of the latter four methods [Tateno, Y. (1978) Dissertation (Univ. Texas, Houston, TX). In this limited comparison, the present method performs considerably better than the UPG method and the F-M method and about equally well as the last two methods. The present method appears to be preferable to the UPG method for analysis of data from populations that have not differentiated much. Indeed, an application of the present method to gene frequency data from some Amerindian populations gives a tree topology far more reasonable than that obtained by the UPG method.

Mathematical formulas are developed for the evolutionary change of restriction cleavage sites in a DNA sequence, allowing unequal rates between transitional and transversional types of nucleotide substitution. Formulas are also developed for the probability of having a particular pattern of site changes among evolutionary lineages, such as parallel gains or losses of sites, and for inferring the presence or absence of a restriction site in an ancestral sequence from data on the present-day sequences. The unordered compatibility method is proposed for inferring the phylogenetic relationships among relatively closely related organisms, treating restriction sites as cladistic characters. Formulas are derived for the probability (P+) of obtaining the correct network for a given number (N) of informative sites for the cases of four and five species. These formulas are applied to evaluate the performance of the method and to estimate the N value required for P+ to be 95% or larger. The method performs well when the branches between ancestral nodes and the branches leading to the two most recent species are more or less equal in length, but performs poorly when the latter two branches are considerably longer than the former.

The branching order of the kingdoms Animalia, Plantae, and Fungi has been a controversial issue. Using the transformed distance method and the maximum parsimony method, we investigated this problem by comparing the sequences of several kinds of macromolecules in organisms spanning all three kingdoms. The analysis was based on the large-subunit and small-subunit ribosomal RNAs, 10 isoacceptor transfer RNA families, and six highly conserved proteins. All three sets of sequences support the same phylogenetic tree: plants and animals are sibling kingdoms that have diverged more recently than the fungi. The ribosomal RNA and protein data sets are large enough so that in both cases the inferred phylogeny is statistically significant. The present report appears to be the first to provide statistically conclusive molecular evidence for the phylogeny of the three kingdoms. The determination of this phylogeny will help us to understand the evolution of various molecular, cellular, and developmental characters shared by any two of the three kingdoms. Noting that the large-subunit rRNA sequences have evolved at similar rates in the three kingdoms, we estimated the ratio of the time since the animal-plant split to the time since the fungal divergence to be 0.90.