Ancestral Origins and Admixture History of Kazakhs

Kazakh people, like many other populations that settled in Central Asia, demonstrate an array of mixed anthropological features of East Eurasian (EEA) and West Eurasian (WEA) populations, indicating a possible scenario of biological admixture between already differentiated EEA and WEA populations. However, their complex biological origin, genomic makeup, and genetic interaction with surrounding populations are not well understood. To decipher their genetic structure and population history, we conducted, to our knowledge, the first whole-genome sequencing study of Kazakhs residing in Xinjiang (KZK). We demonstrated that KZK derived their ancestries from 4 ancestral source populations: East Asian (∼39.7%), West Asian (∼28.6%), Siberian (∼23.6%), and South Asian (∼8.1%). The recognizable interactions of EEA and WEA ancestries in Kazakhs were dated back to the 15th century BCE. Kazakhs were genetically distinctive from the Uyghurs in terms of their overall genomic makeup, although the 2 populations were closely related in genetics, and both showed a substantial admixture of western and eastern peoples. Notably, we identified a considerable sex-biased admixture, with an excess of western males and eastern females contributing to the KZK gene pool. We further identified a set of genes that showed remarkable differentiation in KZK from the surrounding populations, including those associated with skin color (SLC24A5, OCA2), essential hypertension (HLA-DQB1), hypertension (MTHFR, SLC35F3), and neuron development (CNTNAP2). These results advance our understanding of the complex history of contacts between Western and Eastern Eurasians, especially those living or along the old Silk Road.

Multiple Human Population Movements and Cultural Dispersal Events Shaped the Landscape of Chinese Paternal Heritage

Large-scale genomic projects and ancient DNA innovations have ushered in a new paradigm for exploring human evolutionary history. However, the genetic legacy of spatiotemporally diverse ancient Eurasians within Chinese paternal lineages remains unresolved. Here, we report an integrated Y-chromosome genomic database encompassing 15,563 individuals from both modern and ancient Eurasians, including 919 newly reported individuals, to investigate the Chinese paternal genomic diversity. The high-resolution, time-stamped phylogeny reveals multiple diversification events and extensive expansions in the early and middle Neolithic. We identify four major ancient population movements, each associated with technological innovations that have shaped the Chinese paternal landscape. First, the expansion of early East Asians and millet farmers from the Yellow River Basin predominantly carrying O2/D subclades significantly influenced the formation of the Sino-Tibetan people and facilitated the permanent settlement of the Tibetan Plateau. Second, the dispersal of rice farmers from the Yangtze River Valley carrying O1 and certain O2 sublineages reshapes the genetic makeup of southern Han Chinese, as well as the Tai-Kadai, Austronesian, Hmong-Mien, and Austroasiatic people. Third, the Neolithic Siberian Q/C paternal lineages originated and proliferated among hunter-gatherers on the Mongolian Plateau and the Amur River Basin, leaving a significant imprint on the gene pools of northern China. Fourth, the J/G/R paternal lineages derived from western Eurasia, which were initially spread by Yamnaya-related steppe pastoralists, maintain their presence primarily in northwestern China. Overall, our research provides comprehensive genetic evidence elucidating the significant impact of interactions with culturally distinct ancient Eurasians on the patterns of paternal diversity in modern Chinese populations.

Reconstructing complex admixture history using a hierarchical model

Various methods have been proposed to reconstruct admixture histories by analyzing the length of ancestral chromosomal tracts, such as estimating the admixture time and number of admixture events. However, available methods do not explicitly consider the complex admixture structure, which characterizes the joining and mixing patterns of different ancestral populations during the admixture process, and instead assume a simplified one-by-one sequential admixture model. In this study, we proposed a novel approach that considers the non-sequential admixture structure to reconstruct admixture histories. Specifically, we introduced a hierarchical admixture model that incorporated four ancestral populations and developed a new method, called HierarchyMix, which uses the length of ancestral tracts and the number of ancestry switches along genomes to reconstruct the four-way admixture history. By automatically selecting the optimal admixture model using the Bayesian information criterion principles, HierarchyMix effectively estimates the corresponding admixture parameters. Simulation studies confirmed the effectiveness and robustness of HierarchyMix. We also applied HierarchyMix to Uyghurs and Kazakhs, enabling us to reconstruct the admixture histories of Central Asians. Our results highlight the importance of considering complex admixture structures and demonstrate that HierarchyMix is a useful tool for analyzing complex admixture events.

T2T-YAO: A Telomere-to-telomere Assembled Diploid Reference Genome for Han Chinese

Since its initial release in 2001, the human reference genome has undergone continuous improvement in quality, and the recently released telomere-to-telomere (T2T) version - T2T-CHM13 - reaches its highest level of continuity and accuracy after 20 years of effort by working on a simplified, nearly homozygous genome of a hydatidiform mole cell line. Here, to provide an authentic complete diploid human genome reference for the Han Chinese, the largest population in the world, we assembled the genome of a male Han Chinese individual, T2T-YAO, which includes T2T assemblies of all the 22 + X + M and 22 + Y chromosomes in both haploids. The quality of T2T-YAO is much better than those of all currently available diploid assemblies, and its haploid version, T2T-YAO-hp, generated by selecting the better assembly for each autosome, reaches the top quality of fewer than one error per 29.5 Mb, even higher than that of T2T-CHM13. Derived from an individual living in the aboriginal region of the Han population, T2T-YAO shows clear ancestry and potential genetic continuity from the ancient ancestors. Each haplotype of T2T-YAO possesses ∼ 330-Mb exclusive sequences, ∼ 3100 unique genes, and tens of thousands of nucleotide and structural variations as compared with CHM13, highlighting the necessity of a population-stratified reference genome. The construction of T2T-YAO, an accurate and authentic representative of the Chinese population, would enable precise delineation of genomic variations and advance our understandings in the hereditability of diseases and phenotypes, especially within the context of the unique variations of the Chinese population.

Inferring the population history of Kyrgyz in Xinjiang, Northwest China from genome-wide array genotyping

OBJECTIVES

Xinjiang plays a vital role in the trans-Eurasian population migration, language diffusion, and culture and technology exchange. However, the underrepresentation of Xinjiang's genomes has hindered a more comprehensive understanding of Xinjiang's genetic structure and population history.

MATERIALS AND METHODS

We collected and genotyped 70 southern Xinjiang's Kyrgyz (SXJK) individuals and combined the data with modern and ancient Eurasians published. We used allele-frequency methods, including PCA, ADMIXTURE, f-statistics, qpWave/qpAdm, ALDER, Treemix, and haplotype-shared methods including shared-IBD segments, fineSTRUCTURE, and GLOBETROTTER to unveil the fine-scale population structure and reconstruct admixture history.

RESULTS

We identified genetic substructure within the SXJK population with subgroups showing different genetic affinities to West and East Eurasians. All SXJK subgroups were suggested to have close genetic relationships with surrounding Turkic-speaking groups that is, Uyghur, Kyrgyz from north Xinjiang and Tajikistan, and Chinese Kazakh, suggesting a shared ancestry among those populations. Outgroup-f3 and symmetrical f4 statistics showed a high genetic affinity of SXJK to present-day Tungusic, Mongolic-speaking populations and Ancient Northeast Asian (ANA) related groups. Allele sharing and haplotype sharing profiles revealed the east-west admixture pattern of SXJK. The qpAdm-based admixture models showed that SXJK derived ancestry from East Eurasian (ANA and East Asian, 42.7%-83.3%) and West Eurasian (Western Steppe herders and Central Asian, 16.7%-57.3%), the recent east-west admixture event could be traced to 1000 years ago based on ALDER and GLOBETROTTER analysis.

DISCUSSION

The high genetic affinity of SXJK to present-day Tungusic and Mongolic-speaking populations and short-shared IBD segments indicated their shared common ancestry. SXJK harbored a close genetic affinity to ANA-related populations, indicating the Northeast Asian origin of SXJK. The West and East Eurasian admixture models observed in SXJK further provided evidence of the dynamic admixture history in Xinjiang. The east-west admixture pattern and the identified ancestral makeup of SXJK suggested a genetic continuity from some Iron Age Xinjiang populations to present-day SXJK.

Expression profiles of east-west highly differentiated genes in Uyghur genomes

It remains unknown and debatable how European-Asian-differentiated alleles affect individual phenotypes. Here, we made the first effort to analyze the expression profiles of highly differentiated genes with eastern and western origins in 90 Uyghurs using whole-genome (30× to 60×) and transcriptome data. We screened 921 872 east-west highly differentiated genetic variants, of which ∼4.32% were expression quantitative trait loci (eQTLs), ∼0.12% were alternative splicing quantitative trait loci (sQTLs), and ∼0.12% showed allele-specific expression (ASE). The 8305 highly differentiated eQTLs of strong effects appear to have undergone natural selection, associated with immunity and metabolism. European-origin alleles tend to be more biasedly expressed; highly differentiated ASEs were enriched in diabetes-associated genes, likely affecting the diabetes susceptibility in the Uyghurs. We proposed an admixture-induced expression model to dissect the highly differentiated expression profiles. We provide new insights into the genetic basis of phenotypic differentiation between Western and Eastern populations, advancing our understanding of the impact of genetic admixture.

Genomic insights into the genetic structure and population history of Mongolians in Liaoning Province

The Mongolian population exceeds six million and is the largest population among the Mongolic speakers in China. However, the genetic structure and admixture history of the Mongolians are still unclear due to the limited number of samples and lower coverage of single-nucleotide polymorphism (SNP). In this study, we genotyped genome-wide data of over 700,000 SNPs in 38 Mongolian individuals from Fuxin in Liaoning Province to explore the genetic structure and population history based on typical and advanced population genetic analysis methods [principal component analysis (PCA), admixture, F_ST, f₃-statistics, f₄-statistics, qpAdm/qpWave, qpGraph, ALDER, and TreeMix]. We found that Fuxin Mongolians had a close genetic relationship with Han people, northern Mongolians, other Mongolic speakers, and Tungusic speakers in East Asia. Also, we found that Neolithic millet farmers in the Yellow River Basin and West Liao River Basin and Neolithic hunter–gatherers in the Mongolian Plateau and Amur River Basin were the dominant ancestral sources, and there were additional gene flows related to Eurasian Steppe pastoralists and Neolithic Iranian farmers in the gene pool of Fuxin Mongolians. These results shed light on dynamic demographic history, complex population admixture, and multiple sources of genetic diversity in Fuxin Mongolians.

MultiWaverX: modeling latent sex-biased admixture history

Sex-biased gene flow has been common in the demographic history of modern humans. However, the lack of sophisticated methods for delineating the detailed sex-biased admixture process prevents insights into complex admixture history and thus our understanding of the evolutionary mechanisms of genetic diversity. Here, we present a novel algorithm, MultiWaverX, for modeling complex admixture history with sex-biased gene flow. Systematic simulations showed that MultiWaverX is a powerful tool for modeling complex admixture history and inferring sex-biased gene flow. Application of MultiWaverX to empirical data of 17 typical admixed populations in America, Central Asia, and the Middle East revealed sex-biased admixture histories that were largely consistent with the historical records. Notably, fine-scale admixture process reconstruction enabled us to recognize latent sex-biased gene flow in certain populations that would likely be overlooked by much of the routine analysis with commonly used methods. An outstanding example in the real world is the Kazakh population that experienced complex admixture with sex-biased gene flow but in which the overall signature has been canceled due to biased gene flow from an opposite direction.

The genetic echo of the Tarim mummies in modern Central Asians

The diversity of Central Asians has been shaped by multiple migrations and cultural diffusion. Although ancient DNA studies have revealed the demographic changes of the Central Asian since the Bronze Age, the contribution of the ancient populations to the modern Central Asian remains opaque. Herein, we performed high-coverage sequencing of 131 whole genomes of Indo-European-speaking Tajik and Turkic-speaking Kyrgyz populations to explore their genomic diversity and admixture history. By integrating the ancient DNA data, we revealed more details of the origins and admixture history of Central Asians. We found that the major ancestry of present-day Tajik populations can be traced back to the admixture of the Bronze Age Bactria–Margiana Archaeological Complex and Andronovo-related populations. Highland Tajik populations further received additional gene flow from the Tarim mummies, an isolated ancient North Eurasian–related population. The West Eurasian ancestry of Kyrgyz is mainly derived from Historical Era populations in Xinjiang of China. Furthermore, the recent admixture signals detected in both Tajik and Kyrgyz are ascribed to the expansions of Eastern Steppe nomadic pastoralists during the Historical Era.

Bronze and Iron Age population movements underlie Xinjiang population history

The Xinjiang region in northwest China is a historically important geographical passage between East and West Eurasia. By sequencing 201 ancient genomes from 39 archaeological sites, we clarify the complex demographic history of this region. Bronze Age Xinjiang populations are characterized by four major ancestries related to Early Bronze Age cultures from the central and eastern Steppe, Central Asian, and Tarim Basin regions. Admixtures between Middle and Late Bronze Age Steppe cultures continued during the Late Bronze and Iron Ages, along with an inflow of East and Central Asian ancestry. Historical era populations show similar admixed and diverse ancestries as those of present-day Xinjiang populations. These results document the influence that East and West Eurasian populations have had over time in the different regions of Xinjiang.

Genomic diversity and post-admixture adaptation in the Uyghurs

Population admixture results in genome-wide combinations of genetic variants derived from different ancestral populations of distinct ancestry, thus providing a unique opportunity for understanding the genetic determinants of phenotypic variation in humans. Here, we used whole-genome sequencing of 92 individuals with high coverage (30–60×) to systematically investigate genomic diversity in the Uyghurs living in Xinjiang, China (XJU), an admixed population of both European-like and East-Asian-like ancestry. The XJU population shows greater genetic diversity, especially a higher proportion of rare variants, compared with their ancestral source populations, corresponding to greater phenotypic diversity of XJU. Admixture-induced functional variants in EDAR were associated with the diversity of facial morphology in XJU. Interestingly, the interaction of functional variants between SLC24A5 and OCA2 likely influences the diversity of skin pigmentation. Notably, selection has seemingly been relaxed or canceled in several genes with significantly biased ancestry, such as HERC2–OCA2. Moreover, signatures of post-admixture adaptation in XJU were identified, including genes related to metabolism (e.g. CYP2D6), digestion (e.g. COL11A1), olfactory perception (e.g. ANO2) and immunity (e.g. HLA). Our results demonstrated population admixture as a driving force, locally or globally, in shaping human genetic and phenotypic diversity as well as in adaptive evolution.

The Genetic Structure and East-West Population Admixture in Northwest China Inferred From Genome-Wide Array Genotyping

Northwest China is a contacting region for East and West Eurasia and an important center for investigating the migration and admixture history of human populations. However, the comprehensive genetic structure and admixture history of the Altaic speaking populations and Hui group in Northwest China were still not fully characterized due to insufficient sampling and the lack of genome-wide data. Thus, We genotyped genome-wide SNPs for 140 individuals from five Chinese Mongolic, Turkic speaking groups including Dongxiang, Bonan, Yugur, and Salar, as well as the Hui group. Analysis based on allele-sharing and haplotype-sharing were used to elucidate the population history of Northwest Chinese populations, including PCA, ADMIXTURE, pairwise Fst genetic distance, f-statistics, qpWave/qpAdm and ALDER, fineSTRUCTURE and GLOBETROTTER. We observed Dongxiang, Bonan, Yugur, Salar, and Hui people were admixed populations deriving ancestry from both East and West Eurasians, with the proportions of West Eurasian related contributions ranging from 9 to 15%. The genetic admixture was probably driven by male-biased migration- showing a higher frequency of West Eurasian related Y chromosomal lineages than that of mtDNA detected in Northwest China. ALDER-based admixture and haplotype-based GLOBETROTTER showed this observed West Eurasian admixture signal was introduced into East Eurasia approximately 700 ∼1,000 years ago. Generally, our findings provided supporting evidence that the flourish transcontinental communication between East and West Eurasia played a vital role in the genetic formation of northwest Chinese populations.

The Heart of Silk Road "Xinjiang," Its Genetic Portray, and Forensic Parameters Inferred From Autosomal STRs

The Xinjiang Uyghur Autonomous Region of China (XUARC) harbors almost 50 ethnic groups including the Uyghur (UGR: 45.84%), Han (HAN: 40.48%), Kazakh (KZK: 6.50%), Hui (HUI: 4.51%), Kyrgyz (KGZ: 0.86%), Mongol (MGL: 0.81%), Manchu (MCH: 0.11%), and Uzbek (UZK: 0.066%), which make it one of the most colorful regions with abundant cultural and genetic diversities. In our previous study, we established allelic frequency databases for 14 autosomal short tandem repeats (STRs) for four minority populations from XUARC (MCH, KGZ, MGL, and UZK) using the AmpFlSTR® Identifiler PCR Amplification Kit. In this study, we genotyped 2,121 samples using the GoldenEye™ 20A Kit (Beijing PeopleSpot Inc., Beijing, China) amplifying 19 autosomal STR loci for four major ethnic groups (UGR, HAN, KZK, and HUI). These groups make up 97.33% of the total XUARC population. The total number of alleles for all the 19 STRs in these populations ranged from 232 (HAN) to 224 (KZK). We did not observe any departures from the Hardy-Weinberg equilibrium (HWE) in these populations after sequential Bonferroni correction. We did find minimal departure from linkage equilibrium (LE) for a small number of pairwise combinations of loci. The match probabilities for the different populations ranged from 1 in 1.66 × 1023 (HAN) to 6.05 × 1024 (HUI), the combined power of exclusion ranged from 0.999 999 988 (HUI) to 0.999 999 993 (UGR), and the combined power of discrimination ranged from 0.999 999 999 999 999 999 999 983 (HAN) to 0.999 999 999 999 999 999 999 997 (UGR). Genetic distances, principal component analysis (PCA), STRUCTURE analysis, and the phylogenetic tree showed that genetic affinity among studied populations is consistent with linguistic, ethnic, and geographical classifications.

Micro-coevolution of host genetics with gut microbiome in three Chinese ethnic groups

Understanding the micro-coevolution of the human gut microbiome with host genetics is challenging but essential in both evolutionary and medical studies. To gain insight into the interactions between host genetic variation and the gut microbiome, we analyzed both the human genome and gut microbiome collected from a cohort of 190 students in the same boarding college and representing 3 ethnic groups, Uyghur, Kazakh, and Han Chinese. We found that differences in gut microbiome were greater between genetically distinct ethnic groups than those genetically closely related ones in taxonomic composition, functional composition, enterotype stratification, and microbiome genetic differentiation. We also observed considerable correlations between host genetic variants and the abundance of a subset of gut microbial species. Notably, interactions between gut microbiome species and host genetic variants might have coordinated effects on specific human phenotypes. Bacteroides ovatus, previously reported to modulate intestinal immunity, is significantly correlated with the host genetic variant rs12899811 (meta-P = 5.55 × 10^-5), which regulates the VPS33B expression in the colon, acting as a tumor suppressor of colorectal cancer. These results advance our understanding of the micro-coevolution of the human gut microbiome and their interactive effects with host genetic variation on phenotypic diversity.

The genomic origins of the Bronze Age Tarim Basin mummies

The identity of the earliest inhabitants of Xinjiang, in the heart of Inner Asia, and the languages that they spoke have long been debated and remain contentious¹. Here we present genomic data from 5 individuals dating to around 3000-2800 BC from the Dzungarian Basin and 13 individuals dating to around 2100-1700 BC from the Tarim Basin, representing the earliest yet discovered human remains from North and South Xinjiang, respectively. We find that the Early Bronze Age Dzungarian individuals exhibit a predominantly Afanasievo ancestry with an additional local contribution, and the Early-Middle Bronze Age Tarim individuals contain only a local ancestry. The Tarim individuals from the site of Xiaohe further exhibit strong evidence of milk proteins in their dental calculus, indicating a reliance on dairy pastoralism at the site since its founding. Our results do not support previous hypotheses for the origin of the Tarim mummies, who were argued to be Proto-Tocharian-speaking pastoralists descended from the Afanasievo^1,2 or to have originated among the Bactria-Margiana Archaeological Complex³ or Inner Asian Mountain Corridor cultures⁴. Instead, although Tocharian may have been plausibly introduced to the Dzungarian Basin by Afanasievo migrants during the Early Bronze Age, we find that the earliest Tarim Basin cultures appear to have arisen from a genetically isolated local population that adopted neighbouring pastoralist and agriculturalist practices, which allowed them to settle and thrive along the shifting riverine oases of the Taklamakan Desert.

AdmixSim 2: a forward-time simulator for modeling complex population admixture

BACKGROUND

Computer simulations have been widely applied in population genetics and evolutionary studies. A great deal of effort has been made over the past two decades in developing simulation tools. However, there are not many simulation tools suitable for studying population admixture.

RESULTS

We here developed a forward-time simulator, AdmixSim 2, an individual-based tool that can flexibly and efficiently simulate population genomics data under complex evolutionary scenarios. Unlike its previous version, AdmixSim 2 is based on the extended Wright-Fisher model, and it implements many common evolutionary parameters to involve gene flow, natural selection, recombination, and mutation, which allow users to freely design and simulate any complex scenario involving population admixture. AdmixSim 2 can be used to simulate data of dioecious or monoecious populations, autosomes, or sex chromosomes. To our best knowledge, there are no similar tools available for the purpose of simulation of complex population admixture. Using empirical or previously simulated genomic data as input, AdmixSim 2 provides phased haplotype data for the convenience of further admixture-related analyses such as local ancestry inference, association studies, and other applications. We here evaluate the performance of AdmixSim 2 based on simulated data and validated functions via comparative analysis of simulated data and empirical data of African American, Mexican, and Uyghur populations.

CONCLUSIONS

AdmixSim 2 is a flexible simulation tool expected to facilitate the study of complex population admixture in various situations.

Ancient Mitochondrial Genomes Reveal Extensive Genetic Influence of the Steppe Pastoralists in Western Xinjiang

The population prehistory of Xinjiang has been a hot topic among geneticists, linguists, and archaeologists. Current ancient DNA studies in Xinjiang exclusively suggest an admixture model for the populations in Xinjiang since the early Bronze Age. However, almost all of these studies focused on the northern and eastern parts of Xinjiang; the prehistoric demographic processes that occurred in western Xinjiang have been seldomly reported. By analyzing complete mitochondrial sequences from the Xiabandi (XBD) cemetery (3,500–3,300 BP), the up-to-date earliest cemetery excavated in western Xinjiang, we show that all the XBD mitochondrial sequences fall within two different West Eurasian mitochondrial DNA (mtDNA) pools, indicating that the migrants into western Xinjiang from west Eurasians were a consequence of the early expansion of the middle and late Bronze Age steppe pastoralists (Steppe_MLBA), admixed with the indigenous populations from Central Asia. Our study provides genetic links for an early existence of the Indo-Iranian language in southwestern Xinjiang and suggests that the existence of Andronovo culture in western Xinjiang involved not only the dispersal of ideas but also population movement.

Distinctive Gut Microbiota in Patients with Overweight and Obesity with Dyslipidemia and its Responses to Long-term Orlistat and Ezetimibe Intervention: A Randomized Controlled Open-label Trial

This study investigated the gut microbiota and short chain fatty acids (SCFAs) characteristics of subjects with obesity from Xinjiang in northwestern China, a region with a multiethnic culture and characteristic lifestyle, and to explore the potential microbes that respond to a 12-wk medication of orlistat and ezetimibe with a randomized controlled open-label trial manner. The gut microbiota profile of patients with overweight and obesity with dyslipidemia in Xinjiang was distinctive and characterized by enrichment of Lactobacillus and the reduction of the diversity and the depletion of Actinobacteria, Bacteroides, Bifidobacterium, and Bacteroides fragilis. Prevotella-type, Gemmiger-type, and Escherichia/Shigella-type were the gut microbial patterns of the Xinjiang population. However, the fecal SCFAs levels and enterotypes were similar between healthy individuals and patients. These results indicated that the contribution of the gut microbiota to obesity was highly dependent on geography and dietary habits. Waist circumference, total triglyceride (TG), and fasting blood glucose (FBG) were significantly decreased after orlistat therapy, whereas TG, total cholesterol (TC), and low density lipoprotein cholesterol (LDL-C) were significantly decreased by ezetimibe. Overall, the gut microbiota and their SCFAs metabolites were relatively stable after treatment with the two drugs, with alteration of some low-abundant bacteria, i.e., significantly increased Proteobacteria and decreased Alloprevotella after orlistat, and increased Fusobacteria and Fusobacterium after ezetimibe therapy. These results indicated that intestinal malabsorption of dietary fat and cholesterol caused by orlistat and ezetimibe had a limited effect on the overall gut microbial community and their metabolites. Nevertheless, significant correlations between several core microbes that responded to the medications and biochemical data were found; in particular, Actinomyces and Bacteroides were positively correlated with FBG after orlistat intervention, while Clostridium XVIII and Lachnospiracea incertae sedis were negatively correlated with TC and LDL-C after ezetimibe intervention, thus indicating their roles in improving glucolipid metabolism in obesity by acting as potential microbial targets.

A 3,000-year-old, basal S. enterica lineage from Bronze Age Xinjiang suggests spread along the Proto-Silk Road

Salmonella enterica (S. enterica) has infected humans for a long time, but its evolutionary history and geographic spread across Eurasia is still poorly understood. Here, we screened for pathogen DNA in 14 ancient individuals from the Bronze Age Quanergou cemetery (XBQ), Xinjiang, China. In 6 individuals we detected S. enterica. We reconstructed S. enterica genomes from those individuals, which form a previously undetected phylogenetic branch basal to Paratyphi C, Typhisuis and Choleraesuis-the so-called Para C lineage. Based on pseudogene frequency, our analysis suggests that the ancient S. enterica strains were not host adapted. One genome, however, harbors the Salmonella pathogenicity island 7 (SPI-7), which is thought to be involved in (para)typhoid disease in humans. This offers first evidence that SPI-7 was acquired prior to the emergence of human-adapted Paratyphi C around 1,000 years ago. Altogether, our results show that Salmonella enterica infected humans in Eastern Eurasia at least 3,000 years ago, and provide the first ancient DNA evidence for the spread of a pathogen along the Proto-Silk Road.

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 origins and sex-biased admixture of the Huis

The Hui people are unique among Chinese ethnic minorities in that they speak the same language as Han Chinese (HAN) but practice Islam. However, as the second-largest minority group in China numbering well over 10 million, the Huis are under-represented in both global and regional genomic studies. Here, we present the first whole-genome sequencing effort of 234 Hui individuals (NXH) aged over 60 who have been living in Ningxia, where the Huis are mostly concentrated. NXH are genetically more similar to East Asian than to any other global populations. In particular, the genetic differentiation between NXH and HAN (FST = 0.0015) is only slightly larger than that between northern and southern HAN (FST = 0.0010), largely attributed to the western ancestry in NXH (∼10%). Highly differentiated functional variants between NXH and HAN were identified in genes associated with skin pigmentation (e.g., SLC24A5), facial morphology (e.g., EDAR), and lipid metabolism (e.g., ABCG8). The Huis are also distinct from other Muslim groups such as the Uyghurs (FST = 0.0187), especially, NXH derived much less western ancestry (∼10%) compared with the Uyghurs (∼50%). Modeling admixture history indicated that NXH experienced an episode of two-wave admixture. An ancient admixture occurred ∼1,025 years ago, reflecting the intensive west-east contacts during the late Tang Dynasty, and the Five Dynasties and Ten Kingdoms period. A recent admixture occurred ∼500 years ago, corresponding to the Ming Dynasty. Notably, we identified considerable sex-biased admixture, i.e., excess of western males and eastern females contributing to the NXH gene pool. The origins and the genomic diversity of the Hui people imply the complex history of contacts between western and eastern Eurasians.

Genetic analysis of 32 InDels in four ethnic minorities from Chinese Xinjiang

The present study used the previously constructed 32-plex InDels panel to investigated the genetic diversity of four ethnic minorities (Hui, Mongol, Uygur and Kazakh) from Xinjiang, and analyzed the genetic relationships between the four populations and 27 reference populations. No significant deviations were observed from the Hardy-Weinberg equilibrium (HWE) at the 32 InDels for each population. The average observed heterozygosity (H_exp), average polymorphic information content (PIC), combined power of discrimination (CPD) and cumulative probability of exclusion (CPE) for the 32 InDels were all higher than the Qiagen Investigator DIPplex kit in the four populations from Xinjiang. The CPD ranged from 0.999999999999903 (Kazakh) to 0.999999999999952 (Hui) and CPE ranged from 0.9971 (Uygur) to 0.9985 (Hui), which indicated that the 32 InDels were capable for individual identification and could be a supplementary tool in paternity test for these populations. Population genetic analysis by the method of analysis of molecular variance (AMOVA), F_ST, phylogenetic tree, TreeMix-based topology, multi-dimensional scale analysis (MDS), principal components analysis (PCA) and STRUCTURE analysis showed that Xinjiang Hui population has a close relationship with East Asians (EAS), especially Chinese Han, and the populations of Xinjiang Mongol, Uygur and Kazakh showed mixed ancestral components related to EAS and Europeans (EUR).

Ancient Xinjiang mitogenomes reveal intense admixture with high genetic diversity

Xinjiang is a key region in northwestern China, connecting East and West Eurasian populations and cultures for thousands of years. To understand the genetic history of Xinjiang, we sequenced 237 complete ancient human mitochondrial genomes from the Bronze Age through Historical Era (41 archaeological sites). Overall, the Bronze Age Xinjiang populations show high diversity and regional genetic affinities with Steppe and northeastern Asian populations along with a deep ancient Siberian connection for the Tarim Basin Xiaohe individuals. In the Iron Age, in general, Steppe-related and northeastern Asian admixture intensified, with North and East Xinjiang populations showing more affinity with northeastern Asians and South Xinjiang populations showing more affinity with Central Asians. The genetic structure observed in the Historical Era of Xinjiang is similar to that in the Iron Age, demonstrating genetic continuity since the Iron Age with some additional genetic admixture with populations surrounding the Xinjiang region.

Massively parallel sequencing of 165 ancestry-informative SNPs and forensic biogeographical ancestry inference in three southern Chinese Sinitic/Tai-Kadai populations

Ancestry informative markers (AIMs), which are distributed throughout the human genome, harbor significant allele frequency differences among diverse ethnic groups. The use of sets of AIMs to reconstruct population history and genetic relationships is attracting interest in the forensic community, because biogeographic ancestry information for a casework sample can potentially be predicted and used to guide the investigative process. However, subpopulation ancestry inference within East Asia remains in its infancy due to a lack of population reference data collection and incomplete validation work on newly developed or commercial AIM sets. In the present study, 316 Chinese persons, including 85 Sinitic-speaking Haikou Han, 120 Qiongzhong Hlai and 111 Daozhen Gelao individuals belonging to Tai-Kadai-speaking populations, were analyzed using the Precision ID Ancestry Panel (165 AISNPs). Combined with our previous 165-AISNP data (375 individuals from 6 populations), the 1000 Genomes Project and forensic literature, comprehensive population genetic comparisons and ancestry inference were further performed via ADMIXTURE, TreeMix, PCA, f-statistics and N-J tree. Although several nonpolymorphic loci were identified in the three southern Chinese populations, the forensic parameters of this ancestry inference panel were better than those for the 23 STR-based Huaxia Platinum System, which is suitable for use as a robust tool in forensic individual identification and parentage testing. The results based on the ancestry assignment and admixture proportion evaluation revealed that this panel could be used successfully to assign individuals at a continental scale but also possessed obvious limitations in discriminatory power in intercontinental individuals, especially for European-Asian admixed Uyghurs or in populations lacking reference databases. Population genetic analyses further revealed five continental population clusters and three East Asian-focused population subgroups, which is consistent with linguistic affiliations. Ancestry composition and multiple phylogenetic analysis further demonstrated that the geographically isolated Qiongzhong Hlai harbored a close phylogenetic relationship with Austronesian speakers and possessed a homogenous Tai-Kadai-dominant ancestry, which could be used as the ancestral source proxy in population history reconstruction of Tai-Kadai-speaking populations and as one of the representatives for forensic database establishment. In summary, more population-specific AIM sets focused on East Asian subpopulations, comprehensive algorithms and high-coverage population reference data should be developed and validated in the next step.

Population genomics of East Asian ethnic groups

East Asia constitutes one-fifth of the global population and exhibits substantial genetic diversity. However, genetic investigations on populations in this region have been largely under-represented compared with European populations. Nonetheless, the last decade has seen considerable efforts and progress in genome-wide genotyping and whole-genome sequencing of the East-Asian ethnic groups. Here, we review the recent studies in terms of ancestral origin, population relationship, genetic differentiation, and admixture of major East- Asian groups, such as the Chinese, Korean, and Japanese populations. We mainly focus on insights from the whole-genome sequence data and also include the recent progress based on mitochondrial DNA (mtDNA) and Y chromosome data. We further discuss the evolutionary forces driving genetic diversity in East-Asian populations, and provide our perspectives for future directions on population genetics studies, particularly on underrepresented indigenous groups in East Asia.

AdmixSim: A Forward-Time Simulator for Various Complex Scenarios of Population Admixture

Background: Population admixture is a common phenomenon in humans, animals, and plants, and it plays a very important role in shaping individual genetic architecture and population genetic diversity. Inference of population admixture, however, is very challenging and typically relies on in silico simulation. We are aware of the lack of a computerized tool for such a purpose. A simulator capable of generating data under various complex admixture scenarios would facilitate the study of recombination, linkage disequilibrium, ancestry tracing, and admixture dynamics in admixed populations. We described such a simulator here.

Results: We developed a forward-time simulator (AdmixSim) under the standard Wright Fisher model. It can simulate the following admixed populations: (1) multiple ancestral populations; (2) multiple waves of admixture events; (3) fluctuating population size; and (4) admixtures of fluctuating proportions. Analysis of the simulated data by AdmixSim showed that our simulator can quickly and accurately generate data resembling real-world values. We included in AdmixSim all possible parameters that would allow users to modify and simulate any kind of admixture scenario easily, so it is very flexible. AdmixSim records recombination break points and traces of each chromosomal segment from different ancestral populations, with which users can easily perform further analysis and comparative studies with empirical data.

Conclusions:AdmixSim facilitates the study of population admixture by providing a simulation framework with the flexible implementation of various admixture models and parameters.

Forensic Features and Population Genetic Structure of Dong, Yi, Han, and Chuanqing Human Populations in Southwest China Inferred From Insertion/Deletion Markers

Guizhou province in southwest China has abundant genetic and cultural diversities, but the forensic features and genetic structure of Guizhou populations remain poorly understood due to the sparse sampling of present-day populations. Here, we present 30 insertion/deletion polymorphisms (InDels) data of 591 human individuals collected from four populations, Dong, Yi, Han, and Chuanqing residing in Guizhou. We calculated the forensic parameters of 30 InDel loci and found that this panel meets the efficiency of forensic personal identification based on the high combined power of discrimination, but it could only be used as a complementary tool in the parentage testing because of the lower combined probability of exclusion values. The studied populations are genetically closer related to geographically adjacent or linguistically related populations in southern China, such as the Tai-Kadai and Hmong-Mien speaking groups. The unrecognized ethnic Chuanqing people show an additional genetic affinity with Han Chinese, highlighting the role of possible military immigrations in their origin.

Genetic structure and forensic characteristics of the Kyrgyz population from Kizilsu Kirghiz autonomous prefecture based on autosomal DIPs

Living in the heart of Eurasia, the Kyrgyz ethnic minority have a complex human evolutionary and migration history. However, the genetic architecture of the Kyrgyz population has not been fully explored. We studied 526 Kyrgyz samples from Kizilsu Kirghiz Autonomous Prefecture in Xinjiang using the Investigator® DIPplex kit. All loci followed Hardy-Weinberg equilibrium (HWE). The combined power of discrimination (CPD) and combined power of paternity exclusion (CPE) was 0.9999999999988 and 0.9936, respectively. Compared with 90 reference populations, five InDels (HLD99, HLD81, HLD64, HLD118, and HLD111) have the potential to distinguish the Kyrgyz/Uyghur/Kazak population from other East Asian populations. Our results suggested a close genetic relationship between the Kyrgyz population and the Uyghur/Kazak populations, followed by South Asian populations. This was in accordance with the inland migration hypothesis or modern human migration influenced by warfare. Overall, this system can be used as a powerful tool in forensic individual identification and as a complementary tool in paternity cases and biogeographic ancestry analyses.

Forensic genetic polymorphisms and population structure of the Guizhou Bouyei people based on 19 X-STR loci

Background: Guizhou province is located in southwest China with abundant genetic, linguistic and cultural diversity. The Bouyei is one of the 18 officially recognised minority groups in Guizhou, accounting for about 97% of the total Bouyei population in China. However, the genetic history and forensic characterisation of the Bouyei people is largely unknown due to a lack of genetic data.Aim: We aim to investigate genetic polymorphisms and forensic characterisation of the Guizhou Bouyei population, as well as the relationships between the Bouyei and other East Asian populations.Subjects and methods: We genotyped 19 X-STRs in 188 males and 165 females of Guizhou Bouyei using the AGCU X19 STR Kit. We estimated allele frequencies, forensic parameters and genetic distances between the Bouyei and other East Asian populations. We presented the genetic distances in a phylogenetic tree, an MDS plot and a PCA plot.Results: In Guizhou Bouyei individuals, we observed 216 alleles with corresponding frequencies ranging from 0.0019 to 0.6757. All of the six combined powers of PD_m, PD_f, MEC Krüger, MEC Kishida, MEC Desmarais and MEC Desmarais in allele diversity and haplotype diversity are larger than 0.99999995. We found genetic affinities among the Bouyei people and their geographical neighbouring populations in Guizhou, such as the Sui, Miao and Han.Conclusions: The highly polymorphic and informative forensic parameters of the 19 X-STRs in Bouyei people show the powerful potential of those markers in forensic identification and parentage tests. The genetic relationships of the Bouyei with other East Asian populations correspond well with geographic affiliations as well as linguistic classifications.

Forensic features and genetic structure of the Hotan Uyghur inferred from 27 forensic markers

Background: The Uyghur is the fifth largest group in China and also the largest ethnic minority in Xinjiang.Aim: To explore the genetic variations of 27 forensic genetic markers included in the newly developed SureID^® PanGlobal Human DNA Identification System and analyse the genetic relationship between Xinjiang Uyghur and their neighbours.Subjects and methods: We genotyped 27 markers in 2,189 unrelated Uyghur individuals from the Hotan Prefecture in Southwest Xinjiang. Comprehensive population genetic studies among Chinese populations and worldwide populations were conducted via various statistics.Results: The combined power of discrimination (CPD) and the combined power of exclusion (CPE) of the new-generation autosomal STR amplification system in the Hotan Uyghur are 9,9999-E01 and 9,9999-E01, respectively. Population genetic studies indicate that the Hotan Uyghur show a close genetic relationship with geographically different Uyghurs and Kazakhs, while significant genetic differentiation exists between the Hotan Uyghur and some ethnicities from other non-Turkic-speaking populations. The results of population comparisons among the 52 worldwide populations demonstrate that geographically approached intercontinental populations have close genetic relationships.Conclusions: 24 autosomal STRs are highly polymorphic and informative in the Uyghur and this system is suitable for forensic personal identification and paternity testing. Our findings not only reveal that Chinese Uyghur is a homogenous population based on forensic genetic markers, but also indicate that population genetic affinity is closely related to the adjacent populations with common ethnic origin.

MtDNA polymorphism analyses in the Chinese Mongolian group: Efficiency evaluation and further matrilineal genetic structure exploration

BACKGROUND

Profiling of mitochondrial DNA is surely to provide valuable investigative clues for forensic cases involving highly degraded specimens or complex maternal lineage kinship determination. But traditionally used hypervariable region sequencing of mitochondrial DNA is less frequently suggested by the forensic community for insufficient informativeness. Genome-wide sequencing of mitochondrial DNA can provide considerable amount of variant information but can be high cost at the same time.

METHODS

Efficiency of the 60 mitochondrial DNA polymorphic sites dispersing across the control region and coding region of mitochondrial DNA genome was evaluated with 106 Mongolians recruited from the Xinjiang Uyghur Autonomous Region, China, and allele-specific PCR technique was employed for mitochondrial DNA typing.

RESULTS

Altogether 58 haplotypes were observed and the haplotypic diversity, discrimination power and random match probability were calculated to be 0.981, 0.972, and 0.028, respectively. Mitochondrial DNA haplogroup affiliation exhibited an exceeding percentage (12.26%) of west Eurasian lineage (H haplogroup) in the studied Mongolian group, which needed to be further verified with more samples. Furthermore, the genetic relationships between the Xinjiang Mongolian group and the comparison populations were also investigated and the genetic affinity was discovered between the Xinjiang Mongolian group and the Xinjiang Kazak group in this study.

CONCLUSION

It was indicated that the panel was potentially enough to be used as a supplementary tool for forensic applications. And the matrilineal genetic structure analyses based on mitochondrial DNA variants in the Xinjiang Mongolian group could be helpful for subsequent anthropological studies.

PGG.Population: a database for understanding the genomic diversity and genetic ancestry of human populations

There are a growing number of studies focusing on delineating genetic variations that are associated with complex human traits and diseases due to recent advances in next-generation sequencing technologies. However, identifying and prioritizing disease-associated causal variants relies on understanding the distribution of genetic variations within and among populations. The PGG.Population database documents 7122 genomes representing 356 global populations from 107 countries and provides essential information for researchers to understand human genomic diversity and genetic ancestry. These data and information can facilitate the design of research studies and the interpretation of results of both evolutionary and medical studies involving human populations. The database is carefully maintained and constantly updated when new data are available. We included miscellaneous functions and a user-friendly graphical interface for visualization of genomic diversity, population relationships (genetic affinity), ancestral makeup, footprints of natural selection, and population history etc. Moreover, PGG.Population provides a useful feature for users to analyze data and visualize results in a dynamic style via online illustration. The long-term ambition of the PGG.Population, together with the joint efforts from other researchers who contribute their data to our database, is to create a comprehensive depository of geographic and ethnic variation of human genome, as well as a platform bringing influence on future practitioners of medicine and clinical investigators. PGG.Population is available at https://www.pggpopulation.org.

Common variants in SATB2 are associated with schizophrenia in Uygur Chinese population

OBJECTIVE

Schizophrenia is one of the most severe mental disorders and its etiology is supposed to be an interaction between genes and environmental factors. Previous genome-wide association studies of schizophrenia have reported multiple susceptibility loci including rs6704641 in the SATB2 gene. Recently, this locus was further confirmed as a genome-wide significant locus for association with schizophrenia by trans-ancestry meta-analysis of Han Chinese and Caucasian samples. However, there is no report of genetic analysis in Uygur Chinese population, which is considered to have a combined genetic background between eastern Asia and Caucasian. This study is aimed to explore whether SATB2 gene is significantly associated with schizophrenia in Uygur Chinese population, thus providing additional evidence for elucidating the role of SATB2 gene in schizophrenia.

PARTICIPANTS AND METHODS

In this study, we performed a case-control analysis focusing on seven tag single nucleotide polymorphisms located in SATB2 gene among 985 patients with schizophrenia and 1218 healthy controls recruited from the Xinjiang Province of China.

RESULTS

We found that rs6704641 was significantly associated with schizophrenia in both allelic and genotypic distributions (Pallele = 0.008, Pgenotype = 0.028 after correction). In addition, rs16831466 is significantly associated with schizophrenia in allelic distributions (corrected Pallele = 0.041). Besides, several haplotypes of single nucleotide polymorphism are significantly associated with schizophrenia too.

CONCLUSION

Our results suggest that SATB2 is also a susceptibility gene for schizophrenia in Uygur Chinese population, and subsequent functional experiments are necessary to reveal its role in the pathogenesis.

Tai-Kadai-speaking Gelao population: Forensic features, genetic diversity and population structure

Genetic analyses of geographically and ethno-linguistically different populations are essential for understanding population stratification and genomic structure in medical Genome-Wide Association Studies (GWAS) and genetic variation and diversity related to forensic and population genetics studies. Here, we genotyped 30 autosomal insertion/deletion (Indel) markers from 502 Tai-Kadai-speaking Gelao individuals residing in the rugged topographical area in Southeastern China. In addition, two comprehensive population genetic comparisons of 15,327 individuals from 95 worldwide populations and of 6122 individuals from Asia and adjoining populations were conducted based on allele frequency data and raw genotype data, respectively. All studied markers were found to be in Hardy-Weinberg equilibrium. The combined power of discrimination in the Gelao minority group was 0.999999999975, and the combined probability of exclusion was 0.9879. Our results from the forensic statistical parameters indicated that this Indel panel can be independently used as a powerful tool in forensic individual identification but can only be used as a complementary tool in paternity cases involving East Asians. We also found significant allele frequency differences between the Gelao and other continental populations with respect to the markers grouped in clusters ∼Ⅳ, suggesting that these can be used as forensic ancestry informative Indel markers to distinguish the Gelao from other continental populations. Genetic ancestry analyses demonstrated that Tai-Kadai-speaking Gelao share a dominant ancestry component with Hmong-Mien-speaking Miao. Our population genetic results from multidimensional scaling plots, principal component analysis, neighboring-joining tree construction and hierarchical clustering also suggested that the Zunyi Gelao are genetically closer to their linguistically or geographically close populations, such as the Han Chinese, Guizhou Bouyei and the Hubei Tujia, than to Turkic and Tibeto-Burman speakers.

Germline variation in BRCA1/2 is highly ethnic-specific: Evidence from over 30,000 Chinese hereditary breast and ovarian cancer patients

BRCA1 and BRCA2 play essential roles in maintaining the genome stability. Pathogenic germline mutations in these two genes disrupt their function, lead to genome instability and increase the risk of developing breast and ovarian cancers. BRCA mutations have been extensively screened in Caucasian populations, and the resulting information are used globally as the standard reference in clinical diagnosis, treatment and prevention of BRCA-related cancers. Recent studies suggest that BRCA mutations can be ethnic-specific, raising the question whether a Caucasian-based BRCA mutation information can be used as a universal standard worldwide, or whether an ethnicity-based BRCA mutation information system need to be developed for the corresponding ethnic populations. In this study, we used Chinese population as a model to test ethnicity-specific BRCA mutations considering that China has one of the latest numbers of breast cancer patients therefore BRCA mutation carriers. Through comprehensive data mining, standardization and annotation, we collected 1,088 distinct BRCA variants derived from over 30,000 Chinese individuals, one of the largest BRCA data set from a non-Caucasian population covering nearly all known BRCA variants in the Chinese population (https://dbBRCA-Chinese.fhs.umac.mo). Using this data, we performed multi-layered analyses to determine the similarities and differences of BRCA variation between Chinese and non-Chinese ethnic populations. The results show the substantial differences of BRCA data between Chinese and non-Chinese ethnicities. Our study indicates that the current Caucasian population-based BRCA data is not adequate to represent the BRCA status in non-Caucasian populations. Therefore, ethnic-based BRCA standards need to be established to serve for the non-Caucasian populations.

MultiWaver 2.0: modeling discrete and continuous gene flow to reconstruct complex population admixtures

Our goal in developing the MultiWaver software series was to be able to infer population admixture history under various complex scenarios. The earlier version of MultiWaver considered only discrete admixture models. Here, we report a newly developed version, MultiWaver 2.0, that implements a more flexible framework and is capable of inferring multiple-wave admixture histories under both discrete and continuous admixture models. MultiWaver 2.0 can automatically select an optimal admixture model based on the length distribution of ancestral tracks of chromosomes, and the program can estimate the corresponding parameters under the selected model. Specifically, for discrete admixture models, we used a likelihood ratio test (LRT) to determine the optimal discrete model and an expectation-maximization algorithm to estimate the parameters. In addition, according to the principles of the Bayesian Information Criterion (BIC), we compared the optimal discrete model with several continuous admixture models. In MultiWaver 2.0, we also applied a bootstrapping technique to provide levels of support for the chosen model and the confidence interval (CI) of the estimations of admixture time. Simulation studies validated the reliability and effectiveness of our method. Finally, the program performed well when applied to real datasets of typical admixed populations, such as African Americans, Uyghurs, and Hazaras.

AncestryPainter: A Graphic Program for Displaying Ancestry Composition of Populations and Individuals

Ancestry composition of populations and individuals has been extensively investigated in recent years due to advances in the genotyping and sequencing technologies. As the number of populations and individuals used for ancestry inference increases remarkably, say more than 100 populations or 1000 individuals, it is usually challenging to present the ancestry composition in a traditional way using a rectangular graph. To address this issue, we developed a program, AncestryPainter, which can illustrate the ancestry composition of populations and individuals with a rounded and nice-looking graph to save space. Individuals are depicted as length-fixed bars partitioned into colored segments representing different ancestries, and the population of interest can be highlighted as a pie chart in the center of the circle plot. In addition, AncestryPainter can also be applied to display personal ancestry in a way similar to that for displaying population ancestry. AncestryPainter is publicly available at http://www.picb.ac.cn/PGG/resource.php.

Genetic structure and polymorphisms of Gelao ethnicity residing in southwest china revealed by X-chromosomal genetic markers

X-chromosome short tandem repeat markers (X-STRs), due to their special inheritance models, physical location on a single chromosome and the absence of recombination in male meiosis, play an important role in forensic and population genetics. While a series of genetic analyses focusing on the genetic diversity and forensic characteristics of X-STRs are well studied for ethnically/linguistically diverse and demographically large Chinese populations, genetic evidence from Gelao ethnicity is still sparse. Here, we genotyped the first batch of 19 X-STRs in 513 Chinese Gelao individuals (265 females and 248 males), and reported genetic polymorphisms, forensic characteristics based on the single locus and seven linkage groups. DXS10135 with the highest PIC (0.9106) and LG1 (DXS10148-DXS10135-DXS8378) with the largest HD (0.9970) are polymorphic and informative. The CPDs in Gelao males and females are respectively larger than 0.999999999997095 and 0.99999999999999999999918, and the combined MECs are larger than 0.999999975715109. Subsequently, we investigated the population relationships among 14 Chinese populations based on 19 X-STRs and among 23 populations based on 11 overlapped X-STRs. Our results revealed genetic differentiations among Tibeto-Burman, Altaic and other Chinese homogenous populations, and demonstrated that Guizhou Gelao has the genetically closer relationships with Han Chinese and geographically close Guizhou Miao.

Towards broadening Forensic DNA Phenotyping beyond pigmentation: Improving the prediction of head hair shape from DNA

Human head hair shape, commonly classified as straight, wavy, curly or frizzy, is an attractive target for Forensic DNA Phenotyping and other applications of human appearance prediction from DNA such as in paleogenetics. The genetic knowledge underlying head hair shape variation was recently improved by the outcome of a series of genome-wide association and replication studies in a total of 26,964 subjects, highlighting 12 loci of which 8 were novel and introducing a prediction model for Europeans based on 14 SNPs. In the present study, we evaluated the capacity of DNA-based head hair shape prediction by investigating an extended set of candidate SNP predictors and by using an independent set of samples for model validation. Prediction model building was carried out in 9674 subjects (6068 from Europe, 2899 from Asia and 707 of admixed European and Asian ancestries), used previously, by considering a novel list of 90 candidate SNPs. For model validation, genotype and phenotype data were newly collected in 2415 independent subjects (2138 Europeans and 277 non-Europeans) by applying two targeted massively parallel sequencing platforms, Ion Torrent PGM and MiSeq, or the MassARRAY platform. A binomial model was developed to predict straight vs. non-straight hair based on 32 SNPs from 26 genetic loci we identified as significantly contributing to the model. This model achieved prediction accuracies, expressed as AUC, of 0.664 in Europeans and 0.789 in non-Europeans; the statistically significant difference was explained mostly by the effect of one EDAR SNP in non-Europeans. Considering sex and age, in addition to the SNPs, slightly and insignificantly increased the prediction accuracies (AUC of 0.680 and 0.800, respectively). Based on the sample size and candidate DNA markers investigated, this study provides the most robust, validated, and accurate statistical prediction models and SNP predictor marker sets currently available for predicting head hair shape from DNA, providing the next step towards broadening Forensic DNA Phenotyping beyond pigmentation traits.