Evidence that RNA Viruses Drove Adaptive Introgression between Neanderthals and Modern Humans

Tracing the Evolution of Human Immunity Through Ancient DNA

Infections have imposed strong selection pressures throughout human evolution, making the study of natural selection's effects on immunity genes highly complementary to disease-focused research. This review discusses how ancient DNA studies, which have revolutionized evolutionary genetics, increase our understanding of the evolution of human immunity. These studies have shown that interbreeding between modern humans and Neanderthals or Denisovans has influenced present-day immune responses, particularly to viruses. Additionally, ancient genomics enables the tracking of how human immunity has evolved across cultural transitions, highlighting strong selection since the Bronze Age in Europe (<4,500 years) and potential genetic adaptations to epidemics raging during the Middle Ages and the European colonization of the Americas. Furthermore, ancient genomic studies suggest that the genetic risk for noninfectious immune disorders has gradually increased over millennia because alleles associated with increased risk for autoimmunity and inflammation once conferred resistance to infections. The challenge now is to extend these findings to diverse, non-European populations and to provide a more global understanding of the evolution of human immunity.

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.

Making sense of sentinels: wildlife as the One Health bridge for environmental antimicrobial resistance surveillance

Antimicrobial resistance (AMR), arising from decades of imprudent anthropogenic use of antimicrobials in healthcare and agriculture, is considered one of the greatest One Health crises facing healthcare globally. Antimicrobial pollutants released from human-associated sources are intensifying resistance evolution in the environment. Due to various ecological factors, wildlife interact with these polluted ecosystems, acquiring resistant bacteria and genes. Although wildlife are recognized reservoirs and disseminators of AMR in the environment, current AMR surveillance systems still primarily focus on clinical and agricultural settings, neglecting this environmental dimension. Wildlife can serve as valuable sentinels of AMR in the environment, reflecting ecosystem health, and the effectiveness of mitigation strategies. This review explores knowledge gaps surrounding the ecological factors influencing AMR acquisition and dissemination in wildlife, and highlights limitations in current surveillance systems and policy instruments that do not sufficiently address the environmental component of AMR. We discuss the underutilized opportunity of using wildlife as sentinel species in a holistic, One Health-centred AMR surveillance system. By better integrating wildlife into systematic AMR surveillance and policy, and leveraging advances in high-throughput technologies, we can track and predict resistance evolution, assess the ecological impacts, and better understand the complex dynamics of environmental transmission of AMR across ecosystems.

Role of the Neanderthal Genome in Genetic Susceptibility to COVID-19: 3p21.31 Locus in the Spotlight

Since the outbreak of COVID-19, genome-wide association studies have tried to discover the role of genetic predisposition in the clinical variability of this viral infection. The findings of various investigations have led to several loci for COVID-19 genetic susceptibility. Among candidate regions, the 3p21.31 locus has been in the spotlight among scientists, as it can increase the risk of severe COVID-19 by almost two fold. In addition to its substantial association with COVID-19 severity, this locus is related to some common diseases, such as diabetes, malignancies, and coronary artery disease. This locus also harbors evolutionary traces of Neanderthal genomes, which is believed to be the underlying reason for its association with COVID-19 severity. Additionally, the inheritance of this locus from Neanderthals seems to be under positive selection. This review aims to summarize a collection of evidence on the 3p21.31 locus and its impact on COVID-19 outcomes by focusing on the risk variants originated from the Neanderthal genome. Moreover, we discuss candidate genes at this locus and the possible mechanisms by which they influence the progression of COVID-19 symptoms. Better insights into human genetic susceptibility to newly emerging diseases such as COVID-19 and its evolutionary origin can provide fundamentals for risk assessment of different populations as well as the development of personalized prevention and treatments based on genomic medicine.

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.

Extensive longevity and DNA virus-driven adaptation in nearctic Myotis bats

The genus Myotis is one of the largest clades of bats, and exhibits some of the most extreme variation in lifespans among mammals alongside unique adaptations to viral tolerance and immune defense. To study the evolution of longevity-associated traits and infectious disease, we generated near-complete genome assemblies and cell lines for 8 closely related species of Myotis. Using genome-wide screens of positive selection, analyses of structural variation, and functional experiments in primary cell lines, we identify new patterns of adaptation contributing to longevity, cancer resistance, and viral interactions in bats. We find that Myotis bats have some of the most significant variation in cancer risk across mammals and demonstrate a unique DNA damage response in primary cells of the long-lived M. lucifugus. We also find evidence of abundant adaptation in response to DNA viruses - but not RNA viruses - in Myotis and other bats in sharp contrast with other mammals, potentially contributing to the role of bats as reservoirs of zoonoses. Together, our results demonstrate how genomics and primary cells derived from diverse taxa uncover the molecular bases of extreme adaptations in non-model organisms.

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.

Widespread Adaptive Introgression of Major Histocompatibility Complex Genes across Vertebrate Hybrid Zones

Interspecific introgression is a potentially important source of novel variation of adaptive significance. Although multiple cases of adaptive introgression are well documented, broader generalizations about its targets and mechanisms are lacking. Multiallelic balancing selection, particularly when acting through rare allele advantage, is an evolutionary mechanism expected to favor adaptive introgression. This is because introgressed alleles are likely to confer an immediate selective advantage, facilitating their establishment in the recipient species even in the face of strong genomic barriers to introgression. Vertebrate major histocompatibility complex genes are well-established targets of long-term multiallelic balancing selection, so widespread adaptive major histocompatibility complex introgression is expected. Here, we evaluate this hypothesis using data from 29 hybrid zones formed by fish, amphibians, squamates, turtles, birds, and mammals at advanced stages of speciation. The key prediction of more extensive major histocompatibility complex introgression compared to genome-wide introgression was tested with three complementary statistical approaches. We found evidence for widespread adaptive introgression of major histocompatibility complex genes, providing a link between the process of adaptive introgression and an underlying mechanism. Our work identifies major histocompatibility complex introgression as a general mechanism by which species can acquire novel, and possibly regain previously lost, variation that may enhance defense against pathogens and increase adaptive potential.

Ongoing introgression of a secondary sexual plumage trait in a stable avian hybrid zone

Hybrid zones are dynamic systems where natural selection, sexual selection, and other evolutionary forces can act on reshuffled combinations of distinct genomes. The movement of hybrid zones, individual traits, or both are of particular interest for understanding the interplay between selective processes. In a hybrid zone involving two lek-breeding birds, secondary sexual plumage traits of Manacus vitellinus, including bright yellow collar and olive belly color, have introgressed ~50 km asymmetrically across the genomic center of the zone into populations more genetically similar to Manacus candei. Males with yellow collars are preferred by females and are more aggressive than parental M. candei, suggesting that sexual selection was responsible for the introgression of male traits. We assessed the spatial and temporal dynamics of this hybrid zone using historical (1989-1994) and contemporary (2017-2020) transect samples to survey both morphological and genetic variation. Genome-wide single nucleotide polymorphism data and several male phenotypic traits show that the genomic center of the zone has remained spatially stable, whereas the olive belly color of male M. vitellinus has continued to introgress over this time period. Our data suggest that sexual selection can continue to shape phenotypes dynamically, independent of a stable genomic transition between species.

Genetic Signatures of Positive Selection in Human Populations Adapted to High Altitude in Papua New Guinea

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

Local genetic adaptation to habitat in wild chimpanzees

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.

Critically-ill COVID-19 susceptibility gene CCR3 shows natural selection in sub-Saharan Africans

The prevalence of COVID-19 critical illness varies across ethnicities, with recent studies suggesting that genetic factors may contribute to this variation. The aim of this study was to investigate natural selection signals of genes associated with critically-ill COVID-19 in sub-Saharan Africans. Severe COVID-19 SNPs were obtained from the HGI website. Selection signals were assessed in 661 sub-Sahara Africans from 1000 Genomes Project using integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and fixation index (Fst). Allele frequency trajectory analysis of ancient DNA samples were used to validate the existing of selection in sub-Sahara Africans. We also used Mendelian randomization to decipher the correlation between natural selection and critically-ill COVID-19. We identified that CCR3 exhibited significant natural selection signals in sub-Sahara Africans. Within the CCR3 gene, rs17217831-A showed both high iHS (Standardized iHS = 2) and high XP-EHH (Standardized XP-EHH = 2.5) in sub-Sahara Africans. Allele frequency trajectory of CCR3 rs17217831-A revealed natural selection occurring in the recent 1,500 years. Natural selection resulted in increased CCR3 expression in sub-Sahara Africans. Mendelian Randomization provided evidence that increased blood CCR3 expression and eosinophil counts lowered the risk of critically ill COVID-19. Our findings suggest that sub-Saharan Africans are resistant to critically ill COVID-19 due to natural selection and identify CCR3 as a potential novel therapeutic target.

Facilitated introgression from domestic goat into Alpine ibex at immune loci

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1-2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

Archaeology meets environmental genomics: implementing sedaDNA in the study of the human past

Sedimentary ancient DNA (sedaDNA) has become one of the standard applications in the field of paleogenomics in recent years. It has been used for paleoenvironmental reconstructions, detecting the presence of prehistoric species in the absence of macro remains and even investigating the evolutionary history of a few species. However, its application in archaeology has been limited and primarily focused on humans. This article argues that sedaDNA holds significant potential in addressing key archaeological questions concerning the origins, lifestyles, and environments of past human populations. Our aim is to facilitate the integration of sedaDNA into the standard workflows in archaeology as a transformative tool, thereby unleashing its full potential for studying the human past. Ultimately, we not only underscore the challenges inherent in the sedaDNA field but also provide a research agenda for essential enhancements needed for implementing sedaDNA into the archaeological workflow.

G-quadruplex propensity in H. neanderthalensis, H. sapiens and Denisovans mitochondrial genomes

Current methods of processing archaeological samples combined with advances in sequencing methods lead to disclosure of a large part of H. neanderthalensis and Denisovans genetic information. It is hardly surprising that the genome variability between modern humans, Denisovans and H. neanderthalensis is relatively limited. Genomic studies may provide insight on the metabolism of extinct human species or lineages. Detailed analysis of G-quadruplex sequences in H. neanderthalensis and Denisovans mitochondrial DNA showed us interesting features. Relatively similar patterns in mitochondrial DNA are found compared to modern humans, with one notable exception for H. neanderthalensis. An interesting difference between H. neanderthalensis and H. sapiens corresponds to a motif found in the D-loop region of mtDNA, which is responsible for mitochondrial DNA replication. This area is directly responsible for the number of mitochondria and consequently for the efficient energy metabolism of cell. H. neanderthalensis harbor a long uninterrupted run of guanines in this region, which may cause problems for replication, in contrast with H. sapiens, for which this run is generally shorter and interrupted. One may propose that the predominant H. sapiens motif provided a selective advantage for modern humans regarding mtDNA replication and function.

Reconstructing Prehistoric Viral Genomes from Neanderthal Sequencing Data

DNA viruses that produce persistent infections have been proposed as potential causes for the extinction of Neanderthals, and, therefore, the identification of viral genome remnants in Neanderthal sequence reads is an initial step to address this hypothesis. Here, as proof of concept, we searched for viral remnants in sequence reads of Neanderthal genome data by mapping to adenovirus, herpesvirus and papillomavirus, which are double-stranded DNA viruses that may establish lifelong latency and can produce persistent infections. The reconstructed ancient viral genomes of adenovirus, herpesvirus and papillomavirus revealed conserved segments, with nucleotide identity to extant viral genomes and variable regions in coding regions with substantial divergence to extant close relatives. Sequence reads mapped to extant viral genomes showed deamination patterns of ancient DNA, and these ancient viral genomes showed divergence consistent with the age of these samples (≈50,000 years) and viral evolutionary rates (10-5 to 10-8 substitutions/site/year). Analysis of random effects showed that the Neanderthal mapping to genomes of extant persistent viruses is above what is expected by random similarities of short reads. Also, negative control with a nonpersistent DNA virus does not yield statistically significant assemblies. This work demonstrates the feasibility of identifying viral genome remnants in archaeological samples with signal-to-noise assessment.

Shaping immunity: The influence of natural selection on population immune diversity

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Inference of selective forces on house mouse genomes during secondary contact in East Asia

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

An efficient and robust ABC approach to infer the rate and strength of adaptation

Inferring the effects of positive selection on genomes remains a critical step in characterizing the ultimate and proximate causes of adaptation across species, and quantifying positive selection remains a challenge due to the confounding effects of many other evolutionary processes. Robust and efficient approaches for adaptation inference could help characterize the rate and strength of adaptation in nonmodel species for which demographic history, mutational processes, and recombination patterns are not currently well-described. Here, we introduce an efficient and user-friendly extension of the McDonald-Kreitman test (ABC-MK) for quantifying long-term protein adaptation in specific lineages of interest. We characterize the performance of our approach with forward simulations and find that it is robust to many demographic perturbations and positive selection configurations, demonstrating its suitability for applications to nonmodel genomes. We apply ABC-MK to the human proteome and a set of known virus interacting proteins (VIPs) to test the long-term adaptation in genes interacting with viruses. We find substantially stronger signatures of positive selection on RNA-VIPs than DNA-VIPs, suggesting that RNA viruses may be an important driver of human adaptation over deep evolutionary time scales.

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.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

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.

Genetic history of Cambridgeshire before and after the Black Death

The extent of the devastation of the Black Death pandemic (1346-1353) on European populations is known from documentary sources and its bacterial source illuminated by studies of ancient pathogen DNA. What has remained less understood is the effect of the pandemic on human mobility and genetic diversity at the local scale. Here, we report 275 ancient genomes, including 109 with coverage >0.1×, from later medieval and postmedieval Cambridgeshire of individuals buried before and after the Black Death. Consistent with the function of the institutions, we found a lack of close relatives among the friars and the inmates of the hospital in contrast to their abundance in general urban and rural parish communities. While we detect long-term shifts in local genetic ancestry in Cambridgeshire, we find no evidence of major changes in genetic ancestry nor higher differentiation of immune loci between cohorts living before and after the Black Death.

Evolutionary immuno-genetics of endoplasmic reticulum aminopeptidase II (ERAP2)

Endoplasmic reticulum aminopeptidase 2 (ERAP2) is a proteolytic enzyme involved in adaptive immunity. The ERAP2 gene is highly polymorphic and encodes haplotypes that confer resistance against lethal infectious diseases, but also increase the risk for autoimmune disorders. Identifying how ERAP2 influences susceptibility to these traits requires an understanding of the selective pressures that shaped and maintained allelic variation throughout human evolution. Our review discusses the genetic regulation of haplotypes and diversity in naturally occurring ERAP2 allotypes in the global population. We outline how these ERAP2 haplotypes evolved during human history and highlight the presence of Neanderthal DNA sequences in ERAP2 of modern humans. Recent evidence suggests that human adaptation during the last ~10,000 years and historic pandemics left a significant mark on the ERAP2 gene that determines susceptibility to infectious and inflammatory diseases today.

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.

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 evolution: Neanderthal footprints in African genomes

Human and Neanderthal populations met and mixed on multiple occasions over evolutionary time, resulting in the exchange of genetic material. New genomic analyses of diverse African populations reveal a history of bidirectional gene flow and selection acting on introgressed alleles.

Assessing the Presence of Recent Adaptation in the Human Genome With Mixture Density Regression

How much genome differences between species reflect neutral or adaptive evolution is a central question in evolutionary genomics. In humans and other mammals, the presence of adaptive versus neutral genomic evolution has proven particularly difficult to quantify. The difficulty notably stems from the highly heterogeneous organization of mammalian genomes at multiple levels (functional sequence density, recombination, etc.) which complicates the interpretation and distinction of adaptive versus neutral evolution signals. In this study, we introduce mixture density regressions (MDRs) for the study of the determinants of recent adaptation in the human genome. MDRs provide a flexible regression model based on multiple Gaussian distributions. We use MDRs to model the association between recent selection signals and multiple genomic factors likely to affect the occurrence/detection of positive selection, if the latter was present in the first place to generate these associations. We find that an MDR model with two Gaussian distributions provides an excellent fit to the genome-wide distribution of a common sweep summary statistic (integrated haplotype score), with one of the two distributions likely enriched in positive selection. We further find several factors associated with signals of recent adaptation, including the recombination rate, the density of regulatory elements in immune cells, GC content, gene expression in immune cells, the density of mammal-wide conserved elements, and the distance to the nearest virus-interacting gene. These results support the presence of strong positive selection in recent human evolution and highlight MDRs as a powerful tool to make sense of signals of recent genomic adaptation.

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.

An efficient and robust ABC approach to infer the rate and strength of adaptation

Inferring the effects of positive selection on genomes remains a critical step in characterizing the ultimate and proximate causes of adaptation across species, and quantifying positive selection remains a challenge due to the confounding effects of many other evolutionary processes. Robust and efficient approaches for adaptation inference could help characterize the rate and strength of adaptation in non-model species for which demographic history, mutational processes, and recombination patterns are not currently well-described. Here, we introduce an efficient and user-friendly extension of the McDonald-Kreitman test (ABC-MK) for quantifying long-term protein adaptation in specific lineages of interest. We characterize the performance of our approach with forward simulations and find that it is robust to many demographic perturbations and positive selection configurations, demonstrating its suitability for applications to non-model genomes. We apply ABC-MK to the human proteome and a set of known Virus Interacting Proteins (VIPs) to test the long-term adaptation in genes interacting with viruses. We find substantially stronger signatures of positive selection on RNA-VIPs than DNA-VIPs, suggesting that RNA viruses may be an important driver of human adaptation over deep evolutionary time scales.

Dissecting human population variation in single-cell responses to SARS-CoV-2

Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.

The role of genetic selection and climatic factors in the dispersal of anatomically modern humans out of Africa

The evolutionarily recent dispersal of anatomically modern humans (AMH) out of Africa (OoA) and across Eurasia provides a unique opportunity to examine the impacts of genetic selection as humans adapted to multiple new environments. Analysis of ancient Eurasian genomic datasets (~1,000 to 45,000 y old) reveals signatures of strong selection, including at least 57 hard sweeps after the initial AMH movement OoA, which have been obscured in modern populations by extensive admixture during the Holocene. The spatiotemporal patterns of these hard sweeps provide a means to reconstruct early AMH population dispersals OoA. We identify a previously unsuspected extended period of genetic adaptation lasting ~30,000 y, potentially in the Arabian Peninsula area, prior to a major Neandertal genetic introgression and subsequent rapid dispersal across Eurasia as far as Australia. Consistent functional targets of selection initiated during this period, which we term the Arabian Standstill, include loci involved in the regulation of fat storage, neural development, skin physiology, and cilia function. Similar adaptive signatures are also evident in introgressed archaic hominin loci and modern Arctic human groups, and we suggest that this signal represents selection for cold adaptation. Surprisingly, many of the candidate selected loci across these groups appear to directly interact and coordinately regulate biological processes, with a number associated with major modern diseases including the ciliopathies, metabolic syndrome, and neurodegenerative disorders. This expands the potential for ancestral human adaptation to directly impact modern diseases, providing a platform for evolutionary medicine.

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 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.

Possible Causes of Hypertrophic Osteoarthropathy in the La Ferrassie 1 Neanderthal

For over a century, researchers have been perplexed by the unique osteological findings on La Ferrassie 1 (LF1), one of the most complete Neanderthal remains ever found. In 1997, Fennel and Trinkaus proposed that LF1 suffered from hypertrophic osteoarthropathy (HOA), likely secondary to chronic thoracic infection or pulmonary malignancy. This disease process can have many etiologies, and no study has fully explored the possible origin of LF1's HOA. Ultimately, it is most likely that LF1's HOA etiology arose from one of the many infectious diseases that prehistoric Neanderthals were exposed to, specifically a chronic pulmonary RNA virus.

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

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.

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.

Antiviral function and viral antagonism of the rapidly evolving dynein activating adaptor NINL

Viruses interact with the intracellular transport machinery to promote viral replication. Such host-virus interactions can drive host gene adaptation, leaving signatures of pathogen-driven evolution in host genomes. Here, we leverage these genetic signatures to identify the dynein activating adaptor, ninein-like (NINL), as a critical component in the antiviral innate immune response and as a target of viral antagonism. Unique among genes encoding components of active dynein complexes, NINL has evolved under recurrent positive (diversifying) selection, particularly in its carboxy-terminal cargo-binding region. Consistent with a role for NINL in host immunity, we demonstrate that NINL knockout cells exhibit an impaired response to interferon, resulting in increased permissiveness to viral replication. Moreover, we show that proteases encoded by diverse picornaviruses and coronaviruses cleave and disrupt NINL function in a host- and virus-specific manner. Our work reveals the importance of NINL in the antiviral response and the utility of using signatures of host-virus genetic conflicts to uncover new components of antiviral immunity and targets of viral antagonism.

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.

Genetic adaptations to SIV across chimpanzee populations

Central and eastern chimpanzees are infected with Simian Immunodeficiency Virus (SIV) in the wild, typically without developing acute immunodeficiency. Yet the recent zoonotic transmission of chimpanzee SIV to humans, which were naïve to the virus, gave rise to the Human Immunodeficiency Virus (HIV), which causes AIDS and is responsible for one of the deadliest pandemics in human history. Chimpanzees have likely been infected with SIV for tens of thousands of years and have likely evolved to reduce its pathogenicity, becoming semi-natural hosts that largely tolerate the virus. In support of this view, central and eastern chimpanzees show evidence of positive selection in genes involved in SIV/HIV cell entry and immune response to SIV, respectively. We hypothesise that the population first infected by SIV would have experienced the strongest selective pressure to control the lethal potential of zoonotic SIV, and that population genetics will reveal those first critical adaptations. With that aim we used population genetics to investigate signatures of positive selection in the common ancestor of central-eastern chimpanzees. The genes with signatures of positive selection in the ancestral population are significantly enriched in SIV-related genes, especially those involved in the immune response to SIV and those encoding for host genes that physically interact with SIV/HIV (VIPs). This supports a scenario where SIV first infected the central-eastern ancestor and where this population was under strong pressure to adapt to zoonotic SIV. Interestingly, integrating these genes with candidates of positive selection in the two infected subspecies reveals novel patterns of adaptation to SIV. Specifically, we observe evidence of positive selection in numerous steps of the biological pathway responsible for T-helper cell differentiation, including CD4 and multiple genes that SIV/HIV use to infect and control host cells. This pathway is active only in CD4+ cells which SIV/HIV infects, and it plays a crucial role in shaping the immune response so it can efficiently control the virus. Our results confirm the importance of SIV as a selective factor, identify specific genetic changes that may have allowed our closest living relatives to reduce SIV’s pathogenicity, and demonstrate the potential of population genomics to reveal the evolutionary mechanisms used by naïve hosts to reduce the pathogenicity of zoonotic pathogens.

Chimpanzees are at the origin of HIV-1, a virus that generates an incurable disease and that generated a pandemic that has claimed 35 million lives. Chimpanzees have evolved to control the pathogenicity of the virus, which does not typically develop into AIDS in the same way as in humans. Identifying the genetic adaptations responsible for this process provides critical knowledge about SIV and HIV. Our analysis of chimpanzee genetic adaptations identified specific genes and molecular pathways involved in adaptation to SIV, providing important insights into the mechanisms that likely allowed our closest living relatives to control SIV/HIV. Further, we establish SIV as a strong and recurrent selective pressure in central and eastern chimpanzees, two important subspecies of large mammals that are currently endangered.

The genetic and evolutionary determinants of COVID-19 susceptibility

Devastating pandemics, such as that due to COVID-19, can provide strong testimony to our knowledge of the genetic and evolutionary determinants of infectious disease susceptibility and severity. One of the most remarkable aspects of such outbreaks is the stunning interindividual variability observed in the course of infection. In recent decades, enormous progress has been made in the field of the human genetics of infectious diseases, and an increasing number of human genetic factors have been reported to explain, to a great extent, the observed variability for a large number of infectious agents. However, our understanding of the cellular, molecular, and immunological mechanisms underlying such disparities between individuals and ethnic groups, remains very limited. Here, we discuss recent findings relating to human genetic predisposition to infectious disease, from an immunological or population genetic perspective, and show how these and other innovative approaches have been applied to deciphering the genetic basis of human susceptibility to COVID-19 and the severity of this disease. From an evolutionary perspective, we show how past demographic and selection events characterizing the history of our species, including admixture with archaic humans, such as Neanderthals, facilitated modern human adaptation to the threats imposed by ancient pathogens. In the context of emerging infectious diseases, these past episodes of genetic adaptation may contribute to some of the observed population differences in the outcome of SARS-CoV-2 infection and the severity of COVID-19 illness.

Host-Genome Similarity Characterizes the Adaption of SARS-CoV-2 to Humans

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a high mutation rate and many variants have emerged in the last 2 years, including Alpha, Beta, Delta, Gamma and Omicron. Studies showed that the host-genome similarity (HGS) of SARS-CoV-2 is higher than SARS-CoV and the HGS of open reading frame (ORF) in coronavirus genome is closely related to suppression of innate immunity. Many works have shown that ORF 6 and ORF 8 of SARS-CoV-2 play an important role in suppressing IFN-β signaling pathway in vivo. However, the relation between HGS and the adaption of SARS-CoV-2 variants is still not clear. This work investigates HGS of SARS-CoV-2 variants based on a dataset containing more than 40,000 viral genomes. The relation between HGS of viral ORFs and the suppression of antivirus response is studied. The results show that ORF 7b, ORF 6 and ORF 8 are the top 3 genes with the highest HGS. In the past 2 years, the HGS values of ORF 8 and ORF 7B of SARS-CoV-2 have increased greatly. A remarkable correlation is discovered between HGS and inhibition of antivirus response of immune system, which suggests that the similarity between coronavirus and host gnome may be an indicator of the suppression of innate immunity. Among the five variants (Alpha, Beta, Delta, Gamma and Omicron), Delta has the highest HGS and Omicron has the lowest HGS. This finding implies that the high HGS in Delta variant may indicate further suppression of host innate immunity. However, the relatively low HGS of Omicron is still a puzzle. By comparing the mutations in genomes of Alpha, Delta and Omicron variants, a commonly shared mutation ACT > ATT is identified in high-HGS strain populations. The high HGS mutations among the three variants are quite different. This finding strongly suggests that mutations in high HGS strains are different in different variants. Only a few common mutations survive, which may play important role in improving the adaptability of SARS-CoV-2. However, the mechanism for how the mutations help SARS-CoV-2 escape immunity is still unclear. HGS analysis is a new method to study virus−host interaction and may provide a way to understand the rapid mutation and adaption of SARS-CoV-2.

Warthog Genomes Resolve an Evolutionary Conundrum and Reveal Introgression of Disease Resistance Genes

African wild pigs have a contentious evolutionary and biogeographic history. Until recently, desert warthog (Phacochoerus aethiopicus) and common warthog (P. africanus) were considered a single species. Molecular evidence surprisingly suggested they diverged at least 4.4 million years ago, and possibly outside of Africa. We sequenced the first whole-genomes of four desert warthogs and 35 common warthogs from throughout their range. We show that these two species diverged much later than previously estimated, 400,000–1,700,000 years ago depending on assumptions of gene flow. This brings it into agreement with the paleontological record. We found that the common warthog originated in western Africa and subsequently colonized eastern and southern Africa. During this range expansion, the common warthog interbred with the desert warthog, presumably in eastern Africa, underlining this region’s importance in African biogeography. We found that immune system–related genes may have adaptively introgressed into common warthogs, indicating that resistance to novel diseases was one of the most potent drivers of evolution as common warthogs expanded their range. Hence, we solve some of the key controversies surrounding warthog evolution and reveal a complex evolutionary history involving range expansion, introgression, and adaptation to new diseases.

Paleogenomics: The demographic past of prehistoric Europeans

Ancient DNA provides answers to long-standing debates about past human history. New work using demographic modeling on ancient genomes documents the nature and timing of the demographic processes - population size changes, divergences and admixture - that took place in prehistoric Europe.