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     

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     

The Mitochondrial Genome of Chlorogonium elongatum Inferred from the Complete Sequence

The 22,704-bp circular mitochondrial DNA (mtDNA) of the chlamydomonad alga Chlorogonium elongatum was completely cloned and sequenced. The genome encodes seven proteins of the respiratory electron transport chain, subunit 1 of the cytochrome oxidase complex (cox1), apocytochrome b (cob), five subunits of the NADH dehydrogenase complex (nad1, nad2, nad4, nad5, and nad6), a set of three tRNAs (Q, W, M), and the large (LSU)- and small (SSU)-subunit ribosomal RNAs. Six group-I introns were found, two each in the cox1, cob, and nad5 genes. In each intron an open reading frame (ORF) related to maturases or endonucleases was identified. Both the LSU and the SSU rRNA genes are split into fragments intermingled with each other and with other genes. Although the average A + T content is 62.2%, GC-rich clusters were detected in intergenic regions, in variable domains of the rRNA genes, and in introns and intron-encoded ORFs. A comparison of the genome maps reveals that C. elongatum and Chlamydomonas eugametos mtDNAs are more closely related to one another than either is to Chlamydomonas reinhardtii mtDNA. Received: 3 November 1997 / Accepted: 12 January 1998     

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     

A Phylogenetic Study of the Origin of the Domestic Pig Estimated from the Near-Complete mtDNA Genome

The near-complete pig mtDNA genome sequence (15,997 bp) was determined from two domestic pigs (one Chinese Meishan and one Swedish Landrace) and two European wild boars. The sequences were analyzed together with a previously published sequence representing a Swedish domestic pig. The sequences formed three distinct clades, denoted A, E1, and E2, with considerable sequence divergence between them (0.8–1.2%). The results confirm our previous study (based on the sequence of the cytochrome B gene and the control region only) and provide compelling evidence that domestication of pigs must have occurred from both an Asian and a European subspecies of the wild boar. We estimated the time since the divergence of clade A (found in Chinese Meishan pigs) and E1 (found in European domestic pigs) at about 900,000 years before present, long before domestication about 9000 years ago. The pattern of nucleotide substitutions among the sequences was in good agreement with previous interspecific comparisons of mammalian mtDNA; the lowest substitution rates were observed at nonsynonymous sites in protein-coding genes, in the tRNA and rRNA genes, while the highest rates were observed at synonymous sites and in the control region. The presence of Asian clade A in some major European breeds (Large White and Landrace) most likely reflects the documented introgression of Asian germplasm into European stocks during the 18th and 19th centuries. The coexistence of such divergent mtDNA haplotypes for 100+ generations is expected to lead to the presence of recombinant haplotypes if paternal transmission and recombination occur at a low frequency. We found no evidence of such recombination events in the limited sample studied so far. Received: 19 April 2000; Accepted: 15 November 2000     

Monophyletic Origin of the Order Chiroptera and Its Phylogenetic Position Among Mammalia, as Inferred from the Complete Sequence of the Mitochondrial DNA of a Japanese Megabat, the Ryukyu Flying Fox (Pteropus dasymallus)

Complete sequences of mitochondrial DNA (mtDNA) are useful for the reconstruction of phylogenetic trees of mammals and, in particular, for inferring higher-order relationships in mammals. In this study, we determined the complete sequence (16,705 bp) of the mtDNA of a Japanese megabat, the Ryukyu flying fox (Pteropus dasymallus). We analyzed this sequence phylogenetically by comparing it with the complete sequence of mtDNAs of 35 mammals in an effort to reevaluate the enigmatic relationship between Megachiroptera and Microchiroptera and the relationships between them and other mammals. Maximum-likelihood analysis of 12 concatenated mitochondrial proteins from 36 mammals strongly suggested the monophyly of the order Chiroptera and its close relationship to Fereuungulata (Carnivora + Perissodactyla + Cetartiodactyla). We estimated that megabats and microbats diverged approximately 58 MyrBP and discussed the origin and early evolution of Chiroptera based on our findings. Received: 28 January 2000 / Accepted: 30 June 2000     

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     

Analysis of mitochondrial DNA restriction fragment length polymorphisms for examining genetic variability among isolates of Phellinus linteus

 Restriction fragment length polymorphism (RFLP) analysis of mitochondrial DNAs (mtDNAs) from nine Japanese wild isolates of Phellinus linteus was carried out to examine their genetic variability. BamHI and EcoRI digests of mtDNAs from these isolates produced four and five distinct RFLP patterns, respectively. By combining the RFLP patterns obtained with the two endonucleases, mtDNAs from the nine isolates could be assigned to five different genotypes, but no mtDNA variation was detected among the isolates collected from a small area. Distance values calculated among all pairs of mtDNA genotypes, based on the presence or absence of comigrating restriction fragments, were clearly smaller than those among the mtDNA genotypes of Lentinula edodes and Pleurotus ostreatus samples collected worldwide, suggesting the necessity of collecting P. linteus wild isolates for genetic resources from geographically wider areas. Received: June 27, 2002 / Accepted: August 19, 2002 Correspondence to:T. Nakamura     

Recombination of mitochondrial DNAs following transmission of mitochondria among incompatible strains of black Aspergilli

Successful intra- and interspecific mitochondrial transfers were performed by polyethylene glycol (PEG)-induced protoplast fusion among incompatible strains belonging to the Aspergillus niger species aggregate. The mitochondrial DNAs (mtDNAs) of the strains examined were of three main types based on their restriction fragment length polymorphism (RFLP) profiles. mtDNA types 1 and 2 correspond to A. niger and A. tubingensis species, respectively, while type 3 is represented by some Brazilian wild-type isolates (possibly a distinct species or subspecies). mtDNA types 1 and 2 could be further divided into several subgroups (1a–1e and 2a–2f ). All these strains, representing different RFLP groups or subgroups, were fully incompatible with respect to nuclear complementation. The transfer experiments were carried out under selection pressure, using a mitochondrial oligomycin-resistant mutant of mtDNA type 1a as donor. Following fusion mitochondrial oligomycin-resistant progenies were recovered in the presence of oligomycin by selecting for the nuclear phenotypes of the oligomycin-sensitive recipient strains. All attempted transfers were successful, and resulted in different varieties of resistant recombinant mitochondrial progenies at various frequencies. Within the group of strains of mtDNA type 1, the transfer of oligomycin-resistant mitochondria resulted in the appearance of a single recombinant type of RFLP profile in each case. The recombination events were more complex when the transfer of oligomycin resistance occurred between strains representing different species (mtDNA groups 1a→2 and 1a→3). A great variety of recombinant mtDNA RFLP profiles appeared. Explanation for this phenomenon are discussed on the basis of preliminary physical mapping data. Received: 16 July 1996 / Accepted: 2 December 1996     

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     

A Complex Organization of the Gene Encoding Cytochrome Oxidase Subunit 1 in the Mitochondrial Genome of the Dinoflagellate, Crypthecodinium cohnii: Homologous Recombination Generates Two Different cox1 Open Reading Frames

In the course of investigating mitochondrial genome organization in Crypthecodinium cohnii, a non-photosynthetic dinoflagellate, we identified four EcoRI fragments that hybridize to a probe specific for cox1, the gene that encodes subunit 1 of cytochrome oxidase. Cloning and sequence characterization of the four fragments (5.7, 5.1, 4.1, 3.5 kilobase pairs) revealed that cox1 exists in four distinct but related contexts in C. cohnii mtDNA, with a central repeat unit flanked by one of two possible upstream (flanking domain 1 or 2) and downstream (flanking domain 3 or 4) regions. The majority of the cox1 gene is located within the central repeat; however, the C-terminal portion of the open reading frame extends into flanking domains 3 and 4, thereby creating two distinct cox1 coding sequences. The 3′-terminal region of one of the cox1 reading frames can assume an elaborate secondary structure, which potentially could act to stabilize the mature mRNA against nucleolytic degradation. In addition, a high density of small inverted repeats (15–22 base pairs) has been identified at the 5′-end of cox1, further suggesting that hairpin structures could be important for gene regulation. The organization of cox1 in C. cohnii mtDNA appears to reflect homologous recombination events within the central repeat between different cox1 sequence contexts. Such recombining repeats are a characteristic feature of plant (angiosperm) mtDNA, but they have not previously been described in the mitochondrial genomes of protists. Received: 21 December 2000 / Accepted: 30 January 2001     

Polymorphisms of Mitochondrial ATPASE 8/6 Genes and Association with Milk Production Traits in Holstein Cows

The maternal effect has been widely proposed to affect the production traits in domestic animals. However, the sequence polymorphisms of mitochondrial DNA (mtDNA) and association with milk production traits in Holstein cows have remained unclear. In this study, we investigated the single nucleotide polymorphisms (SNPs) of mtDNA ATPase 8/6 genes and association with four milk production traits of interest in 303 Holstein cows. A total of 18 SNPs were detected among the 842 bp fragment of ATPase 8/6 genes, which determined six haplotypes of B. taurus (H1-H4) and B. indicus (H5-H6). The mixed model analysis revealed that there was significant association between haplotype and 305-day milk yield (MY). The highest MY was observed in haplotype H4. However, we did not detect statistically significant differences among haplotypes for the traits of milk fat (MF), milk protein (MP), and somatic cell count (SC). The overall haplotype diversity and nucleotide diversity of ATPase 8/6 genes were 0.563 ± 0.030 and 0.00609 ± 0.00043, respectively. The results suggested that mitochondrial ATPase 8/6 genes could be potentially used as molecular marker to genetically improve milk production in Holstein cows.     

THE EVOLUTIONARY HISTORY OF PARTHENOGENETIC CNEMIDOPHORUS LEMNISCATUS (SAURIA: TEIIDAE). II. MATERNAL ORIGIN AND AGE INFERRED FROM MITOCHONDRIAL DNA ANALYSES

Restriction endonuclease analyses were performed on mitochondrial DNAs (mtDNAs) representing unisexual parthenogenetic (cytotypes A, B, and C) and bisexual (cytotypes D and E) populations of Amazonian lizards presently regarded as Cnemidophorus lemniscatus. The results of mtDNA cleavage map comparisons among these C. lemniscatus indicated that (1) there was no cleavage site variation among the unisexuals, (2) mtDNAs from the bisexual cytotypes D and E differed in sequence from one another by about 13%, and (3) mtDNAs from cytotypes A–C differed from those of cytotype D by about 5% and from those of cytotype E by about 13%. Higher resolution restriction fragment size comparisons confirmed the high degree of similarity among the unisexual mtDNAs, but identified 12 cleavage site variants among the 13 cytotype D mtDNAs examined. Both cladistic and phenetic (UPGMA) analyses of the data indicate that the unisexual and cytotype D mtDNAs form a single clade, suggesting that a female of cytotype D was the maternal progenitor of the unisexuals. The similarity among the unisexual mtDNAs and the variability among those of cytotype D suggest that the three unisexual cytotypes arose recently from a common maternal lineage. The mtDNA variability observed among cytotype D individuals has a strong geographic component, suggesting that the unisexuals arose from one or a few geographically proximal populations. The mtDNA comparisons also support the conclusion, based on allozyme comparisons (Sites et al., 1990, this issue), that cytotypes D and E, although presently allocated to C. lemniscatus, are separate species.     

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     

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     

Mitochondrial genomes of the green macroalga Ulva pertusa (Ulvophyceae,Chlorophyta): novel insights into the evolution of mitogenomes in the Ulvophyceae

To further understand the trends in the evolution of mitochondrial genomes (mitogenomes or mtDNAs) in the Ulvophyceae, the mitogenomes of two separate thalli of Ulva pertusa were sequenced. Two U. pertusa mitogenomes (Up1 and Up2) were 69,333 bp and 64,602 bp in length. These mitogenomes shared two ribosomal RNAs (rRNAs), 28 transfer RNAs (tRNAs), 29 protein‐coding genes, and 12 open reading frames. The 4.7 kb difference in size was attributed to variation in intron content and tandem repeat regions. A total of six introns were present in the smaller U. pertusa mtDNA (Up2), while the larger mtDNA (Up1) had eight. The larger mtDNA had two additional group II introns in two genes (cox1 and cox2) and tandem duplication mutations in noncoding regions. Our results showed the first case of intraspecific variation in chlorophytan mitogenomes and provided further genomic data for the undersampled Ulvophyceae.     

Evolution of the mitochondrial DNA control region in the mbuna (Cichlidae) species flock of lake Malawi, East Africa

Considerable controversy has surrounded the application of mitochondrial DNA data to reconstruction of evolutionary relationships among the endemic cichlids of Lake Malawi. Central to this debate has been the issue of whether lineage sorting is complete, and thus whether these data actually reflect species phylogeny, or simply gene genealogy. Review of all mtDNA control region sequences available for members of one monophyletic subset of this species flock, the Malawi rockfishes, or mbuna, strongly indicates that lineage sorting is incomplete: Character-based analyses of these sequences reconstruct gene, not species, interrelationships. Analysis of the pattern of nucleotide substitutions differentiating these mtDNA alleles suggests that pyrimidine residues undergo transition substitutions more often than do purines. Estimation of the magnitude of derived sequence differentiation in light of the reconstructed gene genealogy suggests that the mbuna may be of considerably more recent vintage than previous molecular characterizations have indicated. Received: 6 April 1996 / Accepted: 3 March 1997     

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     

Variable Substitution Rates of the 18 Domain Sequences in Artemia Hemoglobin

The Artemia hemoglobin is a dimer comprising two nine-domain covalent polymers in quaternary association. Each polymer is encoded by a gene representing nine successive globin domains which have different sequences and are presumed to have been copied originally from a single-domain gene. Two different polymers exist as the result of a complete duplication of the nine-domain gene, allowing the formation of either homodimers or the heterodimer. The total population size of 18 domains comprising nine corresponding pairs, coupled with the probability that they reflect several hundred million years of evolution in the same lineage, provides a unique model in which the process of gene multiplication can be analyzed. The outcome has important implications for the reliability of local molecular clocks. The two polymers differ from each other at 11.7% of amino acid sites; however when corresponding individual domains are compared between polymers, amino acid substitution fluctuates by a factor of 2.7-fold from lowest to highest. This variation is not obvious at the DNA level: Domain pair identity values fluctuate by 1.3-fold. Identity values are, however, uncorrected for multiple substitutions, and both silent and nonsilent changes are pooled. Therefore, to determine the variability in relative substitution rates at the DNA level, we have used the method of Li (1993, J Mol Evol 36:96–99) to determine estimates of nonsynonymous (K _A ) and synonymous (K _S ) substitutions per site for the nine pairs of domains. As expected, the overall level of silent substitutions (K _S of 56.9%) far exceeded nonsilent substitutions (K _A of 6.7%); however, for corresponding domain pairs, K _A fluctuates by 2.3-fold and K _S by 1.7-fold. The large discrepancies reflected in the expressed protein have accrued within a single lineage and the implication is that divergence dates of different genera based on amino acid sequences, even with well-studied proteins of reasonable size, can be wrong by a factor well in excess of 2. Received: 4 June 1997 / Accepted: 17 December 1997