Patterns of coding variation in the complete exomes of three Neandertals

Middle and Late Pleistocene Denisovan subsistence at Baishiya Karst Cave

Genetic and fragmented palaeoanthropological data suggest that Denisovans were once widely distributed across eastern Eurasia1-3. Despite limited archaeological evidence, this indicates that Denisovans were capable of adapting to a highly diverse range of environments. Here we integrate zooarchaeological and proteomic analyses of the late Middle to Late Pleistocene faunal assemblage from Baishiya Karst Cave on the Tibetan Plateau, where a Denisovan mandible and Denisovan sedimentary mitochondrial DNA were found3,4. Using zooarchaeology by mass spectrometry, we identify a new hominin rib specimen that dates to approximately 48-32 thousand years ago (layer 3). Shotgun proteomic analysis taxonomically assigns this specimen to the Denisovan lineage, extending their presence at Baishiya Karst Cave well into the Late Pleistocene. Throughout the stratigraphic sequence, the faunal assemblage is dominated by Caprinae, together with megaherbivores, carnivores, small mammals and birds. The high proportion of anthropogenic modifications on the bone surfaces suggests that Denisovans were the primary agent of faunal accumulation. The chaîne opératoire of carcass processing indicates that animal taxa were exploited for their meat, marrow and hides, while bone was also used as raw material for the production of tools. Our results shed light on the behaviour of Denisovans and their adaptations to the diverse and fluctuating environments of the late Middle and Late Pleistocene of eastern Eurasia.

The evolutionary dynamics of local adaptations under genetic rescue is determined by mutational load and polygenicity

Inbred populations often suffer from increased mutational load and reduced fitness due to lower efficacy of purifying selection in groups with small effective population sizes. Genetic rescue (GR) is a conservation tool that is studied and deployed with the aim of increasing the fitness of such inbred populations by assisted migration of individuals from closely related outbred populations. The success of GR depends on several factors-such as their demographic history and distribution of dominance effects of mutations-that may vary across populations. While we understand the impact of these factors on the dynamics of GR, their impact on local adaptations remains unclear. To this end, we conduct a population genetics simulation study to evaluate the impact of trait complexity (Mendelian vs. polygenic), dominance effects, and demographic history on the efficacy of GR. We find that the impact on local adaptations depends highly on the mutational load at the time of GR, which is in turn shaped dynamically by interactions between demographic history and dominance effects of deleterious variation. Over time local adaptations are generally restored post-GR, though in the short term they are often compromised in the process of purging deleterious variation. We also show that while local adaptations are almost always fully restored, the degree to which ancestral genetic variation affecting the trait is replaced by donor variation can vary drastically and is especially high for complex traits. Our results provide insights on the impact of GR on trait evolution and considerations for the practical implementation of GR.

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

Measuring the Efficiency of Purging by non-random Mating in Human Populations

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

The genetic changes that shaped Neandertals, Denisovans, and modern humans

Modern human ancestors diverged from the ancestors of Neandertals and Denisovans about 600,000 years ago. Until about 40,000 years ago, these three groups existed in parallel, occasionally met, and exchanged genes. A critical question is why modern humans, and not the other two groups, survived, became numerous, and developed complex cultures. Here, we discuss genetic differences among the groups and some of their functional consequences. As more present-day genome sequences become available from diverse groups, we predict that very few, if any, differences will distinguish all modern humans from all Neandertals and Denisovans. We propose that the genetic basis of what constitutes a modern human is best thought of as a combination of genetic features, where perhaps none of them is present in each and every present-day individual.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

Towards correlative archaeology of the human mind

Retracing human cognitive origins started out at the systems level with the top-down interpretation of archaeological records spanning from man-made artifacts to endocasts of ancient skulls. With emerging evolutionary genetics and organoid technologies, it is now possible to deconstruct evolutionary processes on a molecular/cellular level from the bottom-up by functionally testing archaic alleles in experimental models. The current challenge is to complement these approaches with novel strategies that allow a holistic reconstruction of evolutionary patterns across human cognitive domains. We argue that computational neuroarcheology can provide such a critical mesoscale framework at the brain network-level, linking molecular/cellular (bottom-up) to systems (top-down) level data for the correlative archeology of the human mind.

Pharmacogenetic Variation in Neanderthals and Denisovans and Implications for Human Health and Response to Medications

Modern humans carry both Neanderthal and Denisovan (archaic) genome elements that are part of the human gene pool and affect the life and health of living individuals. The impact of archaic DNA may be particularly evident in pharmacogenes-genes responsible for the processing of exogenous substances such as food, pollutants, and medications-as these can relate to changing environmental effects, and beneficial variants may have been retained as modern humans encountered new environments. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remain understudied. Here, we explore 11 key cytochrome P450 genes (CYP450) involved in 75% of all drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these 11 CYP450 genes in archaic individuals and find important predicted phenotypic differences relative to modern human variants. We identify several single nucleotide variants shared between archaic and modern humans in each gene, including some potentially function-altering mutations in archaic CYP450 genes, which may result in altered metabolism in living people carrying these variants. We also identified several variants in the archaic CYP450 genes that are novel and unique to archaic humans as well as one gene, CYP2B6, that shows evidence for a gene duplication found only in Neanderthals and modern Africans. Finally, we highlight CYP2A6, CYP2C9, and CYP2J2, genes which show evidence for archaic introgression into modern humans and posit evolutionary hypotheses that explain their allele frequencies in modern populations.

A Neanderthal/Denisovan GLI3 variant contributes to anatomical variations in mice

Changes in genomic structures underlie phenotypic diversification in organisms. Amino acid-changing mutations affect pleiotropic functions of proteins, although little is known about how mutated proteins are adapted in existing developmental programs. Here we investigate the biological effects of a variant of the GLI3 transcription factor (GLI3R1537C) carried in Neanderthals and Denisovans, which are extinct hominins close to modern humans. R1537C does not compromise protein stability or GLI3 activator-dependent transcriptional activities. In contrast, R1537C affects the regulation of downstream target genes associated with developmental processes. Furthermore, genome-edited mice carrying the Neanderthal/Denisovan GLI3 mutation exhibited various alterations in skeletal morphology. Our data suggest that an extinct hominin-type GLI3 contributes to species-specific anatomical variations, which were tolerated by relaxed constraint in developmental programs during human evolution.

Allelic bias when performing in-solution enrichment of ancient human DNA

In-solution hybridisation enrichment of genetic variation is a valuable methodology in human paleogenomics. It allows enrichment of endogenous DNA by targeting genetic markers that are comparable between sequencing libraries. Many studies have used the 1240k reagent-which enriches 1,237,207 genome-wide SNPs-since 2015, though access was restricted. In 2021, Twist Biosciences and Daicel Arbor Biosciences independently released commercial kits that enabled all researchers to perform enrichments for the same 1240 k SNPs. We used the Daicel Arbor Biosciences Prime Plus kit to enrich 132 ancient samples from three continents. We identified a systematic assay bias that increases genetic similarity between enriched samples and that cannot be explained by batch effects. We present the impact of the bias on population genetics inferences (e.g. Principal Components Analysis, ƒ-statistics) and genetic relatedness (READ). We compare the Prime Plus bias to that previously reported of the legacy 1240k enrichment assay. In ƒ-statistics, we find that all Prime-Plus-generated data exhibit artefactual excess shared drift, such that within-continent relationships cannot be correctly determined. The bias is more subtle in READ, though interpretation of the results can still be misleading in specific contexts. We expect the bias may affect analyses we have not yet tested. Our observations support previously reported concerns for the integration of different data types in paleogenomics. We also caution that technological solutions to generate 1240k data necessitate a thorough validation process before their adoption in the paleogenomic community.

Neanderthal coexistence with Homo sapiens in Europe was affected by herbivore carrying capacity

It has been proposed that climate change and the arrival of modern humans in Europe affected the disappearance of Neanderthals due to their impact on trophic resources; however, it has remained challenging to quantify the effect of these factors. By using Bayesian age models to derive the chronology of the European Middle to Upper Paleolithic transition, followed by a dynamic vegetation model that provides the Net Primary Productivity, and a macroecological model to compute herbivore abundance, we show that in continental regions where the ecosystem productivity was low or unstable, Neanderthals disappeared before or just after the arrival of Homo sapiens. In contrast, regions with high and stable productivity witnessed a prolonged coexistence between both species. The temporal overlap between Neanderthals and H. sapiens is significantly correlated with the carrying capacity of small- and medium-sized herbivores. These results suggest that herbivore abundance released the trophic pressure of the secondary consumers guild, which affected the coexistence likelihood between both human species.

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.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans

Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.

NSD1 gene evolves under episodic selection within primates and mutations of specific exons in humans cause Sotos syndrome

BACKGROUND

Modern human brains and skull shapes differ from other hominids. Brain growth disorders as micro- (ASPM, MCPH1) and macrocephaly (NFIX, GLI3) have been highlighted as relevant for the evolution in humans due to the impact in early brain development. Genes associated with macrocephaly have been reported to cause this change, for example NSD1 which causes Sotos syndrome.

RESULTS

In this study we performed a systematic literature review, located the reported variants associated to Sotos syndrome along the gene domains, compared the sequences with close primates, calculated their similarity, Ka/Ks ratios, nucleotide diversity and selection, and analyzed the sequence and structural conservation with distant primates. We aimed to understand if NSD1 in humans differs from other primates since the evolution of NSD1 has not been analyzed in primates, nor if the localization of the mutations is limited to humans. Our study found that most variations causing Sotos syndrome are in exon 19, 22 and 10. In the primate comparison we did not detect Ka/Ks ratios > 1, but a high nucleotide diversity with non-synonymous variations in exons 10, 5, 9, 11 and 23, and sites under episodic selection in exon 5 and 23, and human, macaque/colobus/tarsier/galago and tarsier/lemur/colobus. Most of the domains are conserved in distant primates with a particular progressive development from a simple PWWP1 in O. garnetti to a complex structure in Human.

CONCLUSION

NSD1 is a chromatin modifier that suggests that the selection could influence brain development during modern human evolution and is not present in other primates; however, nowadays the nucleotide diversity is associated with Sotos syndrome.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

SLITRK1-mediated noradrenergic projection suppression in the neonatal prefrontal cortex

SLITRK1 is an obsessive-compulsive disorder spectrum-disorders-associated gene that encodes a neuronal transmembrane protein. Here we show that SLITRK1 suppresses noradrenergic projections in the neonatal prefrontal cortex, and SLITRK1 functions are impaired by SLITRK1 mutations in patients with schizophrenia (S330A, a revertant of Homo sapiens-specific residue) and bipolar disorder (A444S). Slitrk1-KO newborns exhibit abnormal vocalizations, and their prefrontal cortices show excessive noradrenergic neurites and reduced Semaphorin3A expression, which suppresses noradrenergic neurite outgrowth in vitro. Slitrk1 can bind Dynamin1 and L1 family proteins (Neurofascin and L1CAM), as well as suppress Semaphorin3A-induced endocytosis. Neurofascin-binding kinetics is altered in S330A and A444S mutations. Consistent with the increased obsessive-compulsive disorder prevalence in males in childhood, the prefrontal cortex of male Slitrk1-KO newborns show increased noradrenaline levels, and serotonergic varicosity size. This study further elucidates the role of noradrenaline in controlling the development of the obsessive-compulsive disorder-related neural circuit.

Molecular archaeology of human cognitive traits

The brains and minds of our human ancestors remain inaccessible for experimental exploration. Therefore, we reconstructed human cognitive evolution by projecting nonsynonymous/synonymous rate ratios (ω values) in mammalian phylogeny onto the anatomically modern human (AMH) brain. This atlas retraces human neurogenetic selection and allows imputation of ancestral evolution in task-related functional networks (FNs). Adaptive evolution (high ω values) is associated with excitatory neurons and synaptic function. It shifted from FNs for motor control in anthropoid ancestry (60-41 mya) to attention in ancient hominoids (26-19 mya) and hominids (19-7.4 mya). Selection in FNs for language emerged with an early hominin ancestor (7.4-1.7 mya) and was later accompanied by adaptive evolution in FNs for strategic thinking during recent (0.8 mya-present) speciation of AMHs. This pattern mirrors increasingly complex cognitive demands and suggests that co-selection for language alongside strategic thinking may have separated AMHs from their archaic Denisovan and Neanderthal relatives.

Oxytocin and vasotocin receptor variation and the evolution of human prosociality

Modern human lifestyle strongly depends on complex social traits like empathy, tolerance and cooperation. These diverse facets of social cognition have been associated with variation in the oxytocin receptor (OTR) and its sister genes, the vasotocin/vasopressin receptors (VTR1A/AVPR1A and AVPR1B/VTR1B). Here, we compared the available genomic sequences of these receptors between modern humans, archaic humans, and 12 non-human primate species, and identified sites that show heterozygous variation in modern humans and archaic humans distinct from variation in other primates, and for which we could find association studies with clinical implications. On these sites, we performed a range of analyses (variant clustering, pathogenicity prediction, regulation, linkage disequilibrium frequency), and reviewed the literature on selection data in different modern-human populations. We found five sites with modern human specific variation, where the modern human allele is the major allele in the global population (OTR: rs1042778, rs237885, rs6770632; VTR1A: rs10877969; VTR1B: rs33985287). Among them, variation in the OTR-rs6770632 site was predicted to be the most functional. Two alleles (OTR: rs59190448 and rs237888) present only in modern humans and archaic humans were putatively under positive selection in modern humans, with rs237888 predicted to be a highly functional site. Three sites showed convergent evolution between modern humans and bonobos (OTR: rs2228485 and rs237897; VTR1A: rs1042615), with OTR-rs2228485 ranking highly in terms of functionality and reported to be under balancing selection in modern humans (Schaschl, 2015) [1]. Our findings have implications for understanding hominid prosociality, as well as the similarities between modern human and bonobo social behavior.

Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders

Modern humans exhibit phenotypic traits and molecular events shared with other domesticates that are thought to be by-products of selection for reduced aggression. This is the human self-domestication hypothesis. As one of the first types of responses to a novel environment, epigenetic changes may have also facilitated early self-domestication in humans. Here, we argue that fish species, which have been recently domesticated, can provide model systems to study epigenetic drivers in human self-domestication. To test this, we used in silico approaches to compare genes with epigenetic changes in early domesticates of European sea bass with genes exhibiting methylation changes in anatomically modern humans (comparison 1), and neurodevelopmental cognitive disorders considered to exhibit abnormal self-domestication traits, i.e., schizophrenia, Williams syndrome, and autism spectrum disorders (comparison 2). Overlapping genes in comparison 1 were involved in processes like limb morphogenesis and phenotypes like abnormal jaw morphology and hypopigmentation. Overlapping genes in comparison 2 affected paralogue genes involved in processes such as neural crest differentiation and ectoderm differentiation. These findings pave the way for future studies using fish species as models to investigate epigenetic changes as drivers of human self-domestication and as triggers of cognitive disorders.

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Evolution of genetic networks for human creativity

The genetic basis for the emergence of creativity in modern humans remains a mystery despite sequencing the genomes of chimpanzees and Neanderthals, our closest hominid relatives. Data-driven methods allowed us to uncover networks of genes distinguishing the three major systems of modern human personality and adaptability: emotional reactivity, self-control, and self-awareness. Now we have identified which of these genes are present in chimpanzees and Neanderthals. We replicated our findings in separate analyses of three high-coverage genomes of Neanderthals. We found that Neanderthals had nearly the same genes for emotional reactivity as chimpanzees, and they were intermediate between modern humans and chimpanzees in their numbers of genes for both self-control and self-awareness. 95% of the 267 genes we found only in modern humans were not protein-coding, including many long-non-coding RNAs in the self-awareness network. These genes may have arisen by positive selection for the characteristics of human well-being and behavioral modernity, including creativity, prosocial behavior, and healthy longevity. The genes that cluster in association with those found only in modern humans are over-expressed in brain regions involved in human self-awareness and creativity, including late-myelinating and phylogenetically recent regions of neocortex for autobiographical memory in frontal, parietal, and temporal regions, as well as related components of cortico-thalamo-ponto-cerebellar-cortical and cortico-striato-cortical loops. We conclude that modern humans have more than 200 unique non-protein-coding genes regulating co-expression of many more protein-coding genes in coordinated networks that underlie their capacities for self-awareness, creativity, prosocial behavior, and healthy longevity, which are not found in chimpanzees or Neanderthals.

Whole-exome sequencing of the mummified remains of Cangrande della Scala (1291-1329 CE) indicates the first known case of late-onset Pompe disease

Mummified remains of relevant historical figures are nowadays an important source of information to retrace data concerning their private life and health, especially when historical archives are not available. Next-generation-sequencing was proved to be a valuable tool to unravel the characteristics of these individuals through their genetic heritage. Using the strictest criteria currently available for the validation of ancient DNA sequences, whole-genome and whole-exome sequencing were generated from the mummy remains of an Italian nobleman died almost 700 years ago, Cangrande della Scala. While its genome sequencing could not yield sufficient coverage for in depth investigation, exome sequencing could overcome the limitations of this approach to achieve significantly high coverage on coding regions, thus allowing to perform the first extensive exome analysis of a mummy genome. Similar to a standard "clinical exome analysis" conducted on modern DNA, an in-depth variant annotation, high-quality filtering and interpretation was performed, leading to the identification of a genotype associated with late-onset Pompe disease (glycogen storage disease type II). This genetic diagnosis was concordant with the limited clinical history available for Cangrande della Scala, who likely represents the earliest known case of this autosomal recessive metabolic disorder.

An etiology of human modernity

Following the refutation of the replacement hypothesis, which had proposed that a ‘superior’ hominin species arose in Africa and replaced all other humans existing at the time, the auto-domestication hypothesis remains the only viable explanation for the relatively abrupt change from robust to gracile humans in the Late Pleistocene. It invokes the incidental institution of the domestication syndrome in humans, most probably by newly introduced cultural practices. It also postulates that the induction of exograms compensated for the atrophy of the brain caused by domestication. This new explanation of the origins of modernity in humans elucidates practically all its many aspects, in stark contrast to the superseded replacement hypothesis, which explained virtually nothing. The first results of the domestication syndrome’s genetic exploration have become available in recent years, and they endorse the human self-domestication hypothesis.

Multiple Genomic Events Altering Hominin SIGLEC Biology and Innate Immunity Predated the Common Ancestor of Humans and Archaic Hominins

Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface “self-associated molecular patterns” that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2–3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a “hotspot” in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human–Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal–human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

Deleterious variants in genes regulating mammalian reproduction in Neanderthals, Denisovans and extant humans

STUDY QUESTION

Were Neanderthals and Denisovans (referred here also as extinct hominidae) carrying deleterious variants in genes regulating reproduction?

SUMMARY ANSWER

The majority of extinct hominidae analyzed here, presented a considerable number of deleterious variants per individual in proteins regulating different aspects of reproduction, including gonad and uterine function, and gametogenesis.

WHAT IS KNOWN ALREADY

Neanderthals, Denisovans and extant humans were interfertile and hybridized while occupying geographically overlapping areas in Europe and Asia. This is evidenced by the small archaic genome component (average ∼2%) present in non-African extant humans.

STUDY DESIGN, SIZE, DURATION

The genome of eight extinct hominidae, together with five human genome databases, plus 44 mothers and 48 fathers (fertile controls), were screened to look for deleterious variants in 1734 protein-coding genes regulating reproduction.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Ancient DNA from six Neanderthals and two Denisovans dated between ∼82 000 and 43 000 calibrated years was retrieved from the public European Nucleotide Archive. The hominins analyzed include Altai, Vindija 33.15, 33.19, 33.25 and 33.26, El Sidron 1253, Denisova 3 and 11. Their DNA was analyzed using the CLC Genomics Workbench 12, by mapping overlapping paired-end reads (Illumina, FASTQ files) to the human genome assembly GRCh37 (hg19) (Vindija 33.19, 33.25, 33.26, Denisova 3 and Denisova 11) or by analyzing BAM files (Altai, El Sidron 1253 and Vindija 33.15) (human genome reference, GRCh37 (hg19)). Non-synonymous reproductive variants were classified as deleterious or tolerated (PolyPhen-2 and SIFT analyses) and were compared to deleterious variants obtained from extant human genome databases (Genome Aggregation Database (GnomAD), 1000 Genomes, the Haplotype Map (HapMap), Single Nucleotide Polymorphism Database (dbSNPs)) across different populations. A genetic intersection between extant or extinct DNA variants and other genetic disorders was evaluated by annotating the obtained variants with the Clinical Variant (ClinVar) database.

MAIN RESULTS AND THE ROLE OF CHANCE

Among the eight extinct hominidae analyzed, a total of 9650 non-synonymous variants (only coverage ≥20 reads included; frameshift mutations were excluded) in 1734 reproductive protein-coding genes were found, 24% of which were classified as deleterious. The majority (73%) of the deleterious alleles present in extant humans that are shared between extant humans and extinct hominidae were found to be rare (<1%) in extant human populations. A set of 8044 variants were found uniquely in extinct hominidae. At the single-gene level, no extinct individual was found to be homozygous for deleterious variants in genes necessary for gamete recognition and fusion, and no higher chance of embryo-lethality (calculated by Mendelian Genetics) was found upon simulated mating between extant human and extinct hominidae compared to extant human-extant human. However, three of the eight extinct hominidae were found to be homozygous for 48-69 deleterious variants in 55 genes controlling ovarian and uterine functions, or oogenesis (AKAP1, BUB1B, CCDC141, CDC73, DUSP6, ESR1, ESR2, PATL2, PSMC3IP, SEMA3A, WT1 and WNT4). Moreover, we report the distribution of nine Neanderthal variants in genes associated with a human fertility phenotype found in extant human populations, one of which has been associated with polycystic ovarian syndrome and primary congenital glaucoma.

LIMITATIONS, REASONS FOR CAUTION

While analyzing archaic DNA, stringent filtering criteria were adopted to screen for deleterious variants in Neanderthals and Denisovans, which could result in missing a number of variants. Such restraints preserve the potential for detection of additional deleterious variants in reproductive proteins in extinct hominidae.

WIDER IMPLICATIONS OF THE FINDINGS

This study provides a comprehensive overview of putatively deleterious variants in extant human populations and extinct individuals occurring in 1734 protein-coding genes controlling reproduction and provides the fundaments for future functional studies of extinct variants in human reproduction.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the Department of Biological Science and by the Office of Research and Sponsored Programs at the University of Tulsa (Faculty Research Grant and Faculty Research Summer Fellowship) to M.A. and the University of Tulsa, Tulsa Undergraduate Research Challenge (TURC) program to E.L.; no conflict of interest to declare.

TRIAL REGISTRATION NUMBER

N/A.

Reduced purine biosynthesis in humans after their divergence from Neandertals

We analyze the metabolomes of humans, chimpanzees, and macaques in muscle, kidney and three different regions of the brain. Although several compounds in amino acid metabolism occur at either higher or lower concentrations in humans than in the other primates, metabolites downstream of adenylosuccinate lyase, which catalyzes two reactions in purine synthesis, occur at lower concentrations in humans. This enzyme carries an amino acid substitution that is present in all humans today but absent in Neandertals. By introducing the modern human substitution into the genomes of mice, as well as the ancestral, Neandertal-like substitution into the genomes of human cells, we show that this amino acid substitution contributes to much or all of the reduction of de novo synthesis of purines in humans.

Neandertal introgression and accumulation of hypomorphic mutations in the neuropeptide S (NPS) system promote attenuated functionality

The neuropeptide S (NPS) system plays an important role in fear and fear memory processing but has also been associated with allergic and inflammatory diseases. Genes for NPS and its receptor NPSR1 are found in all tetrapods. Compared to non-human primates, several non-synonymous single-nucleotide polymorphisms (SNPs) occur in both human genes that collectively result in functional attenuation, suggesting adaptive mechanisms in a human context. To investigate historic and geographic origins of these hypomorphic mutations and explore genetic signs of selection, we analyzed ancient genomes and worldwide genotype frequencies of four prototypic SNPs in the NPS system. Neandertal and Denisovan genomes contain exclusively ancestral alleles for NPSR1 while all derived alleles occur in ancient genomes of anatomically modern humans, indicating that they arose in modern Homo sapiens. Worldwide genotype frequencies for three hypomorphic NPSR1 SNPs show significant regional homogeneity but follow a gradient towards increasing derived allele frequencies that supports an out-of-Africa scenario. Increased density of high-frequency polymorphisms around the three NPSR1 loci suggests weak or possibly balancing selection. A hypomorphic mutation in the NPS precursor, however, was detected at high frequency in Eurasian Neandertal genomes and shows genetic signatures indicating that it was introgressed into the human gene pool, particularly in Southern Europe, by interbreeding with Neandertals. We discuss potential evolutionary scenarios including behavior and immune-based natural selection.

An ecological niche shift for Neanderthal populations in Western Europe 70,000 years ago

Middle Paleolithic Neanderthal populations occupied Eurasia for at least 250,000 years prior to the arrival of anatomically modern humans. While a considerable body of archaeological research has focused on Neanderthal material culture and subsistence strategies, little attention has been paid to the relationship between regionally specific cultural trajectories and their associated existing fundamental ecological niches, nor to how the latter varied across periods of climatic variability. We examine the Middle Paleolithic archaeological record of a naturally constrained region of Western Europe between 82,000 and 60,000 years ago using ecological niche modeling methods. Evaluations of ecological niche estimations, in both geographic and environmental dimensions, indicate that 70,000 years ago the range of suitable habitats exploited by these Neanderthal populations contracted and shifted. These ecological niche dynamics are the result of groups continuing to occupy habitual territories that were characterized by new environmental conditions during Marine Isotope Stage 4. The development of original cultural adaptations permitted this territorial stability.

Systematic detection of brain protein-coding genes under positive selection during primate evolution and their roles in cognition

The human brain differs from that of other primates, but the genetic basis of these differences remains unclear. We investigated the evolutionary pressures acting on almost all human protein-coding genes (N = 11,667; 1:1 orthologs in primates) based on their divergence from those of early hominins, such as Neanderthals, and non-human primates. We confirm that genes encoding brain-related proteins are among the most strongly conserved protein-coding genes in the human genome. Combining our evolutionary pressure metrics for the protein-coding genome with recent data sets, we found that this conservation applied to genes functionally associated with the synapse and expressed in brain structures such as the prefrontal cortex and the cerebellum. Conversely, several genes presenting signatures commonly associated with positive selection appear as causing brain diseases or conditions, such as micro/macrocephaly, Joubert syndrome, dyslexia, and autism. Among those, a number of DNA damage response genes associated with microcephaly in humans such as BRCA1, NHEJ1, TOP3A, and RNF168 show strong signs of positive selection and might have played a role in human brain size expansion during primate evolution. We also showed that cerebellum granule neurons express a set of genes also presenting signatures of positive selection and that may have contributed to the emergence of fine motor skills and social cognition in humans. This resource is available online and can be used to estimate evolutionary constraints acting on a set of genes and to explore their relative contributions to human traits.

Theoretical and Methodological Approaches to Ecological Changes, Social Behaviour and Human Intergroup Tolerance 300,000 to 30,000 BP

Archaeological evidence suggests that important shifts were taking place in the character of human social behaviours 300,000 to 30,000 years ago. New artefact types appear and are disseminated with greater frequency. Transfers of both raw materials and finished artefacts take place over increasing distances, implying larger scales of regional mobility and more frequent and friendlier interactions between different communities. Whilst these changes occur during a period of increasing environmental variability, the relationship between ecological changes and transformations in social behaviours is elusive. Here, we explore a possible theoretical approach and methodology for understanding how ecological contexts can influence selection pressures acting on intergroup social behaviours. We focus on the relative advantages and disadvantages of intergroup tolerance in different ecological contexts using agent-based modelling (ABM). We assess the relative costs and benefits of different 'tolerance' levels in between-group interactions on survival and resource exploitation in different environments. The results enable us to infer a potential relationship between ecological changes and proposed changes in between-group behavioural dynamics. We conclude that increasingly harsh environments may have driven changes in hormonal and emotional responses in humans leading to increasing intergroup tolerance, i.e. transformations in social behaviour associated with 'self-domestication'. We argue that changes in intergroup tolerance is a more parsimonious explanation for the emergence of what has been seen as 'modern human behaviour' than changes in hard aspects of cognition or other factors such as cognitive adaptability or population size.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10816-020-09503-5.

The evolutionary history of Neanderthal and Denisovan Y chromosomes

Ancient DNA has provided new insights into many aspects of human history. However, we lack comprehensive studies of the Y chromosomes of Denisovans and Neanderthals because the majority of specimens that have been sequenced to sufficient coverage are female. Sequencing Y chromosomes from two Denisovans and three Neanderthals shows that the Y chromosomes of Denisovans split around 700 thousand years ago from a lineage shared by Neanderthals and modern human Y chromosomes, which diverged from each other around 370 thousand years ago. The phylogenetic relationships of archaic and modern human Y chromosomes differ from the population relationships inferred from the autosomal genomes and mirror mitochondrial DNA phylogenies, indicating replacement of both the mitochondrial and Y chromosomal gene pools in late Neanderthals. This replacement is plausible if the low effective population size of Neanderthals resulted in an increased genetic load in Neanderthals relative to modern humans.

Neanderthal introgression reintroduced functional ancestral alleles lost in Eurasian populations

Neanderthal ancestry remains across modern Eurasian genomes and introgressed sequences influence diverse phenotypes. Here, we demonstrate that introgressed sequences reintroduced thousands of ancestral alleles that were lost in Eurasian populations before introgression. Our simulations and variant effect predictions argue that these reintroduced alleles (RAs) are more likely to be tolerated by modern humans than are introgressed Neanderthal-derived alleles (NDAs) due to their distinct evolutionary histories. Consistent with this, we show enrichment for RAs and depletion for NDAs on introgressed haplotypes with expression quantitative trait loci (eQTL) and phenotype associations. Analysis of available cross-population eQTLs and massively parallel reporter assay data show that RAs commonly influence gene expression independent of linked NDAs. We further validate these independent effects for one RA in vitro. Finally, we demonstrate that NDAs are depleted for regulatory activity compared to RAs, while RAs have activity levels similar to non-introgressed variants. In summary, our study reveals that Neanderthal introgression reintroduced thousands of lost ancestral variants with gene regulatory activity and that these RAs were more tolerated than NDAs. Thus, RAs and their distinct evolutionary histories must be considered when evaluating the effects of introgression.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.

Triangulating Neanderthal cognition: A tale of not seeing the forest for the trees

The inference of Neanderthal cognition, including their cultural and linguistic capabilities, has persisted as a fiercely debated research topic for decades. This lack of consensus is substantially based on inherent uncertainties in reconstructing prehistory out of indirect evidence as well as other methodological limitations. Further factors include systemic difficulties within interdisciplinary discourse, data artifacts, historic research biases, and the sheer scope of the relevant research. Given the degrees of freedom in interpretation ensuing from these complications, any attempt to find approximate answers to the yet unsettled pertinent discourse may not rest on single studies, but instead a careful and comprehensive interdisciplinary synthesis of findings. Triangulating Neanderthals' cognition by considering the plethora of data, diverse perspectives and aforementioned complexities present within the literature constitutes the currently most reliable pathway to tentative conclusions. While some uncertainties remain, such an approach paints the picture of an extensive shared humanity between anatomically modern humans and Neanderthals. This article is categorized under: Cognitive Biology > Evolutionary Roots of Cognition Linguistics > Evolution of Language.

The Neandertal Progesterone Receptor

The hormone progesterone is important for preparing the uterine lining for egg implantation and for maintaining the early stages of pregnancy. The gene encoding the progesterone receptor (PGR) carries introgressed Neandertal haplotypes with two missense substitutions and a mobile Alu element. These Neandertal gene variants have reached nearly 20% frequency in non-Africans and have been associated with preterm birth. Here, we show that one of the missense substitutions appears fixed in Neandertals, while the other substitution as well as the Alu insertion were polymorphic among Neandertals. We show that two Neandertal haplotypes carrying the PGR gene entered the modern human population and that present-day carriers of the Neandertal haplotypes express higher levels of the receptor. In a cohort of present-day Britons, these carriers have more siblings, fewer miscarriages, and less bleeding during early pregnancy suggesting that the Neandertal progesterone receptor alleles promote fertility. This may explain their high frequency in modern human populations.

Selection against archaic hominin genetic variation in regulatory regions

Traces of Neandertal and Denisovan DNA persist in the modern human gene pool, but have been systematically purged by natural selection from genes and other functionally important regions. This implies that many archaic alleles harmed the fitness of hybrid individuals, but the nature of this harm is poorly understood. Here, we show that enhancers contain less Neandertal and Denisovan variation than expected given the background selection they experience, suggesting that selection acted to purge these regions of archaic alleles that disrupted their gene regulatory functions. We infer that selection acted mainly on young archaic variation that arose in Neandertals or Denisovans shortly before their contact with humans; enhancers are not depleted of older variants found in both archaic species. Some types of enhancer appear to have tolerated introgression better than others; compared with tissue-specific enhancers, pleiotropic enhancers show stronger depletion of archaic single-nucleotide polymorphisms. To some extent, evolutionary constraint is predictive of introgression depletion, but certain tissues' enhancers are more depleted of Neandertal and Denisovan alleles than expected given their comparative tolerance to new mutations. Foetal brain and muscle are the tissues whose enhancers show the strongest depletion of archaic alleles, but only brain enhancers show evidence of unusually stringent purifying selection. We conclude that epistatic incompatibilities between human and archaic alleles are needed to explain the degree of archaic variant depletion from foetal muscle enhancers, perhaps due to divergent selection for higher muscle mass in archaic hominins compared with humans.

Genomic analysis of the natural history of attention-deficit/hyperactivity disorder using Neanderthal and ancient Homo sapiens samples

Attention-deficit/hyperactivity disorder (ADHD) is an impairing neurodevelopmental condition highly prevalent in current populations. Several hypotheses have been proposed to explain this paradox, mainly in the context of the Paleolithic versus Neolithic cultural shift but especially within the framework of the mismatch theory. This theory elaborates on how a particular trait once favoured in an ancient environment might become maladaptive upon environmental changes. However, given the lack of genomic data available for ADHD, these theories have not been empirically tested. We took advantage of the largest GWAS meta-analysis available for this disorder consisting of over 20,000 individuals diagnosed with ADHD and 35,000 controls, to assess the evolution of ADHD-associated alleles in European populations using archaic, ancient and modern human samples. We also included Approximate Bayesian computation coupled with deep learning analyses and singleton density scores to detect human adaptation. Our analyses indicate that ADHD-associated alleles are enriched in loss of function intolerant genes, supporting the role of selective pressures in this early-onset phenotype. Furthermore, we observed that the frequency of variants associated with ADHD has steadily decreased since Paleolithic times, particularly in Paleolithic European populations compared to samples from the Neolithic Fertile Crescent. We demonstrate this trend cannot be explained by African admixture nor Neanderthal introgression, since introgressed Neanderthal alleles are enriched in ADHD risk variants. All analyses performed support the presence of long-standing selective pressures acting against ADHD-associated alleles until recent times. Overall, our results are compatible with the mismatch theory for ADHD but suggest a much older time frame for the evolution of ADHD-associated alleles compared to previous hypotheses.

The repertoire of family A-peptide GPCRs in archaic hominins

Given the importance of G-protein coupled receptors in the regulation of many physiological functions, deciphering the relationships between genotype and phenotype in past and present hominin GPCRs is of main interest to understand the evolutionary process that contributed to the present-day variability in human traits and health. Here, we carefully examined the publicly available genomic and protein sequence databases of the archaic hominins (Neanderthal and Denisova) to draw up the catalog of coding variations in GPCRs for peptide ligands, in comparison with living humans. We then searched in the literature the functional changes, phenotypes and risk of disease possibly associated with the detected variants. Our survey suggests that Neanderthal and Denisovan hominins were likely prone to lower risk of obesity, to enhanced platelet aggregation in response to thrombin, to better response to infection, to less anxiety and aggressiveness and to favorable sociability. While some archaic variants were likely advantageous in the past, they might be responsible for maladaptive disorders today in the context of modern life and/or specific regional distribution. For example, an archaic haplotype in the neuromedin receptor 2 is susceptible to confer risk of diabetic nephropathy in type 1 diabetes in present-day Europeans. Paying attention to the pharmacological properties of some of the archaic variants described in this study may be helpful to understand the variability of therapeutic efficacy between individuals or ethnic groups.

Inbreeding, Allee effects and stochasticity might be sufficient to account for Neanderthal extinction

The replacement of Neanderthals by Anatomically Modern Humans has typically been attributed to environmental pressure or a superiority of modern humans with respect to competition for resources. Here we present two independent models that suggest that no such heatedly debated factors might be needed to account for the demise of Neanderthals. Starting from the observation that Neanderthal populations already were small before the arrival of modern humans, the models implement three factors that conservation biology identifies as critical for a small population's persistence, namely inbreeding, Allee effects and stochasticity. Our results indicate that the disappearance of Neanderthals might have resided in the smallness of their population(s) alone: even if they had been identical to modern humans in their cognitive, social and cultural traits, and even in the absence of inter-specific competition, Neanderthals faced a considerable risk of extinction. Furthermore, we suggest that if modern humans contributed to the demise of Neanderthals, that contribution might have had nothing to do with resource competition, but rather with how the incoming populations geographically restructured the resident populations, in a way that reinforced Allee effects, and the effects of inbreeding and stochasticity.

Disease transmission and introgression can explain the long-lasting contact zone of modern humans and Neanderthals

Neanderthals and modern humans both occupied the Levant for tens of thousands of years prior to the spread of modern humans into the rest of Eurasia and their replacement of the Neanderthals. That the inter-species boundary remained geographically localized for so long is a puzzle, particularly in light of the rapidity of its subsequent movement. Here, we propose that infectious-disease dynamics can explain the localization and persistence of the inter-species boundary. We further propose, and support with dynamical-systems models, that introgression-based transmission of alleles related to the immune system would have gradually diminished this barrier to pervasive inter-species interaction, leading to the eventual release of the inter-species boundary from its geographic localization. Asymmetries between the species in the characteristics of their associated ‘pathogen packages’ could have generated feedback that allowed modern humans to overcome disease burden earlier than Neanderthals, giving them an advantage in their subsequent spread into Eurasia.

Modern humans and Neanderthals coexisted in the Levant for tens of thousands of years before modern humans spread and replaced Neanderthals. Here, Greenbaum et al. develop a model showing that transmission of disease and genes can explain the maintenance and then collapse of this contact zone.

Inferred divergent gene regulation in archaic hominins reveals potential phenotypic differences

Sequencing DNA derived from archaic bones has enabled genetic comparison of Neanderthals and anatomically modern humans (AMHs), and revealed that they interbred. However, interpreting what genetic differences imply about their phenotypic differences remains challenging. Here, we introduce an approach for identifying divergent gene regulation between archaic hominins, such as Neanderthals, and AMH sequences, and find 766 genes that are likely to have been divergently regulated (DR) by Neanderthal haplotypes that do not remain in AMHs. DR genes include many involved in phenotypes known to differ between Neanderthals and AMHs, such as the structure of the rib cage and supraorbital ridge development. They are also enriched for genes associated with spontaneous abortion, polycystic ovary syndrome, myocardial infarction and melanoma. Phenotypes associated with modern human variation in these genes' regulation in ~23,000 biobank patients further support their involvement in immune and cardiovascular phenotypes. Comparing DR genes between two Neanderthals and a Denisovan revealed divergence in the immune system and in genes associated with skeletal and dental morphology that are consistent with the archaeological record. These results establish differences in gene regulatory architecture between AMHs and archaic hominins, and provide an avenue for exploring phenotypic differences between archaic groups from genomic information alone.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

A genetic analysis of the Gibraltar Neanderthals

The Forbes' Quarry and Devil's Tower partial crania from Gibraltar are among the first Neanderthal remains ever found. Here, we show that small amounts of ancient DNA are preserved in the petrous bones of the 2 individuals despite unfavorable climatic conditions. However, the endogenous Neanderthal DNA is present among an overwhelming excess of recent human DNA. Using improved DNA library construction methods that enrich for DNA fragments carrying deaminated cytosine residues, we were able to sequence 70 and 0.4 megabase pairs (Mbp) nuclear DNA of the Forbes' Quarry and Devil's Tower specimens, respectively, as well as large parts of the mitochondrial genome of the Forbes' Quarry individual. We confirm that the Forbes' Quarry individual was a female and the Devil's Tower individual a male. We also show that the Forbes' Quarry individual is genetically more similar to the ∼120,000-y-old Neanderthals from Scladina Cave in Belgium (Scladina I-4A) and Hohlenstein-Stadel Cave in Germany, as well as to a ∼60,000- to 70,000-y-old Neanderthal from Russia (Mezmaiskaya 1), than to a ∼49,000-y-old Neanderthal from El Sidrón (El Sidrón 1253) in northern Spain and other younger Neanderthals from Europe and western Asia. This suggests that the Forbes' Quarry fossil predates the latter Neanderthals. The preservation of archaic human DNA in the warm coastal climate of Gibraltar, close to the shores of Africa, raises hopes for the future recovery of archaic human DNA from regions in which climatic conditions are less than optimal for DNA preservation.

Genetic Variation in Pan Species Is Shaped by Demographic History and Harbors Lineage-Specific Functions

Chimpanzees (Pan troglodytes) and bonobos (Pan paniscus) are the closest living relatives of humans, but the two species show distinct behavioral and physiological differences, particularly regarding female reproduction. Despite their recent rapid decline, the demographic histories of the two species have been different during the past 1–2 Myr, likely having an impact on their genomic diversity. Here, we analyze the inferred functional consequences of genetic variation across 69 individuals, making use of the most complete data set of genomes in the Pan clade to date. We test to which extent the demographic history influences the efficacy of purifying selection in these species. We find that small historical effective population sizes (N_e) correlate not only with low levels of genetic diversity but also with a larger number of deleterious alleles in homozygosity and an increased proportion of deleterious changes at low frequencies. To investigate the putative genetic basis for phenotypic differences between chimpanzees and bonobos, we exploit the catalog of putatively deleterious protein-coding changes in each lineage. We show that bonobo-specific nonsynonymous changes are enriched in genes related to age at menarche in humans, suggesting that the prominent physiological differences in the female reproductive system between chimpanzees and bonobos might be explained, in part, by putatively adaptive changes on the bonobo lineage.