The phenotypic legacy of admixture between modern humans and Neandertals

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations

Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.

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

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

Belowground cascading biotic interactions trigger crop diversity benefits

Crop diversification practices offer numerous synergistic benefits. So far, research has traditionally been confined to exploring isolated, unidirectional single-process interactions among plants, soil, and microorganisms. Here, we present a novel and systematic perspective, unveiling the intricate web of plant-soil-microbiome interactions that trigger cascading effects. Applying the principles of cascading interactions can be an alternative way to overcome soil obstacles such as soil compaction and soil pathogen pressure. Finally, we introduce a research framework comprising the design of diversified cropping systems by including commercial varieties and crops with resource-efficient traits, the exploration of cascading effects, and the innovation of field management. We propose that this provides theoretical and methodological insights that can reveal new mechanisms by which crop diversity increases productivity.

Enrichment of a subset of Neanderthal polymorphisms in autistic probands and siblings

Homo sapiens and Neanderthals underwent hybridization during the Middle/Upper Paleolithic age, culminating in retention of small amounts of Neanderthal-derived DNA in the modern human genome. In the current study, we address the potential roles Neanderthal single nucleotide polymorphisms (SNP) may be playing in autism susceptibility in samples of black non-Hispanic, white Hispanic, and white non-Hispanic people using data from the Simons Foundation Powering Autism Research (SPARK), Genotype-Tissue Expression (GTEx), and 1000 Genomes (1000G) databases. We have discovered that rare variants are significantly enriched in autistic probands compared to race-matched controls. In addition, we have identified 25 rare and common SNPs that are significantly enriched in autism on different ethnic backgrounds, some of which show significant clinical associations. We have also identified other SNPs that share more specific genotype-phenotype correlations but which are not necessarily enriched in autism and yet may nevertheless play roles in comorbid phenotype expression (e.g., intellectual disability, epilepsy, and language regression). These results strongly suggest Neanderthal-derived DNA is playing a significant role in autism susceptibility across major populations in the United States.

Neandertal ancestry through time: Insights from genomes of ancient and present-day humans

Gene flow from Neandertals has shaped the landscape of genetic and phenotypic variation in modern humans. We identify the location and size of introgressed Neandertal ancestry segments in more than 300 genomes spanning the last 50,000 years. We study how Neandertal ancestry is shared among individuals to infer the time and duration of the Neandertal gene flow. We find the correlation of Neandertal segment locations across individuals and their divergence to sequenced Neandertals, both support a model of single major Neandertal gene flow. Our catalog of introgressed segments through time confirms that most natural selection-positive and negative-on Neandertal ancestry variants occurred immediately after the gene flow, and provides new insights into how the contact with Neandertals shaped human origins and adaptation.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression

INTRODUCTION

Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear.

MATERIALS AND METHODS

Population genetics analysis and functional genomics work was performed for this locus.

RESULTS

Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human.

CONCLUSION

This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.

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

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

More than a decade of genetic research on the Denisovans

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

Towards correlative archaeology of the human mind

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

Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes

BACKGROUND

Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs.

METHODS

We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis.

RESULTS

Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns.

CONCLUSION

Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Infectious disease in the Pleistocene: Old friends or old foes?

The impact of endemic and epidemic disease on humans has traditionally been seen as a comparatively recent historical phenomenon associated with the Neolithisation of human groups, an increase in population size led by sedentarism, and increasing contact with domesticated animals as well as species occupying opportunistic symbiotic and ectosymbiotic relationships with humans. The orthodox approach is that Neolithisation created the conditions for increasing population size able to support a reservoir of infectious disease sufficient to act as selective pressure. This orthodoxy is the result of an overly simplistic reliance on skeletal data assuming that no skeletal lesions equated to a healthy individual, underpinned by the assumption that hunter-gatherer groups were inherently healthy while agricultural groups acted as infectious disease reservoirs. The work of van Blerkom, Am. J. Phys. Anthropol., vol. suppl 37 (2003), Wolfe et al., Nature, vol. 447 (2007) and Houldcroft and Underdown, Am. J. Phys. Anthropol., vol. 160, (2016) has changed this landscape by arguing that humans and pathogens have long been fellow travelers. The package of infectious diseases experienced by our ancient ancestors may not be as dissimilar to modern infectious diseases as was once believed. The importance of DNA, from ancient and modern sources, to the study of the antiquity of infectious disease, and its role as a selective pressure cannot be overstated. Here we consider evidence of ancient epidemic and endemic infectious diseases with inferences from modern and ancient human and hominin DNA, and from circulating and extinct pathogen genomes. We argue that the pandemics of the past are a vital tool to unlock the weapons needed to fight pandemics of the future.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Exploring the Neandertal legacy of pancreatic ductal adenocarcinoma risk in Eurasians

BACKGROUND

The genomes of present-day non-Africans are composed of 1-3% of Neandertal-derived DNA as a consequence of admixture events between Neandertals and anatomically modern humans about 50-60 thousand years ago. Neandertal-introgressed single nucleotide polymorphisms (aSNPs) have been associated with modern human disease-related traits, which are risk factors for pancreatic ductal adenocarcinoma (PDAC), such as obesity, type 2 diabetes, and inflammation. In this study, we aimed at investigating the role of aSNPs in PDAC in three Eurasian populations.

RESULTS

The high-coverage Vindija Neandertal genome was used to select aSNPs in non-African populations from 1000 Genomes project phase 3 data. Then, the association between aSNPs and PDAC risk was tested independently in Europeans and East Asians, using existing GWAS data on more than 200 000 individuals. We did not find any significant associations between aSNPs and PDAC in samples of European descent, whereas, in East Asians, we observed that the Chr10p12.1-rs117585753-T allele (MAF = 10%) increased the risk to develop PDAC (OR = 1.35, 95%CI 1.19-1.54, P = 3.59 × 10-6), with a P-value close to a threshold that takes into account multiple testing.

CONCLUSIONS

Our results show only a minimal contribution of Neandertal SNPs to PDAC risk.

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

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

The contribution of Neanderthal introgression and natural selection to neurodegenerative diseases

Humans are thought to be more susceptible to neurodegeneration than equivalently-aged primates. It is not known whether this vulnerability is specific to anatomically-modern humans or shared with other hominids. The contribution of introgressed Neanderthal DNA to neurodegenerative disorders remains uncertain. It is also unclear how common variants associated with neurodegenerative disease risk are maintained by natural selection in the population despite their deleterious effects. In this study, we aimed to quantify the genome-wide contribution of Neanderthal introgression and positive selection to the heritability of complex neurodegenerative disorders to address these questions. We used stratified-linkage disequilibrium score regression to investigate the relationship between five SNP-based signatures of natural selection, reflecting different timepoints of evolution, and genome-wide associated variants of the three most prevalent neurodegenerative disorders: Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease. We found no evidence for enrichment of positively-selected SNPs in the heritability of Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease, suggesting that common deleterious disease variants are unlikely to be maintained by positive selection. There was no enrichment of Neanderthal introgression in the SNP-heritability of these disorders, suggesting that Neanderthal admixture is unlikely to have contributed to disease risk. These findings provide insight into the origins of neurodegenerative disorders within the evolution of Homo sapiens and addresses a long-standing debate, showing that Neanderthal admixture is unlikely to have contributed to common genetic risk of neurodegeneration in anatomically-modern humans.

The lingering effects of Neanderthal introgression on human complex traits

The genetic variants introduced into the ancestors of modern humans from interbreeding with Neanderthals have been suggested to contribute an unexpected extent to complex human traits. However, testing this hypothesis has been challenging due to the idiosyncratic population genetic properties of introgressed variants. We developed rigorous methods to assess the contribution of introgressed Neanderthal variants to heritable trait variation and applied these methods to analyze 235,592 introgressed Neanderthal variants and 96 distinct phenotypes measured in about 300,000 unrelated white British individuals in the UK Biobank. Introgressed Neanderthal variants make a significant contribution to trait variation (explaining 0.12% of trait variation on average). However, the contribution of introgressed variants tends to be significantly depleted relative to modern human variants matched for allele frequency and linkage disequilibrium (about 59% depletion on average), consistent with purifying selection on introgressed variants. Different from previous studies (McArthur et al., 2021), we find no evidence for elevated heritability across the phenotypes examined. We identified 348 independent significant associations of introgressed Neanderthal variants with 64 phenotypes. Previous work (Skov et al., 2020) has suggested that a majority of such associations are likely driven by statistical association with nearby modern human variants that are the true causal variants. Applying a customized fine-mapping led us to identify 112 regions across 47 phenotypes containing 4303 unique genetic variants where introgressed variants are highly likely to have a phenotypic effect. Examination of these variants reveals their substantial impact on genes that are important for the immune system, development, and metabolism.

Pathogenic Variants Associated with Rare Monogenic Diseases Established in Ancient Neanderthal and Denisovan Genome-Wide Data

Ancient anatomically modern humans (AMHs) encountered other archaic human species, most notably Neanderthals and Denisovans, when they left Africa and spread across Europe and Asia ~60,000 years ago. They interbred with them, and modern human genomes retain DNA inherited from these interbreeding events. High quality (high coverage) ancient human genomes have recently been sequenced allowing for a direct estimation of individual heterozygosity, which has shown that genetic diversity in these archaic human groups was very low, indicating low population sizes. In this study, we analyze ten ancient human genome-wide data, including four sequenced with high-coverage. We screened these ancient genome-wide data for pathogenic mutations associated with monogenic diseases, and established unusual aggregation of pathogenic mutations in individual subjects, including quadruple homozygous cases of pathogenic variants in the PAH gene associated with the condition phenylketonuria in a ~120,000 years old Neanderthal. Such aggregation of pathogenic mutations is extremely rare in contemporary populations, and their existence in ancient humans could be explained by less significant clinical manifestations coupled with small community sizes, leading to higher inbreeding levels. Our results suggest that pathogenic variants associated with rare diseases might be the result of introgression from other archaic human species, and archaic admixture thus could have influenced disease risk in modern humans.

Long-range regulatory effects of Neandertal DNA in modern humans

The admixture between modern humans and Neandertals has resulted in ∼2% of the genomes of present-day non-Africans being composed of Neandertal DNA. Introgressed Neandertal DNA has been demonstrated to significantly affect the transcriptomic landscape in people today and via this molecular mechanism influence phenotype variation as well. However, little is known about how much of that regulatory impact is mediated through long-range regulatory effects that have been shown to explain ∼20% of expression variation. Here we identified 60 transcription factors (TFs) with their top cis-eQTL SNP in GTEx being of Neandertal ancestry and predicted long-range Neandertal DNA-induced regulatory effects by screening for the predicted target genes of those TFs. We show that the TFs form a significantly connected protein-protein interaction network. Among them are JUN and PRDM5, two brain-expressed TFs that have their predicted target genes enriched in regions devoid of Neandertal DNA. Archaic cis-eQTLs for the 60 TFs include multiple candidates for local adaptation, some of which show significant allele frequency increases over the last ∼10,000 years. A large proportion of the cis-eQTL-associated archaic SNPs have additional associations with various immune traits, schizophrenia, blood cell type composition and anthropometric measures. Finally, we demonstrate that our results are consistent with those of Neandertal DNA-associated empirical trans-eQTLs. Our results suggest that Neandertal DNA significantly influences regulatory networks, that its regulatory reach goes beyond the 40% of genomic sequence it still covers in present-day non-Africans and that via the investigated mechanism Neandertal DNA influences the phenotypic variation in people today.

Evolutionary Origin of Germline Pathogenic MUTYH Variations in Modern Humans

MUTYH plays an essential role in preventing oxidation-caused DNA damage. Pathogenic germline variations in MUTYH damage its function, causing intestinal polyposis and colorectal cancer. Determination of the evolutionary origin of the variation is essential to understanding the etiological relationship between MUTYH variation and cancer development. In this study, we analyzed the origins of pathogenic germline variants in human MUTYH. Using a phylogenic approach, we searched pathogenic MUTYH variants in modern humans in the MUTYH of 99 vertebrates across eight clades. We did not find pathogenic variants shared between modern humans and the non-human vertebrates following the evolutionary tree, ruling out the possibility of cross-species conservation as the origin of human pathogenic variants in MUTYH. We then searched the variants in the MUTYH of 5031 ancient humans and extinct Neanderthals and Denisovans. We identified 24 pathogenic variants in 42 ancient humans dated between 30,570 and 480 years before present (BP), and three pathogenic variants in Neanderthals dated between 65,000 and 38,310 years BP. Data from our study revealed that human pathogenic MUTYH variants mostly arose in recent human history and were partially inherited from Neanderthals.

MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations

Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.

Genetic variations within human gained enhancer elements affect human brain sulcal morphology

The expansion of the cerebral cortex is one of the most distinctive changes in the evolution of the human brain. Cortical expansion and related increases in cortical folding may have contributed to emergence of our capacities for high-order cognitive abilities. Molecular analysis of humans, archaic hominins, and non-human primates has allowed identification of chromosomal regions showing evolutionary changes at different points of our phylogenetic history. In this study, we assessed the contributions of genomic annotations spanning 30 million years to human sulcal morphology measured via MRI in more than 18,000 participants from the UK Biobank. We found that variation within brain-expressed human gained enhancers, regulatory genetic elements that emerged since our last common ancestor with Old World monkeys, explained more trait heritability than expected for the left and right calloso-marginal posterior fissures and the right central sulcus. Intriguingly, these are sulci that have been previously linked to the evolution of locomotion in primates and later on bipedalism in our hominin ancestors.

Moving beyond the adaptationist paradigm for human evolution, and why it matters

The Journal of Human Evolution (JHE) was founded 50 years ago when much of the foundation for how we think about human evolution was in place or being put in place, providing the main framework for how we consider our origins today. Here, we will explore historical developments, including early JHE outputs, as they relate to our understanding of the relationship between phenotypic variation and evolutionary process, and use that as a springboard for considering our current understanding of these links as applied to human evolution. We will focus specifically on how the study of variation itself has shifted us away from taxonomic and adaptationist perspectives toward a richer understanding of the processes shaping human evolutionary history, using literature searches and specific test cases to highlight this. We argue that natural selection, gene exchange, genetic drift, and mutation should not be considered individually when considering the production of hominin diversity. In this context, we offer suggestions for future research directions and reflect on this more complex understanding of human evolution and its broader relevance to society. Finally, we end by considering authorship demographics and practices in the last 50 years within JHE and how a shift in these demographics has the potential to reshape the science of human evolution going forward.

The role of Neanderthal introgression in liver cancer

BACKGROUND

Neanderthal introgressed DNA has been linked to different normal and disease traits including immunity and metabolism-two important functions that are altered in liver cancer. However, there is limited understanding of the relationship between Neanderthal introgression and liver cancer risk. The aim of this study was to investigate the relationship between Neanderthal introgression and liver cancer risk.

METHODS

Using germline and somatic DNA and tumor RNA from liver cancer patients from The Cancer Genome Atlas, along with ancestry-match germline DNA from unaffected individuals from the 1000 Genomes Resource, and allele specific expression data from normal liver tissue from The Genotype-Tissue Expression project we investigated whether Neanderthal introgression impacts cancer etiology. Using a previously generated set of Neanderthal alleles, we identified Neanderthal introgressed haplotypes. We then tested whether somatic mutations are enriched or depleted on Neanderthal introgressed haplotypes compared to modern haplotypes. We also computationally assessed whether somatic mutations have a functional effect or show evidence of regulating expression of Neanderthal haplotypes. Finally, we compared patterns of Neanderthal introgression in liver cancer patients and the general population.

RESULTS

We find Neanderthal introgressed haplotypes exhibit an excess of somatic mutations compared to modern haplotypes. Variant Effect Predictor analysis revealed that most of the somatic mutations on these Neanderthal introgressed haplotypes are not functional. We did observe expression differences of Neanderthal alleles between tumor and normal for four genes that also showed a pattern of enrichment of somatic mutations on Neanderthal haplotypes. However, gene expression was similar between liver cancer patients with modern ancestry and liver cancer patients with Neanderthal ancestry at these genes. Provocatively, when analyzing all genes, we find evidence of Neanderthal introgression regulating expression in tumor from liver cancer patients in two genes, ARK1C4 and OAS1. Finally, we find that most genes do not show a difference in the proportion of Neanderthal introgression between liver cancer patients and the general population.

CONCLUSION

Our results suggest that Neanderthal introgression provides opportunity for somatic mutations to accumulate, and that some Neanderthal introgression may impact liver cancer risk.

Denisovan introgression has shaped the immune system of present-day Papuans

Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.

Denisovan and Neanderthal archaic introgression differentially impacted the genetics of complex traits in modern populations

BACKGROUND

Introgression from extinct Neanderthal and Denisovan human species has been shown to contribute to the genetic pool of modern human populations and their phenotypic spectrum. Evidence of how Neanderthal introgression shaped the genetics of human traits and diseases has been extensively studied in populations of European descent, with signatures of admixture reported for instance in genes associated with pigmentation, immunity, and metabolic traits. However, limited information is currently available about the impact of archaic introgression on other ancestry groups. Additionally, to date, no study has been conducted with respect to the impact of Denisovan introgression on the health and disease of modern populations. Here, we compare the way evolutionary pressures shaped the genetics of complex traits in East Asian and European populations, and provide evidence of the impact of Denisovan introgression on the health of East Asian and Central/South Asian populations.

RESULTS

Leveraging genome-wide association statistics from the Biobank Japan and UK Biobank, we assessed whether Denisovan and Neanderthal introgression together with other evolutionary genomic signatures were enriched for the heritability of physiological and pathological conditions in populations of East Asian and European descent. In EAS, Denisovan-introgressed loci were enriched for coronary artery disease heritability (1.69-fold enrichment, p=0.003). No enrichment for archaic introgression was observed in EUR. We also performed a phenome-wide association study of Denisovan and Neanderthal alleles in six ancestry groups available in the UK Biobank. In EAS, the Denisovan-introgressed SNP rs62391664 in the major histocompatibility complex region was associated with albumin/globulin ratio (beta=-0.17, p=3.57×10^-7). Neanderthal-introgressed alleles were associated with psychiatric and cognitive traits in EAS (e.g., "No Bipolar or Depression"-rs79043717 beta=-1.5, p=1.1×10^-7), and with blood biomarkers (e.g., alkaline phosphatase-rs11244089 beta=0.1, p=3.69×10^-116) and red hair color (rs60733936 beta=-0.86, p=4.49×10^-165) in EUR. In the other ancestry groups, Neanderthal alleles were associated with several traits, also including the use of certain medications (e.g., Central/South East Asia: indapamide - rs732632 beta=-2.38, p=5.22×10^-7).

CONCLUSIONS

Our study provides novel evidence regarding the impact of archaic introgression on the genetics of complex traits in worldwide populations, highlighting the specific contribution of Denisovan introgression in EAS populations.

Discovery of 42 genome-wide significant loci associated with dyslexia

Reading and writing are crucial life skills but roughly one in ten children are affected by dyslexia, which can persist into adulthood. Family studies of dyslexia suggest heritability up to 70%, yet few convincing genetic markers have been found. Here we performed a genome-wide association study of 51,800 adults self-reporting a dyslexia diagnosis and 1,087,070 controls and identified 42 independent genome-wide significant loci: 15 in genes linked to cognitive ability/educational attainment, and 27 new and potentially more specific to dyslexia. We validated 23 loci (13 new) in independent cohorts of Chinese and European ancestry. Genetic etiology of dyslexia was similar between sexes, and genetic covariance with many traits was found, including ambidexterity, but not neuroanatomical measures of language-related circuitry. Dyslexia polygenic scores explained up to 6% of variance in reading traits, and might in future contribute to earlier identification and remediation of dyslexia.

Neandertal introgression partitions the genetic landscape of neuropsychiatric disorders and associated behavioral phenotypes

Despite advances in identifying the genetic basis of psychiatric and neurological disorders, fundamental questions about their evolutionary origins remain elusive. Here, introgressed variants from archaic humans such as Neandertals can serve as an intriguing research paradigm. We compared the number of associations for Neandertal variants to the number of associations of frequency-matched non-archaic variants with regard to human CNS disorders (neurological and psychiatric), nervous system drug prescriptions (as a proxy for disease), and related, non-disease phenotypes in the UK biobank (UKBB). While no enrichment for Neandertal genetic variants were observed in the UKBB for psychiatric or neurological disease categories, we found significant associations with certain behavioral phenotypes including pain, chronotype/sleep, smoking and alcohol consumption. In some instances, the enrichment signal was driven by Neandertal variants that represented the strongest association genome-wide. SNPs within a Neandertal haplotype that was associated with smoking in the UKBB could be replicated in four independent genomics datasets.Our data suggest that evolutionary processes in recent human evolution like admixture with Neandertals significantly contribute to behavioral phenotypes but not psychiatric and neurological diseases. These findings help to link genetic variants in a population to putative past beneficial effects, which likely only indirectly contribute to pathology in modern day humans.

Merging morphological and genetic evidence to assess hybridization in Western Eurasian late Pleistocene hominins

Previous scientific consensus saw human evolution as defined by adaptive differences (behavioural and/or biological) and the emergence of Homo sapiens as the ultimate replacement of non-modern groups by a modern, adaptively more competitive group. However, recent research has shown that the process underlying our origins was considerably more complex. While archaeological and fossil evidence suggests that behavioural complexity may not be confined to the modern human lineage, recent palaeogenomic work shows that gene flow between distinct lineages (for example, Neanderthals, Denisovans, early H. sapiens) occurred repeatedly in the late Pleistocene, probably contributing elements to our genetic make-up that might have been crucial to our success as a diverse, adaptable species. Following these advances, the prevailing human origins model has shifted from one of near-complete replacement to a more nuanced view of partial replacement with considerable reticulation. Here we provide a brief introduction to the current genetic evidence for hybridization among hominins, its prevalence in, and effects on, comparative mammal groups, and especially how it manifests in the skull. We then explore the degree to which cranial variation seen in the fossil record of late Pleistocene hominins from Western Eurasia corresponds with our current genetic and comparative data. We are especially interested in understanding the degree to which skeletal data can reflect admixture. Our findings indicate some correspondence between these different lines of evidence, flag individual fossils as possibly admixed, and suggest that different cranial regions may preserve hybridization signals differentially. We urge further studies of the phenotype to expand our ability to detect the ways in which migration, interaction and genetic exchange have shaped the human past, beyond what is currently visible with the lens of ancient DNA.

The contribution of Neanderthal introgression to modern human traits

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

Selection against admixture and gene regulatory divergence in a long-term primate field study

Genetic admixture is central to primate evolution. We combined 50 years of field observations of immigration and group demography with genomic data from ~9 generations of hybrid baboons to investigate the consequences of admixture in the wild. Despite no obvious fitness costs to hybrids, we found signatures of selection against admixture similar to those described for archaic hominins. These patterns were concentrated near genes where ancestry is strongly associated with gene expression. Our analyses also show that introgression is partially predictable across the genome. This study demonstrates the value of integrating genomic and field data for revealing how "genomic signatures of selection" (e.g., reduced introgression in low-recombination regions) manifest in nature; moreover, it underscores the importance of other primates as living models for human evolution.

Immunoglobulin genes, reproductive isolation and vertebrate speciation

Reproductive isolation drives the formation of new species, and many genes contribute to this through Dobzhansky–Muller incompatibilities (DMIs). These incompatibilities occur when gene divergence affects loci encoding interacting products such as receptors and their ligands. We suggest here that the nature of vertebrate immunoglobulin (IG) genes must make them prone to DMIs. The genes of these complex loci form functional genes through the process of recombination, giving rise to a repertoire of heterodimeric receptors of incredible diversity. This repertoire, within individuals and within species, must defend against pathogens but must also avoid pathogenic self‐reactivity. We suggest that this avoidance of autoimmunity is only achieved through a coordination of evolution between heavy‐ and light‐chain genes, and between these genes and the rest of the genome. Without coordinated evolution, the hybrid offspring of two diverging populations will carry a heavy burden of DMIs, resulting in a loss of fitness. Critical incompatibilities could manifest as incompatibilities between a mother and her divergent offspring. During fetal development, biochemical differences between the parents of hybrid offspring could make their offspring a target of the maternal immune system. This hypothesis was conceived in the light of recent insights into the population genetics of IG genes. This has suggested that antibody genes are probably as susceptible to evolutionary forces as other parts of the genome. Further repertoire studies in human and nonhuman species should now help determine whether antibody genes have been part of the evolutionary forces that drive the development of species.

Legacies at work: plant-soil-microbiome interactions underpinning agricultural sustainability

Agricultural intensification has had long-lasting negative legacies largely because of excessive inputs of agrochemicals (e.g., fertilizers) and simplification of cropping systems (e.g., continuous monocropping). Conventional agricultural management focuses on suppressing these negative legacies. However, there is now increasing attention for creating positive above- and belowground legacies through selecting crop species/genotypes, optimizing temporal and spatial crop combinations, improving nutrient inputs, developing intelligent fertilizers, and applying soil or microbiome inoculations. This can lead to enhanced yields and reduced pest and disease pressure in cropping systems, and can also mitigate greenhouse gas emissions and enhance carbon sequestration in soils. Strengthening positive legacies requires a deeper understanding of plant-soil-microbiome interactions and innovative crop, input, and soil management which can help to achieve agricultural sustainability.

The clinically relevant CYP2C8*3 and CYP2C9*2 haplotype is inherited from Neandertals

Genetic variation in genes encoding cytochrome P450 enzymes influences the metabolism of drugs and endogenous compounds. The locus containing the cytochrome genes CYP2C8 and CYP2C9 on chromosome 10 exhibits linkage disequilibrium between the CYP2C8*3 and CYP2C9*2 alleles, forming a haplotype of ~300 kilobases. This haplotype is associated with altered metabolism of several drugs, most notably reduced metabolism of warfarin and phenytoin, leading to toxicity at otherwise therapeutic doses. Here we show that this haplotype is inherited from Neandertals.

Predicting Archaic Hominin Phenotypes from Genomic Data

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

Detection of Neanderthal Adaptively Introgressed Genetic Variants That Modulate Reporter Gene Expression in Human Immune Cells

Although some variation introgressed from Neanderthals has undergone selective sweeps, little is known about its functional significance. We used a Massively Parallel Reporter Assay (MPRA) to assay 5,353 high-frequency introgressed variants for their ability to modulate the gene expression within 170 bp of endogenous sequence. We identified 2,548 variants in active putative cis-regulatory elements (CREs) and 292 expression-modulating variants (emVars). These emVars are predicted to alter the binding motifs of important immune transcription factors, are enriched for associations with neutrophil and white blood cell count, and are associated with the expression of genes that function in innate immune pathways including inflammatory response and antiviral defense. We combined the MPRA data with other data sets to identify strong candidates to be driver variants of positive selection including an emVar that may contribute to protection against severe COVID-19 response. We endogenously deleted two CREs containing expression-modulation variants linked to immune function, rs11624425 and rs80317430, identifying their primary genic targets as ELMSAN1, and PAN2 and STAT2, respectively, three genes differentially expressed during influenza infection. Overall, we present the first database of experimentally identified expression-modulating Neanderthal-introgressed alleles contributing to potential immune response in modern humans.

pyPheWAS: A Phenome-Disease Association Tool for Electronic Medical Record Analysis

Along with the increasing availability of electronic medical record (EMR) data, phenome-wide association studies (PheWAS) and phenome-disease association studies (PheDAS) have become a prominent, first-line method of analysis for uncovering the secrets of EMR. Despite this recent growth, there is a lack of approachable software tools for conducting these analyses on large-scale EMR cohorts. In this article, we introduce pyPheWAS, an open-source python package for conducting PheDAS and related analyses. This toolkit includes 1) data preparation, such as cohort censoring and age-matching; 2) traditional PheDAS analysis of ICD-9 and ICD-10 billing codes; 3) PheDAS analysis applied to a novel EMR phenotype mapping: current procedural terminology (CPT) codes; and 4) novelty analysis of significant disease-phenotype associations found through PheDAS. The pyPheWAS toolkit is approachable and comprehensive, encapsulating data prep through result visualization all within a simple command-line interface. The toolkit is designed for the ever-growing scale of available EMR data, with the ability to analyze cohorts of 100,000 + patients in less than 2 h. Through a case study of Down Syndrome and other intellectual developmental disabilities, we demonstrate the ability of pyPheWAS to discover both known and potentially novel disease-phenotype associations across different experiment designs and disease groups. The software and user documentation are available in open source at https://github.com/MASILab/pyPheWAS .

Lower promoter activity of the ST8SIA2 gene has been favored in evolving human collective brains

ST8SIA2 is an important molecule regulating expression of the phenotype involved in schizophrenia. Lowered promoter activity of the ST8SIA2 gene is considered to be protective against schizophrenia by conferring tolerance to psychosocial stress. Here, we examined the promoter-type composition of anatomically modern humans (AMHs) and archaic humans (AHs; Neanderthals and Denisovans), and compared the promoter activity at the population level (population promoter activity; PPA) between them. In AMHs, the TCT-type, showing the second lowest promoter activity, was most prevalent in the ancestral population of non-Africans. However, the detection of only the CGT-type from AH samples and recombination tracts in AH sequences showed that the CGT- and TGT-types, exhibiting the two highest promoter activities, were common in AH populations. Furthermore, interspecies gene flow occurred into AMHs from AHs and into Denisovans from Neanderthals, influencing promoter-type compositions independently in both AMHs and AHs. The difference of promoter-type composition makes PPA unique in each population. East and Southeast Asian populations show the lowest PPA. This results from the selective increase of the CGC-type, showing the lowest promoter activity, in these populations. Every non-African population shows significantly lower PPA than African populations, resulting from the TCT-type having the highest prevalence in the ancestral population of non-Africans. In addition, PPA reduction is also found among subpopulations within Africa via a slight increase of the TCT-type. These findings indicate a trend toward lower PPA in the spread of AMHs, interpreted as a continuous adaptation to psychosocial stress arising in migration. This trend is considered as genetic tuning for the evolution of collective brains. The inferred promoter-type composition of AHs differed markedly from that of AMHs, resulting in higher PPA in AHs than in AMHs. This suggests that the trend toward lower PPA is a unique feature in AMH spread.

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.

We lack a comprehensive understanding of how Neanderthal ancestry influences human traits. This study finds that regions with Neanderthal ancestry are broadly depleted of trait-associated variation; yet, introgressed variants likely contributed to human adaptation in a few traits, like skin color and immune response modulation.

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.