Mitochondrial DNA Diversity in Indigenous Populations of the Southern Extent of Siberia, and the Origins of Native American Haplogroups

In search of the ancestors of Native American mitochondrial DNA (mtDNA) haplogroups, we analyzed the mtDNA of 531 individuals from nine indigenous populations in Siberia. All mtDNAs were subjected to high‐resolution RFLP analysis, sequencing of the control‐region hypervariable segment I (HVS‐I), and surveyed for additional polymorphic markers in the coding region. Furthermore, the mtDNAs selected according to haplogroup/subhaplogroup status were completely sequenced. Phylogenetic analyses of the resulting data, combined with those from previously published Siberian arctic and sub‐arctic populations, revealed that remnants of the ancient Siberian gene pool are still evident in Siberian populations, suggesting that the founding haplotypes of the Native American A‐D branches originated in different parts of Siberia. Thus, lineage A complete sequences revealed in the Mansi of the Lower Ob and the Ket of the Lower Yenisei belong to A1, suggesting that A1 mtDNAs occasionally found in the remnants of hunting‐gathering populations of northwestern and northern Siberia belonged to a common gene pool of the Siberian progenitors of Paleoindians. Moreover, lineage B1, which is the most closely related to the American B2, occurred in the Tubalar and Tuvan inhabiting the territory between the upper reaches of the Ob River in the west, to the Upper Yenisei region in the east. Finally, the sequence variants of haplogroups C and D, which are most similar to Native American C1 and D1, were detected in the Ulchi of the Lower Amur. Overall, our data suggest that the immediate ancestors of the Siberian/Beringian migrants who gave rise to ancient (pre‐Clovis) Paleoindians have a common origin with aboriginal people of the area now designated the Altai‐Sayan Upland, as well as the Lower Amur/Sea of Okhotsk region.     

Diversity of Mitochondrial DNA Lineages in South Siberia

To investigate the origin and evolution of aboriginal populations of South Siberia, a comprehensive mitochondrial DNA (mtDNA) analysis (HVR1 sequencing combined with RFLP typing) of 480 individuals, representing seven Altaic‐speaking populations (Altaians, Khakassians, Buryats, Sojots, Tuvinians, Todjins and Tofalars), was performed. Additionally, HVR2 sequence information was obtained for 110 Altaians, providing, in particular, some novel details of the East Asian mtDNA phylogeny. The total sample revealed 81% East Asian (M*, M7, M8, M9, M10, C, D, G, Z, A, B, F, N9a, Y) and 17% West Eurasian (H, U, J, T, I, N1a, X) matrilineal genetic contribution, but with regional differences within South Siberia. The highest influx of West Eurasian mtDNAs was observed in populations from the East Sayan and Altai regions (from 12.5% to 34.5%), whereas in populations from the Baikal region this contribution was markedly lower (less than 10%). The considerable substructure within South Siberian haplogroups B, F, and G, together with the high degree of haplogroup C and D diversity revealed there, allows us to conclude that South Siberians carry the genetic imprint of early‐colonization phase of Eurasia. Statistical analyses revealed that South Siberian populations contain high levels of mtDNA diversity and high heterogeneity of mtDNA sequences among populations (Fst = 5.05%) that might be due to geography but not due to language and anthropological features.     

Northwest Siberian Khanty and Mansi in the junction of West and East Eurasian gene pools as revealed by uniparental markers

Northwest Siberia is geographically remote territory, which has been settled by indigenous human populations probably since the Upper Paleolithic. To investigate the genetic landscape of Northwest Siberians, we have analyzed mitochondrial and Y chromosome DNA polymorphisms of 169 unrelated individuals from Khanty and Mansi ethnic groups in Northwest Siberia. In addition, HVS-I sequences (N = 3522) and Y chromosome SNP data (N = 2175), obtained from the literature, were used to elucidate the genetic relationships among the North Eurasian populations. The results show clinal distributions of mtDNA and Y chromosome haplogroups along East-West axis of Northern Eurasia. In this context, the Ugric-speaking Khanty and Mansi appear as unique intermediate populations carrying Upper Paleolithic and more recent haplotypes typical for both West and East Eurasian gene pools. This admixture indicates that the Khanty and Mansi populations have resided in the contact zone of genetically distinguishable eastern and western Eurasia.     

Mitochondrial haplogroup C in ancient mitochondrial DNA from Ukraine extends the presence of East Eurasian genetic lineages in Neolithic Central and Eastern Europe

Recent studies of ancient mitochondrial DNA (mtDNA) lineages have revealed the presence of East Eurasian mtDNA haplogroups in the Central European Neolithic. Here we report the finding of East Eurasian lineages in ancient mtDNA from two Neolithic cemeteries of the North Pontic Region (NPR) in Ukraine. In our study, comprehensive haplotyping information was obtained for 7 out of 18 specimens. Although the majority of identified mtDNA haplogroups belonged to the traditional West Eurasian lineages of H and U, three specimens were determined to belong to the lineages of mtDNA haplogroup C. This find extends the presence of East Eurasian lineages in Neolithic Europe from the Carpathian Mountains to the northern shores of the Black Sea and provides the first genetic account of Neolithic mtDNA lineages from the NPR.     

Mitochondrial DNA and Y chromosome diversity and the peopling of the Americas: evolutionary and demographic evidence.

A number of important insights into the peopling of the New World have been gained through molecular genetic studies of Siberian and Native American populations. While there is no complete agreement on the interpretation of the mitochondrial DNA (mtDNA) and Y chromosome (NRY) data from these groups, several generalizations can be made. To begin with, the primary migration of ancestral Asians expanded from south-central Siberia into the New World and gave rise to ancestral Amerindians. The initial migration seems to have occurred between 20,000-15,000 calendar years before present (cal BP), i.e., before the emergence of Clovis lithic sites (13,350-12,895 cal BP) in North America. Because an interior route through northern North America was unavailable for human passage until 12,550 cal BP, after the last glacial maximum (LGM), these ancestral groups must have used a coastal route to reach South America by 14,675 cal BP, the date of the Monte Verde site in southern Chile. The initial migration appears to have brought mtDNA haplogroups A-D and NRY haplogroups P-M45a and Q-242/Q-M3 to the New World, with these genetic lineages becoming widespread in the Americas. A second expansion that perhaps coincided with the opening of the ice-free corridor probably brought mtDNA haplogroup X and NRY haplogroups P-M45b, C-M130, and R1a1-M17 to North and Central America. Finally, populations that formerly inhabited Beringia expanded into northern North America after the LGM, and gave rise to Eskimo-Aleuts and Na-Dené Indians.     

Genetic structure of circumpolar populations: a synthesis.

This article summarizes research conducted on the genetic structure of circumpolar populations of Alaska, Siberia, and the Aleutian Archipelago from 1977 to present. Three research programs are described: (1) 1977-1978, on the genetics of small, Yupik-speaking communities of St. Lawrence Island (Savoonga and Gambell) compared with Inupk communities of King Island and Wales, Alaska. This research, based on 25 standard protein and blood markers, demonstrated genetic discontinuity between the populations of the two language groups, patterns of gene flow, and the co-evolution of genes and languages in the Arctic. (2) 1989-1995, on the origins of Native American populations in Siberia. DNA was collected from two Evenki reindeer herding populations: Surinda and Poligus, a Ket community on the Stony Tunguska River, Sulamai, and a cattle-herding village from Gorno Altai, Mendur-Sokhon. Using an assortment of molecular markers, such as mitochondrial DNA, DNA fingerprints, and nonrecombining Y-chromosome markers, this research demonstrated the close genetic affinities among the Siberian and the Native American indigenous groups. (3) 1999 to present, on the origins of populations of the Aleutian Islands and their underlying genetic structure. A total of 11 inhabited islands from the Aleutian Archipelago were sampled and based on mtDNA sequences was shown to cluster tightly with Siberian Eskimo and Chukchi populations. Evidence of genetic drift through founder effect was demonstrated on Bering Island, where the D2 haplotype was fixed. Genetic discontinuity based on mtDNA was shown through SAMOVA between Kamchatkan and Aleut populations. An intimate relationship between geography and genetics through Mantel tests was observed for both Siberia (r = 0.55 P > 0.001) and the Aleutian Islands (r = 0.72 P > 0.000). The genetic structure of the populations of the Aleutian Archipelago most closely approximates the isolation-by-distance model.     

A comprehensive review of HVS-I mitochondrial DNA variation of 19 Iranian populations

Iran is located along the Central Asian corridor, a natural artery that has served as a cross-continental route since the first anatomically modern human populations migrated out of Africa. We compiled and reanalyzed the HVS-I (hypervariable segment-I) of 3840 mitochondrial DNA (mtDNA) sequences from 19 Iranian populations and from 26 groups from adjacent countries to give a comprehensive review of the maternal genetic variation and investigate the impact of historical events and cultural factors on the maternal genetic structure of modern Iranians. We conclude that Iranians have a high level of genetic diversity. Thirty-six haplogroups were observed in Iran's populations, and most of them belong to widespread West-Eurasian haplogroups, such as H, HV, J, N, T, and U. In contrast, the predominant haplogroups observed in most of the adjacent countries studied here are H, M, D, R, U, and C haplogroups. Using principal component analysis, clustering, and genetic distance-based calculations, we estimated moderate genetic relationships between Iranian and other Eurasian groups. Further, analyses of molecular variance and comparing geographic and genetic structures indicate that mtDNA HVS-I sequence diversity does not exhibit any sharp geographic structure in the country. Barring a few from some culturally distinct and naturally separated minorities, most Iranian populations have a homogenous maternal genetic structure.     

Mitochondrial DNA Sequence Diversity in a Sedentary Population from Egypt

The mitochondrial DNA (mtDNA) diversity of 58 individuals from Upper Egypt, more than half (34 individuals) from Gurna, whose population has an ancient cultural history, were studied by sequencing the control‐region and screening diagnostic RFLP markers. This sedentary population presented similarities to the Ethiopian population by the L1 and L2 macrohaplogroup frequency (20.6%), by the West Eurasian component (defined by haplogroups H to K and T to X) and particularly by a high frequency (17.6%) of haplogroup M1. We statistically and phylogenetically analysed and compared the Gurna population with other Egyptian, Near East and sub‐Saharan Africa populations; AMOVA and Minimum Spanning Network analysis showed that the Gurna population was not isolated from neighbouring populations. Our results suggest that the Gurna population has conserved the trace of an ancestral genetic structure from an ancestral East African population, characterized by a high M1 haplogroup frequency. The current structure of the Egyptian population may be the result of further influence of neighbouring populations on this ancestral population.     

Mitochondrial DNA and Y-Chromosome Variation in the Caucasus

We have analyzed mtDNA HVI sequences and Y chromosome haplogroups based on 11 binary markers in 371 individuals, from 11 populations in the Caucasus and the neighbouring countries of Turkey and Iran. Y chromosome haplogroup diversity in the Caucasus was almost as high as in Central Asia and the Near East, and significantly higher than in Europe. More than 27% of the variance in Y‐haplogroups can be attributed to differences between populations, whereas mtDNA showed much lower heterogeneity between populations (less then 5%), suggesting a strong influence of patrilocal social structure. Several groups from the highland region of the Caucasus exhibited low diversity and high differentiation for either or both genetic systems, reflecting enhanced genetic drift in these small, isolated populations. Overall, the Caucasus groups showed greater similarity with West Asian than with European groups for both genetic systems, although this similarity was much more pronounced for the Y chromosome than for mtDNA, suggesting that male‐mediated migrations from West Asia have influenced the genetic structure of Caucasus populations.     

Mitochondrial DNA Variability in Slovaks, with Application to the Roma Origin

To gain insight into the mitochondrial gene pool diversity of European populations, we studied mitochondrial DNA (mtDNA) variability in 207 subjects from western and eastern areas of Slovakia. Sequencing of two hypervariable segments, HVS I and HVS II, in combination with screening of coding region haplogroup-specific RFLP-markers, revealed that the majority of Slovak mtDNAs belong to the common West Eurasian mitochondrial haplogroups (HV, J, T, U, N1, W, and X). However, a few sub-Saharan African (L2a) mtDNAs were detected in a population from eastern part of Slovakia. In addition, about 3% of mtDNAs from eastern Slovakia encompass Roma-specific lineages. By means of complete mtDNA sequencing we demonstrate here that the Roma-specific M-lineages observed in gene pools of different Slavonic populations (Slovaks, Poles and Russians), belong to Indian-specific haplogroups M5a1 and M35. Moreover, we show that haplogroup J lineages found in gene pools of the Roma and some Slavonic populations (Czechs and Slovaks) belong to new subhaplogroup J1a, which is defined by coding region mutation at position 8460.     

Admixture, migrations, and dispersals in Central Asia: evidence from maternal DNA lineages

Mitochondrial DNA (mtDNA) lineages of 232 individuals from 12 Central Asian populations were sequenced for both control region hypervariable segments, and additional informative sites in the coding region were also determined. Most of the mtDNA lineages belong to branches of the haplogroups with an eastern Eurasian (A, B, C, D, F, G, Y, and M haplogroups) or a western Eurasian (HV, JT, UK, I, W, and N haplogroups) origin, with a small fraction of Indian M lineages. This suggests that the extant genetic variation found in Central Asia is the result of admixture of already differentiated populations from eastern and western Eurasia. Nonetheless, two groups of lineages, D4c and G2a, seem to have expanded from Central Asia and might have their Y-chromosome counterpart in lineages belonging to haplotype P(xR1a). The present results suggest that the mtDNA found out of Africa might be the result of a maturation phase, presumably in the Middle East or eastern Africa, that led to haplogroups M and N, and subsequently expanded into Eurasia, yielding a geographically structured group of external branches of these two haplogroups in western and eastern Eurasia, Central Asia being a contact zone between two differentiated groups of peoples.     

High-Resolution mtDNA Studies of the Indian Population: Implications for Palaeolithic Settlement of the Indian Subcontinent

The population of the Indian subcontinent represents a very complex social and cultural structure. Occupying a geographically central position for the early modern human migrations, indications are that the founder group that migrated out of East Africa also reached India. In the present study we used the twin strategy of mapping the whole mitochondrial DNA (mtDNA) using the standard 14 restriction enzymes, and sequencing the non‐transcribed HVSI region, to derive maximum maternal lineages from a sample of non‐tribal Indians. The essential features of the reduced median network of the two datasets were the same. Both showed two demographic expansions of two major haplogroups, ‘M’ and ‘N’. The reduced median network was drawn with inputs from other studies on the Indian population, and correlated with data from other ethnic populations. The coalescence time of expansions and genetic diversity were estimated. A reduced median network was also drawn combining data from studies on Africans, Southeast Asians and West‐Eurasians, tracing the migration of ‘M’ from East Africa to India. A time estimate of the migration of major mtDNA haplogroups from Africa was attempted. The comparison of a set of Indian maternal lineages belonging to different geographical regions of the country, with other populations revealed the in‐situ differentiation and antiquity of the Indian population. Our analysis places the ‘southern route’ migration as the source of haplogroup ‘M’. Multiple migrations might have brought the other major haplogroups, ‘N’ and ‘R’, found in our sample to India. Archaeological evidence of modern humans in the subcontinent supports this mtDNA study.     

Review of Croatian genetic heritage as revealed by mitochondrial DNA and Y chromosomal lineages

The aim of this review is to summarize the existing data collected in high-resolution phylogenetic studies of mitochondrial DNA and Y chromosome variation in mainland and insular Croatian populations. Mitochondrial DNA polymorphisms were explored in 721 individuals by sequencing mtDNA HVS-1 region and screening a selection of 24 restriction fragment length polymorphisms (RFLPs), diagnostic for main Eurasian mtDNA haplogroups. Whereas Y chromosome variation was analyzed in 451 men by using 19 single nucleotide polymorphism (SNP)/indel and 8 short tandem repeat (STR) loci. The phylogeography of mtDNA and Y chromosome variants of Croatians can be adequately explained within typical European maternal and paternal genetic landscape, with the exception of mtDNA haplogroup F and Y-chromosomal haplogroup P* which indicate a connection to Asian populations. Similar to other European and Near Eastern populations, the most frequent mtDNA haplogroups in Croatians were H (41.1%), U5 (10.3%), and J (9.7%). The most frequent Y chromosomal haplogroups in Croatians, I-P37 (41.7%) and R1a-SRY1532 (25%), as well as the observed structuring of Y chromosomal variance reveal a clearly evident Slavic component in the paternal gene pool of contemporary Croatian men. Even though each population and groups of populations are well characterized by maternal and paternal haplogroup distribution, it is important to keep in mind that linking phylogeography of various haplogroups with known historic and prehistoric scenarios should be cautiously performed.     

MtDNA Profile of West Africa Guineans: Towards a Better Understanding of the Senegambia Region

The matrilineal genetic composition of 372 samples from the Republic of Guiné‐Bissau (West African coast) was studied using RFLPs and partial sequencing of the mtDNA control and coding region. The majority of the mtDNA lineages of Guineans (94%) belong to West African specific sub‐clusters of L0‐L3 haplogroups. A new L3 sub‐cluster (L3h) that is found in both eastern and western Africa is present at moderately low frequencies in Guinean populations. A non‐random distribution of haplogroups U5 in the Fula group, the U6 among the “Brame” linguistic family and M1 in the Balanta‐Djola group, suggests a correlation between the genetic and linguistic affiliation of Guinean populations. The presence of M1 in Balanta populations supports the earlier suggestion of their Sudanese origin. Haplogroups U5 and U6, on the other hand, were found to be restricted to populations that are thought to represent the descendants of a southern expansion of Berbers. Particular haplotypes, found almost exclusively in East‐African populations, were found in some ethnic groups with an oral tradition claiming Sudanese origin.     

Mitochondrial DNA haplogroup R predicts survival advantage in severe sepsis in the Han population

PurposeTo determine whether the main mitochondrial DNA (mtDNA) haplogroups of the Han people have an impact on long-term clinical outcome.MethodsWe prospectively studied 181 individuals who were sequentially admitted to the intensive care unit. Demographic and clinical data were recorded along with clinical outcome over 180 days. Follow-up was completed for all study participants. We then determined the mtDNA haplogroups of the patients and 570 healthy, age-matched Han people from Zhejiang province, Southeast China, by analyzing sequences of hypervariable mtDNA segments and testing diagnostic polymorphisms in the mtDNA coding region with DNA probes.ResultThe frequency of the main subhaplogroups of the Han population in the study cohort did not differ significantly from the control group. mtDNA haplogroup R, one of the three main mtDNA haplogroups of the Han people, was a strong independent predictor for the outcome of severe sepsis, conferring a 4.68-fold (95% CI 1.903–10.844, P = 0.001) increased chance of survival at 180 days compared with those without the haplogroup R.ConclusionIn the Han population, mtDNA haplogroup R was a strong independent predictor for the outcome of severe sepsis, conferring an increased chance of long-term survival compared with individuals without the R haplogroup.     

Y-chromosome haplogroup N dispersals from south Siberia to Europe

In order to reconstruct the history of Y-chromosome haplogroup (hg) N dispersals in north Eurasia, we have analyzed the diversity of microsatellite (STR) loci within two major hg N clades, N2 and N3, in a total sample of 1,438 males from 17 ethnic groups, mainly of Siberian and Eastern European origin. Based on STR variance analysis we observed that hg N3a is more diverse in Eastern Europe than in south Siberia. However, analysis of median networks showed that there are two STR subclusters of hg N3a, N3a1 and N3a2, that are characterized by different genetic histories. Age calculation of STR variation within subcluster N3a1 indicated that its first expansion occurred in south Siberia [approximately 10,000 years (ky)] and then this subcluster spread into Eastern Europe where its age is around 8 ky ago. Meanwhile, younger subcluster N3a2 originated in south Siberia (probably in the Baikal region) approximately 4 ky ago. Median network and variance analyses of STR haplotypes suggest that south Siberian N3a2 haplotypes spread further into Volga-Ural region undergoing serial bottlenecks. In addition, median network analysis of STR data demonstrates that haplogroup N2-A is represented by two subclusters, showing recent expansion times. The data obtained allow us to suggest Siberian origin of haplogroups N3 and N2 that are currently widespread in some populations of Eastern Europe.     

Joining the Pillars of Hercules: mtDNA Sequences Show Multidirectional Gene Flow in the Western Mediterranean

Phylogenetic analysis of mitochondrial DNA (mtDNA) performed in Western Mediterranean populations has shown that both shores share a common set of mtDNA haplogroups already found in Europe and the Middle East. Principal co‐ordinates of genetic distances and principal components analyses based on the haplotype frequencies show that the main genetic difference is attributed to the higher frequency of sub‐Saharan L haplogroups in NW Africa, showing some gene flow across the Sahara desert, with a major impact in the southern populations of NW Africa. The AMOVA demonstrates that SW European populations are highly homogeneous whereas NW African populations display a more heterogeneous genetic pattern, due to an east‐west differentiation as a result of gene flow coming from the East. Despite the shared haplogroups found in both areas, the European V and the NW African U6 haplogroups reveal the traces of the Mediterranean Sea permeability to female migrations, and allowed for determination and quantification of the genetic contribution of both shores to the genetic landscape of the geographic area. Comparison of mtDNA data with autosomal markers and Y‐chromosome lineages, analysed in the same populations, shows a congruent pattern, although female‐mediated gene flow seems to have been more intense than male‐mediated gene flow.     

Mitochondrial DNA sequence diversity in two groups of Italian Veneto speakers from Veneto

Although frequencies of mitochondrial DNA (mtDNA) haplogroups in the different European populations are rather homogenous, there are a few European populations or linguistic isolates that show different mtDNA haplogroup distributions; examples are the Saami and Ladin speakers from the eastern Italian Alps.
MtDNA sequence diversity was analysed from subjects from two villages in Veneto. The first, Posina, is situated in the Venetian Alps near Vicenza. The second, Barco di Pravisdomini is a village on the plains near Venice. In spite of their common Veneto dialect, the two group populations have not preserved a genetic homogeneity; particularly, they show differences in T and J haplogroups frequencies. MtDNA diversity in these two groups seems to depend more on their geographic situation.     

Mitochondrial genome diversity at the Bering Strait area highlights prehistoric human migrations from Siberia to northern North America

The patterns of prehistoric migrations across the Bering Land Bridge are far from being completely understood: there still exists a significant gap in our knowledge of the population history of former Beringia. Here, through comprehensive survey of mitochondrial DNA genomes retained in ‘relic'' populations, the Maritime Chukchi, Siberian Eskimos, and Commander Aleuts, we explore genetic contribution of prehistoric Siberians/Asians to northwestern Native Americans. Overall, 201 complete mitochondrial sequences (52 new and 149 published) were selected in the reconstruction of trees encompassing mtDNA lineages that are restricted to Coastal Chukotka and Alaska, the Canadian Arctic, Greenland, and the Aleutian chain. Phylogeography of the resulting mtDNA genomes (mitogenomes) considerably extends the range and intrinsic diversity of haplogroups (eg, A2a, A2b, D2a, and D4b1a2a1) that emerged and diversified in postglacial central Beringia, defining independent origins of Neo-Eskimos versus Paleo-Eskimos, Aleuts, and Tlingit (Na-Dene). Specifically, Neo-Eskimos, ancestral to modern Inuit, not only appear to be of the High Arctic origin but also to harbor Altai/Sayan-related ancestry. The occurrence of the haplogroup D2a1b haplotypes in Chukotka (Sireniki) introduces the possibility that the traces of Paleo-Eskimos have not been fully erased by spread of the Neo-Eskimos or their descendants. Our findings are consistent with the recurrent gene flow model of multiple streams of expansions to northern North America from northeastern Eurasia in late Pleistocene–early Holocene.     

Complete mitogenomes document substantial genetic contribution from the Eurasian Steppe into northern Pakistani Indo-Iranian speakers

To elucidate whether Bronze Age population dispersals from the Eurasian Steppe to South Asia contributed to the gene pool of Indo-Iranian-speaking groups, we analyzed 19,568 mitochondrial DNA (mtDNA) sequences from northern Pakistani and surrounding populations, including 213 newly generated mitochondrial genomes (mitogenomes) from Iranian and Dardic groups, both speakers from the ancient Indo-Iranian branch in northern Pakistan. Our results showed that 23% of mtDNA lineages with west Eurasian origin arose in situ in northern Pakistan since ~5000 years ago (kya), a time depth very close to the documented Indo-European dispersals into South Asia during the Bronze Age. Together with ancient mitogenomes from western Eurasia since the Neolithic, we identified five haplogroups (~8.4% of maternal gene pool) with roots in the Steppe region and subbranches arising (age ~5–2 kya old) in northern Pakistan as genetic legacies of Indo-Iranian speakers. Some of these haplogroups, such as W3a1b that have been found in the ancient samples from the late Bronze Age to the Iron Age period individuals of Swat Valley northern Pakistan, even have sub-lineages (age ~4 kya old) in the southern subcontinent, consistent with the southward spread of Indo-Iranian languages. By showing that substantial genetic components of Indo-Iranian speakers in northern Pakistan can be traced to Bronze Age in the Steppe region, our study suggests a demographic link with the spread of Indo-Iranian languages, and further highlights the corridor role of northern Pakistan in the southward dispersal of Indo-Iranian-speaking groups.Subject terms: Evolutionary biology, Molecular biology