The phenotypic legacy of admixture between modern humans and Neandertals

Microorganisms as Shapers of Human Civilization, from Pandemics to Even Our Genomes: Villains or Friends? A Historical Approach

Universal history is characterized by continuous evolution, in which civilizations are born and die. This evolution is associated with multiple factors, among which the role of microorganisms is often overlooked. Viruses and bacteria have written or decisively contributed to terrible episodes of history, such as the Black Death in 14th century Europe, the annihilation of pre-Columbian American civilizations, and pandemics such as the 1918 Spanish flu or the current COVID-19 pandemic caused by the coronavirus SARS-CoV-2. Nevertheless, it is clear that we could not live in a world without these tiny beings. Endogenous retroviruses have been key to our evolution and for the regulation of gene expression, and the gut microbiota helps us digest compounds that we could not otherwise process. In addition, we have used microorganisms to preserve or prepare food for millennia and more recently to obtain drugs such as antibiotics or to develop recombinant DNA technologies. Due to the enormous importance of microorganisms for our survival, they have significantly influenced the population genetics of different human groups. This paper will review the role of microorganisms as "villains" who have been responsible for tremendous mortality throughout history but also as "friends" who help us survive and evolve.

Populations, Traits, and Their Spatial Structure in Humans

The spatial distribution of genetic variants is jointly determined by geography, past demographic processes, natural selection, and its interplay with environmental variation. A fraction of these genetic variants are "causal alleles" that affect the manifestation of a complex trait. The effect exerted by these causal alleles on complex traits can be independent or dependent on the environment. Understanding the evolutionary processes that shape the spatial structure of causal alleles is key to comprehend the spatial distribution of complex traits. Natural selection, past population size changes, range expansions, consanguinity, assortative mating, archaic introgression, admixture, and the environment can alter the frequencies, effect sizes, and heterozygosities of causal alleles. This provides a genetic axis along which complex traits can vary. However, complex traits also vary along biogeographical and sociocultural axes which are often correlated with genetic axes in complex ways. The purpose of this review is to consider these genetic and environmental axes in concert and examine the ways they can help us decipher the variation in complex traits that is visible in humans today. This initiative necessarily implies a discussion of populations, traits, the ability to infer and interpret "genetic" components of complex traits, and how these have been impacted by adaptive events. In this review, we provide a history-aware discussion on these topics using both the recent and more distant past of our academic discipline and its relevant contexts.

DLX5/6 GABAergic Expression Affects Social Vocalization: Implications for Human Evolution

DLX5 and DLX6 are two closely related transcription factors involved in brain development and in GABAergic differentiation. The DLX5/6 locus is regulated by FoxP2, a gene involved in language evolution and has been associated with neurodevelopmental disorders and mental retardation. Targeted inactivation of Dlx5/6 in mouse GABAergic neurons (Dlx5/6VgatCre mice) results in behavioral and metabolic phenotypes notably increasing lifespan by 33%. Here, we show that Dlx5/6VgatCre mice present a hyper-vocalization and hyper-socialization phenotype. While only 7% of control mice emitted more than 700 vocalizations/10 min, 30% and 56% of heterozygous or homozygous Dlx5/6VgatCre mice emitted more than 700 and up to 1,400 calls/10 min with a higher proportion of complex and modulated calls. Hyper-vocalizing animals were more sociable: the time spent in dynamic interactions with an unknown visitor was more than doubled compared to low-vocalizing individuals. The characters affected by Dlx5/6 in the mouse (sociability, vocalization, skull, and brain shape…) overlap those affected in the "domestication syndrome". We therefore explored the possibility that DLX5/6 played a role in human evolution and "self-domestication" comparing DLX5/6 genomic regions from Neanderthal and modern humans. We identified an introgressed Neanderthal haplotype (DLX5/6-N-Haplotype) present in 12.6% of European individuals that covers DLX5/6 coding and regulatory sequences. The DLX5/6-N-Haplotype includes the binding site for GTF2I, a gene associated with Williams-Beuren syndrome, a hyper-sociability and hyper-vocalization neurodevelopmental disorder. The DLX5/6-N-Haplotype is significantly underrepresented in semi-supercentenarians (>105 years of age), a well-established human model of healthy aging and longevity, suggesting their involvement in the coevolution of longevity, sociability, and speech.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.

Quantifying the contribution of Neanderthal introgression to the heritability of complex traits

Eurasians have ~2% Neanderthal ancestry, but we lack a comprehensive understanding of the genome-wide influence of Neanderthal introgression on modern human diseases and traits. Here, we quantify the contribution of introgressed alleles to the heritability of more than 400 diverse traits. We show that genomic regions in which detectable Neanderthal ancestry remains are depleted of heritability for all traits considered, except those related to skin and hair. Introgressed variants themselves are also depleted for contributions to the heritability of most traits. However, introgressed variants shared across multiple Neanderthal populations are enriched for heritability and have consistent directions of effect on several traits with potential relevance to human adaptation to non-African environments, including hair and skin traits, autoimmunity, chronotype, bone density, lung capacity, and menopause age. Integrating our results, we propose a model in which selection against introgressed functional variation was the dominant trend (especially for cognitive traits); however, for a few traits, introgressed variants provided beneficial variation via uni-directional (e.g., lightening skin color) or bi-directional (e.g., modulating immune response) effects.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

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.

Neanderthal-derived genetic variation in living humans relates to schizophrenia diagnosis, to psychotic symptom severity, and to dopamine synthesis

Schizophrenia has been hypothesized to be a human-specific condition, but experimental approaches to testing this idea have been limited. Because Neanderthals, our closest evolutionary relatives, interbred with modern humans prior to their disappearance from the fossil record, leaving a residual echo that survives in our DNA today, we leveraged new discoveries about ancient hominid DNA to explore this hypothesis in living people in three converging ways. First, in four independent case-control datasets totaling 9,362 individuals, individuals with schizophrenia had less Neanderthal-derived genetic variation than controls (p = .044). Second, in 49 unmedicated inpatients with schizophrenia, having more Neanderthal admixture predicted less severe positive symptoms (p = .046). Finally, using 18 F-fluorodopa PET scanning in 172 healthy individuals, having greater Neanderthal introgression was significantly associated with lower dopamine synthesis capacity in the striatum and pons (p's < 2 × 10-5 ), which is fundamentally important in the pathophysiology and treatment of psychosis. These results may help to elucidate the evolutionary history of a devastating neuropsychiatric disease by supporting the notion of schizophrenia as a human-specific condition. Additionally, the relationship between Neanderthal admixture and dopamine function suggests a potential mechanism whereby Neanderthal admixture may have affected our gene pool to alter schizophrenia risk and/or course.

The history and evolution of the Denisovan-EPAS1 haplotype in Tibetans

Recent studies suggest that admixture with archaic hominins played an important role in facilitating biological adaptations to new environments. For example, interbreeding with Denisovans facilitated the adaptation to high-altitude environments on the Tibetan Plateau. Specifically, the EPAS1 gene, a transcription factor that regulates the response to hypoxia, exhibits strong signatures of both positive selection and introgression from Denisovans in Tibetan individuals. Interestingly, despite being geographically closer to the Denisova Cave, East Asian populations do not harbor as much Denisovan ancestry as populations from Melanesia. Recently, two studies have suggested two independent waves of Denisovan admixture into East Asians, one of which is shared with South Asians and Oceanians. Here, we leverage data from EPAS1 in 78 Tibetan individuals to interrogate which of these two introgression events introduced the EPAS1 beneficial sequence into the ancestral population of Tibetans, and we use the distribution of introgressed segment lengths at this locus to infer the timing of the introgression and selection event. We find that the introgression event unique to East Asians most likely introduced the beneficial haplotype into the ancestral population of Tibetans around 48,700 (16,000-59,500) y ago, and selection started around 9,000 (2,500-42,000) y ago. Our estimates suggest that one of the most convincing examples of adaptive introgression is in fact selection acting on standing archaic variation.

Inferring archaic introgression from hominin genetic data

Questions surrounding the timing, extent, and evolutionary consequences of archaic admixture into human populations have a long history in evolutionary anthropology. More recently, advances in human genetics, particularly in the field of ancient DNA, have shed new light on the question of whether or not Homo sapiens interbred with other hominin groups. By the late 1990s, published genetic work had largely concluded that archaic groups made no lasting genetic contribution to modern humans; less than a decade later, this conclusion was reversed following the successful DNA sequencing of an ancient Neanderthal. This reversal of consensus is noteworthy, but the reasoning behind it is not widely understood across all academic communities. There remains a communication gap between population geneticists and paleoanthropologists. In this review, we endeavor to bridge this gap by outlining how technological advancements, new statistical methods, and notable controversies ultimately led to the current consensus.

Multi-omic analysis elucidates the genetic basis of hydrocephalus

We conducted PrediXcan analysis of hydrocephalus risk in ten neurological tissues and whole blood. Decreased expression of MAEL in the brain was significantly associated (Bonferroni-adjusted p < 0.05) with hydrocephalus. PrediXcan analysis of brain imaging and genomics data in the independent UK Biobank (N = 8,428) revealed that MAEL expression in the frontal cortex is associated with white matter and total brain volumes. Among the top differentially expressed genes in brain, we observed a significant enrichment for gene-level associations with these structural phenotypes, suggesting an effect on disease risk through regulation of brain structure and integrity. We found additional support for these genes through analysis of the choroid plexus transcriptome of a murine model of hydrocephalus. Finally, differential protein expression analysis in patient cerebrospinal fluid recapitulated disease-associated expression changes in neurological tissues, but not in whole blood. Our findings provide convergent evidence highlighting the importance of tissue-specific pathways and mechanisms in the pathophysiology of hydrocephalus.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

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.

Genomic insights into population history and biological adaptation in Oceania

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes1. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

Genome-wide ancestry and introgression in a Zambian baboon hybrid zone

Hybridization in nature offers unique insights into the process of natural selection in incipient species and their hybrids. In order to evaluate the patterns and targets of selection, we examine a recently discovered baboon hybrid zone in the Kafue River Valley of Zambia, where Kinda baboons (Papio kindae) and grey-footed chacma baboons (P. ursinus griseipes) coexist with hybridization. We genotyped baboons at 14,962 variable genome-wide autosomal markers using double-digest RADseq. We compared ancestry patterns from this genome-wide data set to previously reported ancestry from mitochondrial-DNA and Y-chromosome sources. We also fit a Bayesian genomic cline model to scan for genes with extreme patterns of introgression. We show that the Kinda baboon Y chromosome has penetrated the species boundary to a greater extent than either mitochondrial DNA or the autosomal chromosomes. We also find evidence for overall restricted introgression in the JAK/STAT signalling pathway. Echoing results in other species including humans, we find evidence for enhanced and/or directional introgression of immune-related genes or pathways including the toll-like receptor pathway, the blood coagulation pathway, and the LY96 gene. Finally we show enhanced introgression and excess chacma baboon ancestry in the sperm tail gene ODF2. Together, our results elucidate the dynamics of introgressive hybridization in a primate system while identifying genes and pathways possibly under selection.

The Evolutionary History of Common Genetic Variants Influencing Human Cortical Surface Area

Structural brain changes along the lineage leading to modern Homo sapiens contributed to our distinctive cognitive and social abilities. However, the evolutionarily relevant molecular variants impacting key aspects of neuroanatomy are largely unknown. Here, we integrate evolutionary annotations of the genome at diverse timescales with common variant associations from large-scale neuroimaging genetic screens. We find that alleles with evidence of recent positive polygenic selection over the past 2000-3000 years are associated with increased surface area (SA) of the entire cortex, as well as specific regions, including those involved in spoken language and visual processing. Therefore, polygenic selective pressures impact the structure of specific cortical areas even over relatively recent timescales. Moreover, common sequence variation within human gained enhancers active in the prenatal cortex is associated with postnatal global SA. We show that such variation modulates the function of a regulatory element of the developmentally relevant transcription factor HEY2 in human neural progenitor cells and is associated with structural changes in the inferior frontal cortex. These results indicate that non-coding genomic regions active during prenatal cortical development are involved in the evolution of human brain structure and identify novel regulatory elements and genes impacting modern human brain structure.

A genomic region associated with protection against severe COVID-19 is inherited from Neandertals

It was recently shown that the major genetic risk factor associated with becoming severely ill with COVID-19 when infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is inherited from Neandertals. New, larger genetic association studies now allow additional genetic risk factors to be discovered. Using data from the Genetics of Mortality in Critical Care (GenOMICC) consortium, we show that a haplotype at a region on chromosome 12 associated with requiring intensive care when infected with the virus is inherited from Neandertals. This region encodes proteins that activate enzymes that are important during infections with RNA viruses. In contrast to the previously described Neandertal haplotype that increases the risk for severe COVID-19, this Neandertal haplotype is protective against severe disease. It also differs from the risk haplotype in that it has a more moderate effect and occurs at substantial frequencies in all regions of the world outside Africa. Among ancient human genomes in western Eurasia, the frequency of the protective Neandertal haplotype may have increased between 20,000 and 10,000 y ago and again during the past 1,000 y.

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.

The Population-Specific Impact of Neandertal Introgression on Human Disease

Since the discovery of admixture between modern humans and Neandertals, multiple studies investigated the effect of Neandertal-derived DNA on human disease and nondisease phenotypes. These studies have linked Neandertal ancestry to skin- and hair-related phenotypes, immunity, neurological, and behavioral traits. However, these inferences have so far been limited to cohorts with participants of European ancestry. Here, I analyze summary statistics from 40 disease GWAS (genome-wide association study) cohorts of ∼212,000 individuals provided by the Biobank Japan Project for phenotypic effects of Neandertal DNA. I show that Neandertal DNA is associated with autoimmune diseases, prostate cancer and type 2 diabetes. Many of these disease associations are linked to population-specific Neandertal DNA, highlighting the importance of studying a wider range of ancestries to characterize the phenotypic legacy of Neandertals in people today.

Neanderthal-Derived Genetic Variation is Associated with Functional Connectivity in the Brains of Living Humans

Aim: To determine whether Neanderthal-derived genetic variation relates to functional connectivity patterns in the brains of living modern humans. Introduction: Nearly 50,000 years ago, Neanderthals interbred with ancestors of modern humans, imparting a genetic legacy that lives on today. The vestiges of this Neanderthal-derived genetic variation have been previously shown to be enriched in genes coding for neurogenesis and myelination and to alter skull shape and brain structure in living people. Materials and Methods: Using two independent cohorts totaling 553 healthy individuals, we employed multivariate distance matrix regression (MDMR) to determine whether any brain areas exhibited whole-brain functional connectivity patterns that significantly related to the degree of Neanderthal introgression. Identified clusters were then used as regions of interest in follow-up seed-based functional connectivity analyses to determine the connectivity patterns driving the relationships. Results: The MDMR analysis revealed that the percentage of Neanderthal-originating polymorphisms was significantly associated with the functional connectivity patterns of an area of the intraparietal sulcus (IPS) that was nearly identical in both cohorts. Using these IPS clusters as regions of interest in seed-based connectivity analyses, we found, again in both cohorts, that individuals with a higher proportion of Neanderthal-derived genetic variation showed increased IPS functional connectivity with visual processing regions, but decreased IPS connectivity with regions underlying social cognition. Conclusions: These findings demonstrate that the remnants of Neanderthal admixture continue to influence human brain function today, in ways that are consistent with anthropological conceptualizations of Neanderthal phenotypes, including the possibility that Neanderthals may have depended upon visual processing capabilities at the expense of social cognition, and this may have contributed to the extinction of this species through reduced cultural maintenance and inability to cope with fluctuating resources. This and other studies capitalizing on the emerging science surrounding ancient DNA provide a window through which to view an ancient lineage long past.

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.

Canalization and Robustness in Human Genetics and Disease

Canalization refers to the evolution of populations such that the number of individuals who deviate from the optimum trait, or experience disease, is minimized. In the presence of rapid cultural, environmental, or genetic change, the reverse process of decanalization may contribute to observed increases in disease prevalence. This review starts by defining relevant concepts, drawing distinctions between the canalization of populations and robustness of individuals. It then considers evidence pertaining to three continuous traits and six domains of disease. In each case, existing genetic evidence for genotype-by-environment interactions is insufficient to support a strong inference of decanalization, but we argue that the advent of genome-wide polygenic risk assessment now makes an empirical evaluation of the role of canalization in preventing disease possible. Finally, the contributions of both rare and common variants to congenital abnormality and adult onset disease are considered in light of a new kerplunk model of genetic effects.

Patterns and impacts of nonvertical evolution in eukaryotes: a paradigm shift

Evolution of eukaryotic species and their genomes has been traditionally understood as a vertical process in which genetic material is transmitted from parents to offspring along a lineage, and in which genetic exchange is restricted within species boundaries. However, mounting evidence from comparative genomics indicates that this paradigm is often violated. Horizontal gene transfer and mating between diverged lineages blur species boundaries and challenge the reconstruction of evolutionary histories of species and their genomes. Nonvertical evolution might be more restricted in eukaryotes than in prokaryotes, yet it is not negligible and can be common in certain groups. Recognition of such processes brings about the need to incorporate this complexity into our models, as well as to conceptually reframe eukaryotic diversity and evolution. Here, I review the recent work from genomics studies that supports the effects of nonvertical modes of evolution including introgression, hybridization, and horizontal gene transfer in different eukaryotic groups. I then discuss emerging patterns and effects, illustrated by specific examples, that support the conclusion that nonvertical processes are often at the root of important evolutionary transitions and adaptations. I will argue that a paradigm shift is needed to naturally accommodate nonvertical processes in eukaryotic evolution.

Fostering Responsible Research on Ancient DNA

Anticipating and addressing the social implications of scientific work is a fundamental responsibility of all scientists. However, expectations for ethically sound practices can evolve over time as the implications of science come to be better understood. Contemporary researchers who work with ancient human remains, including those who conduct ancient DNA research, face precisely this challenge as it becomes clear that practices such as community engagement are needed to address the important social implications of this work. To foster and promote ethical engagement between researchers and communities, we offer five practical recommendations for ancient DNA researchers: (1) formally consult with communities; (2) address cultural and ethical considerations; (3) engage communities and support capacity building; (4) develop plans to report results and manage data; and (5) develop plans for long-term responsibility and stewardship. Ultimately, every member of a research team has an important role in fostering ethical research on ancient DNA.

Secondary contact and genomic admixture between rhesus and long-tailed macaques in the Indochina Peninsula

Understanding the process and consequences of hybridization is one of the major challenges in evolutionary biology. A growing body of literature has reported evidence of ancient hybridization events or natural hybrid zones in primates, including humans; however, we still have relatively limited knowledge about the pattern and history of admixture because there have been little studies that simultaneously achieved genome-scale analysis and a geographically wide sampling of wild populations. Our study applied double-digest restriction site-associated DNA sequencing to samples from the six localities in and around the provisional hybrid zone of rhesus and long-tailed macaques and evaluated population structure, phylogenetic relationships, demographic history, and geographic clines of morphology and allele frequencies. A latitudinal gradient of genetic components was observed, highlighting the transition from rhesus (north) to long-tailed macaque distribution (south) as well as the presence of one northern population of long-tailed macaques exhibiting unique genetic structure. Interspecific gene flow was estimated to have recently occurred after an isolation period, and the migration rate from rhesus to long-tailed macaques was slightly greater than in the opposite direction. Although some rhesus macaque-biased alleles have widely introgressed into long-tailed macaque populations, the inflection points of allele frequencies have been observed as concentrated around the traditionally recognized interspecific boundary where morphology discontinuously changed; this pattern was more pronounced in the X chromosome than in autosomes. Thus, due to geographic separation before secondary contact, reproductive isolation could have evolved, contributing to the maintenance of an interspecific boundary and species-specific morphological characteristics.

Harmonizing hybridization dissonance in conservation

A dramatic increase in the hybridization between historically allopatric species has been induced by human activities. However, the notion of hybridization seems to lack consistency in two respects. On the one hand, it is inconsistent with the biological species concept, which does not allow for interbreeding between species, and on the other hand, it is considered either as an evolutionary process leading to the emergence of new biodiversity or as a cause of biodiversity loss, with conservation implications. In the first case, we argue that conservation biology should avoid the discussion around the species concept and delimit priorities of conservation units based on the impact on biodiversity if taxa are lost. In the second case, we show that this is not a paradox but an intrinsic property of hybridization, which should be considered in conservation programmes. We propose a novel view of conservation guidelines, in which human-induced hybridization may also be a tool to enhance the likelihood of adaptation to changing environmental conditions or to increase the genetic diversity of taxa affected by inbreeding depression. The conservation guidelines presented here represent a guide for the development of programmes aimed at protecting biodiversity as a dynamic evolutionary system.

Human Stem Cell Resources Are an Inroad to Neandertal DNA Functions

Induced pluripotent stem cells (iPSCs) from diverse humans offer the potential to study human functional variation in controlled culture environments. A portion of this variation originates from an ancient admixture between modern humans and Neandertals, which introduced alleles that left a phenotypic legacy on individual humans today. Here, we show that a large iPSC repository harbors extensive Neandertal DNA, including alleles that contribute to human phenotypes and diseases, encode hundreds of amino acid changes, and alter gene expression in specific tissues. We provide a database of the inferred introgressed Neandertal alleles for each individual iPSC line, together with the annotation of the predicted functional variants. We also show that transcriptomic data from organoids generated from iPSCs can be used to track Neandertal-derived RNA over developmental processes. Human iPSC resources provide an opportunity to experimentally explore Neandertal DNA function and its contribution to present-day phenotypes, and potentially study Neandertal traits.

A view of the Lower to Middle Paleolithic boundary from Northern France, far from the Near East?

Northern France and the Near East play and have played a central role in the debate around the Lower Paleolithic (LP) to Middle Paleolithic (MP) boundary. In the early 1990s, the renewed Saalian record for Northern France began to outline a mosaic model of the LP-to-MP transition-mainly based on Tuffreau's works. It implied the coexistence of Upper Acheulean assemblages (numerous bifaces with few standardized retouched flakes), 'Epi-Acheulean' assemblages (rare bifaces and diversified retouched flakes), and Mousterian assemblages (Levalloisian industries) during the Marine Isotope Stage (MIS) 8-6 period. Since the 2000s, the discovery of new key sites and enhanced field and laboratory methods are challenging this model. We present first a brief historical summary of previous approaches to the LP to MP boundary in Northern France. A large data set of Saalian archaeological units available from previous works has been updated and expanded to include additional sites. This allows us to demonstrate that the current Saalian record from Northern France is both rich and sparse as it is heterogeneously biased through time and space and that these biases limit the accuracy of any attempt to model the LP-to-MP transition. Nevertheless, we describe the differences between pre-MIS 9 and MIS 9 and MIS 8-6 records for lithic industries and discuss whether the current periodization is still relevant considering new data on technological, behavioral, and cultural changes. The comparison between Northern France and Near Eastern records allows regional cultural patterns to be distinguished from global trends in lithic trajectories of change and determination of how they slotted together. Our review of the available data from these two distant regions confirms that the LP-to-MP transition is probably one of the major cultural shifts in human evolution.

Methods for detecting introgressed archaic sequences

Analysis of genome sequences from archaic and modern humans have revealed multiple episodes of admixture between highly-diverged population groups. Statistical methods that attempt to localize DNA segments introduced by these events offer a powerful tool to investigate recent human evolution. We review recent advances in methods for detecting introgressed sequences.

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.

Studying human and nonhuman primate evolutionary biology with powerful in vitro and in vivo functional genomics tools

In recent years, tools for functional genomic studies have become increasingly feasible for use by evolutionary anthropologists. In this review, we provide brief overviews of several exciting in vitro techniques that can be paired with "-omics" approaches (e.g., genomics, epigenomics, transcriptomics, proteomics, and metabolomics) for potentially powerful evolutionary insights. These in vitro techniques include ancestral protein resurrection, cell line experiments using primary, immortalized, and induced pluripotent stem cells, and CRISPR-Cas9 genetic manipulation. We also discuss how several of these methods can be used in vivo, for transgenic organism studies of human and nonhuman primate evolution. Throughout this review, we highlight example studies in which these approaches have already been used to inform our understanding of the evolutionary biology of modern and archaic humans and other primates while simultaneously identifying future opportunities for anthropologists to use this toolkit to help answer additional outstanding questions in evolutionary anthropology.

What have the revelations about Neanderthal DNA revealed about Homo sapiens?

Genetic studies have presented increasing indications about the complexity of the interactions between Homo sapiens, Neanderthals and Denisovans, during Pleistocene. The results indicate potential replacement or admixture of the groups of hominins that lived in the same region at different times. Recently, the time of separation among these hominins in relation to the Last Common Ancestor – LCA has been reasonably well established. Events of mixing with emphasis on the Neanderthal gene flow into H. sapiens outside Africa, Denisovans into H. sapiens ancestors in Oceania and continental Asia, Neanderthals into Denisovans, as well as the origin of some phenotypic features in specific populations such as the color of the skin, eyes, hair and predisposition to develop certain kinds of diseases have also been found. The current information supports the existence of both replacement and interbreeding events, and indicates the need to revise the two main explanatory models, the Multiregional and the Out-of-Africa hypotheses, about the origin and evolution of H. sapiens and its co-relatives. There is definitely no longer the possibility of justifying only one model over the other. This paper aims to provide a brief review and update on the debate around this issue, considering the advances brought about by the recent genetic as well as morphological traits analyses.

Identifying and Interpreting Apparent Neanderthal Ancestry in African Individuals

Admixture has played a prominent role in shaping patterns of human genomic variation, including gene flow with now-extinct hominins like Neanderthals and Denisovans. Here, we describe a novel probabilistic method called IBDmix to identify introgressed hominin sequences, which, unlike existing approaches, does not use a modern reference population. We applied IBDmix to 2,504 individuals from geographically diverse populations to identify and analyze Neanderthal sequences segregating in modern humans. Strikingly, we find that African individuals carry a stronger signal of Neanderthal ancestry than previously thought. We show that this can be explained by genuine Neanderthal ancestry due to migrations back to Africa, predominately from ancestral Europeans, and gene flow into Neanderthals from an early dispersing group of humans out of Africa. Our results refine our understanding of Neanderthal ancestry in African and non-African populations and demonstrate that remnants of Neanderthal genomes survive in every modern human population studied to date.

Drosophila melanogaster Responses against Entomopathogenic Nematodes: Focus on Hemolymph Clots

Several insect innate immune mechanisms are activated in response to infection by entomopathogenic nematodes (EPNs). In this review, we focus on the coagulation of hemolymph, which acts to stop bleeding after injury and prevent access of pathogens to the body cavity. After providing a general overview of invertebrate coagulation systems, we discuss recent findings in Drosophila melanogaster which demonstrate that clots protect against EPN infections. Detailed analysis at the cellular level provided insight into the kinetics of the secretion of Drosophila coagulation factors, including non-classical modes of secretion. Roughly, clot formation can be divided into a primary phase in which crosslinking of clot components depends on the activity of Drosophila transglutaminase and a secondary, phenoloxidase (PO)-dependent phase, characterized by further hardening and melanization of the clot matrix. These two phases appear to play distinct roles in two commonly used EPN infection models, namely Heterorhabditis bacteriophora and Steinernema carpocapsae. Finally, we discuss the implications of the coevolution between parasites such as EPNs and their hosts for the dynamics of coagulation factor evolution.

Polygenic Patterns of Adaptive Introgression in Modern Humans Are Mainly Shaped by Response to Pathogens

Anatomically modern humans carry many introgressed variants from other hominins in their genomes. Some of them affect their phenotype and can thus be negatively or positively selected. Several individual genes have been proposed to be the subject of adaptive introgression, but the possibility of polygenic adaptive introgression has not been extensively investigated yet. In this study, we analyze archaic introgression maps with refined functional enrichment methods to find signals of polygenic adaptation of introgressed variants. We first apply a method to detect sets of connected genes (subnetworks) within biological pathways that present higher-than-expected levels of archaic introgression. We then introduce and apply a new statistical test to distinguish between epistatic and independent selection in gene sets of present-day humans. We identify several known targets of adaptive introgression, and we show that they belong to larger networks of introgressed genes. After correction for genetic linkage, we find that signals of polygenic adaptation are mostly explained by independent and potentially sequential selection episodes. However, we also find some gene sets where introgressed variants present significant signals of epistatic selection. Our results confirm that archaic introgression has facilitated local adaptation, especially in immunity related and metabolic functions and highlight its involvement in a coordinated response to pathogens out of Africa.

Is population structure in the genetic biobank era irrelevant, a challenge, or an opportunity?

Replicable genetic association signals have consistently been found through genome-wide association studies in recent years. The recent dramatic expansion of study sizes improves power of estimation of effect sizes, genomic prediction, causal inference, and polygenic selection, but it simultaneously increases susceptibility of these methods to bias due to subtle population structure. Standard methods using genetic principal components to correct for structure might not always be appropriate and we use a simulation study to illustrate when correction might be ineffective for avoiding biases. New methods such as trans-ethnic modeling and chromosome painting allow for a richer understanding of the relationship between traits and population structure. We illustrate the arguments using real examples (stroke and educational attainment) and provide a more nuanced understanding of population structure, which is set to be revisited as a critical aspect of future analyses in genetic epidemiology. We also make simple recommendations for how problems can be avoided in the future. Our results have particular importance for the implementation of GWAS meta-analysis, for prediction of traits, and for causal inference.

Refinement of evolutionary medicine predictions based on clinical evidence for the manifestations of Mendelian diseases

Prediction methods have become an integral part of biomedical and biotechnological research. However, their clinical interpretations are largely based on biochemical or molecular data, but not clinical data. Here, we focus on improving the reliability and clinical applicability of prediction algorithms. We assembled and curated two large non-overlapping large databases of clinical phenotypes. These phenotypes were caused by missense variations in 44 and 63 genes associated with Mendelian diseases. We used these databases to establish and validate the model, allowing us to improve the predictions obtained from EVmutation, SNAP2 and PoPMuSiC 2.1. The predictions of clinical effects suffered from a lack of specificity, which appears to be the common constraint of all recently used prediction methods, although predictions mediated by these methods are associated with nearly absolute sensitivity. We introduced evidence-based tailoring of the default settings of the prediction methods; this tailoring substantially improved the prediction outcomes. Additionally, the comparisons of the clinically observed and theoretical variations led to the identification of large previously unreported pools of variations that were under negative selection during molecular evolution. The evolutionary variation analysis approach described here is the first to enable the highly specific identification of likely disease-causing missense variations that have not yet been associated with any clinical phenotype.

Occipital hemi-bun development and shape covariation in a longitudinal extant human growth sample

OBJECTIVES

Although the homology of the Neanderthal occipital bun and anatomically modern human "hemi-bun" has long been debated, little is known about the developmental timing and patterning of these two patterns of prominent occipital squama convexity. In this study, occipital hemi-bun ontogeny and cranial shape covariation are assessed in a comparative extant human sample.

MATERIALS AND METHODS

Two-dimensional geometric morphometric methods were used to investigate hemi-bun development in a longitudinal sample of growth study cephalograms representing extant human subjects predominantly of European ancestry. Subjects were each measured at three distinct age points, ranging from 3.0 to 20.4 years, and two-block partial least squares analysis was used to assess patterns of covariation between midsagittal occipital bone morphology and other aspects of craniofacial shape.

RESULTS

Occipital hemi-bun morphology, when present, was found to develop early in ontogeny, in association with anteroposterior elongation of the frontal and parietal bones. No significant pattern of covariation was found between occipital hemi-bun shape and cranial/basicranial breadth, basicranial length, basicranial angle, or midfacial prognathism.

DISCUSSION

This study suggests that the occipital hemi-bun, at least in this extant human population, should not be considered an independent trait, as its development is closely linked to shape variation in the frontal and parietal bones. Importantly, these results suggest that occipital hemi-bun morphology is not significantly influenced by basicranial morphology during development, but instead covaries with changes in midsagittal neurocranial vault shape.

Identification of African-Specific Admixture between Modern and Archaic Humans

Recent work has demonstrated that two archaic human groups (Neanderthals and Denisovans) interbred with modern humans and contributed to the contemporary human gene pool. These findings relied on the availability of high-coverage genomes from both Neanderthals and Denisovans. Here we search for evidence of archaic admixture from a worldwide panel of 1,667 individuals using an approach that does not require the presence of an archaic human reference genome. We find no evidence for archaic admixture in the Andaman Islands, as previously claimed, or on the island of Flores, where Homo floresiensis fossils have been found. However, we do find evidence for at least one archaic admixture event in sub-Saharan Africa, with the strongest signal in Khoesan and Pygmy individuals from Southern and Central Africa. The locations of these putative archaic admixture tracts are weighted against functional regions of the genome, consistent with the long-term effects of purifying selection against introgressed genetic material.

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.

Mapping ICD-10 and ICD-10-CM Codes to Phecodes: Workflow Development and Initial Evaluation

BACKGROUND

The phecode system was built upon the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) for phenome-wide association studies (PheWAS) using the electronic health record (EHR).

OBJECTIVE

The goal of this paper was to develop and perform an initial evaluation of maps from the International Classification of Diseases, 10th Revision (ICD-10) and the International Classification of Diseases, 10th Revision, Clinical Modification (ICD-10-CM) codes to phecodes.

METHODS

We mapped ICD-10 and ICD-10-CM codes to phecodes using a number of methods and resources, such as concept relationships and explicit mappings from the Centers for Medicare & Medicaid Services, the Unified Medical Language System, Observational Health Data Sciences and Informatics, Systematized Nomenclature of Medicine-Clinical Terms, and the National Library of Medicine. We assessed the coverage of the maps in two databases: Vanderbilt University Medical Center (VUMC) using ICD-10-CM and the UK Biobank (UKBB) using ICD-10. We assessed the fidelity of the ICD-10-CM map in comparison to the gold-standard ICD-9-CM phecode map by investigating phenotype reproducibility and conducting a PheWAS.

RESULTS

We mapped >75% of ICD-10 and ICD-10-CM codes to phecodes. Of the unique codes observed in the UKBB (ICD-10) and VUMC (ICD-10-CM) cohorts, >90% were mapped to phecodes. We observed 70-75% reproducibility for chronic diseases and <10% for an acute disease for phenotypes sourced from the ICD-10-CM phecode map. Using the ICD-9-CM and ICD-10-CM maps, we conducted a PheWAS with a Lipoprotein(a) genetic variant, rs10455872, which replicated two known genotype-phenotype associations with similar effect sizes: coronary atherosclerosis (ICD-9-CM: P<.001; odds ratio (OR) 1.60 [95% CI 1.43-1.80] vs ICD-10-CM: P<.001; OR 1.60 [95% CI 1.43-1.80]) and chronic ischemic heart disease (ICD-9-CM: P<.001; OR 1.56 [95% CI 1.35-1.79] vs ICD-10-CM: P<.001; OR 1.47 [95% CI 1.22-1.77]).

CONCLUSIONS

This study introduces the beta versions of ICD-10 and ICD-10-CM to phecode maps that enable researchers to leverage accumulated ICD-10 and ICD-10-CM data for PheWAS in the EHR.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

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.