Showing papers in "Molecular Biology and Evolution in 1992"

Population growth makes waves in the distribution of pairwise genetic differences.

Abstract: Episodes of population growth and decline leave characteristic signatures in the distribution of nucleotide (or restriction) site differences between pairs of individuals. These signatures appear in histograms showing the relative frequencies of pairs of individuals who differ by i sites, where i = 0, 1, .... In this distribution an episode of growth generates a wave that travels to the right, traversing 1 unit of the horizontal axis in each 1/2u generations, where u is the mutation rate. The smaller the initial population, the steeper will be the leading face of the wave. The larger the increase in population size, the smaller will be the distribution's vertical intercept. The implications of continued exponential growth are indistinguishable from those of a sudden burst of population growth Bottlenecks in population size also generate waves similar to those produced by a sudden expansion, but with elevated uppertail probabilities. Reductions in population size initially generate L-shaped distributions with high probability of identity, but these converge rapidly to a new equilibrium. In equilibrium populations the theoretical curves are free of waves. However, computer simulations of such populations generate empirical distributions with many peaks and little resemblance to the theory. On the other hand, agreement is better in the transient (nonequilibrium) case, where simulated empirical distributions typically exhibit waves very similar to those predicted by theory. Thus, waves in empirical distributions may be rich in information about the history of population dynamics.

Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C-content biases.

Abstract: A simple mathematical method is developed to estimate the number of nucleotide substitutions per site between two DNA sequences, by extending Kimura's (1980) two-parameter method to the case where a G+C-content bias exists. This method will be useful when there are strong transition-transversion and G+C-content biases, as in the case of Drosophila mitochondrial DNA.

A Simple Method for Estimating and Testing Minimum-Evolution Trees

Abstract: A simple method for estimating and testing phylogenetic trees under the principle of minimum evolution (ME) is presented. The basic procedure of this method is first to obtain the neighbor-joining (NJ) tree by Saitou and Nei’s method and then to search for a tree with the minimum value of the sum (S) of branch lengths by examining all trees that are closely related to the NJ tree. Once the ME tree is identified, a statistical test is conducted for the difference in S between this tree and other closely related trees. The mathematical method required for conducting this test is developed by using the least-squares approach. Computer simulation has shown that this method identifies the correct tree with a high probability, as long as the number of nucleotides examined is sufficiently large. It has also been shown that the topology of the NJ tree is almost always identical with that of the ME tree. A method for obtaining least-squares estimates (and their standard errors) of branch lengths for a given topology is also presented. This method can be used for testing the reliability of the branching pattern of the ME tree. However, the statistical test of S values is more powerful in rejecting incorrect trees than is the branch-length test or bootstrapping. Furthermore, both a mathematical method for computing the number of trees with a given value of topological difference from the NJ tree and a computer algorithm for identifying all the topologies are developed.

A statistical test for detecting geographic subdivision.

Abstract: A statistical test for detecting genetic differentiation of subpopulations is described that uses molecular variation in samples of DNA sequences from two or more localities. The statistical significance of the test is determined with Monte Carlo simulations. The power of the test to detect genetic differentiation in a selectively neutral Wright-Fisher island model depends on both sample size and the rates of migration, mutation, and recombination. It is found that the power of the test is substantial with samples of size 50, when 4Nm less than 10, where N is the subpopulation size and m is the fraction of migrants in each subpopulation each generation. More powerful tests are obtained with genes with recombination than with genes without recombination.

The number of replications needed for accurate estimation of the bootstrap P value in phylogenetic studies.

Abstract: The bootstrap is a statistical method for obtaining a nonparametric estimate of error (Efron 1979, 1982). Felsenstein ( 1985) was the first to apply this method to phylogeny estimation, and his approach is now widely used. Taxa are held constant, and the characters (for sequence data, nucleotide or amino acid sites) are resampled randomly with replacement. A phylogeny is constructed from each replication of the data, and the frequency of appearance of particular phylogenetic groups (groups of alleles or taxa) among all of the trees constructed by this resampling is the bootstrap confidence limit, or bootstrap P value (BP). The BPS of different nodes within a tree can be used to assess the relative stability of those phylogenetic groups or, if applied strictly, to test their statistical significance (e.g., at the 95% or 99% level). The application of bootstrapping to phylogeny estimation is a tradeoff between the maximum number of replications that can be performed by the researcher in a reasonable amount of time and the minimum number of replications needed for accurate estimation of the BP. The purpose of the present report is to explore the variance (and hence the accuracy) of the phylogenetic BP and to establish guidelines for efficient bootstrap sampling. BP is the proportion of trees containing a particular phylogenetic group. It therefore follows the binomial distribution, which has a variance of o2 = [P( 1 P)/n], where P is the BP and II is the number of replications. Although Li and Gouy ( 1990) recently suggested that more replications are needed for larger numbers of taxa, the accuracy of the BP is a function only of P and n. If the interval containing 95% of the samples (+ 1.96 standard deviations) is used as a measure of accuracy, then the application of the above formula shows that 1,825 replications [ =0.95 X 0.05( 1.96/ 0.0 1 )2] are needed to attain + 1% accuracy at a BP of 0.95 (fig. 1). This is more than an order of magnitude higher than the number of replications (50-100) normally used in phylogenetic analyses. Based on this, a practical guideline for efficient and accurate bootstrap sampling can be made: If one wishes the expectation to be that the 95% confidence range is + 1% of the BP, then one must perform 2,000 bootstrap replications (if BP = 0.95) or 400 replications (if BP = 0.99) in phylogenetic analyses, unless the computational time is prohibitive; additional replications are unnecessary, and fewer replications may sacrifice statistical accuracy. Moreover, statistical testing at the 95% level cannot be made using ~73 replications, even if the group is supported by a BP of 1 .OO. This is because the inaccuracy at a mean BP of 0.95 is greater than f5% (fig. 1) when ~73 replications are used. In other words, a BP > 1 .OO could be obtained when the actual (mean) BP is not significant ( cO.95). Thus, the 50 replications used by Felsenstein ( 1985) in his original example and the 20100 replications used in many subsequent studies (e.g., see Ovenden et al. 1987; Thomas et al. 1989; Jansen et al. 1990; Meyer and Wilson 1990; Douglas et al. 199 1; Irwin et al. 199 1) would appear to be far too few for the intended purpose of statistical testing at the 95% level.

Mitochondrial DNA evolution at a turtle's pace: evidence for low genetic variability and reduced microevolutionary rate in the Testudines.

Abstract: Evidence is compiled suggesting a slowdown in mean microevolutionary rate for turtle mitochondrial DNA (mtDNA). Within each of six species or species complexes of Testudines, representing six genera and three taxonomic families, sequence divergence estimates derived from restriction assays are consistently lower than expectations based on either (a) the dates of particular geographic barriers with which significant mtDNA genetic clades appear associated or (b) the magnitudes of sequence divergence between mtDNA clades in nonturtle species that otherwise exhibit striking phylogeographic concordance with the genetic partitions in turtles. Magnitudes of the inferred rate slowdowns average eightfold relative to the "conventional" mtDNA clock calibration of 2%/Myr sequence divergence between higher animal lineages. Reasons for the postulated deceleration remain unknown, but two intriguing correlates are (a) the exceptionally long generation length most turtles and (b) turtles' low metabolic rate. Both factors have been suspected of influencing evolutionary rates in the DNA sequences of some other vertebrate groups. Uncertainities about the dates of cladogenetic events in these Testudines leave room for alternatives to the slowdown interpretation, but consistency in the direction of the inferred pattern, across several turtle species and evolutionary settings, suggests the need for caution in acceptance of a universal mtDNA-clock calibration for higher animals.

Phylogeny of five fungus-like protoctistan Phytophthora species, inferred from the internal transcribed spacers of ribosomal DNA.

Abstract: Ribosomal DNA variation was used to study evolutionary relationships among five fungal-like protoctistan Phytophthora species. On the basis of morphological and ecological characteristics, four of these species--P. palmivora, P. megakarya, P. capsici, and P. citrophthora--were once thought to be related. Variation within a species was extensively studied in a fifth, outgroup species--P. cinnamomi--known, on the basis of ecological, isozyme, and mitochondrial DNA studies, to be variable. Internal transcribed spacer regions (ITS I, between the 18S and 5.8S rDNAs; and ITS II, between the 5.8S and 25S rDNAs) from 27 isolates of these five species were analyzed by polymerase chain reaction amplification and direct sequencing. Intraspecific variability was undetected or low. Interspecific nucleotide difference was 0.3%-14.6%, and comparisons of variable regions permitted the evaluation of phylogenetic relationships among species. Both neighbor-joining and parsimony analysis of ITS variability support a close relationship between cacao isolates of P. capsici and P. citrophthora and a common lineage for P. palmivora and P. megakarya. Large distance values were estimated between P. cinnamomi and the other species. Inferred relationships based on ITS variability were compared with those based on other characters. The catalog of sequences provides the information necessary to design taxon-specific probes potentially useful in taxonomic, ecological, and population-level studies.

Statistical properties of bootstrap estimation of phylogenetic variability from nucleotide sequences. I. Four taxa with a molecular clock.

Abstract: 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.

Molecular evolution of the fungi: relationship of the Basidiomycetes, Ascomycetes, and Chytridiomycetes.

Abstract: Establishing the phylogeny of fungi and protists often has proved difficult owing to the simple morphologies and convergent characters in these organisms. We used DNA sequences of nuclear small-subunit ribosomal RNA genes to determine phylogenetic relationships among three major classes of organisms considered to be fungi-Basidiomycetes, Ascomycetes and Chytridiomycetes-and to assess the taxonomic position of Neocallimastix, an economically important anaerobic rumen microorganism whose classification is controversial. The Basidiomycetes and Ascomycetes, two classes of nonflagellated fungi, are the most closely related taxa. Chytridiomycetes, though bearing flagella, group with these higher fungi rather than with the protists. Neocallimastix, a eukaryote lacking mitochondria and variously classified as a protist or as a fungus, shows closest molecular affinities with the Chytridiomycete fungi in the order Spizellomycetales.

Two ascomycete classes based on fruiting-body characters and ribosomal DNA sequence.

Abstract: Traditional fruiting body-based classification of ascomycetes has been under attack for 2 decades. Fruiting-body types can converge, and few researchers now assume that either the closed fruiting bodies (cleistothecia) characterizing the class Plectomycetes or the flask-shaped fruiting bodies (perithecia) characterizing the class Pyrenomycetes are stable, unifying characters. Unless we identify characters uniting major ascomycete groups, orders of ascomycetes remain narrowly defined, and supraordinal classification is impossible. We sequenced both strands of 18s rDNA from nine ascomycete fungi, adding three sequences from GenBank into our analysis. The phylogeny, inferred from 162 informative sites in 1,700 bp of DNA sequence data and using yeast as an outgroup, divided the fungi into two groups correlating well both with fruiting-body type and with the traditional classes Plectomycetes and Pyrenomycetes. Each group received strong statistical support. Genera producing cleistothecia, such as Talaromyces (with a Penicillium asexual state) and the human pathogen Ajellomyces capsulatus (causing histoplasmosis), fall within the plectomycete group. Plectomycetes also includes Eremascus albus and the bee pathogen Ascosphaera apis, although both lack typical fruiting bodies. The Dutch elm disease fungus groups with pyrenomycetes such as Neurospora, in spite of its confusing mixture of class-level characters.

Positive selection and interallelic recombination at the merozoite surface antigen-1 (MSA-1) locus of Plasmodium falciparum.

Abstract: DNA sequences of alleles at the merozoite surface antigen-1 (MSA-1) gene locus of the malaria parasite Plasmodium falciparum show evidence of repeated past recombination events between alleles. These include both (1) nonreciprocal recombination events that have homogenized certain gene regions among alleles and (2) reciprocal recombination events that have combined allelic segments with divergent evolutionary histories, thereby enhancing allelic diversity. In three different gene regions, the rate of nonsynonymous nucleotide substitution significantly exceeds that of synonymous nucleotide substitution, implying that positive Darwinian selection has acted to diversify alleles at the amino acid level. The MSA-1 polymorphism seems to be quite ancient; the two major allelic types have been maintained for approximately 35 Myr.

Mitochondrial DNA phylogeny of the Old-World monkey tribe Papionini.

Abstract: The evolution of the Old World monkey tribe Papionini, composed of macaques, baboons, mandrills, drills, and mangabeys, was examined using mitochondrial DNA (mtDNA) sequence data on the cytochrome ox&se subunit II gene. When analyzed cladisticahy, these data support a baboon clade of savannah (Pupio) plus gelada ( Theropithecus) baboons, as well as a clade containing drill (Mandrillus) plus mangabey (Cercocebus) genera. This result stands in opposition to most morphological phylogenies, which break up the baboon clade by placing Papio and Mandrillus as sister taxa and Theropithecus as a more distantly related lineage. Analyses of CO11 gene sequences also suggest that the papionin ancestral stock divided into two lineages, one leading to macaques and the other to the purely African genera. From a molecular evolutionary perspective, the papionin CO11 gene sequences reveal a pattern of amino acid replacements concentrated in the regions spanning the mitochondrial membrane.

Rapid surveying of DNA sequence variation in natural populations.

Abstract: DNA sequencing can be costly and time consuming for population studies because of the relative rarity of variation along exons. These problems can be substantially reduced by the use of the polymerase chain reaction on introns using primers from the exon region. These problems can be further reduced by the use of denaturing gradient gel electrophoresis to identify those alleles in need of sequencing.

Phylogenetic position of phylum Nemertini, inferred from 18S rRNA sequences: molecular data as a test of morphological character homology.

Abstract: Partial 18S rRNA sequence of the nemertine Cerebratulus lacteus was obtained and compared with those of coelomate metazoans and acoelomate platyhelminths to test whether nemertines share a most recent common ancestor with the platyhelminths, as traditionally has been implied, or whether nemertines lie within a protostome coelomate clade, as suggested by more recent morphological analyses. Maximum-parsimony analysis supports the inclusion of the nemertine within a protostome-coelomate clade that falls within a more inclusive coelomate clade. Bootstrap analysis indicates strong support for a monophyletic Coelomata composed of a deuterostome and protostome-coelomate clade. Support for a monophyletic protostome Coelomata is weak. Inference by distance analysis is consistent with that of maximum parsimony. Analysis of down-weighted paired sites by maximum parsimony reveals variation in topology only within the protostome-coelomate clade. The relationships among the protostome coelomates cannot be reliably inferred from the partial sequences, suggesting that coelomate protostomes diversified rapidly. Results with evolutionary parsimony are consistent with the inclusion of the nemertine in a coelomate clade. The molecular inference corroborates recent morphological character analyses that reveal no synapomorphies of nemertines and flatworms but instead suggest that the circulatory system and rhynchocoel of nemertines are homologous to coelomic cavities of protostome coelomates, thus supporting the corresponding hypothesis that nemertines belong within a protostome-coelomate clade. The sequence data provide an independent test of morphological character homology.

The rate and pattern of nucleotide substitution in Drosophila mitochondrial DNA.

Abstract: The nucleotide sequences of a segment of mitochondrial DNA (mtDNA) have been determined for nine species or subspecies of the subgenus Drosophila of the genus Drosophila. This segment contains two complete protein-coding genes (i.e., NADH dehydrogenase subunit 1 and cytochrome b) and a transfer RNA gene (tRNA(ser)). The G+C content at third-codon positions for the two protein-coding genes was 1.5 times higher than that in the D. melanogaster species group, which belongs to the subgenus Sophophora. However, there was a substantial difference between the nucleotide frequencies of G and C. The number of nucleotide substitutions per silent site was more than three times higher than that for nuclear DNA, although it was only 60% of that for mammalian mtDNA. Both parametric and nonparametric analyses revealed a strong transition-transversion bias in nucleotide substitution, as was observed in mammalian mtDNA. Moreover, the rate of substitution of A and T for G and C is higher than that for the opposite direction. This bias seems to be responsible for the extremely A+T-rich base composition of Drosophila mtDNA. It is also noted that the rate of transitional change between A and G is higher than that between T and C.

Evolutionary adaptation to different thermal environments via transcriptional regulation.

Abstract: Populations of the teleost fish Fundu E.C. 1.1.1.27 ). Previously we have shown that the compensating change in the LDH-B enzyme concentration is due to a change in the amount of LDH-B mRNA, but we did not define whether this was due to differences in mRNA stability or to differences in rate of transcription. The results presented here help clarify the molecular mechanism responsible for the variation in Ldh-B gene expression: the rate of transcription from the Ldh-B locus is significantly different between populations, and this difference is responsible for the compensatory change in LDH-B enzyme concentration.

Molecular phylogeny of the volvocine flagellates.

Abstract: Phylogenetic studies of approximately 2,000 bases of sequence from the large and small nuclear-encoded ribosomal RNAs are used to investigate the origins of the genus Volvox. The colonial and multicellular genera currently placed in the family Volvocaceae form a monophyletic group that is significantly closer phylogenetically to Chlamydomonas reinhardtii than it is to the other unicellular green flagellates that were tested, including Chlamydomonas eugametos, Chlorella pyrenoidosa, and Haematococcus lacustris. Statistical analysis of 251 phylogenetically informative nucleotide positions rejects the "volvocine lineage" hypothesis, which postulates a monophyletic evolutionary progression from unicellular organisms (such as Chlamydomonas), through colonial organisms (e.g., Gonium, Pandorina, Eudorina, and Pleodorina) demonstrating increasing size, cell number, and tendency toward cellular differentiation, to multicellular organisms having fully differentiated somatic and reproductive cells (in the genus Volvox). The genus Volvox appears not to be monophyletic. Volvox capensis falls outside a lineage containing other representatives of Volvox (V. aureus, V. carteri, and V. obversus), and both of these Volvox lineages are more closely related to certain colonial genera than they are to each other. This implies either a diphyletic origin of Volvox from different colonial volvocacean ancestors, a phylogenetic derivation of some of the colonial genera from a multicellular (i.e., Volvox) ancestor, or both. Considered together with previously published observations, these results suggest that the different levels of organizational and developmental complexity found in the Volvocaceae represent alternative stable states, among which evolutionary transitions have occurred several times during the phylogenetic history of this group.

Rate and mode differences between nuclear and mitochondrial small-subunit rRNA genes in mushrooms.

Abstract: Sequences from homologous regions of the nuclear and mitochondrial small-subunit rRNA genes from 10 members of the mushroom order Boletales were used to construct evolutionary trees and to compare the rates and modes of evolution. Trees constructed independently for each gene by parsimony and tested by bootstrap analysis have identical topologies in all statistically significant branches. Examination of base substitutions revealed that the nuclear gene is biased toward C-T transitions and that the distribution of transversions in the mitochondrial gene is strongly effected by an A-T bias. When only homologous regions of the two genes were compared, base substitutions per nucleotide were roughly 16-fold greater in the mitochondrial gene. The difference in the frequency of length mutations was at least as great but was impossible to estimate accurately because of their absence in the nuclear gene. Maximum likelihood was used to show that base-substitution rates vary dramatically among the branches. A significant part of the rate inconstancy was caused by an accelerated nuclear rate in one branch and a retarded mitochondrial rate in a different branch. A second part of the rate variability involved a consistent inconstancy: short branches exhibit ratios of mitochondrial to nuclear divergences of less than 1, while longer branches had ratios of approximately 4:1-8:1. This pattern suggests a systematic error in the branch length calculation. The error may be related to the simplicity of the divergence estimates, which assumes that all base positions have an equal probability of change.

Estimation of heterozygosity for single-probe multilocus DNA fingerprints.

Abstract: In spite of the increasing application of DNA fingerprinting to natural populations and to the genetic identification of humans, explicit methods for estimation of basic population genetic parameters from DNA fingerprinting data have not been developed. Contributing to this omission is the inability to determine, for multilocus fingerprinting probes, relatively important genetic information, such as the number of loci, the number of alleles, and the distribution of these alleles into specific loci. One of the most useful genetic parameters that could be derived from such data would be the average heterozygosity, which has traditionally been employed to measure the level of genetic variation within populations and to compare genetic variation among different loci. We derive here explicit formulas for both the estimation of average heterozygosity at multiple hypervariable loci and a maximum value for this estimate. These estimates are based upon the DNA restriction-pattern matrices that are typical for fingerprinting studies of humans and natural populations. For several empirical data sets from our laboratory, estimates of average and maximal heterozygosity are shown to be relatively close to each other. Furthermore, variances of these statistics based on simulation studies are relatively small. These observations, as well as consideration of the effect of missing alleles and alternate numbers of loci, suggest that the average heterozygosity can be accurately estimated using phenotypic DNA fingerprint patterns, because this parameter is relatively insensitive to the lack of certain genetic information.

Molecular evolution of the fungi: human pathogens.

Abstract: The morphological, ecological, and clinical diversity among ascomycete fungi that are pathogenic to humans suggest that the potential for pathogenicity may have arisen multiple times within these higher fungi. We have obtained 18S ribosomal DNA sequences from a diverse group of human pathogenic fungi in order to determine their evolutionary origins. The fungi studied include a skin pathogen that is confined to humans (Trichophyton rubrum) and three systemic, facultative parasites that cause histoplasmosis (Histoplasma capsulatum), blastomycosis (Blastomyces dermatitidis) and coccidioidomycosis (Coccidioides immitis) in humans and other higher animals. Also included in our analysis are representatives of non-pathogenic fungi, as well as two opportunistic pathogens, Pneumocystis carinii and Candida albicans, that cause severe disease in immunocompromised individuals, especially those with AIDS. Two of the fungi we sequenced, T. rubrum and C. immitis, are limited to asexual modes of reproduction and therefore lack the sexual structures that are most useful for evolutionary comparison as well as being essential for classification among the higher fungi. Coccidioides immitis is particularly problematic owing to its contradictory and confusing asexual morphologies, which have caused it to be placed in three fungal classes and the protista. Our analysis shows that the specialized, superficial parasite and the systemic, facultative parasites, including C. immitis, are closely related ascomycetes, which clearly demonstrates the power of molecular characters to compensate for missing or confusing reproductive morphology. Analysis also shows that the opportunistic pathogens are more distantly related, with the likely explanation that pathogenicity has arisen more than once within the Ascomycetes.

Mitochondrial DNA-like sequence in the nuclear genome of an akodontine rodent.

Abstract: Initial amplification and sequencing of a 366-bp fragment of the cytochrome b gene by a conserved primer pair (MVZ 03 and MVZ 04) revealed a nonfunctional copy of the gene with two deletions (one of which is 17 bp in length and the other of which is 3 bp in length) in Chroeomys jelskii, a South American akodontine rodent. By means of an alternative primer to MVZ 03--namely, MVZ 05--from the region of the tRNA for glutamic acid, a functional copy of cytochrome b was subsequently amplified. Both primer pairs amplify functional sequence when applied to purified mitochondrial DNA (mtDNA). Restriction-endonuclease digestion of purified mtDNA from C. jelskii did not reveal any additional sets of bands that would suggest heteroplasmy in the mitochondrial genome. When probed with both functional and nonfunctional gene fragments, MboI restriction digests revealed the same pattern, providing further evidence that the nonfunctional copy must be located in the nucleus. Observed differences in the mitochondrial and nuclear sequences from two populations are consistent with a faster rate of change in mtDNA than in nuclear DNA.

Relative-Rate Test for Nucleotide Substitutions between Two Lineages

Abstract: The molecular-clock hypothesis, i.e., the rate constancy of molecular evolution (Zuckerkandl and Pauling 1965), has been a cornerstone for several issues in molecular evolution, such as the neutral mutation hypothesis (Kimura 1983) and the estimation of phylogenies (Nei 1987, pp. 309-3 13 ) and divergence times (Li and Graur 199 1, pp. 117-l 18). For the purpose of assessing the constancy of molecular evolution, a relative-rate test was developed to determine whether protein replacements are the same in two taxa (Sarich and Wilson 1973), and this test has been extended to nucleotide substitutions ( Wu and Li 1985; Li and Tanimura 1987). In the latter test, the evolutionary distance between taxon 1 and a reference taxon is compared with that between taxon 2 and the reference taxon, so that only two taxa could be simultaneously compared. However, it is of interest in some cases to compare numbers of nucleotide substitutions between two lineages, with each lineage including several taxa. One solution is to perform a series of relative-rate tests between taxa in a pairwise fashion, with one taxon from each of the two lineages. These tests, however, are not independent, so that statistical inferences drawn are not clear. A more rational solution would be to construct a test by using all the sequences from the two lineages. In this note, we extend the relative-rate test of Wu and Li ( 1985) to cases where more than one taxon are sampled for one nucleotide sequence in two lineages compared. Suppose that a phylogenetic tree consists of two lineages, with lineage 1 having n taxa and with lineage 2 having m taxa (fig. 1). First, we present formulas for calculating both the weighted number of substitutions in each lineage as compared with a common reference taxon and the variance of the difference in the weighted number of nucleotide substitutions between the two lineages. Then we assume that the number of nucleotide substitutions follows a Poisson process, so that the standardized normal distribution could be used to test whether the numbers of nucleotide substitutions are the same in the two lineages ( Wu and Li 1985). We define the following terms: N, = the number of nucleotides compared between the ith sequence in lineage 1 and the reference sequence, i = 1, . . . , n. Nj = the number of nucleotides compared between the jth sequence in lineage 2 and the reference sequence, j = 1, . . . , m.

Nucleotide sequence, function, activation, and evolution of the cryptic asc operon of Escherichia coli K12.

Abstract: The cryptic asc (previous called "SAC") operon of Escherichia coli K12 has been completely sequenced. It encodes a repressor (ascG); a PTS enzyme IIasc for the transport of arbutin, salicin, and cellobiose (ascF); and a phospho-beta-glucosidase that hydrolyzes the sugars which are phosphorylated during transport (ascB). ascG and ascFB are transcribed from divergent promoters. The cryptic operon is activated by the insertion of IS186 into the ascG (repressor) gene. The ascFB genes are paralogous to the cryptic bglFB genes, and ascG is paralogous to galR. The duplications that gave rise to these paralogous genes are estimated to have occurred approximately 320 Mya, a time that predates the divergence of E. coli and Salmonella typhimurium.

Nucleotide sequence variation in the mitochondrial 12S rRNA gene and the phylogeny of African mole-rats (Rodentia: Bathyergidae).

Abstract: Mitochondrial DNA (mtDNA) sequence variation was examined in eight taxa of the African rodent family Bathyergidae, as well as in two taxa representative of the Old-World hystricognathid rodent families Petromyidae and Thryonomyidae. A total of 812 bp, constituting domains I-III of the 12S ribosomal rRNA gene, were compared for each taxon. The high levels of intrafamilial mtDNA sequence divergence observed (average 16.8, range 3.5-23.2) support an ancient origin for the five genera, 20-38 Mya. These data do not support the current subfamilial groupings of the Bathyergidae. The eastern African naked mole-rat, Heterocephalus glaber, is the most basal representative of the family, with the silvery mole-rat, Heliophobius, being the next most basal. South African forms [dune, common, and cape mole-rats (Bathyergus, Cryptomys, and Georychus, respectively)] group together. The independent origin of the common mole-rat, relative to the naked mole-rat, suggests that complex social systems evolved in parallel along different bathyergid lineages. The 12S rRNA gene is not evolving at a higher rate within the rodent lineages, relative to that seen for artiodactyls and primates. Bathyergid rodents appear to fall at an extreme end of the spectrum of mammalian variation, with respect to both transition/transversion ratios and divergence, showing much lower transition/transversion ratios than those previously reported for intrafamilial comparisons.

Phylogeny of the Bovidae (Artiodactyla, Mammalia), based on mitochondrial ribosomal DNA sequences.

Abstract: Portions of the 12S and 16S mitochondrial ribosomal genes for 16 species representing nine tribes in the mammal family Bovidae were compared with six previously published orthologous sequences. Phylogenetic analysis of variable nucleotide positions under different constraints and weighting schemes revealed no robust groupings among tribes. Consensus trees support previous hypotheses of monophyly for four clades, including the traditional subfamily Bovinae. However, the basal diversification of bovid tribes, which was largely unresolved by morphological, immunodiffusion, allozyme, and protein sequence data, remains unresolved with the addition of DNA sequence data. The intractability of this systematic problem is consistent with a rapid radiation of the major bovid groups. Several analyses of our data show that monophyly of the Bovidae, which was weakly supported by previous morphological and molecular work, is questionable.

Complete congruence between morphological and rbcL-based molecular phylogenies in birches and related species (Betulaceae).

Abstract: Estimations of phylogenies from morphological and molecular data often show contrasting results. We compared morphological and molecular phylogenies in an ancient family of woody dicots, the Betulaceae (birch family). The phylogeny of the family was estimated from parsimony analysis of morphological characters in the genera Alnus, Betula, Carpinus, Corylus, Ostrya, and Ostryopsis and from parsimony and distance-matrix analyses of DNA sequences of the chloroplast gene encoding the large subunit of ribulose- 1 ,Sbiphosphate carboxylase ( rbcL) in the genera Alnus, Betula, Carpinus, Corylus, and Ostrya and in two outgroups, Quercus and Liquidambar. The topologies obtained by the different methods were completely congruent, and bootstrapping strongly supported the division of the family Betulaceae into two major clades, Betuleae (Alnus and Betula) and Coryleae (other members). Only slightly more homoplasy was present in the rbcL sequence data set than in the morphological set. Relative-rate tests indicated that the Coryleae clade had a faster rate of rbcL evolution than did the Betuleae clade. Heterogeneity of rates of morphological evolution also paralleled those for rbcL.

Freeing phylogenies from artifacts of alignment.

Abstract: Widely used methods for phylogenetic inference, both those that require and those that produce alignments, share certain weaknesses. These weaknesses are discussed, and a method that lacks them is introduced. For each pair of sequences in the data set, the method utilizes both insertion-deletion and amino acid replacement information to estimate a pairwise evolutionary distance. It is also possible to allow regional heterogeneity of replacement rates. Because a likelihood framework is adopted, the standard deviation of each pairwise distance can be estimated. The distance matrix and standard error estimates are used to infer a phylogenetic tree. As an example, this method is used on 10 widely diverged sequences of the second largest RNA polymerase subunit. A pseudo-bootstrap technique is devised to assess the validity of the inferred phylogenetic tree.

The consistency of several phylogeny-inference methods under varying evolutionary rates.

Abstract: A phylogenetic method is a consistent estimator of phylogeny if and only if it is guaranteed to give the correct tree, given that sufficient (possibly infinite) independent data are examined. The following methods are examined for consistency: UPGMA (unweighted pair-group method, averages), NJ (neighbor joining), MF (modified Farris), and P (parsimony). A two-parameter model of nucleotide sequence substitution is used, and the expected distribution of character states is calculated. Without perfect correction for superimposed substitutions, all four methods may be inconsistent if there is but one branch evolving at a faster rate than the other branches. Partial correction of observed distances improves the robustness of the NJ method to rate variation, and perfect correction makes the NJ method a consistent estimator for all combinations of rates that were examined. The sensitivity of all the methods to unequal rates varies over a wide range, so relative-rate tests are unlikely to be a reliable guide for accepting or rejecting phylogenies based on parsimony analysis.

Sequencing errors and molecular evolutionary analysis.

Abstract: Heuristic approaches were used to quantify the influence that sequencing errors have on estimates of nucleotide diversity, substitution rate, and the construction of genealogies. Error rates of less than 1 nucleotide/kb probably have little affect on conclusions about the evolutionary history of highly polymorphic organisms such as Drosophila and Escherichia coli, but organisms with very low nucleotide diversity, such as humans, require greater sequencing accuracy. A scan of GenBank for corrections of previous errors reveals that sequencing errors are highly nonrandom.

Horizontal transfer of hobo transposable elements within the Drosophila melanogaster species complex: evidence from DNA sequencing.

Abstract: The hobo family of transposable elements, one of three transposable-element families that cause hybrid dysgenesis in Drosophila melanogaster, appears to be present in all members of the D. melanogaster species complex: D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Some hobo-hybridizing sequences are also found in the other members of the melanogaster subgroup and in many members of the related montium subgroup. Surveys of older isofemale lines of D. melanogaster suggest that complete hobo elements were absent prior to 50 years ago and that hobo has recently been introduced into the species by horizontal transfer. To test the horizontal transfer hypothesis, the 2.6-kb XhoI fragments of hobo elements from D. melanogaster, D. simulans, and D. mauritiana were cloned and sequenced. The DNA sequences reveal an extremely low level of divergence and support the conclusion that the active hobo element has been horizontally transferred into or among these species in the recent past.