Papuan mitochondrial genomes and the settlement of Sahul

Abrupt onset of intensive human occupation 44,000 years ago on the threshold of Sahul

Archaeological evidence attests multiple early dispersals of Homo sapiens out of Africa, but genetic evidence points to the primacy of a single dispersal 70-40 ka. Laili in Timor-Leste is on the southern dispersal route between Eurasia and Australasia and has the earliest record of human occupation in the eastern Wallacean archipelago. New evidence from the site shows that, unusually in the region, sediment accumulated in the shelter without human occupation, in the window 59-54 ka. This was followed by an abrupt onset of intensive human habitation beginning ~44 ka. The initial occupation is distinctive from overlying layers in the aquatic focus of faunal exploitation, while it has similarities in material culture to other early Homo sapiens sites in Wallacea. We suggest that the intensive early occupation at Laili represents a colonisation phase, which may have overwhelmed previous human dispersals in this part of the world.

Physiography, foraging mobility, and the first peopling of Sahul

The route and speed of migration into Sahul by Homo sapiens remain a major research question in archaeology. Here, we introduce an approach which models the impact of the physical environment on human mobility by combining time-evolving landscapes with Lévy walk foraging patterns, this latter accounting for a combination of short-distance steps and occasional longer moves that hunter-gatherers likely utilised for efficient exploration of new environments. Our results suggest a wave of dispersal radiating across Sahul following riverine corridors and coastlines. Estimated migration speeds, based on archaeological sites and predicted travelled distances, fall within previously reported range from Sahul and other regions. From our mechanistic movement simulations, we then analyse the likelihood of archaeological sites and highlight areas in Australia that hold archaeological potential. Our approach complements existing methods and provides interesting perspectives on the Pleistocene archaeology of Sahul that could be applied to other regions around the world.

Hepatitis B genotypes in Aboriginal and Torres Strait Islander Australians: correlation with clinical course and implications for management

BACKGROUND

The prevalence of chronic hepatitis B (CHB) in Aboriginal and Torres Strait Islander Australians in Far North Queensland (FNQ) is greater than twice that of the general Australian population. CHB is common in Torres Strait Islanders diagnosed with hepatocellular carcinoma (HCC) - and in Aboriginals with HCC living in the Northern Territory - however, Aboriginals diagnosed with HCC in FNQ very rarely have CHB. The explanation for this apparent disparity is uncertain.

AIMS

To determine the HBV genotypes in the FNQ Aboriginal and Torres Strait Islander population and their correlation with clinical phenotype.

METHODS

We determined the HBV genotype of Aboriginal and Torres Strait Islander Australians living with CHB in FNQ and correlated this with demographic and clinical findings.

RESULTS

134/197 (68%) enrolled individuals had a sufficient viral load for genotyping. All 40 people with HBV/D genotype had Aboriginal heritage, whereas 85/93 (91%) with HBV/C had Torres Strait Islander heritage (P < 0.0001). Individuals with HBV/D were younger than those with HBV/C (median (interquartile range) age: 43 (39-48) vs 53 (42-66) years, P = 0.0002). However, they were less likely to be HBeAg positive (1/40 (3%) vs 23/93 (25%), P = 0.001). All three HCCs developed in Torres Strait Islanders; two-thirds were infected with HBV/C14; genotyping was not possible in the other individual. All 10 diagnoses of cirrhosis occurred in Torres Strait Islanders, 6/10 were infected with HBV/C14, genotyping was not possible in the other four individuals.

CONCLUSIONS

HBV genotypes in Aboriginal and Torres Strait Islander Australians in FNQ differ markedly, which could explain the significant differences in the clinical phenotype in the two populations and might be used to inform cost-effective CHB care in the region.

Indigenous Australian genomes show deep structure and rich novel variation

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama-Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000-35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

South Asian maternal and paternal lineages in southern Thailand and the role of sex-biased admixture

Previous genome-wide studies have reported South Asian (SA) ancestry in several Mainland Southeast Asian (MSEA) populations; however, additional details concerning population history, in particular the role of sex-specific aspects of the SA admixture in MSEA populations can be addressed with uniparental markers. Here, we generated ∼2.3 mB sequences of the male-specific portions of the Y chromosome (MSY) of a Tai-Kadai (TK)-speaking Southern Thai group (SouthernThai_TK), and complete mitochondrial (mtDNA) genomes of the SouthernThai_TK and an Austronesian (AN)-speaking Southern Thai (SouthernThai_AN) group. We identified new mtDNA haplogroups, e.g. Q3, E1a1a1, B4a1a and M7c1c3 that have not previously reported in Thai populations, but are frequent in Island Southeast Asia and Oceania, suggesting interactions between MSEA and these regions. SA prevalent mtDNA haplogroups were observed at frequencies of ~35-45% in the Southern Thai groups; both of them showed more genetic relatedness to Austroasiatic (AA) speaking Mon than to any other group. For MSY, SouthernThai_TK had ~35% SA prevalent haplogroups and exhibited closer genetic affinity to Central Thais. We also analyzed published data from other MSEA populations and observed SA ancestry in some additional MSEA populations that also reflects sex-biased admixture; in general, most AA- and AN-speaking groups in MSEA were closer to SA than to TK groups based on mtDNA, but the opposite pattern was observed for the MSY. Overall, our results of new genetic lineages and sex-biased admixture from SA to MSEA groups attest to the additional value that uniparental markers can add to studies of genome-wide variation.

The first maternal genetic study of hunter-gatherers from Vietnam

The current limitation of ancient DNA data from Vietnam led to the controversy surrounding the prehistory of people in this region. The combination of high heat and humidity damaged ancient bones that challenged the study of human evolution, especially when using DNA as study materials. So far, only 4 k years of history have been recorded despite the 65 k years of history of anatomically modern human occupations in Vietnam. Here we report, to our knowledge, the oldest mitogenomes of two hunter-gatherers from Vietnam. We extracted DNA from the femurs of two individuals aged 6.2 k cal BP from the Con Co Ngua (CCN) site in Thanh Hoa, Vietnam. This archeological site is the largest cemetery of the hunter-gatherer population in Southeast Asia (SEA) that was discovered, but their genetics have not been explored until the present. We indicated that the CCN haplotype belongs to a rare haplogroup that was not detected in any present-day Vietnamese individuals. Further matrilineal analysis on CCN mitogenomes showed a close relationship with ancient farmers and present-day populations in SEA. The mitogenomes of hunter-gatherers from Vietnam debate the "two layers" model of peopling history in SEA and provide an alternative solution for studying challenging ancient human samples from Vietnam.

Interspecific introgression patterns reveal the origins of worldwide cultivated bananas in New Guinea

Hybridizations between Musa species and subspecies, enabled by their transport via human migration, were proposed to have played an important role in banana domestication. We exploited sequencing data of 226 Musaceae accessions, including wild and cultivated accessions, to characterize the inter(sub)specific hybridization pattern that gave rise to cultivated bananas. We identified 11 genetic pools that contributed to cultivars, including two contributors of unknown origin. Informative alleles for each of these genetic pools were pinpointed and used to obtain genome ancestry mosaics of accessions. Diploid and triploid cultivars had genome mosaics involving three up to possibly seven contributors. The simplest mosaics were found for some diploid cultivars from New Guinea, combining three contributors, i.e., banksii and zebrina representing Musa acuminata subspecies and, more unexpectedly, the New Guinean species Musa schizocarpa. Breakpoints of M. schizocarpa introgressions were found to be conserved between New Guinea cultivars and the other analyzed diploid and triploid cultivars. This suggests that plants bearing these M. schizocarpa introgressions were transported from New Guinea and gave rise to currently cultivated bananas. Many cultivars showed contrasted mosaics with predominant ancestry from their geographical origin across Southeast Asia to New Guinea. This revealed that further diversification occurred in different Southeast Asian regions through hybridization with other Musa (sub)species, including two unknown ancestors that we propose to be M. acuminata ssp. halabanensis and a yet to be characterized M. acuminata subspecies. These results highlighted a dynamic crop formation process that was initiated in New Guinea, with subsequent diversification throughout Southeast Asia.

Human Genetic Research in Wallacea and Sahul: Recent Findings and Future Prospects

Genomic sequence data from worldwide human populations have provided a range of novel insights into our shared ancestry and the historical migrations that have shaped our global genetic diversity. However, a comprehensive understanding of these fundamental questions has been impeded by the lack of inclusion of many Indigenous populations in genomic surveys, including those from the Wallacean archipelago (which comprises islands of present-day Indonesia located east and west of Wallace's and Lydekker's Lines, respectively) and the former continent of Sahul (which once combined New Guinea and Australia during lower sea levels in the Pleistocene). Notably, these regions have been important areas of human evolution throughout the Late Pleistocene, as documented by diverse fossil and archaeological records which attest to the regional presence of multiple hominin species prior to the arrival of anatomically modern human (AMH) migrants. In this review, we collate and discuss key findings from the past decade of population genetic and phylogeographic literature focussed on the hominin history in Wallacea and Sahul. Specifically, we examine the evidence for the timing and direction of the ancient AMH migratory movements and subsequent hominin mixing events, emphasising several novel but consistent results that have important implications for addressing these questions. Finally, we suggest potentially lucrative directions for future genetic research in this key region of human evolution.

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.

•

239 human genomes from both sides of the Wallacean ecogeographical barriers

•

Signals of selection are dated between -54,000 to -52,000 in modern Oceanian genomes

•

Genes related to immunity and diet were under strong selection

•

Denisovan introgressions participated to the genetic adaptations present in Oceanians

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.

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

We present a two-headed approach called Bayesian Integrated Coalescent Epoch PlotS (BICEPS) for efficient inference of coalescent epoch models. Firstly, we integrate out population size parameters, and secondly, we introduce a set of more powerful Markov chain Monte Carlo (MCMC) proposals for flexing and stretching trees. Even though population sizes are integrated out and not explicitly sampled through MCMC, we are still able to generate samples from the population size posteriors. This allows demographic reconstruction through time and estimating the timing and magnitude of population bottlenecks and full population histories. Altogether, BICEPS can be considered a more muscular version of the popular Bayesian skyline model. We demonstrate its power and correctness by a well-calibrated simulation study. Furthermore, we demonstrate with an application to SARS-CoV-2 genomic data that some analyses that have trouble converging with the traditional Bayesian skyline prior and standard MCMC proposals can do well with the BICEPS approach. BICEPS is available as open-source package for BEAST 2 under GPL license and has a user-friendly graphical user interface.[Bayesian phylogenetics; BEAST 2; BICEPS; coalescent model.].

Is there still evolution in the human population?

It is often claimed that humanity has stopped evolving because modern medicine erased all selection on survival. Even if that would be true, and it is not, there would be other mechanisms of evolution which could still led to changes in allelic frequencies. Here I show, by applying basic evolutionary genetics knowledge, that we expect humanity to evolve. The results from genome sequencing projects have repeatedly affirmed that there are still recent signs of selection in our genomes. I give some examples of such adaptation. Then I briefly discuss what our evolutionary future has in store for us.

Late Pleistocene/Early Holocene sites in the montane forests of New Guinea yield early record of cassowary hunting and egg harvesting

How early human foragers impacted insular forests is a topic with implications across multiple disciplines, including resource management. Paradoxically, terminal Pleistocene and Early Holocene impacts of foraging communities have been characterized as both extreme-as in debates over human-driven faunal extinctions-and minimal compared to later landscape transformations by farmers and herders. We investigated how rainforest hunter-gatherers managed resources in montane New Guinea and present some of the earliest documentation of Late Pleistocene through mid-Holocene exploitation of cassowaries (Aves: Casuariidae). Worldwide, most insular ratites were extirpated by the Late Holocene, following human arrivals, including elephant birds of Madagascar (Aepyornithidae) and moa of Aotearoa/New Zealand (Dinornithiformes)-icons of anthropogenic island devastation. Cassowaries are exceptional, however, with populations persisting in New Guinea and Australia. Little is known of past human exploitation and what factors contributed to their survival. We present a method for inferring past human interaction with mega-avifauna via analysis of microstructural features of archaeological eggshell. We then contextualize cassowary hunting and egg harvesting by montane foragers and discuss the implications of human exploitation. Our data suggest cassowary egg harvesting may have been more common than the harvesting of adults. Furthermore, our analysis of cassowary eggshell microstructural variation reveals a distinct pattern of harvesting eggs in late ontogenetic stages. Harvesting eggs in later stages of embryonic growth may reflect human dietary preferences and foraging seasonality, but the observed pattern also supports the possibility that-as early as the Late Pleistocene-people were collecting eggs in order to hatch and rear cassowary chicks.

Landscape rules predict optimal superhighways for the first peopling of Sahul

Archaeological data and demographic modelling suggest that the peopling of Sahul required substantial populations, occurred rapidly within a few thousand years and encompassed environments ranging from hyper-arid deserts to temperate uplands and tropical rainforests. How this migration occurred and how humans responded to the physical environments they encountered have, however, remained largely speculative. By constructing a high-resolution digital elevation model for Sahul and coupling it with fine-scale viewshed analysis of landscape prominence, least-cost pedestrian travel modelling and high-performance computing, we create over 125 billion potential migratory pathways, whereby the most parsimonious routes traversed emerge. Our analysis revealed several major pathways-superhighways-transecting the continent, that we evaluated using archaeological data. These results suggest that the earliest Australian ancestors adopted a set of fundamental rules shaped by physiological capacity, attraction to visually prominent landscape features and freshwater distribution to maximize survival, even without previous experience of the landscapes they encountered.

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 dataset, 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.

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.

Mitogenomes Reveal Two Major Influxes of Papuan Ancestry across Wallacea Following the Last Glacial Maximum and Austronesian Contact

The tropical archipelago of Wallacea contains thousands of individual islands interspersed between mainland Asia and Near Oceania, and marks the location of a series of ancient oceanic voyages leading to the peopling of Sahul—i.e., the former continent that joined Australia and New Guinea at a time of lowered sea level—by 50,000 years ago. Despite the apparent deep antiquity of human presence in Wallacea, prior population history research in this region has been hampered by patchy archaeological and genetic records and is largely concentrated upon more recent history that follows the arrival of Austronesian seafarers ~3000–4000 years ago (3–4 ka). To shed light on the deeper history of Wallacea and its connections with New Guinea and Australia, we performed phylogeographic analyses on 656 whole mitogenomes from these three regions, including 186 new samples from eight Wallacean islands and three West Papuan populations. Our results point to a surprisingly dynamic population history in Wallacea, marked by two periods of extensive demographic change concentrated around the Last Glacial Maximum ~15 ka and post-Austronesian contact ~3 ka. These changes appear to have greatly diminished genetic signals informative about the original peopling of Sahul, and have important implications for our current understanding of the population history of the region.

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.