Indigenous and foreign Y-chromosomes characterize the Lingayat and Vokkaliga populations of Southwest India

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Forensic and genetic landscape explorations of Chinese Kyrgyz group based on autosomal SNPs, Y-chromosomal SNPs and STRs

To assess the effect of population genetic polymorphism on forensic research, we investigated the genetic polymorphisms of Chinese Kyrgyz group (n = 98) and evaluated forensic application values in Chinese Kyrgyz group and other 26 reference populations at 90 autosomal SNPs, and then combined with 34 SNPs and 37 STRs on Y chromosome to reveal the genetic background of Kyrgyz group in autosomal and Y-chromosomal inheritances, respectively. The 90 autosomal SNPs and 34 Y-chromosomal SNPs were sequenced base on next generation sequencing technology, and 37 Y-chromosomal STRs were analyzed by capillary electrophoresis platform. The results showed that cumulative power of discrimination and cumulative power of exclusion of 90 autosomal SNPs in the panel met the application need of forensic genetics in Kyrgyz group. The forensic effectivenesses of the panel were high in all 27 populations, although there were genetic differences among these populations. The forensic effectiveness of the panel was relatively higher in the European populations, but relatively lower in the African populations. The population genetic results indicated that the Kyrgyz group had the relatively closer genetic relationships with the reference East Asian populations at autosomal SNPs, and there were gene exchanges between the Kyrgyz group and East Asian, European populations based on the analytical results of autosomal SNPs, Y-chromosomal SNPs and STRs.

Genetic insights into the paternal admixture history of Chinese Mongolians via high-resolution customized Y-SNP SNaPshot panels

The Mongolian people, one of the Mongolic-speaking populations, are native to the Mongolian Plateau in North China and southern Siberia. Many ancient DNA studies recently reported extensive population transformations during the Paleolithic to historic periods in this region, while little is known about the paternal genetic legacy of modern geographically different Mongolians. Here, we genotyped 215 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) and 37 Y-chromosomal short tandem repeats (Y-STRs) among 679 Mongolian individuals from Hohhot, Hulunbuir, and Ordos in North China using the AGCU Y37 kit and our developed eight Y-SNP SNaPshot panels (including two panels first reported herein). The C-M130 Y-SNP SNaPshot panel defines 28 subhaplogroups, and the N/O/Q complementary Y-SNP SNaPshot panel defines 30 subhaplogroups of N1b-F2930, N1a1a1a1a3-B197, Q-M242, and O2a2b1a1a1a4a-CTS4658, which improved the resolution our developed Y-SNP SNaPshot panel set and could be applied for dissecting the finer-scale paternal lineages of Mongolic speakers. We found a strong association between Mongolian-prevailing haplogroups and some observed microvariants among the newly generated Y-STR haplotype data, suggesting the possibility of haplogroup prediction based on the distribution of Y-STR haplotypes. We identified three main ancestral sources of the observed Mongolian-dominant haplogroups, including the local lineage of C2*-M217 and incoming lineages from other regions of southern East Asia (O2*-M122, O1b*-P31, and N1*-CTS3750) and western Eurasia (R1*-M173). We also observed DE-M145, D1*-M174, C1*-F3393, G*-M201, I-M170, J*-M304, L-M20, O1a*-M119, and Q*-M242 at relatively low frequencies (< 5.00%), suggesting a complex admixture history between Mongolians and other incoming Eurasians from surrounding regions. Genetic clustering analyses indicated that the studied Mongolians showed close genetic affinities with other Altaic-speaking populations and Sinitic-speaking Hui people. The Y-SNP haplotype/haplogroup-based genetic legacy not only revealed that the stratification among geographically/linguistically/ethnically different Chinese populations was highly consistent with the geographical division and language classification, but also demonstrated that patrilineal genetic materials could provide fine-scale genetic structures among geographically different Mongolian people, suggesting that our developed high-resolution Y-SNP SNaPshot panels have the potential for forensic pedigree searches and biogeographical ancestry inference.

Genetic polymorphisms and phylogenetic analyses of the Ü-Tsang Tibetan from Lhasa based on 30 slowly and moderately mutated Y-STR loci

As a result of the expansion of old Tibet on the Qinghai-Tibet Plateau, Tibetans diverged into three main branches, Ü-Tsang, Amdo, and Kham Tibetan. Ü-Tsang Tibetans are geographically distributed across the wide central and western portions of the Qinghai-Tibet Plateau while Lhasa is the central gathering place for Tibetan culture. The AGCU Y30, a 6-dye fluorescence kit including 30 slowly and moderately mutated Y-STR loci, has been validated for its stability and sensitivity in different biomaterials and diverse Chinese populations (Han and other minorities), and widely used in the practical work of forensic science. However, the 30 Y-STR profiling of Tibetan, especially for Ü-Tsang Tibetan, were insufficient. We utilized the AGCU Y30 to genotype 577 Ü-Tsang Tibetan unrelated males from Lhasa in the Tibet Autonomous Region of China to fill up the full and accurate Y-STR profiles. A total of 552 haplotypes were observed, 536 (97.10%) of which were unique. One hundred and ninety-four alleles were observed at 26 single copy loci and the allelic frequencies ranged from 0.0017 to 0.8180. For the two multi-copy loci DYS385a/b and DYS527a/b, 64 and 36 allelic combinations were observed, respectively. The gene diversity (GD) values ranged from 0.3079 at DYS391 to 0.9142 at DYS385a/b and the overall haplotype diversity (HD) was 0.9998, and its discrimination capacity (DC) was 0.9567. The population genetic analyses demonstrated that Lhasa Ü-Tsang Tibetan had close relationships with other Tibetan populations from Tibet and Qinghai, especially with Ü-Tsang Tibetan. From the perspective of Y haplogroups, the admixture of the southward Qiang people with dominant haplogroup O-M122 and the northward migrations of the initial settlers of East Asia with haplogroup D-M175 hinted the Sino-Tibetan homologous, thus, we could not ignore the gene flows with other Sino-Tibetan populations, especially for Han Chinese, to characterize the forensic genetic landscape of Tibetan.

Forensic characteristics and phylogenetic analyses of one branch of Tai-Kadai language-speaking Hainan Hlai (Ha Hlai) via 23 autosomal STRs included in the Huaxia™ Platinum System

BACKGROUND

Hainan Island, located in the South China Sea and separated from the Leizhou Peninsula by Qiongzhou Strait, is the second largest island after Taiwan in China. With the expansion of Han Chinese and the gradual formation of "South Hlai and North Han", nowadays, Hainan Hlai is the second largest population after Han Chinese in Hainan Island. Ha Hlai, distributed in southwest and southern Hainan Island, is the dominant branch of Hlai and speaks Ha localism.

METHODS

We utilized the Huaxia™ Platinum PCR Amplification System (including 23 autosomal STRs and 2 sex-linked markers) to obtain the first STR profiling batch of 657 Ha Hlai individuals (497 males and 160 females). In order to explore the genetic relationships between the studied Ha Hlai and other reference populations with different language families, population genetic analyses, including PCA, MDS, STRUCTURE, and phylogenetic analysis, were conducted based upon the raw data and allelic frequencies of the polymorphic autosomal STR markers.

RESULTS

In total, 271 distinct alleles were observed at the 23 STR loci. The number of diverse alleles ranged from 7 at TPOX locus to 23 at FGA locus, and the allelic frequencies varied from 0.0008 to 0.5533. In addition, the CPE and CPD were 1-7.39 × 10-10 and 1-3.13 × 10-28 , respectively. The phylogenetic analyses indicated that Ha Hlai is a Tai-Kadai language-speaking and relatively isolated population which has a close genetic and geographical relationship with Hainan Hlai, and M95 is the dominant haplogroup in Ha Hlai (56.18%).

CONCLUSION

The 23 autosomal STR genetic markers were highly polymorphic as well as potentially useful for forensic applications in Hainan Ha Hlai population. The phylogenetic analyses demonstrated that small geographic scale gene flows could not be ignored and the shaping of the unique gene pool for each population was the combination effects of geographic, language, and cultural isolations.

The Geographic Origins of Ethnic Groups in the Indian Subcontinent: Exploring Ancient Footprints with Y-DNA Haplogroups

Several studies have evaluated the movements of large populations to the Indian subcontinent; however, the ancient geographic origins of smaller ethnic communities are not clear. Although historians have attempted to identify the origins of some ethnic groups, the evidence is typically anecdotal and based upon what others have written before. In this study, recent developments in DNA science were assessed to provide a contemporary perspective by analyzing the Y chromosome haplogroups of some key ethnic groups and tracing their ancient geographical origins from genetic markers on the Y-DNA haplogroup tree. A total of 2,504 Y-DNA haplotypes, representing 50 different ethnic groups in the Indian subcontinent, were analyzed. The results identified 14 different haplogroups with 14 geographic origins for these people. Moreover, every ethnic group had representation in more than one haplogroup, indicating multiple geographic origins for these communities. The results also showed that despite their varied languages and cultural differences, most ethnic groups shared some common ancestors because of admixture in the past. These findings provide new insights into the ancient geographic origins of ethnic groups in the Indian subcontinent. With about 2,000 other ethnic groups and tribes in the region, it is expected that more scientific discoveries will follow, providing insights into how, from where, and when the ancestors of these people arrived in the subcontinent to create so many different communities.

Y-STR Haplogroup Diversity in the Jat Population Reveals Several Different Ancient Origins

The Jats represent a large ethnic community that has inhabited the northwest region of India and Pakistan for several thousand years. It is estimated the community has a population of over 123 million people. Many historians and academics have asserted that the Jats are descendants of Aryans, Scythians, or other ancient people that arrived and lived in northern India at one time. Essentially, the specific origin of these people has remained a matter of contention for a long time. This study demonstrated that the origins of Jats can be clarified by identifying their Y-chromosome haplogroups and tracing their genetic markers on the Y-DNA haplogroup tree. A sample of 302 Y-chromosome haplotypes of Jats in India and Pakistan was analyzed. The results showed that the sample population had several different lines of ancestry and emerged from at least nine different geographical regions of the world. It also became evident that the Jats did not have a unique set of genes, but shared an underlying genetic unity with several other ethnic communities in the Indian subcontinent. A startling new assessment of the genetic ancient origins of these people was revealed with DNA science.

High Y-chromosomal Differentiation Among Ethnic Groups of Dir and Swat Districts, Pakistan

The ethnic groups that inhabit the mountainous Dir and Swat districts of northern Pakistan are marked by high levels of cultural and phenotypic diversity. To obtain knowledge of the extent of genetic diversity in this region, we investigated Y-chromosomal diversity in five population samples representing the three main ethnic groups residing within these districts, including Gujars, Pashtuns and Kohistanis. A total of 27 Y-chromosomal short tandem repeats (Y-STRs) and 331 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) were investigated. In the Y-STRs, we observed very high and significant levels of genetic differentiation in nine of the 10 pairwise between-group comparisons (R_ST 0.179-0.746), and the differences were mirrored in the Y-SNP haplogroup frequency distribution. No genetic differences were found between the two Pashtun subethnic groups Tarklanis and Yusafzais (R_ST = 0.000). Utmankhels, also considered Pashtuns culturally, were not closely related to any of the other population samples (R_ST 0.451-0.746). Thus, our findings provide examples of both associations and dissociations between cultural and genetic legacies. When analyzed within a larger continental-scale context, these five ethnic groups fall mostly outside the previously characterized Y-chromosomal gene pools of the Indo-Pakistani subcontinent. Male founder effects, coupled with culturally and topographically based constraints upon marriage and movement, are likely responsible for the high degree of genetic structure in this region.

On the Bantu expansion

Here we report the results of fine resolution Y chromosomal analyses (Y-SNP and Y-STR) of 267 Bantu-speaking males from three populations located in the southeast region of Africa. In an effort to determine the relative Y chromosomal affinities of these three genotyped populations, the findings are interpreted in the context of 74 geographically and ethnically targeted African reference populations representing four major ethno-linguistic groups (Afro-Asiatic, Niger Kordofanin, Khoisan and Pygmoid). In this investigation, we detected a general similarity in the Y chromosome lineages among the geographically dispersed Bantu-speaking populations suggesting a shared heritage and the shallow time depth of the Bantu Expansion. Also, micro-variations in the Bantu Y chromosomal composition across the continent highlight location-specific gene flow patterns with non-Bantu-speaking populations (Khoisan, Pygmy, Afro-Asiatic). Our Y chromosomal results also indicate that the three Bantu-speaking Southeast populations genotyped exhibit unique gene flow patterns involving Eurasian populations but fail to reveal a prevailing genetic affinity to East or Central African Bantu-speaking groups. In addition, the Y-SNP data underscores a longitudinal partitioning in sub-Sahara Africa of two R1b1 subgroups, R1b1-P25* (west) and R1b1a2-M269 (east). No evidence was observed linking the B2a haplogroup detected in the genotyped Southeast African Bantu-speaking populations to gene flow from contemporary Khoisan groups.

Y-STR markers from Ladakh in the Himalayas

A total of 223 samples from the general population of Ladakh in Northwestern India were amplified at 17 Y-STR loci using the AmpFlSTR® Yfiler™ system. The DNA profiles generated were employed to generate allelic frequencies, gene diversity, haplotype diversity and discrimination capacity values as well as number of different haplotypes, fraction of unique haplotypes and Rst pair wise genetic distances. Multidimensional Scaling (MDS) and Correspondence Analysis (CA) were performed with the Rst values and allelic frequencies, respectively. The 17-loci discrimination capacity of Ladakh was found to be 0.8093. Eleven out of the 16 loci have diversity values greater than 0.6, and 13 loci possess values greater than 0.5. Ladakh exhibits no significant genetic difference to seven of the 15 reference forensic databases after Bonferroni correction, three of which are located in South Central Asian and four are from the Himalayan region. Rst genetic distance values before and after Bonferroni corrections illustrate the capacity of the Yfiler system to discriminate among Himalayan populations. The intermediate position of the Ladakh population in the MDS and CA plots likely reflects genetic flow and admixture with neighboring populations. In addition, the longitudinal partition of populations in the MDS and CA plots likely reflect human dispersals such as the silk road migrations.

Northward genetic penetration across the Himalayas viewed from Sherpa people

The Himalayas have been suggested as a natural barrier for human migrations, especially the northward dispersals from the Indian Subcontinent to Tibetan Plateau. However, although the majority of Sherpa have a Tibeto-Burman origin, considerable genetic components from Indian Subcontinent have been observed in Sherpa people living in Tibet. The western Y chromosomal haplogroups R1a1a-M17, J-M304, and F*-M89 comprise almost 17% of Sherpa paternal gene pool. In the maternal side, M5c2, M21d, and U from the west also count up to 8% of Sherpa people. Those lineages with South Asian origin indicate that the Himalayas have been permeable to bidirectional gene flow.

Ladakh, India: the land of high passes and genetic heterogeneity reveals a confluence of migrations

Owing to its geographic location near the longitudinal center of Asia, Ladakh, the land of high passes, has witnessed numerous demographic movements during the past millenniums of occupation. In an effort to view Ladakh's multicultural history from a paternal genetic perspective, we performed a high-resolution Y-chromosomal survey of Ladakh, within the context of Y haplogroup and haplotype distributions of 41 Asian reference populations. The results of this investigation highlight the rich ethnic and genetic diversity of Ladkah which includes genetic contributions from disparate regions of the continent including, West, East, South and Central Asia. The phylogenetic signals from Ladakh are consistent with the Indo-Aryans' occupation during the Neolithic age and its historic connection with Tibet, as well as the East-West gene flow associated with the Silk Road.

Zalloua P, Yepiskoposyan L, Dybo A, Tyler-Smith C, Balanovska E. Deep phylogenetic analysis of haplogroup G1 provides estimates of SNP and STR mutation rates on the human Y-chromosome and reveals migrations of Iranic speakers

Y-chromosomal haplogroup G1 is a minor component of the overall gene pool of South-West and Central Asia but reaches up to 80% frequency in some populations scattered within this area. We have genotyped the G1-defining marker M285 in 27 Eurasian populations (n= 5,346), analyzed 367 M285-positive samples using 17 Y-STRs, and sequenced ~11 Mb of the Y-chromosome in 20 of these samples to an average coverage of 67X. This allowed detailed phylogenetic reconstruction. We identified five branches, all with high geographical specificity: G1-L1323 in Kazakhs, the closely related G1-GG1 in Mongols, G1-GG265 in Armenians and its distant brother clade G1-GG162 in Bashkirs, and G1-GG362 in West Indians. The haplotype diversity, which decreased from West Iran to Central Asia, allows us to hypothesize that this rare haplogroup could have been carried by the expansion of Iranic speakers northwards to the Eurasian steppe and via founder effects became a predominant genetic component of some populations, including the Argyn tribe of the Kazakhs. The remarkable agreement between genetic and genealogical trees of Argyns allowed us to calibrate the molecular clock using a historical date (1405 AD) of the most recent common genealogical ancestor. The mutation rate for Y-chromosomal sequence data obtained was 0.78×10^-9 per bp per year, falling within the range of published rates. The mutation rate for Y-chromosomal STRs was 0.0022 per locus per generation, very close to the so-called genealogical rate. The “clan-based” approach to estimating the mutation rate provides a third, middle way between direct farther-to-son comparisons and using archeologically known migrations, whose dates are subject to revision and of uncertain relationship to genetic events.

Simultaneous analysis of hundreds of Y-chromosomal SNPs for high-resolution paternal lineage classification using targeted semiconductor sequencing

SNPs from the non-recombining part of the human Y chromosome (Y-SNPs) are informative to classify paternal lineages in forensic, genealogical, anthropological, and evolutionary studies. Although thousands of Y-SNPs were identified thus far, previous Y-SNP multiplex tools target only dozens of markers simultaneously, thereby restricting the provided Y-haplogroup resolution and limiting their applications. Here, we overcome this shortcoming by introducing a high-resolution multiplex tool for parallel genotyping-by-sequencing of 530 Y-SNPs using the Ion Torrent PGM platform, which allows classification of 432 worldwide Y haplogroups. Contrary to previous Y-SNP multiplex tools, our approach covers branches of the entire Y tree, thereby maximizing the paternal lineage classification obtainable. We used a default DNA input amount of 10 ng per reaction but preliminary sensitivity testing revealed positive results from as little as 100 pg input DNA. Furthermore, we demonstrate that sample pooling using barcodes is feasible, allowing increased throughput for lower per-sample costs. In addition to the wetlab protocol, we provide a software tool for automated data quality control and haplogroup classification. The unique combination of ultra-high marker density and high sensitivity achievable from low amounts of potentially degraded DNA makes this new multiplex tool suitable for a wide range of Y-chromosome applications.

At the southeast fringe of the Bantu expansion: genetic diversity and phylogenetic relationships to other sub-Saharan tribes

Here, we present 12 loci paternal haplotypes (Y-STR profiles) against the backdrop of the Y-SNP marker system of Bantu males from the Maputo Province of Southeast Africa, a region believed to represent the southeastern fringe of the Bantu expansion. Our Maputo Bantu group was analyzed within the context of 27 geographically relevant reference populations in order to ascertain its genetic relationship to other Bantu and non Bantu (Pygmy, Khoisan and Nilotic) sub-equatorial tribes from West and East Africa. This study entails statistical pair wise comparisons and multidimensional scaling based on YSTR Rst distances, network analyses of Bantu (B2a-M150) and Pygmy (B2b-M112) lineages as well as an assessment of Y-SNP distribution patterns. Several notable findings include the following: 1) the Maputo Province Bantu exhibits a relatively close paternal affinity with both east and west Bantu tribes due to high proportion of Bantu Y chromosomal markers, 2) only traces of Khoisan (1.3%) and Pygmy (1.3%) markers persist in the Maputo Province Bantu gene pool, 3) the occurrence of R1a1a-M17/M198, a member of the Eurasian R1a-M420 branch in the population of the Maputo Province, may represent back migration events and/or recent admixture events, 4) the shared presence of E1b1b1-M35 in all Tanzanian tribes examined, including Bantu and non-Bantu groups, in conjunction with its nearly complete absence in the West African populations indicate that, in addition to a shared linguistic, cultural and genetic heritage, geography (e.g., east vs. west) may have impacted the paternal landscape of sub-Saharan Africa, 5) the admixture and assimilation processes of Bantu elements were both highly complex and region-specific.

•

Maputo Bantus exhibit close affinities with other West and East African Bantus.

•

Traces of Khoisan and Pygmy markers persist in the Maputo Province Bantus.

•

R1a1a-M17/M198 in the Maputo Province may represent back or recent migration.

•

Linguistic, cultural and genetic heritages are reflected in Maputo's gene pool.

•

Admixture and assimilation processes of Bantu elements were region-specific.

The phylogenetic and geographic structure of Y-chromosome haplogroup R1a

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25,000 (95% CI: 21,300-29,000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran.