Major population expansion of East Asians began before neolithic time: evidence of mtDNA genomes

It is a major question in archaeology and anthropology whether human populations started to grow primarily after the advent of agriculture, i.e., the Neolithic time, especially in East Asia, which was one of the centers of ancient agricultural civilization. To answer this question requires an accurate estimation of the time of lineage expansion as well as that of population expansion in a population sample without ascertainment bias. In this study, we analyzed all available mtDNA genomes of East Asians ascertained by random sampling, a total of 367 complete mtDNA sequences generated by the 1000 Genome Project, including 249 Chinese (CHB, CHD, and CHS) and 118 Japanese (JPT). We found that major mtDNA lineages underwent expansions, all of which, except for two JPT-specific lineages, including D4, D4b2b, D4a, D4j, D5a2a, A, N9a, F1a1'4, F2, B4, B4a, G2a1 and M7b1'2'4, occurred before 10 kya, i.e., before the Neolithic time (symbolized by Dadiwan Culture at 7.9 kya) in East Asia. Consistent to this observation, the further analysis showed that the population expansion in East Asia started at 13 kya and lasted until 4 kya. The results suggest that the population growth in East Asia constituted a need for the introduction of agriculture and might be one of the driving forces that led to the further development of agriculture.     

mtDNA variation predicts population size in humans and reveals a major Southern Asian chapter in human prehistory

The relative timing and size of regional human population growth following our expansion from Africa remain unknown. Human mitochondrial DNA (mtDNA) diversity carries a legacy of our population history. Given a set of sequences, we can use coalescent theory to estimate past population size through time and draw inferences about human population history. However, recent work has challenged the validity of using mtDNA diversity to infer species population sizes. Here we use Bayesian coalescent inference methods, together with a global data set of 357 human mtDNA coding-region sequences, to infer human population sizes through time across 8 major geographic regions. Our estimates of relative population sizes show remarkable concordance with the contemporary regional distribution of humans across Africa, Eurasia, and the Americas, indicating that mtDNA diversity is a good predictor of population size in humans. Plots of population size through time show slow growth in sub-Saharan Africa beginning 143-193 kya, followed by a rapid expansion into Eurasia after the emergence of the first non-African mtDNA lineages 50-70 kya. Outside Africa, the earliest and fastest growth is inferred in Southern Asia approximately 52 kya, followed by a succession of growth phases in Northern and Central Asia (approximately 49 kya), Australia (approximately 48 kya), Europe (approximately 42 kya), the Middle East and North Africa (approximately 40 kya), New Guinea (approximately 39 kya), the Americas (approximately 18 kya), and a second expansion in Europe (approximately 10-15 kya). Comparisons of relative regional population sizes through time suggest that between approximately 45 and 20 kya most of humanity lived in Southern Asia. These findings not only support the use of mtDNA data for estimating human population size but also provide a unique picture of human prehistory and demonstrate the importance of Southern Asia to our recent evolutionary past.     

Genetic evidence for unequal effective population sizes of human females and males

The time to the most recent common ancestor (TMRCA) of the human mitochondria (mtDNA) is estimated to be older than that of the nonrecombining portion of the Y chromosome (NRY). Surveys of variation in globally distributed humans typically result in mtDNA TMRCA values just under 200 thousand years ago (kya), whereas those for the NRY range between 46 and 110 kya. A favored hypothesis for this finding is that natural selection has acted on the NRY, leading to a recent selective sweep. An alternate hypothesis is that sex-biased demographic processes are responsible. Here, we re-examine the disparity between NRY and mtDNA TMRCAs using data collected from individual human populations--a sampling strategy that minimizes the confounding influence of population subdivision in global data sets. We survey variation at 782 bp of the mitochondrial cytochrome c oxidase subunit 3 gene as well as at 26.5 kb of noncoding DNA from the NRY in a sample of 25 Khoisan, 24 Mongolians, and 24 Papua New Guineans. Data from both loci in all populations are best described by a model of constant population size, with the exception of Mongolian mtDNA, which appears to be experiencing rapid population growth. Taking these demographic models into account, we estimate the TMRCAs for each locus in each population. A pattern that is remarkably consistent across all three populations is an approximately twofold deeper coalescence for mtDNA than for the NRY. The oldest TMRCAs are observed for the Khoisan (73.6 kya for the NRY and 176.5 kya for mtDNA), whereas those in the non-African populations are consistently lower (averaging 47.7 kya for the NRY and 92.8 kya for mtDNA). Our data do not suggest that differential natural selection is the cause of this difference in TMRCAs. Rather, these results are most consistent with a higher female effective population size.     

Analysis of mitochondrial DNA diversity in the aleuts of the commander islands and its implications for the genetic history of beringia

The Aleuts are aboriginal inhabitants of the Aleutian archipelago, including Bering and Copper (Medny) Islands of the Commanders, and seem to be the survivors of the inhabitants of the southern belt of the Bering Land Bridge that connected Chukotka/Kamchatka and Alaska during the end of the Ice Age. Thirty mtDNA samples collected in the Commanders, as well as seven mtDNA samples from Sireniki Eskimos in Chukotka who belong to the Beringian-specific subhaplogroup D2, were studied through complete sequencing. This analysis has provided evidence that all 37 of these mtDNAs are closely related, since they share the founding haplotype for subhaplogroup D2. We also demonstrated that, unlike the Eskimos and Na-Dene, the Aleuts of the Commanders were founded by a single lineage of haplogroup D2, which had acquired the novel transversion mutation 8910A. The phylogeny of haplogroup D complete sequences showed that (1) the D2 root sequence type originated among the latest inhabitants of Beringia and (2) the Aleut 8910A sublineage of D2 is a part of larger radiation of rooted D2, which gave rise to D2a (Na-Dene), D2b (Aleut), and D2c (Eskimo) sublineages. The geographic specificity and remarkable intrinsic diversity of D2 lineages support the refugial hypothesis, which assumes that the founding population of Eskimo-Aleut originated in Beringan/southwestern Alaskan refugia during the early postglacial period, rather than having reached the shores of Alaska as the result of recent wave of migration from interior Siberia.     

Complete mitochondrial DNA replacement in a Lake Tanganyika cichlid fish

We used nuclear and mitochondrial DNA (mtDNA) sequences from specimens collected throughout Lake Tanganyika to clarify the evolutionary relationship between Lamprologus callipterus and Neolamprologus fasciatus . The nuclear data support the reciprocal monophyly of these two shell-breeding lamprologine cichlids. However, mtDNA sequences show that (i) L. callipterus includes two divergent and geographically disjunct (North–South) mtDNA lineages; and that (ii) N. fasciatus individuals cluster in a lineage sister group to the northern lineage of L. callipterus . The two mtDNA lineages of L. callipterus diverged c . 684 kya to 1.2 Ma, coinciding with a major water level low stand in Lake Tanganyika, which divided the lake into isolated sub-lakes. This suggests that the two mtDNA lineages originated as the result of the separation of L. callipterus populations in different sub-basins. The incongruent phylogenetic position of N. fasciatus can best be explained by an ancient unidirectional introgression from L. callipterus into N. fasciatus. Remarkably, our data indicate that this event resulted in the complete mtDNA replacement in N. fasciatus . Our data suggest that hybridization occurred soon after the divergence of the two L. callipterus mtDNA lineages, probably still during the water level low stand, and that subsequently the invading mtDNA lineage spread throughout the lake.     

The Genetic Legacy of the Pre-Colonial Period in Contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1–2%) of a sub-Saharan component in Bolivians.     

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200–300 thousand years ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~ 50–140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~ 23–50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.     

Revised Timeline and Distribution of the Earliest Diverged Human Maternal Lineages in Southern Africa

The oldest extant human maternal lineages include mitochondrial haplogroups L0d and L0k found in the southern African click-speaking forager peoples broadly classified as Khoesan. Profiling these early mitochondrial lineages allows for better understanding of modern human evolution. In this study, we profile 77 new early-diverged complete mitochondrial genomes and sub-classify another 105 L0d/L0k individuals from southern Africa. We use this data to refine basal phylogenetic divergence, coalescence times and Khoesan prehistory. Our results confirm L0d as the earliest diverged lineage (∼172 kya, 95%CI: 149–199 kya), followed by L0k (∼159 kya, 95%CI: 136–183 kya) and a new lineage we name L0g (∼94 kya, 95%CI: 72–116 kya). We identify two new L0d1 subclades we name L0d1d and L0d1c4/L0d1e, and estimate L0d2 and L0d1 divergence at ∼93 kya (95%CI:76–112 kya). We concur the earliest emerging L0d1’2 sublineage L0d1b (∼49 kya, 95%CI:37–58 kya) is widely distributed across southern Africa. Concomitantly, we find the most recent sublineage L0d2a (∼17 kya, 95%CI:10–27 kya) to be equally common. While we agree that lineages L0d1c and L0k1a are restricted to contemporary inland Khoesan populations, our observed predominance of L0d2a and L0d1a in non-Khoesan populations suggests a once independent coastal Khoesan prehistory. The distribution of early-diverged human maternal lineages within contemporary southern Africans suggests a rich history of human existence prior to any archaeological evidence of migration into the region. For the first time, we provide a genetic-based evidence for significant modern human evolution in southern Africa at the time of the Last Glacial Maximum at between ∼21–17 kya, coinciding with the emergence of major lineages L0d1a, L0d2b, L0d2d and L0d2a.     

Rapid deployment of the five founding Amerind mtDNA haplogroups via coastal and riverine colonization

Numerous studies of variation in mtDNA in Amerindian populations established that four haplogroups are present throughout both North and South America. These four haplogroups (A, B, C, and D) and perhaps a fifth (X) in North America are postulated to be present in the initial founding migration to the Americas. Furthermore, studies of ancient mtDNA in North America suggested long-term regional continuity of the frequencies of these founding haplogroups. Present-day tribal groups possess high frequencies of private mtDNA haplotypes (variants within the major haplogroups), consistent with early establishment of local isolation of regional populations. Clearly these patterns have implications for the mode of colonization of the hemisphere. Recently, the earlier consensus among archaeologists for an initial colonization by Clovis hunters arriving through an ice-free corridor and expanding in a "blitzkrieg " wave was shown to be inconsistent with extensive genetic variability in Native Americans; a coastal migration route avoids this problem. The present paper demonstrates through a computer simulation model how colonization along coasts and rivers could have rapidly spread the founding lineages widely through North America.     

mtDNA variation indicates Mongolia may have been the source for the founding population for the New World.

mtDNA RFLP variation was analyzed in 42 Mongolians from Ulan Bator. All four founding lineage types (A [4.76%], B [2.38%], C [11.9%], and D [19.04%]) identified by Torroni and colleagues were detected. Seven of the nine founding lineage types proposed by Bailliet and colleagues and Merriwether and Ferrell were detected (A2 [4.76%], B [2.38%], C1 [11.9%], D1 [7.14%], D2 [11.9%], X6 [16.7%], and X7 [9.5%]). Sixty-four percent of these 42 individuals had "Amerindian founding lineage" haplotypes. A survey of 24 restriction sites yielded 16 polymorphic sites and 21 different haplotypes. The presence of all four of the founding lineages identified by the Torroni group (and seven of Merriwether and Ferrell's nine founding lineages), combined with Mongolia's location with respect to the Bering Strait, indicates that Mongolia is a potential location for the origin of the founders of the New World. Since lineage B, which is widely distributed in the New World, is absent in Siberia, we conclude that Mongolia or a geographic location common to both contemporary Mongolians and American aboriginals is the more likely origin of the founders of the New World.     

Founding Amerindian mitochondrial DNA lineages in ancient Maya from Xcaret, Quintana Roo

Ancient DNA from the bone remains of 25 out of 28 pre-Columbian individuals from the Late Classic-Postclassic Maya site of Xcaret, Quintana Roo, was recovered, and mitochondrial DNA (mtDNA) was amplified by using the polymerase chain reaction. The presence of the four founding Amerindian mtDNA lineages was investigated by restriction analysis and by direct sequencing in selected individuals. The mtDNA lineages A, B, and C were found in this population. Eighty-four percent of the individuals were lineage A, whereas lineages B and C were present at low frequencies, 4% and 8%, respectively. Lineage D was absent from our sample. One individual did not possess any of the four lineages. Six skeletons out of 7 dated from the Late Classic period were haplotype A, whereas 11 skeletons out of 16 dated from the Postclassic period were also haplotype A. The distribution of mtDNA lineages in the Xcaret population contrasts sharply with that found in ancient Maya from Copán, which lack lineages A and B. On the other hand, our results resemble more closely the frequencies of mtDNA lineages found in contemporary Maya from the Yucatán Peninsula and in other Native American contemporary populations of Mesoamerican origin.     

mtDNA haplogroup X: An ancient link between Europe/Western Asia and North America?

On the basis of comprehensive RFLP analysis, it has been inferred that approximately 97% of Native American mtDNAs belong to one of four major founding mtDNA lineages, designated haplogroups "A"-"D." It has been proposed that a fifth mtDNA haplogroup (haplogroup X) represents a minor founding lineage in Native Americans. Unlike haplogroups A-D, haplogroup X is also found at low frequencies in modern European populations. To investigate the origins, diversity, and continental relationships of this haplogroup, we performed mtDNA high-resolution RFLP and complete control region (CR) sequence analysis on 22 putative Native American haplogroup X and 14 putative European haplogroup X mtDNAs. The results identified a consensus haplogroup X motif that characterizes our European and Native American samples. Among Native Americans, haplogroup X appears to be essentially restricted to northern Amerindian groups, including the Ojibwa, the Nuu-Chah-Nulth, the Sioux, and the Yakima, although we also observed this haplogroup in the Na-Dene-speaking Navajo. Median network analysis indicated that European and Native American haplogroup X mtDNAs, although distinct, nevertheless are distantly related to each other. Time estimates for the arrival of X in North America are 12,000-36,000 years ago, depending on the number of assumed founders, thus supporting the conclusion that the peoples harboring haplogroup X were among the original founders of Native American populations. To date, haplogroup X has not been unambiguously identified in Asia, raising the possibility that some Native American founders were of Caucasian ancestry.     

mtDNA and the origin of the Icelanders: deciphering signals of recent population history

Previous attempts to investigate the origin of the Icelanders have provided estimates of ancestry ranging from a 98% British Isles contribution to an 86% Scandinavian contribution. We generated mitochondrial sequence data for 401 Icelandic individuals and compared these data with >2,500 other European sequences from published sources, to determine the probable origins of women who contributed to Iceland’s settlement. Although the mean number of base-pair differences is high in the Icelandic sequences and they are widely distributed in the overall European mtDNA phylogeny, we find a smaller number of distinct mitochondrial lineages, compared with most other European populations. The frequencies of a number of mtDNA lineages in the Icelanders deviate noticeably from those in neighboring populations, suggesting that founder effects and genetic drift may have had a considerable influence on the Icelandic gene pool. This is in accordance with available demographic evidence about Icelandic population history. A comparison with published mtDNA lineages from European populations indicates that, whereas most founding females probably originated from Scandinavia and the British Isles, lesser contributions from other populations may also have taken place. We present a highly resolved phylogenetic network for the Icelandic data, identifying a number of previously unreported mtDNA lineage clusters and providing a detailed depiction of the evolutionary relationships between European mtDNA clusters. Our findings indicate that European populations contain a large number of closely related mitochondrial lineages, many of which have not yet been sampled in the current comparative data set. Consequently, substantial increases in sample sizes that use mtDNA data will be needed to obtain valid estimates of the diverse ancestral mixtures that ultimately gave rise to contemporary populations.     

Mitochondrial population genomics supports a single pre-Clovis origin with a coastal route for the peopling of the Americas

It is well accepted that the Americas were the last continents reached by modern humans, most likely through Beringia. However, the precise time and mode of the colonization of the New World remain hotly disputed issues. Native American populations exhibit almost exclusively five mitochondrial DNA (mtDNA) haplogroups (A-D and X). Haplogroups A-D are also frequent in Asia, suggesting a northeastern Asian origin of these lineages. However, the differential pattern of distribution and frequency of haplogroup X led some to suggest that it may represent an independent migration to the Americas. Here we show, by using 86 complete mitochondrial genomes, that all Native American haplogroups, including haplogroup X, were part of a single founding population, thereby refuting multiple-migration models. A detailed demographic history of the mtDNA sequences estimated with a Bayesian coalescent method indicates a complex model for the peopling of the Americas, in which the initial differentiation from Asian populations ended with a moderate bottleneck in Beringia during the last glacial maximum (LGM), around approximately 23,000 to approximately 19,000 years ago. Toward the end of the LGM, a strong population expansion started approximately 18,000 and finished approximately 15,000 years ago. These results support a pre-Clovis occupation of the New World, suggesting a rapid settlement of the continent along a Pacific coastal route.     

Phylogeographic analysis of mitochondrial DNA in northern Asian populations

To elucidate the human colonization process of northern Asia and human dispersals to the Americas, a diverse subset of 71 mitochondrial DNA (mtDNA) lineages was chosen for complete genome sequencing from the collection of 1,432 control-region sequences sampled from 18 autochthonous populations of northern, central, eastern, and southwestern Asia. On the basis of complete mtDNA sequencing, we have revised the classification of haplogroups A, D2, G1, M7, and I; identified six new subhaplogroups (I4, N1e, G1c, M7d, M7e, and J1b2a); and fully characterized haplogroups N1a and G1b, which were previously described only by the first hypervariable segment (HVS1) sequencing and coding-region restriction-fragment-length polymorphism analysis. Our findings indicate that the southern Siberian mtDNA pool harbors several lineages associated with the Late Upper Paleolithic and/or early Neolithic dispersals from both eastern Asia and southwestern Asia/southern Caucasus. Moreover, the phylogeography of the D2 lineages suggests that southern Siberia is likely to be a geographical source for the last postglacial maximum spread of this subhaplogroup to northern Siberia and that the expansion of the D2b branch occurred in Beringia ~7,000 years ago. In general, a detailed analysis of mtDNA gene pools of northern Asians provides the additional evidence to rule out the existence of a northern Asian route for the initial human colonization of Asia.     

A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.     

mtDNA analysis of a prehistoric Oneota population: implications for the peopling of the New World.

mtDNA was successfully extracted from 108 individuals from the Norris Farms Oneota, a prehistoric Native American population, to compare the mtDNA diversity from a pre-Columbian population with contemporary Native American and Asian mtDNA lineages and to examine hypotheses about the peopling of the New World. Haplogroup and hypervariable region I sequence data indicate that the lineages from haplogroups A, B, C, and D are the most common among Native Americans but that they were not the only lineages brought into the New World from Asia. The mtDNA evidence does not support the three-wave hypothesis of migration into the New World but rather suggests a single "wave" of people with considerable mtDNA diversity that exhibits a signature of expansion 23,000-37,000 years ago.     

Mapping Bombina mitochondrial genomes: the conundrum of Carpathian Bombina variegata (Anura: Discoglossidae)

Complete mitochondrial genomes (mtDNA) of five individuals representing two haplotypes of Bombina bombina and three of Bombina variegata were compared using restriction site maps. Phylogenetic analyses reveal three ancient mitochondrial lineages: (1) two very similar haplotypes A and B of B. bombina ; (2) almost identical haplotypes D and E of B. variegata ; and (3) haplotype C of B. variegata . Haplotype C is as different from haplotypes D/E as from A/B. These data are strikingly discordant with relationships based on morphology and allozymes. Haplotypes C and D/E represent a pre-Pleistocene mitochondrial divergence within B. variegata , nearly coincident with speciation between B. variegata and B. bombina . Geographical partitioning of the two divergent B. variegata mitochondrial lineages indicates repeated localization of the lineages in separate glacial refugia during the Pleistocene. That nuclear genes do not show a similar divergence, but rather indicate relatively free genetic exchange between populations with divergent mtDNAs, suggests that males dispersed much more widely than females during expansions from glacial refugia. Comparison of Bombina mtDNA maps with a restriction site map of Xenopus laevis mtDNA revealed 16 homologous sites; 12 of these may be nearly invariant across primitive anuran mtDNAs. Two distinct regions of heteroplasmy, representing two regions with variable numbers of sequence repeats [length variable (LV) regions], were characterized. Comparison with the Xenopus map places LV1, present in all five haplotypes, near the 5'-end of the control region, and LV2 present only in B. variegata , near the 3'-end. Although phylogenetic analyses did not group the two major B. variegata lineages together, presence of LV2 in both lineages supports placement of both within B. variegata .     

Geographic Variation in Human Mitochondrial DNA from Papua New Guinea

High resolution mitochondrial DNA (mtDNA) restriction maps, consisting of an average of 370 sites per mtDNA map, were constructed for 119 people from 25 localities in Papua New Guinea (PNG). Comparison of these PNG restriction maps to published maps from Australian, Caucasian, Asian and African mtDNAs reveals that PNG has the lowest amount of mtDNA variation, and that PNG mtDNA lineages originated from Southeast Asia. The statistical significance of geographic structuring of populations with respect to mtDNA was assessed by comparing observed GST values to a distribution of GST values generated by random resampling of the data. These analyses show that there is significant structuring of mtDNA variation among worldwide populations, between highland and coastal PNG populations, and even between two highland PNG populations located approximately 200 km apart. However, coastal PNG populations are essentially panmictic, despite being spread over several hundred kilometers. Highland PNG populations also have more mtDNA variability and more mtDNA types represented per founding lineage than coastal PNG populations. All of these observations are consistent with a more ancient, restricted origin of highland PNG populations, internal isolation of highland PNG populations from one another and from coastal populations, and more recent and extensive population movements through coastal PNG. An apparent linguistic effect on PNG mtDNA variation disappeared when geography was taken into account. The high resolution technique for examining mtDNA variation, coupled with extensive geographic sampling within a single defined area, leads to an enhanced understanding of the influence of geography on mtDNA variation in human populations.