Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors

Pathogens and the diseases they cause have been among the most important selective forces experienced by humans during their evolutionary history. Although adaptive alleles generally arise by mutation, introgression can also be a valuable source of beneficial alleles. Archaic humans, who lived in Europe and Western Asia for more than 200,000 years, were probably well adapted to this environment and its local pathogens. It is therefore conceivable that modern humans entering Europe and Western Asia who admixed with them obtained a substantial immune advantage from the introgression of archaic alleles. Here we document a cluster of three Toll-like receptors (TLR6-TLR1-TLR10) in modern humans that carries three distinct archaic haplotypes, indicating repeated introgression from archaic humans. Two of these haplotypes are most similar to the Neandertal genome, and the third haplotype is most similar to the Denisovan genome. The Toll-like receptors are key components of innate immunity and provide an important first line of immune defense against bacteria, fungi, and parasites. The unusually high allele frequencies and unexpected levels of population differentiation indicate that there has been local positive selection on multiple haplotypes at this locus. We show that the introgressed alleles have clear functional effects in modern humans; archaic-like alleles underlie differences in the expression of the TLR genes and are associated with reduced microbial resistance and increased allergic disease in large cohorts. This provides strong evidence for recurrent adaptive introgression at the TLR6-TLR1-TLR10 locus, resulting in differences in disease phenotypes in modern humans.     

Evolutionary Dynamics of Human Toll-Like Receptors and Their Different Contributions to Host Defense

Infectious diseases have been paramount among the threats to health and survival throughout human evolutionary history. Natural selection is therefore expected to act strongly on host defense genes, particularly on innate immunity genes whose products mediate the direct interaction between the host and the microbial environment. In insects and mammals, the Toll-like receptors (TLRs) appear to play a major role in initiating innate immune responses against microbes. In humans, however, it has been speculated that the set of TLRs could be redundant for protective immunity. We investigated how natural selection has acted upon human TLRs, as an approach to assess their level of biological redundancy. We sequenced the ten human TLRs in a panel of 158 individuals from various populations worldwide and found that the intracellular TLRs—activated by nucleic acids and particularly specialized in viral recognition—have evolved under strong purifying selection, indicating their essential non-redundant role in host survival. Conversely, the selective constraints on the TLRs expressed on the cell surface—activated by compounds other than nucleic acids—have been much more relaxed, with higher rates of damaging nonsynonymous and stop mutations tolerated, suggesting their higher redundancy. Finally, we tested whether TLRs have experienced spatially-varying selection in human populations and found that the region encompassing TLR10-TLR1-TLR6 has been the target of recent positive selection among non-Africans. Our findings indicate that the different TLRs differ in their immunological redundancy, reflecting their distinct contributions to host defense. The insights gained in this study foster new hypotheses to be tested in clinical and epidemiological genetics of infectious disease.     

Specific Gene Expression Responses to Parasite Genotypes Reveal Redundancy of Innate Immunity in Vertebrates

Vertebrate innate immunity is the first line of defense against an invading pathogen and has long been assumed to be largely unspecific with respect to parasite/pathogen species. However, recent phenotypic evidence suggests that immunogenetic variation, i.e. allelic variability in genes associated with the immune system, results in host-parasite genotype-by-genotype interactions and thus specific innate immune responses. Immunogenetic variation is common in all vertebrate taxa and this reflects an effective immunological function in complex environments. However, the underlying variability in host gene expression patterns as response of innate immunity to within-species genetic diversity of macroparasites in vertebrates is unknown. We hypothesized that intra-specific variation among parasite genotypes must be reflected in host gene expression patterns. Here we used high-throughput RNA-sequencing to examine the effect of parasite genotypes on gene expression patterns of a vertebrate host, the three-spined stickleback (Gasterosteus aculeatus). By infecting naïve fish with distinct trematode genotypes of the species Diplostomum pseudospathaceum we show that gene activity of innate immunity in three-spined sticklebacks depended on the identity of an infecting macroparasite genotype. In addition to a suite of genes indicative for a general response against the trematode we also find parasite-strain specific gene expression, in particular in the complement system genes, despite similar infection rates of single clone treatments. The observed discrepancy between infection rates and gene expression indicates the presence of alternative pathways which execute similar functions. This suggests that the innate immune system can induce redundant responses specific to parasite genotypes.     

A Wild C. Elegans Strain Has Enhanced Epithelial Immunity to a Natural Microsporidian Parasite

Microbial pathogens impose selective pressures on their hosts, and combatting these pathogens is fundamental to the propagation of a species. Innate immunity is an ancient system that provides the foundation for pathogen resistance, with epithelial cells in humans increasingly appreciated to play key roles in innate defense. Here, we show that the nematode C. elegans displays genetic variation in epithelial immunity against intestinal infection by its natural pathogen, Nematocida parisii. This pathogen belongs to the microsporidia phylum, which comprises a large phylum of over 1400 species of fungal-related parasites that can infect all animals, including humans, but are poorly understood. Strikingly, we find that a wild C. elegans strain from Hawaii is able to clear intracellular infection by N. parisii, with this ability restricted to young larval animals. Notably, infection of older larvae does not impair progeny production, while infection of younger larvae does. The early-life immunity of Hawaiian larvae enables them to produce more progeny later in life, providing a selective advantage in a laboratory setting—in the presence of parasite it is able to out-compete a susceptible strain in just a few generations. We show that enhanced immunity is dominant to susceptibility, and we use quantitative trait locus mapping to identify four genomic loci associated with resistance. Furthermore, we generate near-isogenic strains to directly demonstrate that two of these loci influence resistance. Thus, our findings show that early-life immunity of C. elegans against microsporidia is a complex trait that enables the host to produce more progeny later in life, likely improving its evolutionary success.     

Sexual antagonism for resistance and tolerance to infection in Drosophila melanogaster

A critical task in evolutionary genetics is to explain the persistence of heritable variation in fitness-related traits such as immunity. Ecological factors can maintain genetic variation in immunity, but less is known about the role of other factors, such as antagonistic pleiotropy, on immunity. Sexually dimorphic immunity—with females often being more immune-competent—may maintain variation in immunity in dioecious populations. Most eco-immunological studies assess host resistance to parasites rather than the host''s ability to maintain fitness during infection (tolerance). Distinguishing between resistance and tolerance is important as they are thought to have markedly different evolutionary and epidemiological outcomes. Few studies have investigated tolerance in animals, and the extent of sexual dimorphism in tolerance is unknown. Using males and females from 50 Drosophila melanogaster genotypes, we investigated possible sources of genetic variation for immunity by assessing both resistance and tolerance to the common bacterial pathogen Pseudomonas aeruginosa. We found evidence of sexual dimorphism and sexual antagonism for resistance and tolerance, and a trade-off between the two traits. Our findings suggest that antagonistic pleiotropy may be a major contributor to variation in immunity, with implications for host–parasite coevolution.     

Searching for signatures of cold adaptations in modern and archaic humans: hints from the brown adipose tissue genes

Adaptation to low temperatures has been reasonably developed in the human species during the colonization of the Eurasian landmass subsequent to Out of Africa migrations of anatomically modern humans. In addition to morphological and cultural changes, also metabolic ones are supposed to have favored human isolation from cold and body heat production and this can be hypothesized also for most Neandertal and at least for some Denisovan populations, which lived in geographical areas that strongly experienced the last glacial period. Modulation of non-shivering thermogenesis, for which adipocytes belonging to the brown adipose tissue are the most specialized cells, might have driven these metabolic adaptations. To perform an exploratory analysis aimed at looking into this hypothesis, variation at 28 genes involved in such functional pathway was investigated in modern populations from different climate zones, as well as in Neandertal and Denisovan genomes. Patterns of variation at the LEPR gene, strongly related to increased heat dissipation by mitochondria, appeared to have been shaped by positive selection in modern East Asians, but not in Europeans. Moreover, a single potentially cold-adapted LEPR allele, different from the supposed adaptive one identified in Homo sapiens, was found also in Neandertal and Denisovan genomes. These findings suggest that independent mechanisms for cold adaptations might have been developed in different non-African human groups, as well as that the evolution of possible enhanced thermal efficiency in Neandertals and in some Denisovan populations has plausibly entailed significant changes also in other functional pathways than in the examined one.     

Enhanced Monocyte Response and Decreased Central Memory T Cells in Children with Invasive Staphylococcus aureus Infections

Staphylococcus aureus has emerged as a significant pathogen causing severe invasive disease in otherwise healthy people. Despite considerable advances in understanding the epidemiology, resistance mechanisms, and virulence factors produced by the bacteria, there is limited knowledge of the in vivo host immune response to acute, invasive S. aureus infections. Herein, we report that peripheral blood mononuclear cells from patients with severe S. aureus infections demonstrate a distinctive and robust gene expression profile which is validated in a distinct group of patients and on a different microarray platform. Application of a systems-wide modular analysis framework reveals significant over-expression of innate immunity genes and under-expression of genes related to adaptive immunity. Simultaneous flow cytometry analyses demonstrated marked alterations in immune cell numbers, with decreased central memory CD4 and CD8 T cells and increased numbers of monocytes. CD14+ monocyte numbers significantly correlated with the gene expression levels of genes related to the innate immune response. These results demonstrate the value of applying a systems biology approach that reveals the significant alterations in the components of circulating blood lymphocytes and monocytes in invasive S. aureus infections.     

Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments—Haplochromis chilotes, H. sauvagei, and Lithochromis rufus—to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation.     

Body temperature variation controls pre-mRNA processing and transcription of antiviral genes and SARS-CoV-2 replication

Antiviral innate immunity represents the first defense against invading viruses and is key to control viral infections, including SARS-CoV-2. Body temperature is an omnipresent variable but was neglected when addressing host defense mechanisms and susceptibility to SARS-CoV-2 infection. Here, we show that increasing temperature in a 1.5°C window, between 36.5 and 38°C, strongly increases the expression of genes in two branches of antiviral immunity, nitric oxide production and type I interferon response. We show that alternative splicing coupled to nonsense-mediated decay decreases STAT2 expression in colder conditions and suggest that increased STAT2 expression at elevated temperature induces the expression of diverse antiviral genes and SARS-CoV-2 restriction factors. This cascade is activated in a remarkably narrow temperature range below febrile temperature, which reflects individual, circadian and age-dependent variation. We suggest that decreased body temperature with aging contributes to reduced expression of antiviral genes in older individuals. Using cell culture and in vivo models, we show that higher body temperature correlates with reduced SARS-CoV-2 replication, which may affect the different vulnerability of children versus seniors toward severe SARS-CoV-2 infection. Altogether, our data connect body temperature and pre-mRNA processing to provide new mechanistic insight into the regulation of antiviral innate immunity.     

The Genetic Basis of Escherichia coli Pathoadaptation to Macrophages

Antagonistic interactions are likely important driving forces of the evolutionary process underlying bacterial genome complexity and diversity. We hypothesized that the ability of evolved bacteria to escape specific components of host innate immunity, such as phagocytosis and killing by macrophages (MΦ), is a critical trait relevant in the acquisition of bacterial virulence. Here, we used a combination of experimental evolution, phenotypic characterization, genome sequencing and mathematical modeling to address how fast, and through how many adaptive steps, a commensal Escherichia coli (E. coli) acquire this virulence trait. We show that when maintained in vitro under the selective pressure of host MΦ commensal E. coli can evolve, in less than 500 generations, virulent clones that escape phagocytosis and MΦ killing in vitro, while increasing their pathogenicity in vivo, as assessed in mice. This pathoadaptive process is driven by a mechanism involving the insertion of a single transposable element into the promoter region of the E. coli yrfF gene. Moreover, transposition of the IS186 element into the promoter of Lon gene, encoding an ATP-dependent serine protease, is likely to accelerate this pathoadaptive process. Competition between clones carrying distinct beneficial mutations dominates the dynamics of the pathoadaptive process, as suggested from a mathematical model, which reproduces the observed experimental dynamics of E. coli evolution towards virulence. In conclusion, we reveal a molecular mechanism explaining how a specific component of host innate immunity can modulate microbial evolution towards pathogenicity.     

Predicting distributions of Wolbachia strains through host ecological contact—Who's manipulating whom?

Reproductive isolation in response to divergent selection is often mediated via third‐party interactions. Under these conditions, speciation is inextricably linked to ecological context. We present a novel framework for understanding arthropod speciation as mediated by Wolbachia, a microbial endosymbiont capable of causing host cytoplasmic incompatibility (CI). We predict that sympatric host sister‐species harbor paraphyletic Wolbachia strains that provide CI, while well‐defined congeners in ecological contact and recently diverged noninteracting congeners are uninfected due to Wolbachia redundancy. We argue that Wolbachia provides an adaptive advantage when coupled with reduced hybrid fitness, facilitating assortative mating between co‐occurring divergent phenotypes—the contact contingency hypothesis. To test this, we applied a predictive algorithm to empirical pollinating fig wasp data, achieving up to 91.60% accuracy. We further postulate that observed temporal decay of Wolbachia incidence results from adaptive host purging—adaptive decay hypothesis—but implementation failed to predict systematic patterns. We then account for post‐zygotic offspring mortality during CI mating, modeling fitness clines across developmental resources—the fecundity trade‐off hypothesis. This model regularly favored CI despite fecundity losses. We demonstrate that a rules‐based algorithm accurately predicts Wolbachia infection status. This has implications among other systems where closely related sympatric species encounter adaptive disadvantage through hybridization.     

Three Prochlorococcus cyanophage genomes: signature features and ecological interpretations

The oceanic cyanobacteria Prochlorococcus are globally important, ecologically diverse primary producers. It is thought that their viruses (phages) mediate population sizes and affect the evolutionary trajectories of their hosts. Here we present an analysis of genomes from three Prochlorococcus phages: a podovirus and two myoviruses. The morphology, overall genome features, and gene content of these phages suggest that they are quite similar to T7-like (P-SSP7) and T4-like (P-SSM2 and P-SSM4) phages. Using the existing phage taxonomic framework as a guideline, we examined genome sequences to establish “core” genes for each phage group. We found the podovirus contained 15 of 26 core T7-like genes and the two myoviruses contained 43 and 42 of 75 core T4-like genes. In addition to these core genes, each genome contains a significant number of “cyanobacterial” genes, i.e., genes with significant best BLAST hits to genes found in cyanobacteria. Some of these, we speculate, represent “signature” cyanophage genes. For example, all three phage genomes contain photosynthetic genes (psbA, hliP) that are thought to help maintain host photosynthetic activity during infection, as well as an aldolase family gene (talC) that could facilitate alternative routes of carbon metabolism during infection. The podovirus genome also contains an integrase gene (int) and other features that suggest it is capable of integrating into its host. If indeed it is, this would be unprecedented among cultured T7-like phages or marine cyanophages and would have significant evolutionary and ecological implications for phage and host. Further, both myoviruses contain phosphate-inducible genes (phoH and pstS) that are likely to be important for phage and host responses to phosphate stress, a commonly limiting nutrient in marine systems. Thus, these marine cyanophages appear to be variations of two well-known phages—T7 and T4—but contain genes that, if functional, reflect adaptations for infection of photosynthetic hosts in low-nutrient oceanic environments.     

Expression Profiling the Temperature-Dependent Amphibian Response to Infection by Batrachochytrium dendrobatidis

Amphibians are experiencing a panzootic of unprecedented proportions caused by the emergence of Batrachochytrium dendrobatidis (Bd). However, all species are not equally at risk of infection, and risk is further modified by environmental variables, specifically temperature. In order to understand how, and when, hosts mount a response to Bd we analysed infection dynamics and patterns of gene expression in the model amphibian species Silurana (Xenopus) tropicalis. Mathematical modelling of infection dynamics demonstrate the existence of a temperature-dependent protective response that is largely independent of the intrinsic growth-rate of Bd. Using temporal expression-profiling by microarrays and qRT-PCR, we characterise this response in the main amphibian lymphoid tissue, the spleen. We demonstrate that clearance of Bd at the host-optimal temperature is not clearly associated with an adaptive immune response, but rather is correlated with the induction of components of host innate immunity including the expression of genes that are associated with the production of the antimicrobial skin peptide preprocareulein (PPCP) as well as inflammatory responses. We find that adaptive immunity appears to be lacking at host-optimal temperatures. This suggests that either Bd does not stimulate, or suppresses, adaptive immunity, or that trade-offs exist between innate and adaptive limbs of the amphibian immune system. At cold temperatures, S. tropicalis loses the ability to mount a PPCP-based innate response, and instead manifests a more pronounced inflammatory reaction that is characterised by the production of proteases and higher pathogen burdens. This study demonstrates the temperature-dependency of the amphibian response to infection by Bd and indicates the influence that changing climates may exert on the ectothermic host response to pathogens.     

Controlled Chaos of Polymorphic Mucins in a Metazoan Parasite (Schistosoma mansoni) Interacting with Its Invertebrate Host (Biomphalaria glabrata)

Invertebrates were long thought to possess only a simple, effective and hence non-adaptive defence system against microbial and parasitic attacks. However, recent studies have shown that invertebrate immunity also relies on immune receptors that diversify (e.g. in echinoderms, insects and mollusks (Biomphalaria glabrata)). Apparently, individual or population-based polymorphism-generating mechanisms exists that permit the survival of invertebrate species exposed to parasites. Consequently, the generally accepted arms race hypothesis predicts that molecular diversity and polymorphism also exist in parasites of invertebrates. We investigated the diversity and polymorphism of parasite molecules (Schistosoma mansoni Polymorphic Mucins, SmPoMucs) that are key factors for the compatibility of schistosomes interacting with their host, the mollusc Biomphalaria glabrata. We have elucidated the complex cascade of mechanisms acting both at the genomic level and during expression that confer polymorphism to SmPoMuc. We show that SmPoMuc is coded by a multi-gene family whose members frequently recombine. We show that these genes are transcribed in an individual-specific manner, and that for each gene, multiple splice variants exist. Finally, we reveal the impact of this polymorphism on the SmPoMuc glycosylation status. Our data support the view that S. mansoni has evolved a complex hierarchical system that efficiently generates a high degree of polymorphism—a “controlled chaos”—based on a relatively low number of genes. This contrasts with protozoan parasites that generate antigenic variation from large sets of genes such as Trypanosoma cruzi, Trypanosoma brucei and Plasmodium falciparum. Our data support the view that the interaction between parasites and their invertebrate hosts are far more complex than previously thought. While most studies in this matter have focused on invertebrate host diversification, we clearly show that diversifying mechanisms also exist on the parasite side of the interaction. Our findings shed new light on how and why invertebrate immunity develops.     

Incipient Balancing Selection through Adaptive Loss of Aquaporins in Natural Saccharomyces cerevisiae Populations

A major goal in evolutionary biology is to understand how adaptive evolution has influenced natural variation, but identifying loci subject to positive selection has been a challenge. Here we present the adaptive loss of a pair of paralogous genes in specific Saccharomyces cerevisiae subpopulations. We mapped natural variation in freeze-thaw tolerance to two water transporters, AQY1 and AQY2, previously implicated in freeze-thaw survival. However, whereas freeze-thaw–tolerant strains harbor functional aquaporin genes, the set of sensitive strains lost aquaporin function at least 6 independent times. Several genomic signatures at AQY1 and/or AQY2 reveal low variation surrounding these loci within strains of the same haplotype, but high variation between strain groups. This is consistent with recent adaptive loss of aquaporins in subgroups of strains, leading to incipient balancing selection. We show that, although aquaporins are critical for surviving freeze-thaw stress, loss of both genes provides a major fitness advantage on high-sugar substrates common to many strains'' natural niche. Strikingly, strains with non-functional alleles have also lost the ancestral requirement for aquaporins during spore formation. Thus, the antagonistic effect of aquaporin function—providing an advantage in freeze-thaw tolerance but a fitness defect for growth in high-sugar environments—contributes to the maintenance of both functional and nonfunctional alleles in S. cerevisiae. This work also shows that gene loss through multiple missense and nonsense mutations, hallmarks of pseudogenization presumed to emerge after loss of constraint, can arise through positive selection.     

Exome sequencing of deer mice on two California Channel Islands identifies potential adaptation to strongly contrasting ecological conditions

Understanding the forces that drive genotypic and phenotypic change in wild populations is a central goal of evolutionary biology. We examined exome variation in populations of deer mice from two of the California Channel Islands: Peromyscus maniculatus elusus from Santa Barbara Island and P. m. santacruzae from Santa Cruz Island exhibit significant differences in olfactory predator recognition, activity timing, aggressive behavior, morphology, prevalence of Sin Nombre virus, and population densities. We characterized variation in protein‐coding regions using exome capture and sequencing of 25 mice from Santa Barbara Island and 22 mice from Santa Cruz Island. We identified and examined 386,256 SNPs using three complementary methods (BayeScan, pcadapt, and LFMM). We found strong differences in molecular variation between the two populations and 710 outlier SNPs in protein‐coding genes that were detected by all three methods. We identified 35 candidate genes from this outlier set that were related to differences in phenotypes between island populations. Enrichment analyses demonstrated that patterns of molecular variation were associated with biological processes related to response to chemical stimuli and regulation of immune processes. Candidate genes associated with olfaction (Gfy, Tlr2, Vmn13r2, numerous olfactory receptor genes), circadian activity (Cry1), anxiety (Brca1), immunity (Cd28, Eif2ak4, Il12a, Syne1), aggression (Cyp19a, Lama2), and body size (Bc16, Syne1) exhibited non‐synonymous mutations predicted to have moderate to large effects. Variation in olfaction‐related genes, including a stop codon in the Santa Barbara Island population, suggests loss of predator‐recognition traits at the molecular level, consistent with a lack of behavioral aversion to fox feces. These findings also suggest that divergent pathogen prevalence and population density may have influenced adaptive immunity and behavioral phenotypes, such as reduced aggression. Overall, our study indicates that ecological differences between islands are associated with signatures of selection in protein‐coding genes underlying phenotypes that promote success in those environments.