Comparison of tet(X4)-containing contigs assembled from metagenomic sequencing data with plasmid sequences of isolates from a cohort of healthy subjects.

Recently discovered tet(X) gene variants have provided new insights into microbial antibiotic resistance mechanisms and their potential consequences for public health. This study focused on detection, analysis, and characterization of Tet(X4)-positive Enterobacterales from the gut microbiota of a healthy cohort of individuals in Singapore using cultivation-dependent and cultivation-independent approaches. Twelve Tet(X4)-positive Enterobacterales strains that were previously obtained from the cohort were fully genome-sequenced and comparatively analyzed. A metagenomic sequencing (MS) data set of the same samples was mined for contigs that harbored the tet(X4) resistance gene. The sequences of tet(X4)-containing contigs and plasmids sequences were compared. The presence of the resistance genes floR and estT (previously annotated as catD) was detected in the same cassette in 10 and 12 out of the 12 tet(X4)-carrying plasmids, respectively. MS detected tet(X4)-containing contigs in 2 out of the 109 subjects, while cultivation-dependent analysis previously reported a prevalence of 10.1%. The tet(X4)-containing sequences assembled from MS data are relatively short (~14 to 33 kb) but show high similarity to the respective plasmid sequences of the isolates. Our findings show that MS can complement efforts in the surveillance of antibiotic resistance genes for clinical samples, while it has a lower sensitivity than a cultivation-based method when the target organism has a low abundance. Further optimization is required if MS is to be utilized in antibiotic resistance surveillance.IMPORTANCEThe global rise in antibiotic resistance makes it necessary to develop and apply new approaches to detect and monitor the prevalence of antibiotic resistance genes in human populations. In this regard, of particular interest are resistances against last-resort antibiotics, such as tigecycline. In this study, we show that metagenomic sequencing can help to detect high abundance of the tigecycline resistance gene tet(X4) in fecal samples from a cohort of healthy human subjects. However, cultivation-based approaches currently remain the most reliable and cost-effective method for detection of antibiotic-resistant bacteria.

Genome-centric analysis of short and long read metagenomes reveals uncharacterized microbiome diversity in Southeast Asians.

Despite extensive efforts to address it, the vastness of uncharacterized 'dark matter' microbial genetic diversity can impact short-read sequencing based metagenomic studies. Population-specific biases in genomic reference databases can further compound this problem. Leveraging advances in hybrid assembly (using short and long reads) and Hi-C technologies in a cross-sectional survey, we deeply characterized 109 gut microbiomes from three ethnicities in Singapore to comprehensively reconstruct 4497 medium and high-quality metagenome assembled genomes, 1708 of which were missing in short-read only analysis and with >28× N50 improvement. Species-level clustering identified 70 (>10% of total) novel gut species out of 685, improved reference genomes for 363 species (53% of total), and discovered 3413 strains unique to these populations. Among the top 10 most abundant gut bacteria in our study, one of the species and >80% of strains were unrepresented in existing databases. Annotation of biosynthetic gene clusters (BGCs) uncovered more than 27,000 BGCs with a large fraction (36-88%) unrepresented in current databases, and with several unique clusters predicted to produce bacteriocins that could significantly alter microbiome community structure. These results reveal significant uncharacterized gut microbial diversity in Southeast Asian populations and highlight the utility of hybrid metagenomic references for bioprospecting and disease-focused studies.

Mutual Exclusion of Methanobrevibacter Species in the Human Gut Microbiota Facilitates Directed Cultivation of a Candidatus Methanobrevibacter Intestini Representative.

Methanogenic Archaea (methanogens) are a phylogenetically diverse group of microorganisms and are considered to be the most abundant archaeal representatives in the human gut. However, the gut methanogen diversity of human populations in many global regions remains poorly investigated. Here, we report the abundance and diversity of gut methanogenic Archaea in a multi-ethnic cohort of healthy Singaporeans by using a concerted approach of metagenomic sequencing, 16S rRNA gene amplicon sequencing, and quantitative PCR. Our results indicate a mutual exclusion of Methanobrevibacter species, i.e., the highly prevalent Methanobrevibacter smithii and the less prevalent Candidatus Methanobrevibacter intestini in more than 80% of the samples when using an amplicon sequencing-based approach. Leveraging on this finding, we were able to select a fecal sample to isolate a representative strain, TLL-48-HuF1, for Candidatus Methanobrevibacter intestini. The analyzed physiological parameters of M. smithii DSM 861T and strain TLL-48-HuF1 suggest high similarity of the two species. Comparative genome analysis and the mutual exclusion of the Methanobrevibacter species indicate potentially different niche adaptation strategies in the human host, which may support the designation of Candidatus M. intestini as a novel species. IMPORTANCE Methanogens are important hydrogen consumers in the gut and are associated with differing host health. Here, we determine the prevalence and abundance of archaeal species in the guts of a multi-ethnic cohort of healthy Singapore residents. While Methanobrevibacter smithii is the most prevalent and abundant methanogen in the human gut of local subjects, the recently proposed Candidatus Methanobrevibacter intestini is the abundant methanogen in a minority of individuals that harbor them. The observed potential mutual exclusion of M. smithii and Ca. M. intestini provides further support to the proposal that the two physiologically similar strains may belong to different Methanobrevibacter species.

The Carbon Footprint of Bioinformatics.

Bioinformatic research relies on large-scale computational infrastructures which have a nonzero carbon footprint but so far, no study has quantified the environmental costs of bioinformatic tools and commonly run analyses. In this work, we estimate the carbon footprint of bioinformatics (in kilograms of CO2 equivalent units, kgCO2e) using the freely available Green Algorithms calculator (www.green-algorithms.org, last accessed 2022). We assessed 1) bioinformatic approaches in genome-wide association studies (GWAS), RNA sequencing, genome assembly, metagenomics, phylogenetics, and molecular simulations, as well as 2) computation strategies, such as parallelization, CPU (central processing unit) versus GPU (graphics processing unit), cloud versus local computing infrastructure, and geography. In particular, we found that biobank-scale GWAS emitted substantial kgCO2e and simple software upgrades could make it greener, for example, upgrading from BOLT-LMM v1 to v2.3 reduced carbon footprint by 73%. Moreover, switching from the average data center to a more efficient one can reduce carbon footprint by approximately 34%. Memory over-allocation can also be a substantial contributor to an algorithm's greenhouse gas emissions. The use of faster processors or greater parallelization reduces running time but can lead to greater carbon footprint. Finally, we provide guidance on how researchers can reduce power consumption and minimize kgCO2e. Overall, this work elucidates the carbon footprint of common analyses in bioinformatics and provides solutions which empower a move toward greener research.

Extended-Spectrum ß-Lactamase-Producing and mcr-1-Positive Escherichia coli from the Gut Microbiota of Healthy Singaporeans.

Multidrug-resistant (MDR) Escherichia coli strains that carry extended-spectrum ß-lactamases (ESBLs) or colistin resistance gene mcr-1 have been identified in the human gut at an increasing incidence worldwide. In this study, we isolated and characterized MDR Enterobacteriaceae from the gut microbiota of healthy Singaporeans and show that the detection rates for ESBL-producing and mcr-positive Enterobacteriaceae are 25.7% (28/109) and 7.3% (8/109), respectively. Whole-genome sequencing analysis of the 37 E. coli isolates assigned them into 25 sequence types and 6 different phylogroups, suggesting that the MDR E. coli gut colonizers are highly diverse. We then analyzed the genetic context of the resistance genes and found that composite transposons played important roles in the cotransfer of blaCTX-M-15/55 and qnrS1, as well as the acquisition of mcr-1. Furthermore, comparative genomic analysis showed that 12 of the 37 MDR E. coli isolates showed high similarity to ESBL-producing E. coli isolates from raw meat products in local markets. By analyzing the core genome single nucleotide polymorphisms (SNPs) shared by these isolates, we identified possible clonal transmission of an MDR E. coli clone between human and raw meat, as well as a group of highly similar IncI2 (Delta) plasmids that might be responsible for the dissemination of mcr-1 in a much wider geographic region. Together, these results suggest that antibiotic resistance may be transmitted between different environmental settings by the expansion of MDR E. coli clones, as well as by the dissemination of resistance plasmids. IMPORTANCE The human gut can harbor both antibiotic-resistant and virulent Escherichia coli which may subsequently cause infections. In this study, we found that multidrug-resistant (MDR) E. coli isolates from the gut of healthy Singaporeans carry a diverse range of antibiotic resistance mechanisms and virulence factor genes and are highly diverse. By comparing their genomes with the extended-spectrum ß-lactamase (ESBL)-producing E. coli isolates from raw meat products that were sampled at a similar time from local markets, we detected an MDR E. coli clone that was possibly transmitted between humans and raw meat products. Furthermore, we also found that a group of resistance plasmids might be responsible for the dissemination of colistin resistance gene mcr-1 in Singapore, Malaysia, and Europe. Our findings call for better countermeasures to block the transmission of antibiotic resistance.

Genome-wide identification of cis DNA methylation quantitative trait loci in three Southeast Asian Populations.

DNA methylation (DNAm) is an epigenetic modification that acts to regulate gene transcription, is essential for cellular processes and plays an important role in complex traits and disease. Variation in DNAm levels is influenced by both genetic and environmental factors. Several studies have examined the extent to which common genetic variation influences DNAm (i.e. mQTLs), however, an improved understanding of mQTLs across diverse human populations is needed to increase their utility in integrative genomic studies in order to further our understanding of complex trait and disease biology. Here, we systematically examine cis-mQTLs in three Southeast Asian populations in the Singapore Integrative Omics (iOmics) Study, comprised of Chinese (n = 93), Indians (n = 83) and Malays (n = 78). A total of 24 851 cis-mQTL probes were associated with at least one SNP in meta- and ethnicity-specific analyses at a stringent significance level. These cis-mQTL probes show significant differences in local SNP heritability between the ethnicities, enrichment in functionally relevant regions using data from the Roadmap Epigenomics Mapping Consortium and are associated with nearby genes and complex traits due to pleiotropy. Importantly, DNAm prediction performance and the replication of cis-mQTLs both within iOmics and between two independent mQTL studies in European and Bangladeshi individuals is best when the genetic distance between the ethnicities is small, with differences in cis-mQTLs likely due to differences in allele frequency and linkage disequilibrium. This study highlights the importance of, and opportunities from, extending investigation of the genetic control of DNAm to Southeast Asian populations.

Paternal gene pool of Malays in Southeast Asia and its applications for the early expansion of Austronesians.

OBJECTIVES: The origin and differentiation of Austronesian populations and their languages have long fascinated linguists, archeologists, and geneticists. However, the founding process of Austronesians and when they separated from their close relatives, such as the Daic and Austro-Asiatic populations in the mainland of Asia, remain unclear. In this study, we explored the paternal origin of Malays in Southeast Asia and the early differentiation of Austronesians. MATERIALS AND METHODS: We generated whole Y-chromosome sequences of 50 Malays and co-analyzed 200 sequences from other Austronesians and related populations. We generated a revised phylogenetic tree with time estimation. RESULTS: We identified six founding paternal lineages among the studied Malays samples. These founding lineages showed a surprisingly coincident expansion age at 5000 to 6000 years ago. We also found numerous mostly close related samples of the founding lineages of Malays among populations from Mainland of Asia. CONCLUSION: Our analyses provided a refined phylogenetic resolution for the dominant paternal lineages of Austronesians found by previous studies. We suggested that the co-expansion of numerous founding paternal lineages corresponds to the initial differentiation of the most recent common ancestor of modern Austronesians. The splitting time and divergence pattern in perspective of paternal Y-chromosome evidence are highly consistent with the previous theories of ethnologists, linguists, and archeologists.

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.

Emergence of tigecycline- and eravacycline-resistant Tet(X4)-producing Enterobacteriaceae in the gut microbiota of healthy Singaporeans.

OBJECTIVES: The recently discovered tigecycline-inactivating enzyme Tet(X4) can confer high-level tigecycline resistance on its hosts, which makes it a public health concern. This study focused on isolation and screening of Tet(X4)-positive Enterobacteriaceae from the gut microbiota of a cohort of healthy individuals in Singapore. METHODS: MinION and Illumina sequencing was performed to obtain the complete genome sequences of Escherichia coli 2EC1-1 and 94EC. Subsequently, 109 human faecal samples were screened retrospectively for eravacycline-resistant Enterobacteriaceae strains, which were further tested for tet(X4) by PCR. The taxonomy of the isolated strains was determined by 16S rRNA gene PCR and Sanger sequencing. RESULTS: Comparative genomic analysis of E. coli 2EC1-1 and 94EC revealed that both carry tet(X4), which is encoded by IncI1-type plasmids p2EC1-1 and p94EC-2, respectively. Retrospective screening of faecal samples collected from 109 healthy individuals showed that the faecal carriage rate of Tet(X4)-producing Enterobacteriaceae is 10.1% (95% CI = 5.1%-17.3%), suggesting that tet(X4) is widely distributed in the gut microbiota of healthy individuals in Singapore. CONCLUSIONS: To the best of our knowledge, this is the first report on the prevalence of tet(X4) in the gut microbiota of a healthy human cohort, as well as the first description of this resistance mechanism outside of China. Our findings suggest that surveillance of tet(X4) in community settings is vital to monitor the spread of this resistance mechanism.

Interethnic analyses of blood pressure loci in populations of East Asian and European descent.

Blood pressure (BP) is a major risk factor for cardiovascular disease and more than 200 genetic loci associated with BP are known. Here, we perform a multi-stage genome-wide association study for BP (max N = 289,038) principally in East Asians and meta-analysis in East Asians and Europeans. We report 19 new genetic loci and ancestry-specific BP variants, conforming to a common ancestry-specific variant association model. At 10 unique loci, distinct non-rare ancestry-specific variants colocalize within the same linkage disequilibrium block despite the significantly discordant effects for the proxy shared variants between the ethnic groups. The genome-wide transethnic correlation of causal-variant effect-sizes is 0.898 and 0.851 for systolic and diastolic BP, respectively. Some of the ancestry-specific association signals are also influenced by a selective sweep. Our results provide new evidence for the role of common ancestry-specific variants and natural selection in ethnic differences in complex traits such as BP.

Dataset on gene expression in the elderly after Mindfulness Awareness Practice or Health Education Program.

It has been reported that relaxation techniques can improve physical health and cognitive function. A number of studies involving different types of relaxation practices showed changes in expression of genes. We investigated the gene expression pattern of a cohort of elderly subjects of Asian descent after weekly (for the first three months) and monthly (for the subsequent six months) intervention. Sixty consenting elderly subjects (aged 60-90 years) with mild cognitive impairment were assigned to either the Mindfulness Awareness Practice (MAP) or Health Education Program (HEP) group in a randomized controlled trial to assess the effectiveness of the programs in preventing further cognitive decline and evaluate the influence on neurological, cellular and biochemical factors. Blood samples were collected before the start of intervention and after nine months for gene expression profiling using Affymetrix Human Genome U133 Plus 2.0 arrays. The dataset is publicly available for further analyses.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error.

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Genetic relatedness of indigenous ethnic groups in northern Borneo to neighboring populations from Southeast Asia, as inferred from genome-wide SNP data.

The region of northern Borneo is home to the current state of Sabah, Malaysia. It is located closest to the southern Philippine islands and may have served as a viaduct for ancient human migration onto or off of Borneo Island. In this study, five indigenous ethnic groups from Sabah were subjected to genome-wide SNP genotyping. These individuals represent the "North Borneo"-speaking group of the great Austronesian family. They have traditionally resided in the inland region of Sabah. The dataset was merged with public datasets, and the genetic relatedness of these groups to neighboring populations from the islands of Southeast Asia, mainland Southeast Asia and southern China was inferred. Genetic structure analysis revealed that these groups formed a genetic cluster that was independent of the clusters of neighboring populations. Additionally, these groups exhibited near-absolute proportions of a genetic component that is also common among Austronesians from Taiwan and the Philippines. They showed no genetic admixture with Austro-Melanesian populations. Furthermore, phylogenetic analysis showed that they are closely related to non-Austro-Melansian Filipinos as well as to Taiwan natives but are distantly related to populations from mainland Southeast Asia. Relatively lower heterozygosity and higher pairwise genetic differentiation index (FST ) values than those of nearby populations indicate that these groups might have experienced genetic drift in the past, resulting in their differentiation from other Austronesians. Subsequent formal testing suggested that these populations have received no gene flow from neighboring populations. Taken together, these results imply that the indigenous ethnic groups of northern Borneo shared a common ancestor with Taiwan natives and non-Austro-Melanesian Filipinos and then isolated themselves on the inland of Sabah. This isolation presumably led to no admixture with other populations, and these individuals therefore underwent strong genetic differentiation. This report contributes to addressing the paucity of genetic data on representatives from this strategic region of ancient human migration event(s).

The fine-scale genetic structure and evolution of the Japanese population.

The contemporary Japanese populations largely consist of three genetically distinct groups-Hondo, Ryukyu and Ainu. By principal-component analysis, while the three groups can be clearly separated, the Hondo people, comprising 99% of the Japanese, form one almost indistinguishable cluster. To understand fine-scale genetic structure, we applied powerful haplotype-based statistical methods to genome-wide single nucleotide polymorphism data from 1600 Japanese individuals, sampled from eight distinct regions in Japan. We then combined the Japanese data with 26 other Asian populations data to analyze the shared ancestry and genetic differentiation. We found that the Japanese could be separated into nine genetic clusters in our dataset, showing a marked concordance with geography; and that major components of ancestry profile of Japanese were from the Korean and Han Chinese clusters. We also detected and dated admixture in the Japanese. While genetic differentiation between Ryukyu and Hondo was suggested to be caused in part by positive selection, genetic differentiation among the Hondo clusters appeared to result principally from genetic drift. Notably, in Asians, we found the possibility that positive selection accentuated genetic differentiation among distant populations but attenuated genetic differentiation among close populations. These findings are significant for studies of human evolution and medical genetics.

Establishing multiple omics baselines for three Southeast Asian populations in the Singapore Integrative Omics Study.

The Singapore Integrative Omics Study provides valuable insights on establishing population reference measurement in 364 Chinese, Malay, and Indian individuals. These measurements include > 2.5 millions genetic variants, 21,649 transcripts expression, 282 lipid species quantification, and 284 clinical, lifestyle, and dietary variables. This concept paper introduces the depth of the data resource, and investigates the extent of ethnic variation at these omics and non-omics biomarkers. It is evident that there are specific biomarkers in each of these platforms to differentiate between the ethnicities, and intra-population analyses suggest that Chinese and Indians are the most biologically homogeneous and heterogeneous, respectively, of the three groups. Consistent patterns of correlations between lipid species also suggest the possibility of lipid tagging to simplify future lipidomics assays. The Singapore Integrative Omics Study is expected to allow the characterization of intra-omic and inter-omic correlations within and across all three ethnic groups through a systems biology approach.The Singapore Genome Variation projects characterized the genetics of Singapore's Chinese, Malay, and Indian populations. The Singapore Integrative Omics Study introduced here goes further in providing multi-omic measurements in individuals from these populations, including genetic, transcriptome, lipidome, and lifestyle data, and will facilitate the study of common diseases in Asian communities.

Phylogeography of Y-chromosome haplogroup O3a2b2-N6 reveals patrilineal traces of Austronesian populations on the eastern coastal regions of Asia.

Austronesian diffusion is considered one of the greatest dispersals in human history; it led to the peopling of an extremely vast region, ranging from Madagascar in the Indian Ocean to Easter Island in Remote Oceania. The Y-chromosome haplogroup O3a2b*-P164(xM134), a predominant paternal lineage of Austronesian populations, is found at high frequencies in Polynesian populations. However, the internal phylogeny of this haplogroup remains poorly investigated. In this study, we analyzed -seventeen Y-chromosome sequences of haplogroup O3a2b*-P164(xM134) and generated a revised phylogenetic tree of this lineage based on 310 non-private Y-chromosome polymorphisms. We discovered that all available O3a2b*-P164(xM134) samples belong to the newly defined haplogroup O3a2b2-N6 and samples from Austronesian populations belong to the sublineage O3a2b2a2-F706. Additionally, we genotyped a series of Y-chromosome polymorphisms in a large collection of samples from China. We confirmed that the sublineage O3a2b2a2b-B451 is unique to Austronesian populations. We found that O3a2b2-N6 samples are widely distributed on the eastern coastal regions of Asia, from Korea to Vietnam. Furthermore, we propose- that the O3a2b2a2b-B451 lineage represents a genetic connection between ancestors of Austronesian populations and ancient populations in North China, where foxtail millet was domesticated about 11,000 years ago. The large number of newly defined Y-chromosome polymorphisms and the revised phylogenetic tree of O3a2b2-N6 will be helpful to explore the origin of proto-Austronesians and the early diffusion process of Austronesian populations.

Characterising private and shared signatures of positive selection in 37 Asian populations.

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

Mapping the genetic diversity of HLA haplotypes in the Japanese populations.

Japan has often been viewed as an Asian country that possesses a genetically homogenous community. The basis for partitioning the country into prefectures has largely been geographical, although cultural and linguistic differences still exist between some of the districts/prefectures, especially between Okinawa and the mainland prefectures. The Major Histocompatibility Complex (MHC) region has consistently emerged as the most polymorphic region in the human genome, harbouring numerous biologically important variants; nevertheless the presence of population-specific long haplotypes hinders the imputation of SNPs and classical HLA alleles. Here, we examined the extent of genetic variation at the MHC between eight Japanese populations sampled from Okinawa, and six other prefectures located in or close to the mainland of Japan, specifically focusing at the haplotypes observed within each population, and what the impact of any variation has on imputation. Our results indicated that Okinawa was genetically farther to the mainland Japanese than were Gujarati Indians from Tamil Indians, while the mainland Japanese from six prefectures were more homogeneous than between northern and southern Han Chinese. The distribution of haplotypes across Japan was similar, although imputation was most accurate for Okinawa and several mainland prefectures when population-specific panels were used as reference.

Dissecting the genetic structure and admixture of four geographical Malay populations.

The Malay people are an important ethnic composition in Southeast Asia, but their genetic make-up and population structure remain poorly studied. Here we conducted a genome-wide study of four geographical Malay populations: Peninsular Malaysian Malay (PMM), Singaporean Malay (SGM), Indonesian Malay (IDM) and Sri Lankan Malay (SLM). All the four Malay populations showed substantial admixture with multiple ancestries. We identified four major ancestral components in Malay populations: Austronesian (17%-62%), Proto-Malay (15%-31%), East Asian (4%-16%) and South Asian (3%-34%). Approximately 34% of the genetic makeup of SLM is of South Asian ancestry, resulting in its distinct genetic pattern compared with the other three Malay populations. Besides, substantial differentiation was observed between the Malay populations from the north and the south, and between those from the west and the east. In summary, this study revealed that the genetic identity of the Malays comprises a mixed entity of multiple ancestries represented by Austronesian, Proto-Malay, East Asian and South Asian, with most of the admixture events estimated to have occurred 175 to 1,500 years ago, which in turn suggests that geographical isolation and independent admixture have significantly shaped the genetic architectures and the diversity of the Malay populations.