Phylogeographic and genome-wide investigations of Vietnam ethnic groups reveal signatures of complex historical demographic movements

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.

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.

The first maternal genetic study of hunter-gatherers from Vietnam

The current limitation of ancient DNA data from Vietnam led to the controversy surrounding the prehistory of people in this region. The combination of high heat and humidity damaged ancient bones that challenged the study of human evolution, especially when using DNA as study materials. So far, only 4 k years of history have been recorded despite the 65 k years of history of anatomically modern human occupations in Vietnam. Here we report, to our knowledge, the oldest mitogenomes of two hunter-gatherers from Vietnam. We extracted DNA from the femurs of two individuals aged 6.2 k cal BP from the Con Co Ngua (CCN) site in Thanh Hoa, Vietnam. This archeological site is the largest cemetery of the hunter-gatherer population in Southeast Asia (SEA) that was discovered, but their genetics have not been explored until the present. We indicated that the CCN haplotype belongs to a rare haplogroup that was not detected in any present-day Vietnamese individuals. Further matrilineal analysis on CCN mitogenomes showed a close relationship with ancient farmers and present-day populations in SEA. The mitogenomes of hunter-gatherers from Vietnam debate the "two layers" model of peopling history in SEA and provide an alternative solution for studying challenging ancient human samples from Vietnam.

The Analysis of Genetic Polymorphism on Mitochondrial Hypervariable Region III in Thai Population

Mitochondrial DNA (mtDNA) analysis is a genetic marker for human identification, especially matrilineal inheritance. Hypervariable regions (HVR) I and II of mtDNA have been currently performed for human identification worldwide. Further examination of HVRIII has been conducted with the aim of enhancing the power of discrimination. The aim of this research is to provide informative data on the polymorphisms of HVRIII in the Thai population in order to establish a national database for human identification. Thai people who were unrelated through the maternal lineage were recruited for blood collections. The mtDNA was extracted by Chelex extraction, amplified by polymerase chain reaction, and analyzed using Sequencing Analysis Software. The most common mutation in HVRIII was base substitution, followed by deletion and insertion. We discovered 40 unique haplotypes, with haplotype 489C being the most frequent. The haplotype diversity, power of discrimination, and random match probability were 0.8014, 0.7987, and 0.2013, respectively. Five-CA repeats were the most frequently observed in nucleotide positions 514-523. Our database can be employed as supplementary markers in addition to nuclear deoxyribonucleic acid (DNA) markers in forensic investigations. Moreover, the data could potentially enhance genetic identification and anthropological genetics research in Thailand.

The endemic Helicobacter pylori population in Southern Vietnam has both South East Asian and European origins

Background

The burden of Helicobacter pylori-induced gastric cancer varies based on predominant H. pylori population in various geographical regions. Vietnam is a high H. pylori burden country with the highest age-standardized incidence rate of gastric cancer (16.3 cases/100,000 for both sexes) in Southeast Asia, despite this data on the H. pylori population is scanty. We examined the global context of the endemic H. pylori population in Vietnam and present a contextual and comparative genomics analysis of 83 H. pylori isolates from patients in Vietnam.

Results

There are at least two major H. pylori populations are circulating in symptomatic Vietnamese patients. The majority of the isolates (~ 80%, 66/83) belong to the hspEastAsia and the remaining belong to hpEurope population (~ 20%, 17/83). In total, 66 isolates (66/83) were cagA positive, 64 were hspEastAsia isolates and two were hpEurope isolates. Examination of the second repeat region revealed that most of the cagA genes were ABD type (63/66; 61 were hspEastAsia isolates and two were hpEurope isolates). The remaining three isolates (all from hspEastAsia isolates) were ABC or ABCC types. We also detected that 4.5% (3/66) cagA gene from hspEastAsia isolates contained EPIYA-like sequences, ESIYA at EPIYA-B segments. Analysis of the vacA allelic type revealed 98.8% (82/83) and 41% (34/83) of the strains harboured the s1 and m1 allelic variant, respectively; 34/83 carried both s1m1 alleles. The most frequent genotypes among the cagA positive isolates were vacA s1m1/cagA + and vacA s1m2/cagA + , accounting for 51.5% (34/66) and 48.5% (32/66) of the isolates, respectively.

Conclusions

There are two predominant lineages of H. pylori circulating in Vietnam; most of the isolates belong to the hspEastAsia population. The hpEurope population is further divided into two smaller clusters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13099-021-00452-2.

Genotypes of informative loci from 1000 Genomes data allude evolution and mixing of human populations

Principal Component Analysis (PCA) projects high-dimensional genotype data into a few components that discern populations. Ancestry Informative Markers (AIMs) are a small subset of SNPs capable of distinguishing populations. We integrate these two approaches by proposing an algorithm to identify necessary informative loci whose removal from the data deteriorates the PCA structure. Unlike classical AIMs, necessary informative loci densely cover the genome, hence can illuminate the evolution and mixing history of populations. We conduct a comprehensive analysis to the genotype data of the 1000 Genomes Project using necessary informative loci. Projections along the top seven principal components demarcate populations at distinct geographic levels. Millions of necessary informative loci along each PC are identified. Population identities along each PC are approximately determined by weighted sums of minor (or major) alleles over the informative loci. Variations of allele frequencies are aligned with the history and direction of population evolution. The population distribution of projections along the top three PCs is recapitulated by a simple demographic model based on several waves of founder population separation and mixing. Informative loci possess locational concentration in the genome and functional enrichment. Genes at two hot spots encompassing dense PC 7 informative loci exhibit differential expressions among European populations. The mosaic of local ancestry in the genome of a mixed descendant from multiple populations can be inferred from partial PCA projections of informative loci. Finally, informative loci derived from the 1000 Genomes data well predict the projections of an independent genotype data of South Asians. These results demonstrate the utility and relevance of informative loci to investigate human evolution.

Extensive Ethnolinguistic Diversity in Vietnam Reflects Multiple Sources of Genetic Diversity

Vietnam features extensive ethnolinguistic diversity and occupies a key position in Mainland Southeast Asia. Yet, the genetic diversity of Vietnam remains relatively unexplored, especially with genome-wide data, because previous studies have focused mainly on the majority Kinh group. Here, we analyze newly generated genome-wide single-nucleotide polymorphism data for the Kinh and 21 additional ethnic groups in Vietnam, encompassing all five major language families in Mainland Southeast Asia. In addition to analyzing the allele and haplotype sharing within the Vietnamese groups, we incorporate published data from both nearby modern populations and ancient samples for comparison. In contrast to previous studies that suggested a largely indigenous origin for Vietnamese genetic diversity, we find that Vietnamese ethnolinguistic groups harbor multiple sources of genetic diversity that likely reflect different sources for the ancestry associated with each language family. However, linguistic diversity does not completely match genetic diversity: There have been extensive interactions between the Hmong-Mien and Tai-Kadai groups; different Austro-Asiatic groups show different affinities with other ethnolinguistic groups; and we identified a likely case of cultural diffusion in which some Austro-Asiatic groups shifted to Austronesian languages during the past 2,500 years. Overall, our results highlight the importance of genome-wide data from dense sampling of ethnolinguistic groups in providing new insights into the genetic diversity and history of an ethnolinguistically diverse region, such as Vietnam.

Population genetic data of 22 autosomal STR loci for the Mong people in Vietnam

This study investigated 22 autosomal short tandem repeat (STR) loci in 156 unrelated individuals from the Mong ethnic minority in Ha Giang Province, Vietnam. Allele frequencies and forensic parameters were calculated, showing the combined Powers of Discrimination reaching 1.000000000000000000000000000000 and the combined Power of Exclusion greater than 0.999999986623. Phylogenetic analysis indicated that the Vietnamese Mong population has close genetic relationships with other Hmong-Mien populations.

SpatialCPie: an R/Bioconductor package for spatial transcriptomics cluster evaluation

Background

Technological developments in the emerging field of spatial transcriptomics have opened up an unexplored landscape where transcript information is put in a spatial context. Clustering commonly constitutes a central component in analyzing this type of data. However, deciding on the number of clusters to use and interpreting their relationships can be difficult.

Results

We introduce SpatialCPie, an R package designed to facilitate cluster evaluation for spatial transcriptomics data. SpatialCPie clusters the data at multiple resolutions. The results are visualized with pie charts that indicate the similarity between spatial regions and clusters and a cluster graph that shows the relationships between clusters at different resolutions. We demonstrate SpatialCPie on several publicly available datasets.

Conclusions

SpatialCPie provides intuitive visualizations of cluster relationships when dealing with Spatial Transcriptomics data.

Vitamin D receptor ApaI polymorphism associated with progression of liver disease in Vietnamese patients chronically infected with hepatitis B virus

Background

Vitamin D derivatives and their receptor (VDR) are potent modulators of immune responses in various diseases including malignancies as well as in metabolic and infectious disorders. The impact of vitamin D receptor polymorphisms on clinical outcomes of hepatitis B virus (HBV) infection is not well understood. This study aims to investigate the potential role of VDR polymorphisms (TaqI, FokI, ApaI, and BsmI) in Vietnamese HBV infected patients and to correlate these polymorphisms with the progression of HBV-related liver disease.

Methods

Four hundred forty-three HBV infected patients of the three clinically well-defined subgroups chronic hepatitis B (CHB, n = 183), liver cirrhosis (LC, n = 89) and hepatocellular carcinoma (HCC, n = 171) and 238 healthy individuals (HC) were enrolled. VDR polymorphisms were genotyped by DNA sequencing and in-house validated ARMS assays. Logistic regression models were applied in order to determine the association of VDR polymorphisms with manifest HBV infection as well as with progression of related liver diseases mulin different genetic models.

Results

The VDR ApaI CA genotype was less frequent in HCC than in CHB patients in different genetic models (codominant model, OR = 0.5, 95%CI = 0.3–0.84, P = 0.004; dominant model, OR = 0.46, 95%CI = 0.27–0.76, P = 0.0023). In the recessive model, the genotype ApaI AA was found more frequently among HCC compared to CHB patients (OR = 2.56, 95%CI = 1.01–6.48, P = 0.04). Similarly, the ApaI CA genotype was less frequent in HCC than in non-HCC group codominant model, OR = 0.6, 95%CI = 0.4–0.98, dominant model, P = 0.04 and OR = 0.6, 95%CI = 0.38–0.90, P = 0.017). The ApaI genotypes CA and AA was significantly associated with higher levels of liver enzymes, bilirubin, and HBV DNA (P < 0.05). No association between TaqI, FokI and BsmI polymorphisms and any clinical outcome as well as liver disease progression was found.

Conclusions

Among the four investigated VDR polymorphisms, ApaI is associated with clinical outcome and liver disease progression in Vietnamese HBV infected patients.

Genotyping of Mycobacterium tuberculosis spreading in Hanoi, Vietnam using conventional and whole genome sequencing methods

Hanoi is the capital of Vietnam, one of the 30 countries with a high tuberculosis (TB) burden. Fundamental data on the molecular epidemiology of the disease is required for future TB management. To identify lineages and genotypes of Mycobacterium tuberculosis (Mtb), conventional genotyping data from clinical isolates of the Hanoi area was compared with whole genome sequencing (WGS) analysis from 332 of 470 samples. It was obtained from lineage-specific single nucleotide variants (SNVs), large sequence polymorphisms, spoligotyping, and variable number of tandem repeats (VNTR) analysis using mycobacterial interspersed repetitive unit (MIRU) and Japan anti-tuberculosis association (JATA) locus sets. This information was directly compared with results obtained from WGS. Mini-satellite repeat unit variants were identified using BLAST search against concatenated short read sequences, the RepUnitTyping tool. WGS analysis revealed that the Mtb strains tested are diverse and classified into lineage (L) 1, 2 and 4 (24.7, 57.2 and 18.1% respectively). The majority of the L2 strains were further divided into ancient and modern Beijing genotypes, and most of the L1 group were EAI4_VNM strains. Although conventional PCR-based genotyping results were mostly consistent with information obtained through WGS analysis, in-depth analysis identified aberrant deletions and spacers that may cause discordance. JATA-VNTR sets, including hypervariable loci, separated large Beijing genotypic clusters generated by MIRU15 into smaller groups. The distribution of repeat unit variants observed within 33 VNTR loci showed clear variation depending on the three lineages. WGS-based pairwise-SNV differences within VNTR-defined genotypic clusters were greater in L1 than in L2 and L4 (P = .001). Direct comparisons between results of PCR-based genotyping and in silico analysis of WGS data would bridge a gap between classical and modern technologies during this transition period, and provide further information on Mtb genotypes in specific geographical areas.

Craniofacial Anthropometric Profile of East Asians

Facial anthropometric data has significant ethnic variation. East Asia, comprised of fourteen countries, represents a significant proportion of the global population. This systematic review presents the facial anthropometric data collected from these countries. The systematic review was conducted in accordance with PRISMA guidelines. An electronic search of the MEDLINE database returned 3054 articles. Twenty articles were considered eligible for inclusion. Nine studies were conducted in China, 1 in Indonesia, 2 in Japan,3 in Korea, 4 in Malaysia, and 1 was a multicentre study conducted in China, Japan, Thailand, and Vietnam. Qualitative and quantitative parameters were extracted from the20 studies. No data was found for the other East Asian countries. There is a paucity of facial anthropometric data for East Asian countries despite their high burden of craniofacial anomalies and a strong demand for cosmetic facial surgery, both of which would benefit from the collection of robust craniofacial norms. It is in the interest of both the craniofacial surgeon and the East Asian patient to collect baseline facial anthropometric data for this population.

Developing pharmacogenetic screening methods for an emergent country: Vietnam

BACKGROUND

The finding of strong associations between certain human leukocyte antigen (HLA) genotypes and the development of severe cutaneous adverse drug reactions (SCARs), [for example, HLA-B*57:01 and abacavir (ABC), HLA-B*15:02 and carbamazepine (CBZ) and HLA-B*58:01 and allopurinol], has led to HLA screening being used to prevent SCARs. Screening has been shown to be of great benefit in a number of studies. Clinical translation from bench to bedside, however, depends upon the development of simple, rapid and cost-effective assays to detect these risk alleles. In highly populated developing countries such as Vietnam, where there is a high prevalence of HLA-B*15:02 and HLA-B*58:01 correlating with a high incidence of CBZ- and allopurinol-induced SCARs, the crucial factor in the implementation of comprehensive screening programs to detect these major risk HLA alleles is the availability of suitable assays.

BODY

We have summarized the role and economic benefits of HLA screening, reviewed published HLA screening methods used currently in pharmacogenetic screening and examined the advantages and disadvantages of assays developed specifically for use in screening for risk alleles in the prevention of HLA-associated SCARs in Vietnam.

CONCLUSION

The optimal approach we propose may serve as a template for the development of screening programs in other emergent countries.

Elucidation of the unique mutation spectrum of severe hearing loss in a Vietnamese pediatric population

The mutational spectrum of deafness in Indochina Peninsula, including Vietnam, remains mostly undetermined. This significantly hampers the progress toward establishing an effective genetic screening method and early customized rehabilitation modalities for hearing loss. In this study, we evaluated the genetic profile of severe-to-profound hearing loss in a Vietnamese pediatric population using a hierarchical genetic analysis protocol that screened 11 known deafness-causing variants, followed by massively parallel sequencing targeting 129 deafness-associated genes. Eighty-seven children with isolated severe-to-profound non-syndromic hearing loss without family history were included. The overall molecular diagnostic yield was estimated to be 31.7%. The mutational spectrum for severe-to-profound non-syndromic hearing loss in our Vietnamese population was unique: The most prevalent variants resided in the MYO15A gene (7.2%), followed by GJB2 (6.9%), MYO7A (5.5%), SLC26A4 (4.6%), TMC1 (1.8%), ESPN (1.8%), POU3F4 (1.8%), MYH14 (1.8%), EYA1 (1.8%), and MR-RNR1 (1.1%). The unique spectrum of causative genes in the Vietnamese deaf population was similar to that in the southern Chinese deaf population. It is our hope that the mutation spectrum provided here could aid in establishing an efficient protocol for genetic analysis of severe-to-profound hearing loss and a customized screening kit for the Vietnamese population.

Complete human mtDNA genome sequences from Vietnam and the phylogeography of Mainland Southeast Asia

Vietnam is an important crossroads within Mainland Southeast Asia (MSEA) and a gateway to Island Southeast Asia, and as such exhibits high levels of ethnolinguistic diversity. However, comparatively few studies have been undertaken of the genetic diversity of Vietnamese populations. In order to gain comprehensive insights into MSEA mtDNA phylogeography, we sequenced 609 complete mtDNA genomes from individuals belonging to five language families (Austroasiatic, Tai-Kadai, Hmong-Mien, Sino-Tibetan and Austronesian) and analyzed them in comparison with sequences from other MSEA countries and Taiwan. Within Vietnam, we identified 399 haplotypes belonging to 135 haplogroups; among the five language families, the sequences from Austronesian groups differ the most from the other groups. Phylogenetic analysis revealed 111 novel Vietnamese mtDNA lineages. Bayesian estimates of coalescence times and associated 95% HPD for these show a peak of mtDNA diversification around 2.5–3 kya, which coincides with the Dong Son culture, and thus may be associated with the agriculturally-driven expansion of this culture. Networks of major MSEA haplogroups emphasize the overall distinctiveness of sequences from Taiwan, in keeping with previous studies that suggested at most a minor impact of the Austronesian expansion from Taiwan on MSEA. We also see evidence for population expansions across MSEA geographic regions and language families.