Neolithic millet farmers contributed to the permanent settlement of the Tibetan Plateau by adopting barley agriculture

Natural selection shaped the protective effect of the mtDNA lineage against obesity in Han Chinese populations

Mitochondria play a key role in lipid metabolism, and mitochondrial DNA (mtDNA) mutations are thus considered to affect obesity susceptibility by altering oxidative phosphorylation and mitochondrial function. In this study, we investigated mtDNA variants that may affect obesity risk in 2877 Han Chinese individuals from three independent populations. The association analysis of 16 basal mtDNA haplogroups with body mass index, waist circumference, and waist-to-hip ratio revealed that only haplogroup M7 was significantly negatively correlated with all three adiposity-related anthropometric traits in the overall cohort, verified by the analysis of a single population, i.e., the Zhengzhou population. Furthermore, subhaplogroup analysis suggested that M7b1a1 was the most likely haplogroup associated with a decreased obesity risk, and the variation T12811C (causing Y159H in ND5) harbored in M7b1a1 may be the most likely candidate for altering the mitochondrial function. Specifically, we found that proportionally more nonsynonymous mutations accumulated in M7b1a1 carriers, indicating that M7b1a1 was either under positive selection or subject to a relaxation of selective constraints. We also found that nuclear variants, especially in DACT2 and PIEZO1, may functionally interact with M7b1a1.

Distinguished biological adaptation architecture aggravated population differentiation of Tibeto-Burman-speaking people

Tibeto-Burman (TB) people have endeavored to adapt to the hypoxic, cold, and high-UV high-altitude environments in the Tibetan Plateau and complex disease exposures in lowland rainforests since the late Paleolithic period. However, the full landscape of genetic history and biological adaptation of geographically diverse TB-speaking people, as well as their interaction mechanism, remain unknown. Here, we generate a whole-genome meta-database of 500 individuals from 39 TB-speaking populations and present a comprehensive landscape of genetic diversity, admixture history, and differentiated adaptative features of geographically different TB-speaking people. We identify genetic differentiation related to geography and language among TB-speaking people, consistent with their differentiated admixture process with incoming or indigenous ancestral source populations. A robust genetic connection between the Tibetan-Yi corridor and the ancient Yellow River people supports their Northern China origin hypothesis. We finally report substructure-related differentiated biological adaptative signatures between highland Tibetans and Loloish speakers. Adaptative signatures associated with the physical pigmentation (EDAR and SLC24A5) and metabolism (ALDH9A1) are identified in Loloish people, which differed from the high-altitude adaptative genetic architecture in Tibetan. TB-related genomic resources provide new insights into the genetic basis of biological adaptation and better reference for the anthropologically informed sampling design in biomedical and genomic cohort research.

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.

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.

Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau

BACKGROUND

Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products.

RESULTS

We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe.

CONCLUSIONS

Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

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.

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.

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.

A Late Pleistocene human genome from Southwest China

Southern East Asia is the dispersal center regarding the prehistoric settlement and migrations of modern humans in Asia-Pacific regions. However, the settlement pattern and population structure of paleolithic humans in this region remain elusive, and ancient DNA can provide direct information. Here, we sequenced the genome of a Late Pleistocene hominin (MZR), dated ∼14.0 thousand years ago from Red Deer Cave located in Southwest China, which was previously reported possessing mosaic features of modern and archaic hominins. MZR is the first Late Pleistocene genome from southern East Asia. Our results indicate that MZR is a modern human who represents an early diversified lineage in East Asia. The mtDNA of MZR belongs to an extinct basal lineage of the M9 haplogroup, reflecting a rich matrilineal diversity in southern East Asia during the Late Pleistocene. Combined with the published data, we detected clear genetic stratification in ancient southern populations of East/Southeast Asia and some degree of south-versus-north divergency during the Late Pleistocene, and MZR was identified as a southern East Asian who exhibits genetic continuity to present day populations. Markedly, MZR is linked deeply to the East Asian ancestry that contributed to First Americans.

Vegetation History and Survival Patterns of the Earliest Village on the Qinghai-Tibetan Plateau

The upper Yellow River valley in the northeastern Qinghai-Tibetan Plateau (QTP) is an important corridor for prehistoric migration to the hinterland plateau. However, most studies have focused on the Neolithic Age, with limited evidence for earlier periods. The Shalongka (SLK) site on the northeastern QTP spans the Epipaleolithic to Bronze Age and contains cultural deposits, so provides a good basis for unraveling the evolutionary history of the human-land relationship. In this study, we sampled the 420-cm-thick section T1406E at the SLK site and undertook lithologic stratigraphic description and analysis of grain size, redness, magnetic susceptibility, geochemical elements, pollen and charcoal. Dating control was provided by accelerated mass spectrometry 14C and optically stimulated luminescence methods. Results show that SLK site was affected by the local fluvial sedimentary environment. The absolute dating results of the SLK site have revealed that humans occupied the site during the Epipaleolithic (8.5-7.3 cal ka BP), Yangshao culture (5.9-5.1 ka) and Qijia Culture (4.1-3.9 cal ka BP). Pollen analysis showed that the humans lived in a landscape that was predominated by forest-steppe. Consolidating with multidisciplinary evidence, we learned that Epipaleolithic sites were occupied by microlithic hunter-gatherers and comprised by relatively fixed seasonal central campsites, and their mobility was significantly decreased from the early to late period. Subsequently, farmers of the Yangshao culture migrated from the low elevation (Chinese Loess Plateau) to the upper Yellow River valleys on the QTP and founded the earliest settlement villages (~5.9 ka) on the QTP. People of the Qijia culture adopted diversified survival strategies under the settled lifestyle. In all, we infered that SLK site may play an important role in the communication and integration between different people and cultures.

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.

The genomic history of southwestern Chinese populations demonstrated massive population migration and admixture among proto-Hmong-Mien speakers and incoming migrants

Southwest China was the crossroad for the initial settler people of East Asia, which shows the highest diversity in languages and genetics. This region played a significant role in the formation of the genetic makeup of the proto-Hmong-Mien-speaking people and in the north-to-south human expansion during the Neolithic-to-historic transformation. Their genetic history covering migration events and the admixture processes still needs to be further explored. Therefore, in the current study, we have generated genome-wide data from three genomic aspects covering autosomal, mitochondrial and Y-chromosomal regions in 260 Hmong-Mien, Tibeto-Burman, and Sinitic people from 29 different southwestern Chinese groups, and further analyzed them with 2676 published modern and ancient Eurasian genomes. Here, we have noticed a new southwestern East Asian genetic cline composed of the Hmong-Mien-specific ancestry enriched in modern Hmong and Pathen. This newly identified southern inland East Asian lineage contributed to a great extent of the gene pool in the modern southern East Asians. We also have observed genetic substructure among Hmong-Mien-speaking populations. The southern Hmong-Mien-speaking people showed more genetic affinity with modern Tai-Kadai/Austroasiatic people, while the northern Hmong-Mien speakers expressed a closer genetic connection with the Neolithic-to-modern northern East Asians. Moreover, southwestern Sinitic populations had a strong genomic affinity with the adjacent Hmong-Mien-speaking populations and the lowlander Tibeto-Burman-speaking populations, which suggested the large-scale genetic admixture occurred between them. Allele-sharing-based qpAdm/qpGraph results further confirmed that all included southwestern Chinese populations could be modeled as a mixed result of the major ancestry component from the northern millet farmers in the Yellow River basin and the minor ancestry component from the southern rice farmers in the Yangtze River basin. Usually, this newly identified Hmong-Mien-associated southern East Asian ancestry could improve our understanding of the full-scale genetic landscape of the evolutionary and admixture history of southwestern East Asians. Further ancient genomic studies from southeastern China are required to shed deeper light on our established phylogeny context.

Genomic history and forensic characteristics of Sherpa highlanders on the Tibetan Plateau inferred from high-resolution InDel panel and genome-wide SNPs

Sherpa people, one of the high-altitude hypoxic adaptive populations, mainly reside in Nepal and the southern Tibet Autonomous Region. The genetic origin and detailed evolutionary profiles of Sherpas remain to be further explored and comprehensively characterized. Here we analyzed the newly-generated InDel genotype data from 628 Dingjie Sherpas by merging with 4222 worldwide InDel profiles and collected genome-wide SNP data (approximately 600K SNPs) from 1612 individuals in 191 modern and ancient populations to explore and reconstruct the fine-scale genetic structure of Sherpas and their relationships with nearby modern and ancient East Asians based on the shared alleles and haplotypes. The forensic parameters of 57 autosomal InDels (A-InDels) included in our used new-generation InDel amplification system showed that this focused InDel panel is informative and polymorphic in Dingjie Sherpas, suggesting that it can be used as the supplementary tool for forensic personal identification and parentage testing in Dingjie Sherpas. Descriptive findings from the PCA, ADMIXTURE, and TreeMix-based phylogenies suggested that studied Nepal Sherpas showed excess allele sharing with neighboring Tibeto-Burman Tibetans. Furthermore, patterns of allele sharing in f-statistics demonstrated that Nepal Sherpas had a different evolutionary history compared with their neighbors from Nepal (Newar and Gurung) but showed genetic similarity with 2700-year-old Chokhopani and modern Tibet Tibetans. QpAdm/qpGraph-based admixture sources and models further showed that Sherpas, core Tibetans, and Chokhopani formed one clade, which could be fitted as having the main ancestry from late Neolithic Qijia millet farmers and other deep ancestries from early Asians. Chromosome painting profiles and shared IBD fragments inferred from fineSTRUCTURE and ChromoPainter not only confirmed the abovementioned genomic affinity patterns but also revealed the fine-scale genetic microstructures among Sino-Tibetan speakers. Finally, natural-selection signals revealed via iHS, nSL and iHH12 showed natural selection signatures associated with disease susceptibility in Sherpas. Generally, we provided the comprehensive landscape of admixture and evolutionary history of Sherpa people based on the shared alleles and haplotypes from the InDel-based genotype data and high-density genome-wide SNP data. The more detailed genetic landscape of Sherpa people should be further confirmed and characterized via ancient genomes or single-molecule real-time sequencing technology.

Mitogenomics of modern Mongolic-speaking populations

Here, we present a comprehensive data set of 489 complete mitogenomes (211 of which are new) from four Mongolic-speaking populations (Mongols, Barghuts, Khamnigans, and Buryats) to investigate their matrilineal genetic structure, ancestry and relationship with other ethnic groups. We show that along with very high levels of genetic diversity and lack of genetic differentiation, Mongolic-speaking populations exhibit strong genetic resemblance to East Asian populations of Chinese, Japanese, and Uyghurs. Phylogeographic analysis of complete mitogenomes reveals the presence of different components in the gene pools of modern Mongolic-speaking populations-the main East Eurasian component is represented by mtDNA lineages of East Asian, Siberian and autochthonous (the Baikal region/Mongolian) ancestry, whereas the less pronounced West Eurasian component can be ascribed to Europe and West Asia/Caucasus. We also observed that up to one third of the mtDNA subhaplogroups identified in Mongolic-speaking populations can be considered as Mongolic-specific with the coalescence age of most of them not exceeding 1.7 kya. This coincides well with the population size growth which started around 1.1 kya and is detectable only in the Bayesian Skyline Plot constructed based on Mongolic-specific mitogenomes. Our data suggest that the genetic structure established during the Mongol empire is still retained in present-day Mongolic-speaking populations.

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.

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.

A Mitochondrial DNA Variant Elevates the Risk of Gallstone Disease by Altering Mitochondrial Function

BACKGROUND AND AIMS

Gallstone disease (cholelithiasis) is a cholesterol-related metabolic disorders with strong familial predisposition. Mitochondrial DNA (mtDNA) variants accumulated during human evolution are associated with some metabolic disorders related to modified mitochondrial function. The mechanistic links between mtDNA variants and gallstone formation need further exploration.

METHODS

In this study, we explored the possible associations of mtDNA variants with gallstone disease by comparing 104 probands and 300 controls in a Chinese population. We constructed corresponding cybrids using trans-mitochondrial technology to investigate the underlying mechanisms of these associations. Mitochondrial respiratory chain complex activity and function and cholesterol metabolism were assessed in the trans-mitochondrial cell models.

RESULTS

Here, we found a significant association of mtDNA 827A>G with an increased risk of familial gallstone disease in a Chinese population (odds ratio [OR]: 4.5, 95% confidence interval [CI]: 2.1-9.4, P=1.2×10-4). Compared with 827A cybrids (haplogroups B4a and B4c), 827G cybrids (haplogroups B4b and B4d) had impaired mitochondrial respiratory chain complex activity and function and activated JNK and AMPK signaling pathways. Additionally, the 827G cybrids showed disturbances in cholesterol transport and accelerated development of gallstones. Specifically, cholesterol transport through the transporter ABCG5/8 was increased via activation of the AMPK signaling pathway in 827G cybrids.

CONCLUSIONS

Our findings reveal that mtDNA 827A>G induces aberrant mitochondrial function and abnormal cholesterol transport, resulting in increased occurrence of gallstones. The results provide an important biological basis for the clinical diagnosis and prevention of gallstone disease 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.

Ancient genomes reveal tropical bovid species in the Tibetan Plateau contributed to the prevalence of hunting game until the late Neolithic

Local wild bovids have been determined to be important prey on the northeastern Tibetan Plateau (NETP), where hunting game was a major subsistence strategy until the late Neolithic, when farming lifestyles dominated in the neighboring Loess Plateau. However, the species affiliation and population ecology of these prehistoric wild bovids in the prehistoric NETP remain unknown. Ancient DNA (aDNA) analysis is highly informative in decoding this puzzle. Here, we applied aDNA analysis to fragmented bovid and rhinoceros specimens dating ∼5,200 y B.P. from the Neolithic site of Shannashuzha located in the marginal area of the NETP. Utilizing both whole genomes and mitochondrial DNA, our results demonstrate that the range of the present-day tropical gaur (Bos gaurus) extended as far north as the margins of the NETP during the late Neolithic from ∼29°N to ∼34°N. Furthermore, comparative analysis with zooarchaeological and paleoclimatic evidence indicated that a high summer temperature in the late Neolithic might have facilitated the northward expansion of tropical animals (at least gaur and Sumatran-like rhinoceros) to the NETP. This enriched the diversity of wildlife, thus providing abundant hunting resources for humans and facilitating the exploration of the Tibetan Plateau as one of the last habitats for hunting game in East Asia.

Population genetic analysis of Shaanxi male Han Chinese population reveals genetic differentiation and homogenization of East Asians

BACKGROUND

Shaanxi province, located in the upper Yellow River, has been evidenced as the geographic origin of Chinese civilization, Sino-Tibetan-speaking language, and foxtail or broomcorn millet farmers via the linguistic phylogenetic spectrum, archeological documents, and genetic evidence. Nowadays, Han Chinese is the dominant population in this area. The formation process of modern Shaanxi Han population reconstructed via the ancient DNA is on the way, however, the patterns of genetic relationships of modern Shaanxi Han, allele frequency distributions of high mutated short tandem repeats (STRs) and corresponding forensic parameters are remained to be explored.

METHODS

Here, we successfully genotyped 23 autosomal STRs in 630 unrelated Shaanxi male Han individuals using the recently updated Huaxia Platinum PCR amplification system. Forensic allele frequency and parameters of all autosomal STRs were assessed. And comprehensive population genetic structure was explored via various typical statistical technologies.

RESULTS

Population genetic analysis based on the raw-genotype dataset among 15,803 Eurasian individuals and frequency datasets among 56 populations generally illustrated that linguistic stratification is significantly associated with the genetic substructure of the East Asian population. Principal component analysis, multidimensional scaling plots and phylogenetic tree further demonstrated that Shaanxi Han has a close genetic relationship with geographically close Shanxi Han, and showed that Han Chinese is a homogeneous population during the historic and recent admixture from the STR variations. Except for Sinitic-speaking populations, Shaanxi Han harbored more alleles sharing with Tibeto-Burman-speaking populations than with other reference populations. Focused on the allele frequency correlation and forensic parameters, all loci are in accordance with the minimum requirements of HWE and LD. The observed combined probability of discrimination of 8.2201E-28 and the cumulative power of exclusion of 0.9999999995 in Shaanxi Han demonstrated that the studied STR loci are informative and polymorphic, and this system can be used as a powerful routine forensic tool in personal identification and parentage testing.

CONCLUSION

Both the geographical and linguistic divisions have shaped the genetic structure of modern East Asian. And more forensic reference data should be obtained for ethnically, culturally, geographically and linguistically different populations for better routine forensic practice and population genetic studies.

Massively parallel sequencing of mitogenome sequences reveals the forensic features and maternal diversity of tai-kadai-speaking hlai islanders

As a single maternally inherited locus, human mitochondrial DNA (mtDNA) is geographically arranged and plays a key role in forensic applications. Hlai population has been evidenced as the most typical and unmixed representative of the Tai-Kadai-speaking populations via genome-wide analyses. However, forensic features and maternal diversity of the complete mitogenomes in this Tai-Kadai ancestrally related population are scarce. Thus, we sequenced the complete mitogenomes in 127 Hainan Hlais and found 109 distinct haplotypes belonging to 43 terminal haplogroups resulting in the haplotype diversity of 0.9970. Our results of comprehensive population comparisons showed that Hlai islanders had a close genetic affinity with Tai-Kadai-speaking populations from Southeast Asia, which is consistent with the back-migration of Chinese Neolithic farmers into this region via the inland route. Besides, maternally genetic evidence further revealed a close genetic relationship between Tai-Kadai-speaking and Austronesian-speaking populations when only East Asian dataset was considered, which is consistent with the common origin from Yangtze rice farmers and then spread southward along the inland and coastal routes, respectively. In the reconstructed phylogenetic tree and median-joining networks, the vast majority of Hlais were clustered in exclusive clades, which demonstrated that Hlai people probably had undergone founder effect or genetic bottleneck in their history, and remained genetically isolated for a long time. Collectively, Hainan Hlai did not exhibit detectable maternal gene flow from surrounding or incoming populations. Mitogenome information generated in this study is a contribution in mitigating the underrepresentation of Chinese data in forensic mitogenetics and will assist geography-, metapopulation-, as well as phylogeny-based queries.