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.

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.

Severity of gastric intestinal metaplasia predicts the risk of gastric cancer: a prospective multicentre cohort study (GCEP).

OBJECTIVE: To investigate the incidence of gastric cancer (GC) attributed to gastric intestinal metaplasia (IM), and validate the Operative Link on Gastric Intestinal Metaplasia (OLGIM) for targeted endoscopic surveillance in regions with low-intermediate incidence of GC. METHODS: A prospective, longitudinal and multicentre study was carried out in Singapore. The study participants comprised 2980 patients undergoing screening gastroscopy with standardised gastric mucosal sampling, from January 2004 and December 2010, with scheduled surveillance endoscopies at year 3 and 5. Participants were also matched against the National Registry of Diseases Office for missed diagnoses of early gastric neoplasia (EGN). RESULTS: There were 21 participants diagnosed with EGN. IM was a significant risk factor for EGN (adjusted-HR 5.36; 95% CI 1.51 to 19.0; p<0.01). The age-adjusted EGN incidence rates for patients with and without IM were 133.9 and 12.5 per 100 000 person-years. Participants with OLGIM stages III-IV were at greatest risk (adjusted-HR 20.7; 95% CI 5.04 to 85.6; p<0.01). More than half of the EGNs (n=4/7) attributed to baseline OLGIM III-IV developed within 2 years (range: 12.7-44.8 months). Serum trefoil factor 3 distinguishes (Area Under the Receiver Operating Characteristics 0.749) patients with OLGIM III-IV if they are negative for H. pylori. Participants with OLGIM II were also at significant risk of EGN (adjusted-HR 7.34; 95% CI 1.60 to 33.7; p=0.02). A significant smoking history further increases the risk of EGN among patients with OLGIM stages II-IV. CONCLUSIONS: We suggest a risk-stratified approach and recommend that high-risk patients (OLGIM III-IV) have endoscopic surveillance in 2 years, intermediate-risk patients (OLGIM II) in 5 years.

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.

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.

Lessons learnt from easing COVID-19 restrictions: an analysis of countries and regions in Asia Pacific and Europe.

The COVID-19 pandemic is an unprecedented global crisis. Many countries have implemented restrictions on population movement to slow the spread of severe acute respiratory syndrome coronavirus 2 and prevent health systems from becoming overwhelmed; some have instituted full or partial lockdowns. However, lockdowns and other extreme restrictions cannot be sustained for the long term in the hope that there will be an effective vaccine or treatment for COVID-19. Governments worldwide now face the common challenge of easing lockdowns and restrictions while balancing various health, social, and economic concerns. To facilitate cross-country learning, this Health Policy paper uses an adapted framework to examine the approaches taken by nine high-income countries and regions that have started to ease COVID-19 restrictions: five in the Asia Pacific region (ie, Hong Kong [Special Administrative Region], Japan, New Zealand, Singapore, and South Korea) and four in Europe (ie, Germany, Norway, Spain, and the UK). This comparative analysis presents important lessons to be learnt from the experiences of these countries and regions. Although the future of the virus is unknown at present, countries should continue to share their experiences, shield populations who are at risk, and suppress transmission to save lives.

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.

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.

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.

Evaluating the contribution of rare variants to type 2 diabetes and related traits using pedigrees.

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis-expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants.

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.

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.

Genomic prediction of tuberculosis drug-resistance: benchmarking existing databases and prediction algorithms.

BACKGROUND: It is possible to predict whether a tuberculosis (TB) patient will fail to respond to specific antibiotics by sequencing the genome of the infecting Mycobacterium tuberculosis (Mtb) and observing whether the pathogen carries specific mutations at drug-resistance sites. This advancement has led to the collation of TB databases such as PATRIC and ReSeqTB that possess both whole genome sequences and drug resistance phenotypes of infecting Mtb isolates. Bioinformatics tools have also been developed to predict drug resistance from whole genome sequencing (WGS) data. Here, we evaluate the performance of four popular tools (TBProfiler, MyKrobe, KvarQ, PhyResSE) with 6746 isolates compiled from publicly available databases, and subsequently identify highly probable phenotyping errors in the databases by genetically predicting the drug phenotypes using all four software. RESULTS: Our results show that these bioinformatics tools generally perform well in predicting the resistance status for two key first-line agents (isoniazid, rifampicin), but the accuracy is lower for second-line injectables and fluoroquinolones. The error rates in the databases are also non-trivial, reaching as high as 31.1% for prothionamide, and that phenotypes from ReSeqTB are more susceptible to errors. CONCLUSIONS: The good performance of the automated software for drug resistance prediction from TB WGS data shown in this study further substantiates the usefulness and promise of utilising genetic data to accurately profile TB drug resistance, thereby reducing misdiagnoses arising from error-prone culture-based drug susceptibility testing.

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.

Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis.

Despite the progress in human leukocyte antigen (HLA) causal variant mapping, independent localization of major histocompatibility complex (MHC) risk from classical HLA genes is challenging. Here, we conducted a large-scale MHC fine-mapping analysis of rheumatoid arthritis (RA) in a Japanese population (6,244 RA cases and 23,731 controls) population by using HLA imputation, followed by a multi-ethnic validation study including east Asian and European populations (n = 7,097 and 23,149, respectively). Our study identified an independent risk of a synonymous mutation at HLA-DOA, a non-classical HLA gene, on anti-citrullinated protein autoantibody (ACPA)-positive RA risk (p = 1.4 × 10(-9)), which demonstrated a cis-expression quantitative trait loci (cis-eQTL) effect on HLA-DOA expression. Trans-ethnic comparison revealed different linkage disequilibrium (LD) patterns in HLA-DOA and HLA-DRB1, explaining the observed HLA-DOA variant risk heterogeneity among ethnicities, which was most evident in the Japanese population. Although previous HLA fine-mapping studies have identified amino acid polymorphisms of the classical HLA genes as driving genetic susceptibility to disease, our study additionally identifies the dosage contribution of a non-classical HLA gene to disease etiology. Our study contributes to the understanding of HLA immunology in human diseases and suggests the value of incorporating additional ancestry in MHC fine-mapping.

Genome-wide association study identified copy number variants associated with sporadic colorectal cancer risk.

BACKGROUND: Multiple single nucleotide polymorphisms (SNPs) have been associated with colorectal cancer (CRC) risk. The role of structural or copy number variants (CNV) in CRC, however, remained unclear. We investigated the role of CNVs in patients with sporadic CRC. METHODS: A genome-wide association study (GWAS) was performed on 1000 Singapore Chinese patients aged 50 years or more with no family history of CRC and 1000 ethnicity-matched, age-matched and gender-matched healthy controls using the Affymetrix SNP 6 platform. After 16 principal component corrections, univariate and multivariate segmentations followed by association testing were performed on 1830 samples that passed quality assurance tests. RESULTS: A rare CNV region (CNVR) at chromosome 14q11 (OR=1.92 (95% CI 1.59 to 2.32), p=2.7e-12) encompassing CHD8, and common CNVR at chromosomes 3q13.12 (OR=1.54 (95% CI 1.33 to 1.77), p=2.9e-9) and 12p12.3 (OR=1.69 (95% CI 1.41 to 2.01), p=2.8e-9) encompassing CD47 and RERG/ARHGDIB, respectively, were significantly associated with CRC risk. CNV loci were validated in an independent replication panel using an optimised copy number assay. Whole-genome expression data in matched tumours of a subset of cases demonstrated that copy number loss at CHD8 was significantly associated with dysregulation of several genes that perturb the Wnt, TP53 and inflammatory pathways. CONCLUSIONS: A rare CNVR at 14q11 encompassing the chromatin modifier CHD8 was significantly associated with sporadic CRC risk. Copy number loss at CHD8 altered expressions of genes implicated in colorectal tumourigenesis.

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).

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

Southeast Asia (SEA) is enriched with a complex history of peopling. Malaysia, which is located at the crossroads of SEA, has been recognized as one of the hubs for early human migration. To unravel the genomic complexity of the native inhabitants of Malaysia, we sequenced 12 samples from 3 indigenous populations from Peninsular Malaysia and 4 native populations from North Borneo to a high coverage of 28-37×. We showed that the Negritos from Peninsular Malaysia shared a common ancestor with the East Asians, but exhibited some level of gene flow from South Asia, while the North Borneo populations exhibited closer genetic affinity towards East Asians than the Malays. The analysis of time of divergence suggested that ancestors of Negrito were the earliest settlers in the Malay Peninsula, whom first separated from the Papuans ~ 50-33 thousand years ago (kya), followed by East Asian (~ 40-15 kya), while the divergence time frame between North Borneo and East Asia populations predates the Austronesian expansion period implies a possible pre-Neolithic colonization. Substantial Neanderthal ancestry was confirmed in our genomes, as was observed in other East Asians. However, no significant difference was observed, in terms of the proportion of Denisovan gene flow into these native inhabitants from Malaysia. Judging from the similar amount of introgression in the Southeast Asians and East Asians, our findings suggest that the Denisovan gene flow may have occurred before the divergence of these populations and that the shared similarities are likely an ancestral component.