Genetic evidence of paleolithic colonization and neolithic expansion of modern humans on the tibetan plateau

Evolutionary profiles and complex admixture landscape in East Asia: New insights from modern and ancient Y chromosome variation perspectives

Human Y-chromosomes are characterized by nonrecombination and uniparental inheritance, carrying traces of human history evolution and admixture. Large-scale population-specific genomic sources based on advanced sequencing technologies have revolutionized our understanding of human Y chromosome diversity and its anthropological and forensic applications. Here, we reviewed and meta-analyzed the Y chromosome genetic diversity of modern and ancient people from China and summarized the patterns of founding lineages of spatiotemporally different populations associated with their origin, expansion, and admixture. We emphasized the strong association between our identified founding lineages and language-related human dispersal events correlated with the Sino-Tibetan, Altaic, and southern Chinese multiple-language families related to the Hmong-Mien, Tai-Kadai, Austronesian, and Austro-Asiatic languages. We subsequently summarize the recent advances in translational applications in forensic and anthropological science, including paternal biogeographical ancestry inference (PBGAI), surname investigation, and paternal history reconstruction. Whole-Y sequencing or high-resolution panels with high coverage of terminal Y chromosome lineages are essential for capturing the genomic diversity of ethnolinguistically diverse East Asians. Generally, we emphasized the importance of including more ethnolinguistically diverse, underrepresented modern and spatiotemporally different ancient East Asians in human genetic research for a comprehensive understanding of the paternal genetic landscape of East Asians with a detailed time series and for the reconstruction of a reference database in the PBGAI, even including new technology innovations of Telomere-to-Telomere (T2T) for new genetic variation discovery.

Development of an ARMS multiplex real-time PCR assay for the detection of HLA-B*13:01 genotype by detecting highly specific SNPs

OBJECTIVES

HLA-B*13:01 was strongly associated with Dapsone Hypersensitivity Syndrome (DHS). This study aimed to develop and validate a rapid and economical method for HLA-B*13:01 genotyping.

METHODS

Two tubes multiplex real-time PCR detection system comprising amplification refractory mutation system primers and TaqMan probes was established for HLA-B*13:01 genotyping. Sequence-based typing was applied to validate the accuracy of the assay.

RESULTS

The accuracy of the assay was 100% for HLA-B*13:01 genotyping. The detection limit of the new method was 0.025 ng DNA. The positive rate of HLA-B*13:01 in the Bouyei (20%, n = 50) populations was significantly higher than that in the Uighur population (4%, n = 100), Han (4.5%, n = 200), and Tibetan (1%, n = 100) ( P  < 0.05).

CONCLUSION

The proposed method is rapid and reliable for HLA-B*13:01 screening in a clinical setting.

Differentiated genomic footprints suggest isolation and long-distance migration of Hmong-Mien populations

BACKGROUND

The underrepresentation of Hmong-Mien (HM) people in Asian genomic studies has hindered our comprehensive understanding of the full landscape of their evolutionary history and complex trait architecture. South China is a multi-ethnic region and indigenously settled by ethnolinguistically diverse HM, Austroasiatic (AA), Tai-Kadai (TK), Austronesian (AN), and Sino-Tibetan (ST) people, which is regarded as East Asia's initial cradle of biodiversity. However, previous fragmented genetic studies have only presented a fraction of the landscape of genetic diversity in this region, especially the lack of haplotype-based genomic resources. The deep characterization of demographic history and natural-selection-relevant genetic architecture of HM people was necessary.

RESULTS

We reported one HM-specific genomic resource and comprehensively explored the fine-scale genetic structure and adaptative features inferred from the genome-wide SNP data of 440 HM individuals from 33 ethnolinguistic populations, including previously unreported She. We identified solid genetic differentiation between HM people and Han Chinese at 7.64‒15.86 years ago (kya) and split events between southern Chinese inland (Miao/Yao) and coastal (She) HM people in the middle Bronze Age period and the latter obtained more gene flow from Ancient Northern East Asians. Multiple admixture models further confirmed that extensive gene flow from surrounding ST, TK, and AN people entangled in forming the gene pool of Chinese coastal HM people. Genetic findings of isolated shared unique ancestral components based on the sharing alleles and haplotypes deconstructed that HM people from the Yungui Plateau carried the breadth of previously unknown genomic diversity. We identified a direct and recent genetic connection between Chinese inland and Southeast Asian HM people as they shared the most extended identity-by-descent fragments, supporting the long-distance migration hypothesis. Uniparental phylogenetic topology and network-based phylogenetic relationship reconstruction found ancient uniparental founding lineages in southwestern HM people. Finally, the population-specific biological adaptation study identified the shared and differentiated natural selection signatures among inland and coastal HM people associated with physical features and immune functions. The allele frequency spectrum of cancer susceptibility alleles and pharmacogenomic genes showed significant differences between HM and northern Chinese people.

CONCLUSIONS

Our extensive genetic evidence combined with the historical documents supported the view that ancient HM people originated from the Yungui regions associated with ancient "Three-Miao tribes" descended from the ancient Daxi-Qujialing-Shijiahe people. Then, some have recently migrated rapidly to Southeast Asia, and some have migrated eastward and mixed respectively with Southeast Asian indigenes, Liangzhu-related coastal ancient populations, and incoming southward ST people. Generally, complex population migration, admixture, and adaptation history contributed to the complicated patterns of population structure of geographically diverse HM people.

A Tibetan group from Ngawa Tibetan and Qiang Autonomous Prefecture, southwest China, is rich in genetic polymorphisms at 36 autosomal STR loci and shares a complex genetic structure with other Chinese populations

The Tibetan people are ancient and populous, constituting the seventh-largest of the fifty-five ethnic minority groups in China. The Ngawa Tibetan and Qiang Autonomous Prefecture (NTQAP), situated on the border of northwest and southwest China, has its distinct group relationships. Short tandem repeat (STR) is extremely polymorphic and extensively used in the application of forensic medicine and population genetics. However, it is not clear the genetic information including linkage disequilibrium (LD) by 36 autosomal STR (A-STR) markers in the Tibetan group from NTQAP. The Tibetan population from NTQAP of southwest China was examined for 36 A-STR loci in the research. Every marker across the 36 A-STR loci was consistent with Hardy-Weinberg equilibrium (HWE). The results of the calculation revealed that the total discrimination power (TDP) is 1-2.2552 × 10-42 and the cumulative probability of exclusion (CPE) is 1-1.3031 × 10-16. Subsequently, a total of 345 alleles with allelic frequencies ranging from 0.00382 to 0.55343 were identified, and the allelic numbers varied from 5 in both the TH01 and TPOX markers to 28 in the SE33 locus. The Ngawa Tibetan population, along with other Chinese populations, exhibited influences from historical factors and regional distribution, as indicated by the results of population genetics analysis. We thus first explored the genetic characteristics and correlated forensic parameters of the 36 A-STR markers in NTQAP to fill the gap in the Tibetan population. It was discovered that these 36 autosomal STR markers supplemented forensic STR databases and offered extremely valuable polymorphisms for Chinese forensic applications, such as parentage testing and personal identification. Moreover, the study would contribute additional information regarding the substructure and diversity in the Chinese population.

Whole mitochondrial genome analysis in highland Tibetans: further matrilineal genetic structure exploration

Introduction: The Qinghai-Tibet Plateau is one of the last terrestrial environments conquered by modern humans. Tibetans are among the few high-altitude settlers in the world, and understanding the genetic profile of Tibetans plays a pivotal role in studies of anthropology, genetics, and archaeology. Methods: In this study, we investigated the maternal genetic landscape of Tibetans based on the whole mitochondrial genome collected from 145 unrelated native Lhasa Tibetans. Molecular diversity indices, haplotype diversity (HD), Tajima's D and Fu's Fs were calculated and the Bayesian Skyline Plot was obtained to determining the genetic profile and population fluctuation of Lhasa Tibetans. To further explore the genetic structure of Lhasa Tibetans, we collected 107 East Asian reference populations to perform principal component analysis (PCA), multidimensional scaling (MDS), calculated Fst values and constructed phylogenetic tree. Results: The maternal genetic landscape of Tibetans showed obvious East Asian characteristics, M9a (28.28%), R (11.03%), F1 (12.41%), D4 (9.66%), N (6.21%), and M62 (4.14%) were the dominant haplogroups. The results of PCA, MDS, Fst and phylogenetic tree were consistent: Lhasa Tibetans clustered with other highland Tibeto-Burman speakers, there was obvious genetic homogeneity of Tibetans in Xizang, and genetic similarity between Tibetans and northern Han people and geographically adjacent populations was found. In addition, specific maternal lineages of Tibetans also be determined in this study. Discussion: In general, this study further shed light on long-time matrilineal continuity on the Tibetan Plateau and the genetic connection between Tibetans and millet famers in the Yellow River Basin, and further revealed that multiple waves of population interaction and admixture during different historical periods between lowland and highland populations shaped the maternal genetic profile of Tibetans.

Unveiling the genetic landscape of high-altitude adaptive ethnic groups with polymorphic markers: Implications of comprehensive forensic appraisals and population genetic investigations

Background

Deletion/insertion polymorphisms (DIPs), a novel class of biomarker, have been widely utilized in forensic areas for individual identification, paternity tests, and ancestral origin inference due to its applicability to degraded samples and low mutation rates. Despite the availability of a well-established commercial kit, the Investigator® DIPplex kit (Qiagen), certain loci exhibit limited levels of polymorphisms in East Asian populations, particularly in Chinese populations.

Objective

This dissertation seeks to undertake a comprehensive evaluation about the forensic efficiency of a self-developed multiplex amplification system in high-altitude adaptive ethnic groups of China. Healthy unrelated Tibetan individuals residing in Tibet Autonomous Region and Qinghai Province were genotyped using previously reported 43 deletion/insertion polymorphism loci. Forensic statistical analyses including allele frequencies and forensic parameters were conducted in the two Tibetan groups, and the genetic relatedness of the studied groups with reference populations from the 1000 Genomes Project Phase 3 were investigated.

Results

Forensic statistical results showed that the polymorphism information content values of the 43 deletion/insertion polymorphism loci in the two Tibetan groups exceeded 0.35. Moreover, the combined power of discrimination using the 43 deletion/insertion polymorphism loci was calculated to be 0.9999999999999999984 in the Qinghai Tibetan group and 0.9999999999999999921 in the Tibet Tibetan group. The cumulative power of exclusion using the 43 deletion/insertion polymorphism loci was calculated to be 0.999782512 in the Qinghai Tibetan group and 0.999886205 in the Tibet Tibetan group. Analysis of population genetics demonstrated that the two studied Tibetan groups shared close genetic relationships with East Asia populations.

Conclusion

The set of 43 deletion/insertion polymorphism loci exhibited remarkable forensic efficacy, rendering it a promising tool for forensic practice. Population genetic analyses indicated that the two Tibetan groups had closer genetic affinities to East Asian populations.

Developmental validation of a high-resolution panel genotyping 639 Y-chromosome SNP and InDel markers and its evolutionary features in Chinese populations

Uniparental-inherited haploid genetic marker of Y-chromosome single nucleotide polymorphisms (Y-SNP) have the power to provide a deep understanding of the human evolutionary past, forensic pedigree, and bio-geographical ancestry information. Several international cross-continental or regional Y-panels instead of Y-whole sequencing have recently been developed to promote Y-tools in forensic practice. However, panels based on next-generation sequencing (NGS) explicitly developed for Chinese populations are insufficient to represent the Chinese Y-chromosome genetic diversity and complex population structures, especially for Chinese-predominant haplogroup O. We developed and validated a 639-plex panel including 633 Y-SNPs and 6 Y-Insertion/deletions, which covered 573 Y haplogroups on the Y-DNA haplogroup tree. In this panel, subgroups from haplogroup O accounted for 64.4% of total inferable haplogroups. We reported the sequencing metrics of 354 libraries sequenced with this panel, with the average sequencing depth among 226 individuals being 3,741×. We illuminated the high level of concordance, accuracy, reproducibility, and specificity of the 639-plex panel and found that 610 loci were genotyped with as little as 0.03 ng of genomic DNA in the sensitivity test. 94.05% of the 639 loci were detectable in male-female mixed DNA samples with a mix ratio of 1:500. Nearly all of the loci were genotyped correctly when no more than 25 ng/μL tannic acid, 20 ng/μL humic acid, or 37.5 μM hematin was added to the amplification mixture. More than 80% of genotypes were obtained from degraded DNA samples with a degradation index of 11.76. Individuals from the same pedigree shared identical genotypes in 11 male pedigrees. Finally, we presented the complex evolutionary history of 183 northern Chinese Hans and six other Chinese populations, and found multiple founding lineages that contributed to the northern Han Chinese gene pool. The 639-plex panel proved an efficient tool for Chinese paternal studies and forensic applications.

Polygenic adaptation leads to a higher reproductive fitness of native Tibetans at high altitude

The adaptation of Tibetans to high-altitude environments has been studied extensively. However, the direct assessment of evolutionary adaptation, i.e., the reproductive fitness of Tibetans and its genetic basis, remains elusive. Here, we conduct systematic phenotyping and genome-wide association analysis of 2,252 mother-newborn pairs of indigenous Tibetans, covering 12 reproductive traits and 76 maternal physiological traits. Compared with the lowland immigrants living at high altitudes, indigenous Tibetans show better reproductive outcomes, reflected by their lower abortion rate, higher birth weight, and better fetal development. The results of genome-wide association analyses indicate a polygenic adaptation of reproduction in Tibetans, attributed to the genomic backgrounds of both the mothers and the newborns. Furthermore, the EPAS1-edited mice display higher reproductive fitness under chronic hypoxia, mirroring the situation in Tibetans. Collectively, these results shed new light on the phenotypic pattern and the genetic mechanism of human reproductive fitness in extreme environments.

Maternal genetic history of ancient Tibetans over the past 4000 years

The settlement of the Tibetan Plateau epitomizes human adaptation to a high-altitude environment that poses great challenges to human activity. Here, we reconstruct a 4000-year maternal genetic history of Tibetans using 128 ancient mitochondrial genome data from 37 sites in Tibet. The phylogeny of haplotypes M9a1a, M9a1b, D4g2, G2a'c, and D4i show that ancient Tibetans share the most recent common ancestor with ancient Middle and Upper Yellow River populations around the Early and Middle Holocene. In addition, the connections between Tibetans and Northeastern Asians vary over the past 4000 years, with a stronger matrilineal connection between the two during 4000 BP-3000 BP, and a weakened connection after 3000 BP, that are coincident with climate change, followed by a reinforced connection after the Tubo period (1400 BP-1100 BP). Besides, an over 4000-year matrilineal continuity is observed in some of the maternal lineages. We also find the maternal genetic structure of ancient Tibetans is correlated to the geography and interactions between ancient Tibetans and ancient Nepal and Pakistan populations. Overall, the maternal genetic history of Tibetans can be characterized as a long-term matrilineal continuity with frequent internal and external population interactions that are dynamically shaped by geography, climate changes, as well as historical events.

Deep phenotyping of 11,880 highlanders reveals novel adaptive traits in native Tibetans

Tibetans are the ideal population to study genetic adaptation in extreme environments. Here, we performed systematic phenotyping of 11,880 highlanders, covering 133 quantitative traits of 13 organ systems. We provided a comprehensive phenotypic atlas by comparing altitude adaptation and altitude acclimatization. We found the differences between adaptation and acclimatization are quantitative rather than qualitative, with a whole-system "blunted effect" seen in the adapted Tibetans. We characterized twelve different functional changes between adaptation and acclimatization. More importantly, we established a landscape of adaptive phenotypes of indigenous Tibetans, including 45 newly identified Tibetan adaptation-nominated traits, involving specific changes of Tibetans in internal organ state, metabolism, eye morphology, and skin pigmentation. In addition, we observed a sex-biased pattern between altitude acclimatization and adaptation. The generated atlas of phenotypic landscape provides new insights into understanding of human adaptation to high-altitude environments, and it serves as a valuable blueprint for future medical and physiological studies.

Associations between high-altitude adaptation and risk of cardiovascular diseases: a bidirectional Mendelian randomization study

High-altitude adaptation (HAA) was reported to be significantly associated with reduced risks for multiple cardiovascular diseases (CVDs). However, the causality and direction of the associations are largely uncharacterized. We aimed to examine the potential causal relationships between HAA and six types of CVD, including coronary artery disease (CAD), cerebral aneurysm, ischemic stroke, peripheral artery disease, arrhythmia and atrial fibrillation. We obtained the summary data from largest available genome-wide association study of HAA and six types of CVD. Two-sample bidirectional Mendelian randomization (MR) analyses were performed to infer the causality between them. In the sensitivity analyses, MR-Egger regression analyses and MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) global analyses were used to assess the pleiotropic effects; Cochran's Q tests were used to test the heterogeneity by inverse variance-weighted (IVW) and MR-Egger methods; and the leave-one-out analyses were used to examine whether some single nucleotide polymorphisms (SNPs) could influence the results independently. The MR main analyses showed that the genetically instrumented HAA was significantly causally associated with the reduced risks of CAD (odds ratio [OR] = 0.029; 95% confidence interval [CI] = 0.004-0.234; P = 8.64 × 10-4). In contrast, there was no statistically significant relationship between CVDs and HAA. Our findings provide evidence for the causal effects of HAA on the reduced risks of CAD. However, there is no causality of CVDs on HAA. These findings might be helpful in developing the prevention and intervention strategies for CAD.

Searching for molecular hypoxia sensors among oxygen-dependent enzymes

The ability to sense and respond to changes in cellular oxygen levels is critical for aerobic organisms and requires a molecular oxygen sensor. The prototypical sensor is the oxygen-dependent enzyme PHD: hypoxia inhibits its ability to hydroxylate the transcription factor HIF, causing HIF to accumulate and trigger the classic HIF-dependent hypoxia response. A small handful of other oxygen sensors are known, all of which are oxygen-dependent enzymes. However, hundreds of oxygen-dependent enzymes exist among aerobic organisms, raising the possibility that additional sensors remain to be discovered. This review summarizes known and potential hypoxia sensors among human O2-dependent enzymes and highlights their possible roles in hypoxia-related adaptation and diseases.

Large-scale genome sequencing redefines the genetic footprints of high-altitude adaptation in Tibetans

BACKGROUND

Tibetans are genetically adapted to high-altitude environments. Though many studies have been conducted, the genetic basis of the adaptation remains elusive due to the poor reproducibility for detecting selective signatures in the Tibetan genomes.

RESULTS

Here, we present whole-genome sequencing (WGS) data of 1001 indigenous Tibetans, covering the major populated areas of the Qinghai-Tibetan Plateau in China. We identify 35 million variants, and more than one-third of them are novel variants. Utilizing the large-scale WGS data, we construct a comprehensive map of allele frequency and linkage disequilibrium and provide a population-specific genome reference panel, referred to as 1KTGP. Moreover, with the use of a combined approach, we redefine the signatures of Darwinian-positive selection in the Tibetan genomes, and we characterize a high-confidence list of 4320 variants and 192 genes that have undergone selection in Tibetans. In particular, we discover four new genes, TMEM132C, ATP13A3, SANBR, and KHDRBS2, with strong signals of selection, and they may account for the adaptation of cardio-pulmonary functions in Tibetans. Functional annotation and enrichment analysis indicate that the 192 genes with selective signatures are likely involved in multiple organs and physiological systems, suggesting polygenic and pleiotropic effects.

CONCLUSIONS

Overall, the large-scale Tibetan WGS data and the identified adaptive variants/genes can serve as a valuable resource for future genetic and medical studies of high-altitude populations.

The matrilineal ancestry of Nepali populations

The Tibetan plateau and high mountain ranges of Nepal are one of the challenging geographical regions inhabited by modern humans. While much of the ethnographic and population-based genetic studies were carried out to investigate the Tibetan and Sherpa highlanders, little is known about the demographic processes that enabled the colonization of the hilly areas of Nepal. Thus, the present study aimed to investigate the past demographic events that shaped the extant Nepalese genetic diversity using mitochondrial DNA (mtDNA) variations from ethnic Nepalese groups. We have analyzed mtDNA sequences of 999 Nepalese and compared data with 38,622 published mtDNA sequences from rest of the world. Our analysis revealed that the genomic landscapes of prehistoric Himalayan settlers of Nepal were similar to that of the low-altitude extant Nepalese (LAN), especially Newar and Magar population groups, but differ from contemporary high-altitude Sherpas. LAN might have derived their East Eurasian ancestry mainly from low-altitude Tibeto-Burmans, who likely have migrated from East Asia and assimilated across the Eastern Himalayas extended from the Eastern Nepal to the North-East of India, Bhutan, Tibet and Northern Myanmar. We also identified a clear genetic sub-structure across different ethnic groups of Nepal based on mtDNA haplogroups and ectodysplasin-A receptor (EDAR) gene polymorphism. Our comprehensive high-resolution mtDNA-based genetic study of Tibeto-Burman communities reconstructs the maternal origins of prehistoric Himalayan populations and sheds light on migration events that have brought most of the East Eurasian ancestry to the present-day Nepalese population.

A single-tube multiplex real-time PCR for HLA-B*38:02 genotype by detecting highly specific SNPs

Purpose: HLA-B*38:02 is closely related to carbimazole/methimazole-induced agranulocytosis. This study aimed to develop and validate a rapid and economical method for HLA-B*38:02 genotyping. Methods: A single-tube multiplex real-time PCR detection system comprising amplification refractory mutation system primers and TaqMan probes was established for HLA-B*38:02 genotyping. Sequence-based typing was applied to validate the accuracy of the assay. Results: The accuracy of the assay was 100% for HLA-B*38:02 genotyping. The detection limit of the new method was 0.05 ng DNA. The positive rate of HLA-B*38:02 in the Han (8%, n = 100), Bouyei (17.8%, n = 90) and Tibetan (12.7%, n = 110) populations was significantly higher than that in the Uighur population (1%, n = 100) (p < 0.05). Conclusion: The proposed method is rapid and reliable for HLA-B*38:02 screening in a clinical setting.

Comparison of Cardiorespiratory Fitness of Chinese Tibetan Adolescents with Their Han Counterparts: A Cross-Sectional Retrospective Study

Cardiorespiratory fitness (CRF) is a core element of healthy physical fitness. Foreign attention to CRF in adolescents at different altitudes is high, while less research has been conducted on Chinese adolescents. In order to compare the CRF of Chinese Tibetan adolescents with their Han counterparts born and raised at high altitude and Chinese Han adolescents at sea level. A total of 2748 participants, including Chinese Tibetan adolescents, Chinese Han adolescents born and raised at high altitudes, and Chinese Han adolescents at sea level aged 12-18 years old, were obtained using convenience sampling and random cluster sampling. The method of the 20 m shuttle run test (20 m SRT) test was used to derive VO_2max by equation. One-way ANOVA and LSD methods were conducted, and effect sizes were calculated to compare the CRF of the three types of adolescents. Regression analysis was used to analyze the relationship between altitude and VO_2max. The VO_2max scores of Chinese Tibetan adolescents and Chinese Han adolescents at sea level were higher than Chinese Han adolescents born and raised at high altitudes. For both boys and girls, the VO_2max scores of Chinese Tibetan adolescents exceeded Chinese Han adolescents at sea level after the age of 16 years old. Regression analysis showed that altitude was inversely associated with VO_2max. The pace of lung growth may distinguish Chinese Tibetan adolescents from Chinese Han adolescents born and raised at high altitudes. The results of the study suggest that we should focus on the changes in CRF in adolescents at different altitudes and should adopt different CRF interventions for adolescents at different altitudes.

Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia

As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.

Cultural and demic co-diffusion of Tubo Empire on Tibetan Plateau

A high point of Tibetan Plateau (TP) civilization, the expansive Tubo Empire (618-842 AD) wielded great influence across ancient western China. However, whether the Tubo expansion was cultural or demic remains unclear due to sparse ancient DNA sampling. Here, we reported ten ancient genomes at 0.017- to 0.867-fold coverages from the Dulan site with typical Tubo archaeological culture dating to 1308-1130 BP. Nine individuals from three different grave types have close relationship with previously reported ancient highlanders from the southwestern Himalayas and modern core-Tibetan populations. A Dulan-related Tubo ancestry contributed overwhelmingly (95%-100%) to the formation of modern Tibetans. A genetic outlier with dominant Eurasian steppe-related ancestry suggesting a potential population movement into the Tubo-controlled regions from Central Asia. Together with archeological evidence from burial styles and customs, our study suggested the impact of the Tubo empire on the northeast edge of the TP involved both cultural and demic diffusion.

Genetic adaptation of skin pigmentation in highland Tibetans

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Ancient DNA from Tubo Kingdom-related tombs in northeastern Tibetan Plateau revealed their genetic affinity to both Tibeto-Burman and Altaic populations

The rise of the Tubo Kingdom is considered as the key period for the formation of modern groups on the Tibetan Plateau. The ethnic origin of the residents of the Tubo Kingdom is quite complex, and their genetic structure remains unclear. The tombs of the Tubo Kingdom period in Dulan County, Qinghai Province, dating back to the seventh century, are considered to be the remains left by Tubo conquerors or the Tuyuhun people dominated by the Tubo Kingdom. The human remains of these tombs are ideal materials for studying the population dynamics in the Tubo Kingdom. In this paper, we analyzed the genome-wide data of eight remains from these tombs by shotgun sequencing and multiplex PCR panels and compared the results with data of available ancient and modern populations across East Asia. Genetic continuity between ancient Dulan people with ancient Xianbei tribes in Northeast Asia, ancient settlers on the Tibetan Plateau, and modern Tibeto-Burman populations was found. Surprisingly, one out of eight individuals showed typical genetic features of populations from Central Asia. In summary, the genetic diversity of ancient Dulan people and their affiliations with other populations provide an example of the complex origin of the residents in the Tubo Kingdom and their long-distance connection with populations in a vast geographic region across ancient Asia.

Genomic insight into the population history and biological adaptations of high-altitude Tibetan highlanders in Nagqu

Tibetan, one of the largest indigenous populations living in the high-altitude region of the Tibetan Plateau (TP), has developed a suite of physiological adaptation strategies to cope with the extreme highland environment in TP. Here, we reported genome-wide SNP data from 48 Kham-speaking Nagqu Tibetans and analyzed it with published data from 1,067 individuals in 167 modern and ancient populations to characterize the detailed Tibetan subgroup history and population substructure. Overall, the patterns of allele sharing and haplotype sharing suggested (1) the relatively genetic homogeny between the studied Nagqu Tibetans and ancient Nepalese as well as present-day core Tibetans from Lhasa, Nagqu, and Shigatse; and (2) the close relationship between our studied Kham-speaking Nagqu Tibetans and Kham-speaking Chamdo Tibetans. The fitted qpAdm models showed that the studied Nagqu Tibetans could be fitted as having the main ancestry from late Neolithic upper Yellow River millet farmers and deeply diverged lineages from Southern East Asians (represented by Upper Paleolithic Guangxi_Longlin and Laos_Hoabinhian), and a non-neglectable western Steppe herder-related ancestry (∼3%). We further scanned the candidate genomic regions of natural selection for our newly generated Nagqu Tibetans and the published core Tibetans via FST, iHS, and XP-EHH tests. The genes overlapping with these regions were associated with essential human biological functions such as immune response, enzyme activity, signal transduction, skin development, and energy metabolism. Together, our results shed light on the admixture and evolutionary history of Nagqu Tibetan populations.

Genetic evidence of tri-genealogy hypothesis on the origin of ethnic minorities in Yunnan

BACKGROUND

Yunnan is located in Southwest China and consists of great cultural, linguistic, and genetic diversity. However, the genomic diversity of ethnic minorities in Yunnan is largely under-investigated. To gain insights into population history and local adaptation of Yunnan minorities, we analyzed 242 whole-exome sequencing data with high coverage (~ 100-150 ×) of Yunnan minorities representing Achang, Jingpo, Dai, and Deang, who were linguistically assumed to be derived from three ancient lineages (the tri-genealogy hypothesis), i.e., Di-Qiang, Bai-Yue, and Bai-Pu.

RESULTS

Yunnan minorities show considerable genetic differences. Di-Qiang populations likely migrated from the Tibetan area about 6700 years ago. Genetic divergence between Bai-Yue and Di-Qiang was estimated to be 7000 years, and that between Bai-Yue and Bai-Pu was estimated to be 5500 years. Bai-Pu is relatively isolated, but gene flow from surrounding Di-Qiang and Bai-Yue populations was also found. Furthermore, we identified genetic variants that are differentiated within Yunnan minorities possibly due to the living circumstances and habits. Notably, we found that adaptive variants related to malaria and glucose metabolism suggest the adaptation to thalassemia and G6PD deficiency resulting from malaria resistance in the Dai population.

CONCLUSIONS

We provided genetic evidence of the tri-genealogy hypothesis as well as new insights into the genetic history and local adaptation of the Yunnan minorities.

Uniparental Genetic Analyses Reveal Multi-Ethnic Background of Dunhuang Foyemiaowan Population (220–907 CE) With Typical Han Chinese Archaological Culture

The relationship between archeological culture and ethnicity is invariably complex. This is especially the case for periods of national division and rapid inter-ethnic exchange, such as China’s Sixteen Kingdoms (304–439 CE) and Northern and Southern Dynasties (420–589 CE). Going by tomb shape and grave goods, the Foyemiaowan cemetery at Dunhuang exhibits a typical third–tenth century Han style. Despite this, the ethnic makeup of the Foyemiaowan population has remained unclear. We therefore analyzed 485 Y-chromosomal SNPs and entire mitochondrial genomes of 34 Foyemiaowan samples. Our study yielded the following discoveries: (1) principal component analysis revealed that the Foyemiaowan population was closely clustered with Tibeto-Burman populations on the paternal side and close to Mongolic-speaking populations on the maternal side; (2) lineage comparisons at the individual level showed that the Foyemiaowan population consisted of primarily Tibeto-Burman and Han Chinese related lineages (Oα-M117, 25%;Oβ-F46, 18.75%), partially Altaic speaking North Eurasian lineages (N-F1206, 18.75%) and a slight admixture of southern East Asian lineages (O1b1a2-Page59, 6.25%; O1b1a1-PK4, 3.13%). Similarly, the maternal gene pool of Foyemiaowan contained northern East Asian (A, 4.17%; CZ, 16.67%; D, 20.83%; G, 4.17%; M9, 4.17%), southern East Asian (B, 12.51%; F, 20.83%) and western Eurasian (H, 4.17%; J, 4.17%) related lineages; (3) we discovered a relatively high genetic diversity among the Foyemiaowan population (0.891) in our ancient reference populations, indicating a complex history of population admixture. Archeological findings, stable isotope analysis and historical documents further corroborated our results. Although in this period China’s central government had relinquished control of the Hexi Corridor and regional non-Han regimes became the dominant regional power, Foyemiaowan’s inhabitants remained strongly influenced by Han culture.

Regulatory Role of N6-Methyladenosine in Longissimus Dorsi Development in Yak

N6-methyladenine (m6A) RNA undergoes epigenetic modification, which is the most extensive intermediate chemical modification in mRNA. Although this modification occurs in all living organisms, it is the most widely studied among mammals. However, to date, no study has investigated the m6A transcriptome-wide map of yak and its potential biological functions in muscle development. In this study, the differences of m6A methylation and gene expression in yak muscle development belonging to three age groups, namely 3 years (group A), 6 months (group M), and 90-day-old fetuses (group E), were determined by using methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and RNA sequencing (RNA-Seq). In these three groups, a total of 6,278 (A), 9,298 (E), and 9,584 (M) m6A peaks were identified, with average densities between 1.02 and 2.01. m6A peaks were mostly enriched in the stop codon, 3′ untranslated region (UTR) region, and inner long exon region with consensus motifs of UGACA. In all the three stages, the m6A peak enrichment level was significantly negatively correlated with mRNA abundance (Pearson's correlation coefficient r = −0.22 to −0.32, p < 10⁻¹⁶). The functional enrichment of genes consistently modified by m6A methylation, particularly those genes that regulate cell differentiation as well as muscle growth and development, was observed at all three stages. Moreover, m6A abundance was negatively associated with gene expression levels, indicating that m6A might play a vital role in modulating gene expression during yak muscle development. This comprehensive map thus provides a solid foundation for determining the potential functional role of m6A RNA modification in yak muscle growth.

Ancient genomes from the Himalayas illuminate the genetic history of Tibetans and their Tibeto-Burman speaking neighbors

Present-day Tibetans have adapted both genetically and culturally to the high altitude environment of the Tibetan Plateau, but fundamental questions about their origins remain unanswered. Recent archaeological and genetic research suggests the presence of an early population on the Plateau within the past 40 thousand years, followed by the arrival of subsequent groups within the past 10 thousand years. Here, we obtain new genome-wide data for 33 ancient individuals from high elevation sites on the southern fringe of the Tibetan Plateau in Nepal, who we show are most closely related to present-day Tibetans. They derive most of their ancestry from groups related to Late Neolithic populations at the northeastern edge of the Tibetan Plateau but also harbor a minor genetic component from a distinct and deep Paleolithic Eurasian ancestry. In contrast to their Tibetan neighbors, present-day non-Tibetan Tibeto-Burman speakers living at mid-elevations along the southern and eastern margins of the Plateau form a genetic cline that reflects a distinct genetic history. Finally, a comparison between ancient and present-day highlanders confirms ongoing positive selection of high altitude adaptive alleles.

Denisovans and Homo sapiens on the Tibetan Plateau: dispersals and adaptations

Recent archaeological discoveries suggest that both archaic Denisovans and Homo sapiens occupied the Tibetan Plateau earlier than expected. Genetic studies show that a pulse of Denisovan introgression was involved in the adaptation of Tibetan populations to high-altitude hypoxia. These findings challenge the traditional view that the plateau was one of the last places on earth colonized by H. sapiens and warrant a reappraisal of the population history of this highland. Here, we integrate archaeological and genomic evidence relevant to human dispersal, settlement, and adaptation in the region. We propose two testable models to address the peopling of the plateau in the broader context of H. sapiens dispersal and their encounters with Denisovans in Asia.

Same total normal forms sperm counts of males from Lhasa and Shanghai, China

Male infertility may be caused by genetic and/or environmental factors that impair spermatogenesis and sperm maturation. High-altitude (HA) hypoxic environments represent one of the most serious challenges faced by humans that reside in these areas. To assess the influence of the plateau environment on semen parameters, 2,798 males, including 1,111 native Tibetans and 1,687 Han Chinese individuals living in the plains (HCILP) who underwent pre-pregnancy checkups, were enrolled in this study. The semen samples of males were evaluated to determine conventional sperm parameters, sperm morphology, and sperm movement. Reproductive endocrine hormones (REHs) were detected in 474 males, including 221 Tibetans and 253 HCILP. Due to recurrent abortions in partners, the DNA fragmentation index (DFI) of 133 native Tibetans and 393 HCILP individuals was further compared. Luteinizing hormone (LH) (4.94 ± 2.12 vs. 3.29 ± 1.43 U/L), prolactin (11.34 ± 3.87 vs. 8.97 ± 3.48 nmol/L), E2/T (0.22 ± 0.11 vs 0.11 ± 0.05), median total sperm motility (61.20% vs. 51.56%), and DFI (23.11% vs. 7.22%) were higher in males from plateau areas while median progressive motility (PR) (35.60% vs. 41.12%), total number of PR sperms (51.61 vs. 59.63 mil/ejaculate), percentage of normal form sperms (3.70% vs. 6.00%), curvilinear velocity (36.10 vs. 48.97 μm/s), straight-line (rectilinear) velocity (14.70 vs. 31.52 μm/s), estradiol (103.82 ± 45.92 vs. 146.01 ± 39.73 pmol/L), progesterone (0.29 ± 0.27 vs. 2.22 ± 0.84 nmol/L), testosterone (4.90 ± 1.96 vs. 14.36 ± 5.24 nmol/L), and testosterone secretion index (ratio of testosterone to LH) (33.45 ± 22.86 vs 145.78 ± 73.41) were lower than those in males from the plains. There was no difference in median total sperm number (157.76 vs. 151.65 mil/mL), sperm concentration (52.40 vs. 51.79 mil/mL), volume (3.10 vs. 3.10 mL), total normal form sperms (5.91 vs. 6.58 mil/ejaculate, p50), and follicle-stimulating hormone (FSH) levels (4.13 ± 2.55 U/L vs 3.82 ± 2.35 U/L) between the two groups of males. The REH and sperm parameters of males from HA hypoxic environments were adaptively altered. Although the total number of PR sperm decreased and DFI increased, the Tibetan population that lives at HAs has been found to grown continuously and rapidly. These results supplement prior findings regarding the impact of HA on male reproductive function.

Genetic characterization of the highlander Tibetan population from Qinghai-Tibet Plateau revealed by X chromosomal STRs

Tibetans are considered an East Asian ethnic group and primarily live in the high Tibetan plateau, the western Sichuan and Yunnan mountains of central and southern China, and areas throughout the Himalayas and around the Tibetan plateau. These people exhibit rare molecular machinery that allows them to adapt to hypoxic environments in the Qinghai-Tibet Plateau and make them a potential candidate for providing insights related to medical genetic, molecular medicine and human population studies. In the current study, we have genotyped 549 individuals with Investigator Argus X-12 Kit. For 12 X-STRs, a total of 174 unique alleles were found, among them DXS10134 and DXS10135 were the most polymorphic loci. All of the loci were in Hardy-Weinberg Equilibrium (HWE). The numbers of observed haplotypes in Highlander Tibetans males were 161,112, 96 and 108, respectively, whereas haplotype diversities (HD) were 0.9959, 0.9880, 0.9809 and 0.9873, respectively. The combined discrimination power for males (PDm) was 0.999 999 99701 and for females (PDf) was 0.999 999 999 999 9958. This study represents an extensive report on X chromosomal STR markers variation in the Highlander Tibetans population for forensic applications and population genetic studies.

Improving the regional Y-STR haplotype resolution utilizing haplogroup-determining Y-SNPs and the application of machine learning in Y-SNP haplogroup prediction in a forensic Y-STR database: A pilot study on male Chinese Yunnan Zhaoyang Han population

Improving the resolution of the current widely used Y-chromosomal short tandem repeat (Y-STR) dataset is of great importance for forensic investigators, and the current approach is limited, except for the addition of more Y-STR loci. In this research, a regional Y-DNA database was investigated to improve the Y-STR haplotype resolution utilizing a Y-SNP Pedigree Tagging System that includes 24 Y-chromosomal single nucleotide polymorphism (Y-SNP) loci. This pilot study was conducted in the Chinese Yunnan Zhaoyang Han population, and 3473 unrelated male individuals were enrolled. Based on data on the male haplogroups under different panels, the matched or near-matching (NM) Y-STR haplotype pairs from different haplogroups indicated the critical roles of haplogroups in improving the regional Y-STR haplotype resolution. A classic median-joining network analysis was performed using Y-STR or Y-STR/Y-SNP data to reconstruct population substructures, which revealed the ability of Y-SNPs to correct misclassifications from Y-STRs. Additionally, population substructures were reconstructed using multiple unsupervised or supervised dimensionality reduction methods, which indicated the potential of Y-STR haplotypes in predicting Y-SNP haplogroups. Haplogroup prediction models were built based on nine publicly accessible machine-learning (ML) approaches. The results showed that the best prediction accuracy score could reach 99.71% for major haplogroups and 98.54% for detailed haplogroups. Potential influences on prediction accuracy were assessed by adjusting the Y-STR locus numbers, selecting Y-STR loci with various mutabilities, and performing data processing. ML-based predictors generally presented a better prediction accuracy than two available predictors (Nevgen and EA-YPredictor). Three tree models were developed based on the Yfiler Plus panel with unprocessed input data, which showed their strong generalization ability in classifying various Chinese Han subgroups (validation dataset). In conclusion, this study revealed the significance and application prospects of Y-SNP haplogroups in improving regional Y-STR databases. Y-SNP haplogroups can be used to discriminate NM Y-STR haplotype pairs, and it is important for forensic Y-STR databases to develop haplogroup prediction tools to improve the accuracy of biogeographic ancestry inferences.

The impact of prehistoric human dispersals on the presence of tobacco-related oral cancer in Northeast India

BACKGROUND

Northeast (NE) India is a subject of debate for predicting its involvement in prehistoric anatomically modern human (AMH) dispersal. The unique lifestyle and genetic characteristics of native ethnic groups in this region are believed to be responsible for their susceptibility to tobacco-related oral cancer (TrOC). The present study assessed mitochondrial macro-haplogroup (mHG) diversity and TrOC susceptibility autosomal loci to evaluate the impact of prehistoric AMH dispersal on the present day's high TrOC prevalence in major NE Indian ethnics.

METHODS

We considered 175 unrelated individuals from 35 ethnic groups and previously published 374 sequences for sequencing-based assessment of mtDNA-based marker by subsequent analyses like haplogroup diversity, phylogenetic, genetic structure by AMOVA, and MDS, descriptive statistics of demographic parameters, and migration analysis. Besides, we selected prolonged tobacco-chewing 124 case-control individuals from similar ethnic backgrounds for genotyping 115 autosomal loci in Sequenom iPLEX MassARRAY™ platform and mined 1000genome data (n = 398) for consequent global admixture and ancestry-specific allele frequencies-based analyses.

RESULTS

Our mtDNA-based findings suggested that NE populations were distinct from other Indian populations, owing to the first wave of migration from ancient southern China (∼54kya) and two successive spatial expansion events at ∼45kya and ∼43kya. Consequently, it probably acted as another source for prehistoric AMH dispersal in N/NE Asia. Besides, the second wave of back-migration from SE Asia (∼40kya) probably replaced the mitochondrial footprints of survivors from the first migrants and introduced the TrOC susceptibility traits in this region. Afterward, the autosomal marker-based observations on the transition of the disease-associated admixture component 'K6' from SE Asia reconfirmed these results. Moreover, we also observed that the mitochondrial mHG 'R' is significantly associated with the risk of TrOC (OR > 9.5) in NE India. Furthermore, the possible onset of the phenotypic expression of those traits was predicted at ∼4kya, thus, contributing to present-day's TrOC prevalence.

CONCLUSIONS

This study reflects its uniqueness by revealing an updated AMH dispersal route for the peopling in and out of NE India, which probably introduced the disease-causing traits in the ancestral NE Indian population. Those traits were then imprinted in their genome to get transferred through their respective generations, forming the present-day's TrOC-prevalent NE Indian population.

A comparative analysis of differentially expressed mRNAs, miRNAs and circRNAs provides insights into the key genes involved in the high-altitude adaptation of yaks

BACKGROUND

Yaks that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from cattle and have adapted well to the extreme conditions on the plateau. However, the mechanisms used by these animals for the regulation of gene expression at high altitude are not fully understood.

RESULTS

Here, we sequenced nine lung transcriptomes of yaks at altitudes of 3400, 4200 and 5000 m, and low-altitude Zaosheng cattle, which is a closely related species, served as controls. The analysis identified 21,764 mRNAs, 1377 circRNAs and 1209 miRNAs. By comparing yaks and cattle, 4975 mRNAs, 252 circRNAs and 75 miRNAs were identified differentially expressed. By comparing yaks at different altitudes, we identified 756 mRNAs, 64 circRNAs and 83 miRNAs that were differentially expressed (fold change ≥2 and P-value < 0.05). The pathways enriched in the mRNAs, circRNAs and miRNAs identified from the comparison of yaks and cattle were mainly associated with metabolism, including 'glycosaminoglycan degradation', 'pentose and glucuronate interconversions' and 'flavone and flavonol biosynthesis', and the mRNAs, circRNAs and miRNAs identified from the comparison of yaks at different altitude gradients were significantly enriched in metabolic pathways and immune and genetic information processing pathways. The core RNAs were identified from the mRNA-miRNA-circRNA networks constructed using the predominant differentially expressed RNAs. The core genes specific to the difference between yaks and cattle were associated with the endoplasmic reticulum and fat deposition, but those identified from the comparison among yaks at different altitude gradients were associated with maintenance of the normal biological functions of cells.

CONCLUSIONS

This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in yaks and might contribute to improvements in the understanding and prevention of hypoxia-related diseases.

Fine-Scale Genetic Structure and Natural Selection Signatures of Southwestern Hans Inferred From Patterns of Genome-Wide Allele, Haplotype, and Haplogroup Lineages

The evolutionary and admixture history of Han Chinese have been widely discussed via traditional autosomal and uniparental genetic markers [e.g., short tandem repeats, low-density single nucleotide polymorphisms). However, their fine-scale genetic landscapes (admixture scenarios and natural selection signatures) based on the high-density allele/haplotype sharing patterns have not been deeply characterized. Here, we collected and generated genome-wide data of 50 Han Chinese individuals from four populations in Guizhou Province, one of the most ethnolinguistically diverse regions, and merged it with over 3,000 publicly available modern and ancient Eurasians to describe the genetic origin and population admixture history of Guizhou Hans and their neighbors. PCA and ADMIXTURE results showed that the studied four populations were homogeneous and grouped closely to central East Asians. Genetic homogeneity within Guizhou populations was further confirmed via the observed strong genetic affinity with inland Hmong-Mien people through the observed genetic clade in Fst and outgroup f₃/f₄-statistics. qpGraph-based phylogenies and f₄-based demographic models illuminated that Guizhou Hans were well fitted via the admixture of ancient Yellow River Millet farmers related to Lajia people and southern Yangtze River farmers related to Hanben people. Further ChromoPainter-based chromosome painting profiles and GLOBETROTTER-based admixture signatures confirmed the two best source matches for southwestern Hans, respectively, from northern Shaanxi Hans and southern indigenes with variable mixture proportions in the historical period. Further three-way admixture models revealed larger genetic contributions from coastal southern East Asians into Guizhou Hans compared with the proposed inland ancient source from mainland Southeast Asia. We also identified candidate loci (e.g., MTUS2, NOTCH4, EDAR, ADH1B, and ABCG2) with strong natural selection signatures in Guizhou Hans via iHS, nSL, and ihh, which were associated with the susceptibility of the multiple complex diseases, morphology formation, alcohol and lipid metabolism. Generally, we provided a case and ideal strategy to reconstruct the detailed demographic evolutionary history of Guizhou Hans, which provided new insights into the fine-scale genomic formation of one ethnolinguistically specific targeted population from the comprehensive perspectives of the shared unlinked alleles, linked haplotypes, and paternal and maternal lineages.

Peopling History of the Tibetan Plateau and Multiple Waves of Admixture of Tibetans Inferred From Both Ancient and Modern Genome-Wide Data

Archeologically attested human occupation on the Tibetan Plateau (TP) can be traced back to 160 thousand years ago (kya) via the archaic Xiahe people and 30∼40 kya via the Nwya Devu anatomically modern human. However, the history of the Tibetan populations and their migration inferred from the ancient and modern DNA remains unclear. Here, we performed the first ancient and modern genomic meta-analysis among 3,017 Paleolithic to present-day Eastern Eurasian genomes (2,444 modern individuals from 183 populations and 573 ancient individuals). We identified a close genetic connection between the ancient-modern highland Tibetans and lowland island/coastal Neolithic Northern East Asians (NEA). This observed genetic affinity reflected the primary ancestry of high-altitude Tibeto-Burman speakers originated from the Neolithic farming populations in the Yellow River Basin. The identified pattern was consistent with the proposed common north-China origin hypothesis of the Sino-Tibetan languages and dispersal patterns of the northern millet farmers. We also observed the genetic differentiation between the highlanders and lowland NEAs. The former harbored more deeply diverged Hoabinhian/Onge-related ancestry and the latter possessed more Neolithic southern East Asian (SEA) or Siberian-related ancestry. Our reconstructed qpAdm and qpGraph models suggested the co-existence of Paleolithic and Neolithic ancestries in the Neolithic to modern East Asian highlanders. Additionally, we found that Tibetans from Ü-Tsang/Ando/Kham regions showed a strong population stratification consistent with their cultural background and geographic terrain. Ü-Tsang Tibetans possessed a stronger Chokhopani-affinity, Ando Tibetans had more Western Eurasian related ancestry and Kham Tibetans harbored greater Neolithic southern EA ancestry. Generally, ancient and modern genomes documented multiple waves of human migrations in the TP's past. The first layer of local hunter-gatherers mixed with incoming millet farmers and arose the Chokhopani-associated Proto-Tibetan-Burman highlanders, which further respectively mixed with additional genetic contributors from the western Eurasian Steppe, Yellow River and Yangtze River and finally gave rise to the modern Ando, Ü-Tsang and Kham Tibetans.

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.

The History and Driving Force for Prehistoric Human Expansion Upward to the Hinterland of the Tibetan Plateau Post–Last Glacial Maximum

The timing and motivation of prehistoric human expansion into the hinterland of the Tibetan Plateau (TP) is a widely debated scientific issue. Recent archaeological studies have brought forward predictions of the earliest human occupation of the TP to the late–Middle Pleistocene. However, massive human occupation of the TP did not appear until the termination of the Last Glacial Maximum (LGM). The spatio-temporal distribution of prehistoric hunter-gatherers on the TP varies significantly before the permanent occupation after 3600 BP (before present). Here, we report on environmental-archaeological evidence from the Canxionggashuo (CXGS) site in Yushu Prefecture, which provides information that is key to understanding the dynamics of post-LGM human occupation on the TP. Radiocarbon dating has revealed two occupation periods of the CXGS site at 8600–7100 cal (calibrated years) BP and 2400–2100 cal BP. The charcoal concentration in cultural layers correlates well with paleo–human activities. Hunter-gatherers expanded westwards from the northeastern margin of the TP to the hinterland of the TP during the warming period of the early–middle Holocene (~11,500–6000 BP). However, these groups retreated during the middle–late Holocene (~6000–3600 BP) under a cooling-drying climate. Prehistoric humans finally occupied the hinterland of the TP permanently after 3600 BP, with an enhanced cold-adaptive lifestyle, although the climate was still deteriorating.

Characterization of structural variation in Tibetans reveals new evidence of high-altitude adaptation and introgression

BACKGROUND

Structural variation (SV) acts as an essential mutational force shaping the evolution and function of the human genome. However, few studies have examined the role of SVs in high-altitude adaptation and little is known of adaptive introgressed SVs in Tibetans so far.

RESULTS

Here, we generate a comprehensive catalog of SVs in a Chinese Tibetan (n = 15) and Han (n = 10) population using nanopore sequencing technology. Among a total of 38,216 unique SVs in the catalog, 27% are sequence-resolved for the first time. We systematically assess the distribution of these SVs across repeat sequences and functional genomic regions. Through genotyping in additional 276 genomes, we identify 69 Tibetan-Han stratified SVs and 80 candidate adaptive genes. We also discover a few adaptive introgressed SV candidates and provide evidence for a deletion of 335 base pairs at 1p36.32.

CONCLUSIONS

Overall, our results highlight the important role of SVs in the evolutionary processes of Tibetans' adaptation to the Qinghai-Tibet Plateau and provide a valuable resource for future high-altitude adaptation studies.

Local adaptation of Mycobacterium tuberculosis on the Tibetan Plateau

During its global dispersal, Mycobacterium tuberculosis (Mtb) has encountered varied geographic environments and host populations. Although local adaptation seems to be a plausible model for describing long-term host-pathogen interactions, genetic evidence for this model is lacking. Here, we analyzed 576 whole-genome sequences of Mtb strains sampled from different regions of high-altitude Tibet. Our results show that, after sequential introduction of a few ancestral strains, the Tibetan Mtb population diversified locally while maintaining strict separation from the Mtb populations on the lower altitude plain regions of China. The current population structure and estimated past population dynamics suggest that the modern Beijing sublineage strains, which expanded over most of China and other global regions, did not show an expansion advantage in Tibet. The mutations in the Tibetan strains showed a higher proportion of A > G/T > C transitions than strains from the plain regions, and genes encoding DNA repair enzymes showed evidence of positive selection. Moreover, the long-term Tibetan exclusive selection for truncating mutations in the thiol-oxidoreductase encoding sseA gene suggests that Mtb was subjected to local selective pressures associated with oxidative stress. Collectively, the population genomics of Mtb strains in the relatively isolated population of Tibet provides genetic evidence that Mtb has adapted to local environments.

Combined Low-/High-Density Modern and Ancient Genome-Wide Data Document Genomic Admixture History of High-Altitude East Asians

The Tibetan Plateau (TP) is considered to be one of the last terrestrial environments conquered by the anatomically modern human. Understanding of the genetic background of highland Tibetans plays a pivotal role in archeology, anthropology, genetics, and forensic investigations. Here, we genotyped 22 forensic genetic markers in 1,089 Tibetans residing in Nagqu Prefecture and collected 1,233,013 single nucleotide polymorphisms (SNPs) in the highland East Asians (Sherpa and Tibetan) from the Simons Genome Diversity Project and ancient Tibetans from Nepal and Neolithic farmers from northeastern Qinghai-Tibetan Plateau from public databases. We subsequently merged our two datasets with other worldwide reference populations or eastern ancient Eurasians to gain new insights into the genetic diversity, population movements, and admixtures of high-altitude East Asians via comprehensive population genetic statistical tools [principal component analysis (PCA), multidimensional scaling plot (MDS), STRUCTURE/ADMIXTURE, f3 , f4 , qpWave/qpAdm, and qpGraph]. Besides, we also explored their forensic characteristics and extended the Chinese National Database based on STR data. We identified 231 alleles with the corresponding allele frequencies spanning from 0.0005 to 0.5624 in the forensic low-density dataset, in which the combined powers of discrimination and the probability of exclusion were 1-1.22E-24 and 0.999999998, respectively. Additionally, comprehensive population comparisons in our low-density data among 57 worldwide populations via the Nei's genetic distance, PCA, MDS, NJ tree, and STRUCTURE analysis indicated that the highland Tibeto-Burman speakers kept the close genetic relationship with ethnically close populations. Findings from the 1240K high-density dataset not only confirmed the close genetic connection between modern Highlanders, Nepal ancients (Samdzong, Mebrak, and Chokhopani), and the upper Yellow River Qijia people, suggesting the northeastern edge of the TP served as a geographical corridor for ancient population migrations and interactions between highland and lowland regions, but also evidenced that late Neolithic farmers permanently colonized into the TP by adopting cold-tolerant barley agriculture that was mediated via the acculturation of idea via the millet farmer and not via the movement of barley agriculturalist as no obvious western Eurasian admixture signals were identified in our analyzed modern and ancient populations. Besides, results from the qpAdm-based admixture proportion estimation and qpGraph-based phylogenetic relationship reconstruction consistently demonstrated that all ancient and modern highland East Asians harbored and shared the deeply diverged Onge/Hoabinhian-related eastern Eurasian lineage, suggesting a common Paleolithic genetic legacy existed in high-altitude East Asians as the first layer of their gene pool.

Shared paternal ancestry of Han, Tai-Kadai-speaking, and Austronesian-speaking populations as revealed by the high resolution phylogeny of O1a-M119 and distribution of its sub-lineages within China

OBJECTIVES

The aim of this research was to explore the origin, diversification, and demographic history of O1a-M119 over the past 10,000 years, as well as its role during the formation of East Asian and Southeast Asian populations, particularly the Han, Tai-Kadai-speaking, and Austronesian-speaking populations.

MATERIALS AND METHODS

Y-chromosome sequences (n = 141) of the O1a-M119 lineage, including 17 newly generated in this study, were used to reconstruct a revised phylogenetic tree with age estimates, and identify sub-lineages. The geographic distribution of 12 O1a-M119 sub-lineages was summarized, based on 7325 O1a-M119 individuals identified among 60,009 Chinese males.

RESULTS

A revised phylogenetic tree, age estimation, and distribution maps indicated continuous expansion of haplogroup O1a-M119 over the past 10,000 years, and differences in demographic history across geographic regions. We propose several sub-lineages of O1a-M119 as founding paternal lineages of Han, Tai-Kadai-speaking, and Austronesian-speaking populations. The sharing of several young O1a-M119 sub-lineages with expansion times less than 6000 years between these three population groups supports a partial common ancestry for them in the Neolithic Age; however, the paternal genetic divergence pattern is much more complex than previous hypotheses based on ethnology, archeology, and linguistics.

DISCUSSION

Our analyses contribute to a better understanding of the demographic history of O1a-M119 sub-lineages over the past 10,000 years during the emergence of Han, Austronesians, Tai-Kadai-speaking populations. The data described in this study will assist in understanding of the history of Han, Tai-Kadai-speaking, and Austronesian-speaking populations from ethnology, archeology, and linguistic perspectives in the future.

Dated phylogeny suggests early Neolithic origin of Sino-Tibetan languages

An accurate reconstruction of Sino-Tibetan language evolution would greatly advance our understanding of East Asian population history. Two recent phylogenetic studies attempted to do so but several of their conclusions are different from each other. Here we reconstruct the phylogeny of the Sino-Tibetan language family, using Bayesian computational methods applied to a larger and linguistically more diverse sample. Our results confirm previous work in finding that the ancestral Sino-Tibetans first split into Sinitic and Tibeto-Burman clades, and support the existence of key internal relationships. But we find that the initial divergence of this group occurred earlier than previously suggested, at approximately 8000 years before the present, coinciding with the onset of millet-based agriculture and significant environmental changes in the Yellow River region. Our findings illustrate that key aspects of phylogenetic history can be replicated in this complex language family, and calls for a more nuanced understanding of the first Sino-Tibetan speakers in relation to the "early farming dispersal" theory of language evolution.

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.

Chromatin accessibility landscape and regulatory network of high-altitude hypoxia adaptation

High-altitude adaptation of Tibetans represents a remarkable case of natural selection during recent human evolution. Previous genome-wide scans found many non-coding variants under selection, suggesting a pressing need to understand the functional role of non-coding regulatory elements (REs). Here, we generate time courses of paired ATAC-seq and RNA-seq data on cultured HUVECs under hypoxic and normoxic conditions. We further develop a variant interpretation methodology (vPECA) to identify active selected REs (ASREs) and associated regulatory network. We discover three causal SNPs of EPAS1, the key adaptive gene for Tibetans. These SNPs decrease the accessibility of ASREs with weakened binding strength of relevant TFs, and cooperatively down-regulate EPAS1 expression. We further construct the downstream network of EPAS1, elucidating its roles in hypoxic response and angiogenesis. Collectively, we provide a systematic approach to interpret phenotype-associated noncoding variants in proper cell types and relevant dynamic conditions, to model their impact on gene regulation.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Haptoglobin polymorphisms in Latin American populations

Several genetic polymorphisms of the haptoglobin gene (HP) or haptoglobin-related gene (HPR) were reported to show a population-specific distribution and to be associated with not only serum haptoglobin (HP) but also cholesterol levels. For such association studies, it is important to know the distribution of polymorphisms or their haplotypes in the populations concerned. However, no comprehensive genetic studies have explored this in Latin Americans, and not every human variation or genotype is available in a database. In this study, we determined the genotypes of common HP (HP¹ and HP²), HP^del, rs5471, rs5472, and rs2000999 in several Latin American populations. Haplotypes of rs5472-common HP-rs2000999 polymorphisms were estimated. We did not encounter any HP^del, and the frequencies of rs5471 A, rs5472 A, HP¹, and rs2000999 G were higher than their counterpart alleles in studied populations. All of the alleles with higher frequency in the Latin Americans are associated with higher serum HP and lower cholesterol levels. Both A-1-G (probably HP^1S) and G-1-G (probably HP^1F) haplotypes were higher in Latin American populations than those in other geographic regions. In addition, the genetic influx from populations of other continents into Peruvians seems to be relatively lower than that into other Latin Americans.

Cross-Species Insights Into Genomic Adaptations to Hypoxia

Over millions of years, vertebrate species populated vast environments spanning the globe. Among the most challenging habitats encountered were those with limited availability of oxygen, yet many animal and human populations inhabit and perform life cycle functions and/or daily activities in varying degrees of hypoxia today. Of particular interest are species that inhabit high-altitude niches, which experience chronic hypobaric hypoxia throughout their lives. Physiological and molecular aspects of adaptation to hypoxia have long been the focus of high-altitude populations and, within the past decade, genomic information has become increasingly accessible. These data provide an opportunity to search for common genetic signatures of selection across uniquely informative populations and thereby augment our understanding of the mechanisms underlying adaptations to hypoxia. In this review, we synthesize the available genomic findings across hypoxia-tolerant species to provide a comprehensive view of putatively hypoxia-adaptive genes and pathways. In many cases, adaptive signatures across species converge on the same genetic pathways or on genes themselves [i.e., the hypoxia inducible factor (HIF) pathway). However, specific variants thought to underlie function are distinct between species and populations, and, in most cases, the precise functional role of these genomic differences remains unknown. Efforts to standardize these findings and explore relationships between genotype and phenotype will provide important clues into the evolutionary and mechanistic bases of physiological adaptations to environmental hypoxia.

Why and when was lactase persistence selected for? Insights from Central Asian herders and ancient DNA

The genetic adaptation of humans to the consumption of milk from dairying animals is one of the most emblematic cases of recent human evolution. While the phenotypic change under selection, lactase persistence (LP), is known, the evolutionary advantage conferred to persistent individuals remains obscure. One informative but underappreciated observation is that not all populations whose ancestors had access to milk genetically adapted to become lactase persistent. Indeed, Central Asian herders are mostly lactase nonpersistent, despite their significant dietary reliance on dairy products. Investigating the temporal dynamic of the -13.910:C>T Eurasian mutation associated with LP, we found that, after its emergence in Ukraine 5,960 before present (BP), the T allele spread between 4,000 BP and 3,500 BP throughout Eurasia, from Spain to Kazakhstan. The timing and geographical progression of the mutation coincides well with the migration of steppe populations across and outside of Europe. After 3,000 BP, the mutation strongly increased in frequency in Europe, but not in Asia. We propose that Central Asian herders have adapted to milk consumption culturally, by fermentation, and/or by colonic adaptation, rather than genetically. Given the possibility of a nongenetic adaptation to avoid intestinal symptoms when consuming dairy products, the puzzle then becomes this: why has LP been selected for at all?