Analyses of mitochondrial DNA nest ratite birds within the Neognathae: supporting a neotenous origin of ratite morphological characters

It is commonly acknowledged that the basal divergence among extant birds is between Palaeognathae and Neognathae. However, recent analyses of mitochondrial sequences have challenged that notion. In order to investigate this fundamental phylogenetic question, the complete mitochondrial DNA (mtDNA) molecule of the rook Corvus frugilegus (Passeriformes) was sequenced and included in phylogenetic analyses with the previously reported complete mtDNAs of the chicken Gallus gallus (Galliformes) and two ratite species, the ostrich Struthio camelus and the rhea Rhea americana (Struthioniformes). The analyses reconstructed a split between Passeriformes and a branch including Galliformes and Struthioniformes. Thus, the result is inconsistent with the traditional understanding of a basal avian divergence between Palaeognathae and Neognathae. The findings suggest that the morphological characteristics of the ratites are secondarily acquired, probably through neoteny and that the ratites are descendants of flying, neognathous ancestors.     

The mtDNA sequence of the ostrich and the divergence between paleognathous and neognathous birds

The complete mitochondrial DNA (mtDNA) molecule of the ostrich, Struthio camelus, was sequenced. The size of the molecule is 16,591 nucleotides. Since the ostrich represents the paleognathous birds, comparison with the mtDNA of the neognathous chicken, the only avian species reported so far in databases, made it possible to identify common and, probably, general avian mtDNA characteristics. Relative to other vertebrates, the avian NADH6 and tRNA-Glu genes are positioned upstream of the control region rather than the cytochrome b gene. The NADH3 gene of the ostrich is terminated by a stop codon at position 207. Thus, the gene is about 140 nucleotides shorter than in other vertebrates. The sequence for L-strand origin of replication is missing in both birds, and four transfer RNA genes of the two avian mtDNAs deviate from common characteristics of tRNAs of vertebrate mtDNAs by having an adenine (and not a thymidine) at position 8. Due to the absence of suitable fossils, most paleontological datings of avian divergences are conjectural. Molecular dating of the divergence between the ostrich and the chicken indicates that these two avian lineages separated 80-90 MYA. Phylogenetic analysis of complete cytochrome b genes of six avian orders showed that Passeriformes represent the earliest divergence among recent birds, contradicting the commonly accepted notion of a basal position of the Palaeognathae among recent birds.      

Mitochondrial DNA of the coral sarcophyton glaucum contains a gene for a homologue of bacterial muts: A possible case of gene transfer from the nucleus to the mitochondrion

The nucleotide sequences of two segments of 6,737 ntp and 258 ntp of the 18.4-kb circular mitochondrial (mt) DNA molecule of the soft coral Sarcophyton glaucum (phylum Cnidaria, class Anthozoa, subclass Octocorallia, order Alcyonacea) have been determined. The larger segment contains the 3′ 191 ntp of the gene for subunit 1 of the respiratory chain NADH dehydrogenase (ND1), complete genes for cytochrome b (Cyt b), ND6, ND3, ND4L, and a bacterial MutS homologue (MSH), and the 5′ terminal 1,124 ntp of the gene for the large subunit rRNA (l-rRNA). These genes are arranged in the order given and all are transcribed from the same strand of the molecule. The smaller segment contains the 3′ terminal 134 ntp of the ND4 gene and a complete tRNA^f-Met gene, and these genes are transcribed in opposite directions. As in the hexacorallian anthozoan, Metridium senile, the mt-genetic code of S. glaucum is near standard: that is, in contrast to the situation in mt-genetic codes of other invertebrate phyla, AGA and AGG specify arginine, and ATA specifies isoleucine. However, as appears to be universal for metazoan mt-genetic codes, TGA specifies tryptophan rather than termination. Also, as in M. senile the mt-tRNA^f-Met gene has primary and secondary structural features resembling those of Escherichia coli initiator tRNA, including standard dihydrouridine and TψC loop sequences, and a mismatched nucleotide pair at the top of the amino-acyl stem. The presence of a mutS gene homologue, which has not been reported to occur in any other known mtDNA, suggests that there is mismatch repair activity in S. glaucum mitochondria. In support of this, phylogenetic analysis of MutS family protein sequences indicates that the S. glaucum mtMSH protein is more closely related to the nuclear DNA-encoded mitochondrial mismatch repair protein (MSH1) of the yeast Saccharomyces cerevisiae than to eukaryotic homologues involved in nuclear function, or to bacterial homologues. Regarding the possible origin of the S. glaucum mtMSH gene, the phylogenetic analysis results, together with comparative base composition considerations, and the absence of an MSH gene in any other known mtDNA best support the hypothesis that S. glaucum mtDNA acquired the mtMSH gene from nuclear DNA early in the evolution of octocorals. The presence of mismatch repair activity in S. glaucum mitochondria might be expected to influence the rate of evolution of this organism's mtDNA. Received: 13 January 1997 / Accepted: 23 September 1997     

Molecular Timing of Primate Divergences as Estimated by Two Nonprimate Calibration Points

The complete mitochondrial DNA (mtDNA) molecule of the hamadryas baboon, Papio hamadryas, was sequenced and included in a molecular analysis of 24 complete mammalian mtDNAs. The particular aim of the study was to time the divergence between Cercopithecoidea and Hominoidea. That divergence, set at 30 million years before present (MYBP) was a fundamental reference for the original proposal of recent hominoid divergences, according to which the split among gorilla, chimpanzee, and Homo took place 5 MYBP. In the present study the validity of the postulated 30 MYBP dating of the Cercopithecoidea/Hominoidea divergence was examined by applying two independent nonprimate molecular references, the divergence between artiodactyls and cetaceans set at 60 MYBP and that between Equidae and Rhinocerotidae set at 50 MYBP. After calibration for differences in evolutionary rates, application of the two references suggested that the Cercopithecoidea/Hominoidea divergence took place >50 MYBP. Consistent with the marked shift in the dating of the Cercopithecoidea/Hominoidea split, all hominoid divergences receive a much earlier dating. Thus the estimated date of the divergence between Pan (chimpanzee) and Homo is 10–13 MYBP and that between Gorilla and the Pan/Homo linage ≈17 MYBP. The same datings were obtained in an analysis of clocklike evolving genes. The findings show that recalculation is necessary of all molecular datings based directly or indirectly on a Cercopithecoidea/Hominoidea split 30 MYBP. Received: 1 April 1998 / Accepted: 1 July 1998     

Evolution of the Mitochondrial Cytochrome Oxidase II Gene in Collembola

The sequence of the mitochondrial COII gene has been widely used to estimate phylogenetic relationships at different taxomonic levels across insects. We investigated the molecular evolution of the COII gene and its usefulness for reconstructing phylogenetic relationships within and among four collembolan families. The collembolan COII gene showed the lowest A + T content of all insects so far examined, confirming that the well-known A + T bias in insect mitochondrial genes tends to increase from the basal to apical orders. Fifty-seven percent of all nucleotide positions were variable and most of the third codon positions appeared free to vary. Values of genetic distance between congeneric species and between families were remarkably high; in some cases the latter were higher than divergence values between other orders of insects. The remarkably high divergence levels observed here provide evidence that collembolan taxa are quite old; divergence levels among collembolan families equaled or exceeded divergences among pterygote insect orders. Once the saturated third-codon positions (which violated stationarity of base frequencies) were removed, the COII sequences contained phylogenetic information, but the extent of that information was overestimated by parsimony methods relative to likelihood methods. In the phylogenetic analysis, consistent statistical support was obtained for the monophyly of all four genera examined, but relationships among genera/families were not well supported. Within the genus Orchesella, relationships were well resolved and agreed with allozyme data. Within the genus Isotomurus, although three pairs of populations were consistently identified, these appeared to have arisen in a burst of evolution from an earlier ancestor. Isotomurus italicus always appeared as basal and I. palustris appeared to harbor a cryptic species, corroborating allozyme data. Received: 12 January 1996 / Accepted: 10 August 1996     

The Mitochondrial Genome of the Sperm Whale and a New Molecular Reference for Estimating Eutherian Divergence Dates

Extant cetaceans are systematically divided into two suborders: Mysticeti (baleen whales) and Odontoceti (toothed whales). In this study, we have sequenced the complete mitochondrial (mt) genome of an odontocete, the sperm whale (Physeter macrocephalus), and included it in phylogenetic analyses together with the previously sequenced complete mtDNAs of two mysticetes (the fin and blue whales) and a number of other mammals, including five artiodactyls (the hippopotamus, cow, sheep, alpaca, and pig). The most strongly supported cetartiodactyl relationship was: outgroup,((pig, alpaca),((cow, sheep),(hippopotamus,(sperm whale,(baleen whales))))). As in previous analyses of complete mtDNAs, the sister-group relationship between the hippopotamus and the whales received strong support, making both Artiodactyla and Suiformes (pigs, peccaries, and hippopotamuses) paraphyletic. In addition, the analyses identified a sister-group relationship between Suina (the pig) and Tylopoda (the alpaca), although this relationship was not strongly supported. The paleontological records of both mysticetes and odontocetes extend into the Oligocene, suggesting that the mysticete and odontocete lineages diverged 32–34 million years before present (MYBP). Use of this divergence date and the complete mtDNAs of the sperm whale and the two baleen whales allowed the establishment of a new molecular reference, O/M-33, for dating other eutherian divergences. There was a general consistency between O/M-33 and the two previously established eutherian references, A/C-60 and E/R-50. Cetacean (whale) origin, i.e., the divergence between the hippopotamus and the cetaceans, was dated to ≈55 MYBP, while basal artiodactyl divergences were dated to ≥65 MYBP. Molecular estimates of Tertiary eutherian divergences were consistent with the fossil record. Received: 12 July 1999 / Accepted: 28 February 2000     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

Male-Driven Evolution Among Eoaves? A Test of the Replicative Division Hypothesis in a Heterogametic Female (ZW) System

Because avian females are heterogametic, the reverse of mammals, avian sex chromosomes undergo significantly different patterns and numbers of DNA replications than do those in mammals. This makes the W (female-specific) and the Z chromosomes an excellent model system for the study of the replicative division hypothesis, which purports that DNA substitution rate is determined by the number of germline replications. The sex-specific chromosome in birds (the W) is predicted to change at the slowest rate of all avian chromosomes because it undergoes the fewest rounds of replication per unit of evolutionary time. Using published data on gametogenesis from a variety of sources, we estimated the ratio of male-to-female germline replications (c) in galliforms and anseriforms to be approximately 4.4. The value of c should predict the value of the ratio of male-to-female mutation rates (α_m) if the replicative division hypothesis is true. Homologous DNA sequences including an intron and parts of two exons of the CHD gene were obtained from the W and the Z chromosomes in ostrich, sage grouse, canvasback duck, tundra swan, and snow goose. The exons show significantly different nucleotide composition from the introns, and the W-linked exons show evidence of relaxed constraint. The Z-linked intron is diverging ≈ 3.1 times faster than the W-linked intron. From this, α_m was calculated to be approximately 4.1, with a confidence interval of 3.1 to 5.1. The data support the idea that the number of replicative divisions is a major determinant of substitution rate in the Eoavian genome. Received: 19 January 1999 / Accepted: 8 June 1999     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

The mitochondrial DNA molecule of sumatran orangutan and a molecular proposal for two (Bornean and Sumatran) species of orangutan

The complete mitochondrial DNA (mtDNA) molecule of Sumatran orangutan, plus the complete mitochondrial control region of another Sumatran specimen and the control regions and five protein-coding genes of two specimens of Bornean orangutan were sequenced and compared with a previously reported complete mtDNA of Bornean orangutan. The two orangutans are presently separated at the subspecies level. Comparison with five different species pairs—namely, harbor seal/grey seal, horse/donkey, fin whale/blue whale, common chimpanzee/pygmy chimpanzee, and Homo/common chimpanzee—showed that the molecular difference between Sumatran and Bornean orangutan is much greater than that between the seals, and greater than that between the two chimpanzees, but similar to that between the horse and the donkey and the fin and blue whales. Considering their limited morphological distinction the comparison revealed unexpectedly great molecular difference between the two orangutans. The nucleotide difference between the orangutans is about 75% of that between Homo and the common chimpanzee, whereas the amino acid difference exceeds that between Homo and the common chimpanzee. On the basis of their molecular distinction we propose that the two orangutans should be recognized as different species, Pongo pygmaeus, Bornean orangutan, and P. abelii, Sumatran orangutan. Received: 15 May 1996 / Accepted: 21 June 1996     

The Role of Yeast VDAC Genes on the Permeability of the Mitochondrial Outer Membrane

In addition to the POR1 gene, which encodes the well-characterized voltage dependent anion-selective channel (YVDAC1) of the mitochondrial outer membrane, the yeast Saccharomyces cerevisiae contains a second gene (POR2) encoding a protein (YVDAC2) with 50% sequence identity to YVDAC1. Mitochondria isolated from yeast cells deleted for the POR1 gene (Δpor1) had a profoundly reduced outer membrane permeability as measured by the ability of an intermembrane space dehydrogenase to oxidize exogenously added NADH. Mitochondria missing either YVDAC1 or both YVDAC1 and YVDAC2 showed a 2-fold increase in the rate of NADH oxidation when the outer membrane was deliberately damaged. Mitochondria from parental cells showed only a 10% increase indicating that the outer membrane is highly permeable to NADH. In the absence of YVDAC1, we calculate that the outer membrane permeability to NADH is reduced 20-fold. The low NADH permeability in the presence of YVDAC2 was not due to the low levels of YVDAC2 expression as mitochondria from cells expressing levels of YVDAC2 comparable to those of YVDAC1 in parental cells showed no substantial increase in NADH permeability, indicating a minimal role of YVDAC2 in this permeability. The residual permeability may be due to other pathways because cells missing both genes can still grow on nonfermentable carbon sources. However, YVDAC1 is clearly the major pathway for NADH flux through the outer membrane in these mitochondria. Received: 23 May 1997/Revised: 3 October 1997     

Comparative Genomics of Mitochondrial DNA in Members of the Drosophila melanogaster Subgroup

In this study, a comparative genomics approach is employed to investigate the forces that shape evolutionary change in the mitochondrial DNA (mtDNA) of members of the Drosophila melanogaster subgroup. This approach facilitates differentiation of the patterns of variation resulting from processes acting at a higher level from those acting on a single gene. The mitochondrial genomes of three isofemale lines of D. simulans (siI, -II, and -III), two of D. melanogaster (Oregon R and a line from Zimbabwe), and D. mauritiana (maI and -II), and one of D. sechellia were sequenced and compared with that derived from D. yakuba. Data presented here indicate that at least three broad mechanisms shape the evolutionary dynamics of mtDNA in these taxa. The first set of mechanisms is intrinsic to the molecule. Dominant processes may be interpreted as selection for an increased rate of replication of the mtDNA molecule, biases in DNA repair, and differences in the pattern of nucleotide substitution among strands. In the genes encoded on the major strand (62% of the coding DNA) changes to or from C predominate, whereas on the minor changes to or from G predominate. The second set of mechanisms affects distinct lineages. There are evolutionary rate differences among lineages, possibly owing to population demographic changes or changes in mutational biases. This is supported by the heterogeneity found in synonymous, nonsynonymous, and silent substitutions. The third set of mechanisms differentially affects distinct genes. A maximum-likelihood sliding-window analysis detected four disjunct regions that have a significantly different nucleotide substitution process from that derived from the complete sequence. These data show the potential for comparative genomics to tease apart subtle forces that shape the evolution of DNA. Received: 30 July 1999 / Accepted: 16 March 2000     

The Complete Mitochondrial Genome of the Hagfish Myxine glutinosa: Unique Features of the Control Region

The complete sequence of the mitochondrial DNA of the hagfish Myxine glutinosa has been determined. The hagfish mtDNA (18,909 bp) is the longest vertebrate mtDNA determined so far. The gene arrangement conforms to the consensus vertebrate type and differs from that of lampreys. The exceptionally long (3628-bp) control region of the hagfish contains the typical conserved elements found in other vertebrate mtDNAs but is characterized by a large number of putative hairpins, which can potentially fold into a highly compact secondary structure that appears to be unique to hagfish. The comparison of the mtDNAs of two M. glutinosa specimens, excluding the control region, shows a 0.6% divergence at the nucleotide level as a sample of intraspecies polymorphism. Received: 21 August 2000 / Accepted: 2 March 2001     

Evolution of a Metal-Binding Cluster in the NH2-Terminal Variable Region of Avian Fast Skeletal Muscle Troponin T: Functional Divergence on the Basis of Tolerance to Structural Drifting

This study investigated the evolution of a transition metal ion-binding cluster ([H–X–X–X–H] _n ; Tx) in the alternatively spliced NH₂-terminal variable region of avian pectoral muscle troponin T (TnT). Encoded by avian fast skeletal muscle TnT-specific P exons, Tx-like structures were expressed in the breast muscle TnT's of almost all birds examined. Their presence results in the developmentally up-regulated high molecular weight pectoral muscle TnT. Sequence analysis and metal affinity chromatography revealed that in Galliformes and Craciformes, the Tx structure evolved into multiple H–X–X–X–H pairs with a high-affinity metal-binding capacity. Turkey, chicken, quail, and curassow breast muscle TnT's contain nine, seven, four, and three consecutive or closely located metal-binding sites, respectively, in the NH₂-terminal region. The metal-binding affinity of the Tx element increased as the number of His pairs increased due to the duplication of P exons and the conversion of other exon sequences. The data show two related components of avian pectoral muscle TnT evolution: a larger, more acidic NH₂-terminal segment and a cluster of transition metal-binding sites, both of which may have functional significance for their selection value. The evolution of the Tx segment in the NH₂-terminal variable region of avian pectoral muscle TnT demonstrates a functional divergence on the basis of tolerance to structural drifting. Received: 2 May 2000 / Accepted: 5 September 2000     

The complete mitochondrial DNA (mtDNA) of the donkey and mtDNA comparisons among four closely related mammalian species-pairs

The nucleotide sequence of the complete mitochondrial genome of the donkey, Equus asinus, was determined. The length of the molecule is 16,670 bp. The length, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of two types of repetitive motifs in the control region. The sequence of the repeats is (a) 5′-CACACCCA and (b) 5′-TGCGCGCA, respectively. The order of (a) and (b) can be expressed as {n[2(a)+(b)]+m(a)}. In 32 different clones analyzed the number of n and m ranged from 0 to 9 and 1 to 7. The two rRNA genes, the 13 peptide-coding genes, and the 22 tRNA genes of the donkey and the horse, Equus caballus, were compared in detail. Total nucleotide difference outside the control region was 6.9%. Nucleotide difference between peptide-coding genes ranged from 6.4% to 9.4% with a mean of 8.0%. In the inferred protein sequences of the 13 peptide-coding genes the amino acid difference was 0.2–8.8%, and the mean for the 13 concatenated amino acid sequences was 1.9%. In the 22 tRNA genes, the mean difference was 3.5%, and that in the two rRNA genes was 4.1%. The mtDNA differences between the donkey and the horse suggest that the evolutionary separation of the two species occurred ≈9 million years ago. Analyses of differences among the mtDNAs of three other species-pairs, harbor seal/grey seal, fin whale/blue whale, and Homo/common chimpanzee, showed that the relative evolutionary rate of individual peptide-coding genes varies among different species-pairs and modes of comparison. The findings show that the superimposition of sequence data of one lineage for resolving and dating evolutionary divergences of other lineages should be performed with caution unless based on comprehensive data. Received: 15 October 1995 / Accepted: 15 April 1996     

An Estimate of Divergence Time of Parazoa and Eumetazoa and That of Cephalochordata and Vertebrata by Aldolase and Triose Phosphate Isomerase Clocks

Previously we suggested that four proteins including aldolase and triose phosphate isomerase (TPI) evolved with approximately constant rates over long periods covering the whole animal phyla. The constant rates of aldolase and TPI evolution were reexamined based on three different models for estimating evolutionary distances. It was shown that the evolutionary rates remain essentially unchanged in comparisons not only between different classes of vertebrates but also between vertebrates and arthropods and even between animals and plants, irrespective of the models used. Thus these enzymes might be useful molecular clocks for inferring divergence times of animal phyla. To know the divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata, the aldolase cDNAs from Ephydatia fluviatilis, a freshwater sponge, and the TPI cDNAs from Ephydatia fluviatilis and Branchiostoma belcheri, an amphioxus, have been cloned and sequenced. Comparisons of the deduced amino acid sequences of aldolase and TPI from the freshwater sponge with known sequences revealed that the Parazoa–Eumetazoa split occurred about 940 million years ago (Ma) as determined by the average of two proteins and three models. Similarly, the aldolase and TPI clocks suggest that vertebrates and amphioxus last shared a common ancestor around 700 Ma and they possibly diverged shortly after the divergence of deuterostomes and protostomes.     

Ancient Origin of the Null Allele se 428 of the Human ABO-Secretor Locus (FUT2)

In human populations, a null allele having several nucleotide differences from the wild-type allele is segregating at the FUT2 locus (the ABO-Secretor locus) encoding α(1,2)fucosyltransferase. To estimate the age of the most recent common ancestor (MRCA) of these two alleles, we sequenced FUT2 homologues from chimpanzee, gorilla, orangutan, and green monkey. Since we did not detect acceleration or any heterogeneity in the substitution rate at this locus among these species, the age of the MRCA was estimated to be around 3 MYA, assuming the divergence time of human and chimpanzee to be 5 MYA. We developed a simple test to examine whether or not the old age of the MRCA of the FUT2 is consistent with that expected for two divergent neutral alleles sampled from a random mating population. An application of the test to the data at FUT2 indicated that the age of the MRCA is too old to be explained by the simple neutral assumptions, although our test depends on accurate estimation of the divergence time of human and chimpanzee in units of twice the human population size. Various possibilities including balancing selection are discussed to explain this old age of the MRCA. Received: 9 May 1999 / Accepted: 20 September 1999     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998