Antiquity and diversity of aboriginal Australian Y-chromosomes

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.

A custom hybridisation enrichment forensic intelligence panel to infer biogeographic ancestry, hair and eye colour, and Y chromosome lineage

Massively parallel sequencing can provide genetic data for hundreds to thousands of loci in a single assay for various types of forensic testing. However, available commercial kits require an initial PCR amplification of short-to-medium sized targets which limits their application for highly degraded DNA. Development and optimisation of large PCR multiplexes also prevents creation of custom panels that target different suites of markers for identity, biogeographic ancestry, phenotype, and lineage markers (Y-chromosome and mtDNA). Hybridisation enrichment, an alternative approach for target enrichment prior to sequencing, uses biotinylated probes to bind to target DNA and has proven successful on degraded and ancient DNA. We developed a customisable hybridisation capture method, that uses individually mixed baits to allow tailored and targeted enrichment to specific forensic questions of interest. To allow collection of forensic intelligence data, we assembled and tested a custom panel of hybridisation baits to infer biogeographic ancestry, hair and eye colour, and paternal lineage (and sex) on modern male and female samples with a range of self-declared ancestries and hair/eye colour combinations. The panel correctly estimated biogeographic ancestry in 9/12 samples (75%) but detected European admixture in three individuals from regions with admixed demographic history. Hair and eye colour were predicted correctly in 83% and 92% of samples respectively, where intermediate eye colour and blond hair were problematic to predict. Analysis of Y-chromosome SNPs correctly assigned sex and paternal haplogroups, the latter complementing and supporting biogeographic ancestry predictions. Overall, we demonstrate the utility of this hybridisation enrichment approach to forensic intelligence testing using a combined suite of biogeographic ancestry, phenotype, and Y-chromosome SNPs for comprehensive biological profiling.

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.

Unveiling forensically relevant biogeographic, phenotype and Y-chromosome SNP variation in Pakistani ethnic groups using a customized hybridisation enrichment forensic intelligence panel

Massively parallel sequencing following hybridisation enrichment provides new opportunities to obtain genetic data for various types of forensic testing and has proven successful on modern as well as degraded and ancient DNA. A customisable forensic intelligence panel that targeted 124 SNP markers (67 ancestry informative markers, 23 phenotype markers from the HIrisplex panel, and 35 Y-chromosome SNPs) was used to examine biogeographic ancestry, phenotype and sex and Y-lineage in samples from different ethnic populations of Pakistan including Pothwari, Gilgit, Baloach, Pathan, Kashmiri and Siraiki. Targeted sequencing and computational data analysis pipeline allowed filtering of variants across the targeted loci. Study samples showed an admixture between East Asian and European ancestry. Eye colour was predicted accurately based on the highest p-value giving overall prediction accuracy of 92.8%. Predictions were consistent with reported hair colour for all samples, using the combined highest p-value approach and step-wise model incorporating probability thresholds for light or dark shade. Y-SNPs were successfully recovered only from male samples which indicates the ability of this method to identify biological sex and allow inference of Y-haplogroup. Our results demonstrate practicality of using hybridisation enrichment and MPS to aid in human intelligence gathering and will open many insights into forensic research in South Asia.

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.

Exploitation of a novel slowly mutating Y-STRs set and evaluation of slowly mutating Y-STRs plus Y-SNPs typing strategy in forensic genetics and evolutionary research

The Y-chromosome short tandem repeats (Y-STRs) loci with different mutation rates existing in the Y chromosome non-recombination region (NRY) allow to be applied in human forensics, genealogical researches, historical investigations and evolutionary studies. Currently, there is a high demand for pedigree search to narrow the scope of crime investigations. However, the commonly used Y-STRs kits generally contain Y-STRs with high mutation rates that could cause individuals from the same pedigree to display different haplotypes. Herein, we put forward a new strategy of Slowly Mutating (SM) Y-STRs plus Y-SNPs typing, which could not only improve the resolution and accuracy of pedigree search, but also be applicable to evolutionary research. First, we developed a nine SM Y-STRs assay by evaluating their mutation rates in 210 pedigrees. Then the gene diversity and efficiency of the SM Y-STRs and 172 Y-SNPs sets were investigated by 2304 unrelated males from 24 populations. Furthermore, network and time estimation analyses were performed to evaluate the new strategy's capability to reconstruct phylogenetic tree and reliability to infer the time to the most recent common ancestor (TMRCA). The nine SM Y-STRs assay even had a higher resolution and a comparable capacity of revealing population genetic differentiation compared to 172 Y-SNPs system. This new strategy could optimize the phylogenetic tree generated by commonly used Y-STR panels and obtain a quite consistent time estimations with the published dating.

Comparison of Markov Chain Monte Carlo Software for the Evolutionary Analysis of Y-Chromosomal Microsatellite Data

The evolutionary analysis of genetic data is an important subject of modern bioscience, with practical applications in diverse fields. Parameters of interest in this context include effective population sizes, mutation rates, population growth rates and the times to most recent common ancestors. Studying Y-chromosomal microsatellite data, in particular, has proven useful to unravel the recent patrilineal history of Homo sapiens populations. We compared the individual analysis options and technical details of four software tools that are widely used for this purpose, namely BATWING, BEAST, IMa2 and LAMARC, all of which use Bayesian coalescent-based Markov chain Monte Carlo (MCMC) methods for parameter estimation. More specifically, we simulated datasets for either eight or 20 hypothetical Y-chromosomal microsatellites, assuming a mutation rate of 0.0030 per generation and a constant or exponentially increasing population size, and used these data to evaluate the parameter estimation capacity of each tool. The datasets comprised between 100 and 1000 samples. In addition to runtime, the practical utility of the tools of interest can also be expected to depend critically upon the convergence behavior of the actual MCMC implementation. In fact, we found that runtime increased, and convergence rate decreased, with increasing sample size as expected. BATWING performed best with respect to runtime and convergence behavior, but only supports simple evolutionary models. As regards the spectrum of evolutionary models covered, and also in terms of cross-platform usability, BEAST provided the greatest flexibility. Finally, IMa2 and LAMARC turned out best to incorporate elaborate migration models in the analysis process.

Analysis of the South Australian Aboriginal population using the Global AIMs Nano ancestry test

We investigate the ability of the 31 SNP loci in the Global AIMs Nano set to distinguish self-declared Australian Aboriginal individuals from European, Oceanic, African, Native American and East Asian populations. Human evolution suggests that Australian Aboriginal individuals came to Australia approximately 50 000 years ago, during the time it made up part of Sahul. Since then the colonisation of Australia by Europeans has meant significant admixture within the Australian Aboriginal population. These two events present themselves in our study with the Aboriginal population creating a continuous genetic cline between the Oceanic and European groups. We also assigned the Aboriginal individuals into their traditional regional groups to determine whether there was any ability to distinguish these from each other. We found similar results to studies using other markers, namely that the more remote regions (that have been less affected by admixture) diverged from the rest. Overall, we found the ability of the GNano system to differentiate self-declared Australian Aboriginal individuals was reasonable but had limitations that need to be recognised if these assignments are applied to unknown individuals.

Ancient nuclear genomes enable repatriation of Indigenous human remains

After European colonization, the ancestral remains of Indigenous people were often collected for scientific research or display in museum collections. For many decades, Indigenous people, including Native Americans and Aboriginal Australians, have fought for their return. However, many of these remains have no recorded provenance, making their repatriation very difficult or impossible. To determine whether DNA-based methods could resolve this important problem, we sequenced 10 nuclear genomes and 27 mitogenomes from ancient pre-European Aboriginal Australians (up to 1540 years before the present) of known provenance and compared them to 100 high-coverage contemporary Aboriginal Australian genomes, also of known provenance. We report substantial ancient population structure showing strong genetic affinities between ancient and contemporary Aboriginal Australian individuals from the same geographic location. Our findings demonstrate the feasibility of successfully identifying the origins of unprovenanced ancestral remains using genomic methods.

A standardized framework for representation of ancestry data in genomics studies, with application to the NHGRI-EBI GWAS Catalog

The accurate description of ancestry is essential to interpret, access, and integrate human genomics data, and to ensure that these benefit individuals from all ancestral backgrounds. However, there are no established guidelines for the representation of ancestry information. Here we describe a framework for the accurate and standardized description of sample ancestry, and validate it by application to the NHGRI-EBI GWAS Catalog. We confirm known biases and gaps in diversity, and find that African and Hispanic or Latin American ancestry populations contribute a disproportionately high number of associations. It is our hope that widespread adoption of this framework will lead to improved analysis, interpretation, and integration of human genomics data.

Aboriginal Australian mitochondrial genome variation - an increased understanding of population antiquity and diversity

Aboriginal Australians represent one of the oldest continuous cultures outside Africa, with evidence indicating that their ancestors arrived in the ancient landmass of Sahul (present-day New Guinea and Australia) ~55 thousand years ago. Genetic studies, though limited, have demonstrated both the uniqueness and antiquity of Aboriginal Australian genomes. We have further resolved known Aboriginal Australian mitochondrial haplogroups and discovered novel indigenous lineages by sequencing the mitogenomes of 127 contemporary Aboriginal Australians. In particular, the more common haplogroups observed in our dataset included M42a, M42c, S, P5 and P12, followed by rarer haplogroups M15, M16, N13, O, P3, P6 and P8. We propose some major phylogenetic rearrangements, such as in haplogroup P where we delinked P4a and P4b and redefined them as P4 (New Guinean) and P11 (Australian), respectively. Haplogroup P2b was identified as a novel clade potentially restricted to Torres Strait Islanders. Nearly all Aboriginal Australian mitochondrial haplogroups detected appear to be ancient, with no evidence of later introgression during the Holocene. Our findings greatly increase knowledge about the geographic distribution and phylogenetic structure of mitochondrial lineages that have survived in contemporary descendants of Australia’s first settlers.

Aboriginal mitogenomes reveal 50,000 years of regionalism in Australia

Aboriginal Australians represent one of the longest continuous cultural complexes known. Archaeological evidence indicates that Australia and New Guinea were initially settled approximately 50 thousand years ago (ka); however, little is known about the processes underlying the enormous linguistic and phenotypic diversity within Australia. Here we report 111 mitochondrial genomes (mitogenomes) from historical Aboriginal Australian hair samples, whose origins enable us to reconstruct Australian phylogeographic history before European settlement. Marked geographic patterns and deep splits across the major mitochondrial haplogroups imply that the settlement of Australia comprised a single, rapid migration along the east and west coasts that reached southern Australia by 49-45 ka. After continent-wide colonization, strong regional patterns developed and these have survived despite substantial climatic and cultural change during the late Pleistocene and Holocene epochs. Remarkably, we find evidence for the continuous presence of populations in discrete geographic areas dating back to around 50 ka, in agreement with the notable Aboriginal Australian cultural attachment to their country.

Mitochondrial DNA diversity of present-day Aboriginal Australians and implications for human evolution in Oceania

Aboriginal Australians are one of the more poorly studied populations from the standpoint of human evolution and genetic diversity. Thus, to investigate their genetic diversity, the possible date of their ancestors' arrival and their relationships with neighboring populations, we analyzed mitochondrial DNA (mtDNA) diversity in a large sample of Aboriginal Australians. Selected mtDNA single-nucleotide polymorphisms and the hypervariable segment haplotypes were analyzed in 594 Aboriginal Australians drawn from locations across the continent, chiefly from regions not previously sampled. Most (~78%) samples could be assigned to mtDNA haplogroups indigenous to Australia. The indigenous haplogroups were all ancient (with estimated ages >40 000 years) and geographically widespread across the continent. The most common haplogroup was P (44%) followed by S (23%) and M42a (9%). There was some geographic structure at the haplotype level. The estimated ages of the indigenous haplogroups range from 39 000 to 55 000 years, dates that fit well with the estimated date of colonization of Australia based on archeological evidence (~47 000 years ago). The distribution of mtDNA haplogroups in Australia and New Guinea supports the hypothesis that the ancestors of Aboriginal Australians entered Sahul through at least two entry points. The mtDNA data give no support to the hypothesis of secondary gene flow into Australia during the Holocene, but instead suggest long-term isolation of the continent.

Deep Roots for Aboriginal Australian Y Chromosomes

Australia was one of the earliest regions outside Africa to be colonized by fully modern humans, with archaeological evidence for human presence by 47,000 years ago (47 kya) widely accepted [1, 2]. However, the extent of subsequent human entry before the European colonial age is less clear. The dingo reached Australia about 4 kya, indirectly implying human contact, which some have linked to changes in language and stone tool technology to suggest substantial cultural changes at the same time [3]. Genetic data of two kinds have been proposed to support gene flow from the Indian subcontinent to Australia at this time, as well: first, signs of South Asian admixture in Aboriginal Australian genomes have been reported on the basis of genome-wide SNP data [4]; and second, a Y chromosome lineage designated haplogroup C^∗, present in both India and Australia, was estimated to have a most recent common ancestor around 5 kya and to have entered Australia from India [5]. Here, we sequence 13 Aboriginal Australian Y chromosomes to re-investigate their divergence times from Y chromosomes in other continents, including a comparison of Aboriginal Australian and South Asian haplogroup C chromosomes. We find divergence times dating back to ∼50 kya, thus excluding the Y chromosome as providing evidence for recent gene flow from India into Australia.

•

We have sequenced 13 Aboriginal Australian Y chromosomes

•

These diverged from Y chromosomes in other continents around 50,000 years ago

•

They diverged from Papua New Guinean Y chromosomes soon after this

•

We find no evidence for Holocene male gene flow to Australia from South Asia