Genomic imputation of ancient Asian populations contrasts local adaptation in pre- and post-agricultural Japan.

Early modern humans lived as hunter-gatherers for millennia before agriculture, yet the genetic adaptations of these populations remain a mystery. Here, we investigate selection in the ancient hunter-gatherer-fisher Jomon and contrast pre- and post-agricultural adaptation in the Japanese archipelago. Building on the successful validation of imputation with ancient Asian genomes, we identify selection signatures in the Jomon, particularly robust signals from KITLG variants, which may have influenced dark pigmentation evolution. The Jomon lacks well-known adaptive variants (EDAR, ADH1B, and ALDH2), marking their emergence after the advent of farming in the archipelago. Notably, the EDAR and ADH1B variants were prevalent in the archipelago 1,300 years ago, whereas the ALDH2 variant could have emerged later due to its absence in other ancient genomes. Overall, our study underpins local adaptation unique to the Jomon population, which in turn sheds light on post-farming selection that continues to shape contemporary Asian populations.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation.

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

Ancient Genomes From Bronze Age Remains Reveal Deep Diversity and Recent Adaptive Episodes for Human Oral Pathobionts.

Ancient microbial genomes can illuminate pathobiont evolution across millenia, with teeth providing a rich substrate. However, the characterization of prehistoric oral pathobiont diversity is limited. In Europe, only preagricultural genomes have been subject to phylogenetic analysis, with none compared to more recent archaeological periods. Here, we report well-preserved microbiomes from two 4,000-year-old teeth from an Irish limestone cave. These contained bacteria implicated in periodontitis, as well as Streptococcus mutans, the major cause of caries and rare in the ancient genomic record. Despite deriving from the same individual, these teeth produced divergent Tannerella forsythia genomes, indicating higher levels of strain diversity in prehistoric populations. We find evidence of microbiome dysbiosis, with a disproportionate quantity of S. mutans sequences relative to other oral streptococci. This high abundance allowed for metagenomic assembly, resulting in its first reported ancient genome. Phylogenetic analysis indicates major postmedieval population expansions for both species, highlighting the inordinate impact of recent dietary changes. In T. forsythia, this expansion is associated with the replacement of older lineages, possibly reflecting a genome-wide selective sweep. Accordingly, we see dramatic changes in T. forsythia's virulence repertoire across this period. S. mutans shows a contrasting pattern, with deeply divergent lineages persisting in modern populations. This may be due to its highly recombining nature, allowing for maintenance of diversity through selective episodes. Nonetheless, an explosion in recent coalescences and significantly shorter branch lengths separating bacteriocin-carrying strains indicate major changes in S. mutans demography and function coinciding with sugar popularization during the industrial period.

Genomic insights into a tripartite ancestry in the Southern Ryukyu Islands.

A tripartite structure for the genetic origin of Japanese populations states that present-day populations are descended from three main ancestors: (1) the indigenous Jomon hunter-gatherers; (2) a Northeast Asian component that arrived during the agrarian Yayoi period; and (3) a major influx of East Asian ancestry in the imperial Kofun period. However, the genetic heterogeneity observed in different regions of the Japanese archipelago highlights the need to assess the applicability and suitability of this model. Here, we analyse historic genomes from the southern Ryukyu Islands, which have unique cultural and historical backgrounds compared with other parts of Japan. Our analysis supports the tripartite structure as the best fit in this region, with significantly higher estimated proportions of Jomon ancestry than mainland Japanese. Unlike the main islands, where each continental ancestor was directly brought by immigrants from the continent, those who already possessed the tripartite ancestor migrated to the southern Ryukyu Islands and admixed with the prehistoric people around the eleventh century AD, coinciding with the emergence of the Gusuku period. These results reaffirm the tripartite model in the southernmost extremes of the Japanese archipelago and show variability in how the structure emerged in diverse geographic regions.

Millennium-old pathogenic Mendelian mutation discovery for multiple osteochondromas from a Gaelic Medieval graveyard.

Only a limited number of genetic diseases are diagnosable in archaeological individuals and none have had causal mutations identified in genome-wide screens. Two individuals from the Gaelic Irish Medieval burial ground of Ballyhanna, Co. Donegal, showed evidence of bone tumours consistent with the autosomal dominant condition multiple osteochondromas. Genome sequencing of the earlier individual uncovered a missense mutation in the second exon of EXT1, a specific lesion that has been identified in several modern patients. The later individual lacked this but displayed a novel frameshift mutation leading to a premature stop codon and loss of function in the same gene. These molecular confirmations of a paleopathological diagnosis within a single rural ancient context are surprisingly disjunct, given the observation of clusters of this disease in modern isolated populations and a de novo mutation rate of only 10%.

Ancient Maltese genomes and the genetic geography of Neolithic Europe.

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

Ancient genomics reveals tripartite origins of Japanese populations.

Prehistoric Japan underwent rapid transformations in the past 3000 years, first from foraging to wet rice farming and then to state formation. A long-standing hypothesis posits that mainland Japanese populations derive dual ancestry from indigenous Jomon hunter-gatherer-fishers and succeeding Yayoi farmers. However, the genomic impact of agricultural migration and subsequent sociocultural changes remains unclear. We report 12 ancient Japanese genomes from pre- and postfarming periods. Our analysis finds that the Jomon maintained a small effective population size of ~1000 over several millennia, with a deep divergence from continental populations dated to 20,000 to 15,000 years ago, a period that saw the insularization of Japan through rising sea levels. Rice cultivation was introduced by people with Northeast Asian ancestry. Unexpectedly, we identify a later influx of East Asian ancestry during the imperial Kofun period. These three ancestral components continue to characterize present-day populations, supporting a tripartite model of Japanese genomic origins.

A curated dataset of modern and ancient high-coverage shotgun human genomes.

Over the last few years, genome-wide data for a large number of ancient human samples have been collected. Whilst datasets of captured SNPs have been collated, high coverage shotgun genomes (which are relatively few but allow certain types of analyses not possible with ascertained captured SNPs) have to be reprocessed by individual groups from raw reads. This task is computationally intensive. Here, we release a dataset including 35 whole-genome sequenced samples, previously published and distributed worldwide, together with the genetic pipeline used to process them. The dataset contains 72,041,355 sites called across 19 ancient and 16 modern individuals and includes sequence data from four previously published ancient samples which we sequenced to higher coverage (10-18x). Such a resource will allow researchers to analyse their new samples with the same genetic pipeline and directly compare them to the reference dataset without re-processing published samples. Moreover, this dataset can be easily expanded to increase the sample distribution both across time and space.

Inferring Population Histories for Ancient Genomes Using Genome-Wide Genealogies.

Ancient genomes anchor genealogies in directly observed historical genetic variation and contextualize ancestral lineages with archaeological insights into their geography and cultural associations. However, the majority of ancient genomes are of lower coverage and cannot be directly built into genealogies. Here, we present a fast and scalable method, Colate, the first approach for inferring ancestral relationships through time between low-coverage genomes without requiring phasing or imputation. Our approach leverages sharing patterns of mutations dated using a genealogy to infer coalescence rates. For deeply sequenced ancient genomes, we additionally introduce an extension of the Relate algorithm for joint inference of genealogies incorporating such genomes. Application to 278 present-day and 430 ancient DNA samples of >0.5x mean coverage allows us to identify dynamic population structure and directional gene flow between early farmer and European hunter-gatherer groups. We further show that the previously reported, but still unexplained, increase in the TCC/TTC mutation rate, which is strongest in West Eurasia today, was already present at similar strength and widespread in the Late Glacial Period ~10k-15k years ago, but is not observed in samples >30k years old. It is strongest in Neolithic farmers, and highly correlated with recent coalescence rates between other genomes and a 10,000-year-old Anatolian hunter-gatherer. This suggests gene-flow among ancient peoples postdating the last glacial maximum as widespread and localizes the driver of this mutational signal in both time and geography in that region. Our approach should be widely applicable in future for addressing other evolutionary questions, and in other species.

Evaluating genotype imputation pipeline for ultra-low coverage ancient genomes.

Although ancient DNA data have become increasingly more important in studies about past populations, it is often not feasible or practical to obtain high coverage genomes from poorly preserved samples. While methods of accurate genotype imputation from > 1 × coverage data have recently become a routine, a large proportion of ancient samples remain unusable for downstream analyses due to their low coverage. Here, we evaluate a two-step pipeline for the imputation of common variants in ancient genomes at 0.05-1 × coverage. We use the genotype likelihood input mode in Beagle and filter for confident genotypes as the input to impute missing genotypes. This procedure, when tested on ancient genomes, outperforms a single-step imputation from genotype likelihoods, suggesting that current genotype callers do not fully account for errors in ancient sequences and additional quality controls can be beneficial. We compared the effect of various genotype likelihood calling methods, post-calling, pre-imputation and post-imputation filters, different reference panels, as well as different imputation tools. In a Neolithic Hungarian genome, we obtain ~ 90% imputation accuracy for heterozygous common variants at coverage 0.05 × and > 97% accuracy at coverage 0.5 ×. We show that imputation can mitigate, though not eliminate reference bias in ultra-low coverage ancient genomes.

Population genomics of the Viking world.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent.

A dynastic elite in monumental Neolithic society.

The nature and distribution of political power in Europe during the Neolithic era remains poorly understood1. During this period, many societies began to invest heavily in building monuments, which suggests an increase in social organization. The scale and sophistication of megalithic architecture along the Atlantic seaboard, culminating in the great passage tomb complexes, is particularly impressive2. Although co-operative ideology has often been emphasised as a driver of megalith construction1, the human expenditure required to erect the largest monuments has led some researchers to emphasize hierarchy3-of which the most extreme case is a small elite marshalling the labour of the masses. Here we present evidence that a social stratum of this type was established during the Neolithic period in Ireland. We sampled 44 whole genomes, among which we identify the adult son of a first-degree incestuous union from remains that were discovered within the most elaborate recess of the Newgrange passage tomb. Socially sanctioned matings of this nature are very rare, and are documented almost exclusively among politico-religious elites4-specifically within polygynous and patrilineal royal families that are headed by god-kings5,6. We identify relatives of this individual within two other major complexes of passage tombs 150 km to the west of Newgrange, as well as dietary differences and fine-scale haplotypic structure (which is unprecedented in resolution for a prehistoric population) between passage tomb samples and the larger dataset, which together imply hierarchy. This elite emerged against a backdrop of rapid maritime colonization that displaced a unique Mesolithic isolate population, although we also detected rare Irish hunter-gatherer introgression within the Neolithic population.

Origins and genetic legacies of the Caribbean Taino.

The Caribbean was one of the last parts of the Americas to be settled by humans, but how and when the islands were first occupied remains a matter of debate. Ancient DNA can help answering these questions, but the work has been hampered by poor DNA preservation. We report the genome sequence of a 1,000-year-old Lucayan Taino individual recovered from the site of Preacher's Cave in the Bahamas. We sequenced her genome to 12.4-fold coverage and show that she is genetically most closely related to present-day Arawakan speakers from northern South America, suggesting that the ancestors of the Lucayans originated there. Further, we find no evidence for recent inbreeding or isolation in the ancient genome, suggesting that the Lucayans had a relatively large effective population size. Finally, we show that the native American components in some present-day Caribbean genomes are closely related to the ancient Taino, demonstrating an element of continuity between precontact populations and present-day Latino populations in the Caribbean.

Insular Celtic population structure and genomic footprints of migration.

Previous studies of the genetic landscape of Ireland have suggested homogeneity, with population substructure undetectable using single-marker methods. Here we have harnessed the haplotype-based method fineSTRUCTURE in an Irish genome-wide SNP dataset, identifying 23 discrete genetic clusters which segregate with geographical provenance. Cluster diversity is pronounced in the west of Ireland but reduced in the east where older structure has been eroded by historical migrations. Accordingly, when populations from the neighbouring island of Britain are included, a west-east cline of Celtic-British ancestry is revealed along with a particularly striking correlation between haplotypes and geography across both islands. A strong relationship is revealed between subsets of Northern Irish and Scottish populations, where discordant genetic and geographic affinities reflect major migrations in recent centuries. Additionally, Irish genetic proximity of all Scottish samples likely reflects older strata of communication across the narrowest inter-island crossing. Using GLOBETROTTER we detected Irish admixture signals from Britain and Europe and estimated dates for events consistent with the historical migrations of the Norse-Vikings, the Anglo-Normans and the British Plantations. The influence of the former is greater than previously estimated from Y chromosome haplotypes. In all, we paint a new picture of the genetic landscape of Ireland, revealing structure which should be considered in the design of studies examining rare genetic variation and its association with traits.

The population genomics of archaeological transition in west Iberia: Investigation of ancient substructure using imputation and haplotype-based methods.

We analyse new genomic data (0.05-2.95x) from 14 ancient individuals from Portugal distributed from the Middle Neolithic (4200-3500 BC) to the Middle Bronze Age (1740-1430 BC) and impute genomewide diploid genotypes in these together with published ancient Eurasians. While discontinuity is evident in the transition to agriculture across the region, sensitive haplotype-based analyses suggest a significant degree of local hunter-gatherer contribution to later Iberian Neolithic populations. A more subtle genetic influx is also apparent in the Bronze Age, detectable from analyses including haplotype sharing with both ancient and modern genomes, D-statistics and Y-chromosome lineages. However, the limited nature of this introgression contrasts with the major Steppe migration turnovers within third Millennium northern Europe and echoes the survival of non-Indo-European language in Iberia. Changes in genomic estimates of individual height across Europe are also associated with these major cultural transitions, and ancestral components continue to correlate with modern differences in stature.

Capturing goats: documenting two hundred years of mitochondrial DNA diversity among goat populations from Britain and Ireland.

The domestic goat (Capra hircus) plays a key role in global agriculture, being especially prized in regions of marginal pasture. However, the advent of industrialized breeding has seen a dramatic reduction in genetic diversity within commercial populations, while high extinction rates among feral herds have further depleted the reservoir of genetic variation available. Here, we present the first survey of whole mitochondrial genomic variation among the modern and historical goat populations of Britain and Ireland using a combination of mtDNA enrichment and high throughput sequencing. Fifteen historical taxidermy samples, representing the indigenous 'Old Goat' populations of the islands, were sequenced alongside five modern Irish dairy goats and four feral samples from endangered populations in western Ireland. Phylogenetic and network analyses of European mitochondrial variation revealed distinct groupings dominated by historical British and Irish samples, which demonstrate a degree of maternal genetic structure between the goats of insular and continental Europe. Several Irish modern feral samples also fall within these clusters, suggesting continuity between these dwindling populations and the ancestral 'Old Goats' of Ireland and Britain.