Ancient genomes from Iceland reveal the making of a human population

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

Stable population structure in Europe since the Iron Age, despite high mobility

Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000-3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

Detection of chromosomal aneuploidy in ancient genomes

Ancient DNA is a valuable tool for investigating genetic and evolutionary history that can also provide detailed profiles of the lives of ancient individuals. In this study, we develop a generalised computational approach to detect aneuploidies (atypical autosomal and sex chromosome karyotypes) in the ancient genetic record and distinguish such karyotypes from contamination. We confirm that aneuploidies can be detected even in low-coverage genomes ( ~ 0.0001-fold), common in ancient DNA. We apply this method to ancient skeletal remains from Britain to document the first instance of mosaic Turner syndrome (45,X0/46,XX) in the ancient genetic record in an Iron Age individual sequenced to average 9-fold coverage, the earliest known incidence of an individual with a 47,XYY karyotype from the Early Medieval period, as well as individuals with Klinefelter (47,XXY) and Down syndrome (47,XY, + 21). Overall, our approach provides an accessible and automated framework allowing for the detection of individuals with aneuploidies, which extends previous binary approaches. This tool can facilitate the interpretation of burial context and living conditions, as well as elucidate past perceptions of biological sex and people with diverse biological traits.

Genetic Load and Adaptive Potential of a Recovered Avian Species that Narrowly Avoided Extinction

High genetic diversity is a good predictor of long-term population viability, yet some species persevere despite having low genetic diversity. Here we study the genomic erosion of the Seychelles paradise flycatcher (Terpsiphone corvina), a species that narrowly avoided extinction after having declined to 28 individuals in the 1960s. The species recovered unassisted to over 250 individuals in the 1990s and was downlisted from Critically Endangered to Vulnerable in the International Union for the Conservation of Nature Red List in 2020. By comparing historical, prebottleneck (130+ years old) and modern genomes, we uncovered a 10-fold loss of genetic diversity. Highly deleterious mutations were partly purged during the bottleneck, but mildly deleterious mutations accumulated. The genome shows signs of historical inbreeding during the bottleneck in the 1960s, but low levels of recent inbreeding after demographic recovery. Computer simulations suggest that the species long-term small Ne reduced the masked genetic load and made the species more resilient to inbreeding and extinction. However, the reduction in genetic diversity due to the chronically small Ne and the severe bottleneck is likely to have reduced the species adaptive potential to face environmental change, which together with a higher load, compromises its long-term population viability. Thus, small ancestral Ne offers short-term bottleneck resilience but hampers long-term adaptability to environmental shifts. In light of rapid global rates of population decline, our work shows that species can continue to suffer the effect of their decline even after recovery, highlighting the importance of considering genomic erosion and computer modeling in conservation assessments.

Genetic diversity and historical demography of underutilised goat breeds in North-Western Europe

In the last decade, several studies aimed at dissecting the genetic architecture of local small ruminant breeds to discover which variations are involved in the process of adaptation to environmental conditions, a topic that has acquired priority due to climate change. Considering that traditional breeds are a reservoir of such important genetic variation, improving the current knowledge about their genetic diversity and origin is the first step forward in designing sound conservation guidelines. The genetic composition of North-Western European archetypical goat breeds is still poorly exploited. In this study we aimed to fill this gap investigating goat breeds across Ireland and Scandinavia, including also some other potential continental sources of introgression. The PCA and Admixture analyses suggest a well-defined cluster that includes Norwegian and Swedish breeds, while the crossbred Danish landrace is far apart, and there appears to be a close relationship between the Irish and Saanen goats. In addition, both graph representation of historical relationships among populations and f4-ratio statistics suggest a certain degree of gene flow between the Norse and Atlantic landraces. Furthermore, we identify signs of ancient admixture events of Scandinavian origin in the Irish and in the Icelandic goats. The time when these migrations, and consequently the introgression, of Scandinavian-like alleles occurred, can be traced back to the Viking colonisation of these two isles during the Viking Age (793-1066 CE). The demographic analysis indicates a complicated history of these traditional breeds with signatures of bottleneck, inbreeding and crossbreeding with the improved breeds. Despite these recent demographic changes and the historical genetic background shaped by centuries of human-mediated gene flow, most of them maintained their genetic identity, becoming an irreplaceable genetic resource as well as a cultural heritage.

The genetic legacy of African Americans from Catoctin Furnace

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

Imputation of ancient human genomes

Due to postmortem DNA degradation and microbial colonization, most ancient genomes have low depth of coverage, hindering genotype calling. Genotype imputation can improve genotyping accuracy for low-coverage genomes. However, it is unknown how accurate ancient DNA imputation is and whether imputation introduces bias to downstream analyses. Here we re-sequence an ancient trio (mother, father, son) and downsample and impute a total of 43 ancient genomes, including 42 high-coverage (above 10x) genomes. We assess imputation accuracy across ancestries, time, depth of coverage, and sequencing technology. We find that ancient and modern DNA imputation accuracies are comparable. When downsampled at 1x, 36 of the 42 genomes are imputed with low error rates (below 5%) while African genomes have higher error rates. We validate imputation and phasing results using the ancient trio data and an orthogonal approach based on Mendel's rules of inheritance. We further compare the downstream analysis results between imputed and high-coverage genomes, notably principal component analysis, genetic clustering, and runs of homozygosity, observing similar results starting from 0.5x coverage, except for the African genomes. These results suggest that, for most populations and depths of coverage as low as 0.5x, imputation is a reliable method that can improve ancient DNA studies.

Across the North Sea and beyond

Increasing numbers of ancient genomes from the Viking period retrieved across the North Sea regions are revealing a complex layer of genetic ancestries and a past cosmopolitanism that was triggered by different mobility patterns.

The genetic history of Scandinavia from the Roman Iron Age to the present

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

Founder vs. non-founder BRCA1/2 pathogenic alleles: the analysis of Belarusian breast and ovarian cancer patients and review of other studies on ethnically homogenous populations

The spectrum of BRCA1/2 mutations demonstrates significant interethnic variations. We analyzed for the first time the entire BRCA1/2 coding region in 340 Belarusian cancer patients with clinical signs of BRCA1/2-related disease, including 168 women with bilateral and/or early-onset breast cancer (BC), 104 patients with ovarian cancer and 68 subjects with multiple primary malignancies involving BC and/or OC. BRCA1/2 pathogenic alleles were detected in 98 (29%) women, with 67 (68%) of these being represented by founder alleles. Systematic comparison with other relevant studies revealed that the founder effect observed in Belarus is among the highest estimates observed worldwide. These findings are surprising, given that the population of Belarus did not experience geographic or cultural isolation throughout history.

The population genomic legacy of the second plague pandemic

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.¹ It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).² Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17^th-19^th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large F_ST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

The genetic history of the Southern Arc: A bridge between West Asia and Europe

By sequencing 727 ancient individuals from the Southern Arc (Anatolia and its neighbors in Southeastern Europe and West Asia) over 10,000 years, we contextualize its Chalcolithic period and Bronze Age (about 5000 to 1000 BCE), when extensive gene flow entangled it with the Eurasian steppe. Two streams of migration transmitted Caucasus and Anatolian/Levantine ancestry northward, and the Yamnaya pastoralists, formed on the steppe, then spread southward into the Balkans and across the Caucasus into Armenia, where they left numerous patrilineal descendants. Anatolia was transformed by intra-West Asian gene flow, with negligible impact of the later Yamnaya migrations. This contrasts with all other regions where Indo-European languages were spoken, suggesting that the homeland of the Indo-Anatolian language family was in West Asia, with only secondary dispersals of non-Anatolian Indo-Europeans from the steppe.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

Functional Studies of Genetic Variants Associated with Human Diseases in Notch Signaling-Related Genes Using Drosophila

Rare variants in the many genes related to Notch signaling cause diverse Mendelian diseases that affect myriad organ systems. In addition, genome- and exome-wide association studies have linked common and rare variants in Notch-related genes to common diseases and phenotypic traits. Moreover, somatic mutations in these genes have been observed in many types of cancer, some of which are classified as oncogenic and others as tumor suppressive. While functional characterization of some of these variants has been performed through experimental studies, the number of "variants of unknown significance" identified in patients with diverse conditions keeps increasing as high-throughput sequencing technologies become more commonly used in the clinic. Furthermore, as disease gene discovery efforts identify rare variants in human genes that have yet to be linked to a disease, the demand for functional characterization of variants in these "genes of unknown significance" continues to increase. In this chapter, we describe a workflow to functionally characterize a rare variant in a Notch signaling related gene that was found to be associated with late-onset Alzheimer's disease. This pipeline involves informatic analysis of the variant of interest using diverse human and model organism databases, followed by in vivo experiments in the fruit fly Drosophila melanogaster. The protocol described here can be used to study variants that affect amino acids that are not conserved between human and fly. By "humanizing" the almondex gene in Drosophila with mutant alleles and heterologous genomic rescue constructs, a missense variant in TM2D3 (TM2 Domain Containing 3) was shown to be functionally damaging. This, and similar approaches, greatly facilitate functional interpretations of genetic variants in the human genome and propel personalized medicine.

Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia

Large-scale movement of organisms across their habitable range, or migration, is an important evolutionary process that can shape genetic diversity and influence the adaptive spread of alleles. Although human migrations have been studied in great detail with modern and ancient genomes, recent anthropogenic influence on reducing the biogeographical constraints on the migration of nonnative species has presented opportunities in several study systems to ask the questions about how repeated introductions shape genetic diversity in the introduced range. We present an extensive overview of population structure of North American Arabidopsis thaliana by studying a set of 500 whole-genome sequenced and over 2,800 RAD-seq genotyped individuals in the context of global diversity represented by Afro-Eurasian genomes. We use methods based on haplotype and rare-allele sharing as well as phylogenetic modeling to identify likely sources of introductions of extant N. American A. thaliana from the native range in Africa and Eurasia. We find evidence of admixture among the introduced lineages having increased haplotype diversity and reduced mutational load. We also detect signals of selection in immune-system-related genes that may impart qualitative disease resistance to pathogens of bacterial and oomycete origin. We conclude that multiple introductions to a nonnative range can rapidly enhance the adaptive potential of a colonizing species by increasing haplotypic diversity through admixture. Our results lay the foundation for further investigations into the functional significance of admixture.

Genomic transformation and social organization during the Copper Age-Bronze Age transition in southern Iberia

The emerging Bronze Age (BA) of southeastern Iberia saw marked social changes. Late Copper Age (CA) settlements were abandoned in favor of hilltop sites, and collective graves were largely replaced by single or double burials with often distinctive grave goods indirectly reflecting a hierarchical social organization, as exemplified by the BA El Argar group. We explored this transition from a genomic viewpoint by tripling the amount of data available for this period. Concomitant with the rise of El Argar starting ~2200 cal BCE, we observe a complete turnover of Y-chromosome lineages along with the arrival of steppe-related ancestry. This pattern is consistent with a founder effect in male lineages, supported by our finding that males shared more relatives at sites than females. However, simple two-source models do not find support in some El Argar groups, suggesting additional genetic contributions from the Mediterranean that could predate the BA.

Fitness and Power: The Contribution of Genetics to the History of Differential Reproduction

Textual evidence from pre-modern societies supports the prediction that status differences among men translate to variance in reproductive success. In recent years, analysis of genetic data has opened up new ways of studying this relationship. By investigating cases that range over several millennia, these analyses repeatedly document the replacement of local men by newcomers and reveal instances of exceptional reproductive success of specific male lineages. These findings suggest that violent population transfers and conquests could generate considerable reproductive advantages for male dominants. At the same time, this does not always seem to have been the case. Moreover, it is difficult to link such outcomes to particular historical characters or events, or to identify status-biased reproductive inequalities within dominant groups. The proximate factors that mediated implied imbalances in reproductive success often remain unclear. A better understanding of the complex interplay between social power and genetic fitness will only arise from sustained transdisciplinary engagement.

The Genetic Analyses of French Canadians of Quebec Facilitate the Characterization of New Cancer Predisposing Genes Implicated in Hereditary Breast and/or Ovarian Cancer Syndrome Families

The French Canadian population of the province of Quebec has been recognized for its contribution to research in medical genetics, especially in defining the role of heritable pathogenic variants in cancer predisposing genes. Multiple carriers of a limited number of pathogenic variants in BRCA1 and BRCA2, the major risk genes for hereditary breast and/or ovarian cancer syndrome families, have been identified in French Canadians, which is in stark contrast to the array of over 2000 different pathogenic variants reported in each of these genes in other populations. As not all such cancer syndrome families are explained by BRCA1 and BRCA2, newly proposed gene candidates identified in other populations have been investigated for their role in conferring risk in French Canadian cancer families. For example, multiple carriers of distinct variants were identified in PALB2 and RAD51D. The unique genetic architecture of French Canadians has been attributed to shared ancestry due to common ancestors of early settlers of this population with origins mainly from France. In this review, we discuss the merits of genetically characterizing cancer predisposing genes in French Canadians of Quebec. We focused on genes that have been implicated in hereditary breast and/or ovarian cancer syndrome families as they have been the most thoroughly characterized cancer syndromes in this population. We describe how genetic analyses of French Canadians have facilitated: (i) the classification of variants in BRCA1 and BRCA2; (ii) the identification and classification of variants in newly proposed breast and/or ovarian cancer predisposing genes; and (iii) the identification of a new breast cancer predisposing gene candidate, RECQL. The genetic architecture of French Canadians provides a unique opportunity to evaluate new candidate cancer predisposing genes regardless of the population in which they were identified.

Genomics of Endometriosis: From Genome Wide Association Studies to Exome Sequencing

This review aims at better understanding the genetics of endometriosis. Endometriosis is a frequent feminine disease, affecting up to 10% of women, and characterized by pain and infertility. In the most accepted hypothesis, endometriosis is caused by the implantation of uterine tissue at ectopic abdominal places, originating from retrograde menses. Despite the obvious genetic complexity of the disease, analysis of sibs has allowed heritability estimation of endometriosis at ~50%. From 2010, large Genome Wide Association Studies (GWAS), aimed at identifying the genes and loci underlying this genetic determinism. Some of these loci were confirmed in other populations and replication studies, some new loci were also found through meta-analyses using pooled samples. For two loci on chromosomes 1 (near CCD42) and chromosome 9 (near CDKN2A), functional explanations of the SNP (Single Nucleotide Polymorphism) effects have been more thoroughly studied. While a handful of chromosome regions and genes have clearly been identified and statistically demonstrated as at-risk for the disease, only a small part of the heritability is explained (missing heritability). Some attempts of exome sequencing started to identify additional genes from families or populations, but are still scarce. The solution may reside inside a combined effort: increasing the size of the GWAS designs, better categorize the clinical forms of the disease before analyzing genome-wide polymorphisms, and generalizing exome sequencing ventures. We try here to provide a vision of what we have and what we should obtain to completely elucidate the genetics of this complex disease.

The genetic structure of Norway

The aim of the present study was to describe the genetic structure of the Norwegian population using genotypes from 6369 unrelated individuals with detailed information about places of residence. Using standard single marker- and haplotype-based approaches, we report evidence of two regions with distinctive patterns of genetic variation, one in the far northeast, and another in the south of Norway, as indicated by fixation indices, haplotype sharing, homozygosity, and effective population size. We detect and quantify a component of Uralic Sami ancestry that is enriched in the North. On a finer scale, we find that rates of migration have been affected by topography like mountain ridges. In the broader Scandinavian context, we detect elevated relatedness between the mid- and northern border areas towards Sweden. The main finding of this study is that despite Norway's long maritime history and as a former Danish territory, the region closest to mainland Europe in the south appears to have been an isolated region in Norway, highlighting the open sea as a barrier to gene flow into Norway.

Molecular genetics of inherited retinal degenerations in Icelandic patients

The study objective was to delineate the genetics of inherited retinal degenerations (IRDs) in Iceland, a small nation of 364.000 and a genetic isolate. Benefits include delineating novel pathogenic genetic variants and defining genetically homogenous patients as potential investigative molecular therapy candidates. The study sample comprised patients with IRD in Iceland ascertained through national centralized genetic and ophthalmological services at Landspitali, a national social support institute, and the Icelandic patient association. Information on patients' disease, syndrome, and genetic testing was collected in a clinical registry. Variants were reevaluated according to ACMG/AMP guidelines. Overall, 140 IRD patients were identified (point prevalence of 1/2.600), of which 70 patients had a genetic evaluation where two-thirds had an identified genetic cause. Thirteen disease genes were found in patients with retinitis pigmentosa, with the RLBP1 gene most common (n = 4). The c.1073 + 5G > A variant in the PRPF31 gene was homozygous in two RP patients. All tested patients with X-linked retinoschisis (XLRS) had the same possibly unique RS1 pathogenic variant, c.441G > A (p.Trp147X). Pathologic variants and genes for IRDs in Iceland did not resemble those described in ancestral North-Western European nations. Four variants were reclassified as likely pathogenic. One novel pathogenic variant defined a genetically homogenous XLRS patient group.

The composition of the founding population of Iceland: A new perspective from 3D analyses of basicranial shape

The settlement of Iceland in the Viking Age has been the focus of much research, but the composition of the founding population remains the subject of debate. Some lines of evidence suggest that almost all the founding population were Scandinavian, while others indicate a mix of Scandinavians and people of Scottish and Irish ancestry. To explore this issue further, we used three-dimensional techniques to compare the basicrania of skeletons from archaeological sites in Iceland, Scandinavia, and the British Isles. Our analyses yielded two main results. One was that the founding population likely consisted of roughly equal numbers of Scandinavians and people from the British Isles. The other was that the immigrants who originated from the British Isles included individuals of southern British ancestry as well as individuals of Scottish and Irish ancestry. The first of these findings is consistent with the results of recent analyses of modern and ancient DNA, while the second is novel. Our study, therefore, strengthens the idea that the founding population was a mix of Scandinavians and people from the British Isles, but also raises a new possibility regarding the regions from which the settlers originated.

Whole genome sequencing data of multiple individuals of Pakistani descent

Here we report whole genome sequencing of four individuals (H3, H4, H5, and H6) from a family of Pakistani descent. Whole genome sequencing yielded 1084.92, 894.73, 1068.62, and 1005.77 million mapped reads corresponding to 162.73, 134.21, 160.29, and 150.86 Gb sequence data and 52.49x, 43.29x, 51.70x, and 48.66x average coverage for H3, H4, H5, and H6, respectively. We identified 3,529,659, 3,478,495, 3,407,895, and 3,426,862 variants in the genomes of H3, H4, H5, and H6, respectively, including 1,668,024 variants common in the four genomes. Further, we identified 42,422, 39,824, 28,599, and 35,206 novel variants in the genomes of H3, H4, H5, and H6, respectively. A major fraction of the variants identified in the four genomes reside within the intergenic regions of the genome. Single nucleotide polymorphism (SNP) genotype based comparative analysis with ethnic populations of 1000 Genomes database linked the ancestry of all four genomes with the South Asian populations, which was further supported by mitochondria based haplogroup analysis. In conclusion, we report whole genome sequencing of four individuals of Pakistani descent.

Measurement(s)	SNV • genome
Technology Type(s)	whole genome sequencing • DNA sequencing
Factor Type(s)	individual
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	Pakistan

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12642761

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 history^1,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.

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 history^1,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.

The mitogenome portrait of Umbria in Central Italy as depicted by contemporary inhabitants and pre-Roman remains

Umbria is located in Central Italy and took the name from its ancient inhabitants, the Umbri, whose origins are still debated. Here, we investigated the mitochondrial DNA (mtDNA) variation of 545 present-day Umbrians (with 198 entire mitogenomes) and 28 pre-Roman individuals (obtaining 19 ancient mtDNAs) excavated from the necropolis of Plestia. We found a rather homogeneous distribution of western Eurasian lineages across the region, with few notable exceptions. Contemporary inhabitants of the eastern part, delimited by the Tiber River and the Apennine Mountains, manifest a peculiar mitochondrial proximity to central-eastern Europeans, mainly due to haplogroups U4 and U5a, and an overrepresentation of J (30%) similar to the pre-Roman remains, also excavated in East Umbria. Local genetic continuities are further attested to by six terminal branches (H1e1, J1c3, J2b1, U2e2a, U8b1b1 and K1a4a) shared between ancient and modern mitogenomes. Eventually, we identified multiple inputs from various population sources that likely shaped the mitochondrial gene pool of ancient Umbri over time, since early Neolithic, including gene flows with central-eastern Europe. This diachronic mtDNA portrait of Umbria fits well with the genome-wide population structure identified on the entire peninsula and with historical sources that list the Umbri among the most ancient Italic populations.

Latest trends in archaeogenetic research of west Eurasians

During the past ten years, archaeogenetic research has exponentially grown to study the genetic history of human populations, using genome-wide data from large numbers of ancient individuals. Of the entire globe, Europe and the Near East are the regions where ancient DNA data is by far most abundant with over 2500 genomes published at present. In this review, we focus on archaeological contexts that have received less attention in the literature, specifically those associated with west Eurasian hunter-gatherers as well as populations from the Iron Age and later historical periods. In addition, we emphasize a recent shift from continent-wide to regional and even site-specific studies, which is starting to provide novel insights into sociocultural aspects of past societies.

Ancient genome-wide DNA from France highlights the complexity of interactions between Mesolithic hunter-gatherers and Neolithic farmers

Starting from 12,000 years ago in the Middle East, the Neolithic lifestyle spread across Europe via separate continental and Mediterranean routes. Genomes from early European farmers have shown a clear Near Eastern/Anatolian genetic affinity with limited contribution from hunter-gatherers. However, no genomic data are available from modern-day France, where both routes converged, as evidenced by a mosaic cultural pattern. Here, we present genome-wide data from 101 individuals from 12 sites covering today's France and Germany from the Mesolithic (N = 3) to the Neolithic (N = 98) (7000-3000 BCE). Using the genetic substructure observed in European hunter-gatherers, we characterize diverse patterns of admixture in different regions, consistent with both routes of expansion. Early western European farmers show a higher proportion of distinctly western hunter-gatherer ancestry compared to central/southeastern farmers. Our data highlight the complexity of the biological interactions during the Neolithic expansion by revealing major regional variations.

Mapping co-ancestry connections between the genome of a Medieval individual and modern Europeans

Historical genetic links among similar populations can be difficult to establish. Identity by descent (IBD) analyses find genomic blocks that represent direct genealogical relationships among individuals. However, this method has rarely been applied to ancient genomes because IBD stretches are progressively fragmented by recombination and thus not recognizable after few tens of generations. To explore such genealogical relationships, we estimated long IBD blocks among modern Europeans, generating networks to uncover the genetic structures. We found that Basques, Sardinians, Icelanders and Orcadians form, each of them, highly intraconnected sub-clusters in a European network, indicating dense genealogical links within small, isolated populations. We also exposed individual genealogical links -such as the connection between one Basque and one Icelandic individual- that cannot be uncovered with other, widely used population genetics methods such as PCA or ADMIXTURE. Moreover, using ancient DNA technology we sequenced a Late Medieval individual (Barcelona, Spain) to high genomic coverage and identified IBD blocks shared between her and modern Europeans. The Medieval IBD blocks are statistically overrepresented only in modern Spaniards, which is the geographically closest population. This approach can be used to produce a fine-scale reflection of shared ancestry across different populations of the world, offering a direct genetic link from the past to the present.

The GenomeAsia 100K Project enables genetic discoveries across Asia

The underrepresentation of non-Europeans in human genetic studies so far has limited the diversity of individuals in genomic datasets and led to reduced medical relevance for a large proportion of the world’s population. Population-specific reference genome datasets as well as genome-wide association studies in diverse populations are needed to address this issue. Here we describe the pilot phase of the GenomeAsia 100K Project. This includes a whole-genome sequencing reference dataset from 1,739 individuals of 219 population groups and 64 countries across Asia. We catalogue genetic variation, population structure, disease associations and founder effects. We also explore the use of this dataset in imputation, to facilitate genetic studies in populations across Asia and worldwide.

Using whole-genome sequencing data from 1,739 individuals, the GenomeAsia 100K Project catalogues genetic variation, population structure and disease associations to facilitate genetic studies in Asian populations and increase representation in genetics studies worldwide.

The genetic landscape of Scotland and the Isles

Britain and Ireland are known to show population genetic structure; however, large swathes of Scotland, in particular, have yet to be described. Delineating the structure and ancestry of these populations will allow variant discovery efforts to focus efficiently on areas not represented in existing cohorts. Thus, we assembled genotype data for 2,554 individuals from across the entire archipelago with geographically restricted ancestry, and performed population structure analyses and comparisons to ancient DNA. Extensive geographic structuring is revealed, from broad scales such as a NE to SW divide in mainland Scotland, through to the finest scale observed to date: across 3 km in the Northern Isles. Many genetic boundaries are consistent with Dark Age kingdoms of Gaels, Picts, Britons, and Norse. Populations in the Hebrides, the Highlands, Argyll, Donegal, and the Isle of Man show characteristics of isolation. We document a pole of Norwegian ancestry in the north of the archipelago (reaching 23 to 28% in Shetland) which complements previously described poles of Germanic ancestry in the east, and "Celtic" to the west. This modern genetic structure suggests a northwestern British or Irish source population for the ancient Gaels that contributed to the founding of Iceland. As rarer variants, often with larger effect sizes, become the focus of complex trait genetics, more diverse rural cohorts may be required to optimize discoveries in British and Irish populations and their considerable global diaspora.

Sequence variation at ANAPC1 accounts for 24% of the variability in corneal endothelial cell density

The corneal endothelium is vital for transparency and proper hydration of the cornea. Here, we conduct a genome-wide association study of corneal endothelial cell density (cells/mm2), coefficient of cell size variation (CV), percentage of hexagonal cells (HEX) and central corneal thickness (CCT) in 6,125 Icelanders and find associations at 10 loci, including 7 novel. We assess the effects of these variants on various ocular biomechanics such as corneal hysteresis (CH), as well as eye diseases such as glaucoma and corneal dystrophies. Most notably, an intergenic variant close to ANAPC1 (rs78658973[A], frequency = 28.3%) strongly associates with decreased cell density and accounts for 24% of the population variance in cell density (β = -0.77 SD, P = 1.8 × 10-314) and associates with increased CH (β = 0.19 SD, P = 2.6 × 10-19) without affecting risk of corneal diseases and glaucoma. Our findings indicate that despite correlations between cell density and eye diseases, low cell density does not increase the risk of disease.