Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries

Methodological opportunities in genomic data analysis to advance health equity

The causes and consequences of inequities in genomic research and medicine are complex and widespread. However, it is widely acknowledged that underrepresentation of diverse populations in human genetics research risks exacerbating existing health disparities. Efforts to improve diversity are ongoing, but an often-overlooked source of inequity is the choice of analytical methods used to process, analyse and interpret genomic data. This choice can influence all areas of genomic research, from genome-wide association studies and polygenic score development to variant prioritization and functional genomics. New statistical and machine learning techniques to understand, quantify and correct for the impact of biases in genomic data are emerging within the wider genomic research and genomic medicine ecosystems. At this crucial time point, it is important to clarify where improvements in methods and practices can, or cannot, have a role in improving equity in genomics. Here, we review existing approaches to promote equity and fairness in statistical analysis for genomics, and propose future methodological developments that are likely to yield the most impact for equity.

Archaic hominin admixture and its consequences for modern humans

As anatomically modern humans dispersed out of Africa, they encountered and mated with now extinct hominins, including Neanderthals and Denisovans. It is now well established that all non-African individuals derive approximately 2% of their genome from Neanderthal ancestors and individuals of Melanesian and Australian aboriginal ancestry inherited an additional 2%-5% of their genomes from Denisovan ancestors. Attention has started to shift from documenting amounts of archaic admixture and identifying introgressed segments to understanding their molecular, phenotypic, and evolutionary consequences and refining models of human history. Here, we review recent insights into admixture between modern and archaic humans, emphasizing methodological innovations and the functional and phenotypic effects Neanderthal and Denisovan sequences have in contemporary individuals.

A scoping review of statistical methods to investigate colocalization between genetic associations and microRNA expression in osteoarthritis

Background

Genetic colocalization analysis is a statistical method that evaluates whether two traits (e.g., osteoarthritis [OA] risk and microRNA [miRNA] expression levels) share the same or distinct genetic association signals in a locus typically identified in genome-wide association studies (GWAS). This method is useful for providing insights into the biological relevance of genetic association signals, particularly in intergenic regions, which can help to elucidate disease mechanisms in OA and other complex traits.

Objectives

To review the existing literature on genetic colocalization methods, assess their suitability for studying OA, and investigate their capacity to integrate miRNA data, while bearing in view their statistical assumptions.

Design

We followed scoping review methodology and used Covidence software for data management. Search terms for colocalization, GWAS, and genetic or statistical models were used in the databases MEDLINE and EMBASE, searched till March 4, 2024.

Results

Our search returned 546 peer-reviewed papers, of which 96 were included following title/abstract and full-text screening. Based on both cumulative and annual publication counts, the most cited method for colocalization analysis was coloc. Four papers examined OA-related phenotypes, and none examined miRNA. An approach to colocalization analysis using miRNA was postulated based on further hand-searching.

Conclusions

Colocalization analysis is a largely unexplored method in OA. Many of the approaches to colocalization analysis identified in this review, including the integration of GWAS and miRNA data, may help to elucidate genetic and epigenetic factors implicated in OA and other complex traits.

Profiling genetically driven alternative splicing across the Indonesian archipelago

One of the regulatory mechanisms influencing the functional capacity of genes is alternative splicing (AS). Previous studies exploring the splicing landscape of human tissues have shown that AS has contributed to human biology, especially in disease progression and the immune response. Nonetheless, this phenomenon remains poorly characterized across human populations, and it is unclear how genetic and environmental variation contribute to AS. Here, we examine a set of 115 Indonesian samples from three traditional island populations spanning the genetic ancestry cline that characterizes Island Southeast Asia. We conduct a global AS analysis between islands to ascertain the degree of functionally significant AS events and their consequences. Using an event-based statistical model, we detected over 1,500 significant differential AS events across all comparisons. Additionally, we identify over 6,000 genetic variants associated with changes in splicing (splicing quantitative trait loci [sQTLs]), some of which are driven by Papuan-like genetic ancestry, and only show partial overlap with other publicly available sQTL datasets derived from other populations. Computational predictions of RNA binding activity reveal that a fraction of these sQTLs directly modulate the binding propensity of proteins involved in the splicing regulation of immune genes. Overall, these results contribute toward elucidating the role of genetic variation in shaping gene regulation in one of the most diverse regions in the world.

Exploring Genetic Variants and Platinum Chemotherapy Response in Indonesian Non-Small Cell Lung Cancer Patients: Insights from ERCC2 rs13181

Purpose

Individual responses to platinum-based treatment for Non-Small Cell Lung Cancer (NSCLC) are influenced by genetic polymorphisms, including Single Nucleotide Polymorphisms (SNPs). This study aimed to explore the role of ERCC2 in the Nucleotide Excision Repair (NER) pathway for platinum-based chemotherapy in NSCLC. While ERCC2 is widely studied, data for Southeast Asian populations are lacking. Addressing this gap could improve personalized treatment strategies for NSCLC in this demographic.

Patients and Methods

This study recruited 82 NSCLC patients with wildtype mutations of EGFR at Dr. H.A. Rotinsulu Lung Hospital, Bandung, and Dharmais Cancer Hospital, Jakarta. Data were collected prospectively from whole blood samples and medical records, while the effectiveness of chemotherapy was assessed by evaluating the response using RECIST 1.1 criteria on fourth cycle of chemotherapy.

Results

The results of this study showed the presence of genotype variation among the subjects, with frequency distribution as follows: AA genotype (82.9%), AC genotype (15.9%), and CC genotype (1.2%). The analysis of the association between ERCC2 rs13181 CC + AC versus AA with RECIST 1.1 yielded an odds ratio (OR) of 1.042 (95% CI: 0.292-3.715; p=0.950). A multivariate analysis that included cancer stage and chemotherapy regimen as additional variables produced an adjusted odds ratio (aOR) of 0.970 (95% CI: 0.263-3.568; p=0.963).

Conclusion

This study did not find statistically significant associations between ERCC2 rs13181 polymorphisms and chemotherapy responses. However, this research highlights the presence of genetic variation within the Indonesian population, with the AA genotype being the most prevalent, which may influence chemotherapy responses. The results provided preliminary data and lay the foundation for future comprehensive cohort observational investigations.

Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations

Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population

BACKGROUND

Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource.

RESULTS

We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue.

CONCLUSIONS

By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Synthetic biology in Indonesia: Potential and projection in a country with mega biodiversity

Synthetic biology has gained many interest around the globe in the last two decades, not only due to its rapid development but also the potential to provide addressable solutions using standardized design of biological systems. Considering its huge population, biodiversity, and natural resources, Indonesia could play an important role in shaping the future of synthetic biology towards a sustainable bio-circular economy. Here, we provide an overview of synthetic biology development in Indonesia, especially on exploring the potential of our biodiversity. We also discuss some potentials of synthetic biology in solving national issues. Furthermore, we also provide the projection and future landscape of synthetic biology development in Indonesia. In addition, we briefly explain the potential challenges that may arise during the development.

Denisovan introgression has shaped the immune system of present-day Papuans

Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.