The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Mitochondrial DNA Phylogeny of the Family Cichlidae: Monophyly and Fast Molecular Evolution of the Neotropical Assemblage

A mitochondrial DNA (mtDNA) phylogeny of cichlid fish is presented for the most taxonomically inclusive data set compiled to date (64 taxa). 16S rDNA data establish with confidence relationships among major lineages of cichlids, with a general pattern congruent with previous morphological studies and less inclusive molecular phylogenies based on nuclear genes. Cichlids from Madagascar and India are the most basal groups of the family Cichlidae and sister to African–Neotropical cichlids. The cichlid phylogeny suggests drift-vicariance events, consistent with the fragmentation of Gondwana, to explain current biogeographic distributions. Important phylogenetic findings include the placement of the controversial genus Heterochromis basal among African cichlids, the South American genus Retroculus as the most basal taxon of the Neotropical cichlid assemblage, and the close relationship of the Neotropical genera Cichla with Astronotus rather than with the crenicichlines. Based on a large number of South American genera, the Neotropical cichlids are defined as a monophyletic assemblage and shown to harbor significantly higher levels of genetic variation than their African counterparts. Relative rate tests suggest that Neotropical cichlids have experienced accelerated rates of molecular evolution. But these high evolutionary rates were significantly higher among geophagine cichlids. Received: 18 September 1998 / Accepted: 16 December 1998     

Characterization and Evolution of the Mitochondrial DNA Control Region in Hornbills (Bucerotiformes)

We determined the mitochondrial DNA control region sequences of six Bucerotiformes. Hornbills have the typical avian gene order and their control region is similar to other avian control regions in that it is partitioned into three domains: two variable domains that flank a central conserved domain. Two characteristics of the hornbill control region sequence differ from that of other birds. First, domain I is AT rich as opposed to AC rich, and second, the control region is approximately 500 bp longer than that of other birds. Both these deviations from typical avian control region sequence are explainable on the basis of repeat motifs in domain I of the hornbill control region. The repeat motifs probably originated from a duplication of CSB-1 as has been determined in chicken, quail, and snowgoose. Furthermore, the hornbill repeat motifs probably arose before the divergence of hornbills from each other but after the divergence of hornbills from other avian taxa. The mitochondrial control region of hornbills is suitable for both phylogenetic and population studies, with domains I and II probably more suited to population and phylogenetic analyses, respectively.     

Ancient Mitochondrial DNA Reveals the Origin of Sus scrofa from Rebun Island, Japan

The Kabukai A site (5 to 8C A.D.) of the Okhotsk cultural area is on Rebun Island, a small island near the coast, north–northwest of Hokkaido, Japan. Specimens of Sus scrofa, called the Sakhalin pig, were discovered in five cultural layers at the Kabukai A site. Ancient DNA was extracted from the remains of 42 Sakhalin pig bones. Thirty-nine nucleotide sequences of the 574-bp mitochondrial DNA control region, estimated to have originated from at least 21 individuals, were amplified and analyzed phylogenetically. Nine distinct haplotypes (A1, A2, A3, B1, B2, C1, C2, D1, and D2) from this site were classified into four haplotype groups (A, B, C, and D) by parsimonious network analysis. Phylogenetic analysis of 9 ancient and 55 modern haplotypes indicated that the population of Sakhalin pigs at the Kabukai A site belonged to two distinct clusters; haplotype groups A and B formed a cluster comprised only of themselves, and haplotype groups C and D belonged to the cluster of one of the two genetic groups of Japanese wild boars uniquely distributed in the western part of Japan, including one northeast Mongolian wild boar. Analysis of the haplotype distribution among three archaeological sites and their historical transitions among the five layers reflecting the cultural periods at the Kabukai A site suggests that the Sakhalin pig populations were introduced from Sakhalin island and the Amur River basin in the northeastern Eurasian continent together with some cultural influences. Received: 18 April 2000 / Accepted: 24 November 2000     

The Molecular Evolution of the Vertebrate Trypsinogens

We expand the already large number of known trypsinogen nucleotide and amino acid sequences by presenting additional trypsinogen sequences from the tunicate (Boltenia villosa), the lamprey (Petromyzon marinus), the pufferfish (Fugu rubripes), and the frog (Xenopus laevis). The current array of known trypsinogen sequences now spans the entire vertebrate phylogeny. Phylogenetic analysis is made difficult by the presence of multiple isozymes within species and rates of evolution that vary highly between both species and isozymes. We nevertheless present a Fitch-Margoliash phylogeny constructed from pairwise distances. We employ this phylogeny as a vehicle for speculation on the evolution of the trypsinogen gene family as well as the general modes of evolution of multigene families. Unique attributes of the lamprey and tunicate trypsinogens are noted. Received: 12 July 1997     

Phylogeographic Patterns and Evolution of the Mitochondrial DNA Control Region in Two Neotropical Cats (Mammalia, Felidae)

The ocelot (Leopardus pardalis) and margay (L. wiedii) are sister-species of Neotropical cats which evolved from a lineage that migrated into South America during the formation of the Panamanian land bridge 3–5 million years ago. Patterns of population genetic divergence of each species were studied by phylogenetic analyses of mitochondrial DNA (mtDNA) control region sequences in individuals sampled across the distribution of these taxa. Abundant genetic diversity and remarkably concordant phylogeographic partitions for both species were observed, identifying parallel geographic regions which likely reflect historical faunal barriers. Inferred aspects of phylogeography, population genetic structure, and demographic history were used to formulate conservation recommendations for these species. In addition, observed patterns of sequence variation provided insight into the molecular evolution of the mtDNA control region in closely related felids. Received: 26 January 1998 / Accepted: 14 May 1998     

Prehistoric Introduction of Domestic Pigs onto the Okinawa Islands: Ancient Mitochondrial DNA Evidence

Ancient DNAs of Sus scrofa specimens excavated from archaeological sites on the Okinawa islands were examined to clarify the genetic relationships among prehistoric Sus scrofa, modern wild boars and domestic pigs inhabiting the Ryukyu archipelago, the Japanese islands, and the Asian continent. We extracted remain DNA from 161 bone specimens excavated from 12 archaeological sites on the Okinawa islands and successfully amplified mitochondrial DNA control region fragments from 33 of 161 specimens. Pairwise difference between prehistoric and modern S. scrofa nucleotide sequences showed that haplotypes of the East Asian domestic pig lineage were found from archaeological specimens together with Ryukyu wild boars native to the Ryukyu archipelago. Phylogenetic analysis of 14 ancient sequences (11 haplotypes; 574 bp) indicated that S. scrofa specimens from two Yayoi-Heian sites (Kitahara and Ara shellmiddens) and two Recent Times sites (Wakuta Kiln and Kiyuna sites) are grouped with modern East Asian domestic pigs. Sus scrofa specimens from Shimizu shellmidden (Yayoi-Heian Period) were very closely related to modern Sus scrofa riukiuanus but had a unique nucleotide insertion, indicating that the population is genetically distinct from the lineage of modern Ryukyu wild boars. This genetic evidence suggests that domestic pigs from the Asian continent were introduced to the Okinawa islands in the early Yayoi-Heian period (1700–2000 BP), or earlier.     

Molecular Phylogenetic Inference of the Woolly Mammoth Mammuthus primigenius, Based on Complete Sequences of Mitochondrial Cytochrome b and 12S Ribosomal RNA Genes

Complete sequences of cytochrome b (1,137 bases) and 12S ribosomal RNA (961 bases) genes in mitochondrial DNA were successfully determined from the woolly mammoth (Mammuthus primigenius), African elephant (Loxodonta africana), and Asian elephant (Elephas maximus). From these sequence data, phylogenetic relationships among three genera were examined. Molecular phylogenetic trees reconstructed by the neighbor-joining and the maximum parsimony methods provided an identical topology both for cytochrome b and 12S rRNA genes. These results support the ``Mammuthus-Loxodonta' clade, which is contrary to some previous morphological reports that Mammuthus is more closely related to Elephas than to Loxodonta. Received: 8 April 1997 / Accepted: 23 July 1997     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

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     

Circular Permutations in the Molecular Evolution of DNA Methyltransferases

Circular permutations of genes during molecular evolution often are regarded as elusive, although a simple model can explain these rearrangements. The model assumes that first a gene duplication of the precursor gene occurs in such a way that both genes become fused in frame, leading to a tandem protein. After generation of a new start codon within the 5′ part of the tandem gene and a stop at an equivalent position in the 3′ part of the gene, a protein is encoded that represents a perfect circular permutation of the precursor gene product. The model is illustrated here by the molecular evolution of adenine-N⁶ DNA methyltransferases. β- and γ-type enzymes of this family can be interconverted by a single circular permutation event. Interestingly, tandem proteins, proposed as evolutionary intermediates during circular permutation, can be directly observed in the case of adenine methyltransferases, because some enzymes belonging to type IIS, like the FokI methyltransferase, are built up by two fused enzymes, both of which are active independently of each other. The mechanism for circular permutation illustrated here is very easy and applicable to every protein. Thus, circular permutation can be regarded as a normal process in molecular evolution and a changed order of conserved amino acid motifs should not be interpreted to argue against divergent evolution. Received: 17 November 1998 / Accepted: 19 February 1999     

Genetic Structure of the Ancestral Population of Modern Humans

Neutral DNA polymorphisms from an 8-kb segment of the dystrophin gene, previously ascertained in a worldwide sample (n= 250 chromosomes), were used to characterize the population ancestral to the present-day human groups. The ancestral state of each polymorphic site was determined by comparing human variants with their orthologous sites in the great apes. The ``age before fixation' of the underlying mutations was estimated from the frequencies of the new alleles and analyzed in the context of these polymorphisms' distribution among 13 populations from Africa, Europe, Asia, New Guinea, and the Americas (n= 860 chromosomes in total). Seventeen polymorphisms older tan 100,000–200,000 years, which contributed ∼90% to the overall nucleotide diversity, were common to all human groups. Polymorphisms endemic to human groups or continentally restricted were younger than 100,000–200,000 years. Africans (six populations) with 13 such sites stood out from the rest of the world (seven populations), where only 2 population-specific variants were observed. The similarity of the frequencies of the old polymorphisms in Africans and non-Africans suggested a similar profile of genetic variability in the population before the modern human's divergence. This ancestral population was characterized by an effective size of about 10,000 as estimated from the nucleotide diversity; this size may describe the number of breeding individuals over a long time during the Middle Pleistocene or reflect a speciation bottleneck from an initially larger population at the end of this period. Received: 3 February 1998 / Accepted: 9 February 1998     

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     

The Complete Mitochondrial DNA Sequence of the Pig (Sus scrofa)

The complete mitochondrial genome sequence of the pig, Sus scrofa, was determined. The length of the sequence presented is 16,679 nucleotides. This figure is not absolute, however, due to pronounced heteroplasmy caused by variable numbers of the motif GTACACGTGC in the control region of different molecules. A phylogenetic study was performed on the concatenated amino acid and nucleotide sequences of 12 protein-coding genes of the mitochondrial genome. The analysis identified the pig (Suiformes) as a sister group of a cow/whale clade, making Artiodactyla paraphyletic. The split between pig and cow/whale was molecularly dated at 65 million years before present. Received: 2 December 1997 / Accepted: 20 February 1998     

Evolution of 28S Ribosomal DNA in Chaetognaths: Duplicate Genes and Molecular Phylogeny

The chaetognaths are an extraordinarily homogeneous phylum of animals at the morphological level, with a bauplan that can be traced back to the Cambrian. Despite the attention of zoologists for over two centuries, there is little agreement on classification within the phylum. We have used a molecular biological approach to investigate the phylogeny of extant chaetognaths. A rapidly evolving expansion segment toward the 5′ end of 28S ribosomal DNA (rDNA) was amplified using the polymerase chain reaction (PCR), cloned, and sequenced from 26 chaetognath samples representing 18 species. An unusual finding was the presence of two distinct classes of 28S rDNA gene in chaetognaths; our analyses suggest these arose by a gene (or gene cluster) duplication in a common ancestor of extant chaetognaths. The two classes of chaetognath 28S rDNA have been subject to different rates of molecular evolution; we present evidence that both are expressed and functional. In phylogenetic reconstructions, the two classes of 28S rDNA yield trees that root each other; these clearly demonstrate that the Aphragmophora and Phragmophora are natural groups. Within the Aphragmophora, we find good support for the groupings denoted Solidosagitta, Parasagitta, and Pseudosagitta. The relationships between several well-supported groups within the Aphragmophora are uncertain; we suggest this reflects rapid, recent radiation during chaetognath evolution. Received: 19 March 1996 / Accepted: 5 August 1996     

Evolution of Human Polyomavirus JC: Implications for the Population History of Humans

The polyomavirus JC virus (JCV), the etiological agent of progressive multifocal leukoencephalopathy, is ubiquitous in the human population, infecting children asymptomatically, then persisting in the kidney. The main mode of transmission of JCV is from parents to children through long-term cohabitation. Twelve JCV subtypes that occupy unique domains in Europe, Africa, and Asia have been identified. Here, we attempted to elucidate the evolutionary relationships among JCV strains worldwide using the whole-genome approach with which a highly reliable phylogeny of JCV strains can be reconstructed. Sixty-five complete JCV DNA sequences, derived from various geographical regions and belonging to 11 of the 12 known subtypes, were subjected to phylogenetic analysis using three independent methods: the neighbor-joining, maximum parsimony, and maximum likelihood methods. The trees obtained with these methods consistently indicated that ancestral JCVs were divided into three superclusters, designated as Types A, B, and C. A split in Type A generated two subtypes, EU-a and -b, mainly containing European and Mediterranean strains. The first split in Type B generated Af2 (the major African subtype). Subsequent splits in Type B generated B1-c (a minor European subtype) and all seven Asian subtypes (B1-a, -b, -d, B2, MY, CY, and SC). Type C generated a single subtype (Af1), consisting of strains derived from western Africa. While the present findings provided a basis on which to classify JCV into types or subtypes, they have several implications for the divergence and migration of human populations. Received: 4 April 2001 / Accepted: 31 July 2001     

Mutation and Recombination in Cattle Satellite DNA: A Feedback Model for the Evolution of Satellite DNA Repeats

The cattle genome contains several distinct centromeric satellites with interrelated evolutionary histories. We compared these satellites in Bovini species that diverged 0.2 to about 5 Myr ago. Quantification of hybridization signals by phosphor imaging revealed a large variation in the relative amounts of the major satellites. In the genome of water buffalo this has led to the complete deletion of satellite III. Comparative sequencing and PCR-RFLP analysis of satellites IV, 1.711a, and 1.711b from the related Bos and Bison species revealed heterogeneities in 0.5 to 2% of the positions, again with variations in the relative amounts of sequence variants. Restriction patterns generated by double digestions suggested a recombination of sequence variants. Our results are compatible with a model of the life history of satellites during which homogeneity of interacting repeat units is both cause and consequence of the rapid turnover of satellite DNA. Initially, a positive feedback loop leads to a rapid saltatory amplification of homogeneous repeat units. In the second phase, mutations inhibit the interaction of repeat units and coexisting sequence variants amplify independently. Homogenization by the spreading of one of the variants is prevented by recombination and the satellite is eventually outcompeted by another, more homogeneous tandem repeat sequence. Received: 21 July 2000 / Accepted: 30 October 2000