Reconstructing the Human Genetic History of Mainland Southeast Asia: Insights from Genome-Wide Data from Thailand and Laos

Genome diversity and signatures of natural selection in mainland Southeast Asia

Mainland Southeast Asia (MSEA) has rich ethnic and cultural diversity with a population of nearly 300 million1,2. However, people from MSEA are underrepresented in the current human genomic databases. Here we present the SEA3K genome dataset (phase I), generated by deep short-read whole-genome sequencing of 3,023 individuals from 30 MSEA populations, and long-read whole-genome sequencing of 37 representative individuals. We identified 79.59 million small variants and 96,384 structural variants, among which 22.83 million small variants and 24,622 structural variants are unique to this dataset. We observed a high genetic heterogeneity across MSEA populations, reflected by the varied combinations of genetic components. We identified 44 genomic regions with strong signatures of Darwinian positive selection, covering 89 genes involved in varied physiological systems such as physical traits and immune response. Furthermore, we observed varied patterns of archaic Denisovan introgression in MSEA populations, supporting the proposal of at least two distinct instances of Denisovan admixture into modern humans in Asia3. We also detected genomic regions that suggest adaptive archaic introgressions in MSEA populations. The large number of novel genomic variants in MSEA populations highlight the necessity of studying regional populations that can help answer key questions related to prehistory, genetic adaptation and complex diseases.

Prediction of Skin Color Using Forensic DNA Phenotyping in Asian Populations: A Focus on Thailand

Forensic DNA phenotyping (FDP) has emerged as an essential tool in criminal investigations, enabling the prediction of physical traits based on genetic information. This review explores the genetic factors influencing skin pigmentation, particularly within Asian populations, with a focus on Thailand. Key genes such as Oculocutaneous Albinism II (OCA2), Dopachrome Tautomerase (DCT), KIT Ligand (KITLG), and Solute Carrier Family 24 Member 2 (SLC24A2) are examined for their roles in melanin production and variations that lead to different skin tones. The OCA2 gene is highlighted for its role in transporting ions that help stabilize melanosomes, while specific variants in the DCT gene, including single nucleotide polymorphisms (SNPs) rs2031526 and rs3782974, are discussed for their potential effects on pigmentation in Asian groups. The KITLG gene, crucial for developing melanocytes, includes the SNP rs642742, which is linked to lighter skin in East Asians. Additionally, recent findings on the SLC24A2 gene are presented, emphasizing its connection to pigmentation through calcium regulation in melanin production. Finally, the review addresses the ethical considerations of using FDP in Thailand, where advances in genetic profiling raise concerns about privacy, consent, and discrimination. Establishing clear guidelines is vital to balancing the benefits of forensic DNA applications with the protection of individual rights.

Paternal genetic landscape of contemporary Thai populations in the borderland provinces of Thailand and Myanmar

This study investigates the genetic structure and demographic history of contemporary Thai populations from Tak and Ranong, two border provinces between Thailand and Myanmar with complex immigration histories. We analyzed 20 Y-SNP markers (327 samples) and 24 Y-STR profiles (241 samples), along with published datasets, and observed significant genetic differentiation between these regions. Tak is dominated by O2a2b1a1a (44.53%) and O1b1a1a (21.17%). The presence of NO1 (14.6%) suggests genetic contributions from nearby populations such as the Karen, Lawa, and Khuen. Notably, Tak shows low gene diversity (GD) at the DYS391 (0.2208), contrasting with Ranong (0.5247) and the overall Thailand population (0.4708). Ranong exhibits a more diverse haplogroup distribution, with prominent frequencies of O1b1a1a (20.75%), F (16.98%) and O2a2b1a1a (13.21%). The presence of haplogroups K (11.32%) and R* (7.55%) reflects a genetic connection with the Maniq or other mainland Negrito groups, as well as historical gene flow from South Asia. These findings enhance our understanding of the genetic landscape shaped by migration and localized admixture. We also expanded the Thailand Y-chromosome database with novel Y-STR haplotypes, updated forensic parameters and Y-haplogroup frequencies, providing valuable resources for paternal ancestry inference and kinship verification in both forensic science and population genetics.

Under the name of "Lua": revisiting genetic heterogeneity and population ancestry of Austroasiatic speakers in northern Thailand through genomic analysis

BACKGROUND

Austroasiatic (AA)-speaking populations in northern Thailand are of significant interest due to their status as indigenous descendants and their location at the crossroads of AA prehistoric distribution across Southern China, the Indian Subcontinent, and Mainland Southeast Asia. However, the complexity of ethnic identification can result in inaccuracies regarding the origin and migration history of these populations. To address this, we have conducted a genome-wide SNP analysis of 89 individuals from two Lavue and three Lwa-endonym populations. We then combined our outcomes with previously published data to elucidate the genetic diversity and clustering of AA groups in northern Thailand.

RESULTS

Our findings align with existing linguistic classifications, revealing different genetic compositions among the three branches of the Mon-Khmer subfamily within the AA family: Monic, Khmuic, and Palaungic. Although the term "Lua" ethnicity is confusingly used to identify ethnic groups belonging to both Khmuic and Palaungic branches, our genomic data indicate that the Khmuic-speaking Lua living on the eastern side of the region are relatively distant from the Palaungic-speaking Lavue and Lwa populations living on the western side. The Lavue populations, primarily inhabiting mountainous areas, exhibit a genetic makeup unique to the AA family, with a close genetic relationship to the Karenic subgroup of the Sino-Tibetan language family. Conversely, the Lwa and Blang populations, residing in lowland river valleys, display genetic signatures resulting from admixture with Tai-Kadai-speaking ethnic groups.

CONCLUSION

Utilizing genome-wide SNP markers, our findings indicate genetic heterogeneity among the Lua, Lavue, and Lwa ethnic groups. The intricate interplay of genetics, cultural heritage, and historical influences has shaped these ethnic communities. Our study underscores the importance of accurate ethnic classifications, emphasizing the use of self-identified endonyms, names created and used by the ethnic groups themselves. This approach respects the AA communities in northern Thailand and acknowledges their significant contributions to advancing our understanding of genetic anthropology.

A genotype imputation reference panel specific for native Southeast Asian populations

We report the development of a "Southeast Asian Specific (SEA-specific) Reference Panel" through a "Cross-panel Imputation" approach, consisting of 2550 samples derived from the GA100K, SG10K, and the Peninsular Malaysia Orang Asli (OA) datasets, covering 113,851,450 variants. The SEA-specific panel produced more high confidence variants than 1000 Genomes Project (1KGP) when imputing the OA (8.9 million SEA-specific vs 8.1 million 1KGP) and the Singapore Genome Variation Project (SGVP) (12.5 million SEA-specific vs 11.8 million 1KGP) genotyping datasets. Further, the SEA-specific panel imputed SNPs with better estimated quality scores (INFO, DR2 and R2) on the OA genotyping dataset when comparing with TOPMED and the Human Genome Diversity Project, but performed similarly on SGVP dataset. This panel also exhibited higher recall and non-reference disconcordance rates, indicating the influence of ancestry closeness of the reference panel. However, we note that the imputation accuracy may be compromised by the size of the reference panel.

Genetic variability of 23 autosomal STRs in Austroasiatic-speaking populations from Thailand

Austroasiatic (AA) speakers constitute around 4% of the population of Thailand, while the majority (89.4%) speak Kra-Dai (KD) languages. Previous forensic and population genetic studies in various Thai populations have employed a limited number of short tandem repeats (STRs). This study aims to expand the investigation of the genetic makeup of AA populations in Thailand and their relationship to KD populations using a larger number of autosomal STRs with the VeriFiler™ Plus PCR Amplification Kit. We generated 593 new genotypes from AA-speaking groups and combined them with previously reported data from AA and KD groups. A total of 1,129 genotypes across 23 STR loci were used to construct the largest allelic frequency profile for Thai and Lao populations. However, several loci deviated from Hardy-Weinberg equilibrium, likely due to the reduced genetic diversity in some highland populations, which should be considered in forensic investigations. Beyond forensic applications, our findings reveal genetic differences between AA-speaking groups in Northern and Northeastern Thailand. The AA groups from Northeastern Thailand exhibit greater genetic homogeneity and diversity, likely due to population interactions. In contrast, reduced diversity and increased heterogeneity in AA groups from Northern Thailand are possibly driven by genetic drift and cultural and geographic isolation. In conclusion, we emphasize the usefulness of increasing the number of autosomal STRs in forensic and anthropological genetic studies. Additional Y-STR and X-STR data from various AA-speaking groups in Thailand would further enhance and strengthen forensic STR databases in the region.

Tracing the fine-scale demographic history and recent admixture in Hmong-Mien speakers

The linguistic, historical, and subsistent uniqueness of Hmong-Mien (HM) speakers offers a wonderful opportunity to investigate how these factors impact the genetic structure. The genetic differentiation among HM speakers and their population history are not well characterized. Here, we generate genome-wide data from 65 Yao ethnicity samples and analyze them with published data, particularly by leveraging haplotype-based methods. We determined that the fine-scale genetic substructure of HM speakers corresponds better with linguistic classification than with geography. Particularly, parallels between serial founder events and language differentiations can be observed in West Hmongic speakers. Multiple lines of evidence indicate that ~500-year-old GaoHuaHua individuals are most closely related to West Hmongic-speaking Bunu. The strong genetic bottleneck of some HM-speaking groups, especially Bunu, could potentially be associated with their long-term practice of swidden agriculture to some degree. The inferred admixture dates for most of the HM speakers overlap with the reign of the Ming dynasty (1368-1644 CE). Besides a common genetic origin for HM speakers, their genetic ancestry is shared primarily with neighboring Han Chinese and Tai-Kadai speakers in south China. In conclusion, our analyses reveal that recent isolation and admixture events have contributed to the genetic population history of present-day HM speakers.

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.

Malaria resistance-related biological adaptation and complex evolutionary footprints inferred from one integrative Tai-Kadai-related genomic resource

Pathogen‒host adaptative interactions and complex population demographical processes, including admixture, drift, and Darwen selection, have considerably shaped the Neolithic-to-Modern Western Eurasian population structure and genetic susceptibility to modern human diseases. However, the genetic footprints of evolutionary events in East Asia remain unknown due to the underrepresentation of genomic diversity and the design of large-scale population studies. We reported one aggregated database of genome-wide SNP variations from 796 Tai-Kadai (TK) genomes, including that of Bouyei first reported here, to explore the genetic history, population structure, and biological adaptative features of TK people from southern China and Southeast Asia. We found geography-related population substructure among TK people using the state-of-the-art population genetic structure reconstruction techniques based on the allele frequency spectrum and haplotype-resolved phased fragments. We found that the northern TK people from Guizhou harbored one TK-dominant ancestry maximized in the Bouyei people, and the southern TK people from Thailand were more influenced by Southeast Asians and indigenous people. We reconstructed fitted admixture models and demographic graphs, which showed that TK people received gene flow from ancient southern rice farmer-related lineages related to the Hmong-Mien and Austroasiatic people and from northern millet farmers associated with the Sino-Tibetan people. Biological adaptation focused on our identified unique TK lineages related to Bouyei, which showed many adaptive signatures conferring Malaria resistance and low-rate lipid metabolism. Further gene enrichment, the allele frequency distribution of derived alleles, and their correlation with the incidence of Malaria further confirmed that CR1 played an essential role in the resistance of Malaria in the ancient "Baiyue" tribes.

Forensic efficiency and genetic polymorphisms of 12 X-chromosomal STR loci in Northeastern Thai populations

Northeastern Thailand comprises one-third of the country and is home to various populations, with Lao Isan constituting the majority, while others are considered minority groups. Previous studies on forensic short tandem repeats (STRs) in Thailand predominantly focused on autosomal STRs but there was a paucity of X-STRs, exclusively reported from the North and Central regions of the country. In this study, we have newly established a 12 X-STRs from a total of 896 samples from Northeastern Thailand, encompassing Lao Isan as the major group in the region, alongside nine minor populations (Khmer, Mon, Nyahkur, Bru, Kuy, Phutai, Kalueang, Nyaw, and Saek). Across all ten populations, the combined powers of discrimination in both genders were high and the combined mean exclusion chance (MEC) indices calculated for deficiency, normal trio and duo cases were also high (> 0.99999). DXS10148 emerged as the most informative marker, while DXS7423 was identified as the least informative. Genetic comparison based on X-STRs frequency supported genetic distinction of cerain minor groups such as Kuy, Saek and Nyahkur from other northeastern Thai groups as well as genetic differences according to the geographic region of Thai groups (Northeast, North and Central). In sum, the overall results on population genetics are in agreement with earlier reports on other genetic systems, indicating the informativeness of X-STRs for use in anthropological genetics studies. From a forensic perspective, despite the limitations of small sample sizes for minority groups, the present results contribute to filling the gap in the reference X-STRs database of the major group Lao Isan, providing valuable frequency data for forensic applications in Thailand and neighboring countries.

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.

Genomic portrait and relatedness patterns of the Iron Age Log Coffin culture in northwestern Thailand

The Iron Age of highland Pang Mapha, northwestern Thailand, is characterised by a mortuary practice known as Log Coffin culture. Dating between 2300 and 1000 years ago, large coffins carved from individual teak trees have been discovered in over 40 caves and rock shelters. While previous studies focussed on the cultural development of the Log Coffin-associated sites, the origins of the practice, connections with other wooden coffin-using groups in Southeast Asia, and social structure within the region remain understudied. Here, we present genome-wide data from 33 individuals from five Log Coffin culture sites to study genetic ancestry profiles and genetic interconnectedness. The Log Coffin-associated genomes can be modelled as an admixture between Hòabìnhian hunter-gatherer-, Yangtze River farmer-, and Yellow River farmer-related ancestry. This indicates different influence spheres from Bronze and Iron Age individuals from northeastern Thailand as reflected by cultural practices. Our analyses also identify close genetic relationships within the sites and more distant connections between sites in the same and different river valleys. In combination with high mitochondrial haplogroup diversity and genome-wide homogeneity, the Log Coffin-associated groups from northwestern Thailand seem to have been a large, well-connected community, where genetic relatedness played a significant role in the mortuary ritual.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

Genetic diversity and ancestry of the Khmuic-speaking ethnic groups in Thailand: a genome-wide perspective

The Khmuic-speaking populations are believed to be the descendants of one of the earliest groups to settle in Mainland Southeast Asia. In Thailand, there are two agricultural Khmuic-speaking ethnic groups, the Khamu and Lua (Htin). These peoples primarily reside in scattered locations along the mountainous Thailand-Laos border in Nan province. In this study, we conducted genome-wide SNP analysis on 81 individuals from three Khamu and two Lua villages in northern Thailand. Our findings revealed that both the Khamu and Lua groups possess genetic structures that are distinct from other ethnicities in Southeast Asia, indicating a unique history of migration and settlement. Within the Khmuic group, the Khamu populations living in different locations exhibited similar genetic structures and displayed genetic affinities only with some hill-tribes and Tai-Kadai (Kra-Dai)-speaking groups in Thailand, suggesting potential intermixing or cultural exchange. Furthermore, the Lua people displayed a distinctive population structure, which could be attributed to the founder effect and endogamous marriage practices. Additionally, we discovered a relationship between the Khmuic-speaking populations in Thailand and a Neolithic ancient sample obtained from the Tham Pha Ling archaeological site in Laos. This study provides new insight into genetic substructure within the Khmuic-speaking people and their potential relationship to the indigenous inhabitants of Mainland Southeast Asia.

South Asian maternal and paternal lineages in southern Thailand and the role of sex-biased admixture

Previous genome-wide studies have reported South Asian (SA) ancestry in several Mainland Southeast Asian (MSEA) populations; however, additional details concerning population history, in particular the role of sex-specific aspects of the SA admixture in MSEA populations can be addressed with uniparental markers. Here, we generated ∼2.3 mB sequences of the male-specific portions of the Y chromosome (MSY) of a Tai-Kadai (TK)-speaking Southern Thai group (SouthernThai_TK), and complete mitochondrial (mtDNA) genomes of the SouthernThai_TK and an Austronesian (AN)-speaking Southern Thai (SouthernThai_AN) group. We identified new mtDNA haplogroups, e.g. Q3, E1a1a1, B4a1a and M7c1c3 that have not previously reported in Thai populations, but are frequent in Island Southeast Asia and Oceania, suggesting interactions between MSEA and these regions. SA prevalent mtDNA haplogroups were observed at frequencies of ~35-45% in the Southern Thai groups; both of them showed more genetic relatedness to Austroasiatic (AA) speaking Mon than to any other group. For MSY, SouthernThai_TK had ~35% SA prevalent haplogroups and exhibited closer genetic affinity to Central Thais. We also analyzed published data from other MSEA populations and observed SA ancestry in some additional MSEA populations that also reflects sex-biased admixture; in general, most AA- and AN-speaking groups in MSEA were closer to SA than to TK groups based on mtDNA, but the opposite pattern was observed for the MSY. Overall, our results of new genetic lineages and sex-biased admixture from SA to MSEA groups attest to the additional value that uniparental markers can add to studies of genome-wide variation.

Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations

BACKGROUND

Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear.

RESULTS

Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci.

CONCLUSIONS

Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.

Reanalyzing the genetic history of Kra-Dai speakers from Thailand and new insights into their genetic interactions beyond Mainland Southeast Asia

Thailand is a country where over 60 languages from five language families (Austroasiatic, Austronesian, Hmong-Mien, Kra-Dai, and Sino-Tibetan) are spoken. The Kra-Dai language family is the most prevalent, and Thai, the official language of the country, belongs to it. Previous genome-wide studies on Thailand populations revealed a complex population structure and put some hypotheses forward concerning the population history of the country. However, many published populations have not been co-analyzed, and some aspects of population history were not explored adequately. In this study, we employ new methods to re-analyze published genome-wide genetic data on Thailand populations, with a focus on 14 Kra-Dai-speaking groups. Our analyses reveal South Asian ancestry in Kra-Dai-speaking Lao Isan and Khonmueang, and in Austroasiatic-speaking Palaung, in contrast to a previous study in which the data were generated. We support the admixture scenario for the formation of Kra-Dai-speaking groups from Thailand who harbor both Austroasiatic-related ancestry and Kra-Dai-related ancestry from outside of Thailand. We also provide evidence of bidirectional admixture between Southern Thai and Nayu, an Austronesian-speaking group from Southern Thailand. Challenging some previously reported genetic analyses, we reveal a close genetic relationship between Nayu and Austronesian-speaking groups from Island Southeast Asia (ISEA).

The genomic diversity of Taiwanese Austronesian groups: Implications for the "Into- and Out-of-Taiwan" models

The origin and dispersal of the Austronesian language family, one of the largest and most widespread in the world, have long attracted the attention of linguists, archaeologists, and geneticists. Even though there is a growing consensus that Taiwan is the source of the spread of Austronesian languages, little is known about the migration patterns of the early Austronesians who settled in and left Taiwan, i.e. the "Into-Taiwan" and "out-of-Taiwan" events. In particular, the genetic diversity and structure within Taiwan and how this relates to the into-/out-of-Taiwan events are largely unexplored, primarily because most genomic studies have largely utilized data from just two of the 16 recognized Highland Austronesian groups in Taiwan. In this study, we generated the largest genome-wide data set of Taiwanese Austronesians to date, including six Highland groups and one Lowland group from across the island and two Taiwanese Han groups. We identified fine-scale genomic structure in Taiwan, inferred the ancestry profile of the ancestors of Austronesians, and found that the southern Taiwanese Austronesians show excess genetic affinities with the Austronesians outside of Taiwan. Our findings thus shed new light on the Into- and Out-of-Taiwan dispersals.

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.

Interpreting mismatches between linguistic and genetic patterns among speakers of Tanimuka (Eastern Tukanoan) and Yukuna (Arawakan)

Northwestern Amazonia is home to a great degree of linguistic diversity, and the human societies in that region are part of complex networks of interaction that predate the arrival of Europeans. This study investigates the population and language contact dynamics between two languages found within this region, Yukuna and Tanimuka, which belong to the Arawakan and Tukanoan language families, respectively. We use evidence from linguistics, ethnohistory, ethnography and population genetics to provide new insights into the contact dynamics between these and other human groups in NWA. Our results show that the interaction between these groups intensified in the last 500 years, to the point that it is difficult to differentiate between them genetically. However, this close interaction has led to more substantial contact-induced language changes in Tanimuka than in Yukuna, consistent with a scenario of language shift and asymmetrical power relations.

Genomic footprints of Kam Sweet Rice domestication indicate possible migration routes of the Dong people in China and provide resources for future rice breeding

The Dong people are one of China's 55 recognized ethnic minorities, but there has been a long-standing debate about their origins. In this study, we performed whole-genome resequencing of Kam Sweet Rice (KSR), a valuable, rare, and ancient rice landrace unique to the Dong people. Through comparative genomic analyses of KSR and other rice landraces from south of the Yangtze River Basin in China, we provide evidence that the ancestors of the Dong people likely originated from the southeast coast of China at least 1000 years ago. Alien introgression and admixture in KSR demonstrated multiple migration events in the history of the Dong people. Genomic footprints of domestication demonstrated characteristics of KSR that arose from artificial selection and geographical adaptation by the Dong people. The key genes GS3, Hd1, and DPS1 (related to agronomic traits) and LTG1 and MYBS3 (related to cold tolerance) were identified as domestication targets, reflecting crop improvement and changes in the geographical environment of the Dong people during migration. A genome-wide association study revealed a candidate yield-associated gene, Os01g0923300, a specific haplotype in KSR that is important for regulating grain number per panicle. RNA-sequencing and quantitative reverse transcription-PCR results showed that this gene was more highly expressed in KSR than in ancestral populations, indicating that it may have great value in increasing yield potential in other rice accessions. In summary, our work develops a novel approach for studying human civilization and migration patterns and provides valuable genomic datasets and resources for future breeding of high-yield and climate-resilient rice varieties.

The social lives of isolates (and small language families): the case of the Northwest Amazon

The Americas are home to patches of extraordinary linguistic (genealogical) diversity. These high-diversity areas are particularly unexpected given the recent population of the Americas. In this paper, we zoom in on one such area, the Northwest Amazon, and address the question of how the diversity in this area has persisted to the present. We contrast two hypotheses that claim opposite mechanisms for the maintenance of diversity: the isolation hypothesis suggests that isolation facilitates the preservation of diversity, while the integration hypothesis proposes that conscious identity preservation in combination with contact drives diversity maintenance. We test predictions for both hypotheses across four disciplines: biogeography, cultural anthropology, population genetics and linguistics. Our results show signs of both isolation and integration, but they mainly suggest considerable diversity in how groups of speakers have interacted with their surroundings.

Genome-wide allele and haplotype-sharing patterns suggested one unique Hmong-Mein-related lineage and biological adaptation history in Southwest China

BACKGROUND

Fine-scale genetic structure of ethnolinguistically diverse Chinese populations can fill the gap in the missing diversity and evolutionary landscape of East Asians, particularly for anthropologically informed Chinese minorities. Hmong-Mien (HM) people were one of the most significant indigenous populations in South China and Southeast Asia, which were suggested to be the descendants of the ancient Yangtze rice farmers based on linguistic and archeological evidence. However, their deep population history and biological adaptative features remained to be fully characterized.

OBJECTIVES

To explore the evolutionary and adaptive characteristics of the Miao people, we genotyped genome-wide SNP data in Guizhou HM-speaking populations and merged it with modern and ancient reference populations via a comprehensive population genetic analysis and evolutionary admixture modeling.

RESULTS

The overall genetic admixture landscape of Guizhou Miao showed genetic differentiation between them and other linguistically diverse Guizhou populations. Admixture models further confirmed that Miao people derived their primary ancestry from geographically close Guangxi Gaohuahua people. The estimated identity by descent and effective population size confirmed a plausible population bottleneck, contributing to their unique genetic diversity and population structure patterns. We finally identified several natural selection candidate genes associated with several biological pathways.

CONCLUSIONS

Guizhou Miao possessed a specific genetic structure and harbored a close genetic relationship with geographically close southern Chinese indigenous populations and Guangxi historical people. Miao people derived their major ancestry from geographically close Guangxi Gaohuahua people and experienced a plausible population bottleneck which contributed to the unique pattern of their genetic diversity and structure. Future ancient DNA from Shijiahe and Qujialing will provide new insights into the origin of the Miao people.

Development of molecular diagnostic platform for α0 -thalassemia 44.6 kb (Chiang Rai, --CR ) deletion in individuals with microcytic red blood cells across Thailand

INTRODUCTION

The α⁰ -thalassemia 44.6 kb or Chiang Rai (--^CR ) deletion has been reported in northern Thailand and is capable of causing hemoglobin (Hb) H disease and a lethal α-thalassemia genotype, Hb Bart's hydrops fetalis, in this region. However, there are no current data regarding the frequency of --^CR nationwide due to a lack of effective diagnostic assay. Therefore, this study aimed to develop a reliable platform for simultaneous genotyping of --^CR and two common α⁰ -thalassemias in Thailand (--^SEA and --^THAI ) and investigate the frequency of --^CR across Thailand.

METHODS

Multiplex gap-PCR assay and five renewable plasmid DNA controls for --^CR , --^SEA , --^THAI , α2-globin (HBA2), and β-actin (ACTB) were newly developed and validated with reference methods. The developed assay was further tested on 1046 unrelated individuals with a reduced mean corpuscular volume (MCV) of less than 75 fl for investigating genotypic and allelic spectrum of --^CR .

RESULTS

Our developed assay showed 100% concordance with reference methods. The results were valid and reproducible throughout hundreds of reactions. Comparison of the genotypic and allelic spectra revealed that heterozygous --^SEA (--^SEA /αα) and --^SEA alleles were dominant with the frequency of 22.85% (239/1046) and 13.34% (279/2092), respectively. Of these, --^THAI and --^CR were relatively rare in this population and comparable to each other with the allelic frequency of 0.14% (3/2092).

CONCLUSION

This study successfully established a reliable molecular diagnostic platform for genotyping of --^CR , --^SEA , and --^THAI in a single reaction. Additionally, we demonstrated the frequency of --^CR in Thailand for the first time and provided knowledge basis for the planning of severe α-thalassemia prevention and control programs in Thailand, where thalassemia is endemic.

Genomic perspectives on human dispersals during the Holocene

Nearly 20 y ago, Jared Diamond and Peter Bellwood reviewed the evidence for the associated spread of farming and large language families by the demographic expansions of farmers. Since then, advances in obtaining and analyzing genomic data from modern and ancient populations have transformed our knowledge of human dispersals during the Holocene. Here, we provide an overview of Holocene dispersals in the light of genomic evidence and conclude that they have a complex history. Even when there is a demonstrated connection between a demographic expansion of people, the spread of agriculture, and the spread of a particular language family, the outcome in the results of contact between expanding and resident groups is highly variable. Further research is needed to identify the factors and social circumstances that have influenced this variation and complex history.

Ancient DNA from Protohistoric Period Cambodia indicates that South Asians admixed with local populations as early as 1st-3rd centuries CE

Indian cultural influence is remarkable in present-day Mainland Southeast Asia (MSEA), and it may have stimulated early state formation in the region. Various present-day populations in MSEA harbor a low level of South Asian ancestry, but previous studies failed to detect such ancestry in any ancient individual from MSEA. In this study, we discovered a substantial level of South Asian admixture (ca. 40-50%) in a Protohistoric individual from the Vat Komnou cemetery at the Angkor Borei site in Cambodia. The location and direct radiocarbon dating result on the human bone (95% confidence interval is 78-234 calCE) indicate that this individual lived during the early period of Funan, one of the earliest states in MSEA, which shows that the South Asian gene flow to Cambodia started about a millennium earlier than indicated by previous published results of genetic dating relying on present-day populations. Plausible proxies for the South Asian ancestry source in this individual are present-day populations in Southern India, and the individual shares more genetic drift with present-day Cambodians than with most present-day East and Southeast Asian populations.

Genetic admixture history and forensic characteristics of Tibeto-Burman-speaking Qiang people explored via the newly developed Y-STR panel and genome-wide SNP data

Fine-scale patterns of population genetic structure and diversity of ethnolinguistically diverse populations are important for biogeographical ancestry inference, kinship testing, and development and validation of new kits focused on forensic personal identification. Analyses focused on forensic markers and genome-wide single nucleotide polymorphism (SNP) data can provide new insights into the origin, admixture processes, and forensic characteristics of targeted populations. Qiang people had a large sample size among Tibeto-Burmanspeaking populations, which widely resided in the middle latitude of the Tibetan Plateau. However, their genetic structure and forensic features have remained uncharacterized because of the paucity of comprehensive genetic analyses. Here, we first developed and validated the forensic performance of the AGCU-Y30 Y-short tandem repeats (STR) panel, which contains slowly and moderately mutating Y-STRs, and then we conducted comprehensive population genetic analyses based on Y-STRs and genome-wide SNPs to explore the admixture history of Qiang people and their neighbors. The validated results of this panel showed that the new Y-STR kit was sensitive and robust enough for forensic applications. Haplotype diversity (HD) ranging from 0.9932 to 0.9996 and allelic frequencies ranging from 0.001946 to 0.8326 in 514 Qiang people demonstrated that all included markers were highly polymorphic in Tibeto-Burman people. Population genetic analyses based on Y-STRs [RST, FST, multidimensional scaling (MDS) analysis, neighboring-joining (NJ) tree, principal component analysis (PCA), and median-joining network (MJN)] revealed that the Qiang people harbored a paternally close relationship with lowland Tibetan-Yi corridor populations. Furthermore, we conducted a comprehensive population admixture analysis among modern and ancient Eurasian populations based on genome-wide shared SNPs. We found that the Qiang people were a genetically admixed population and showed closest relationship with Tibetan and Neolithic Yellow River farmers. Admixture modeling showed that Qiang people shared the primary ancestry related to Tibetan, supporting the hypothesis of common origin between Tibetan and Qiang people from North China.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

Genetic substructure of Guizhou Tai-Kadai-speaking people inferred from genome-wide single nucleotide polymorphisms data

The genome-wide characteristics and admixture history of the Tai-Kadai-speaking populations are essential for understanding the population genetic diversity in southern China. We genotyped about 700,000 single nucleotide polymorphisms (SNPs) of 239 individuals from six Tai-Kadai-speaking populations residing in the mountainous Guizhou Province of southwestern China. We merged the genome-wide data with available populations and ancients in East and Southeast Asia to infer Tai-Kadai-speaking populations’ admixture history and genetic structure. We observed a genetic substructure within the studied six populations in the PCA, ADMIXTURE, ChromoPainter, GLOBETROTTER, f-statistics, and qpWave analysis. The Dong, Zhuang, and Bouyei people had a strong genetic affinity with other Tai-Kadai-speaking and Austronesian groups in the surrounding area. However, Gelao showed an affinity to Sino-Tibetan groups, and Mulao people were genetically close to Hmong-Mien populations. qpAdm further illuminated that Gelao and Dong_Tongren composited more Han-related ancestry than Dong, Zhuang, Bouyei, and Mulao people. Meanwhile, we observed high frequencies of Y-chromosome haplogroup O in studied Tai-Kadai-speaking groups except for Gelao people with a high haplogroup N frequency. From the maternal side, haplogroup M7 was frequent in studied populations except for Tongren Dong, who had a high frequency of haplogroup B5. Our newly reported data are helpful for further exploring population dynamics in southern China.

Genomic Insights Into the Demographic History of the Southern Chinese

Southern China is the birthplace of rice-cultivating agriculture and different language families and has also witnessed various human migrations that facilitated cultural diffusions. The fine-scale demographic history in situ that forms present-day local populations, however, remains unclear. To comprehensively cover the genetic diversity in East and Southeast Asia, we generated genome-wide SNP data from 211 present-day Southern Chinese and co-analyzed them with ∼1,200 ancient and modern genomes. In Southern China, language classification is significantly associated with genetic variation but with a different extent of predictability, and there is strong evidence for recent shared genetic history particularly in Hmong–Mien and Austronesian speakers. A geography-related genetic sub-structure that represents the major genetic variation in Southern East Asians is established pre-Holocene and its extremes are represented by Neolithic Fujianese and First Farmers in Mainland Southeast Asia. This sub-structure is largely reduced by admixture in ancient Southern Chinese since > ∼2,000 BP, which forms a “Southern Chinese Cluster” with a high level of genetic homogeneity. Further admixture characterizes the demographic history of the majority of Hmong–Mien speakers and some Kra-Dai speakers in Southwest China happened ∼1,500–1,000 BP, coeval to the reigns of local chiefdoms. In Yellow River Basin, we identify a connection of local populations to genetic sub-structure in Southern China with geographical correspondence appearing > ∼9,000 BP, while the gene flow likely closely related to “Southern Chinese Cluster” since the Longshan period (∼5,000–4,000 BP) forms ancestry profile of Han Chinese Cline.

Ancient genomes from the last three millennia support multiple human dispersals into Wallacea

Previous research indicates that human genetic diversity in Wallacea-islands in present-day Eastern Indonesia and Timor-Leste that were never part of the Sunda or Sahul continental shelves-has been shaped by complex interactions between migrating Austronesian farmers and indigenous hunter-gatherer communities. Yet, inferences based on present-day groups proved insufficient to disentangle this region's demographic movements and admixture timings. Here, we investigate the spatio-temporal patterns of variation in Wallacea based on genome-wide data from 16 ancient individuals (2600-250 years BP) from the North Moluccas, Sulawesi and East Nusa Tenggara. While ancestry in the northern islands primarily reflects contact between Austronesian- and Papuan-related groups, ancestry in the southern islands reveals additional contributions from Mainland Southeast Asia that seem to predate the arrival of Austronesians. Admixture time estimates further support multiple and/or continuous admixture involving Papuan- and Asian-related groups throughout Wallacea. Our results clarify previously debated times of admixture and suggest that the Neolithic dispersals into Island Southeast Asia are associated with the spread of multiple genetic ancestries.

Fully exploiting SNP arrays: a systematic review on the tools to extract underlying genomic structure

Single nucleotide polymorphisms (SNPs) are the most abundant type of genomic variation and the most accessible to genotype in large cohorts. However, they individually explain a small proportion of phenotypic differences between individuals. Ancestry, collective SNP effects, structural variants, somatic mutations or even differences in historic recombination can potentially explain a high percentage of genomic divergence. These genetic differences can be infrequent or laborious to characterize; however, many of them leave distinctive marks on the SNPs across the genome allowing their study in large population samples. Consequently, several methods have been developed over the last decade to detect and analyze different genomic structures using SNP arrays, to complement genome-wide association studies and determine the contribution of these structures to explain the phenotypic differences between individuals. We present an up-to-date collection of available bioinformatics tools that can be used to extract relevant genomic information from SNP array data including population structure and ancestry; polygenic risk scores; identity-by-descent fragments; linkage disequilibrium; heritability and structural variants such as inversions, copy number variants, genetic mosaicisms and recombination histories. From a systematic review of recently published applications of the methods, we describe the main characteristics of R packages, command-line tools and desktop applications, both free and commercial, to help make the most of a large amount of publicly available SNP data.

Fine-Scale Population Admixture Landscape of Tai–Kadai-Speaking Maonan in Southwest China Inferred From Genome-Wide SNP Data

Guizhou Province harbors extensive ethnolinguistic and cultural diversity with Sino-Tibetan-, Hmong–Mien-, and Tai–Kadai-speaking populations. However, previous genetic analyses mainly focused on the genetic admixture history of the former two linguistic groups. The admixture history of Tai–Kadai-speaking populations in Guizhou needed to be characterized further. Thus, we genotyped genome-wide SNP data from 41 Tai–Kadai-speaking Maonan people and made a comprehensive population genetic analysis to explore their genetic origin and admixture history based on the pattern of the sharing alleles and haplotypes. We found a genetic affinity among geographically different Tai–Kadai-speaking populations, especially for Guizhou Maonan people and reference Maonan from Guangxi. Furthermore, formal tests based on the f₃/f₄-statistics further identified an adjacent connection between Maonan and geographically adjacent Hmong–Mien and Sino-Tibetan people, which was consistent with their historically documented shared material culture (Zhang et al., iScience, 2020, 23, 101032). Fitted qpAdm-based two-way admixture models with ancestral sources from northern and southern East Asians demonstrated that Maonan people were an admixed population with primary ancestry related to Guangxi historical people and a minor proportion of ancestry from Northeast Asians, consistent with their linguistically supported southern China origin. Here, we presented the landscape of genetic structure and diversity of Maonan people and a simple demographic model for their evolutionary process. Further whole-genome-sequence–based projects can be presented with more detailed information about the population history and adaptative history of the Guizhou Maonan people.

Indian genetic heritage in Southeast Asian populations

The great ethnolinguistic diversity found today in mainland Southeast Asia (MSEA) reflects multiple migration waves of people in the past. Maritime trading between MSEA and India was established at the latest 300 BCE, and the formation of early states in Southeast Asia during the first millennium CE was strongly influenced by Indian culture, a cultural influence that is still prominent today. Several ancient Indian-influenced states were located in present-day Thailand, and various populations in the country are likely to be descendants of people from those states. To systematically explore Indian genetic heritage in MSEA populations, we generated genome-wide SNP data (using the Affymetrix Human Origins array) for 119 present-day individuals belonging to 10 ethnic groups from Thailand and co-analyzed them with published data using PCA, ADMIXTURE, and methods relying on f-statistics and on autosomal haplotypes. We found low levels of South Asian admixture in various MSEA populations for whom there is evidence of historical connections with the ancient Indian-influenced states but failed to find this genetic component in present-day hunter-gatherer groups and relatively isolated groups from the highlands of Northern Thailand. The results suggest that migration of Indian populations to MSEA may have been responsible for the spread of Indian culture in the region. Our results also support close genetic affinity between Kra-Dai-speaking (also known as Tai-Kadai) and Austronesian-speaking populations, which fits a linguistic hypothesis suggesting cladality of the two language families.

Genomic Insights Into the Population History and Biological Adaptation of Southwestern Chinese Hmong-Mien People

Hmong-Mien (HM) -speaking populations, widely distributed in South China, the north of Thailand, Laos, and Vietnam, have experienced different settlement environments, dietary habits, and pathogenic exposure. However, their specific biological adaptation remained largely uncharacterized, which is important in the population evolutionary genetics and Trans-Omics for regional Precision Medicine. Besides, the origin and genetic diversity of HM people and their phylogenetic relationship with surrounding modern and ancient populations are also unknown. Here, we reported genome-wide SNPs in 52 representative Miao people and combined them with 144 HM people from 13 geographically representative populations to characterize the full genetic admixture and adaptive landscape of HM speakers. We found that obvious genetic substructures existed in geographically different HM populations; one localized in the HM clines, and others possessed affinity with Han Chinese. We also identified one new ancestral lineage specifically existed in HM people, which spatially distributed from Sichuan and Guizhou in the north to Thailand in the south. The sharing patterns of the newly identified homogenous ancestry component combined the estimated admixture times via the decay of linkage disequilibrium and haplotype sharing in GLOBETROTTER suggested that the modern HM-speaking populations originated from Southwest China and migrated southward in the historic period, which is consistent with the reconstructed phenomena of linguistic and archeological documents. Additionally, we identified specific adaptive signatures associated with several important human nervous system biological functions. Our pilot work emphasized the importance of anthropologically informed sampling and deeply genetic structure reconstruction via whole-genome sequencing in the next step in the deep Chinese Population Genomic Diversity Project (CPGDP), especially in the regions with rich ethnolinguistic diversity.

Forensic and genetic characterizations of diverse southern Thai populations based on 15 autosomal STRs

Southern Thailand is home to various populations; the Moklen, Moken and Urak Lawoi’ sea nomads and Maniq negrito are the minority, while the southern Thai groups (Buddhist and Muslim) are the majority. Although previous studies have generated forensic STR dataset for major groups, such data of the southern Thai minority have not been included; here we generated a regional forensic database of southern Thailand. We newly genotyped common 15 autosomal STRs in 184 unrelated southern Thais, including all minorities and majorities. When combined with previously published data of major southern Thais, this provides a total of 334 southern Thai samples. The forensic parameter results show appropriate values for personal identification and paternity testing; the probability of excluding paternity is 0.99999622, and the combined discrimination power is 0.999999999999999. Probably driven by genetic drift and/or isolation with small census size, we found genetic distinction of the Maniq and sea nomads from the major groups, which were closer to the Malay and central Thais than the other Thai groups. The allelic frequency results can strength the regional forensic database in southern Thailand and also provide useful information for anthropological perspective.

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.

Genomic Insight Into the Population Structure and Admixture History of Tai-Kadai-Speaking Sui People in Southwest China

Sui people, which belong to the Tai-Kadai-speaking family, remain poorly characterized due to a lack of genome-wide data. To infer the fine-scale population genetic structure and putative genetic sources of the Sui people, we genotyped 498,655 genome-wide single-nucleotide polymorphisms (SNPs) using SNP arrays in 68 Sui individuals from seven indigenous populations in Guizhou province and Guangxi Zhuang Autonomous Region in Southwest China and co-analyzed with available East Asians via a series of population genetic methods including principal component analysis (PCA), ADMIXTURE, pairwise Fst genetic distance, f-statistics, qpWave, and qpAdm. Our results revealed that Guangxi and Guizhou Sui people showed a strong genetic affinity with populations from southern China and Southeast Asia, especially Tai-Kadai- and Hmong-Mien-speaking populations as well as ancient Iron Age Taiwan Hanben, Gongguan individuals supporting the hypothesis that Sui people came from southern China originally. The indigenous Tai-Kadai-related ancestry (represented by Li), Northern East Asian-related ancestry, and Hmong-Mien-related lineage contributed to the formation processes of the Sui people. We identified the genetic substructure within Sui groups: Guizhou Sui people were relatively homogeneous and possessed similar genetic profiles with neighboring Tai-Kadai-related populations, such as Maonan. While Sui people in Yizhou and Huanjiang of Guangxi might receive unique, additional gene flow from Hmong-Mien-speaking populations and Northern East Asians, respectively, after the divergence within other Sui populations. Sui people could be modeled as the admixture of ancient Yellow River Basin farmer-related ancestry (36.2-54.7%) and ancient coastal Southeast Asian-related ancestry (45.3-63.8%). We also identified the potential positive selection signals related to the disease susceptibility in Sui people via integrated haplotype score (iHS) and number of segregating sites by length (nSL) scores. These genomic findings provided new insights into the demographic history of Tai-Kadai-speaking Sui people and their interaction with neighboring populations in Southern China.

Genomic Insight Into the Population Admixture History of Tungusic-Speaking Manchu People in Northeast China

Manchu is the third-largest ethnic minority in China and has the largest population size among the Tungusic-speaking groups. However, the genetic origin and admixture history of the Manchu people are far from clear due to the sparse sampling and a limited number of markers genotyped. Here, we provided the first batch of genome-wide data of genotyping approximate 700,000 single-nucleotide polymorphisms (SNPs) in 93 Manchu individuals collected from northeast China. We merged the newly generated data with data of publicly available modern and ancient East Asians to comprehensively characterize the genetic diversity and fine-scale population structure, as well as explore the genetic origin and admixture history of northern Chinese Manchus. We applied both descriptive methods of ADMIXTURE, fineSTRUCTURE, F ST , TreeMix, identity by decedent (IBD), principal component analysis (PCA), and qualitative f-statistics (f 3, f 4, qpAdm, and qpWave). We found that Liaoning Manchus have a close genetic relationship and significant admixture signal with northern Han Chinese, which is in line with the cluster patterns in the haplotype-based results. Additionally, the qpAdm-based admixture models showed that modern Manchu people were formed as major ancestry related to Yellow River farmers and minor ancestry linked to ancient populations from Amur River Bain, or others. In summary, the northeastern Chinese Manchu people in Liaoning were an exception to the coherent genetic structure of Tungusic-speaking populations, probably due to the large-scale population migrations and genetic admixtures in the past few hundred years.