The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Medieval DNA from Soqotra points to Eurasian origins of an isolated population at the crossroads of Africa and Arabia

Soqotra, an island situated at the mouth of the Gulf of Aden in the northwest Indian Ocean between Africa and Arabia, is home to ~60,000 people subsisting through fishing and semi-nomadic pastoralism who speak a Modern South Arabian language. Most of what is known about Soqotri history derives from writings of foreign travellers who provided little detail about local people, and the geographic origins and genetic affinities of early Soqotri people has not yet been investigated directly. Here we report genome-wide data from 39 individuals who lived between ~650 and 1750 CE at six locations across the island and document strong genetic connections between Soqotra and the similarly isolated Hadramawt region of coastal South Arabia that likely reflects a source for the peopling of Soqotra. Medieval Soqotri can be modelled as deriving ~86% of their ancestry from a population such as that found in the Hadramawt today, with the remaining ~14% best proxied by an Iranian-related source with up to 2% ancestry from the Indian sub-continent, possibly reflecting genetic exchanges that occurred along with archaeologically documented trade from these regions. In contrast to all other genotyped populations of the Arabian Peninsula, genome-level analysis of the medieval Soqotri is consistent with no sub-Saharan African admixture dating to the Holocene. The deep ancestry of people from medieval Soqotra and the Hadramawt is also unique in deriving less from early Holocene Levantine farmers and more from groups such as Late Pleistocene hunter-gatherers from the Levant (Natufians) than other mainland Arabians. This attests to migrations by early farmers having less impact in southernmost Arabia and Soqotra and provides compelling evidence that there has not been complete population replacement between the Pleistocene and Holocene throughout the Arabian Peninsula. Medieval Soqotra harboured a small population that showed qualitatively different marriage practices from modern Soqotri, with first-cousin unions occurring significantly less frequently than today.

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.

Accurate detection of identity-by-descent segments in human ancient DNA

Long DNA segments shared between two individuals, known as identity-by-descent (IBD), reveal recent genealogical connections. Here we introduce ancIBD, a method for identifying IBD segments in ancient human DNA (aDNA) using a hidden Markov model and imputed genotype probabilities. We demonstrate that ancIBD accurately identifies IBD segments >8 cM for aDNA data with an average depth of >0.25× for whole-genome sequencing or >1× for 1240k single nucleotide polymorphism capture data. Applying ancIBD to 4,248 ancient Eurasian individuals, we identify relatives up to the sixth degree and genealogical connections between archaeological groups. Notably, we reveal long IBD sharing between Corded Ware and Yamnaya groups, indicating that the Yamnaya herders of the Pontic-Caspian Steppe and the Steppe-related ancestry in various European Corded Ware groups share substantial co-ancestry within only a few hundred years. These results show that detecting IBD segments can generate powerful insights into the growing aDNA record, both on a small scale relevant to life stories and on a large scale relevant to major cultural-historical events.

A genetic history of the Balkans from Roman frontier to Slavic migrations

The rise and fall of the Roman Empire was a socio-political process with enormous ramifications for human history. The Middle Danube was a crucial frontier and a crossroads for population and cultural movement. Here, we present genome-wide data from 136 Balkan individuals dated to the 1st millennium CE. Despite extensive militarization and cultural influence, we find little ancestry contribution from peoples of Italic descent. However, we trace a large-scale influx of people of Anatolian ancestry during the Imperial period. Between ∼250 and 550 CE, we detect migrants with ancestry from Central/Northern Europe and the Steppe, confirming that "barbarian" migrations were propelled by ethnically diverse confederations. Following the end of Roman control, we detect the large-scale arrival of individuals who were genetically similar to modern Eastern European Slavic-speaking populations, who contributed 30%-60% of the ancestry of Balkan people, representing one of the largest permanent demographic changes anywhere in Europe during the Migration Period.

Genetic continuity and change among the Indigenous peoples of California

Before the colonial period, California harboured more language variation than all of Europe, and linguistic and archaeological analyses have led to many hypotheses to explain this diversity1. We report genome-wide data from 79 ancient individuals from California and 40 ancient individuals from Northern Mexico dating to 7,400-200 years before present (BP). Our analyses document long-term genetic continuity between people living on the Northern Channel Islands of California and the adjacent Santa Barbara mainland coast from 7,400 years BP to modern Chumash groups represented by individuals who lived around 200 years BP. The distinctive genetic lineages that characterize present-day and ancient people from Northwest Mexico increased in frequency in Southern and Central California by 5,200 years BP, providing evidence for northward migrations that are candidates for spreading Uto-Aztecan languages before the dispersal of maize agriculture from Mexico2-4. Individuals from Baja California share more alleles with the earliest individual from Central California in the dataset than with later individuals from Central California, potentially reflecting an earlier linguistic substrate, whose impact on local ancestry was diluted by later migrations from inland regions1,5. After 1,600 years BP, ancient individuals from the Channel Islands lived in communities with effective sizes similar to those in pre-agricultural Caribbean and Patagonia, and smaller than those on the California mainland and in sampled regions of Mexico.

Prospective Cohort Study to Evaluate Toxicity and Cosmesis between Interstitial APBI and Ultra-Hypo-Fractionated Whole Breast Irradiation in Patients with Early Breast Cancer

PURPOSE/OBJECTIVE(S)

Both interstitial multi-catheter brachytherapy based accelerated partial breast irradiation (APBI) and ultra-hypo-fractionated whole breast irradiation (WBI) with 5 fractions can complete treatment in a short time and result in lower treatment-related costs. Although data are available on their toxicity compared with conventional and hypo-fractionated regimens, there are few data directly comparing these two techniques.

MATERIALS/METHODS

This prospective cohort study enrolled early breast cancer patients with negative ink margin, T1 or T2 stage, and tumor size of 3 cm or less, with 0 to 3 positive lymph nodes after adequate axillary dissection or negative sentinel lymph node biopsy (N0-1), regardless of receptor status (ER, PR, HER2neu). Of the total 31 patients, 15 received APBI with multi-catheter interstitial brachytherapy (35 Gy in10 fractions over 5 days with a minimum interval of 6 hours between fractions) and 16 patients received WBI (27 Gy in 5 fractions over 5 days) with either 3DCRT or VMAT technique. Acute toxicity was assessed using the RTOG toxicity grading and breast cosmesis was assessed using the Harvard breast cosmesis score. Toxicity and cosmesis were assessed immediately after treatment, after one month, and after three months of treatment. Results were compared between groups using the Mann Whitney U test.

RESULTS

All patients completed a 3-month follow-up after completion of treatment. There was no statistically significant difference between the 2 groups in terms of highest grade toxicity (p = 0.97) and acute toxicity immediately after RT (p = 1.0) and after 1 month of RT (p = 0.91). However, at 3 months, more patients in the WBI group had a residual skin reaction than in the APBI group (p = 0.03). There was no statistically significant difference in breast cosmesis scores at any time point. 53.3%, 6.6% had grade 1 and 2 skin toxicities respectively immediately after RT in the APBI group and 62.5% of patients had grade 1 skin toxicity in WBI group. After 1 month 26.6%, 26.6% & 13.3% patients had grade 2, grade 1 and grade 3 toxicities respectively. In the WBI group, at 1 month, 50%, 18.7% and 12.5% patients had grade 1, 3 and 2 skin toxicities respectively. At 3 months, 42.8% and 7.1% of patients in the APBI group had grade 1 and grade 2 toxicities, respectively. In the WBI group, 78.5% and 7.1% of patients had grade 1 and 2 skin toxicities respectively.

CONCLUSION

In this study, no difference in acute toxicity or breast cosmesis was demonstrated between APBI-based interstitial brachytherapy and ultra-hypo-fractionated WBI. Since both regimens can be administered within 5 days, the decision between them is challenging for both physicians and patients. This study can be considered hypothesis-generating and used for further planning of equivalence or noninferiority studies.

Clinical Outcome and Patterns of Failure in Patients with Solitary Plasmacytoma Treated with Radical Radiotherapy

PURPOSE/OBJECTIVE(S)

We intended to assess the patterns of care and failure in patients with solitary plasmacytoma (SP), treated at our institute from 2011-21.

MATERIALS/METHODS

Data pertaining to the management of patients with SP was abstracted by retrospective chart review. Progression free survival (PFS) and overall survival (OS) were defined as the time intervals from the date of diagnosis to the date of progression (local, systemic or both) and death or last follow up, respectively. PFS and OS were analyzed by Kaplan Meier method.

RESULTS

Atotal of 72 patients were diagnosed with SP (61-bone plasmacytoma; 11-extramedullary plasmacytoma) in between 2011 and 2021. The median age at diagnosis was 52 years. The male: female ratio was 2.27. The ECOG PS were 0, 1, 2, 3, and 4 in 1, 46, 11, 10, and 3 patients, respectively. The median diameter of the lesion was 6.7 cm (range 2.3-18.4 cm). The median M protein was 0.375 g/dl (range 0.63-5.3 g/dl). On bone marrow examination, 8 and 64 patients had no and minimal marrow involvement, respectively (median bone marrow plasma cells 4%). Radical radiotherapy (RT) was administered in 71 patients, the techniques being 2D-conventional in 35, 3D-conformal in 14, IMRT in 8 and VMAT in 5 patients (details not available in 9). The median RT dose was 45 Gy in 25 fractions over 5 weeks (range 45-50 Gy). In addition, 5(7%) patients underwent surgery prior to RT and 6(8.5%) patients received systemic chemotherapy. After RT, the median quantum of symptom relief was 75% (range 0-100%). The best response to RT was CR, PR, SD and PD in 41(57.7%),18(25.4%),3(4.2%) and 4(5.6%) patients, respectively (response assessment could not be done in 5). Local control with primary RT was achieved in 58(81.7%) patients. No RT related grade 3/4 acute or late toxicity was reported. The median follow up was 48.3 months (range 1.3-141.6 months). A total of 28(39.4%) patients experienced disease progression (systemic in 15, local in 7 and both in 6), out of which 13 could be successfully salvaged with further systemic chemotherapy (mostly bortezomib, lenalidomide/pomalidomide and dexamethasone based) with or without local treatment. On progression to multiple myeloma, 9 patients underwent high dose chemotherapy(melphalan) followed by autologous stem cell transplant. Nine (12.7%) patients died, the causes being COVID-19 infection in 1, surgical complications in 1, PD in 2, second malignancies in 3 and unknown in 2. The actuarial rates of PFS & OS were 58.4% & 92.5% at 3 years and 50.3% & 89.4% at 5 years, respectively. Metachronous second malignancies were noted in 4 patients (B-ALL-1, lung cancer-1, metastatic carcinoma cervix-1, metastatic carcinoma gall bladder-1).

CONCLUSION

In patients with SP, radical RT with a modest dose (45-50 Gy) leads to satisfactory symptom relief, response and local control. Close to 40% patients eventually progress to multiple myeloma, out of which 50% can be successfully salvaged with further systemic chemotherapy. The 5-year OS approaches 90% in this cohort.

Optimum Radiation Dose for Palliation in Head and Neck SCC: A Phase III Randomized Controlled Trial (OPRAH)

PURPOSE/OBJECTIVE(S)

Palliative radiotherapy has been used to alleviate symptoms in locally advanced (T4b) head and neck Squamous cell carcinoma, but optimum dose is still debated. The trial has been designed to determine if high palliative radiation dose (30 Gy/ 5#/ 5 days) improves progression free survival compared to low palliative dose (20 Gy/ 5# 5 days). The secondary objective was to compare toxicity rates profile between the two radiation doses (CTCAE version 5) and compare Quality of life between the two arms with QOL - HN35 questionnaire.

MATERIALS/METHODS

This is a prospective randomized phase III randomized controlled trial. We hypothesized that by increasing the dose per fraction and total dose there would be an improvement in PFS without difference in late toxicity and QOL. The present study will also perform an exploratory analysis of effect of palliative radiotherapy on immune function and its effect on survival. Considering median PFS 2.5 months in 20 Gy arm we assumed that PFS will improve by 1.25 months (total PFS 3.75 months) in 30 Gy arm. To reach the hypothesis with 90% power and alpha error of 0.05, attrition rate of 10% sample size calculated was 286. Patients with cT4b disease of larynx/ Hypopharynx/ Oropharynx/ Oral cavity not amenable to surgery, Carotid encasement > 180 degree, histologically proven squamous cell carcinoma, Age > 18 years and < 75 years ECOG PS 0-2 were eligible.

RESULTS

Until December 2022, we completed accrual of 239 patients. Median age: 53 years (Range: 29-70 years); Gender: Male: Female = 9: 1 (214: 25). 145 patients had history of smoking, 87 had additional history of tobacco chewing. 67 patients had history of alcohol consumption, 13 patients had history of smoking, tobacco and alcoholism. All patients had stage IVA or Stage IVB disease. Primary site was oral cavity in 58%, Oropharynx in 19%, larynx in 10% and hypopharynx in 8% cases. 72.3% patients had a PS of 1, 21.7% had PS 2 at presentation. Prior to RT requirement of analgesic was: No pain medicine: 5; WHO step I:173, WHO step II: 52, WHO step III:14. All patients received conformal radiation with 6 MV linear accelerator. Target volume was GTV plus 1 cm margin to form CTV and 5 mm expansion to generate PTV. Out of 239 patients 67 patients had a disease progression. Mucositis is the commonest side effect affecting most of the patients and improve with best supportive care.

CONCLUSION

High dose palliation appears to be equally well tolerated. However, the impact of such dose escalation will be established if final analysis reaches the pre-defined endpoint.

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 since it first was made available and crossed the threshold of >10,000 ancient individuals with genome-wide data at the end of 2022. This note is intended as a citable description of the AADR.

ancIBD - Screening for identity by descent segments in human ancient DNA

Long DNA sequences shared between two individuals, known as Identical by descent (IBD) segments, are a powerful signal for identifying close and distant biological relatives because they only arise when the pair shares a recent common ancestor. Existing methods to call IBD segments between present-day genomes cannot be straightforwardly applied to ancient DNA data (aDNA) due to typically low coverage and high genotyping error rates. We present ancIBD, a method to identify IBD segments for human aDNA data implemented as a Python package. Our approach is based on a Hidden Markov Model, using as input genotype probabilities imputed based on a modern reference panel of genomic variation. Through simulation and downsampling experiments, we demonstrate that ancIBD robustly identifies IBD segments longer than 8 centimorgan for aDNA data with at least either 0.25x average whole-genome sequencing (WGS) coverage depth or at least 1x average depth for in-solution enrichment experiments targeting a widely used aDNA SNP set ('1240k'). This application range allows us to screen a substantial fraction of the aDNA record for IBD segments and we showcase two downstream applications. First, leveraging the fact that biological relatives up to the sixth degree are expected to share multiple long IBD segments, we identify relatives between 10,156 ancient Eurasian individuals and document evidence of long-distance migration, for example by identifying a pair of two approximately fifth-degree relatives who were buried 1410km apart in Central Asia 5000 years ago. Second, by applying ancIBD, we reveal new details regarding the spread of ancestry related to Steppe pastoralists into Europe starting 5000 years ago. We find that the first individuals in Central and Northern Europe carrying high amounts of Steppe-ancestry, associated with the Corded Ware culture, share high rates of long IBD (12-25 cM) with Yamnaya herders of the Pontic-Caspian steppe, signaling a strong bottleneck and a recent biological connection on the order of only few hundred years, providing evidence that the Yamnaya themselves are a main source of Steppe ancestry in Corded Ware people. We also detect elevated sharing of long IBD segments between Corded Ware individuals and people associated with the Globular Amphora culture (GAC) from Poland and Ukraine, who were Copper Age farmers not yet carrying Steppe-like ancestry. These IBD links appear for all Corded Ware groups in our analysis, indicating that individuals related to GAC contexts must have had a major demographic impact early on in the genetic admixtures giving rise to various Corded Ware groups across Europe. These results show that detecting IBD segments in aDNA can generate new insights both on a small scale, relevant to understanding the life stories of people, and on the macroscale, relevant to large-scale cultural-historical events.

Entwined African and Asian genetic roots of medieval peoples of the Swahili coast

The urban peoples of the Swahili coast traded across eastern Africa and the Indian Ocean and were among the first practitioners of Islam among sub-Saharan people1,2. The extent to which these early interactions between Africans and non-Africans were accompanied by genetic exchange remains unknown. Here we report ancient DNA data for 80 individuals from 6 medieval and early modern (AD 1250-1800) coastal towns and an inland town after AD 1650. More than half of the DNA of many of the individuals from coastal towns originates from primarily female ancestors from Africa, with a large proportion-and occasionally more than half-of the DNA coming from Asian ancestors. The Asian ancestry includes components associated with Persia and India, with 80-90% of the Asian DNA originating from Persian men. Peoples of African and Asian origins began to mix by about AD 1000, coinciding with the large-scale adoption of Islam. Before about AD 1500, the Southwest Asian ancestry was mainly Persian-related, consistent with the narrative of the Kilwa Chronicle, the oldest history told by people of the Swahili coast3. After this time, the sources of DNA became increasingly Arabian, consistent with evidence of growing interactions with southern Arabia4. Subsequent interactions with Asian and African people further changed the ancestry of present-day people of the Swahili coast in relation to the medieval individuals whose DNA we sequenced.

Ancient DNA from Mesopotamia suggests distinct Pre-Pottery and Pottery Neolithic migrations into Anatolia

We present the first ancient DNA data from the Pre-Pottery Neolithic of Mesopotamia (Southeastern Turkey and Northern Iraq), Cyprus, and the Northwestern Zagros, along with the first data from Neolithic Armenia. We show that these and neighboring populations were formed through admixture of pre-Neolithic sources related to Anatolian, Caucasus, and Levantine hunter-gatherers, forming a Neolithic continuum of ancestry mirroring the geography of West Asia. By analyzing Pre-Pottery and Pottery Neolithic populations of Anatolia, we show that the former were derived from admixture between Mesopotamian-related and local Epipaleolithic-related sources, but the latter experienced additional Levantine-related gene flow, thus documenting at least two pulses of migration from the Fertile Crescent heartland to the early farmers of Anatolia.

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 genetic probe into the ancient and medieval history of Southern Europe and West Asia

Literary and archaeological sources have preserved a rich history of Southern Europe and West Asia since the Bronze Age that can be complemented by genetics. Mycenaean period elites in Greece did not differ from the general population and included both people with some steppe ancestry and others, like the Griffin Warrior, without it. Similarly, people in the central area of the Urartian Kingdom around Lake Van lacked the steppe ancestry characteristic of the kingdom's northern provinces. Anatolia exhibited extraordinary continuity down to the Roman and Byzantine periods, with its people serving as the demographic core of much of the Roman Empire, including the city of Rome itself. During medieval times, migrations associated with Slavic and Turkic speakers profoundly affected the region.

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.

Three assays for in-solution enrichment of ancient human DNA at more than a million SNPs

The strategy of in-solution enrichment for hundreds of thousands of single-nucleotide polymorphisms (SNPs) has been used to analyze >70% of individuals with genome-scale ancient DNA published to date. This approach makes it economical to study ancient samples with low proportions of human DNA and increases the rate of conversion of sampled remains into interpretable data. So far, nearly all such data have been generated using a set of bait sequences targeting about 1.24 million SNPs (the "1240k reagent"), but synthesis of the reagent has been cost-effective for only a few laboratories. In 2021, two companies, Daicel Arbor Biosciences and Twist Bioscience, made available assays that target the same core set of SNPs along with supplementary content. We test all three assays on a common set of 27 ancient DNA libraries and show that all three are effective at enriching many hundreds of thousands of SNPs. For all assays, one round of enrichment produces data that are as useful as two. In our testing, the "Twist Ancient DNA" assay produces the highest coverages, greatest uniformity on targeted positions, and almost no bias toward enriching one allele more than another relative to shotgun sequencing. We also identify hundreds of thousands of targeted SNPs for which there is minimal allelic bias when comparing 1240k data to either shotgun or Twist data. This facilitates coanalysis of the large data sets that have been generated using 1240k and Twist capture, as well as shotgun sequencing approaches.

Social stratification without genetic differentiation at the site of Kulubnarti in Christian Period Nubia

Relatively little is known about Nubia's genetic landscape prior to the influence of the Islamic migrations that began in the late 1st millennium CE. Here, we increase the number of ancient individuals with genome-level data from the Nile Valley from three to 69, reporting data for 66 individuals from two cemeteries at the Christian Period (~650-1000 CE) site of Kulubnarti, where multiple lines of evidence suggest social stratification. The Kulubnarti Nubians had ~43% Nilotic-related ancestry (individual variation between ~36-54%) with the remaining ancestry consistent with being  introduced through Egypt and ultimately deriving from an ancestry pool like that found in the Bronze and Iron Age Levant. The Kulubnarti gene pool - shaped over a millennium - harbors disproportionately female-associated West Eurasian-related ancestry. Genetic similarity among individuals from the two cemeteries supports a hypothesis of social division without genetic distinction. Seven pairs of inter-cemetery relatives suggest fluidity between cemetery groups. Present-day Nubians are not directly descended from the Kulubnarti Nubians, attesting to additional genetic input since the Christian Period.

Ancient genomes reveal long-range influence of the pre-Columbian culture and site of Tiwanaku

Tiwanaku civilization flourished in the Lake Titicaca basin between 500 and 1000 CE and at its apogee influenced wide areas across the southern Andes. Despite a considerable amount of archaeological data, little is known about the Tiwanaku population. We analyzed 17 low-coverage genomes from individuals dated between 300 and 1500 CE and demonstrated genetic continuity in the Lake Titicaca basin throughout this period, which indicates that the substantial cultural and political changes in the region were not accompanied by large-scale population movements. Conversely, the ritual center of Tiwanaku revealed high diversity, including individuals with primarily local genetic ancestry and those with foreign admixture or provenance from as far as the Amazon. Nonetheless, most human offerings associated with the Akapana platform exhibited pure Titicaca basin ancestry and dated to ca. 950 CE—the onset of Tiwanaku’s decline as a sociopolitical center. Our results strengthen the view of Tiwanaku as a complex and far-reaching polity.

Genomic insights into the formation of human populations in East Asia

The deep population history of East Asia remains poorly understood owing to a lack of ancient DNA data and sparse sampling of present-day people1,2. Here we report genome-wide data from 166 East Asian individuals dating to between 6000 BC and AD 1000 and 46 present-day groups. Hunter-gatherers from Japan, the Amur River Basin, and people of Neolithic and Iron Age Taiwan and the Tibetan Plateau are linked by a deeply splitting lineage that probably reflects a coastal migration during the Late Pleistocene epoch. We also follow expansions during the subsequent Holocene epoch from four regions. First, hunter-gatherers from Mongolia and the Amur River Basin have ancestry shared by individuals who speak Mongolic and Tungusic languages, but do not carry ancestry characteristic of farmers from the West Liao River region (around 3000 BC), which contradicts theories that the expansion of these farmers spread the Mongolic and Tungusic proto-languages. Second, farmers from the Yellow River Basin (around 3000 BC) probably spread Sino-Tibetan languages, as their ancestry dispersed both to Tibet-where it forms approximately 84% of the gene pool in some groups-and to the Central Plain, where it has contributed around 59-84% to modern Han Chinese groups. Third, people from Taiwan from around 1300 BC to AD 800 derived approximately 75% of their ancestry from a lineage that is widespread in modern individuals who speak Austronesian, Tai-Kadai and Austroasiatic languages, and that we hypothesize derives from farmers of the Yangtze River Valley. Ancient people from Taiwan also derived about 25% of their ancestry from a northern lineage that is related to, but different from, farmers of the Yellow River Basin, which suggests an additional north-to-south expansion. Fourth, ancestry from Yamnaya Steppe pastoralists arrived in western Mongolia after around 3000 BC but was displaced by previously established lineages even while it persisted in western China, as would be expected if this ancestry was associated with the spread of proto-Tocharian Indo-European languages. Two later gene flows affected western Mongolia: migrants after around 2000 BC with Yamnaya and European farmer ancestry, and episodic influences of later groups with ancestry from Turan.

A genetic history of the pre-contact Caribbean

Humans settled the Caribbean about 6,000 years ago, and ceramic use and intensified agriculture mark a shift from the Archaic to the Ceramic Age at around 2,500 years ago1-3. Here we report genome-wide data from 174 ancient individuals from The Bahamas, Haiti and the Dominican Republic (collectively, Hispaniola), Puerto Rico, Curaçao and Venezuela, which we co-analysed with 89 previously published ancient individuals. Stone-tool-using Caribbean people, who first entered the Caribbean during the Archaic Age, derive from a deeply divergent population that is closest to Central and northern South American individuals; contrary to previous work4, we find no support for ancestry contributed by a population related to North American individuals. Archaic-related lineages were >98% replaced by a genetically homogeneous ceramic-using population related to speakers of languages in the Arawak family from northeast South America; these people moved through the Lesser Antilles and into the Greater Antilles at least 1,700 years ago, introducing ancestry that is still present. Ancient Caribbean people avoided close kin unions despite limited mate pools that reflect small effective population sizes, which we estimate to be a minimum of 500-1,500 and a maximum of 1,530-8,150 individuals on the combined islands of Puerto Rico and Hispaniola in the dozens of generations before the individuals who we analysed lived. Census sizes are unlikely to be more than tenfold larger than effective population sizes, so previous pan-Caribbean estimates of hundreds of thousands of people are too large5,6. Confirming a small and interconnected Ceramic Age population7, we detect 19 pairs of cross-island cousins, close relatives buried around 75 km apart in Hispaniola and low genetic differentiation across islands. Genetic continuity across transitions in pottery styles reveals that cultural changes during the Ceramic Age were not driven by migration of genetically differentiated groups from the mainland, but instead reflected interactions within an interconnected Caribbean world1,8.

Mitochondrial genome diversity on the Central Siberian Plateau with particular reference to the prehistory of northernmost Eurasia

The Central Siberian Plateau was the last geographic area in Eurasia to become habitable by modern humans after the Last Glacial Maximum (LGM). Through a comprehensive dataset of mitochondrial DNA (mtDNA) genomes retained in the remnats of earlier ("Old") Siberians, primarily the Ket, Tofalar, and Todzhi, we explored genetic links between the Yenisei-Sayan region and Northeast Eurasia (best represented by the Yukaghir) over the last 10,000 years. We generated 218 new complete mtDNA sequences and placed them into compound phylogenies with 7 newly obtained and 70 published ancient mitochondrial genomes. We have considerably extended the mtDNA sequence diversity (at the entire mtDNA genome level) of autochthonous Siberians, which remain poorly sampled, and these new data may have a broad impact on the study of human migration. We compared present-day mtDNA diversity in these groups with complete mitochondrial genomes from ancient samples from the region and placed the samples into combined genealogical trees. The resulting components were used to clarify the origins and expansion history of mtDNA lineages that evolved in the refugia of south-central Siberia and beyond, as well as multiple phases of connection between this region and distant parts of Eurasia.

ContamLD: estimation of ancient nuclear DNA contamination using breakdown of linkage disequilibrium

We report a method called ContamLD for estimating autosomal ancient DNA (aDNA) contamination by measuring the breakdown of linkage disequilibrium in a sequenced individual due to the introduction of contaminant DNA. ContamLD leverages the idea that contaminants should have haplotypes uncorrelated to those of the studied individual. Using simulated data, we confirm that ContamLD accurately infers contamination rates with low standard errors: for example, less than 1.5% standard error in cases with less than 10% contamination and 500,000 sequences covering SNPs. This method is optimized for application to aDNA, taking advantage of characteristic aDNA damage patterns to provide calibrated contamination estimates, and is available at https://github.com/nathan-nakatsuka/ContamLD .

Ancient genomes in South Patagonia reveal population movements associated with technological shifts and geography

Archaeological research documents major technological shifts among people who have lived in the southern tip of South America (South Patagonia) during the last thirteen millennia, including the development of marine-based economies and changes in tools and raw materials. It has been proposed that movements of people spreading culture and technology propelled some of these shifts, but these hypotheses have not been tested with ancient DNA. Here we report genome-wide data from 20 ancient individuals, and co-analyze it with previously reported data. We reveal that immigration does not explain the appearance of marine adaptations in South Patagonia. We describe partial genetic continuity since ~6600 BP and two later gene flows correlated with technological changes: one between 4700-2000 BP that affected primarily marine-based groups, and a later one impacting all <2000 BP groups. From ~2200-1200 BP, mixture among neighbors resulted in a cline correlated to geographic ordering along the coast.

A Paleogenomic Reconstruction of the Deep Population History of the Andes

There are many unanswered questions about the population history of the Central and South Central Andes, particularly regarding the impact of large-scale societies, such as the Moche, Wari, Tiwanaku, and Inca. We assembled genome-wide data on 89 individuals dating from ∼9,000-500 years ago (BP), with a particular focus on the period of the rise and fall of state societies. Today's genetic structure began to develop by 5,800 BP, followed by bi-directional gene flow between the North and South Highlands, and between the Highlands and Coast. We detect minimal admixture among neighboring groups between ∼2,000-500 BP, although we do detect cosmopolitanism (people of diverse ancestries living side-by-side) in the heartlands of the Tiwanaku and Inca polities. We also highlight cases of long-range mobility connecting the Andes to Argentina and the Northwest Andes to the Amazon Basin. VIDEO ABSTRACT.

Author Correction: The spread of steppe and Iranian-related ancestry in the islands of the western Mediterranean

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Insights into human genetic variation and population history from 929 diverse genomes

Genome sequences from diverse human groups are needed to understand the structure of genetic variation in our species and the history of, and relationships between, different populations. We present 929 high-coverage genome sequences from 54 diverse human populations, 26 of which are physically phased using linked-read sequencing. Analyses of these genomes reveal an excess of previously undocumented common genetic variation private to southern Africa, central Africa, Oceania, and the Americas, but an absence of such variants fixed between major geographical regions. We also find deep and gradual population separations within Africa, contrasting population size histories between hunter-gatherer and agriculturalist groups in the past 10,000 years, and a contrast between single Neanderthal but multiple Denisovan source populations contributing to present-day human populations.

Differential DNA methylation of vocal and facial anatomy genes in modern humans

Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated genes (SOX9, ACAN, COL2A1, NFIX and XYLT1). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Ancient West African foragers in the context of African population history

Our knowledge of ancient human population structure in sub-Saharan Africa, particularly prior to the advent of food production, remains limited. Here we report genome-wide DNA data from four children-two of whom were buried approximately 8,000 years ago and two 3,000 years ago-from Shum Laka (Cameroon), one of the earliest known archaeological sites within the probable homeland of the Bantu language group1-11. One individual carried the deeply divergent Y chromosome haplogroup A00, which today is found almost exclusively in the same region12,13. However, the genome-wide ancestry profiles of all four individuals are most similar to those of present-day hunter-gatherers from western Central Africa, which implies that populations in western Cameroon today-as well as speakers of Bantu languages from across the continent-are not descended substantially from the population represented by these four people. We infer an Africa-wide phylogeny that features widespread admixture and three prominent radiations, including one that gave rise to at least four major lineages deep in the history of modern humans.

The formation of human populations in South and Central Asia

By sequencing 523 ancient humans, we show that the primary source of ancestry in modern South Asians is a prehistoric genetic gradient between people related to early hunter-gatherers of Iran and Southeast Asia. After the Indus Valley Civilization's decline, its people mixed with individuals in the southeast to form one of the two main ancestral populations of South Asia, whose direct descendants live in southern India. Simultaneously, they mixed with descendants of Steppe pastoralists who, starting around 4000 years ago, spread via Central Asia to form the other main ancestral population. The Steppe ancestry in South Asia has the same profile as that in Bronze Age Eastern Europe, tracking a movement of people that affected both regions and that likely spread the distinctive features shared between Indo-Iranian and Balto-Slavic languages.

An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers

We report an ancient genome from the Indus Valley Civilization (IVC). The individual we sequenced fits as a mixture of people related to ancient Iranians (the largest component) and Southeast Asian hunter-gatherers, a unique profile that matches ancient DNA from 11 genetic outliers from sites in Iran and Turkmenistan in cultural communication with the IVC. These individuals had little if any Steppe pastoralist-derived ancestry, showing that it was not ubiquitous in northwest South Asia during the IVC as it is today. The Iranian-related ancestry in the IVC derives from a lineage leading to early Iranian farmers, herders, and hunter-gatherers before their ancestors separated, contradicting the hypothesis that the shared ancestry between early Iranians and South Asians reflects a large-scale spread of western Iranian farmers east. Instead, sampled ancient genomes from the Iranian plateau and IVC descend from different groups of hunter-gatherers who began farming without being connected by substantial movement of people.

Ancient DNA from the skeletons of Roopkund Lake reveals Mediterranean migrants in India

Situated at over 5,000 meters above sea level in the Himalayan Mountains, Roopkund Lake is home to the scattered skeletal remains of several hundred individuals of unknown origin. We report genome-wide ancient DNA for 38 skeletons from Roopkund Lake, and find that they cluster into three distinct groups. A group of 23 individuals have ancestry that falls within the range of variation of present-day South Asians. A further 14 have ancestry typical of the eastern Mediterranean. We also identify one individual with Southeast Asian-related ancestry. Radiocarbon dating indicates that these remains were not deposited simultaneously. Instead, all of the individuals with South Asian-related ancestry date to ~800 CE (but with evidence of being deposited in more than one event), while all other individuals date to ~1800 CE. These differences are also reflected in stable isotope measurements, which reveal a distinct dietary profile for the two main groups.

Remains of several hundred humans are scattered around Roopkund Lake, situated over 5,000 meters above sea level in the Himalayan Mountains. Here the authors analyze genome-wide data from 38 skeletons and find 3 clusters with different ancestries and dates, showing the people were desposited in multiple catastrophic events.

The genomic history of the Iberian Peninsula over the past 8000 years

We assembled genome-wide data from 271 ancient Iberians, of whom 176 are from the largely unsampled period after 2000 BCE, thereby providing a high-resolution time transect of the Iberian Peninsula. We document high genetic substructure between northwestern and southeastern hunter-gatherers before the spread of farming. We reveal sporadic contacts between Iberia and North Africa by ~2500 BCE and, by ~2000 BCE, the replacement of 40% of Iberia's ancestry and nearly 100% of its Y-chromosomes by people with Steppe ancestry. We show that, in the Iron Age, Steppe ancestry had spread not only into Indo-European-speaking regions but also into non-Indo-European-speaking ones, and we reveal that present-day Basques are best described as a typical Iron Age population without the admixture events that later affected the rest of Iberia. Additionally, we document how, beginning at least in the Roman period, the ancestry of the peninsula was transformed by gene flow from North Africa and the eastern Mediterranean.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.

Author Correction: Ancient genomes indicate population replacement in Early Neolithic Britain

In the version of this Article originally published, there were errors in the colour ordering of the legend in Fig. 5b, and in the positions of the target and surrogate populations in Fig. 5c. This has now been corrected. The conclusions of the study are in no way affected. The errors have been corrected in the HTML and PDF versions of the article.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations

BACKGROUND

Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations.

RESULTS

Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations.

CONCLUSIONS

We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.