Genetic characterization of the Makrani people of Pakistan from mitochondrial DNA control-region data

Ethnogenetic analysis reveals the Bronze Age genetic affiliation of Yashkuns with West Eurasians

Focusing on the Yashkun population of Gilgit-Baltistan, an administrative territory in northern Pakistan, our study investigated mtDNA haplotypes as indicators of ancient gene flow and genetic diversity. Genomic DNA was extracted and evaluated for quality using agarose gel electrophoresis. The complete control region of mtDNA (nt 16024-576) was amplified via PCR, and sequencing was performed using the Big Dye Terminator Kit on an Applied Biosystems Genetic Analyzer. Sequence alignment and analysis were conducted using Geneious software and the Cambridge Reference Sequence (CRS). Haplotypes were classified into phylogenetic haplogroups with HaploGrep 2 and the Mito Tool platform, referencing PhyloTree version 17. The study identified considerable genetic diversity within the Yashkun population, with major macro-haplogroups H (29.47%), T (13.68%), and M (13.68%) representing West Eurasian, South Asian, and East & Southeast Asian lineages. The highest genetic affiliation was with West Eurasian haplogroups (73.68%), indicating significant gene flow from the West. The population's genetic distinctiveness was highlighted by the calculation of forensic genetic characteristics, such as power of discrimination (PD = 0.9975), genetic diversity (GD = 0.9865), and random match probability (RMP = 0.01349). The knowledge of the genetic makeup of the Yashkun population was further improved by the inference of phylogenetic links made possible by comparisons with other indigenous groups in Pakistan. These results attest to the reliability of mtDNA data for forensic applications such as ancestry tracking, individual identification, and additions to national forensic databases. The inclusion of mtDNA analysis in forensic science is strongly supported by this groundbreaking study on the Yashkun community, especially in areas like Pakistan with diverse genetic histories.

Genetic diversity and forensic application of Y-filer STRs in four major ethnic groups of Pakistan

17 Y-chromosomal STRs which are part of the Yfiler Amplification Kit were investigated in 493 unrelated Pakistani individuals belonging to the Punjabi, Sindhi, Baloch, and Pathan ethnic groups. We have assessed the forensic parameters and population genetic structure for each group. Among the 493 unrelated individuals from four ethnic groups (128 Baloch, 122 Pathan, 108 Punjabi, and 135 Sindhi), 82 haplotypes were observed with haplotype diversity (HD) of 0.9906 in Baloch, 102 haplotypes with HD value of 0.9957 in Pathans, 80 haplotypes with HD value of 0.9924 in Punjabi, and 105 haplotypes with HD value of 0.9945 in the Sindhi population. The overall gene diversity for Baloch, Pathan, Punjabi, and Sindhi populations was 0.6367, 0.6479, 0.6657, and 0.6112, respectively. The results had shown us that Pakistani populations do not have a unique set of genes but share the genetic affinity with regional (Central Asia and Northern India) populations. The observed low gene diversity (heterozygosity) values may be because of endogamy trends and this observation is equally supported by the results of forensic parameters which are mostly static across 4 combinations (minimal STRs, extended 11 Y-STRs, Powerplex 12 Y System, and Yfiler 17 Y-STRs) of STRs in these four populations.

Contrasting maternal and paternal genetic histories among five ethnic groups from Khyber Pakhtunkhwa, Pakistan

Northwest Pakistan has served as a point of entry to South Asia for different populations since ancient times. However, relatively little is known about the population genetic history of the people residing within this region. To better understand human dispersal in the region within the broader history of the subcontinent, we analyzed mtDNA diversity in 659 and Y-chromosome diversity in 678 individuals, respectively, from five ethnic groups (Gujars, Jadoons, Syeds, Tanolis and Yousafzais), from Swabi and Buner Districts, Khyber Pakhtunkhwa Province, Pakistan. The mtDNAs of all individuals were subject to control region sequencing and SNP genotyping, while Y-chromosomes were analyzed using 54 SNPs and 19 STR loci. The majority of the mtDNAs belonged to West Eurasian haplogroups, with the rest belonging to either South or East Asian lineages. Four of the five Pakistani populations (Gujars, Jadoons, Syeds, Yousafzais) possessed strong maternal genetic affinities with other Pakistani and Central Asian populations, whereas one (Tanolis) did not. Four haplogroups (R1a, R1b, O3, L) among the 11 Y-chromosome lineages observed among these five ethnic groups contributed substantially to their paternal genetic makeup. Gujars, Syeds and Yousafzais showed strong paternal genetic affinities with other Pakistani and Central Asian populations, whereas Jadoons and Tanolis had close affinities with Turkmen populations from Central Asia and ethnic groups from northeast India. We evaluate these genetic data in the context of historical and archeological evidence to test different hypotheses concerning their origins and biological relationships.

Current pool of ultimate collection of mitochondrial DNA from remnants of Kalash

The mitochondrial DNA (mtDNA) complete control region coverage of 111 individuals from Kalash population of Pakistan has been presented for forensic applications and to infer their genetic parameters. We detected in total 14 different haplotypes with only five unique and nine shared by more than one individual. This population has come up with quite lower haplotype diversity (0.8393) and very higher random match probability (0.1682), and ultimately lower power of discrimination (0.832). Additionally, haplogroup distribution reveals the genetic ancestry of Kalash, mainly from West Eurasia (98.8%) and very little from South Asia (0.9%). Neither African lineages nor East Asian genetic segments were detected among these Kalash. This study will contribute to the database development for forensic applications as well as to track the evolutionary highlights of this ethnic group.

Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch

Skin pigmentation is one of the most prominent and variable phenotypes in humans. We compared the alleles of 163 SNPs and indels from the Human Pigmentation (HuPi) AmpliSeq™ Custom panel, and biogeographic ancestry with the quantitative skin pigmentation levels on the upper arm, lower arm, and forehead of 299 Pakistani individuals from three subpopulations: Baloch, Pashtun, and Punjabi. The biogeographic ancestry of each individual was estimated using the Precision ID Ancestry Panel. All individuals were mainly of mixed South-Central Asian and European ancestry. However, the Baloch individuals also had an average proportion of Sub-Saharan African ancestry of approximately 10%, whereas it was <1% in the Punjabi and Pashtun individuals. The pairwise genetic distances between the Pashtun, Punjabi, and Baloch subpopulations based on the ancestry markers were statistically significantly different. Individuals from the Pashtun subpopulation had statistically significantly lower skin pigmentation than individuals from the Punjabi and Baloch subpopulations (p < 0.05). The proportions of European and Sub-Saharan African ancestry and five SNPs (rs1042602, rs10831496, rs1426654, rs16891982, and rs12913832) were statistically significantly associated with skin pigmentation at either the upper arm, lower arm or forehead in the Pakistani population after correction for multiple testing (p < 10⁻³). A model based on four of these SNPs (rs1426654, rs1042602, rs16891982, and rs12913832) explained 33% of the upper arm skin pigmentation. The four SNPs and the proportions of European and Sub-Saharan African ancestry explained 37% of the upper arm skin pigmentation. Our results indicate that the four likely causative SNPs, rs1426654, rs1042602, rs16891982, and rs12913832 located in SLC24A5, TYR, SLC45A2, and HERC2, respectively, are essential for skin color variation in the admixed Pakistani subpopulations.

Mitochondrial DNA diversity in the Khattak and Kheshgi of the Peshawar Valley, Pakistan

The strategic location of Pakistan and its presence at the crossroads of Asia has resulted in it playing a central role in both prehistoric and historic human migratory events, thereby linking and facilitating contacts between the inhabitants of the Middle East, Central Asia, China and South Asia. Despite the importance of this region and its inhabitants for our understanding of modern human origins and population dispersals, the nature of mitochondrial DNA (mtDNA) variation among members of the myriad populations of this area has largely been unexplored. Here, we report mtDNA control region sequences in 58 individuals from the Khattak and the Kheshgi, two major Pakhtun tribes residing within the Peshawar Valley of northwestern Pakistan. The results reveal that these ethnic groups are genetically heterogeneous, having 55.7% West Eurasian, 33.9% South Asian and 10.2% East Asian haplogroups. The genetic diversity observed for the Kheshgi was somewhat higher than that of the Khattak. A multidimensional scaling plot based on haplogroup frequencies for the Khattak, Kheshgi and neighboring populations indicates that the Khattak have close affinities with Baluch, Uzbek and Kazak populations but are only distantly related to the Kheshgi and other Pakistani populations. By contrast, the Kheshgi cluster closely with other Pakhtun or Pathan populations of Pakistan, suggesting a possible common maternal gene pool shared amongst them. These mtDNA data allow us to begin reconstructing the origins of the Khattak and Kheshgi and describe their complex interactions with populations from the surrounding regions.

Mitochondrial control region diversity in Sindhi ethnic group of Pakistan

The entire mitochondrial DNA control region (nt 16024-576) of 88 unrelated individuals of Sindhi ethnic group residing in different parts of Sindh province of Pakistan was sequenced. Out of 66 different observed haplotypes 50 were unique and 16 were shared by more than one individual. Results showed admixture of mtDNA pool constituting the haplogroups derived mainly from South Asia (47.6%) and West Eurasian (35.7%) whereas the contribution of the African haplogroup was very small (2.4%). High values of genetic diversity (0.992), power of discrimination (0.981) and low value of random match probability (0.018) indicates that mtDNA analysis for this population can effectively be used for forensic casework. The results are valuable contribution towards building mtDNA population variation database for this particular ethnic group from Pakistan.

Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia

Background

From a mtDNA dominant perspective, the exit from Africa of modern humans to colonize Eurasia occurred once, around 60 kya, following a southern coastal route across Arabia and India to reach Australia short after. These pioneers carried with them the currently dominant Eurasian lineages M and N. Based also on mtDNA phylogenetic and phylogeographic grounds, some authors have proposed the coeval existence of a northern route across the Levant that brought mtDNA macrohaplogroup N to Australia. To contrast both hypothesis, here we reanalyzed the phylogeography and respective ages of mtDNA haplogroups belonging to macrohaplogroup M in different regions of Eurasia and Australasia.

Results

The macrohaplogroup M has a historical implantation in West Eurasia, including the Arabian Peninsula. Founder ages of M lineages in India are significantly younger than those in East Asia, Southeast Asia and Near Oceania. Moreover, there is a significant positive correlation between the age of the M haplogroups and its longitudinal geographical distribution. These results point to a colonization of the Indian subcontinent by modern humans carrying M lineages from the east instead the west side.

Conclusions

The existence of a northern route, previously proposed for the mtDNA macrohaplogroup N, is confirmed here for the macrohaplogroup M. Both mtDNA macrolineages seem to have differentiated in South East Asia from ancestral L3 lineages. Taking this genetic evidence and those reported by other disciplines we have constructed a new and more conciliatory model to explain the history of modern humans out of Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0816-8) contains supplementary material, which is available to authorized users.

EMPOP-quality mtDNA control region sequences from Kashmiri of Azad Jammu & Kashmir, Pakistan

The mitochondrial DNA (mtDNA) control region (nucleotide position 16024-576) sequences were generated through Sanger sequencing method for 317 self-identified Kashmiris from all districts of Azad Jammu & Kashmir Pakistan. The population sample set showed a total of 251 haplotypes, with a relatively high haplotype diversity (0.9977) and a low random match probability (0.54%). The containing matrilineal lineages belonging to three different phylogeographic origins of Western Eurasian (48.9%), South Asian (47.0%) and East Asian (4.1%). The present study was compared to previous data from Pakistan and other worldwide populations (Central Asia, Western Asia, and East & Southeast Asia). The dataset is made available through EMPOP under accession number EMP00679 and will serve as an mtDNA reference database in forensic casework in Pakistan.

Genetic analysis of mitochondrial DNA control region variations in four tribes of Khyber Pakhtunkhwa, Pakistan

Due to its geo strategic position at the crossroad of Asia, Pakistan has gained crucial importance of playing its pivotal role in subsequent human migratory events, both prehistoric and historic. This human movement became possible through an ancient overland network of trails called "The Silk Route" linking Asia Minor, Middle East China, Central Asia and Southeast Asia. This study was conducted to analyze complete mitochondrial control region samples of 100 individuals of four major Pashtun tribes namely, Bangash, Khattak, Mahsuds and Orakzai in the province of Khyber Pakhtunkhwa, Pakistan. All Pashtun tribes revealed high genetic diversity which is comparable to the other Central Asian, Southeast Asian and European populations. The configuration of genetic variation and heterogeneity further unveiled through Multidimensional Scaling, Principal Component Analysis and phylogenetic analysis. The results revealed that Pashtun are the composite mosaic of West Eurasian ancestry of numerous geographic origin. They received substantial gene flow during different invasive movements and have a high element of the Western provenance. The most common haplogroups reported in this study are: South Asian haplogroups M (28%) and R (8%); whereas, West Asians haplogroups are present, albeit in high frequencies (67%) and widespread over all; HV (15%), U (17%), H (9%), J (8%), K (8%), W (4%), N (3%) and T (3%). Moreover, we linked the unexplored genetic connection between Ashkenazi Jews and Pashtun. The presence of specific haplotypes J1b (4%) and K1a1b1a (5%) pointed to a genetic connection of Jewish conglomeration in Khattak tribe. This was a result of an ancient genetic influx in the early Neolithic period that led to the formation of a diverse genetic substratum in present day Pashtun.