Neanderthal Introgression Research Articles

The timing of human adaptation from Neanderthal introgression.

Admixture has the potential to facilitate adaptation by providing alleles that are immediately adaptive in a new environment or by simply increasing the long-term reservoir of genetic diversity for future adaptation. A growing number of cases of adaptive introgression are being identified in species across the tree of life, however the timing of selection, and therefore the importance of the different evolutionary roles of admixture, is typically unknown. Here, we investigate the spatio-temporal history of selection favoring Neanderthal-introgressed alleles in modern human populations. Using both ancient and present-day samples of modern humans, we integrate the known demographic history of populations, namely population divergence and migration, with tests for selection. We model how a sweep placed along different branches of an admixture graph acts to modify the variance and covariance in neutral allele frequencies among populations at linked loci. Using a method based on this model of allele frequencies, we study previously identified cases of adaptive Neanderthal introgression. From these, we identify cases in which Neanderthal-introgressed alleles were quickly beneficial and other cases in which they persisted at low frequency for some time. For some of the alleles that persisted at low frequency, we show that selection likely independently favored them later on in geographically separated populations. Our work highlights how admixture with ancient hominins has contributed to modern human adaptation and contextualizes observed levels of Neanderthal ancestry in present-day and ancient samples.

Bioarchaeology-related studies in the Arabian Gulf: potentialities and shortcomings

Archaeological studies provide a powerful tool to understand the prehistoric societies, especially when combined to cutting-edge morphological and molecular anthropological analyses, allowing reconstructing past population dynamics, admixture events, and socio-cultural changes. Despite the advances achieved in the last decades by archaeological studies worldwide, several regions of the World have been spared from this scientific improvement due to various reasons. The Arabian Gulf represents a unique ground to investigate, being the passageway for human migrations and one of the hypothesized areas in which Neanderthal introgression occurred. A number of archaeological sites are currently present in the Arabian Gulf and have witnessed the antiquity and the intensiveness of the human settlements in the region. Nevertheless, the archaeological and anthropological investigation in the Gulf is still in its infancy. Data collected through archaeological studies in the area have the potential to help answering adamant questions of human history from the beginning of the structuring of genetic diversity in human species to the Neolithisation process. This review aims at providing an overview of the archaeological studies in the Arabian Gulf with special focus to Qatar, highlighting potentialities and shortcomings.

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.

The Spatial Signature of Introgression After a Biological Invasion With Hybridization

The accumulation of genome-wide molecular data has emphasized the important role of hybridization in the evolution of many organisms, which may carry introgressed genomic segments resulting from past admixture events with other taxa. Despite a number of examples of hybridization occurring during biological invasions, the resulting spatial patterns of genomic introgression remain poorly understood. Preliminary simulation studies have suggested a heterogeneous spatial level of introgression for invasive taxa after range expansion. We investigated in detail the robustness of this pattern and its persistence over time for both invasive and local organisms. Using spatially explicit simulations, we explored the spatial distribution of introgression across the area of colonization of an invasive taxon hybridizing with a local taxon. The general pattern for neutral loci supported by our results is an increasing introgression of local genes into the invasive taxon with the increase in the distance from the source of the invasion and a decreasing introgression of invasive genes into the local taxon. However, we also show there is some variation in this general trend depending on the scenario investigated. Spatial heterogeneity of introgression within a given taxon is thus an expected neutral pattern in structured populations after a biological invasion with a low to moderate amount of hybridization. We further show that this pattern is consistent with published empirical observations. Using additional simulations, we argue that the spatial pattern of Neanderthal introgression in modern humans, which has been documented to be higher in Asia than in Europe, can be explained by a model of hybridization with Neanderthals in Eurasia during the range expansion of modern humans from Africa. Our results support the view that weak hybridization during range expansion may explain spatially heterogeneous introgression patterns without the need to invoke selection.

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 prior to introgression. Our simulations and variant effect predictions argue that these reintroduced alleles (RAs) are more likely to be tolerated by modern humans than 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 (MPRA) 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.ONE SENTENCE SUMMARYNeanderthal interbreeding with anatomically modern humans restored thousands of ancient alleles that were previously lost in Eurasian populations.

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 nature of Neanderthal introgression revealed by 27,566 Icelandic genomes.

Human evolutionary history is rich with the interbreeding of divergent populations. Most humans outside of Africa trace about 2% of their genomes to admixture from Neanderthals, which occurred 50-60 thousand years ago1. Here we examine the effect of this event using 14.4million putative archaic chromosome fragments that were detected in fully phased whole-genome sequences from 27,566Icelanders, corresponding to a range of 56,388-112,709 unique archaic fragments that cover 38.0-48.2% of the callable genome. On the basis of the similarity with known archaic genomes, we assign 84.5% of fragments to an Altai or Vindija Neanderthal origin and 3.3% to Denisovan origin; 12.2% of fragments are of unknown origin. We find that Icelanders have more Denisovan-like fragments than expected through incomplete lineage sorting. This is best explained by Denisovan gene flow, either into ancestors of the introgressing Neanderthals or directly into humans. A within-individual, paired comparison of archaic fragments with syntenic non-archaic fragments revealed that, although the overall rate of mutation was similar in humans and Neanderthals during the 500thousand years that their lineages were separate, there were differences in the relative frequencies of mutation types-perhaps due to different generation intervals for males and females. Finally, we assessed 271phenotypes, report 5 associations driven by variants in archaic fragments and show that the majority of previously reported associations are better explained by non-archaic variants.

An evolutionary recent IFN/IL-6/CEBP axis is linked to monocyte expansion and tuberculosis severity in humans.

Monocyte counts are increased during human tuberculosis (TB) but it has not been determined whether Mycobacterium tuberculosis (Mtb) directly regulates myeloid commitment. We demonstrated that exposure to Mtb directs primary human CD34+ cells to differentiate into monocytes/macrophages. In vitro myeloid conversion did not require type I or type II IFN signaling. In contrast, Mtb enhanced IL-6 responses by CD34+ cell cultures and IL-6R neutralization inhibited myeloid differentiation and decreased mycobacterial growth in vitro. Integrated systems biology analysis of transcriptomic, proteomic and genomic data of large data sets of healthy controls and TB patients established the existence of a myeloid IL-6/IL6R/CEBP gene module associated with disease severity. Furthermore, genetic and functional analysis revealed the IL6/IL6R/CEBP gene module has undergone recent evolutionary selection, including Neanderthal introgression and human pathogen adaptation, connected to systemic monocyte counts. These results suggest Mtb co-opts an evolutionary recent IFN-IL6-CEBP feed-forward loop, increasing myeloid differentiation linked to severe TB in humans.

Impact and Evolutionary Determinants of Neanderthal Introgression on Transcriptional and Post-Transcriptional Regulation

Archaic admixture is increasingly recognized as an important source of diversity in modern humans, and Neanderthal haplotypes cover 1%-3% of the genome of present-day Eurasians. Recent work has shown that archaic introgression has contributed to human phenotypic diversity, mostly through the regulation of gene expression. Yet the mechanisms through which archaic variants alter gene expression and the forces driving the introgression landscape at regulatory regions remain elusive. Here, we explored the impact of archaic introgression on transcriptional and post-transcriptional regulation. We focused on promoters and enhancers across 127 different tissues as well as on microRNA (miRNA)-mediated regulation. Although miRNAs themselves harbor few archaic variants, we found that some of these variants may have a strong impact on miRNA-mediated gene regulation. Enhancers were by far the regulatory elements most affected by archaic introgression: up to one-third of the tissues we tested presented significant enrichments. Specifically, we found strong enrichments of archaic variants in adipose-related tissues and primary Tcells, even after accounting for various genomic and evolutionary confounders such as recombination rate and background selection. Interestingly, we identified signatures of adaptive introgression at enhancers of some key regulators of adipogenesis, raising the interesting hypothesis of a possible adaptation of early Eurasians to colder climates. Collectively, this study sheds new light on the mechanisms through which archaic admixture has impacted gene regulation in Eurasians and, more generally, increases our understanding of the contribution of Neanderthals to the regulation of acquired immunity and adipose homeostasis in modern humans.

Defining the genetic and evolutionary architecture of alternative splicing in response to infection

Host and environmental factors contribute to variation in human immune responses, yet the genetic and evolutionary drivers of alternative splicing in response to infection remain largely uncharacterised. Leveraging 970 RNA-sequencing profiles of resting and stimulated monocytes from 200 individuals of African- and European-descent, we show that immune activation elicits a marked remodelling of the isoform repertoire, while increasing the levels of erroneous splicing. We identify 1,464 loci associated with variation in isoform usage (sQTLs), 9% of them being stimulation-specific, which are enriched in disease-related loci. Furthermore, we detect a longstanding increased plasticity of immune gene splicing, and show that positive selection and Neanderthal introgression have both contributed to diversify the splicing landscape of human populations. Together, these findings suggest that differential isoform usage has been an important substrate of innovation in the long-term evolution of immune responses and a more recent vehicle of population local adaptation.

Analyses of Neanderthal introgression suggest that Levantine and southern Arabian populations have a shared population history.

Modern humans are thought to have interbred with Neanderthals in the Near East soon after modern humans dispersed out of Africa. This introgression event likely took place in either the Levant or southern Arabia depending on the dispersal route out of Africa that was followed. In this study, we compare Neanderthal introgression in contemporary Levantine and southern Arabian populations to investigate Neanderthal introgression and to study Near Eastern population history. We analyzed genotyping data on >400,000 autosomal SNPs from seven Levantine and five southern Arabian populations and compared these data to those from populations from around the world including Neanderthal and Denisovan genomes. We used f4 and D statistics to estimate and compare levels of Neanderthal introgression between Levantine, southern Arabian, and comparative global populations. We also identified 1,581 putative Neanderthal-introgressed SNPs within our dataset and analyzed their allele frequencies as a means to compare introgression patterns in Levantine and southern Arabian genomes. We find that Levantine and southern Arabian populations have similar levels of Neanderthal introgression to each other but lower levels than other non-Africans. Furthermore, we find that introgressed SNPs have very similar allele frequencies in the Levant and southern Arabia, which indicates that Neanderthal introgression is similarly distributed in Levantine and southern Arabian genomes. We infer that the ancestors of contemporary Levantine and southern Arabian populations received Neanderthal introgression prior to separating from each other and that there has been extensive gene flow between these populations.

SU5 - IDENTIFICATION OF POLYGENIC ADAPTATION IN ATTENTION-DEFICIT/HYPERACTIVITY DISORDER USING GWAS DATA

Background During the human Diaspora out of Africa, Anatomically Modern Humans (AMH) adapted to the new environments they encountered. On top, adaptation process could have been fueled by interbreeding with homo archaic species -such as Neanderthal and Denisovan- present at that time in Eurasia. In this context of natural adaptation, Attention-Deficit/Hyperactivity Disorder (ADHD) poses an interesting evolutionary paradox: despite being associated with significant impairment, morbidity and increased rates of mortality, it is a highly prevalent phenotype. Being a highly heritable trait, several selective hypotheses have been proposed to explain the prevalence of this phenotype. In the present study, we have analysed the signatures of recent polygenic selection and the effect of archaic inbreeding in ADHD. Methods Databases: i) The effect size estimates of millions of SNPs reported at a GWAS meta-analysis on ADHD conducted by the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) and the Psychiatric Genomics Consortium (PGC) on approximately 20 K patients and 35 K controls. ii) The SDS statistic for detecting recent selection estimated in the UK10K dataset. iii) Genotypes reported at archaic European anatomically modern humans individuals (ranging from ~2 kilo years ago (kya) to ~45 kya. iv) The sequenced Altai Neanderthal genome. v) Map of Neanderthal introgression haplotypes in current humans. Analyses: The presence of signatures of recent polygenic adaptation of ADHD in UK10K samples was estimated by means of the SDS statistic. For ancestral human samples, regression analysis was conducted between the time of each sample and the average ADHD genetic score (GS) by SNP. For the Altai Neanderthal sample, average ADHD GS was compared with the one observed with ancestral human samples. The map of introgressed Neanderthal alleles was used for estimating the effect of introgressed Neanderthal alleles in ADHD in current humans. Results Our results suggest that current North European populations show signatures of recent ( Discussion Overall, our results are compatible with a thrifty-gene-like hypothesis such as the hunter-farmer theory, where some phenotypes related to ADHD (e.g. hyperactivity, impulsivity) may have been beneficial in the past as hunters but the change of lifestyle during the last 45 kya towards farmers would have provided them as detrimental in the new environment.

Simultaneous Viral Exposure and Protection from Neanderthal Introgression

In this issue, Enard and Petrov present intriguing results on the possibility of genetic traces left behind in our genomes from adaptation to past viral epidemics that may have been initiated by interaction with Neanderthal archaic hominins. The work highlights how powerful infectious agents can act as a selective force to shape our genetic makeup.

Model-based detection and analysis of introgressed Neanderthal ancestry in modern humans.

Genetic evidence has revealed that the ancestors of modern human populations outside Africa and their hominin sister groups, notably Neanderthals, exchanged genetic material in the past. The distribution of these introgressed sequence tracts along modern-day human genomes provides insight into the selective forces acting on them and the role of introgression in the evolutionary history of hominins. Studying introgression patterns on the X-chromosome is of particular interest, as sex chromosomes are thought to play a special role in speciation. Recent studies have developed methods to localize introgressed ancestries, reporting long regions that are depleted of Neanderthal introgression and enriched in genes, suggesting negative selection against the Neanderthal variants. On the other hand, enriched Neanderthal ancestry in hair- and skin-related genes suggests that some introgressed variants facilitated adaptation to new environments. Here, we present a model-based introgression detection method called dical-admix. We demonstrate its efficiency and accuracy through extensive simulations and apply it to detect tracts of Neanderthal introgression in modern human individuals from the 1000 Genomes Project. Our findings are largely concordant with previous studies, consistent with weak selection against Neanderthal ancestry. We find evidence that selection against Neanderthal ancestry was due to higher genetic load in Neanderthals resulting from small effective population size, rather than widespread Dobzhansky-Müller incompatibilities (DMIs) that could contribute to reproductive isolation. Moreover, we confirm the previously reported low level of introgression on the X-chromosome, but find little evidence that DMIs contributed to this pattern.

Variation and Functional Impact of Neanderthal Ancestry in Western Asia.

Neanderthals contributed genetic material to modern humans via multiple admixture events. Initial admixture events presumably occurred in Western Asia shortly after humans migrated out of Africa. Despite being a focal point of admixture, earlier studies indicate lower Neanderthal introgression rates in some Western Asian populations as compared with other Eurasian populations. To better understand the genome-wide and phenotypic impact of Neanderthal introgression in the region, we sequenced whole genomes of nine present-day Europeans, Africans, and the Western Asian Druze at high depth, and analyzed available whole genome data from various other populations, including 16 genomes from present-day Turkey. Our results confirmed previous observations that contemporary Western Asian populations, on an average, have lower levels of Neanderthal-introgressed DNA relative to other Eurasian populations. Modern Western Asians also show comparatively high variability in Neanderthal ancestry, which may be attributed to the complex demographic history of the region. We further replicated the previously described depletion of putatively functional sequences among Neanderthal-introgressed haplotypes. Still, we find dozens of common Neanderthal-introgressed haplotypes in the Turkish sample associated with human phenotypes, including anthropometric and metabolic traits, as well as the immune response. One of these haplotypes is unusually long and harbors variants that affect the expression of members of the CCR gene family and are associated with celiac disease. Overall, our results paint a complex first picture of the genomic impact of Neanderthal introgression in the Western Asian populations.

Distinct selective forces and Neanderthal introgression shaped genetic diversity at genes involved in neurodevelopmental disorders

In addition to high intelligence, humans evolved specialized social-cognitive skills, which are specifically affected in children with autism spectrum disorder (ASD). Genes affected in ASD represent suitable candidates to study the evolution of human social cognition. We performed an evolutionary analysis on 68 genes associated to neurodevelopmental disorders; our data indicate that genetic diversity was shaped by distinct selective forces, including natural selection and introgression from archaic hominins. We discuss the possibility that segregation distortion during spermatogenesis accounts for a subset of ASD mutations. Finally, we detected modern-human-specific alleles in DYRK1A and TCF4. These variants are located within regions that display chromatin features typical of transcriptional enhancers in several brain areas, strongly suggesting a regulatory role. These SNPs thus represent candidates for association with neurodevelopmental disorders, and await experimental validation in future studies.

The Strength of Selection against Neanderthal Introgression.

Hybridization between humans and Neanderthals has resulted in a low level of Neanderthal ancestry scattered across the genomes of many modern-day humans. After hybridization, on average, selection appears to have removed Neanderthal alleles from the human population. Quantifying the strength and causes of this selection against Neanderthal ancestry is key to understanding our relationship to Neanderthals and, more broadly, how populations remain distinct after secondary contact. Here, we develop a novel method for estimating the genome-wide average strength of selection and the density of selected sites using estimates of Neanderthal allele frequency along the genomes of modern-day humans. We confirm that East Asians had somewhat higher initial levels of Neanderthal ancestry than Europeans even after accounting for selection. We find that the bulk of purifying selection against Neanderthal ancestry is best understood as acting on many weakly deleterious alleles. We propose that the majority of these alleles were effectively neutral—and segregating at high frequency—in Neanderthals, but became selected against after entering human populations of much larger effective size. While individually of small effect, these alleles potentially imposed a heavy genetic load on the early-generation human–Neanderthal hybrids. This work suggests that differences in effective population size may play a far more important role in shaping levels of introgression than previously thought.

The Genetic Cost of Neanderthal Introgression.

Approximately 2-4% of genetic material in human populations outside Africa is derived from Neanderthals who interbred with anatomically modern humans. Recent studies have shown that this Neanderthal DNA is depleted around functional genomic regions; this has been suggested to be a consequence of harmful epistatic interactions between human and Neanderthal alleles. However, using published estimates of Neanderthal inbreeding and the distribution of mutational fitness effects, we infer that Neanderthals had at least 40% lower fitness than humans on average; this increased load predicts the reduction in Neanderthal introgression around genes without the need to invoke epistasis. We also predict a residual Neanderthal mutational load in non-Africans, leading to a fitness reduction of at least 0.5%. This effect of Neanderthal admixture has been left out of previous debate on mutation load differences between Africans and non-Africans. We also show that if many deleterious mutations are recessive, the Neanderthal admixture fraction could increase over time due to the protective effect of Neanderthal haplotypes against deleterious alleles that arose recently in the human population. This might partially explain why so many organisms retain gene flow from other species and appear to derive adaptive benefits from introgression.

Strong Selective Sweeps on the X Chromosome in the Human-Chimpanzee Ancestor Explain Its Low Divergence.

The human and chimpanzee X chromosomes are less divergent than expected based on autosomal divergence. We study incomplete lineage sorting patterns between humans, chimpanzees and gorillas to show that this low divergence can be entirely explained by megabase-sized regions comprising one-third of the X chromosome, where polymorphism in the human-chimpanzee ancestral species was severely reduced. We show that background selection can explain at most 10% of this reduction of diversity in the ancestor. Instead, we show that several strong selective sweeps in the ancestral species can explain it. We also report evidence of population specific sweeps in extant humans that overlap the regions of low diversity in the ancestral species. These regions further correspond to chromosomal sections shown to be devoid of Neanderthal introgression into modern humans. This suggests that the same X-linked regions that undergo selective sweeps are among the first to form reproductive barriers between diverging species. We hypothesize that meiotic drive is the underlying mechanism causing these two observations.

Neanderthal origin of the haplotypes carrying the functional variant Val92Met in the MC1R in modern humans.

Skin color is one of the most visible and important phenotypes of modern humans. Melanocyte-stimulating hormone and its receptor played an important role in regulating skin color. In this article, we present evidence of Neanderthal introgression encompassing the melanocyte-stimulating hormone receptor gene MC1R. The haplotypes from Neanderthal introgression diverged with the Altai Neanderthal 103.3 ka, which postdates the anatomically modern human-Neanderthal divergence. We further discovered that all of the putative Neanderthal introgressive haplotypes carry the Val92Met variant, a loss-of-function variant in MC1R that is associated with multiple dermatological traits including skin color and photoaging. Frequency of this Neanderthal introgression is low in Europeans (∼5%), moderate in continental East Asians (∼30%), and high in Taiwanese aborigines (60-70%). As the putative Neanderthal introgressive haplotypes carry a loss-of-function variant that could alter the function of MC1R and is associated with multiple traits related to skin color, we speculate that the Neanderthal introgression may have played an important role in the local adaptation of Eurasians to sunlight intensity.