Search

Summary Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations., Graphical abstract

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity on the region. In a comparison with a European dataset, we identify 166 uniquely Indonesia-specific eQTLs, highlighting the benefits of performing association studies on non-European populations. By combining local ancestry and archaic introgression inference eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits. Our findings illustrate how local ancestry and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations.

Indonesia is the world’s fourth most populous country, host to striking levels of human diversity, regional patterns of admixture, and varying degrees of introgression from both Neanderthals and Denisovans. However, it has been largely excluded from the human genomics sequencing boom of the last decade. To serve as a benchmark dataset of molecular phenotypes across the region, we generated genome-wide CpG methylation and gene expression measurements in over 100 individuals from three locations that capture the major genomic and geographical axes of diversity across the Indonesian archipelago. Investigating between- and within-island differences, we find up to 10.55% of tested genes are differentially expressed between the islands of Sumba and New Guinea. Variation in gene expression is closely associated with DNA methylation, with expression levels of 9.80% of genes correlating with nearby promoter CpG methylation, and many of these genes being differentially expressed between islands. Genes identified in our differential expression and methylation analyses are enriched in pathways involved in immunity, highlighting Indonesia's tropical role as a source of infectious disease diversity and the strong selective pressures these diseases have exerted on humans. Finally, we identify robust within-island variation in DNA methylation and gene expression, likely driven by fine-scale environmental differences across sampling sites. Together, these results strongly suggest complex relationships between DNA methylation, transcription, archaic hominin introgression and immunity, all jointly shaped by the environment. This has implications for the application of genomic medicine, both in critically understudied Indonesia and globally, and will allow a better understanding of the interacting roles of genomic and environmental factors shaping molecular and complex phenotypes., Author summary Understanding how gene expression varies across individuals and environments is fundamental for the success of molecular medicine. However, the vast majority of studies of gene expression, particularly in healthy individuals, have primarily involved populations of European descent. It is unclear whether these findings are generalizable to the rest of humanity. We examined gene expression and DNA methylation in 116 individuals across Indonesia, the world's fourth most populous country. Our three populations–Mentawai off the western coast of Sumatra, Sumba in central Indonesia, and a group of hunter-gatherers, the Korowai, living in the province of Indonesian New Guinea–span the main axes of genetic and geographic diversity in the region. Up to 10% of genes show differences in expression and methylation patterns between islands, with many involved in immune function. Variation between villages within each location can likely be attributed to small-scale environmental differences. These findings emphasise the need to consider a broad range of genetic and environmental backgrounds when examining how molecular patterns and processes are shared across human populations.

International audience; Genome sequences are known for two archaic hominins-Neanderthals and Denisovans-which interbred with anatomically modern humans as they dispersed out of Africa. We identified high-confidence archaic haplotypes in 161 new genomes spanning 14 island groups in Island Southeast Asia and New Guinea and found large stretches of DNA that are inconsistent with a single introgressing Denisovan origin. Instead, modern Papuans carry hundreds of gene variants from two deeply divergent Denisovan lineages that separated over 350 thousand years ago. Spatial and temporal structure among these lineages suggest that introgression from one of these Denisovan groups predominantly took place east of the Wallace line and continued until near the end of the Pleistocene. A third Denisovan lineage occurs in modern East Asians. This regional mosaic suggests considerable complexity in archaic contact, with modern humans interbreeding with multiple Denisovan groups that were geographically isolated from each other over deep evolutionary time.

The Austronesian dispersal across the Indonesian Ocean to Madagascar and the Comoros has been well documented, but in an unexplained anomaly, few to no traces have been found of the Austronesian expansion in East Africa or the Arabian Peninsula. To revisit this peculiarity, we surveyed the Western Indian Ocean rim populations to identify potential Austronesian genetic ancestry. We generated full mitochondrial DNA genomes and genome-wide genotyping data for these individuals and compared them with the Banjar, the Indonesian source population of the westward Austronesian dispersal. We find strong support for Asian genetic contributions to maternal lineages and autosomal variation in modern day Somalia and Yemen. Surprisingly, this input reveals two apparently different geographic origins and timings of admixture for the Austronesian contact; one at a very early phase (likely associated with the early Austronesian dispersals), and a later movement dating to the end of nineteenth century. These Austronesian gene flows come, respectively, from Madagascar and directly from an unidentified location in Island Southeast Asia. This result reveals a far more complex dynamic of Austronesian dispersals through the Western Indian Ocean than has previously been understood and suggests that Austronesian movements within the Indian Ocean may have been part of a lengthy process, probably continuing well into the modern era., GenoToul bioinformatics facility of Genopole Toulouse Midi-Pyrénées, France. We thank the CNRGH production team for genotyping data generation. This work was supported by French Ministry of Research grant ANR-14-CE31-0013-01 (OCEOADAPTO) to F.-X.R, the French Ministry of Foreign and European Affairs (French Archaeological Mission in Borneo (MAFBO) to F.-X.R.), and the French Embassy in Indonesia through its Cultural and Cooperation Services (Institut Français en Indonésie); a fellowship from the Alexander von Humboldt Foundation to M.P.C.; funding from the Max Planck Society to M.S.; and grants from COMPETE 2020 and a Fundação para a Ciência e a Tecnologia-funded project (POCI-01-0145-FEDER-016609) to V.F

At the dawn of the second millennium, the expansion of the Indian Ocean trading network aligned with the emergence of an outward-oriented community along the East African coast to create a cosmopolitan cultural and trading zone known as the Swahili Corridor. On the basis of analyses of new genome-wide genotyping data and uniparental data in 276 individuals from coastal Kenya and the Comoros islands, along with large-scale genetic datasets from the Indian Ocean rim, we reconstruct historical population dynamics to show that the Swahili Corridor is largely an eastern Bantu genetic continuum. Limited gene flows from the Middle East can be seen in Swahili and Comorian populations at dates corresponding to historically documented contacts. However, the main admixture event in southern insular populations, particularly Comorian and Malagasy groups, occurred with individuals from Island Southeast Asia as early as the 8th century, reflecting an earlier dispersal from this region. Remarkably, our results support recent archaeological and linguistic evidence-based suggestions that the Comoros archipelago was the earliest location of contact between Austronesian and African populations in the Swahili Corridor. We acknowledge support from the GenoToul bioinformatics facility of Genopole Toulouse Midi-Pyre´ne´es, France. This research was supported by French Ministry of Research grant ANR-14-CE31-0013-01 (OCEOADAPTO) to F.-X.R, a Rutherford Fellowship from the Royal Society of New Zealand (RDF-10- MAU-001), a fellowship from the Alexander von Humboldt Foundation to M.P.C., and grants from COMPETE 2020 and a Fundac¸a˜o para a Cieˆncia e a Tecnologia-funded project (POCI-01-0145-FEDER-016609) to V.F. We thank Kamil Merito for his help collecting East African samples. We thank all the local Swahili and Comorian communities who participated in this study.

The Indian Ocean has long been a hub of interacting human populations. Following land- and sea-based routes, trade drove cultural contacts between far-distant ethnic groups in Asia, India, the Middle East and Africa, creating one of the world’s first proto-globalized environments. However, the extent to which population mixing was mediated by trade is poorly understood. Reconstructing admixture times from genomic data in 3,006 individuals from 187 regional populations reveals a close association between bouts of human migration and trade volumes during the last 2,000 years across the Indian Ocean trading system. Temporal oscillations in trading activity match phases of contraction and expansion in migration, with high water marks following the expansion of the Silk Roads in the 5th century AD, the rise of maritime routes in the 11th century and a drastic restructuring of the trade network following the arrival of Europeans in the 16th century. The economic fluxes of the Indian Ocean trade network therefore directly shaped exchanges of genes, in addition to goods and concepts.

International audience; Although situated ∼400 km from the east coast of Africa, Madagascar exhibits cultural, linguistic, and genetic traits from both Southeast Asia and Eastern Africa. The settlement history remains contentious; we therefore used a grid-based approach to sample at high resolution the genomic diversity (including maternal lineages, paternal lineages, and genome-wide data) across 257 villages and 2,704 Malagasy individuals. We find a common Bantu and Austronesian descent for all Malagasy individuals with a limited paternal contribution from Europe and the Middle East. Admixture and demographic growth happened recently, suggesting a rapid settlement of Madagascar during the last millennium. However, the distribution of African and Asian ancestry across the island reveals that the admixture was sex biased and happened heterogeneously across Madagascar, suggesting independent colonization of Madagascar from Africa and Asia rather than settlement by an already admixed population. In addition, there are geographic influences on the present genomic diversity, independent of the admixture, showing that a few centuries is sufficient to produce detectable genetic structure in human populations.

International audience; The history of human settlement in Southeast Asia has been complex and involved several distinct dispersal events. Here, we report the analyses of 1825 individuals from Southeast Asia including new genome-wide genotype data for 146 individuals from three Mainland Southeast Asian (Burmese, Malay and Vietnamese) and four Island Southeast Asian (Dusun, Filipino, Kankanaey and Murut) populations. While confirming the presence of previously recognised major ancestry components in the Southeast Asian population structure, we highlight the Kankanaey Igorots from the highlands of the Philippine Mountain Province as likely the closest living representatives of the source population that may have given rise to the Austronesian expansion. This conclusion rests on independent evidence from various analyses of autosomal data and uniparental markers. Given the extensive presence of trade goods, cultural and linguistic evidence of Indian influence in Southeast Asia starting from 2.5 kya, we also detect traces of a South Asian signature in different populations in the region dating to the last couple of thousand years.

Western Eurasia, notably the Near East and South Asia (Indian sub-continent), has interacted with Indonesia through Indian Ocean trade (the Maritime Silk Route) for more than 2000 years. The Indianization, and later Islamization, of Indonesia was enacted largely through trading activities, but also spread with help from the many Indianized and Islamic kingdoms that reigned over parts of the Indonesian archipelago during this time. Western Eurasian interaction left behind not only imported trade goods and cultural features, but also genetic traces. To locate the primary areas of Western Eurasian genetic influence in Indonesia, we have assembled published uniparental genetic data from ∼2900 Indonesian individuals. Frequency distributions show that Western Eurasian paternal lineages are found more commonly than Western Eurasian maternal lineages. Furthermore, the origins of these paternal lineages are more diverse than the corresponding maternal lineages, predominantly tracing back to South West and South Asia, and the Indian sub-continent, respectively. Indianized kingdoms in the Indonesian archipelago likely played a major role in dispersing Western Eurasian lineages, as these kingdoms overlap geographically with the current distribution of individuals carrying Western Eurasian genetic markers. Our data highlight the important role of these Western Eurasian migrants in contributing to the complexity of genetic diversity across the Indonesian archipelago today.

The Austronesian expansion, one of the last major human migrations, influenced regions as distant as tropical Asia, Remote Oceania and Madagascar, off the east coast of Africa. The identity of the Asian groups that settled Madagascar is particularly mysterious. While language connects Madagascar to the Ma’anyan of southern Borneo, haploid genetic data are more ambiguous. Here, we screened genome-wide diversity in 211 individuals from the Ma’anyan and surrounding groups in southern Borneo. Surprisingly, the Ma’anyan are characterized by a distinct, high frequency genomic component that is not found in Malagasy. This novel genetic layer occurs at low levels across Island Southeast Asia and hints at a more complex model for the Austronesian expansion in this region. In contrast, Malagasy show genomic links to a range of Island Southeast Asian groups, particularly from southern Borneo, but do not have a clear genetic connection with the Ma’anyan despite the obvious linguistic association.

International audience; The appearance of people associated with the Lapita culture in the South Pacific around 3,000 years ago marked the beginning of the last major human dispersal to unpopulated lands. However, the relationship of these pioneers to the long-established Papuan people of the New Guinea region is unclear. Here we present genome-wide ancient DNA data from three individuals from Vanuatu (about 3,100-2,700 years before present) and one from Tonga (about 2,700-2,300 years before present), and analyse them with data from 778 present-day East Asians and Oceanians. Today, indigenous people of the South Pacific harbour a mixture of ancestry from Papuans and a population of East Asian origin that no longer exists in unmixed form, but is a match to the ancient individuals. Most analyses have interpreted the minimum of twenty-five per cent Papuan ancestry in the region today as evidence that the first humans to reach Remote Oceania, including Polynesia, were derived from population mixtures near New Guinea, before their further expansion into Remote Oceania. However, our finding that the ancient individuals had little to no Papuan ancestry implies that later human population movements spread Papuan ancestry through the South Pacific after the first peopling of the islands.

It is commonly thought that human genetic diversity in non-African populations was primarily shaped by a recent out-of-Africa dispersal. Here we present a large geographical Y chromosome study using 459 high coverage sequences, including 302 newly reported here. We date the Y chromosomal Most Recent Common Ancestor (MRCA) in Africa at 254 (95% CI 192- 307) kya and differentiation of African and non-African lineages 52-121 kya. The age estimates for major non-African founder haplogroups cluster closely in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. We find that the highest basal Y chromosome diversity outside Africa has been preserved in South and Southeast Asians while extant Y chromosome diversity in Europe and the Near East stems from a small number of mid-Holocene founders. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck, followed by a fast recovery in Old World populations dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Pathogens found within local environments are a major cause of morbidity and mortality. This is particularly true in Indonesia, where infectious diseases such as malaria or dengue are a significant part of the disease burden within the country. One way to strengthen the control of infectious diseases is through better surveillance, however unequal investment in medical funding throughout Indonesia, particularly in rural areas, has resulted in under-reporting of cases. Here, we use transcriptome data from 117 healthy individuals living on the islands of Mentawai, Sumba, and the Indonesian side of New Guinea Island to explore which pathogens are present within whole blood. We are able to detect a broad range of taxa within RNA-sequencing data generated from whole blood, including bacteria, viruses, archaea, and eukaryotes. Using independent component analysis, we find that two of these pathogens—Flaviviridae and Plasmodium—have the most noticeable effects on expression profiles. We also identify specific genes linked with Plasmodium and Flavivirus abundance and find that both of these infections are most pronounced in the easternmost island within our Indonesian dataset. This study provides a framework for novel applications of RNA-seq as surveillance and a better understanding of environmental contributors affecting gene expression within Indonesia.