Resurrecting surviving Neandertal lineages from modern human genomes

Exploring Mitonuclear Discordance: Ghost Introgression From an Ancient Extinction Lineage in the Odorrana swinhoana Complex

Mitonuclear discordance, the incongruence between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA), is a well-documented phenomenon with various potential explanations. One emerging hypothesis, ghost introgression, refers to the genetic contribution of an ancient, extinct or unsampled lineage and can now be tested using modern genomic data and demographic models. In this study, we investigated the evolutionary history of the Odorrana swinhoana complex (Anura: Ranidae), which includes O. swinhoana, O. utsunomiyaorum and an unidentified population with highly divergent mtDNA. While mitochondrial phylogeny suggested this population as a basal lineage, nuclear data from ddRADseq revealed it as a mixture of the most derived O. swinhoana nuclear sequences combined with ancient mtDNA. Demographic modelling further supported ghost introgression, as all models incorporating a ghost population outperformed those without it. These findings suggest that an eastward expansion of western O. swinhoana replaced an ancient Odorrana lineage, leaving only its mtDNA and fragments of its nuclear genome in the hybrid population. Our results provide one of the first documented cases of ghost introgression in amphibians and highlight its potential as a widespread evolutionary process. This study also underscores the risks of relying solely on mtDNA for phylogenetic reconstruction and species delimitation.

Reconciling Gene Tree Discordance and Biogeography in European Crows

Reconstructing the evolutionary history of young lineages diverging with gene flow is challenging due to factors like incomplete lineage sorting, introgression, and selection causing gene tree discordance. The European crow hybrid zone between all-black carrion crows and grey-coated hooded crows exemplifies this challenge. Most of the genome in Western and Central European carrion crow populations is near-identical to hooded crows, but differs substantially from their Iberian congeners. A notable exception is a single major-effect colour-locus under sexual selection aligning with the 'species' tree. To understand the underlying evolutionary processes, we reconstructed the biogeographic history of the species complex. During the Pleistocene carrion and hooded crows took refuge in the Iberian Peninsula and the Middle East, respectively. Allele-sharing of all-black Western European populations with likewise black Iberian crows at the colour-locus represents the last trace of carrion crow ancestry, resisting gene flow from expanding hooded crow populations that have homogenised most of the genome. A model of colour-locus introgression from an Iberian ancestor into hooded crow populations near the Pyrenées was significantly less supported. We found no positive relationship between introgression and recombination rate consistent with the absence of genome-wide, polygenic barriers in this young species complex. Overall, this study portrays a scenario where few large-effect loci, subject to divergent sexual selection, resist rampant and asymmetric gene exchange. This study underscores the importance of integrating population demography and biogeography to accurately interpret patterns of gene tree discordance following population divergence.

Neanderthal adaptive introgression shaped LCT enhancer region diversity without linking to lactase persistence in East Asian populations

Positive selection at the 2q21.3 enhancer region for lactase gene (LCT) expression in Europeans and Africans has long been attributed to selection for lactase persistence (LP), the capacity of adults to digest lactose in milk, presumably because of the benefits associated with milk consumption. While considered a classic example of gene-culture coevolution, recently doubts have been raised about the link between selection at 2q21.3 and LP. Analysis of additional populations could shed further light; here, we demonstrate that a haplotype spanning ~467 kb at the 2q21.3 locus has risen to high frequency in East Asians (~25%) but is absent from Africans and Europeans. This haplotype likely derived from Neanderthals and has been under positive selection in East Asians. The East Asian-specific haplotype is associated with alterations in LCT expression and promoter methylation in certain cell types, similar to what is observed with LP-associated haplotypes in Europeans. Moreover, its frequency is comparable to that of LP in East Asians, suggesting a potential association with LP in East Asians. However, it is highly unlikely that selection in East Asians was related to milk-drinking habits. We find that this haplotype impacts the expression of UBXN4, DARS1, and DARS1-AS1 in immune cells and is associated with neutrophil and white blood cell counts. Hence, the selection might be linked to certain aspects of immune function. This implies that selection on 2q21.3 has thus either occurred for different reasons in different populations or the selection observed in other populations is also not due to LP.

Archaic introgression and the distribution of shared variation under stabilizing selection

Many phenotypic traits are under stabilizing selection, which maintains a population's mean phenotypic value near some optimum. The dynamics of traits and trait architectures under stabilizing selection have been extensively studied for single populations at steady state. However, natural populations are seldom at steady state and are often structured in some way. Admixture and introgression events may be common, including over human evolutionary history. Because stabilizing selection results in selection against the minor allele at a trait-affecting locus, alleles from the minor parental ancestry will be selected against after admixture. We show that the site-frequency spectrum can be used to model the genetic architecture of such traits, allowing for the study of trait architecture dynamics in complex multi-population settings. We use a simple deterministic two-locus model to predict the reduction of introgressed ancestry around trait-contributing loci. From this and individual-based simulations, we show that introgressed-ancestry is depleted around such loci. When introgression between two diverged populations occurs in both directions, as has been inferred between humans and Neanderthals, the locations of such regions with depleted introgressed ancestry will tend to be shared across populations. We argue that stabilizing selection for shared phenotypic optima may explain recent observations in which regions of depleted human-introgressed ancestry in the Neanderthal genome overlap with Neanderthal-ancestry deserts in humans.

A refined analysis of Neanderthal-introgressed sequences in modern humans with a complete reference genome

BACKGROUND

Leveraging long-read sequencing technologies, the first complete human reference genome, T2T-CHM13, corrects assembly errors in previous references and resolves the remaining 8% of the genome. While studies on archaic admixture in modern humans have so far relied on the GRCh37 reference due to the availability of archaic genome data, the impact of T2T-CHM13 in this field remains unexplored.

RESULTS

We remap the sequencing reads of the high-quality Altai Neanderthal and Denisovan genomes onto GRCh38 and T2T-CHM13. Compared to GRCh37, we find that T2T-CHM13 significantly improves read mapping quality in archaic samples. We then apply IBDmix to identify Neanderthal-introgressed sequences in 2504 individuals from 26 geographically diverse populations using different reference genomes. We observe that commonly used pre-phasing filtering strategies in public datasets substantially influence archaic ancestry determination, underscoring the need for careful filter selection. Our analysis identifies approximately 51 Mb of Neanderthal sequences unique to T2T-CHM13, predominantly in genomic regions where GRCh38 and T2T-CHM13 assemblies diverge. Additionally, we uncover novel instances of population-specific archaic introgression in diverse populations, spanning genes involved in metabolism, olfaction, and ion-channel function. Finally, to facilitate the exploration of archaic alleles and adaptive signals in human genomics and evolutionary research, we integrate these introgressed sequences and adaptive signals across all reference genomes into a visualization database, ASH ( www.arcseqhub.com ).

CONCLUSIONS

Our study enhances the detection of archaic variations in modern humans, highlights the importance of utilizing the T2T-CHM13 reference, and provides novel insights into the functional consequences of archaic hominin admixture.

Human dispersal into East Eurasia: ancient genome insights and the need for research on physiological adaptations

Humans have long pondered their genesis. The answer to the great question of where Homo sapiens come from has evolved in conjunction with biotechnologies that have allowed us to more brightly illuminate our distant past. The "Multiregional Evolution" model was once the hegemonic theory of Homo sapiens origins, but in the last 30 years, it has been supplanted by the "Out of Africa" model. Here, we review the major findings that have resulted in this paradigmatic shift. These include hominin brain expansion, classical insight from the mitochondrial genome (mtDNA) regarding the timing of the divergence point between Africans and non-Africans, and next-generation sequencing (NGS) of the Neanderthal and Denisovan genomes. These findings largely bolstered the "Out of Africa" model, although they also revealed a small degree of introgression of the Neanderthal and Denisovan genomes into those of non-African Homo sapiens. We also review paleogenomic studies for which migration route, north or south, early migrants to East Eurasia most likely traversed. Whichever route was taken, the migrants moved to higher latitudes, which necessitated adaptation for lower light conditions, colder clines, and pro-adipogenic mechanisms to counteract food scarcity. Further genetic and epigenetic investigations of these physiological adaptations constitute an integral aspect of the story of human origins and human migration to East Asia.

The Chimpanzee Brainnetome Atlas reveals distinct connectivity and gene expression profiles relative to humans

Chimpanzees (Pan troglodytes) are one of humans' closest living relatives, making them the most directly relevant comparison point for understanding human brain evolution. Zeroing in on the differences in brain connectivity between humans and chimpanzees can provide key insights into the specific evolutionary changes that might have occurred along the human lineage. However, such comparisons are hindered by the absence of cross-species brain atlases established within the same framework. To address this gap, we developed the Chimpanzee Brainnetome Atlas (ChimpBNA) using a connectivity-based parcellation framework. Leveraging this new resource, we found substantial divergence in connectivity patterns between the two species across most association cortices, notably in the lateral temporal and dorsolateral prefrontal cortex. These differences deviate sharply from the pattern of cortical expansion observed when comparing humans to chimpanzees, highlighting more complex and nuanced connectivity changes in brain evolution than previously recognized. Additionally, we identified regions displaying connectional asymmetries that differed between species, likely resulting from evolutionary divergence. Genes highly expressed in regions of divergent connectivities were enriched in cell types crucial for cortical projection circuits and synapse formation, whose pronounced differences in expression patterns hint at genetic influences on neural circuit development, function, and evolution. Our study provides a fine-scale chimpanzee brain atlas and highlights the chimpanzee-human connectivity divergence in a rigorous and comparative manner. In addition, these results suggest potential gene expression correlates for species-specific differences by linking neuroimaging and genetic data, offering insights into the evolution of human-unique cognitive capabilities.

The shared genetic architecture and evolution of human language and musical rhythm

This study aimed to test theoretical predictions over biological underpinnings of previously documented phenotypic correlations between human language-related and musical rhythm traits. Here, after identifying significant genetic correlations between rhythm, dyslexia and various language-related traits, we adapted multivariate methods to capture genetic signals common to genome-wide association studies of rhythm (N = 606,825) and dyslexia (N = 1,138,870). The results revealed 16 pleiotropic loci (P < 5 × 10-8) jointly associated with rhythm impairment and dyslexia, and intricate shared genetic and neurobiological architectures. The joint genetic signal was enriched for foetal and adult brain cell-specific regulatory regions, highlighting complex cellular composition in their shared underpinnings. Local genetic correlation with a key white matter tract (the left superior longitudinal fasciculus-I) substantiated hypotheses about auditory-motor connectivity as a genetically influenced, evolutionarily relevant neural endophenotype common to rhythm and language processing. Overall, we provide empirical evidence of multiple aspects of shared biology linking language and musical rhythm, contributing novel insight into the evolutionary relationships between human musicality and linguistic communication traits.

Archaic hominin admixture and its consequences for modern humans

As anatomically modern humans dispersed out of Africa, they encountered and mated with now extinct hominins, including Neanderthals and Denisovans. It is now well established that all non-African individuals derive approximately 2% of their genome from Neanderthal ancestors and individuals of Melanesian and Australian aboriginal ancestry inherited an additional 2%-5% of their genomes from Denisovan ancestors. Attention has started to shift from documenting amounts of archaic admixture and identifying introgressed segments to understanding their molecular, phenotypic, and evolutionary consequences and refining models of human history. Here, we review recent insights into admixture between modern and archaic humans, emphasizing methodological innovations and the functional and phenotypic effects Neanderthal and Denisovan sequences have in contemporary individuals.

Ignoring population structure in hominin evolutionary models can lead to the inference of spurious admixture events

Genomic and ancient DNA data have revolutionized palaeoanthropology and our vision of human evolution, with indisputable landmarks like the sequencing of Neanderthal and Denisovan genomes. Yet, using genetic data to identify, date and quantify evolutionary events-such as ancient bottlenecks or admixture-is not straightforward, as inferences may depend on model assumptions. In the last two decades, the idea that Neanderthals and members of the Homo sapiens lineage interbred has gained momentum. From the status of unlikely theory, it has reached consensus among human evolutionary biologists. This theory is mainly supported by statistical approaches that depend on demographic models minimizing or ignoring population structure, despite its widespread occurrence and the fact that, when ignored, population structure can lead to the inference of spurious demographic events. We simulated genomic data under a structured and admixture-free model of human evolution, and found that all the tested admixture approaches identified long Neanderthal fragments in our simulated genomes and an admixture event that never took place. We also observed that several published admixture models failed to predict important empirical diversity or admixture statistics, and that we could identify several scenarios from our structured model that better predicted these statistics jointly. Using a simulated time series of ancient DNA, the structured scenarios could also predict the trajectory of the empirical D statistics. Our results suggest that models accounting for population structure are fundamental to improve our understanding of human evolution, and that admixture between Neanderthals and H. sapiens needs to be re-evaluated in the light of structured models. Beyond the Neanderthal case, we argue that ancient hybridization events, which are increasingly documented in many species, including with other hominins, may also benefit from such re-evaluation.

PITX2 expression and Neanderthal introgression in HS3ST3A1 contribute to variation in tooth dimensions in modern humans

Dental morphology varies greatly throughout evolution, including in the human lineage, but little is known about the biology of this variation. Here, we use multiomics analyses to examine the genetics of variation in tooth crown dimensions. In a human cohort with mixed continental ancestry, we detected genome-wide significant associations at 18 genome regions. One region includes EDAR, a gene known to impact dental features in East Asians. Furthermore, we find that EDAR variants increase the mesiodistal diameter of all teeth, following an anterior-posterior gradient of decreasing strength. Among the 17 novel-associated regions, we replicate 7/13 in an independent human cohort and find that 4/12 orthologous regions affect molar size in mice. Two association signals point to compelling candidate genes. One is ∼61 kb from PITX2, a major determinant of tooth development. Another overlaps HS3ST3A1, a paralogous neighbor of HS3ST3B1, a tooth enamel knot factor. We document the expression of Pitx2 and Hs3st3a1 in enamel knot and dental epithelial cells of developing mouse incisors. Furthermore, associated SNPs in PITX2 and HS3ST3A1 overlap enhancers active in these cells, suggesting a role for these SNPs in gene regulation during dental development. In addition, we document that Pitx2 and Hs3st3a1/Hs3st3b1 knockout mice show alterations in dental morphology. Finally, we find that associated SNPs in HS3ST3A1 are in a DNA tract introgressed from Neanderthals, consistent with an involvement of HS3ST3A1 in tooth size variation during human evolution.

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

Gene flow from Neanderthals has shaped genetic and phenotypic variation in modern humans. We generated a catalog of Neanderthal ancestry segments in more than 300 genomes spanning the past 50,000 years. We examined how Neanderthal ancestry is shared among individuals over time. Our analysis revealed that the vast majority of Neanderthal gene flow is attributable to a single, shared extended period of gene flow that occurred between 50,500 to 43,500 years ago, as evidenced by ancestry correlation, colocalization of Neanderthal segments across individuals, and divergence from the sequenced Neanderthals. Most natural selection-positive and negative-on Neanderthal variants occurred rapidly after the gene flow. Our findings provide new insights into how contact with Neanderthals shaped modern human origins and adaptation.

The MUC19 gene in Denisovans, Neanderthals, and Modern Humans: An Evolutionary History of Recurrent Introgression and Natural Selection

We study the gene MUC19, for which modern humans carry a Denisovan-like haplotype. MUC19 is a mucin, a glycoprotein that forms gels with various biological functions. We find the diagnostic variants for the Denisovan-like MUC19 haplotype at high frequencies in admixed Latin American individuals among global populations, and at highest frequency in 23 ancient Indigenous American individuals, all predating population admixture with Europeans and Africans. We find that the Denisovan-like MUC19 haplotype carries a higher copy number of a 30 base-pair variable number tandem repeat, and that copy numbers of this repeat are exceedingly high in American populations and are under positive selection. This study provides the first example of positive selection acting on archaic alleles at coding sites and VNTRs. Finally, we find that some Neanderthals carry the Denisovan-like MUC19 haplotype, and that it was likely introgressed into human populations through Neanderthal introgression rather than Denisovan introgression.

Computational Genomics and Its Applications to Anthropological Questions

The advent of affordable genome sequencing and the development of new computational tools have established a new era of genomic knowledge. Sequenced human genomes number in the tens of thousands, including thousands of ancient human genomes. The abundance of data has been met with new analysis tools that can be used to understand populations' demographic and evolutionary histories. Thus, a variety of computational methods now exist that can be leveraged to answer anthropological questions. This includes novel likelihood and Bayesian methods, machine learning techniques, and a vast array of population simulators. These computational tools provide powerful insights gained from genomic datasets, although they are generally inaccessible to those with less computational experience. Here, we outline the theoretical workings behind computational genomics methods, limitations and other considerations when applying these computational methods, and examples of how computational methods have already been applied to anthropological questions. We hope this review will empower other anthropologists to utilize these powerful tools in their own research.

Haplotype-Based Approach Represents Locus Specificity in the Genomic Diversification Process in Humans (Homo sapiens)

BACKGROUND/OBJECTIVES

Recent progress in evolutionary genomics on human (Homo sapiens) populations has revealed complex demographic events and genomic changes. These include population expansion with complicated migration, substantial population structure, and ancient introgression from other hominins, as well as human characteristics selections. Nevertheless, the genomic regions in which such evolutionary events took place have remained unclear.

METHODS

Here, we focused on eight loci containing the haplotypes that were previously presented as atypical for the mutation pattern in sequence and/or geographic distribution pattern with the model of recent African origin, which constitute two major clusters: African only, and global. This was the consensus model before information regarding introgression from Neanderthal (Homo neanderthalensis) was available. We compared diversity in identical datasets of the modern human population genome, with the 1000 Genomes project among them.

RESULTS/CONCLUSIONS

This study identified representative genomic regions that show traces of various demographic events and genomic changes that modern humans have undergone by categorizing the relationships in sequence similarity and in worldwide geographic distribution among haplotypes.

Molecular characterization of the archaic HLA-B*73:01 allele reveals presentation of a unique peptidome and skewed engagement by KIR2DL2

HLA class I alleles of archaic origin may have been retained in modern humans because they provide immunity against diseases to which archaic humans had evolved resistance. According to this model, archaic introgressed alleles were somehow distinct from those that evolved in African populations. Here we show that HLA-B*73:01, a rare allotype with putative archaic origins, has a relatively rare peptide binding motif with an unusually long-tailed peptide length distribution. We also find that HLA-B*73:01 combines a restricted and unique peptidome with high-cell surface expression, characteristics that make it well-suited to combat one or a number of closely-related pathogens. Furthermore, a crystal structure of HLA-B*73:01 in complex with KIR2DL2 highlights differences from previously solved structures with HLA-C molecules. These molecular characteristics distinguish HLA-B*73:01 from other HLA class I alleles previously investigated and may have provided early modern human migrants that inherited this allele with a selective advantage as they colonized Europe and Asia.

Archaic introgression contributed to shape the adaptive modulation of angiogenesis and cardiovascular traits in human high-altitude populations from the Himalayas

It is well established that several Homo sapiens populations experienced admixture with extinct human species during their evolutionary history. Sometimes, such a gene flow could have played a role in modulating their capability to cope with a variety of selective pressures, thus resulting in archaic adaptive introgression events. A paradigmatic example of this evolutionary mechanism is offered by the EPAS1 gene, whose most frequent haplotype in Himalayan highlanders was proved to reduce their susceptibility to chronic mountain sickness and to be introduced in the gene pool of their ancestors by admixture with Denisovans. In this study, we aimed at further expanding the investigation of the impact of archaic introgression on more complex adaptive responses to hypobaric hypoxia evolved by populations of Tibetan/Sherpa ancestry, which have been plausibly mediated by soft selective sweeps and/or polygenic adaptations rather than by hard selective sweeps. For this purpose, we used a combination of composite-likelihood and gene network-based methods to detect adaptive loci in introgressed chromosomal segments from Tibetan WGS data and to shortlist those presenting Denisovan-like derived alleles that participate to the same functional pathways and are absent in populations of African ancestry, which are supposed to do not have experienced Denisovan admixture. According to this approach, we identified multiple genes putatively involved in archaic introgression events and that, especially as regards TBC1D1, RASGRF2, PRKAG2, and KRAS, have plausibly contributed to shape the adaptive modulation of angiogenesis and of certain cardiovascular traits in high-altitude Himalayan peoples. These findings provided unprecedented evidence about the complexity of the adaptive phenotype evolved by these human groups to cope with challenges imposed by hypobaric hypoxia, offering new insights into the tangled interplay of genetic determinants that mediates the physiological adjustments crucial for human adaptation to the high-altitude environment.

Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

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.

Adaptive Evolution of Two Distinct Adaptive Haplotypes of Neanderthal Origin at the Immunoglobulin Heavy-chain Locus in East Asian and European Populations

Immunoglobulins (Igs) have a crucial role in humoral immunity. Two recent studies have reported a high-frequency Neanderthal-introgressed haplotype throughout Eurasia and a high-frequency Neanderthal-introgressed haplotype specific to southern East Asia at the immunoglobulin heavy-chain (IGH) gene locus on chromosome 14q32.33. Surprisingly, we found the previously reported high-frequency Neanderthal-introgressed haplotype does not exist throughout Eurasia. Instead, our study identified two distinct high-frequency haplotypes of putative Neanderthal origin in East Asia and Europe, although they shared introgressed alleles. Notably, the alleles of putative Neanderthal origin reduced the expression of IGHG1 and increased the expression of IGHG2 and IGHG3 in various tissues. These putatively introgressed alleles also affected the production of IgG1 upon antigen stimulation and increased the risk of systemic lupus erythematosus. Additionally, the greatest genetic differentiation across the whole genome between southern and northern East Asians was observed for the East Asian haplotype of putative Neanderthal origin. The frequency decreased from southern to northern East Asia and correlated positively with the genome-wide proportion of southern East Asian ancestry, indicating that this putative positive selection likely occurred in the common ancestor of southern East Asian populations before the admixture with northern East Asian populations.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

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.

Selection Shapes the Genomic Landscape of Introgressed Ancestry in a Pair of Sympatric Sea Urchin Species

A growing number of recent studies have demonstrated that introgression is common across the tree of life. However, we still have a limited understanding of the fate and fitness consequence of introgressed variation at the whole-genome scale across diverse taxonomic groups. Here, we implemented a phylogenetic hidden Markov model to identify and characterize introgressed genomic regions in a pair of well-diverged, nonsister sea urchin species: Strongylocentrotus pallidus and Strongylocentrotus droebachiensis. Despite the old age of introgression, a sizable fraction of the genome (1% to 5%) exhibited introgressed ancestry, including numerous genes showing signals of historical positive selection that may represent cases of adaptive introgression. One striking result was the overrepresentation of hyalin genes in the identified introgressed regions despite observing considerable overall evidence of selection against introgression. There was a negative correlation between introgression and chromosome gene density, and two chromosomes were observed with considerably reduced introgression. Relative to the nonintrogressed genome-wide background, introgressed regions had significantly reduced nucleotide divergence (dXY) and overlapped fewer protein-coding genes, coding bases, and genes with a history of positive selection. Additionally, genes residing within introgressed regions showed slower rates of evolution (dN, dS, dN/dS) than random samples of genes without introgressed ancestry. Overall, our findings are consistent with widespread selection against introgressed ancestry across the genome and suggest that slowly evolving, low-divergence genomic regions are more likely to move between species and avoid negative selection following hybridization and introgression.

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.

A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

The temporal and genomic scale of selection following hybridization

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.

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.

Neuroimaging genomics as a window into the evolution of human sulcal organization

Primate brain evolution has involved prominent expansions of the cerebral cortex, with largest effects observed in the human lineage. Such expansions were accompanied by fine-grained anatomical alterations, including increased cortical folding. However, the molecular bases of evolutionary alterations in human sulcal organization are not yet well understood. Here, we integrated data from recently completed large-scale neuroimaging genetic analyses with annotations of the human genome relevant to various periods and events in our evolutionary history. These analyses identified single-nucleotide polymorphism (SNP) heritability enrichments in fetal brain human-gained enhancer (HGE) elements for a number of sulcal structures, including the central sulcus, which is implicated in human hand dexterity. We zeroed in on a genomic region that harbors DNA variants associated with left central sulcus shape, an HGE element, and genetic loci involved in neurogenesis including ZIC4, to illustrate the value of this approach for probing the complex factors contributing to human sulcal evolution.

Leveraging shared ancestral variation to detect local introgression

Introgression is a common evolutionary phenomenon that results in shared genetic material across non-sister taxa. Existing statistical methods such as Patterson's D statistic can detect introgression by measuring an excess of shared derived alleles between populations. The D statistic is effective to detect genome-wide patterns of introgression but can give spurious inferences of introgression when applied to local regions. We propose a new statistic, D+, that leverages both shared ancestral and derived alleles to infer local introgressed regions. Incorporating both shared derived and ancestral alleles increases the number of informative sites per region, improving our ability to identify local introgression. We use a coalescent framework to derive the expected value of this statistic as a function of different demographic parameters under an instantaneous admixture model and use coalescent simulations to compute the power and precision of D+. While the power of D and D+ is comparable, D+ has better precision than D. We apply D+ to empirical data from the 1000 Genome Project and Heliconius butterflies to infer local targets of introgression in humans and in butterflies.

Genome evolution is surprisingly predictable after initial hybridization

Over the past two decades, evolutionary biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common - not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past ~40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni×X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

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

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

Archaic Introgression Shaped Human Circadian Traits

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Disentangling archaic introgression and genomic signatures of selection at human immunity genes

Pathogens and infectious diseases have imposed exceptionally strong selective pressure on ancient and modern human genomes and contributed to the current variation in many genes. There is evidence that modern humans acquired immune variants through interbreeding with ancient hominins, but the impact of such variants on human traits is not fully understood. The main objectives of this research were to infer the genetic signatures of positive selection that may be involved in adaptation to infectious diseases and to investigate the function of Neanderthal alleles identified within a set of 50 Lithuanian genomes. Introgressed regions were identified using the machine learning tool ArchIE. Recent positive selection signatures were analysed using iHS. We detected high-scoring signals of positive selection at innate immunity genes (EMB, PARP8, HLAC, and CDSN) and evaluated their interactions with the structural proteins of pathogens. Interactions with human immunodeficiency virus (HIV) 1 and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were identified. Overall, genomic regions introgressed from Neanderthals were shown to be enriched in genes related to immunity, keratinocyte differentiation, and sensory perception.

Diverse African genomes reveal selection on ancient modern human introgressions in Neanderthals

Comparisons of Neanderthal genomes to anatomically modern human (AMH) genomes show a history of Neanderthal-to-AMH introgression stemming from interbreeding after the migration of AMHs from Africa to Eurasia. All non-sub-Saharan African AMHs have genomic regions genetically similar to Neanderthals that descend from this introgression. Regions of the genome with Neanderthal similarities have also been identified in sub-Saharan African populations, but their origins have been unclear. To better understand how these regions are distributed across sub-Saharan Africa, the source of their origin, and what their distribution within the genome tells us about early AMH and Neanderthal evolution, we analyzed a dataset of high-coverage, whole-genome sequences from 180 individuals from 12 diverse sub-Saharan African populations. In sub-Saharan African populations with non-sub-Saharan African ancestry, as much as 1% of their genomes can be attributed to Neanderthal sequence introduced by recent migration, and subsequent admixture, of AMH populations originating from the Levant and North Africa. However, most Neanderthal homologous regions in sub-Saharan African populations originate from migration of AMH populations from Africa to Eurasia ∼250 kya, and subsequent admixture with Neanderthals, resulting in ∼6% AMH ancestry in Neanderthals. These results indicate that there have been multiple migration events of AMHs out of Africa and that Neanderthal and AMH gene flow has been bi-directional. Observing that genomic regions where AMHs show a depletion of Neanderthal introgression are also regions where Neanderthal genomes show a depletion of AMH introgression points to deleterious interactions between introgressed variants and background genomes in both groups-a hallmark of incipient speciation.

The shared genetic architecture and evolution of human language and musical rhythm

Rhythm and language-related traits are phenotypically correlated, but their genetic overlap is largely unknown. Here, we leveraged two large-scale genome-wide association studies performed to shed light on the shared genetics of rhythm (N=606,825) and dyslexia (N=1,138,870). Our results reveal an intricate shared genetic and neurobiological architecture, and lay groundwork for resolving longstanding debates about the potential co-evolution of human language and musical traits.

Past human expansions shaped the spatial pattern of Neanderthal ancestry

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

A Catalog of Coding Sequence Variations in Salivary Proteins' Genes Occurring during Recent Human Evolution

Saliva houses over 2000 proteins and peptides with poorly clarified functions, including proline-rich proteins, statherin, P-B peptides, histatins, cystatins, and amylases. Their genes are poorly conserved across related species, reflecting an evolutionary adaptation. We searched the nucleotide substitutions fixed in these salivary proteins' gene loci in modern humans compared with ancient hominins. We mapped 3472 sequence variants/nucleotide substitutions in coding, noncoding, and 5'-3' untranslated regions. Despite most of the detected variations being within noncoding regions, the frequency of coding variations was far higher than the general rate found throughout the genome. Among the various missense substitutions, specific substitutions detected in PRB1 and PRB2 genes were responsible for the introduction/abrogation of consensus sequences recognized by convertase enzymes that cleave the protein precursors. Overall, these changes that occurred during the recent human evolution might have generated novel functional features and/or different expression ratios among the various components of the salivary proteome. This may have influenced the homeostasis of the oral cavity environment, possibly conditioning the eating habits of modern humans. However, fixed nucleotide changes in modern humans represented only 7.3% of all the substitutions reported in this study, and no signs of evolutionary pressure or adaptative introgression from archaic hominins were found on the tested genes.

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.

Archaic Introgression Shaped Human Circadian Traits

Introduction

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultra-violet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology, and whether archaic introgression adaptively contributed to human chronotypes remains unknown.

Results

Here we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, RORC), and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among eQTLs for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency.

Conclusions

These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

TP53 germline pathogenic variants in modern humans were likely originated during recent human history

TP53 is crucial for maintaining genome stability and preventing oncogenesis. Germline pathogenic variation in TP53 damages its function, causing genome instability and increased cancer risk. Despite extensive study in TP53, the evolutionary origin of the human TP53 germline pathogenic variants remains largely unclear. In this study, we applied phylogenetic and archaeological approaches to identify the evolutionary origin of TP53 germline pathogenic variants in modern humans. In the phylogenic analysis, we searched 406 human TP53 germline pathogenic variants in 99 vertebrates distributed in eight clades of Primate, Euarchontoglires, Laurasiatheria, Afrotheria, Mammal, Aves, Sarcopterygii and Fish, but we observed no direct evidence for the cross-species conservation as the origin; in the archaeological analysis, we searched the variants in 5031 ancient human genomes dated between 45045 and 100 years before present, and identified 45 pathogenic variants in 62 ancient humans dated mostly within the last 8000 years; we also identified 6 pathogenic variants in 3 Neanderthals dated 44000 to 38515 years before present and 1 Denisovan dated 158 550 years before present. Our study reveals that TP53 germline pathogenic variants in modern humans were likely originated in recent human history and partially inherited from the extinct Neanderthals and Denisovans.

Ghost admixture in eastern gorillas

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic 'ghost' lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species.

Engineering Human Pluripotent Stem Cell Lines to Evade Xenogeneic Transplantation Barriers

Allogeneic human pluripotent stem cell (hPSC)-derived cells and tissues for therapeutic transplantation must necessarily overcome immunological rejection by the recipient. To define these barriers and to create cells capable of evading rejection for preclinical testing in immunocompetent mouse models, we genetically ablated β2m, Tap1, Ciita, Cd74, Mica, and Micb to limit expression of HLA-I, HLA-II, and natural killer cell activating ligands in hPSCs. Though these and even unedited hPSCs readily formed teratomas in cord blood-humanized immunodeficient mice, grafts were rapidly rejected by immunocompetent wild-type mice. Transplantation of these cells that also expressed covalent single chain trimers of Qa1 and H2-Kb to inhibit natural killer cells and CD55, Crry, and CD59 to inhibit complement deposition led to persistent teratomas in wild-type mice. Expression of additional inhibitory factors such as CD24, CD47, and/or PD-L1 had no discernible impact on teratoma growth or persistence. Transplantation of HLA-deficient hPSCs into mice genetically deficient in complement and depleted of natural killer cells also led to persistent teratomas. Thus, T cell, NK cell, and complement evasion are necessary to prevent immunological rejection of hPSCs and their progeny. These cells and versions expressing human orthologs of immune evasion factors can be used to refine tissue- and cell type-specific immune barriers, and to conduct preclinical testing in immunocompetent mouse models.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

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

The temporal and genomic scale of selection following hybridization

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the dynamics underlying these patterns within hybrid populations have been lacking. Here, we develop methods based on the Wavelet Transform to understand the spatial genomic scale of local ancestry variation and its association with recombination rates. We present theory and use simulations to show how wavelet-based decompositions of ancestry variance along the genome and the correlation between ancestry and recombination reflect the joint effects of recombination, genetic drift, and genome-wide selection against introgressed alleles. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of local ancestry variation at varying spatial genomic scales through time. Using wavelet approaches to identify the genomic scale of variance in ancestry and its correlates, we show that these methods can detect temporally localized effects of drift and selection. We apply these methods to previously published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio), and to inferred Neanderthal introgression in modern humans. Across systems, we find that upwards of 20% of the variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. We also see signals of selection at fine genomic scales and much longer time scales. However, we show that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available, and can help shed light on generalities of the genomic consequences of interspecific hybridization.

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 role of geography, ecology, and hybridization in the evolutionary history of Canary Island Descurainia

PREMISE

Oceanic islands offer the opportunity to understand evolutionary processes underlying rapid diversification. Along with geographic isolation and ecological shifts, a growing body of genomic evidence has suggested that hybridization can play an important role in island evolution. Here we use genotyping-by-sequencing (GBS) to understand the roles of hybridization, ecology, and geographic isolation in the radiation of Canary Island Descurainia (Brassicaceae).

METHODS

We carried out GBS for multiple individuals of all Canary Island species and two outgroups. Phylogenetic analyses of the GBS data were performed using both supermatrix and gene tree approaches and hybridization events were examined using D-statistics and Approximate Bayesian Computation. Climatic data were analyzed to examine the relationship between ecology and diversification.

RESULTS

Analysis of the supermatrix data set resulted in a fully resolved phylogeny. Species networks suggest a hybridization event has occurred for D. gilva, with these results being supported by Approximate Bayesian Computation analysis. Strong phylogenetic signals for temperature and precipitation indicate one major ecological shift within Canary Island Descurainia.

CONCLUSIONS

Inter-island dispersal played a significant role in the diversification of Descurainia, with evidence of only one major shift in climate preferences. Despite weak reproductive barriers and the occurrence of hybrids, hybridization appears to have played only a limited role in the diversification of the group with a single instance detected. The results highlight the need to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups prone to hybridization; patterns that might otherwise be obscured in species trees.

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.

Modern Japanese ancestry-derived variants reveal the formation process of the current Japanese regional gradations

Modern Japanese people have two major ancestral populations: indigenous Jomon hunter-gatherers and continental East Asian farmers. To determine the formation process of the current Japanese population, we developed a detection method for variants derived from ancestral populations using a summary statistic, the ancestry marker index (AMI). We applied AMI to modern Japanese population samples and identified 208,648 single nucleotide polymorphisms (SNPs) that were likely derived from the Jomon people (Jomon-derived variants). Analysis of Jomon-derived variants in 10,842 modern Japanese individuals recruited from all over Japan revealed that the admixture proportions of the Jomon people varied between prefectures, probably owing to the prehistoric population size difference. The estimated allele frequencies of genome-wide SNPs in the ancestral populations of the modern Japanese suggested their adaptive phenotypic characteristics to their respective livelihoods. Based on our findings, we propose a formation model for the genotypic and phenotypic gradations of the current Japanese archipelago populations.

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.