Erratum: The Beaker phenomenon and the genomic transformation of northwest Europe

This corrects the article DOI: 10.1038/nature25738.

Reconstructing the genetic history of late Neanderthals

Although it is known that Neandertals contributed DNA to modern humans1,2, not much is known about the genetic diversity of Neandertals or the relationship between late Neandertal populations at the time when their last interactions with early modern humans occurred and before they eventually disappeared. Our ability to retrieve DNA from a larger number of Neandertal individuals has been limited by poor preservation of endogenous DNA3 and large amounts of microbial and present-day human DNA that contaminate Neandertal skeletal remains3–5. Here we use hypochlorite treatment6 of as little as 9 mg of bone or tooth powder to generate between 1- and 2.7-fold genomic coverage of five 39,000- to 47,000-year-old Neandertals (i.e. late Neandertals), thereby doubling the number of Neandertals for which genome sequences are available. Genetic similarity among late Neandertals is well predicted by their geographical location, and comparison to the genome of an older Neandertal from the Caucasus2,7 indicates that a population turnover is likely to have occurred, either in the Caucasus or throughout Europe, towards the end of Neandertal history. We find that the bulk of Neandertal gene flow into early modern humans originated from one or more source populations that diverged from the Neandertals studied here at least 70,000 years ago, but after they split from a previously sequenced Neandertal from Siberia2 ~150,000 years ago. Although four of these Neandertals post-date the putative arrival of early modern humans into Europe, we do not detect any recent gene flow from early modern humans in their ancestry.

A comprehensive genomic history of extinct and living elephants

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

The Beaker phenomenon and the genomic transformation of northwest Europe

From around 2750 to 2500 bc, Bell Beaker pottery became widespread across western and central Europe, before it disappeared between 2200 and 1800 bc. The forces that propelled its expansion are a matter of long-standing debate, and there is support for both cultural diffusion and migration having a role in this process. Here we present genome-wide data from 400 Neolithic, Copper Age and Bronze Age Europeans, including 226 individuals associated with Beaker-complex artefacts. We detected limited genetic affinity between Beaker-complex-associated individuals from Iberia and central Europe, and thus exclude migration as an important mechanism of spread between these two regions. However, migration had a key role in the further dissemination of the Beaker complex. We document this phenomenon most clearly in Britain, where the spread of the Beaker complex introduced high levels of steppe-related ancestry and was associated with the replacement of approximately 90% of Britain's gene pool within a few hundred years, continuing the east-to-west expansion that had brought steppe-related ancestry into central and northern Europe over the previous centuries.

Extracorporeal cellular therapy (ELAD) in severe alcoholic hepatitis: A multinational, prospective, controlled, randomized trial

Severe alcoholic hepatitis (sAH) is associated with a poor prognosis. There is no proven effective treatment for sAH, which is why early transplantation has been increasingly discussed. Hepatoblastoma-derived C3A cells express anti-inflammatory proteins and growth factors and were tested in an extracorporeal cellular therapy (ELAD) study to establish their effect on survival for subjects with sAH. Adults with sAH, bilirubin ≥8 mg/dL, Maddrey's discriminant function ≥ 32, and Model for End-Stage Liver Disease (MELD) score ≤ 35 were randomized to receive standard of care (SOC) only or 3-5 days of continuous ELAD treatment plus SOC. After a minimum follow-up of 91 days, overall survival (OS) was assessed by using a Kaplan-Meier survival analysis. A total of 203 subjects were enrolled (96 ELAD and 107 SOC) at 40 sites worldwide. Comparison of baseline characteristics showed no significant differences between groups and within subgroups. There was no significant difference in serious adverse events between the 2 groups. In an analysis of the intent-to-treat population, there was no difference in OS (51.0% versus 49.5%). The study failed its primary and secondary end point in a population with sAH and with a MELD ranging from 18 to 35 and no upper age limit. In the prespecified analysis of subjects with MELD < 28 (n = 120), ELAD was associated with a trend toward higher OS at 91 days (68.6% versus 53.6%; P = .08). Regression analysis identified high creatinine and international normalized ratio, but not bilirubin, as the MELD components predicting negative outcomes with ELAD. A new trial investigating a potential benefit of ELAD in younger subjects with sufficient renal function and less severe coagulopathy has been initiated. Liver Transplantation 24 380-393 2018 AASLD.

Biological Sexing of a 4000-Year-Old Egyptian Mummy Head to Assess the Potential of Nuclear DNA Recovery from the Most Damaged and Limited Forensic Specimens

High throughput sequencing (HTS) has been used for a number of years in the field of paleogenomics to facilitate the recovery of small DNA fragments from ancient specimens. Recently, these techniques have also been applied in forensics, where they have been used for the recovery of mitochondrial DNA sequences from samples where traditional PCR-based assays fail because of the very short length of endogenous DNA molecules. Here, we describe the biological sexing of a ~4000-year-old Egyptian mummy using shotgun sequencing and two established methods of biological sex determination (R_X and R_Y), by way of mitochondrial genome analysis as a means of sequence data authentication. This particular case of historical interest increases the potential utility of HTS techniques for forensic purposes by demonstrating that data from the more discriminatory nuclear genome can be recovered from the most damaged specimens, even in cases where mitochondrial DNA cannot be recovered with current PCR-based forensic technologies. Although additional work remains to be done before nuclear DNA recovered via these methods can be used routinely in operational casework for individual identification purposes, these results indicate substantial promise for the retrieval of probative individually identifying DNA data from the most limited and degraded forensic specimens.

Population Turnover in Remote Oceania Shortly after Initial Settlement

Ancient DNA from Vanuatu and Tonga dating to about 2,900-2,600 years ago (before present, BP) has revealed that the "First Remote Oceanians" associated with the Lapita archaeological culture were directly descended from the population that, beginning around 5000 BP, spread Austronesian languages from Taiwan to the Philippines, western Melanesia, and eventually Remote Oceania. Thus, ancestors of the First Remote Oceanians must have passed by the Papuan-ancestry populations they encountered in New Guinea, the Bismarck Archipelago, and the Solomon Islands with minimal admixture [1]. However, all present-day populations in Near and Remote Oceania harbor >25% Papuan ancestry, implying that additional eastward migration must have occurred. We generated genome-wide data for 14 ancient individuals from Efate and Epi Islands in Vanuatu from 2900-150 BP, as well as 185 present-day individuals from 18 islands. We find that people of almost entirely Papuan ancestry arrived in Vanuatu by around 2300 BP, most likely reflecting migrations a few hundred years earlier at the end of the Lapita period, when there is also evidence of changes in skeletal morphology and cessation of long-distance trade between Near and Remote Oceania [2, 3]. Papuan ancestry was subsequently diluted through admixture but remains at least 80%-90% in most islands. Through a fine-grained analysis of ancestry profiles, we show that the Papuan ancestry in Vanuatu derives from the Bismarck Archipelago rather than the geographically closer Solomon Islands. However, the Papuan ancestry in Polynesia-the most remote Pacific islands-derives from different sources, documenting a third stream of migration from Near to Remote Oceania.

The genomic history of southeastern Europe

Farming was first introduced to Europe in the mid-7^th millennium BCE–associated with migrants from Anatolia who settled in the Southeast before spreading throughout Europe. To understand the dynamics of this process, we analyzed genome-wide ancient DNA data from 225 individuals who lived in southeastern Europe and surrounding regions between 12,000 and 500 BCE. We document a West-East cline of ancestry in indigenous hunter-gatherers and–in far-eastern Europe–early stages in the formation of Bronze Age Steppe ancestry. We show that the first farmers of northern and western Europe passed through southeastern Europe with limited hunter-gatherer admixture, but that some groups that remained mixed extensively, without the male-biased hunter-gatherer admixture that prevailed later in the North and West. Southeastern Europe continued to be a nexus between East and West, with intermittent genetic contact with the Steppe up to 2000 years before the migrations that replaced much of northern Europe’s population.

Meiotic Genes Are Enriched in Regions of Reduced Archaic Ancestry

About 1-6% of the genetic ancestry of modern humans today originates from admixture with archaic humans. It has recently been shown that autosomal genomic regions with a reduced proportion of Neanderthal and Denisovan ancestries (NA and DA) are significantly enriched in genes that are more expressed in testis than in other tissues. To determine whether a cellular segregation pattern would exist, we combined maps of archaic introgression with a cross-analysis of three transcriptomic datasets deciphering the transcriptional landscape of human gonadal cell types. We reveal that the regions deficient in both NA and DA contain a significant enrichment of genes transcribed in meiotic germ cells. The interbreeding of anatomically modern humans with archaic humans may have introduced archaic-derived alleles that contributed to genetic incompatibilities affecting meiosis that were subsequently purged by natural selection.

The genetic prehistory of the Baltic Sea region

While the series of events that shaped the transition between foraging societies and food producers are well described for Central and Southern Europe, genetic evidence from Northern Europe surrounding the Baltic Sea is still sparse. Here, we report genome-wide DNA data from 38 ancient North Europeans ranging from ~9500 to 2200 years before present. Our analysis provides genetic evidence that hunter-gatherers settled Scandinavia via two routes. We reveal that the first Scandinavian farmers derive their ancestry from Anatolia 1000 years earlier than previously demonstrated. The range of Mesolithic Western hunter-gatherers extended to the east of the Baltic Sea, where these populations persisted without gene-flow from Central European farmers during the Early and Middle Neolithic. The arrival of steppe pastoralists in the Late Neolithic introduced a major shift in economy and mediated the spread of a new ancestry associated with the Corded Ware Complex in Northern Europe.

Genome diversity in the Neolithic Globular Amphorae culture and the spread of Indo-European languages

It is unclear whether Indo-European languages in Europe spread from the Pontic steppes in the late Neolithic, or from Anatolia in the Early Neolithic. Under the former hypothesis, people of the Globular Amphorae culture (GAC) would be descended from Eastern ancestors, likely representing the Yamnaya culture. However, nuclear (six individuals typed for 597 573 SNPs) and mitochondrial (11 complete sequences) DNA from the GAC appear closer to those of earlier Neolithic groups than to the DNA of all other populations related to the Pontic steppe migration. Explicit comparisons of alternative demographic models via approximate Bayesian computation confirmed this pattern. These results are not in contrast to Late Neolithic gene flow from the Pontic steppes into Central Europe. However, they add nuance to this model, showing that the eastern affinities of the GAC in the archaeological record reflect cultural influences from other groups from the East, rather than the movement of people.

A high-coverage Neandertal genome from Vindija Cave in Croatia

To date, the only Neandertal genome that has been sequenced to high quality is from an individual found in Southern Siberia. We sequenced the genome of a female Neandertal from ~50,000 years ago from Vindija Cave, Croatia, to ~30-fold genomic coverage. She carried 1.6 differences per 10,000 base pairs between the two copies of her genome, fewer than present-day humans, suggesting that Neandertal populations were of small size. Our analyses indicate that she was more closely related to the Neandertals that mixed with the ancestors of present-day humans living outside of sub-Saharan Africa than the previously sequenced Neandertal from Siberia, allowing 10 to 20% more Neandertal DNA to be identified in present-day humans, including variants involved in low-density lipoprotein cholesterol concentrations, schizophrenia, and other diseases.

Erratum to: A multi-stage genome-wide association study of uterine fibroids in African Americans

The article "A multi-stage genome-wide association study of uterine fibroids in African Americans", written by Jacklyn N. Hellwege, was originally published Online First without open access. After publication in volume 136, issue 10, page 1363-1373 the author decided to opt for Open Choice and to make the article an open access publication. Therefore, the copyright of the article has been changed to

Interpreting short tandem repeat variations in humans using mutational constraint

Identifying regions of the genome that are depleted of mutations can distinguish potentially deleterious variants. Short tandem repeats (STRs), also known as microsatellites, are among the largest contributors of de novo mutations in humans. However, per-locus studies of STR mutations have been limited to highly ascertained panels of several dozen loci. Here we harnessed bioinformatics tools and a novel analytical framework to estimate mutation parameters for each STR in the human genome by correlating STR genotypes with local sequence heterozygosity. We applied our method to obtain robust estimates of the impact of local sequence features on mutation parameters and used these estimates to create a framework for measuring constraint at STRs by comparing observed versus expected mutation rates. Constraint scores identified known pathogenic variants with early-onset effects. Our metric will provide a valuable tool for prioritizing pathogenic STRs in medical genetics studies.

A multi-stage genome-wide association study of uterine fibroids in African Americans

Uterine fibroids are benign tumors of the uterus affecting up to 77% of women by menopause. They are the leading indication for hysterectomy, and account for $34 billion annually in the United States. Race/ethnicity and age are the strongest known risk factors. African American (AA) women have higher prevalence, earlier onset, and larger and more numerous fibroids than European American women. We conducted a multi-stage genome-wide association study (GWAS) of fibroid risk among AA women followed by in silico genetically predicted gene expression profiling of top hits. In Stage 1, cases and controls were confirmed by pelvic imaging, genotyped and imputed to 1000 Genomes. Stage 2 used self-reported fibroid and GWAS data from 23andMe, Inc. and the Black Women's Health Study. Associations with fibroid risk were modeled using logistic regression adjusted for principal components, followed by meta-analysis of results. We observed a significant association among 3399 AA cases and 4764 AA controls at rs739187 (risk-allele frequency = 0.27) in CYTH4 (OR (95% confidence interval) = 1.23 (1.16-1.30), p value = 7.82 × 10-9). Evaluation of the genetic association results with MetaXcan identified lower predicted gene expression of CYTH4 in thyroid tissue as significantly associated with fibroid risk (p value = 5.86 × 10-8). In this first multi-stage GWAS for fibroids among AA women, we identified a novel risk locus for fibroids within CYTH4 that impacts gene expression in thyroid and has potential biological relevance for fibroids.

Reconstructing Prehistoric African Population Structure

We assembled genome-wide data from 16 prehistoric Africans. We show that the anciently divergent lineage that comprises the primary ancestry of the southern African San had a wider distribution in the past, contributing approximately two-thirds of the ancestry of Malawi hunter-gatherers ∼8,100-2,500 years ago and approximately one-third of the ancestry of Tanzanian hunter-gatherers ∼1,400 years ago. We document how the spread of farmers from western Africa involved complete replacement of local hunter-gatherers in some regions, and we track the spread of herders by showing that the population of a ∼3,100-year-old pastoralist from Tanzania contributed ancestry to people from northeastern to southern Africa, including a ∼1,200-year-old southern African pastoralist. The deepest diversifications of African lineages were complex, involving either repeated gene flow among geographically disparate groups or a lineage more deeply diverging than that of the San contributing more to some western African populations than to others. We finally leverage ancient genomes to document episodes of natural selection in southern African populations. PAPERCLIP.

The promise of discovering population-specific disease-associated genes in South Asia

The more than 1.5 billion people who live in South Asia are correctly viewed not as a single large population but as many small endogamous groups. We assembled genome-wide data from over 2,800 individuals from over 260 distinct South Asian groups. We identified 81 unique groups, 14 of which had estimated census sizes of more than 1 million, that descend from founder events more extreme than those in Ashkenazi Jews and Finns, both of which have high rates of recessive disease due to founder events. We identified multiple examples of recessive diseases in South Asia that are the result of such founder events. This study highlights an underappreciated opportunity for decreasing disease burden among South Asians through discovery of and testing for recessive disease-associated genes.

Genetic origins of the Minoans and Mycenaeans

The origins of the Bronze Age Minoan and Mycenaean cultures have puzzled archaeologists for more than a century. We assembled genome-wide data from nineteen ancient individuals, including Minoans from Crete, Mycenaeans from mainland Greece, and their eastern neighbours from southwestern Anatolia. We show that Minoans and Mycenaeans were genetically similar, having at least three quarters of their ancestry from the first Neolithic farmers of western Anatolia and the Aegean^1,2, and most of the remainder from ancient populations like those of the Caucasus³ and Iran^4,5. However, the Mycenaeans differed from Minoans in deriving additional ancestry from an ultimate source related to the hunter-gatherers of eastern Europe and Siberia^6–8, introduced via a proximal source related to either the inhabitants of either the Eurasian steppe^1,6,9 or Armenia^4,9. Modern Greeks resemble the Mycenaeans, but with some additional dilution of the early Neolithic ancestry. Our results support the idea of continuity but not isolation in the history of populations of the Aegean, before and after the time of its earliest civilizations.

A Case Report of Cystic Pheochromocytoma

BACKGROUND Pheochromocytoma is a rare catecholamine-producing tumor with an estimated incidence of less than 0.1% in the global population. We present a case of cystic pheochromocytoma that was diagnosed as an incidental finding. The patient presented with abdominal pain and had a history of hypertension. CASE REPORT A 64-year-old man with hypertension presented with a clinical history of intermittent abdominal pain for one year. He denied sweating, palpitations, headache or back pain. He was found to have an elevated blood pressure of 170/90 and no palpable abdominal mass. Contrast-enhanced computed tomography (CT) imaging of the abdomen and pelvis were performed that showed cystic mass measuring 9 cm in diameter arising from the left adrenal gland with contrast-enhancing mural nodules. Magnetic resonance imaging (MRI) confirmed the cystic nature of the mass. Laboratory analysis showed an elevated plasma normetanephrine (NMN) of 1,087 pg/ml and metanephrine (MN) of 372 pg/ml; 24-hour urine showed elevated levels of NMN and MN, 3,002 mg/24 h and 1,596 mg/24 h, respectively. Given the laboratory and radiologic findings, a diagnosis of cystic pheochromocytoma was made. After controlling blood pressure with the alpha-blocker, doxazosin, the patient was hydrated and scheduled for an elective adrenalectomy. The histopathology of the excised adrenal gland was consistent with a cystic pheochromocytoma. CONCLUSIONS Cystic pheochromocytoma is a very rare tumor that may present without symptoms. The clinical course of cystic pheochromocytoma is similar to that of solid pheochromocytoma. Early surgical intervention is recommended, following blood pressure control with an alpha-blocker, and adequate hydration.

Palaeogenomes of Eurasian straight-tusked elephants challenge the current view of elephant evolution

The straight-tusked elephants Palaeoloxodon spp. were widespread across Eurasia during the Pleistocene. Phylogenetic reconstructions using morphological traits have grouped them with Asian elephants (Elephas maximus), and many paleontologists place Palaeoloxodon within Elephas. Here, we report the recovery of full mitochondrial genomes from four and partial nuclear genomes from two P. antiquus fossils. These fossils were collected at two sites in Germany, Neumark-Nord and Weimar-Ehringsdorf, and likely date to interglacial periods ~120 and ~244 thousand years ago, respectively. Unexpectedly, nuclear and mitochondrial DNA analyses suggest that P. antiquus was a close relative of extant African forest elephants (Loxodonta cyclotis). Species previously referred to Palaeoloxodon are thus most parsimoniously explained as having diverged from the lineage of Loxodonta, indicating that Loxodonta has not been constrained to Africa. Our results demonstrate that the current picture of elephant evolution is in need of substantial revision.

A Working Model of the Deep Relationships of Diverse Modern Human Genetic Lineages Outside of Africa

A major topic of interest in human prehistory is how the large-scale genetic structure of modern populations outside of Africa was established. Demographic models have been developed that capture the relationships among small numbers of populations or within particular geographical regions, but constructing a phylogenetic tree with gene flow events for a wide diversity of non-Africans remains a difficult problem. Here, we report a model that provides a good statistical fit to allele-frequency correlation patterns among East Asians, Australasians, Native Americans, and ancient western and northern Eurasians, together with archaic human groups. The model features a primary eastern/western bifurcation dating to at least 45,000 years ago, with Australasians nested inside the eastern clade, and a parsimonious set of admixture events. While our results still represent a simplified picture, they provide a useful summary of deep Eurasian population history that can serve as a null model for future studies and a baseline for further discoveries.

Ancestry and demography and descendants of Iron Age nomads of the Eurasian Steppe

During the 1st millennium before the Common Era (BCE), nomadic tribes associated with the Iron Age Scythian culture spread over the Eurasian Steppe, covering a territory of more than 3,500 km in breadth. To understand the demographic processes behind the spread of the Scythian culture, we analysed genomic data from eight individuals and a mitochondrial dataset of 96 individuals originating in eastern and western parts of the Eurasian Steppe. Genomic inference reveals that Scythians in the east and the west of the steppe zone can best be described as a mixture of Yamnaya-related ancestry and an East Asian component. Demographic modelling suggests independent origins for eastern and western groups with ongoing gene-flow between them, plausibly explaining the striking uniformity of their material culture. We also find evidence that significant gene-flow from east to west Eurasia must have occurred early during the Iron Age.

Mitochondrial DNA analysis of eneolithic trypillians from Ukraine reveals neolithic farming genetic roots

The agricultural revolution in Eastern Europe began in the Eneolithic with the Cucuteni-Trypillia culture complex. In Ukraine, the Trypillian culture (TC) existed for over two millennia (ca. 5,400-2,700 BCE) and left a wealth of artifacts. Yet, their burial rituals remain a mystery and to date almost nothing is known about the genetic composition of the TC population. One of the very few TC sites where human remains can be found is a cave called Verteba in western Ukraine. This report presents four partial and four complete mitochondrial genomes from nine TC individuals uncovered in the cave. The results of this analysis, combined with the data from previous reports, indicate that the Trypillian population at Verteba carried, for the most part, a typical Neolithic farmer package of mitochondrial DNA (mtDNA) lineages traced to Anatolian farmers and Neolithic farming groups of central Europe. At the same time, the find of two specimens belonging to haplogroup U8b1 at Verteba can be viewed as a connection of TC with the Upper Paleolithic European populations. At the level of mtDNA haplogroup frequencies, the TC population from Verteba demonstrates a close genetic relationship with population groups of the Funnel Beaker/ Trichterbecker cultural complex from central and northern Europe (ca. 3,950-2,500 BCE).

Differences in the rare variant spectrum among human populations

Mutations occur at vastly different rates across the genome, and populations, leading to differences in the spectrum of segregating polymorphisms. Here, we investigate variation in the rare variant spectrum in a sample of human genomes representing all major world populations. We find at least two distinct signatures of variation. One, consistent with a previously reported signature is characterized by an increased rate of TCC>TTC mutations in people from Western Eurasia and South Asia, likely related to differences in the rate, or efficiency of repair, of damage due to deamination of methylated guanine. We describe the geographic extent of this signature and show that it is detectable in the genomes of ancient, but not archaic humans. The second signature is private to certain Native American populations, and is concentrated at CpG sites. We show that this signature is not driven by differences in the CpG mutation rate, but is a result of the fact that highly mutable CpG sites are more likely to undergo multiple independent mutations across human populations, and the spectrum of such mutations is highly sensitive to recent demography. Both of these effects dramatically affect the spectrum of rare variants across human populations, and should be taken into account when using mutational clocks to make inference about demography.