Reconstruction of Y-chromosome phylogeny reveals two neolithic expansions of Tibeto-Burman populations

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.

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.

Forager-farmer transition at the crossroads of East and Southeast Asia 4900 years ago

The southward expansion of East Asian farmers profoundly influenced the social evolution of Southeast Asia by introducing cereal agriculture. However, the timing and routes of cereal expansion in key regions are unclear due to limited empirical evidence. Here we report macrofossil, microfossil, multiple isotopic (C/N/Sr/O) and paleoproteomic data directly from radiocarbon-dated human samples, which were unearthed from a site in Xingyi in central Yunnan and which date between 7000 and 3300 a BP. Dietary isotopes reveal the earliest arrival of millet ca. 4900 a BP, and greater reliance on plant and animal agriculture was indicated between 3800 and 3300 a BP. The dietary differences between hunter-gatherer and agricultural groups are also evident in the metabolic and immune system proteins analysed from their skeletal remains. The results of paleoproteomic analysis indicate that humans had divergent biological adaptations, with and without farming. The combined application of isotopes, archaeobotanical data and proteomics provides a new approach to documenting dietary and health changes across major subsistence transitions.

Ancient genomes reveal millet farming-related demic diffusion from the Yellow River into southwest China

The study of southwest China is vital for understanding the dispersal and development of farming because of the coexistence of millet and rice in this region since the Neolithic period.1,2 However, the process of the Neolithic transition in southwest China is largely unknown, mainly due to the lack of ancient DNA from the Neolithic period. Here, we report genome-wide data from 11 human samples from the Gaoshan and Haimenkou sites with mixed farming of millet and rice dating to between 4,500 and 3,000 years before present in southwest China. The two ancient groups derived approximately 90% of their ancestry from the Neolithic Yellow River farmers, suggesting a demic diffusion of millet farming to southwest China. We inferred their remaining ancestry to be derived from a Hòabìnhian-related hunter-gatherer lineage. We did not detect rice farmer-related ancestry in the two ancient groups, which indicates that they likely adopted rice farming without genetic assimilation. We, however, observed rice farmer-related ancestry in the formation of some present-day Tibeto-Burman populations. Our results suggested the occurrence of both demic and cultural diffusion in the development of Neolithic mixed farming in some parts of southwest China.

Genomic formation of Tibeto-Burman speaking populations in Guizhou, Southwest China

Sino-Tibetan is the most prominent language family in East Asia. Previous genetic studies mainly focused on the Tibetan and Han Chinese populations. However, due to the sparse sampling, the genetic structure and admixture history of Tibeto-Burman-speaking populations in the low-altitude region of Southwest China still need to be clarified. We collected DNA from 157 individuals from four Tibeto-Burman-speaking groups from the Guizhou province in Southwest China. We genotyped the samples at about 700,000 genome-wide single nucleotide polymorphisms. Our results indicate that the genetic variation of the four Tibeto-Burman-speaking groups in Guizhou is at the intermediate position in the modern Tibetan-Tai-Kadai/Austronesian genetic cline. This suggests that the formation of Tibetan-Burman groups involved a large-scale gene flow from lowland southern Chinese. The southern ancestry could be further modelled as deriving from Vietnam's Late Neolithic-related inland Southeast Asia agricultural populations and Taiwan's Iron Age-related coastal rice-farming populations. Compared to the Tibeto-Burman speakers in the Tibetan-Yi Corridor reported previously, the Tibeto-Burman groups in the Guizhou region received additional gene flow from the southeast coastal area of China. We show a difference between the genetic profiles of the Tibeto-Burman speakers of the Tibetan-Yi Corridor and the Guizhou province. Vast mountain ranges and rivers in Southwest China may have decelerated the westward expansion of the southeast coastal East Asians. Our results demonstrate the complex genetic profile in the Guizhou region in Southwest China and support the multiple waves of human migration in the southern area of East Asia.

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.

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.

Ancient Mitochondrial Genomes Provide New Clues to the Origin of Domestic Cattle in China

Cattle are one of the six livestock species that have occupied an important place in Chinese history. Previous ancient DNA studies have indicated that Chinese taurine cattle (Bos taurus taurus) are exotic, but the exact route and diffusion by which they were introduced to China is unknown. In this study, we extracted the mitochondrial genomes of 34 cases of ancient taurine cattle (from the late Neolithic to Qin and Han dynasties) excavated from sites in northern China and the eastern Eurasian steppe, and successfully obtained 14 mitochondrial genome sequences. The results of ancient DNA analysis reveal that with cultural exchange and trade, there was close genetic exchange between domestic taurine cattle in different regions. The haplotypes shared by domestic cattle have genetic continuity, reflecting the strong cultural influence of the large capital city sites such as Taosi, Shimao and Erlitou on the surrounding areas. This study suggests that ancient northern Chinese taurine cattle may have accompanied the westward transmission of agricultural or painted pottery culture and thus had a maternal genetic contribution to modern Tibetan cattle.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BACKGROUND

Non-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.

RESULTS

Here, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.

CONCLUSIONS

Our findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

Forensic characteristics and genetic substructure analysis of the Handan Han population, Northern China

We analysed the forensic characteristics and substructure of the Handan Han population based on 36 Y-STR (short tandem repeat) and Y-SNP (single nucleotide polymorphism) markers. The two most dominant haplogroups in Handan Han, O2a2b1a1a1-F8 (17.95%) and O2a2b1a2a1a (21.51%), and their abundant downstream branches, reflected the strong expansion of the precursor of the Hans in Handan. The present results enrich the forensic database and explore the genetic relationships between Handan Han and other neighbouring and/or linguistically close populations, which suggests that the current concise overview of the Han intricate substructure remains oversimplified.

Multiple migrations from East Asia led to linguistic transformation in NorthEast India and mainland Southeast Asia

NorthEast India, with its unique geographic location in the midst of the Himalayas and Bay of Bengal, has served as a passage for the movement of modern humans across the Indian subcontinent and East/Southeast Asia. In this study we look into the population genetics of a unique population called the Khasi, speaking a language (also known as the Khasi language) belonging to the Austroasiatic language family and residing amidst the Tibeto-Burman speakers as an isolated population. The Khasi language belongs to one of the three major broad classifications or phyla of the Austroasiatic language and the speakers of the three sub-groups are separated from each other by large geographical distances. The Khasi speakers are separated from their nearest Austroasiatic language-speaking sub-groups: the “Mundari” sub-family from East and peninsular India and the “Mon-Khmers” in Mainland Southeast Asia. We found the Khasi population to be genetically distinct from other Austroasiatic speakers, i.e. Mundaris and Mon-Khmers, but relatively similar to the geographically proximal Tibeto Burmans. The possible reasons for this genetic-linguistic discordance lie in the admixture history of different migration events that originated from East Asia and proceeded possibly towards Southeast Asia. We found at least two distinct migration events from East Asia. While the ancestors of today’s Tibeto-Burman speakers were affected by both, the ancestors of Khasis were insulated from the second migration event. Correlating the linguistic similarity of Tibeto-Burman and Sino-Tibetan languages of today’s East Asians, we infer that the second wave of migration resulted in a linguistic transition while the Khasis could preserve their linguistic identity.

An initial investigation of dental morphology variation among three southern Naga ethnic groups of Northeast India

OBJECTIVES

This study examines dental morphology trait prevalence among three southern Naga groups and compares them to 10 ethnic groups from other regions of South Asia to accomplish two objectives: assess the biological relationship of these Tibeto-Burman-speakers to speakers of non-Tibeto-Burman languages in other South Asian regions, and determine which traits distinguish northeast Indians from other South Asians.

METHODS

Dental morphology traits were scored with the Arizona State University Dental Anthropology System. Tooth-trait combinations were evaluated for significant inter-trait correlation and intra-trait correspondence within dental fields. Comparisons were based on simple trait prevalence and with Smith's MMD. Affinities based on the former were accomplished with correspondence analysis and principal components analysis. Affinities based on the latter were undertaken with neighbor-joining cluster analysis and multidimensional scaling.

RESULTS

After elimination due to inter-trait correlations and uniform prevalence, biodistances based on the remaining 17 tooth-trait combinations identify significant differences between northeast Indians and other South Asian ethnic groups due to high frequencies of shoveling on the maxillary incisors and Cusp 6 on the mandibular molars coupled with low frequencies of Carabelli's trait and Cusp 5 on UM1 and UM2, respectively.

CONCLUSIONS

Patterns of biodistances obtained from dental morphology are consilient with those obtained from DNA indicating statistically significant differences between northeast Indians from members of ethnic groups of other regions of South Asia. Researchers should explore the sex-specific patterns. Biodistances should not be limited to "key" teeth within dental fields, for in almost every case traits present on mesial and distal teeth yield non-redundant information.

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.

Paternal genetic structure of the Qiang ethnic group in China revealed by high-resolution Y-chromosome STRs and SNPs

The Qiang population mainly lived in Beichuan Qiang Autonomous County of Sichuan Province. It is one of the nomads in China, distributed along the Minjiang River. The Qiang population was assumed to have great affinity with the Han, the largest ethnic group in China, when it refers to the genetic origin. Whereas, it is deeply understudied, especially from the Y chromosome. Here in this study, we used validated high-resolution Y-chromosome single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs) panels to study the Qiang ethnic group to unravel their paternal genetic, forensic and phylogenetic characteristics. A total of 422 male samples of the Qiang ethnic group were genotyped by 233 Y-SNPs and 29 Y-STRs. Haplogroup O-M175 (N = 312) was the most predominant haplogroup in the Qiang ethnic group, followed by D-M174 (N = 32) and C-M130 (N = 32), N-M231 (N = 27), and Q-M242 (N = 15). After further subdivision, O2a-M324 (N = 213) accounted for the majority of haplogroup O. Haplogroup C2b-Z1338 (N = 29), D1a-CTS11577 (N = 30). O2a2b1a1a1-F42 (N = 48), O2a1b1a1a1a-F11 (N = 35), and O2a2b1a1-M117 (N = 21) represented other large terminal haplogroups. The results unveiled that Qiang ethnic group was a population with a high percentage of haplogroup O2a2b1a1a1-F42 (48/422) and O2a1b1a1a1a-F11 (35/422), and O2a2b1a1-M117 (21/422), which has never been reported. Its haplogroup distribution pattern was different from any of the Han populations, implying that the Qiang ethnic group had its unique genetic pattern. Mismatch analysis indicated that the biggest mismatch number in haplogroup O2a2b1a1a1-F42 was 21, while that of haplogroup O2a1b1a1a1a-F11 was 20. The haplotype diversity of the Qiang ethnic group equaled 0.999788, with 392 haplotypes observed, of which 367 haplotypes were unique. The haplogroup diversity of the Qiang ethnic group reached 0.9767, and 53 terminal haplogroups were observed (The haplogroup diversity of the Qiang ethnic group was the highest among Qiang and all Han subgroups, indicating the larger genetic diversity of the Qiang ethnic group.). Haplogroup O2a2b1a1a1-F42 was the most predominant haplogroup, including 11.37 % of the Qiang individuals. Median-joining trees showed gene flow between the Qiang and Han individuals. Our results indicated that 1) the highest genetic diversity was observed in the Qiang ethnic group compared to any of the former studied Chinese population, suggesting that the Qiang might be an older paternal branch; 2) the haplogroup D-M174 individuals of Qiang, Tibetans and Japanese distributed in three different subclades, which was unable to identify through low-resolution Y-SNP panel; and 3) the Qiang had lower proportion of haplogroup D compared to Yi and Tibetan ethnic groups, showing that the Qiang had less genetic communication with them than with Han Chinese.

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.

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.

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.

Origin of ethnic groups, linguistic families, and civilizations in China viewed from the Y chromosome

East Asia, geographically extending to the Pamir Plateau in the west, to the Himalayan Mountains in the southwest, to Lake Baikal in the north and to the South China Sea in the south, harbors a variety of people, cultures, and languages. To reconstruct the natural history of East Asians is a mission of multiple disciplines, including genetics, archaeology, linguistics, and ethnology. Geneticists confirm the recent African origin of modern East Asians. Anatomically modern humans arose in Africa and immigrated into East Asia via a southern route approximately 50,000 years ago. Following the end of the Last Glacial Maximum approximately 12,000 years ago, rice and millet were domesticated in the south and north of East Asia, respectively, which allowed human populations to expand and linguistic families and ethnic groups to develop. These Neolithic populations produced a strong relation between the present genetic structures and linguistic families. The expansion of the Hongshan people from northeastern China relocated most of the ethnic populations on a large scale approximately 5300 years ago. Most of the ethnic groups migrated to remote regions, producing genetic structure differences between the edge and center of East Asia. In central China, pronounced population admixture occurred and accelerated over time, which subsequently formed the Han Chinese population and eventually the Chinese civilization. Population migration between the north and the south throughout history has left a smooth gradient in north-south changes in genetic structure. Observation of the process of shaping the genetic structure of East Asians may help in understanding the global natural history of modern humans.

Genetic reconstruction and phylogenetic analysis by 193 Y-SNPs and 27 Y-STRs in a Chinese Yi ethnic group

Yi is the seventh-largest ethnic group in China and features mountainous regional characteristics. The Liangshan Yi Autonomous Prefecture is the largest Yi agglomeration with isolated geographical conditions, profoundly impeding genetic communication. Here, we investigated 427 unrelated males of Liangshan from 193 Y-chromosome single nucleotide polymorphisms (Y-SNPs) and 27 Y-chromosome short tandem repeats (Y-STRs) to reveal the genetic structure and paternal phylogeny of the group. The haplogroup diversity reached 0.9169 with 46 different subhaplogroups by 193 Y-SNPs, and the haplotype diversity reached 0.9999 by 27 Y-STR loci. Multidimensional scaling (MDS), N-J tree, and Network were constructed to decipher and visualize the genetic relations between Yi and worldwide groups. Our results revealed: (1) the Network by Y-STRs and Y-SNPs showed the haplogroup D1a1a-M15 in the Liangshan Yi population was a ramification of Tibetan groups' expansion from west to east on the plateau; (2) the haplogroup distribution and the mismatch mutation analysis indicated the haplogroup O2a2b1a1a1a4a2-Z25929 of Liangshan Yi derived from manifold Southeast Asian immigrants; (3) a high-resolution Y-SNPs panel is vital to depict accurate paternal derivations and build an integrated and refining genetic structure of ethnic groups.

Genomic insights into the formation of human populations in East Asia

The deep population history of East Asia remains poorly understood owing to a lack of ancient DNA data and sparse sampling of present-day people1,2. Here we report genome-wide data from 166 East Asian individuals dating to between 6000 BC and AD 1000 and 46 present-day groups. Hunter-gatherers from Japan, the Amur River Basin, and people of Neolithic and Iron Age Taiwan and the Tibetan Plateau are linked by a deeply splitting lineage that probably reflects a coastal migration during the Late Pleistocene epoch. We also follow expansions during the subsequent Holocene epoch from four regions. First, hunter-gatherers from Mongolia and the Amur River Basin have ancestry shared by individuals who speak Mongolic and Tungusic languages, but do not carry ancestry characteristic of farmers from the West Liao River region (around 3000 BC), which contradicts theories that the expansion of these farmers spread the Mongolic and Tungusic proto-languages. Second, farmers from the Yellow River Basin (around 3000 BC) probably spread Sino-Tibetan languages, as their ancestry dispersed both to Tibet-where it forms approximately 84% of the gene pool in some groups-and to the Central Plain, where it has contributed around 59-84% to modern Han Chinese groups. Third, people from Taiwan from around 1300 BC to AD 800 derived approximately 75% of their ancestry from a lineage that is widespread in modern individuals who speak Austronesian, Tai-Kadai and Austroasiatic languages, and that we hypothesize derives from farmers of the Yangtze River Valley. Ancient people from Taiwan also derived about 25% of their ancestry from a northern lineage that is related to, but different from, farmers of the Yellow River Basin, which suggests an additional north-to-south expansion. Fourth, ancestry from Yamnaya Steppe pastoralists arrived in western Mongolia after around 3000 BC but was displaced by previously established lineages even while it persisted in western China, as would be expected if this ancestry was associated with the spread of proto-Tocharian Indo-European languages. Two later gene flows affected western Mongolia: migrants after around 2000 BC with Yamnaya and European farmer ancestry, and episodic influences of later groups with ancestry from Turan.

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.

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.

Paternal heritage of the Han Chinese in Henan province (Central China): high diversity and evidence of in situ Neolithic expansions

Background: Due to their long history, complex admixture processes and large population sizes, more research is required to explore the fine genetic structure of Han populations from different geographic locations of China.Aim: To characterise the paternal genetic structure of the Han Chinese in Henan province, which was once the central living region of the ancient Huaxia population, the precursors of the Han Chinese.Subjects and methods: We sequenced Y chromosomes of 60 males from Zhengzhou, Henan Province, and reconstructed a phylogenetic tree for these samples with age estimation.Results: We observed high diversity of paternal lineages in our collection. We found that the in situ Neolithic expansion of the "Major lineages" contributed to a large portion of the paternal gene pool of the Han population in Henan Province. We also detected a large number of "Minor lineages" that diverged in the Palaeolithic Age.Conclusion: We suggest that the high genetic diversity in the paternal gene pool of modern Han populations is mainly attributed to the reservation of a larger number of lineages that diverged in the Palaeolithic Age, while the recent expansion of limited lineages contributed to the majority of the gene pool of modern Han populations. We propose that such a structure is a basal characteristic for the genetic structure of modern Han populations.

Genetic diversity and phylogenetic structure of four Tibeto-Burman-speaking populations in Tibetan-Yi corridor revealed by insertion/deletion polymorphisms

BACKGROUND

Insertion/deletion polymorphisms (InDels), combined with all the desirable features of both short tandem repeat and single nucleotide polymorphism, have been used in archaeological and anthropological research, population genetics and forensic application.

METHODS

Thirty InDels in 530 individuals residing in the Tibetan-Yi corridor (142 Dujiangyan Tibetans, 164 Muli Tibetans, 187 Xichang Yis, and 37 Yanyuan Mosuos) were genotyped using the Investigator DIPplex. Forensic parameters and allele frequency spectrum were calculated. Genetic relationships between the investigated populations and worldwide and nationwide populations were assessed based on both the allele frequency distribution and genotype data.

RESULTS

The combined powers of exclusion were 0.9807 (Dujiangyan Tibetan), 0.9880 (Muli Tibetan), 0.9852 (Xichang Yi) and 0.9892 (Yanyuan Mosuo). The combined powers of discrimination were 0.999999999983 (Dujiangyan Tibetan), 0.999999999942 (Muli Tibetan), 0.999999999982 (Xichang Yi) and 0.999999999962 (Yanyuan Mosuo), respectively. The comprehensive population comparisons among worldwide and nationwide populations uniformly illustrated that the investigated populations have a genetically closer relationship with Tibeto-Burman-speaking populations and geographically adjacent populations.

CONCLUSION

These 30 loci can be regarded as an efficient genetic tool in forensic individual identification and as a supplementary tool in paternity testing in Dujiangyan Tibetan, Muli Tibetan, Xichang Yi, and Yanyuan Mosuo. The genetic proximity between the four populations in the Tibetan-Yi corridor and other populations is strongly correlated with the linguistic origin and geographical distance.

Inferring the population history of Tai-Kadai-speaking people and southernmost Han Chinese on Hainan Island by genome-wide array genotyping

Hainan Island, located between East Asia and Southeast Asia, represents an ideal region for the study of the genetic architecture of geographically isolated populations. However, the genetic structure and demographic history of the indigenous Tai-Kadai-speaking Hlai people and recent expanded southernmost Han Chinese on this island are poorly characterized due to a lack of genetic data. Thus, we collected and genotyped 36 Qiongzhong Hlai and 48 Haikou Han individuals at 497,637 single nucleotide polymorphisms (SNPs). We applied principal component analysis, ADMIXTURE, symmetrical D-statistics, admixture-f₃ statistics, qpWave, and qpAdm analysis to infer the population history. Our results revealed the East Asian populations are characterized by a north-south genetic cline with Hlai at the southernmost end. We have not detected recent gene flow from neighboring populations into Hlai, therefore, we used Hlai as an unadmixed proxy to model the admixture history of mainland Tai-Kadai-speaking populations and southern Han Chinese. The mainland Tai-Kadai-speaking populations are suggested deriving a larger number of their ancestry from Hlai-related lineage, but also having admixture from South Asian-related or other neighboring populations. The Hlai group is also suggested to contribute about half of the ancestry to Han Chinese in Hainan. The complex patterns of genetic structure in East Asia were shaped via language categories, geographical boundaries, and large southward population movements with language dispersal and agriculture propagation.

Paternal Y chromosomal genotyping reveals multiple large-scale admixtures in the formation of Lolo-Burmese-speaking populations in southwest China

Background: Bai and Yi people are two Tibeto-Burman speaking ethnic groups in Yunnan, southwest China. The genetic structure and history of these two groups are largely unknown due to a lack of available genetic data.Aim: To investigate the paternal genetic structure and population relationship of the Yi and Bai people.Subjects and methods: We collected samples from 278 Bai individuals and 283 Yi individuals from Yunnan and subsequently genotyped 43 phylogenetically relevant Y-SNPs in those samples. We estimated haplogroup frequencies and merged our data with a reference database including 46 representative worldwide populations to infer genetic relationships.Results: Y chromosomal haplogroup O-M175 is the dominant lineage in both Bai and Yi people. The Bai and Yi show a close genetic relationship with other Tibeto-Burman-speaking populations with high frequencies of haplogroup O2a2b1a1-Page23, which is also confirmed by PCA. The frequencies of the Tai-Kadai specific lineage O1a-M119, the southern China widespread lineage O1b-P31 and the eastern China enriched lineage O2a1b-002611, are also relatively high in our studied populations.Conclusions: The paternal Y chromosomal affinity of the Bai and Yi with Tibeto-Burman groups is consistent with the language classification. During the formation of the Bai and Yi populations, there were multiple large-scale admixtures, including the expansion of Neolithic farming populations from northern China, the assimilation of Tai-Kadai-speaking populations in southwest China, the demographic expansion driven by Neolithic agricultural revolution from southern China, and the admixture with populations of military immigration from northern and eastern China.

The paternal genetic structure of Jingpo and Dai in southwest China

Yunnan province harbours substantial genetic, cultural and linguistic diversity, with the largest number of Aborigines in China, but the relationship among these Aborigines remains enigmatic. This study genotyped 45 Y chromosomal single nucleotide polymorphisms (SNPs) of 500 males from two aboriginal cross-border populations, Jingpo and Dai, from Dehong, Yunnan. It is reported that Haplogroup O2a2b1a1-M117 is the dominant lineage in both Jingpo and Dai. The Jingpo people show affinity with Tibeto-Burman speaking populations with a relatively high frequency of Haplogroup D-M174, and the Dai people are generally genetically similar with Tai-Kadai speaking populations with high frequencies of Haplogroup O1a-M119 and O1b1a1a-M95, which is consistent with their language classification.

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

The permanent human settlement of the Tibetan Plateau (TP) has been suggested to have been facilitated by the introduction of barley agriculture ∼3.6 kilo-years ago (ka). However, how barley agriculture spread onto the TP remains unknown. Given that the lower altitudes in the northeastern TP were occupied by millet cultivators from 5.2 ka, who also adopted barley farming ∼4 ka, it is highly possible that it was millet farmers who brought barley agriculture onto the TP ∼3.6 ka. To test this hypothesis, we analyzed mitochondrial DNA (mtDNA) from 8277 Tibetans and 58 514 individuals from surrounding populations, including 682 newly sequenced whole mitogenomes. Multiple lines of evidence, together with radiocarbon dating of cereal remains at different elevations, supports the scenario that two haplogroups (M9a1a1c1b1a and A11a1a), which are common in contemporary Tibetans (20.9%) and were probably even more common (40–50%) in early Tibetans prior to historical immigrations to the TP, represent the genetic legacy of the Neolithic millet farmers. Both haplogroups originated in northern China between 10.0–6.0 ka and differentiated in the ancestors of modern Tibetans ∼5.2–4.0 ka, matching the dispersal history of millet farming. By showing that substantial genetic components in contemporary Tibetans can trace their ancestry back to the Neolithic millet farmers, our study reveals that millet farmers adopted and brought barley agriculture to the TP ∼3.6–3.3 ka, and made an important contribution to the Tibetan gene pool.

Genetic structure and forensic characteristics of Tibeto-Burman-speaking Ü-Tsang and Kham Tibetan Highlanders revealed by 27 Y-chromosomal STRs

Culturally diverse Tibetans (Ü-Tsang, Kham and Ando) harboring a unique molecular mechanism that allows them to successfully adapt to hypoxic environments in the Qinghai-Tibet Plateau have been a subject of great interest in medical genetics, linguistics, archeology and forensic science. However, forensic characteristics and genetic variations of the Y-chromosomal 27-marker haplotype included in the Yfiler Plus system in the Ü-Tsang and Kham Tibeto-Burman-speaking Tibetans remain unexplored. Thus, we genotyped 27 Y-STRs in 230 Shigatse Ü-Tsang Tibetans (SUT) and 172 Chamdo Kham Tibetans (CKT) to investigate the forensic characterization and genetic affinity of Chinese Tibetan Highlanders. The haplotype diversities were 0.999962028 in SUT and 0.999796002 in CKT. Forensic diversity measures indicated that this 27-Y-STR amplification system is appropriate for routine forensic applications, such as identifying and separating unrelated males in deficiency paternity cases, male disaster victims and missing person identification and determining male components in sexual assault cases. Moreover, the genetic relationships among 63 worldwide populations (16,282 individuals), 16 Asian populations, and 21 Chinese populations were analyzed and reconstructed using principal component analysis, multidimensional scaling plots and a phylogenetic tree. Considerable genetic differences were observed between Tibetan populations and other geographically/ethnically diverse populations (Han Chinese). Our studied SUT and CKT have a genetically closer relationship with Gansu Ando Tibetans than with other Asians. In total, our analyses indicated that subpopulation structures exist among Asian and Chinese populations, and population-specific reference databases should be established for forensic applications.

Dated language phylogenies shed light on the ancestry of Sino-Tibetan

Given its size and geographical extension, Sino-Tibetan is of the highest importance for understanding the prehistory of East Asia, and of neighboring language families. Based on a dataset of 50 Sino-Tibetan languages, we infer phylogenies that date the origin of the language family to around 7200 B.P., linking the origin of the language family with the late Cishan and the early Yangshao cultures.

The Sino-Tibetan language family is one of the world’s largest and most prominent families, spoken by nearly 1.4 billion people. Despite the importance of the Sino-Tibetan languages, their prehistory remains controversial, with ongoing debate about when and where they originated. To shed light on this debate we develop a database of comparative linguistic data, and apply the linguistic comparative method to identify sound correspondences and establish cognates. We then use phylogenetic methods to infer the relationships among these languages and estimate the age of their origin and homeland. Our findings point to Sino-Tibetan originating with north Chinese millet farmers around 7200 B.P. and suggest a link to the late Cishan and the early Yangshao cultures.

Population Genetic Diversity and Phylogenetic Characteristics for High-Altitude Adaptive Kham Tibetan Revealed by DNATyperTM 19 Amplification System

Tibetans residing in the high-altitude inhospitable environment have undergone significant natural selection of their genetic architecture. Recently, highly mutational autosomal short tandem repeats were widely used not only in the anthropology and population genetics to investigate the genetic structure and relationships, but also in the medical genetics to explore the pathogenesis of multiple genetic diseases and in the forensic science to identify individual and parentage relatedness. However, genetic variants and forensic efficiency of DNATyper^TM 19 amplification system and genetic background of Kham Tibetan remain uncharacterized. Thus, we genotyped 19 forensic genetic markers in 11,402 Kham Tibetans to gain insight into the genetic diversity of Chinese high-altitude adaptive population. Highly discriminating and polymorphic forensic measures were observed, which indicated that this new-developed DNATyper 19 PCR amplification is suitable for routine forensic identification purposes and Chinese national DNA database establishment. Pairwise genetic distances among the comprehensive population comparisons suggested that this high-altitude adaptive Kham Tibetan has genetically closer relationships with lowlanders of Tibeto-Burman-speaking populations (Chengdu Tibetan, Liangshan Tibetan, and Liangshan Yi). Genetic substructure analyses via phylogenetic reconstruction, principal component analysis, and multidimensional scaling analysis in both nationwide and worldwide contexts suggested that the genetic proximity exists along the linguistic, ethnic, and continental geographical boundary. Further studies with whole-genome sequencing of modern or archaic Kham Tibetans would be useful in reconstructing the Tibetan population history.