Caries trajectories from childhood to adolescence: Analysis of data from a nationwide school dental service.

OBJECTIVES: The aim of the study was to assess patterns of longitudinal changes in caries status among school-going children in Singapore. METHODS: Dental records for a single cohort of students who received dental examinations in six standard examination years between 2009 and 2017 were analysed (n = 24 699). Group-based trajectory modelling with a zero-inflated Poisson distribution was carried out to determine dental caries trajectories in the permanent dentition. Associations between sociodemographic factors and trajectory group membership were assessed using multinomial logistic regression. RESULTS: The predicted population distribution across the four caries trajectory groups identified was 65.0% ('none'), 16.8% ('low'), 14.8% ('medium') and 3.4% ('high'). The 'none' trajectory group had a decayed, missing and filled teeth (DMFT) score of 0 throughout the 8 years. Higher baseline DMFT counts and nonlinear increases in DMFT scores were noted for the 'low', 'medium' and 'high' trajectory groups. The correlation coefficient between DMFT counts in years 6 and 8 was 0.91, as compared to 0.77 between baseline and year 1. Factors associated with the 'high' caries trajectory include lower socio-economic status, female gender, Chinese race (compared to the Indian race), enrolment in primary schools in the Eastern and Western regions of Singapore, and enrolment in public secondary schools. CONCLUSIONS: Under a nationwide school dental service, four trajectory patterns of caries counts in the permanent dentition were identified over 8 years. Among students in the 'low', 'medium' and 'high' trajectory groups, greater caries increment was noted during the transition from primary to secondary school. The correlation between DMFT counts in successive examinations was stronger in older than younger ages.

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.

Whole genome sequencing reveals hidden transmission of carbapenemase-producing Enterobacterales.

Carbapenemase-producing Enterobacterales (CPE) infection control practices are based on the paradigm that detected carriers in the hospital transmit to other patients who stay in the same ward. The role of plasmid-mediated transmission at population level remains largely unknown. In this retrospective cohort study over 4.7 years involving all multi-disciplinary public hospitals in Singapore, we analysed 779 patients who acquired CPE (1215 CPE isolates) detected by clinical or surveillance cultures. 42.0% met putative clonal transmission criteria, 44.8% met putative plasmid-mediated transmission criteria and 13.2% were unlinked. Only putative clonal transmissions associated with direct ward contact decreased in the second half of the study. Both putative clonal and plasmid-mediated transmission associated with indirect (no temporal overlap in patients' admission period) ward and hospital contact did not decrease during the study period. Indirect ward and hospital contact were identified as independent risk factors associated with clonal transmission. In conclusion, undetected CPE reservoirs continue to evade hospital infection prevention measures. New measures are needed to address plasmid-mediated transmission, which accounted for 50% of CPE dissemination.

Estimating the economic cost of carbapenem resistant Enterobacterales healthcare associated infections in Singapore acute-care hospitals.

Quantifying the costs of hospital associated infections (HAIs) caused by carbapenem-resistant Enterobacterales (CRE) can aid hospital decision makers in infection prevention and control decisions. We estimate the costs of a CRE HAI by infection type and the annual costs of CRE HAIs to acute-care hospitals in Singapore. We used tree diagrams to estimate the costs (in Singapore dollar) of different CRE HAI types from the health service perspective and compared them to the costs of carbapenem-susceptible HAIs. We used two approaches to estimate costs-direct costs of consumables for infection prevention and treatment; and costs associated with lost bed days. Cost of a HAI were extrapolated to annual CRE HAI incidence in Singapore acute-care hospitals to estimate the annual cost to the hospitals. We found that the cost of a CRE HAI based on direct cost and lost bed days are SGD$9,913 (95% CI, SGD$9,431-10,395) and SGD$10,044 (95% CI, SGD$9,789-10,300) respectively. CRE HAIs are markedly higher than the carbapenem-susceptible HAIs for all infection types. In both approaches, CRE pneumonia was the costliest infection. Based on a CRE HAI incidence of 233 per 100,000 inpatient admissions, CRE HAIs costed SGD$12.16M (95% CI, SGD$11.84-12.48M) annually based on direct costs, and SGD$12.33M (95% CI, SGD$12.01-12.64M) annually based on lost bed days. In conclusion, we described the cost of CRE HAIs in Singapore hospitals and identified infections with the highest costs. The findings may be useful in informing future economic evaluations of competing CRE HAI prevention and treatment programmes.

Diagnostics for COVID-19: moving from pandemic response to control.

Diagnostics have proven to be crucial to the COVID-19 pandemic response. There are three major methods for the detection of SARS-CoV-2 infection and their role has evolved during the course of the pandemic. Molecular tests such as PCR are highly sensitive and specific at detecting viral RNA, and are recommended by WHO for confirming diagnosis in individuals who are symptomatic and for activating public health measures. Antigen rapid detection tests detect viral proteins and, although they are less sensitive than molecular tests, have the advantages of being easier to do, giving a faster time to result, of being lower cost, and able to detect infection in those who are most likely to be at risk of transmitting the virus to others. Antigen rapid detection tests can be used as a public health tool for screening individuals at enhanced risk of infection, to protect people who are clinically vulnerable, to ensure safe travel and the resumption of schooling and social activities, and to enable economic recovery. With vaccine roll-out, antibody tests (which detect the host's response to infection or vaccination) can be useful surveillance tools to inform public policy, but should not be used to provide proof of immunity, as the correlates of protection remain unclear. All three types of COVID-19 test continue to have a crucial role in the transition from pandemic response to pandemic control.

Prioritising positively selected variants in whole-genome sequencing data using FineMAV.

BACKGROUND: In population genomics, polymorphisms that are highly differentiated between geographically separated populations are often suggestive of Darwinian positive selection. Genomic scans have highlighted several such regions in African and non-African populations, but only a handful of these have functional data that clearly associates candidate variations driving the selection process. Fine-Mapping of Adaptive Variation (FineMAV) was developed to address this in a high-throughput manner using population based whole-genome sequences generated by the 1000 Genomes Project. It pinpoints positively selected genetic variants in sequencing data by prioritizing high frequency, population-specific and functional derived alleles. RESULTS: We developed a stand-alone software that implements the FineMAV statistic. To graphically visualise the FineMAV scores, it outputs the statistics as bigWig files, which is a common file format supported by many genome browsers. It is available as a command-line and graphical user interface. The software was tested by replicating the FineMAV scores obtained using 1000 Genomes Project African, European, East and South Asian populations and subsequently applied to whole-genome sequencing datasets from Singapore and China to highlight population specific variants that can be subsequently modelled. The software tool is publicly available at https://github.com/fadilla-wahyudi/finemav . CONCLUSIONS: The software tool described here determines genome-wide FineMAV scores, using low or high-coverage whole-genome sequencing datasets, that can be used to prioritize a list of population specific, highly differentiated candidate variants for in vitro or in vivo functional screens. The tool displays these scores on the human genome browsers for easy visualisation, annotation and comparison between different genomic regions in worldwide human populations.

Impact of BMI and waist circumference on epigenome-wide DNA methylation and identification of epigenetic biomarkers in blood: an EWAS in multi-ethnic Asian individuals.

BACKGROUND: The prevalence of obesity and its related chronic diseases have been increasing especially in Asian countries. Obesity-related genetic variants have been identified, but these explain little of the variation in BMI. Recent studies reported associations between DNA methylation and obesity, mostly in non-Asian populations. METHODS: We performed an epigenome-wide association study (EWAS) on general adiposity (body mass index, BMI) and abdominal adiposity (waist circumference, WC) in 409 multi-ethnic Asian individuals and replicated BMI and waist-associated DNA methylation CpGs identified in other populations. The cross-lagged panel model and Mendelian randomization were used to assess the temporal relationship between methylation and BMI. The temporal relationship between the identified CpGs and inflammation and metabolic markers was also examined. RESULTS: EWAS identified 116 DNA methylation CpGs independently associated with BMI and eight independently associated with WC at false discovery rate PFDR < 0.05 in 409 Asian samples. We replicated 110 BMI-associated CpGs previously reported in Europeans and identified six novel BMI-associated CpGs and two novel WC-associated CpGs. We observed high consistency in association direction of effect compared to studies in other populations. Causal relationship analyses indicated that BMI was more likely to be the cause of DNA methylation alteration, rather than the consequence. The causal analyses using BMI-associated methylation risk score also suggested that higher levels of the inflammation marker IL-6 were likely the consequence of methylation change. CONCLUSION: Our study provides evidence of an association between obesity and DNA methylation in multi-ethnic Asians and suggests that obesity can drive methylation change. The results also suggested possible causal influence that obesity-related methylation changes might have on inflammation and lipoprotein levels.

New Polygenic Risk Score to Predict High Myopia in Singapore Chinese Children.

Purpose: The purpose of this study was to develop an Asian polygenic risk score (PRS) to predict high myopia (HM) in Chinese children in the Singapore Cohort of Risk factors for Myopia (SCORM) cohort. Methods: We included children followed from 6 to 11 years old until teenage years (12-18 years old). Cycloplegic autorefraction, ultrasound biometry, Illumina HumanHap 550, or 550 Duo Beadarrays, demographics, and environmental factors data were obtained. The PRS was generated from the Consortium for Refractive Error and Myopia genomewide association study (n = 542,934) and the Strabismus, Amblyopia, and Refractive Error in Singapore children Study (n = 500). The Growing Up in Singapore Towards healthy Outcomes Cohort study (n = 339) was the replication cohort. The outcome was teenage HM (≤ -5.00 D) with predictive performance assessed using the area under the curve (AUC). Results: Mean baseline age ± SD was 7.85 ± 0.84 (n = 1004) and 571 attended the teenage visit; 23.3% had HM. In multivariate analysis, the PRS was associated with a myopic spherical equivalent with an incremental R2 of 0.041 (95% confidence interval [CI] = 0.010, 0.073; P < 0.001). AUC for HM (0.77 [95% CI = 0.71-0.83]) performed better (P = 0.02) with the PRS compared with a model without (0.72 [95% CI = 0.65, 0.78]). Children at the top 25% PRS risk had a 2.34-fold-greater risk of HM (95% CI = 1.53, 3.55; P < 0.001). Conclusions: The new Asian PRS improved the predictive performance to detect children at risk of HM. Translational Relevance: Clinicians may use the PRS with other predictive factors to identify high risk children and guide interventions to reduce the risk of HM later in life.

Whole-genome analysis of drug-resistant Mycobacterium tuberculosis reveals novel mutations associated with fluoroquinolone resistance.

Multidrug-resistant and extensively drug-resistant tuberculosis (M/XDR-TB) remains a global public-health challenge. Known mutations in quinolone resistance-determination regions cannot fully explain phenotypic fluoroquinolone (FQ) resistance in Mycobacterium tuberculosis (Mtb). The aim of this study was to look for novel mutations in Mtb associated with resistance to FQ drugs using whole-genome sequencing analysis. Whole-genome sequences of 659 Mtb strains, including 214 with phenotypic FQ resistance and 445 pan-susceptible isolates, were explored for mutations associated with FQ resistance overall and with resistance to individual FQ drugs (ofloxacin, levofloxacin, moxifloxacin and gatifloxacin). Three novel genes (recC, Rv2005c and PPE59) associated with FQ resistance were identified (P < 0.00001 based on screening analysis and absence of relevant mutations in a pan-susceptible validation set of 360 strains). Nine novel single nucleotide polymorphisms (SNPs), including in gyrB (G5383A and G6773A), gyrA (G7892A), recC (G725900C and G726857T/C), Rv2005c (C2251373G, G2251420C and C2251725T) and PPE59 (C3847269T), were used for diagnostic performance analysis. Enhancing the known SNP set with five of these novel SNPs, including gyrA [G7892A (Leu247Leu)], recC [G725900C (Leu893Leu) and G726857T/C (Arg484Arg)], Rv2005c [G2251420C (Pro205Arg)] and PPE59 [C3847269T (Asn35Asn)] increased the sensitivity of detection of FQ-resistant Mtb from 83.2% (178/214) to 86.9% (186/214) while maintaining 100% specificity (360/360). No specific mutation associated with resistance to only a single drug (ofloxacin, levofloxacin, moxifloxacin or gatifloxacin) was found. In conclusion, this study reports possible additional mutations associated with FQ resistance in Mtb.

Cost effectiveness analysis of a polygenic risk tailored breast cancer screening programme in Singapore.

BACKGROUND: This study aimed to evaluate the cost-effectiveness of a breast cancer screening programme that incorporates genetic testing using breast cancer associated single nucleotide polymorphisms (SNPs), against the current biennial mammogram-only screening programme to aid in its implementation into the current programme in Singapore. METHODS: A Markov model was used to compare the costs and health outcomes of the current screening programme, against a polygenic risk-tailored screening programme, which can advise a long-term screening strategy depending on the individual's polygenic risk. The model took the perspective of the healthcare system, with a time horizon of 40 years, following women from the age of 35 to 74. Epidemiological and cost data were taken from Asian studies, and an annual discount rate of 3% was used. The model outcome was the incremental cost-effectiveness ratio (ICER), calculated from the difference in costs per quality-adjusted life year (QALY). Scenarios with varying risk thresholds for each polygenic risk group were examined. One-way and probabilistic sensitivity analyses were performed to assess parameter uncertainty. RESULTS: The ICER for a polygenic risk-tailored breast cancer screening programme, compared with the current biennial mammogram-only screening programme, was - 3713.80 SGD/QALY, with incremental costs < 0 and incremental effects > 0. The scenario analysis of different polygenic risk cutoffs showed that the ICERs remain negative, with all ICERs falling within the south-east quadrant of the cost-effectiveness plane, indicating that tailored screening is more cost effective than mammogram-only screening, with lower costs and higher QALYs to be gained. This suggests that a polygenic risk-tailored breast cancer screening programme is cost effective, entailing lower cost than the current mammogram-only programme, while causing no additional harm to women. CONCLUSION: Results from this cost-effectiveness analysis show that polygenic risk-tailored screening is cost effective with an ICER of - 3713.80 SGD/QALY. Tailored screening remains cost effective even across varying percentile cutoffs for each risk group. While the results look promising for incorporating polygenic risk into the current breast cancer screening programme, further studies should be conducted to address various limitations.

Childhood Factors and Dental Caries in the Permanent Dentition: Findings of an 8-Year Study Under a Nationwide School Dental Service.

OBJECTIVE: This study, using a retrospective analysis of nationwide dental records of school-going children over 8 years, aims to identify childhood factors associated with dental caries incidence in the permanent dentition up to adolescence. METHODS: Students studying in primary and secondary schools in Singapore are eligible for free basic dental care under the nationwide School Dental Service (SDS). All available dental records, general health records, and sociodemographic data from 2009 to 2017 were extracted as anonymised records, for a full cohort of 29,617 students that were enrolled in Primary 1 in 2009. Multiple logistic regression and modified Poisson regression were applied to identify risk factors for caries incidence in the permanent dentition over 8 years. RESULTS: Caries occurred in the permanent dentition of 9389 (31.7%) students in the 8 years. Risk of caries incidence in the permanent dentition was associated with baseline caries in the primary (risk ratio [RR]: 1.88; 95% CI: 1.81, 1.95) and permanent (RR: 1.54; 95% CI: 1.47, 1.61) dentition, tooth hypoplasia (RR: 1.58; 95% CI: 1.49, 1.68), and poor baseline oral hygiene (RR: 1.07; 95% CI: 1.03, 1.12). Sociodemographic predictors of caries incidence include a lower socioeconomic status, Chinese ethnicity, female gender, and enrolment in nonmainstream schools or schools in the Eastern and Western regions of Singapore. CONCLUSIONS: Both clinical and sociodemographic factors in childhood are associated with caries incidence in the permanent dentition and can be used for structuring dental service provision and identifying caries-susceptible individuals and groups for early prevention and intervention.

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.

Exploring Factors Underlying Ethnic Difference in Age-related Macular Degeneration Prevalence.

AIMS: To assess contributions of dietary and genetic factors to ethnic differences in AMD prevalence. DESIGN: Population-based analytical study. METHODS: In the Blue Mountains Eye Study, Australia (European ancestry n = 2826) and Multi-Ethnic Cohort Study, Singapore (Asian ancestry, n = 1900), AMD was assessed from retinal photographs. Patterns of dietary composition and scores of the Alternative Healthy Eating Index were computed using food frequency questionnaire data. Genetic susceptibility to AMD was determined using either single nucleotide polymorphisms (SNPs) of the complement factor H and age-related maculopathy susceptibility 2 genes, or combined odds-weighted genetic risk scores of 24 AMD-associated SNPs. Associations of AMD with ethnicity, diet, and genetics were assessed using logistic regression. Six potential mediators covering genetic, diet and lifestyle factors were assessed for their contributions to AMD risk difference between the two samples using mediation analyses. RESULTS: Age-standardized prevalence of any (early or late) AMD was higher in the European (16%) compared to Asian samples (9%, p < .01). Mean AMD-related genetic risk scores were also higher in European (33.3 ± 4.4) than Asian (Chinese) samples (31.7 ± 3.7, p < .001). In a model simultaneously adjusting for age, ethnicity, genetic susceptibility and Alternative Healthy Eating Index scores, only age and genetic susceptibility were significantly associated with AMD. Genetic risk scores contributed 19% of AMD risk difference between the two samples while intake of polyunsaturated fatty acids contributed 7.2%. CONCLUSION: Genetic susceptibility to AMD was higher in European compared to Chinese samples and explained more of the AMD risk difference between the two samples than the dietary factors investigated.

Genome-wide meta-analysis associates GPSM1 with type 2 diabetes, a plausible gene involved in skeletal muscle function.

Genome-wide association studies (GWASs) have identified many genetic variations associated with type 2 diabetes mellitus (T2DM) in Asians, but understanding the functional genetic variants that influence traits is often a complex process. In this study, fine mapping and other analytical strategies were performed to investigate the effects of G protein signaling modulator 1 (GPSM1) on insulin resistance in skeletal muscle. A total of 128 single-nucleotide polymorphisms (SNPs) within GPSM1 were analysed in 21,897 T2DM cases and 32,710 healthy controls from seven GWASs. The SNP rs28539249 in intron 9 of GPSM1 showed a nominally significant association with T2DM in Asians (OR = 1.07, 95% CI = 1.04-1.10, P < 10-4). The GPSM1 mRNA was increased in skeletal muscle and correlated with T2DM traits across obese mice model. An eQTL for the cis-acting regulation of GPSM1 expression in human skeletal muscle was identified for rs28539249, and the increased GPSM1 expression related with T2DM traits within GEO datasets. Another independent Asian cohort showed that rs28539249 is associated with the skeletal muscle expression of CACFD1, GTF3C5, SARDH, and FAM163B genes, which are functionally enriched for endoplasmic reticulum stress (ERS) and unfolded protein response (UPR) pathways. Moreover, rs28539249 locus was predicted to disrupt regulatory regions in human skeletal muscle with enriched epigenetic marks and binding affinity for CTCF. Supershift EMSA assays followed luciferase assays demonstrated the CTCF specifically binding to rs28539249-C allele leading to decreased transcriptional activity. Thus, the post-GWAS annotation confirmed the Asian-specific association of genetic variant in GPSM1 with T2DM, suggesting a role for the variant in the regulation in skeletal muscle.

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.

A trans-ancestral meta-analysis of genome-wide association studies reveals loci associated with childhood obesity.

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.