Genomic structure of the native inhabitants of Peninsular Malaysia and North Borneo suggests complex human population history in Southeast Asia

High-coverage whole-genome sequencing of a Jakun individual from the "Orang Asli" Proto-Malay subtribe from Peninsular Malaysia

Jakun, a Proto-Malay subtribe from Peninsular Malaysia, is believed to have inhabited the Malay Archipelago during the period of agricultural expansion approximately 4 thousand years ago (kya). However, their genetic structure and population history remain inconclusive. In this study, we report the genome structure of a Jakun female, based on whole-genome sequencing, which yielded an average coverage of 35.97-fold. We identified approximately 3.6 million single-nucleotide variations (SNVs) and 517,784 small insertions/deletions (indels). Of these, 39,916 SNVs were novel (referencing dbSNP151), and 10,167 were nonsynonymous (nsSNVs), spanning 5674 genes. Principal Component Analysis (PCA) revealed that the Jakun genome sequence closely clustered with the genomes of the Cambodians (CAM) and the Metropolitan Malays from Singapore (SG_MAS). The ADMIXTURE analysis further revealed potential admixture from the EA and North Borneo populations, as corroborated by the results from the F3, F4, and TreeMix analyses. Mitochondrial DNA analysis revealed that the Jakun genome carried the N21a haplogroup (estimated to have occurred ~19 kya), which is commonly found among Malays from Malaysia and Indonesia. From the whole-genome sequence data, we identified 825 damaging and deleterious nonsynonymous single-nucleotide polymorphisms (nsSNVs) affecting 720 genes. Some of these variants are associated with age-related macular degeneration, atrial fibrillation, and HDL cholesterol level. Additionally, we located a total of 3310 variants on 32 core adsorption, distribution, metabolism, and elimination (ADME) genes. Of these, 193 variants are listed in PharmGKB, and 21 are nsSNVs. In summary, the genetic structure identified in the Jakun individual could enhance the mapping of genetic variants for disease-based population studies and further our understanding of the human migration history in Southeast Asia.

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.

Sequence analyses of Malaysian Indigenous communities reveal historical admixture between Hoabinhian hunter-gatherers and Neolithic farmers

Southeast Asia comprises 11 countries that span mainland Asia across to numerous islands that stretch from the Andaman Sea to the South China Sea and Indian Ocean. This region harbors an impressive diversity of history, culture, religion and biology. Indigenous people of Malaysia display substantial phenotypic, linguistic, and anthropological diversity. Despite this remarkable diversity which has been documented for centuries, the genetic history and structure of indigenous Malaysians remain under-studied. To have a better understanding about the genetic history of these people, especially Malaysian Negritos, we sequenced whole genomes of 15 individuals belonging to five indigenous groups from Peninsular Malaysia and one from North Borneo to high coverage (30X). Our results demonstrate that indigenous populations of Malaysia are genetically close to East Asian populations. We show that present-day Malaysian Negritos can be modeled as an admixture of ancient Hoabinhian hunter-gatherers and Neolithic farmers. We observe gene flow from South Asian populations into the Malaysian indigenous groups, but not into Dusun of North Borneo. Our study proposes that Malaysian indigenous people originated from at least three distinct ancestral populations related to the Hoabinhian hunter-gatherers, Neolithic farmers and Austronesian speakers.

The Peopling and Migration History of the Natives in Peninsular Malaysia and Borneo: A Glimpse on the Studies Over the Past 100 years

Southeast Asia (SEA) has one of the longest records of modern human habitation out-of-Africa. Located at the crossroad of the mainland and islands of SEA, Peninsular Malaysia is an important piece of puzzle to the map of peopling and migration history in Asia, a question that is of interest to many anthropologists, archeologists, and population geneticists. This review aims to revisit our understanding to the population genetics of the natives from Peninsular Malaysia and Borneo over the past century based on the chronology of the technology advancement: 1) Anthropological and Physical Characterization; 2) Blood Group Markers; 3) Protein Markers; 4) Mitochondrial and Autosomal DNA Markers; and 5) Whole Genome Analysis. Subsequently some missing gaps of the study are identified. In the later part of this review, challenges of studying the population genetics of natives will be elaborated. Finally, we conclude our review by reiterating the importance of unveiling migration history and genetic diversity of the indigenous populations as a steppingstone towards comprehending disease evolution and etiology.

Limits and convergence properties of the sequentially Markovian coalescent

Several methods based on the sequentially Markovian coalescent (SMC) make use of full genome sequence data from samples to infer population demographic history including past changes in population size, admixture, migration events and population structure. More recently, the original theoretical framework has been extended to allow the simultaneous estimation of population size changes along with other life history traits such as selfing or seed banking. The latter developments enhance the applicability of SMC methods to nonmodel species. Although convergence proofs have been given using simulated data in a few specific cases, an in-depth investigation of the limitations of SMC methods is lacking. In order to explore such limits, we first develop a tool inferring the best case convergence of SMC methods assuming the true underlying coalescent genealogies are known. This tool can be used to quantify the amount and type of information that can be confidently retrieved from given data sets prior to the analysis of the real data. Second, we assess the inference accuracy when the assumptions of SMC approaches are violated due to departures from the model, namely the presence of transposable elements, variable recombination and mutation rates along the sequence, and SNP calling errors. Third, we deliver a new interpretation of SMC methods by highlighting the importance of the transition matrix, which we argue can be used as a set of summary statistics in other statistical inference methods, uncoupling the SMC from hidden Markov models (HMMs). We finally offer recommendations to better apply SMC methods and build adequate data sets under budget constraints.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

Origin of the pork tapeworm Taenia solium in Bali and Papua, Indonesia

Global distributions of zoonotic pathogens have been strongly affected by the history of human dispersal and domestication of livestock. The pork tapeworm Taenia solium is distributed worldwide as the cause of neurocysticercosis, one of the most serious neglected tropical diseases. T. solium has been reported in Indonesia but only endemic to restricted areas such as Bali and Papua. Previous studies indicated the distinctiveness of a mitochondrial haplotype confirmed in Papua, but only one isolate has been examined to date. In this study, genetic characterization of T. solium and pigs in Bali and Papua was conducted to clarify the distributional history of the parasite. Mitochondrial haplotype network analysis clearly showed that Indonesian T. solium comprises a unique haplogroup which was the first to diverge among Asian genotypes, indicating its single origin and the fact that it was not introduced in the recent past from other area in Asia in which it is endemic. Although phylogenetic analysis based on the mitochondrial D-loop revealed multiple origins of pigs in Bali and Papua, the majority of pigs belonged to the Pacific Clade, which is widely dispersed throughout the Island Southeast Asia (ISEA) and Oceania due to Neolithic human dispersal. Given the results of our network analysis, it is likely that the Pacific Clade pigs played a key role in the dispersal of T. solium. The data suggest that T. solium was introduced from mainland Asia into Western Indonesia, including Bali, by modern humans in the late Pleistocene, or in the early to middle Holocene along with the Pacific Clade pigs. Introduction into New Guinea most likely occurred in the late Holocene through the spread of Pacific Clade pigs. Over time, T. solium has been eradicated from most of Indonesia through the middle to modern ages owing to religious and cultural practices.

Shared Signature of Recent Positive Selection on the TSBP1-BTNL2-HLA-DRA Genes in Five Native Populations from North Borneo

North Borneo (NB) is home to more than 40 native populations. These natives are believed to have undergone local adaptation in response to environmental challenges such as the mosquito-abundant tropical rainforest. We attempted to trace the footprints of natural selection from the genomic data of NB native populations using a panel of ∼2.2 million genome-wide single nucleotide polymorphisms. As a result, an ∼13-kb haplotype in the Major Histocompatibility Complex Class II region encompassing candidate genes TSBP1–BTNL2–HLA-DRA was identified to be undergoing natural selection. This putative signature of positive selection is shared among the five NB populations and is estimated to have arisen ∼5.5 thousand years (∼220 generations) ago, which coincides with the period of Austronesian expansion. Owing to the long history of endemic malaria in NB, the putative signature of positive selection is postulated to be driven by Plasmodium parasite infection. The findings of this study imply that despite high levels of genetic differentiation, the NB populations might have experienced similar local genetic adaptation resulting from stresses of the shared environment.

Autosomal STR Profiling and Databanking in Malaysia: Current Status and Future Prospects

Science and technology are extensively used in criminal investigation. From the mid- to late-1980s, one of the scientific discoveries that has had a particularly remarkable impact on this field has been the use of highly variable DNA sequence regions (minisatellites) in the human genome for individual identification. The technique was initially referred to as DNA fingerprinting, but is now more widely referred to as DNA profiling. Since then, many new developments have occurred within this area of science. These include the introduction of new genetic markers (microsatellites also known as short tandem repeats/STRs), the use of the polymerase chain reaction for target amplification, the development of DNA databases (databanking), and the advancement and/or improvement of genotyping protocols and technologies. In 2019, we described the progress of DNA profiling and DNA databanking in Malaysia for the first time. This report included information on DNA analysis regulations and legislation, STR genotyping protocols, database management, and accreditation status. Here, we provide an update on the performance of our DNA databank (numbers of DNA profiles and hits) plus the technical issues associated with correctly assigning the weight of evidence for DNA profiles in an ethnically diverse population, and the potential application of rapid DNA testing in the country. A total of 116,534 DNA profiles were obtained and stored in the Forensic DNA Databank of Malaysia (FDDM) by 2019, having increased from 70,570 in 2017. The number of hits increased by more than three-fold in just two years, where 17 and 69 hits between the DNA profiles stored in the FDDM and those from crime scenes, suspects, detainees, drug users, convicts, missing persons, or volunteers were recorded in 2017 and 2019, respectively. Forensic DNA analysis and databanking are thus progressing well in Malaysia and have already contributed to many criminal investigations. However, several other issues are discussed here, including the need for STR population data for uncharacterized population groups, and pilot trials for adopting rapid DNA profiling technology. These aspects should be considered by policy makers and law enforcement agencies in order to increase the reliability and efficiency of DNA profiling in criminal cases and in kinship analysis in Malaysia.

Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data

Several methods based on the Sequential Markovian coalescence (SMC) have been developed that make use of genome sequence data to uncover population demographic history, which is of interest in its own right and is a key requirement to generate a null model for selection tests. While these methods can be applied to all possible kind of species, the underlying assumptions are sexual reproduction in each generation and non-overlapping generations. However, in many plants, invertebrates, fungi and other taxa, those assumptions are often violated due to different ecological and life history traits, such as self-fertilization or long term dormant structures (seed or egg-banking). We develop a novel SMC-based method to infer 1) the rates/parameters of dormancy and of self-fertilization, and 2) the populations' past demographic history. Using simulated data sets, we demonstrate the accuracy of our method for a wide range of demographic scenarios and for sequence lengths from one to 30 Mb using four sampled genomes. Finally, we apply our method to a Swedish and a German population of Arabidopsis thaliana demonstrating a selfing rate of ca. 0.87 and the absence of any detectable seed-bank. In contrast, we show that the water flea Daphnia pulex exhibits a long lived egg-bank of three to 18 generations. In conclusion, we here present a novel method to infer accurate demographies and life-history traits for species with selfing and/or seed/egg-banks. Finally, we provide recommendations for the use of SMC-based methods for non-model organisms, highlighting the importance of the per site and the effective ratios of recombination over mutation.

Analysis of five deep-sequenced trio-genomes of the Peninsular Malaysia Orang Asli and North Borneo populations

Background

Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated.

Results

We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10^− 8 – 1.33 × 10^− 8, 1.0 × 10^− 9 – 2.9 × 10^− 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples.

Conclusion

Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.