Different selective pressures shape the evolution of Toll-like receptors in human and African great ape populations

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

Humans with inherited MyD88 and IRAK-4 deficiencies are predisposed to hypoxemic COVID-19 pneumonia

X-linked recessive deficiency of TLR7, a MyD88- and IRAK-4-dependent endosomal ssRNA sensor, impairs SARS-CoV-2 recognition and type I IFN production in plasmacytoid dendritic cells (pDCs), thereby underlying hypoxemic COVID-19 pneumonia with high penetrance. We report 22 unvaccinated patients with autosomal recessive MyD88 or IRAK-4 deficiency infected with SARS-CoV-2 (mean age: 10.9 yr; 2 mo to 24 yr), originating from 17 kindreds from eight countries on three continents. 16 patients were hospitalized: six with moderate, four with severe, and six with critical pneumonia, one of whom died. The risk of hypoxemic pneumonia increased with age. The risk of invasive mechanical ventilation was also much greater than in age-matched controls from the general population (OR: 74.7, 95% CI: 26.8-207.8, P < 0.001). The patients' susceptibility to SARS-CoV-2 can be attributed to impaired TLR7-dependent type I IFN production by pDCs, which do not sense SARS-CoV-2 correctly. Patients with inherited MyD88 or IRAK-4 deficiency were long thought to be selectively vulnerable to pyogenic bacteria, but also have a high risk of hypoxemic COVID-19 pneumonia.

Unlocking life-threatening COVID-19 through two types of inborn errors of type I IFNs

Since 2003, rare inborn errors of human type I IFN immunity have been discovered, each underlying a few severe viral illnesses. Autoantibodies neutralizing type I IFNs due to rare inborn errors of autoimmune regulator (AIRE)-driven T cell tolerance were discovered in 2006, but not initially linked to any viral disease. These two lines of clinical investigation converged in 2020, with the discovery that inherited and/or autoimmune deficiencies of type I IFN immunity accounted for approximately 15%-20% of cases of critical COVID-19 pneumonia in unvaccinated individuals. Thus, insufficient type I IFN immunity at the onset of SARS-CoV-2 infection may be a general determinant of life-threatening COVID-19. These findings illustrate the unpredictable, but considerable, contribution of the study of rare human genetic diseases to basic biology and public health.

Evolution of toll-like receptor gene family in amphibians

The evolutionary position and lifestyle of amphibians highlights the important roles of the immune system in adaptive radiation and their adaptation to a complex pathogenic environment. Toll-like receptors (TLRs) are membrane-like sensors that recognize and bind conserved molecular motifs in pathogens to initiate downstream immune responses. To understand the evolutionary patterns of TLRs in amphibians, we analyzed TLR genes from the genomes and transcriptomes of 102 amphibian species. Phylogenetic results showed that 578 intact amphibian TLR sequences belonged to 16 TLR genes and were divided into seven subfamilies. The TLR4 subfamily was only identified in the Anura. Purification selection plays a leading role in amphibian TLR evolution and mean ω (dN/dS) values ranged from 0.252 for TLR7 to 0.381 for TLR19. Furthermore, the ω values of different domains were significantly different. We found positive selection patterns for 141 of 12,690 codons (1.1%) in all amphibian TLRs, most of which were located in leucine-rich repeats (LRRs). We also observed low to moderate levels of single-nucleotide polymorphisms (SNPs) in Pelophylax nigromaculatus and Bombina orientalis. This study provided critical primers, meaningful information regarding TLR gene family evolution in amphibians, and insights into the complex evolutionary patterns and implications of TLR polymorphisms.

Variations within Toll-like receptor (TLR) and TLR signalling pathway-related genes and their synergistic effects on the risk of Guillain-Barré Syndrome

Guillain-Barré Syndrome (GBS) is the commonest postinfectious polyradiculopathy. Though genetic background of the host seems to play an important role in the susceptibility to GBS, genes conferring major risk are not yet known. Dysregulation of Toll-like receptor (TLR) molecules exacerbates immune-inflammatory responses and the genetic variations within TLR pathway-related genes contribute to differential risk to infection. To delineate the impact of genetic variations within TLR2, TLR3, and TLR4 genes and TLR signaling pathway-related genes such as MyD88, TRIF, TRAF3, TRAF6, IRF3, NFκβ1, and IκBα on risk of developing GBS. Fourteen polymorphisms located within TLR2 (rs3804099; rs111200466), TLR3 (rs3775290; rs3775291), TLR4 (rs1927911, rs11536891), MyD88 (rs7744, rs4988453), TRIF (rs8120 TRAF3 (rs12147254), TRAF6 (rs4755453), IRF3 (rs2304204), NFκβ1 (rs28362491) and IκBα (rs696) genes were genotyped in 150 GBS patients and 150 healthy subjects either by PCR-RFLP or TaqMan Allelic Discrimination Assay. Genotypes of two polymorphic variants, Del/Del of rs111200466 Insertion and Deletion (INDEL) polymorphism of TLR2 gene and TT of rs3775290 single nucleotide polymorphism (SNP) of TLR3 gene had significantly higher frequencies among GBS patients, while the frequencies of TT genotype of rs3804099 of TLR2 gene and TT genotype of rs11536891 SNP of TLR4 gene were significantly higher in controls. Gene-gene interaction study by Multifactor Dimensionality Reduction (MDR) analysis also suggested a significant combined effect of TLR2, and NFκβ1 genes on the risk of GBS. The SNPs in the IκBα and IRF3 genes correlated with severity of GBS. The genes encoding TLRs and TLR signalling pathway-related molecules could serve as crucial genetic markers of susceptibility and severity of GBS. This article is protected by copyright. All rights reserved.

Evolutionary History of the Toll-Like Receptor Gene Family across Vertebrates

Adaptation to a wide range of pathogenic environments is a major aspect of the ecological adaptations of vertebrates during evolution. Toll-like receptors (TLRs) are ancient membrane-bound sensors in animals and are best known for their roles in detecting and defense against invading pathogenic microorganisms. To understand the evolutionary history of the vertebrate TLR gene family, we first traced the origin of single-cysteine cluster TLRs that share the same protein architecture with vertebrate TLRs in early-branching animals and then analyzed all members of the TLR family in over 200 species covering all major vertebrate clades. Our results indicate that although the emergence of single-cysteine cluster TLRs predates the separation of bilaterians and cnidarians, most vertebrate TLR members originated shortly after vertebrate emergence. Phylogenetic analyses divided 1,726 vertebrate TLRs into 8 subfamilies, and TLR3 may represent the most ancient subfamily that emerged before the branching of deuterostomes. Our analysis reveals that purifying selection predominated in the evolution of all vertebrate TLRs, with mean dN/dS (ω) values ranging from 0.082 for TLR21 in birds to 0.434 for TLR11 in mammals. However, we did observe patterns of positive selection acting on specific codons (527 of 60,294 codons across all vertebrate TLRs, 8.7‰), which are significantly concentrated in ligand-binding extracellular domains and suggest host–pathogen coevolutionary interactions. Additionally, we found stronger positive selection acting on nonviral compared with viral TLRs, indicating the more essential nonredundant function of viral TLRs in host immunity. Taken together, our findings provide comprehensive insight into the complex evolutionary processes of the vertebrate TLR gene family, involving gene duplication, pseudogenization, purification, and positive selection.

Presence of Genetic Variants Among Young Men With Severe COVID-19

IMPORTANCE

Severe coronavirus disease 2019 (COVID-19) can occur in younger, predominantly male, patients without preexisting medical conditions. Some individuals may have primary immunodeficiencies that predispose to severe infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).

OBJECTIVE

To explore the presence of genetic variants associated with primary immunodeficiencies among young patients with COVID-19.

DESIGN, SETTING, AND PARTICIPANTS

Case series of pairs of brothers without medical history meeting the selection criteria of young (age <35 years) brother pairs admitted to the intensive care unit (ICU) due to severe COVID-19. Four men from 2 unrelated families were admitted to the ICUs of 4 hospitals in the Netherlands between March 23 and April 12, 2020. The final date of follow-up was May 16, 2020. Available family members were included for genetic variant segregation analysis and as controls for functional experiments.

EXPOSURE

Severe COVID-19.

MAIN OUTCOME AND MEASURES

Results of rapid clinical whole-exome sequencing, performed to identify a potential monogenic cause. Subsequently, basic genetic and immunological tests were performed in primary immune cells isolated from the patients and family members to characterize any immune defects.

RESULTS

The 4 male patients had a mean age of 26 years (range, 21-32), with no history of major chronic disease. They were previously well before developing respiratory insufficiency due to severe COVID-19, requiring mechanical ventilation in the ICU. The mean duration of ventilatory support was 10 days (range, 9-11); the mean duration of ICU stay was 13 days (range, 10-16). One patient died. Rapid clinical whole-exome sequencing of the patients and segregation in available family members identified loss-of-function variants of the X-chromosomal TLR7. In members of family 1, a maternally inherited 4-nucleotide deletion was identified (c.2129_2132del; p.[Gln710Argfs*18]); the affected members of family 2 carried a missense variant (c.2383G>T; p.[Val795Phe]). In primary peripheral blood mononuclear cells from the patients, downstream type I interferon (IFN) signaling was transcriptionally downregulated, as measured by significantly decreased mRNA expression of IRF7, IFNB1, and ISG15 on stimulation with the TLR7 agonist imiquimod as compared with family members and controls. The production of IFN-γ, a type II IFN, was decreased in patients in response to stimulation with imiquimod.

CONCLUSIONS AND RELEVANCE

In this case series of 4 young male patients with severe COVID-19, rare putative loss-of-function variants of X-chromosomal TLR7 were identified that were associated with impaired type I and II IFN responses. These preliminary findings provide insights into the pathogenesis of COVID-19.

Comparing the selective landscape of TLR7 and TLR8 across primates reveals unique sites under positive selection in Alouatta

Among primates, susceptibility to yellow fever (YFV), a single-stranded (ss) RNA virus, ranges from complete resistance to high susceptibility. Howler monkeys (genus Alouatta) are the most susceptible to YFV. In order to identify Alouatta-specific genetic factors that may be responsible for their susceptibility, we collected skin samples from howler monkey museum specimens of the species A. caraya and A. guariba clamitans. We compared the rate of nonsynonymous to synonymous (dN/dS) changes of Toll-like receptor (TLR) 7 and TLR8, the two genes responsible for detecting all ssRNA viruses, across the Primate order. Overall, we found that the TLR7 gene is under stronger purifying selection in howler monkeys compared to other New World and Old World primates, but TLR8 is under the same selective pressure. When we evaluated dN/dS at each codon, we found six codons under positive selection in Alouatta TLR8 and two codons under positive selection in TLR7. The changes in TLR7 are unique to A. guariba clamitans and are found in functionally important regions likely to affect detection of ssRNA viruses by TLR7/TLR8, as well as downstream signaling. These amino acid differences in A. guariba clamitans may play a role in YFV susceptibility. These results have implications for identifying genetic factors affecting YFV susceptibility in primates.

Positively selected variants in functionally important regions of TLR7 in Alouatta guariba clamitans with yellow fever virus exposure in Northern Argentina

BACKGROUND

In 2007-2009, a major yellow fever virus (YFV) outbreak in Northern Argentina decimated the local howler monkey (Alouatta) population.

AIMS

To evaluate whether the surviving howler monkeys possess advantageous genetic variants inherited from monkeys alive prior to the YFV outbreak, we explored the relationship between Toll-like receptor (TLR) 7 and TLR8 gene variation and YFV susceptibility.

METHODS

We used samples from Alouatta individuals in Misiones, Argentina alive before the YFV outbreak, individuals that died during the outbreak, and individuals that survived the outbreak and are alive today. We measured genetic divergence between Alouatta YFV exposure groups and evaluated Alouatta-specific substitutions for functional consequences.

RESULTS

We did not find different allele frequencies in the post-YFV exposure Alouatta group compared to the pre-exposure group. We identified three nonsynonymous variants in TLR7 in Alouatta guariba clamitans. Two of these substitutions are under positive selection in functionally important regions of the gene.

DISCUSSION AND CONCLUSIONS

Our results did not indicate that surviving howler monkey spossess advantageous genetic variants at greater frequency than those alive before the YFV outbreak. However, the positively selected unique coding differences in A. guariba clamitans are in the region important in pathogen detection which may affect YFV resistance. Morework is necessary to fully explore this hypothesis.

Common homozygosity for predicted loss-of-function variants reveals both redundant and advantageous effects of dispensable human genes

Humans homozygous or hemizygous for variants predicted to cause a loss of function (LoF) of the corresponding protein do not necessarily present with overt clinical phenotypes. We report here 190 autosomal genes with 207 predicted LoF variants, for which the frequency of homozygous individuals exceeds 1% in at least one human population from five major ancestry groups. No such genes were identified on the X and Y chromosomes. Manual curation revealed that 28 variants (15%) had been misannotated as LoF. Of the 179 remaining variants in 166 genes, only 11 alleles in 11 genes had previously been confirmed experimentally to be LoF. The set of 166 dispensable genes was enriched in olfactory receptor genes (41 genes). The 41 dispensable olfactory receptor genes displayed a relaxation of selective constraints similar to that observed for other olfactory receptor genes. The 125 dispensable nonolfactory receptor genes also displayed a relaxation of selective constraints consistent with greater redundancy. Sixty-two of these 125 genes were found to be dispensable in at least three human populations, suggesting possible evolution toward pseudogenes. Of the 179 LoF variants, 68 could be tested for two neutrality statistics, and 8 displayed robust signals of positive selection. These latter variants included a known FUT2 variant that confers resistance to intestinal viruses, and an APOL3 variant involved in resistance to parasitic infections. Overall, the identification of 166 genes for which a sizeable proportion of humans are homozygous for predicted LoF alleles reveals both redundancies and advantages of such deficiencies for human survival.

Intrinsically disordered regions are abundant in simplexvirus proteomes and display signatures of positive selection

Whereas the majority of herpesviruses co-speciated with their mammalian hosts, human herpes simplex virus 2 (HSV-2, genus Simplexvirus) most likely originated from the cross-species transmission of chimpanzee herpesvirus 1 to an ancestor of modern humans. We exploited the peculiar evolutionary history of HSV-2 to investigate the selective events that drove herpesvirus adaptation to a new host. We show that HSV-2 intrinsically disordered regions (IDRs)-that is, protein domains that do not adopt compact three-dimensional structures-are strongly enriched in positive selection signals. Analysis of viral proteomes indicated that a significantly higher portion of simplexvirus proteins is disordered compared with the proteins of other human herpesviruses. IDR abundance in simplexvirus proteomes was not a consequence of the base composition of their genomes (high G + C content). Conversely, protein function determines the IDR fraction, which is significantly higher in viral proteins that interact with human factors. We also found that the average extent of disorder in herpesvirus proteins tends to parallel that of their human interactors. These data suggest that viruses that interact with fast-evolving, disordered human proteins, in turn, evolve disordered viral interactors poised for innovation. We propose that the high IDR fraction present in simplexvirus proteomes contributes to their wider host range compared with other herpesviruses.

Genetic defects in fungal recognition and susceptibility to invasive pulmonary aspergillosis

The interindividual variability in the onset and clinical course of invasive pulmonary aspergillosis (IPA) raises fundamental questions about its actual pathogenesis. Clinical and epidemiological studies have reported only a few examples of monogenic defects, however an expanding number of common polymorphisms associated with IPA has been identified. Understanding how genetic variation regulates the immune response to Aspergillus provides critical insights into the human immunobiology of IPA by pinpointing directly relevant immune molecules and interacting pathways. Most of the genetic defects reported to increase susceptibility to infection were described or suggested to impair fungal recognition by the innate immune system. In this review, we discuss the contribution of host genetic variation in pattern recognition receptors to the development of IPA. An improved understanding of the molecular and cellular processes that regulate human susceptibility to IPA is ultimately expected to pave the way toward personalized medical interventions based on host-directed risk stratification and individualized immunotherapy.

Selective constraint acting on TLR2 and TLR4 genes of Japanese Rana frogs

Toll-like receptors (TLRs) are an important component of innate immunity, the first line of pathogen defence. One of the major roles of TLRs includes recognition of pathogen-associated molecular patterns. Amphibians are currently facing population declines and even extinction due to chytridiomycosis caused by the Batrachochytrium dendrobatidis (Bd) fungus. Evidence from other vertebrates shows that TLR2 and TLR4 are involved in innate immunity against various fungi. Such genes therefore may play a functional role in amphibian-chytridiomycosis dynamics. Frogs from East Asia appear to be tolerant to Bd, so we examined the genetic diversity that underlies TLR2 and TLR4 from three Japanese Ranidae frog species, Rana japonica, R. ornativentris and R. tagoi tagoi (n = 5 per species). We isolated 27 TLR2 and 20 TLR4 alleles and found that these genes are evolutionarily conserved, with overall evidence supporting purifying selection. In contrast, site-by-site analysis of selection identified several specific codon sites under positive selection, some of which were located in the variable leucine rich repeat domains. In addition, preliminary expression levels of TLR2 and TLR4 from transcriptome data showed overall low expression. Although it remains unclear whether infectious pathogens are a selective force acting on TLRs of Japanese frogs, our results support that certain sites in TLRs of these species may have experienced pathogen-mediated selection.

Human genetics of infectious diseases: Unique insights into immunological redundancy

For almost any given human-tropic virus, bacterium, fungus, or parasite, the clinical outcome of primary infection is enormously variable, ranging from asymptomatic to lethal infection. This variability has long been thought to be largely determined by the germline genetics of the human host, and this is increasingly being demonstrated to be the case. The number and diversity of known inborn errors of immunity is continually increasing, and we focus here on autosomal and X-linked recessive traits underlying complete deficiencies of the encoded protein. Schematically, four types of infectious phenotype have been observed in individuals with such deficiencies, each providing information about the redundancy of the corresponding human gene, in terms of host defense in natural conditions. The lack of a protein can confer vulnerability to a broad range of microbes in most, if not all patients, through the disruption of a key immunological component. In such cases, the gene concerned is of low redundancy. However, the lack of a protein may also confer vulnerability to a narrow range of microbes, sometimes a single pathogen, and not necessarily in all patients. In such cases, the gene concerned is highly redundant. Conversely, the deficiency may be apparently neutral, conferring no detectable predisposition to infection in any individual. In such cases, the gene concerned is completely redundant. Finally, the lack of a protein may, paradoxically, be advantageous to the host, conferring resistance to one or more infections. In such cases, the gene is considered to display beneficial redundancy. These findings reflect the current state of evolution of humans and microbes, and should not be considered predictive of redundancy, or of a lack of redundancy, in the distant future. Nevertheless, these observations are of potential interest to present-day biologists testing immunological hypotheses experimentally and physicians managing patients with immunological or infectious conditions.

Host Genetic Signatures of Susceptibility to Fungal Disease

Our relative inability to predict the development of fungal disease and its clinical outcome raises fundamental questions about its actual pathogenesis. Several clinical risk factors are described to predispose to fungal disease, particularly in immunocompromised and severely ill patients. However, these alone do not entirely explain why, under comparable clinical conditions, only some patients develop infection. Recent clinical and epidemiological studies have reported an expanding number of monogenic defects and common polymorphisms associated with fungal disease. By directly implicating genetic variation in the functional regulation of immune mediators and interacting pathways, these studies have provided critical insights into the human immunobiology of fungal disease. Most of the common genetic defects reported were described or suggested to impair fungal recognition by the innate immune system. Here, we review common genetic variation in pattern recognition receptors and its impact on the immune response against the two major fungal pathogens Candida albicans and Aspergillus fumigatus. In addition, we discuss potential strategies and opportunities for the clinical translation of genetic information in the field of medical mycology. These approaches are expected to transfigure current clinical practice by unleashing an unprecedented ability to personalize prophylaxis, therapy and monitoring for fungal disease.

Rapidly Evolving Toll-3/4 Genes Encode Male-Specific Toll-Like Receptors in Drosophila

Animal Toll-like receptors (TLRs) have evolved through a pattern of duplication and divergence. Whereas mammalian TLRs directly recognize microbial ligands, Drosophila Tolls bind endogenous ligands downstream of both developmental and immune signaling cascades. Here, we find that most Toll genes in Drosophila evolve slowly with little gene turnover (gains/losses), consistent with their important roles in development and indirect roles in microbial recognition. In contrast, we find that the Toll-3/4 genes have experienced an unusually rapid rate of gene gains and losses, resulting in lineage-specific Toll-3/4s and vastly different gene repertoires among Drosophila species, from zero copies (e.g., D. mojavensis) to nineteen copies (e.g., D. willistoni). In D. willistoni, we find strong evidence for positive selection in Toll-3/4 genes, localized specifically to an extracellular region predicted to overlap with the binding site of Spätzle, the only known ligand of insect Tolls. However, because Spätzle genes are not experiencing similar selective pressures, we hypothesize that Toll-3/4s may be rapidly evolving because they bind to a different ligand, akin to TLRs outside of insects. We further find that most Drosophila Toll-3/4 genes are either weakly expressed or expressed exclusively in males, specifically in the germline. Unlike other Toll genes in D. melanogaster, Toll-3, and Toll-4 have apparently escaped from essential developmental roles, as knockdowns have no substantial effects on viability or male fertility. Based on these findings, we propose that the Toll-3/4 genes represent an exceptionally rapidly evolving lineage of Drosophila Toll genes, which play an unusual, as-yet-undiscovered role in the male germline.

Evolutionary analysis of Old World arenaviruses reveals a major adaptive contribution of the viral polymerase

The Old World (OW) arenavirus complex includes several species of rodent-borne viruses, some of which (i.e., Lassa virus, LASV and Lymphocytic choriomeningitis virus, LCMV) cause human diseases. Most LCMV and LASV infections are caused by rodent-to-human transmissions. Thus, viral evolution is largely determined by events that occur in the wildlife reservoirs. We used a set of human- and rodent-derived viral sequences to investigate the evolutionary history underlying OW arenavirus speciation, as well as the more recent selective events that accompanied LASV spread in West Africa. We show that the viral RNA polymerase (L protein) was a major positive selection target in OW arenaviruses and during LASV out-of-Nigeria migration. No evidence of selection was observed for the glycoprotein, whereas positive selection acted on the nucleoprotein (NP) during LCMV speciation. Positively selected sites in L and NP are surrounded by highly conserved residues, and the bulk of the viral genome evolves under purifying selection. Several positively selected sites are likely to modulate viral replication/transcription. In both L and NP, structural features (solvent exposed surface area) are important determinants of site-wise evolutionary rate variation. By incorporating several rodent-derived sequences, we also performed an analysis of OW arenavirus codon adaptation to the human host. Results do not support a previously hypothesized role of codon adaptation in disease severity for non-Nigerian strains. In conclusion, L and NP represent the major selection targets and possible determinants of disease presentation; these results suggest that field surveys and experimental studies should primarily focus on these proteins.

Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus)

Background

The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense.

Results

The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d_N/d_S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition.

Conclusions

This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3983-0) contains supplementary material, which is available to authorized users.

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

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

Evolution of toll-like receptors in the context of terrestrial ungulates and cetaceans diversification

Background

Toll-like receptors (TLRs) are the frontline actors in the innate immune response to various pathogens and are expected to be targets of natural selection in species adapted to habitats with contrasting pathogen burdens. The recent publication of genome sequences of giraffe and okapi together afforded the opportunity to examine the evolution of selected TLRs in broad range of terrestrial ungulates and cetaceans during their complex habitat diversification. Through direct sequence comparisons and standard evolutionary approaches, the extent of nucleotide and protein sequence diversity in seven Toll-like receptors (TLR2, TLR3, TLR4, TLR5, TLR7, TLR9 and TLR10) between giraffe and closely related species was determined. In addition, comparison of the patterning of key TLR motifs and domains between giraffe and related species was performed. The quantification of selection pressure and divergence on TLRs among terrestrial ungulates and cetaceans was also performed.

Results

Sequence analysis shows that giraffe has 94–99% nucleotide identity with okapi and cattle for all TLRs analyzed. Variations in the number of Leucine-rich repeats were observed in some of TLRs between giraffe, okapi and cattle. Patterning of key TLR domains did not reveal any significant differences in the domain architecture among giraffe, okapi and cattle. Molecular evolutionary analysis for selection pressure identifies positive selection on key sites for all TLRs examined suggesting that pervasive evolutionary pressure has taken place during the evolution of terrestrial ungulates and cetaceans. Analysis of positively selected sites showed some site to be part of Leucine-rich motifs suggesting functional relevance in species-specific recognition of pathogen associated molecular patterns. Notably, clade analysis reveals significant selection divergence between terrestrial ungulates and cetaceans in viral sensing TLR3. Mapping of giraffe TLR3 key substitutions to the structure of the receptor indicates that at least one of giraffe altered sites coincides with TLR3 residue known to play a critical role in receptor signaling activity.

Conclusion

There is overall structural conservation in TLRs among giraffe, okapi and cattle indicating that the mechanism for innate immune response utilizing TLR pathways may not have changed very much during the evolution of these species. However, a broader phylogenetic analysis revealed signatures of adaptive evolution among terrestrial ungulates and cetaceans, including the observed selection divergence in TLR3. This suggests that long term ecological dynamics has led to species-specific innovation and functional variation in the mechanisms mediating innate immunity in terrestrial ungulates and cetaceans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0901-7) contains supplementary material, which is available to authorized users.

Population genetics of mouse lemur vomeronasal receptors: current versus past selection and demographic inference

BACKGROUND

A major effort is underway to use population genetic approaches to identify loci involved in adaptation. One issue that has so far received limited attention is whether loci that show a phylogenetic signal of positive selection in the past also show evidence of ongoing positive selection at the population level. We address this issue using vomeronasal receptors (VRs), a diverse gene family in mammals involved in intraspecific communication and predator detection. In mouse lemurs, we previously demonstrated that both subfamilies of VRs (V1Rs and V2Rs) show a strong signal of directional selection in interspecific analyses. We predicted that ongoing sexual selection and/or co-evolution with predators may lead to current directional or balancing selection on VRs. Here, we re-sequence 17 VRs and perform a suite of selection and demographic analyses in sympatric populations of two species of mouse lemurs (Microcebus murinus and M. ravelobensis) in northwestern Madagascar.

RESULTS

M. ravelobensis had consistently higher genetic diversity at VRs than M. murinus. In general, we find little evidence for positive selection, with most loci evolving under purifying selection and one locus even showing evidence of functional loss in M. ravelobensis. However, a few loci in M. ravelobensis show potential evidence of positive selection. Using mismatch distributions and expansion models, we infer a more recent colonisation of the habitat by M. murinus than by M. ravelobensis, which most likely speciated in this region earlier on.

CONCLUSIONS

These findings suggest that the analysis of VR variation is useful in inferring demographic and phylogeographic history of mouse lemurs. In conclusion, this study reveals a substantial heterogeneity over time in selection on VR loci, suggesting that VR evolution is episodic.

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.

The evolutionary history of genes involved in spoken and written language: beyond FOXP2

Humans possess a communication system based on spoken and written language. Other animals can learn vocalization by imitation, but this is not equivalent to human language. Many genes were described to be implicated in language impairment (LI) and developmental dyslexia (DD), but their evolutionary history has not been thoroughly analyzed. Herein we analyzed the evolution of ten genes involved in DD and LI. Results show that the evolutionary history of LI genes for mammals and aves was comparable in vocal-learner species and non-learners. For the human lineage, several sites showing evidence of positive selection were identified in KIAA0319 and were already present in Neanderthals and Denisovans, suggesting that any phenotypic change they entailed was shared with archaic hominins. Conversely, in FOXP2, ROBO1, ROBO2, and CNTNAP2 non-coding changes rose to high frequency after the separation from archaic hominins. These variants are promising candidates for association studies in LI and DD.

Models for preclinical studies in aging-related disorders: One is not for all

Preclinical studies are essentially based on animal models of a particular disease. The primary purpose of preclinical efficacy studies is to support generalization of treatment-effect relationships to human subjects. Researchers aim to demonstrate a causal relationship between an investigational agent and a disease-related phenotype in such models. Numerous factors can muddle reliable inferences about such cause-effect relationships, including biased outcome assessment due to experimenter expectations. For instance, responses in a particular inbred mouse might be specific to the strain, limiting generalizability. Selecting well-justified and widely acknowledged model systems represents the best start in designing preclinical studies, especially to overcome any potential bias related to the model itself. This is particularly true in the research that focuses on aging, which carries unique challenges, mainly attributable to the fact that our already long lifespan makes designing experiments that use people as subjects extremely difficult and largely impractical.

Diverse selective regimes shape genetic diversity at ADAR genes and at their coding targets

A-to-I RNA editing operated by ADAR enzymes is extremely common in mammals. Several editing events in coding regions have pivotal physiological roles and affect protein sequence (recoding events) or function. We analyzed the evolutionary history of the 3 ADAR family genes and of their coding targets. Evolutionary analysis indicated that ADAR evolved adaptively in primates, with the strongest selection in the unique N-terminal domain of the interferon-inducible isoform. Positively selected residues in the human lineage were also detected in the ADAR deaminase domain and in the RNA binding domains of ADARB1 and ADARB2. During the recent history of human populations distinct variants in the 3 genes increased in frequency as a result of local selective pressures. Most selected variants are located within regulatory regions and some are in linkage disequilibrium with eQTLs in monocytes. Finally, analysis of conservation scores of coding editing sites indicated that editing events are counter-selected within regions that are poorly tolerant to change. Nevertheless, a minority of recoding events occurs at highly conserved positions and possibly represents the functional fraction. These events are enriched in pathways related to HIV-1 infection and to epidermis/hair development. Thus, both ADAR genes and their targets evolved under variable selective regimes, including purifying and positive selection. Pressures related to immune response likely represented major drivers of evolution for ADAR genes. As for their coding targets, we suggest that most editing events are slightly deleterious, although a minority may be beneficial and contribute to antiviral response and skin homeostasis.

Study of genetic diversity of KIR and TLR in the Rabhas, an endogamous primitive tribe of India

The Rabha tribe is a little known small endogamous population belonging to Indo-mongoloid group of north-eastern India. We have analyzed 16 KIR and 5 TLR gene polymorphisms in the Rabha population of northern West Bengal, India for the first time. The observed frequencies of the KIR genes (except framework and pseudogene loci) ranged between 0.26 (KIR2DS3) and 0.96 (KIR2DL1). Comparisons based on KIR polymorphism have revealed that although the Rabhas are of Indian origin the presence of mongoloid component in their gene pool cannot be denied. The frequencies of the 5 TLR genes ranged between 0.90 (TLR4) and 0.46 (TLR5). TLR variations found in the Rabhas may play a synergistic role in fighting against the bacterial invasions. Our results may contribute to the understanding of (1) genetic background and extent of genetic admixture in the Rabhas, (2) population migration events and (3) KIR-disease-TLR interactions.

Positive selection underlies the species-specific binding of Plasmodium falciparum RH5 to human basigin

Plasmodium falciparum, the causative agent of the deadliest form of malaria, is a member of the Laverania subgenus, which includes ape-infecting parasites. P. falciparum is thought to have originated in gorillas, although infection is now restricted to humans. Laverania parasites display remarkable host-specificity, which is partially mediated by the interaction between parasite ligands and host receptors. We analyse the evolution of BSG (basigin) and GYPA (glycophorin A) in primates/hominins, as well as of their Plasmodium-encoded ligands, PfRH5 and PfEBA175. We show that, in primates, positive selection targeted two sites in BSG (F27 and H102), both involved in PfRH5 binding. A population genetics-phylogenetics approach detected the strongest selection for the gorilla lineage: one of the positively selected sites (K191) is a major determinant of PfRH5 binding affinity. Analysis of RH5 genes indicated episodic selection on the P. falciparum branch; the positively selected W447 site is known to stabilize the interaction with human basigin. Conversely, we detect no selection in the receptor-binding region of EBA175 in the P. falciparum lineage. Its host receptor, GYPA, shows evidence of positive selection in all hominid lineages; selected codons include glycosylation sites that modulate PfEBA175 binding affinity. Data herein provide an evolutionary explanation for species-specific binding of the PfRH5-BSG ligand-receptor pair and support the hypothesis that positive selection at these genes drove the host shift leading to the emergence of P. falciparum as a human pathogen.

Causes and consequences of coagulation activation in sepsis: an evolutionary medicine perspective

BACKGROUND

Coagulation and innate immunity have been linked together for at least 450 million years of evolution. Sepsis, one of the world's leading causes of death, is probably the condition in which this evolutionary link is more evident. However, the biological and the clinical relevance of this association have only recently gained the attention of the scientific community.

DISCUSSION

During sepsis, the host response to a pathogen is invariably associated with coagulation activation. For several years, coagulation activation has been solely regarded as a mechanism of tissue damage, a concept that led to several clinical trials of anticoagulant agents for sepsis. More recently, this paradigm has been challenged by the failure of these clinical trials, and by a growing bulk of evidence supporting the concept that coagulation activation is beneficial for pathogen clearance. In this article we discuss recent basic and clinical data that point to a more balanced view of the detrimental and beneficial consequences of coagulation activation in sepsis. Reappraisal of the association between coagulation and immune activation from an evolutionary medicine perspective offers a unique opportunity to gain new insights about the pathogenesis of sepsis, paving the way to more successful approaches in both basic and clinical research in this field.

OASes and STING: adaptive evolution in concert

OAS (2′–5′-oligoadenylate synthases) proteins and cyclic GMP–AMP synthase (cGAS, gene symbol: MB21D1) patrol the cytoplasm for the presence of foreign nucleic acids. Upon binding to double-stranded RNA or double-stranded DNA, OAS proteins and cGAS produce nucleotide second messengers to activate RNase L and STING (stimulator of interferon genes, gene symbol: TMEM173), respectively; this leads to the initiation of antiviral responses. We analyzed the evolutionary history of the MB21D1–TMEM173 and OAS–RNASEL axes in primates and bats and found evidence of widespread positive selection in both orders. In TMEM173, residue 230, a major determinant of response to natural ligands and to mimetic drugs (e.g., DMXAA), was positively selected in Primates and Chiroptera. In both orders, selection also targeted an α-helix/loop element in RNase L that modulates the enzyme preference for single-stranded RNA versus stem loops. Analysis of positively selected sites in OAS1, OAS2, and MB21D1 revealed parallel evolution, with the corresponding residues being selected in different genes. As this cannot result from gene conversion, these data suggest that selective pressure acting on OAS and MB21D1 genes is related to nucleic acid recognition and to the specific mechanism of enzyme activation, which requires a conformational change. Finally, a population genetics-phylogenetics analysis in humans, chimpanzees, and gorillas detected several positively selected sites in most genes. Data herein shed light into species-specific differences in infection susceptibility and in response to synthetic compounds, with relevance for the design of synthetic compounds as vaccine adjuvants.

Interleukin-37 gene variants segregated anciently coexist during hominid evolution

IL37 is a member of IL-1 cytokine family but conveys anti-inflammatory functions. The biological characteristic and genetic heterogeneity of IL37 are not fully understood yet. Here using the whole-genome sequencing data from 1000 Genomes Project, we performed population and evolutionary genetic analysis of human IL37 gene. First, 2184 IL37 gene sequences from different human populations were retrieved. The IL37 protein sequences were inferred from the coding DNA sequences and multiple species alignment was made. Then, the phylogenetic tree of IL37 was built and dN/dS ratios were calculated for each evolutionary branch, the classic McDonald and Kreitman test was also performed. Next, we conducted intraspecific evolutionary genetic analysis and built the genealogy network of 116 unique IL37 haplotypes through median-joining network analysis. Finally, we compared IL37 sequences between the modern and archaic humans. Our results for the first time provide solid evidence that common IL37 variants other than NCBI reference sequence are present worldwide. Our data also supports that IL37 variants are shaped and maintained by selection instead of neutral evolution. We further identified that human IL37 variants consist of two major haplogroups and their presence in archaic humans corroborates its ancient origin in hominid evolution. In conclusion, these data indicate that common IL37 variants are maintained among human populations by selective force, suggesting their potential involvements in immune regulation and human diseases. In addition, the ancient history of IL37 variants reveals interesting insight into the complicated human evolutionary history.

Constraint and adaptation in newt toll-like receptor genes

Acute die-offs of amphibian populations worldwide have been linked to the emergence of viral and fungal diseases. Inter and intraspecific immunogenetic differences may influence the outcome of infection. Toll-like receptors (TLRs) are an essential component of innate immunity and also prime acquired defenses. We report the first comprehensive assessment of TLR gene variation for urodele amphibians. The Lissotriton newt TLR repertoire includes representatives of 13 families and is compositionally most similar to that of the anuran Xenopus. Both ancient and recent gene duplications have occurred in urodeles, bringing the total number of TLR genes to at least 21. Purifying selection has predominated the evolution of newt TLRs in both long (∼70 Ma) and medium (∼18 Ma) timescales. However, we find evidence for both purifying and positive selection acting on TLRs in two recently diverged (2-5 Ma) allopatric evolutionary lineages (Lissotriton montandoni and L. vulgaris graecus). Overall, both forms of selection have been stronger in L. v. graecus, while constraint on most TLR genes in L. montandoni appears relaxed. The differences in selection regimes are unlikely to be biased by demographic effects because these were controlled by means of a historical demographic model derived from an independent data set of 62 loci. We infer that TLR genes undergo distinct trajectories of adaptive evolution in closely related amphibian lineages, highlight the potential of TLRs to capture the signatures of different assemblages of pathogenic microorganisms, and suggest differences between lineages in the relative roles of innate and acquired immunity.

Albuminoid genes: evolving at the interface of dispensability and selection

The albuminoid gene family comprises vitamin D-binding protein (GC), alpha-fetoprotein (AFP), afamin (AFM), and albumin (ALB). Albumin is the most abundant human serum protein, and, as the other family members, acts as a transporter of endogenous and exogenous substances including thyroxine, fatty acids, and drugs. Instead, the major cargo of GC is 25-hydroxyvitamin D. We performed an evolutionary study of albuminoid genes and we show that ALB evolved adaptively in mammals. Most positively selected sites are located within albumin-binding sites for fatty acids and thyroxine, as well as at the contact surface with neonatal Fc receptor. Positive selection was also detected for residues forming the prostaglandin-binding pocket. Adaptation to hibernation/torpor might explain the signatures of episodic positive selection we detected for few mammalian lineages. Application of a population genetics-phylogenetics approach showed that purifying selection represented a major force acting on albuminoid genes in both humans and chimpanzees, with the strongest constraint observed for human GC. Population genetic analysis revealed that GC was also the target of locally exerted selective pressure, which drove the frequency increase of different haplotypes in distinct human populations. A search for known variants that modulate GC and 25-hydroxyvitamin D concentrations revealed linkage disequilibrium with positively selected variants, although European and Asian major GC haplotypes carry alleles with reported opposite effect on GC concentration. Data herein indicate that albumin, an extremely abundant housekeeping protein, was the target of pervasive and episodic selection in mammals, whereas GC represented a selection target during the recent evolution of human populations.

Human genome variability, natural selection and infectious diseases

The recent availability of large-scale sequencing DNA data allowed researchers to investigate how genomic variation is distributed among populations. While demographic factors explain genome-wide population genetic diversity levels, scans for signatures of natural selection pinpointed several regions under non-neutral evolution. Recent studies found an enrichment of immune-related genes subjected to natural selection, suggesting that pathogens and infectious diseases have imposed a strong selective pressure throughout human history. Pathogen-mediated selection often targeted regulatory sites of genes belonging to the same biological pathway. Results from these studies have the potential to identify mutations that modulate infection susceptibility by integrating a population genomic approach with molecular immunology data and large-scale functional annotations.