Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels.

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

First insight into oral microbiome diversity in Papua New Guineans reveals a specific regional signature.

The oral microbiota is a highly complex and diversified part of the human microbiome. Being located at the interface between the human body and the exterior environment, this microbiota can deepen our understanding of the environmental impacts on the global status of human health. This research topic has been well addressed in Westernized populations, but these populations only represent a fraction of human diversity. Papua New Guinea hosts very diverse environments and one of the most unique human biological diversities worldwide. In this study we performed the first known characterization of the oral microbiome in 85 Papua New Guinean individuals living in different environments, using a qualitative and quantitative approach. We found a significant geographical structure of the Papua New Guineans oral microbiome, especially in the groups most isolated from urban spaces. In comparison to other global populations, two bacterial genera related to iron absorption were significantly more abundant in Papua New Guineans and Aboriginal Australians, which suggests a shared oral microbiome signature. Further studies will be needed to confirm and explore this possible regional-specific oral microbiome profile.

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.

Chronology of natural selection in Oceanian genomes.

As human populations left Asia to first settle in Oceania around 50,000 years ago, they entered a territory ecologically separated from the Old World for millions of years. We analyzed genomic data of 239 modern Oceanian individuals to detect and date signals of selection specific to this region. Combining both relative and absolute dating approaches, we identified a strong selection pattern between 52,000 and 54,000 years ago in the genomes of descendants of the first settlers of Sahul. This strikingly corresponds to the dates of initial settlement as inferred from archaeological evidence. Loci under selection during this period, some showing enrichment in Denisovan ancestry, overlap genes involved in the immune response and diet, especially based on plants. Pathogens and natural resources, especially from endemic plants, therefore appear to have acted as strong selective pressures on the genomes of the first settlers of Sahul.

Episodes of Diversification and Isolation in Island Southeast Asian and Near Oceanian Male Lineages.

Island Southeast Asia (ISEA) and Oceania host one of the world's richest assemblages of human phenotypic, linguistic, and cultural diversity. Despite this, the region's male genetic lineages are globally among the last to remain unresolved. We compiled ∼9.7 Mb of Y chromosome (chrY) sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this data set allows us to fully resolve and date the regional chrY phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region's initial settlement ∼50 kya, and extensive expansions <6 kya. Notably, ∼40-25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in ISEA, New Guinea, and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ∼25 kya, at the time of the last glacial maximum. This improved chrY topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and ISEA, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the chrY presented here thus enables a detailed exploration of past isolation, interaction, and change in one of the world's least understood regions.

Papua New Guinean Genomes Reveal the Complex Settlement of North Sahul.

The settlement of Sahul, the lost continent of Oceania, remains one of the most ancient and debated human migrations. Modern New Guineans inherited a unique genetic diversity tracing back 50,000 years, and yet there is currently no model reconstructing their past population dynamics. We generated 58 new whole-genome sequences from Papua New Guinea, filling geographical gaps in previous sampling, specifically to address alternative scenarios of the initial migration to Sahul and the settlement of New Guinea. Here, we present the first genomic models for the settlement of northeast Sahul considering one or two migrations from Wallacea. Both models fit our data set, reinforcing the idea that ancestral groups to New Guinean and Indigenous Australians split early, potentially during their migration in Wallacea where the northern route could have been favored. The earliest period of human presence in Sahul was an era of interactions and gene flow between related but already differentiated groups, from whom all modern New Guineans, Bismarck islanders, and Indigenous Australians descend. The settlement of New Guinea was probably initiated from its southeast region, where the oldest archaeological sites have been found. This was followed by two migrations into the south and north lowlands that ultimately reached the west and east highlands. We also identify ancient gene flows between populations in New Guinea, Australia, East Indonesia, and the Bismarck Archipelago, emphasizing the fact that the anthropological landscape during the early period of Sahul settlement was highly dynamic rather than the traditional view of extensive isolation.

Genome skims analysis of betel palms (Areca spp., Arecaceae) and development of a profiling method to assess their plastome diversity.

The betel nut (Areca catechu L., Arecaceae) is a monoecious cultivated palm tree that is widespread in tropical regions. It is mainly cultivated for producing areca nuts, from which seeds are extracted and chewed by local populations principally in the Indo-Pacific region. Seeds contain alkaloids which are central nervous system stimulants and are highly addictive. Wild relatives of the betel nut are distributed in South Asia and Australasia, with ca. 40-50 Areca species currently recognized. The geographic origin(s) of the betel nut and its subsequent diffusion and diversification remains poorly documented. Here, a genome skimming approach was applied to screen nucleotidic variation in the most abundant genomic regions. Low coverage sequencing data allowed us to assemble full plastomes, mitochondrial regions (either full mitogenomes or the full set of mitochondrial genes) and the nuclear ribosomal DNA cluster for nine representatives of the Areca genus collected in the field and herbarium collections (including a 182-years old specimen collected during the Dumont d'Urville's expedition). These three genomic compartments provided similar phylogenetic signals, and revealed very low genomic diversity in our sample of cultivated betel nut. We finally developed a genotyping method targeting 34 plastid DNA microsatellites. This plastome profiling approach is useful for tracing the spread of matrilineages, and in combination with nuclear genomic data, can resolve the history of the betel nut. Our method also proves to be efficient for analyzing herbarium specimens, even those collected >100 years ago.

Phenotypic differences between highlanders and lowlanders in Papua New Guinea.

OBJECTIVES: Altitude is one of the most demanding environmental pressures for human populations. Highlanders from Asia, America and Africa have been shown to exhibit different biological adaptations, but Oceanian populations remain understudied [Woolcock et al., 1972; Cotes et al., 1974; Senn et al., 2010]. We tested the hypothesis that highlanders phenotypically differ from lowlanders in Papua New Guinea, as a result of inhabiting the highest mountains in Oceania for at least 20,000 years. MATERIALS AND METHODS: We collected data for 13 different phenotypes related to altitude for 162 Papua New Guineans living at high altitude (Mont Wilhelm, 2,300-2,700 m above sea level (a.s.l.) and low altitude (Daru, <100m a.s.l.). Multilinear regressions were performed to detect differences between highlanders and lowlanders for phenotypic measurements related to body proportions, pulmonary function, and the circulatory system. RESULTS: Six phenotypes were significantly different between Papua New Guinean highlanders and lowlanders. Highlanders show shorter height (p-value = 0.001), smaller waist circumference (p-value = 0.002), larger Forced Vital Capacity (FVC) (p-value = 0.008), larger maximal (p-value = 3.20e -4) and minimal chest depth (p-value = 2.37e -5) and higher haemoglobin concentration (p-value = 3.36e -4). DISCUSSION: Our study reports specific phenotypes in Papua New Guinean highlanders potentially related to altitude adaptation. Similar to other human groups adapted to high altitude, the evolutionary history of Papua New Guineans appears to have also followed an adaptive biological strategy for altitude.

Stochastic models support rapid peopling of Late Pleistocene Sahul.

The peopling of Sahul (the combined continent of Australia and New Guinea) represents the earliest continental migration and settlement event of solely anatomically modern humans, but its patterns and ecological drivers remain largely conceptual in the current literature. We present an advanced stochastic-ecological model to test the relative support for scenarios describing where and when the first humans entered Sahul, and their most probable routes of early settlement. The model supports a dominant entry via the northwest Sahul Shelf first, potentially followed by a second entry through New Guinea, with initial entry most consistent with 50,000 or 75,000 years ago based on comparison with bias-corrected archaeological map layers. The model's emergent properties predict that peopling of the entire continent occurred rapidly across all ecological environments within 156-208 human generations (4368-5599 years) and at a plausible rate of 0.71-0.92 km year-1. More broadly, our methods and approaches can readily inform other global migration debates, with results supporting an exit of anatomically modern humans from Africa 63,000-90,000 years ago, and the peopling of Eurasia in as little as 12,000-15,000 years via inland routes.

Papuan mitochondrial genomes and the settlement of Sahul.

New Guineans represent one of the oldest locally continuous populations outside Africa, harboring among the greatest linguistic and genetic diversity on the planet. Archeological and genetic evidence suggest that their ancestors reached Sahul (present day New Guinea and Australia) by at least 55,000 years ago (kya). However, little is known about this early settlement phase or subsequent dispersal and population structuring over the subsequent period of time. Here we report 379 complete Papuan mitochondrial genomes from across Papua New Guinea, which allow us to reconstruct the phylogenetic and phylogeographic history of northern Sahul. Our results support the arrival of two groups of settlers in Sahul within the same broad time window (50-65 kya), each carrying a different set of maternal lineages and settling Northern and Southern Sahul separately. Strong geographic structure in northern Sahul remains visible today, indicating limited dispersal over time despite major climatic, cultural, and historical changes. However, following a period of isolation lasting nearly 20 ky after initial settlement, environmental changes postdating the Last Glacial Maximum stimulated diversification of mtDNA lineages and greater interactions within and beyond Northern Sahul, to Southern Sahul, Wallacea and beyond. Later, in the Holocene, populations from New Guinea, in contrast to those of Australia, participated in early interactions with incoming Asian populations from Island Southeast Asia and continuing into Oceania.

Testing for Betel Nut Alkaloids in Hair of Papuans Abusers using UPLC-MS/MS and UPLC-Q-Tof-MS.

Betel nut is the fruit of Areca palm, growing in Papua New Guinea. Mixed with limestone and stick mustard, arecoline and guvacoline, which are present in betel nut, are hydrolyzed into arecaidine and guvacine, respectively. As part of the study on dietary habits of Papuans residents, our laboratory was asked to analyze the four alkaloids in hair to document long-term exposure. Hair samples were collected from 19 adult subjects (males = 11; females = 8), by some of the authors, and were sent to the laboratory for analysis. The four alkaloids have very similar chemical structures. In order to accurately identify the drugs, two methods were developed. First, the compounds were identified using an ultra-high-performance liquid chromatography system coupled to time-of-flight mass spectrometry. Then, they were quantified by an ultra-high-performance liquid chromatography system coupled to tandem mass spectrometry. After decontamination with dichloromethane, hair samples were cut into very small segments and 20 mg were incubated in methanol for 2 h 30 min in an ultrasound bath. After cooling, the methanol was evaporated to dryness in presence of 20-µL octanol to prevent volatilization. Nicotine-d4 was used as an internal standard. Linearity was observed for concentrations ranging from the limit of quantification to 20 ng/mg for arecoline, arecaidine, guvacine and guvacoline. Measured concentrations were in the range 60 pg/mg to 18 ng/mg for arecoline (n = 19), 14 pg/mg to 2.5 ng/mg for guvacoline (n = 11), 63 pg/mg to 3.8 ng/mg for arecaidine (n = 11) and 100 pg/mg to 3.2 ng/mg for guvacine (n = 6). There was no correlation between concentrations of arecoline and arecaidine (ratio from 0.01 to 0.18) and guvacoline and guvacine (ratio from 0.06 to 3.50). However, the identification of these substances in hair is a good marker of consumption of betel nut and allows us to document a local practice that remains difficult to evaluate just by questioning.

Disentangling Immediate Adaptive Introgression from Selection on Standing Introgressed Variation in Humans.

Recent studies have reported evidence suggesting that portions of contemporary human genomes introgressed from archaic hominin populations went to high frequencies due to positive selection. However, no study to date has specifically addressed the postintrogression population dynamics of these putative cases of adaptive introgression. Here, for the first time, we specifically define cases of immediate adaptive introgression (iAI) in which archaic haplotypes rose to high frequencies in humans as a result of a selective sweep that occurred shortly after the introgression event. We define these cases as distinct from instances of selection on standing introgressed variation (SI), in which an introgressed haplotype initially segregated neutrally and subsequently underwent positive selection. Using a geographically diverse data set, we report novel cases of selection on introgressed variation in living humans and shortlist among these cases those whose selective sweeps are more consistent with having been the product of iAI rather than SI. Many of these novel inferred iAI haplotypes have potential biological relevance, including three that contain immune-related genes in West Siberians, South Asians, and West Eurasians. Overall, our results suggest that iAI may not represent the full picture of positive selection on archaically introgressed haplotypes in humans and that more work needs to be done to analyze the role of SI in the archaic introgression landscape of living humans.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Evolution of the pygmy phenotype: evidence of positive selection fro genome-wide scans in African, Asian, and Melanesian pygmies.

Human pygmy populations inhabit different regions of the world, from Africa to Melanesia. In Asia, short-statured populations are often referred to as "negritos." Their short stature has been interpreted as a consequence of thermoregulatory, nutritional, and/or locomotory adaptations to life in tropical forests. A more recent hypothesis proposes that their stature is the outcome of a life history trade-off in high-mortality environments, where early reproduction is favored and, consequently, early sexual maturation and early growth cessation have coevolved. Some serological evidence of deficiencies in the growth hormone/insulin-like growth factor axis have been previously associated with pygmies' short stature. Using genome-wide single-nucleotide polymorphism genotype data, we first tested whether different negrito groups living in the Philippines and Papua New Guinea are closely related and then investigated genomic signals of recent positive selection in African, Asian, and Papuan pygmy populations. We found that negritos in the Philippines and Papua New Guinea are genetically more similar to their nonpygmy neighbors than to one another and have experienced positive selection at different genes. These results indicate that geographically distant pygmy groups are likely to have evolved their short stature independently. We also found that selection on common height variants is unlikely to explain their short stature and that different genes associated with growth, thyroid function, and sexual development are under selection in different pygmy groups.

Human adaptation and plant use in highland New Guinea 49,000 to 44,000 years ago.

After their emergence by 200,000 years before the present in Africa, modern humans colonized the globe, reaching Australia and New Guinea by 40,000 to 50,000 years ago. Understanding how humans lived and adapted to the range of environments in these areas has been difficult because well-preserved settlements are scarce. Data from the New Guinea Highlands (at an elevation of ~2000 meters) demonstrate the exploitation of the endemic nut Pandanus and yams in archaeological sites dated to 49,000 to 36,000 years ago, which are among the oldest human sites in this region. The sites also contain stone tools thought to be used to remove trees, which suggests that the early inhabitants cleared forest patches to promote the growth of useful plants.