Theodore G. Schurr

The structure of human mitochondrial DNA variation

SummaryRestriction analysis of mitochondrial DNA (mtDNA) of 3065 humans from 62 geographic samples identified 149 haplotypes and 81 polymorphic sites. These data were used to test several aspects of the evolutionary past of the human species. A dendrogram depicting the genetic relatedness of all haplotypes shows that the native African populations have the greatest diversity and, consistent with evidence from a variety of sources, suggests an African origin for our species. The data also indicate that two individuals drawn, at random from the entire sample will differ at approximately 0.4% of their mtDNA nucleotide sites, which is somewhat higher than previous estimates. Human mtDNA also exhibits more interpopulation heterogeneity (GST=0.351±0.025) than does nuclear DNA (GST=0.12). Moreover, the virtual absence of intermediate levels of linkage disequilibrium between pairs of sites is consistent with the absence of genetic recombination and places constraints on the rate of mutation. Tests of the selective neutrality of mtDNA variation, including the Ewens-Watterson and Tajima tests, indicate a departure in the direction consistent with purifying selection, but this departure is more likely due to the rapid growth of the human population and the geographic heterogeneity of the variation. The lack of a good fit to neutrality poses problems for the estimation of times of coalescence from human mtDNA data.

mtDNA Variation in the South African Kung and Khwe—and Their Genetic Relationships to Other African Populations

The mtDNA variation of 74 Khoisan-speaking individuals (Kung and Khwe) from Schmidtsdrift, in the Northern Cape Province of South Africa, was examined by high-resolution RFLP analysis and control region (CR) sequencing. The resulting data were combined with published RFLP haplotype and CR sequence data from sub-Saharan African populations and then were subjected to phylogenetic analysis to deduce the evolutionary relationships among them. More than 77% of the Kung and Khwe mtDNA samples were found to belong to the major mtDNA lineage, macrohaplogroup L* (defined by a HpaI site at nucleotide position 3592), which is prevalent in sub-Saharan African populations. Additional sets of RFLPs subdivided macrohaplogroup L* into two extended haplogroups-L1 and L2-both of which appeared in the Kung and Khwe. Besides revealing the significant substructure of macrohaplogroup L* in African populations, these data showed that the Biaka Pygmies have one of the most ancient RFLP sublineages observed in African mtDNA and, thus, that they could represent one of the oldest human populations. In addition, the Kung exhibited a set of related haplotypes that were positioned closest to the root of the human mtDNA phylogeny, suggesting that they, too, represent one of the most ancient African populations. Comparison of Kung and Khwe CR sequences with those from other African populations confirmed the genetic association of the Kung with other Khoisan-speaking peoples, whereas the Khwe were more closely linked to non-Khoisan-speaking (Bantu) populations. Finally, the overall sequence divergence of 214 African RFLP haplotypes defined in both this and an earlier study was 0.364%, giving an estimated age, for all African mtDNAs, of 125,500-165,500 years before the present, a date that is concordant with all previous estimates derived from mtDNA and other genetic data, for the time of origin of modern humans in Africa.

The Dual Origin and Siberian Affinities of Native American Y Chromosomes

The Y chromosomes of 549 individuals from Siberia and the Americas were analyzed for 12 biallelic markers, which defined 15 haplogroups. The addition of four microsatellite markers increased the number of haplotypes to 111. The major Native American founding lineage, haplogroup M3, accounted for 66% of male Y chromosomes and was defined by the biallelic markers M89, M9, M45, and M3. The founder haplotype also harbored the microsatellite alleles DYS19 (10 repeats), DYS388 (11 repeats), DYS390 (11 repeats), and DYS391 (10 repeats). In Siberia, the M3 haplogroup was confined to the Chukotka peninsula, adjacent to Alaska. The second major group of Native American Y chromosomes, haplogroup M45, accounted for about one-quarter of male lineages. M45 was subdivided by the biallelic marker M173 and by the four microsatellite loci alleles into two major subdivisions: M45a, which is found throughout the Americas, and M45b, which incorporates the M173 variant and is concentrated in North and Central America. In Siberia, M45a haplotypes, including the direct ancestor of haplogroup M3, are concentrated in Middle Siberia, whereas M45b haplotypes are found in the Lower Amur River and Sea of Okhotsk regions of eastern Siberia. Among the remaining 5% of Native American Y chromosomes is haplogroup RPS4Y-T, found in North America. In Siberia, this haplogroup, along with haplogroup M45b, is concentrated in the Lower Amur River/Sea of Okhotsk region. These data suggest that Native American male lineages were derived from two major Siberian migrations. The first migration originated in southern Middle Siberia with the founding haplotype M45a (10-11-11-10). In Beringia, this gave rise to the predominant Native American lineage, M3 (10-11-11-10), which crossed into the New World. A later migration came from the Lower Amur/Sea of Okhkotsk region, bringing haplogroup RPS4Y-T and subhaplogroup M45b, with its associated M173 variant. This migration event contributed to the modern genetic pool of the Na-Dene and Amerinds of North and Central America.

Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation

Understanding the distribution of human genetic variation is an important foundation for research into the genetics of common diseases. Some of the alleles that modify common disease risk are themselves likely to be common and, thus, amenable to identification using gene-association methods. A problem with this approach is that the large sample sizes required for sufficient statistical power to detect alleles with moderate effect make gene-association studies susceptible to false-positive findings as the result of population stratification [1, 2]. Such type I errors can be eliminated by using either family-based association tests or methods that sufficiently adjust for population stratification [3–5]. These methods require the availability of genetic markers that can detect and, thus, control for sources of genetic stratification among populations. In an effort to investigate population stratification and identify appropriate marker panels, we have analysed 11,555 single nucleotide polymorphisms in 203 individuals from 12 diverse human populations. Individuals in each population cluster to the exclusion of individuals from other populations using two clustering methods. Higher-order branching and clustering of the populations are consistent with the geographic origins of populations and with previously published genetic analyses. These data provide a valuable resource for the definition of marker panels to detect and control for population stratification in population-based gene identification studies. Using three US resident populations (European-American, African-American and Puerto Rican), we demonstrate how such studies can proceed, quantifying proportional ancestry levels and detecting significant admixture structure in each of these populations.

Melanesian mtDNA Complexity

Melanesian populations are known for their diversity, but it has been hard to grasp the pattern of the variation or its underlying dynamic. Using 1,223 mitochondrial DNA (mtDNA) sequences from hypervariable regions 1 and 2 (HVR1 and HVR2) from 32 populations, we found the among-group variation is structured by island, island size, and also by language affiliation. The more isolated inland Papuan-speaking groups on the largest islands have the greatest distinctions, while shore dwelling populations are considerably less diverse (at the same time, within-group haplotype diversity is less in the most isolated groups). Persistent differences between shore and inland groups in effective population sizes and marital migration rates probably cause these differences. We also add 16 whole sequences to the Melanesian mtDNA phylogenies. We identify the likely origins of a number of the haplogroups and ancient branches in specific islands, point to some ancient mtDNA connections between Near Oceania and Australia, and show additional Holocene connections between Island Southeast Asia/Taiwan and Island Melanesia with branches of haplogroup E. Coalescence estimates based on synonymous transitions in the coding region suggest an initial settlement and expansion in the region at ∼30–50,000 years before present (YBP), and a second important expansion from Island Southeast Asia/Taiwan during the interval ∼3,500–8,000 YBP. However, there are some important variance components in molecular dating that have been overlooked, and the specific nature of ancestral (maternal) Austronesian influence in this region remains unresolved.

Parallel Evolution of Genes and Languages in the Caucasus Region

We analyzed 40 single nucleotide polymorphism and 19 short tandem repeat Y-chromosomal markers in a large sample of 1,525 indigenous individuals from 14 populations in the Caucasus and 254 additional individuals representing potential source populations. We also employed a lexicostatistical approach to reconstruct the history of the languages of the North Caucasian family spoken by the Caucasus populations. We found a different major haplogroup to be prevalent in each of four sets of populations that occupy distinct geographic regions and belong to different linguistic branches. The haplogroup frequencies correlated with geography and, even more strongly, with language. Within haplogroups, a number of haplotype clusters were shown to be specific to individual populations and languages. The data suggested a direct origin of Caucasus male lineages from the Near East, followed by high levels of isolation, differentiation, and genetic drift in situ. Comparison of genetic and linguistic reconstructions covering the last few millennia showed striking correspondences between the topology and dates of the respective gene and language trees and with documented historical events. Overall, in the Caucasus region, unmatched levels of gene-language coevolution occurred within geographically isolated populations, probably due to its mountainous terrain.

Mitochondrial DNA and Y Chromosome Variation Provides Evidence for a Recent Common Ancestry between Native Americans and Indigenous Altaians

The Altai region of southern Siberia has played a critical role in the peopling of northern Asia as an entry point into Siberia and a possible homeland for ancestral Native Americans. It has an old and rich history because humans have inhabited this area since the Paleolithic. Today, the Altai region is home to numerous Turkic-speaking ethnic groups, which have been divided into northern and southern clusters based on linguistic, cultural, and anthropological traits. To untangle Altaian genetic histories, we analyzed mtDNA and Y chromosome variation in northern and southern Altaian populations. All mtDNAs were assayed by PCR-RFLP analysis and control region sequencing, and the nonrecombining portion of the Y chromosome was scored for more than 100 biallelic markers and 17 Y-STRs. Based on these data, we noted differences in the origin and population history of Altaian ethnic groups, with northern Altaians appearing more like Yeniseian, Ugric, and Samoyedic speakers to the north, and southern Altaians having greater affinities to other Turkic speaking populations of southern Siberia and Central Asia. Moreover, high-resolution analysis of Y chromosome haplogroup Q has allowed us to reshape the phylogeny of this branch, making connections between populations of the New World and Old World more apparent and demonstrating that southern Altaians and Native Americans share a recent common ancestor. These results greatly enhance our understanding of the peopling of Siberia and the Americas.

Y chromosome polymorphisms in native American and Siberian populations: identification of native American Y chromosome haplotypes.

Abstract We have initiated a study of ancient male migrations from Siberia to the Americas using Y chromosome polymorphisms. The first polymorphism examined, a C→T transition at nucleotide position 181 of the DYS199 locus, was previously reported only in Native American populations. To investigate the origin of this DYS199 polymorphism, we screened Y chromosomes from a number of Siberian, Asian, and Native American populations for this and other markers. This survey detected the T allele in all five Native American populations studied at an average frequency of 61%, and in two of nine native Siberian populations, the Siberian Eskimo (21%) and the Chukchi (17%). This finding suggested that the DYS199 T allele may have originated in Beringia and was then spread throughout the New World by the founding populations of the major subgroups of modern Native Americans. We further characterized Native American Y chromosome variation by analyzing two additional Y chromosome polymorphisms, the DYS287 Y Alu polymorphic (YAP) element insertion and a YAP-associated A→G transition at DYS271, both commonly found in Africans. We found neither African allele associated with the DYS199 T allele in any of the Native American or native Siberian populations. However, we did find DYS287 YAP+ individuals who harbored the DYS199 C allele in one Native American population, the Mixe, and in one Asian group, the Tibetans. A correlation of these Y chromosome alleles in Native Americans with those of the DYS1 locus, as detected by the p49a/p49f (p49a,f) probes on TaqI-digested genomic DNA, revealed a complete association of DYS1 alleles (p49a,f haplotypes) 13, 18, 66, 67 and 69 with the DYS199 T allele, while DYS1 alleles 8 and 63 were associated with both the DYS199 C and T allele.

A private allele ubiquitous in the Americas

The three-wave migration hypothesis of Greenberg et al. has permeated the genetic literature on the peopling of the Americas. Greenberg et al. proposed that Na-Dene, Aleut-Eskimo and Amerind are language phyla which represent separate migrations from Asia to the Americas. We show that a unique allele at autosomal microsatellite locus D9S1120 is present in all sampled North and South American populations, including the Na-Dene and Aleut-Eskimo, and in related Western Beringian groups, at an average frequency of 31.7%. This allele was not observed in any sampled putative Asian source populations or in other worldwide populations. Neither selection nor admixture explains the distribution of this regionally specific marker. The simplest explanation for the ubiquity of this allele across the Americas is that the same founding population contributed a large fraction of ancestry to all modern Native American populations.

Geographic population structure analysis of worldwide human populations infers their biogeographical origins

The search for a method that utilizes biological information to predict humans’ place of origin has occupied scientists for millennia. Over the past four decades, scientists have employed genetic data in an effort to achieve this goal but with limited success. While biogeographical algorithms using next-generation sequencing data have achieved an accuracy of 700 km in Europe, they were inaccurate elsewhere. Here we describe the Geographic Population Structure (GPS) algorithm and demonstrate its accuracy with three data sets using 40,000–130,000 SNPs. GPS placed 83% of worldwide individuals in their country of origin. Applied to over 200 Sardinians villagers, GPS placed a quarter of them in their villages and most of the rest within 50 km of their villages. GPS’s accuracy and power to infer the biogeography of worldwide individuals down to their country or, in some cases, village, of origin, underscores the promise of admixture-based methods for biogeography and has ramifications for genetic ancestry testing.