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.     

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     

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     

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     

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     

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     

Sociality and the Rate of rDNA Sequence Evolution in Wasps (Vespidae) and Honeybees (Apis)

Sequence data of mitochondrial 16S ribosomal DNA (mt-rDNA) and nuclear 28S ribosomal DNA (nuc-rDNA) were compared in two honeybee species (Apis mellifera and Apis dorsata) and a selection of 22 wasp species (Vespidae) with different levels of sociality. The averge substitution rates in mt-rDNA and nuc-rDNA were almost-equal in solitary species. In species with larger nests, however, the difference between the nuclear and the mitochondrial substitution rate significantly increased. The average substitution ratio, ψ (nucleotide substitutions in mt-rDNA/nucleotide substitutions in nuc-rDNA) was 1.48 ± 0.12 (SE) among the solitary Eumeninae, 3.70 ± 0.15 among five primitive social Stenogastrinae species, 3.24 ± 0.20 among five Polistinae species, 5.76 ± 0.33 among nine highly eusocial Vespinae, and 12.7 in the two Apis species. The high egg-laying rate and the effective population size skew between the sexes may contribute to the rise of the substitution ratio in the highly eusocial species. Drift and bottleneck effects in the mitochondrial DNA pool during speciation events as well as polyandry may further enhance this phenomenon. Received: 12 January 1998 / Accepted: 28 April 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.     

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     

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     

A New Method for Locating Changes in a Tree Reveals Distinct Nucleotide Polymorphism vs Divergence Patterns in Mouse Mitochondrial Control Region

A new, model-based method was devised to locate nucleotide changes in a given phylogenetic tree. For each site, the posterior probability of any possible change in each branch of the tree is computed. This probabilistic method is a valuable alternative to the maximum parsimony method when base composition is skewed (i.e., different from 25% A, 25% C, 25% G, 25% T): computer simulations showed that parsimony misses more rare → common than common → rare changes, resulting in biased inferred change matrices, whereas the new method appeared unbiased. The probabilistic method was applied to the analysis of the mutation and substitution processes in the mitochondrial control region of mouse. Distinct change patterns were found at the polymorphism (within species) and divergence (between species) levels, rejecting the hypothesis of a neutral evolution of base composition in mitochondrial DNA. Received: 15 March 1999 / Accepted: 7 October 1999     

Estimating the Transition/Transversion Ratio from Independent Pairwise Comparisons with an Assumed Phylogeny

A method is presented for estimating the transition/transversion ratio (TI/TV), based on phylogenetically independent comparisons. TI/TV is a parameter of some models used in phylogeny estimation intended to reflect the fact that nucleotide substitutions are not all equally likely. Previous attempts to estimate TI/TV have commonly faced three problems: (1) few taxa; (2) nonindependence among pairwise comparisons; and (3) multiple hits make the apparent TI/TV between two sequences decrease over time since their divergence, giving a misleading impression of relative substitution probabilities. We have made use of the time dependency, modeling how the observed TI/TV changes over time and extrapolating to estimate the ``instantaneous' TI/TV—the relevant parameter for phylogenetic inference. To illustrate our method, TI/TV was estimated for two mammalian mitochondrial genes. For 26 pairs of cytochrome b sequences, the estimate of TI/TV was 5.5; 16 pairs of 12s rRNA yielded an estimate of 9.5. These estimates are higher than those given by the maximum likelihood method and than those obtained by averaging all possible pairwise comparisons (with or without a two-parameter correction for multiple substitutions). We discuss strengths, weaknesses, and further uses of our method. Received: 22 August 1995 / Accepted: 26 July 1996     

Lineage Specificity of the Evolutionary Dynamics of the mtDNA D-Loop Region in Rodents

This paper reports an intraorder study on the D-loop-containing region of the mitochondrial DNA in rodents. A complete multialignment of this region is not feasible with the exception of some conserved regions. The comparative analysis of 25 complete rodent sequences from 23 species plus one lagomorph has revealed that only the central domain (CD), a conserved region of about 80 bp in the extended termination-associated sequences (ETAS) domain, adjacent to the CD, the ETAS1, and conserved sequence block (CSB) 1 blocks are present in all rodent species, whereas the presence of CSB2 and CSB3 is erratic within the order. We have also found a conserved region of 90 bp located between tRNAPro and ETAS1 present in fat dormouse, squirrel, guinea pig, and rabbit. Repeated sequences are present in both the ETAS and the CSB domain, but the repeats differ in length, copy number, and base composition in different species. The potential use of the D-loop for evolutionary studies has been investigated; the presence/absence of conserved blocks and/or repeated sequences cannot be used as a reliable phylogenetic marker, since in some cases they may be shared by distantly related organisms but not by close ones, while in other ones a relationship between tree topology and presence/absence of such motifs is observed. Better results can be obtained by the use of the CD, which, however, due to its reduced size, when used for tracing a phylogenetic tree, shows some nodes with low statistical support. Received: 26 February 2001 / Accepted: 6 June 2001     

鲹科鱼类线粒体DNA控制区结构及系统发育关系

采用PCR技术获得了9种鲹科鱼类的线粒体DNA控制区全序列,并结合从GenBank中下载的3种鲹科鱼类的相应序列采用ClustalW排序后,对控制区结构进行分析,识别了其终止序列区、中央保守区和保守序列区3个区域,指出了终止相关序列的主体是TACAT与其反向互补序列ATGTA以及一系列保守序列(CSB-F、CSB-E、CSB-D和CSB-1、CSB-2、CSB-3),并给出了它们的一般形式,同时在康氏似鲹控制区的5′和3′两端发现重复序列。以尖吻鲈作为外类群,应用邻接法构建的分子系统树表明:鲹科鱼类分为鲹亚科、鰤亚科、鲳鲹亚科和鰆鲹亚科4个亚科,各自形成单系群;鲹亚科与鰤亚科形成姐妹群,鲳鲹亚科再与他们聚在一起,鰆鲹亚科处于鲹科鱼类的基部,与前面3个亚科聚在一起。     

Determination of the evolutionary relationships in Rattus sensu lato (Rodentia : Muridae) using L1 (LINE-1) amplification events

We determined ∼215 bp of DNA sequence from the 3′-untranslated region (UTR) of 240 cloned L1 (LINE-1) elements isolated from 22 species of Rattus sensu lato and Rattus sensu stricto murine rodents. The sequences were sorted into different L1 subfamilies, and oligonucleotides cognate to them were hybridized to genomic DNA of various taxa. From the distribution of the L1 subfamilies in the various species, we inferred the partial phylogeny of Rattus sensu lato. The four Maxomys species comprise a well-defined clade separate from a monophyletic cluster that contains the two Leopoldamys and four Niviventer species. The Niviventer/Leopoldamys clade, in turn, shares a node with the clade that contains Berylmys, Sundamys, Bandicota, and Rattus sensu stricto. The evolutionary relationships that we deduced agree with and significantly extend the phylogeny of Rattus sensu lato established by other molecular criteria. Furthermore, the L1 amplification events scored here produced a unique phylogenetic tree, that is, in no case did a character (a given L1 amplification event) appear on more than one branch. The lack of homoplasy found in this study supports the robustness of L1 amplification events as phylogenetic markers for the study of mammalian evolution. Received: 8 November 1996 / Accepted: 11 April 1997