Diego Forni

A Common Polymorphism in TLR3 Confers Natural Resistance to HIV-1 Infection

TLR3 recognizes dsRNA and activates antiviral immune responses through the production of inflammatory cytokines and type I IFNs. Genetic association studies have provided evidence concerning the role of a polymorphism in TLR3 (rs3775291, Leu412Phe) in viral infection susceptibility. We genotyped rs3775291 in a population of Spanish HIV-1–exposed seronegative (HESN) individuals who remain HIV seronegative despite repeated exposure through i.v. injection drug use (IDU-HESN individuals) as witnessed by their hepatitis C virus seropositivity. The frequency of individuals carrying at least one 412Phe allele was significantly higher in IDU-HESN individuals compared with that of a matched control sample (odds ratio for a dominant model = 1.87; 95% confidence interval, 1.06–3.34; p = 0.023). To replicate this finding, we analyzed a cohort of Italian, sexually HESN individuals. Similar results were obtained: the frequency of individuals carrying at least one 412Phe allele was significantly higher compared with that of a matched control sample (odds ratio, 1.79; 95% confidence interval, 1.05–3.08; p = 0.029). In vitro infection assays showed that in PBMCs carrying the 412Phe allele, HIV-1Ba-L replication was significantly reduced (p = 0.025) compared with that of Leu/Leu homozygous samples and was associated with a higher expression of factors suggestive of a state of immune activation (IL-6, CCL3, CD69). Similarly, stimulation of PBMCs with a TLR3 agonist indicated that the presence of the 412Phe allele results in a significantly increased expression of CD69 and higher production of proinflammatory cytokines including IL-6 and CCL3. The data of this study indicate that a common TLR3 allele confers immunologically mediated protection from HIV-1 and suggest the potential use of TLR3 triggering in HIV-1 immunotherapy.

Evolutionary insights into host–pathogen interactions from mammalian sequence data

Infections are one of the major selective pressures acting on humans, and host-pathogen interactions contribute to shaping the genetic diversity of both organisms. Evolutionary genomic studies take advantage of experiments that natural selection has been performing over millennia. In particular, inter-species comparative genomic analyses can highlight the genetic determinants of infection susceptibility or severity. Recent examples show how evolution-guided approaches can provide new insights into host–pathogen interactions, ultimately clarifying the basis of host range and explaining the emergence of different diseases. We describe the latest developments in comparative immunology and evolutionary genetics, showing their relevance for understanding the molecular determinants of infection susceptibility in mammals.

A functional variant in ERAP1 predisposes to multiple sclerosis.

The ERAP1 gene encodes an aminopeptidase involved in antigen processing. A functional polymorphism in the gene (rs30187, Arg528Lys) associates with susceptibility to ankylosying spondylitis (AS), whereas a SNP in the interacting ERAP2 gene increases susceptibility to another inflammatory autoimmune disorder, Crohns disease (CD). We analysed rs30187 in 572 Italian patients with CD and in 517 subjects suffering from multiple sclerosis (MS); for each cohort, an independent sex- and age-matched control group was genotyped. The frequency of the 528Arg allele was significantly higher in both disease cohorts compared to the respective control population (for CD, OR = 1.20 95%CI: 1.01–1.43, p = 0.036; for RRMS, OR = 1.26; 95%CI: 1.04–1.51, p = 0.01). Meta-analysis with the Wellcome Trust Cases Control Consortium GWAS data confirmed the association with MS (pmeta = 0.005), but not with CD. In AS, the rs30187 variant has a predisposing effect only in an HLA-B27 allelic background. It remains to be evaluated whether interaction between ERAP1 and distinct HLA class I alleles also affects the predisposition to MS, and explains the failure to provide definitive evidence for a role of rs30187 in CD. Results herein support the emerging concept that a subset of master-regulatory genes underlay the pathogenesis of autoimmunity.

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.

Nonstructural Proteins Are Preferential Positive Selection Targets in Zika Virus and Related Flaviviruses

The Flavivirus genus comprises several human pathogens such as dengue virus (DENV), Japanese encephalitis virus (JEV), and Zika virus (ZIKV). Although ZIKV usually causes mild symptoms, growing evidence is linking it to congenital birth defects and to increased risk of Guillain-Barré syndrome. ZIKV encodes a polyprotein that is processed to produce three structural and seven nonstructural (NS) proteins. We investigated the evolution of the viral polyprotein in ZIKV and in related flaviviruses (DENV, Spondweni virus, and Kedougou virus). After accounting for saturation issues, alignment uncertainties, and recombination, we found evidence of episodic positive selection on the branch that separates DENV from the other flaviviruses. NS1 emerged as the major selection target, and selected sites were located in immune epitopes or in functionally important protein regions. Three of these sites are located in an NS1 region that interacts with structural proteins and is essential for virion biogenesis. Analysis of the more recent evolutionary history of ZIKV lineages indicated that positive selection acted on NS5 and NS4B, this latter representing the preferential target. All selected sites were located in the N-terminal portion of NS4B, which inhibits interferon response. One of the positively selected sites (26M/I/T/V) in ZIKV also represents a selection target in sylvatic DENV2 isolates, and a nearby residue evolves adaptively in JEV. Two additional positively selected sites are within a protein region that interacts with host (e.g. STING) and viral (i.e. NS1, NS4A) proteins. Notably, mutations in the NS4B region of other flaviviruses modulate neurovirulence and/or neuroinvasiveness. These results suggest that the positively selected sites we identified modulate viral replication and contribute to immune evasion. These sites should be prioritized in future experimental studies. However, analyses herein detected no selective events associated to the spread of the Asian/American ZIKV lineage.

A Nonsense Polymorphism (R392X) in TLR5 Protects from Obesity but Predisposes to Diabetes

The TLR5 gene encodes an innate immunity receptor. Mice lacking Tlr5 (T5KO) develop insulin resistance and increased adiposity. Owing to the segregation of a dominant nonsense polymorphism (R392X, rs5744168), a portion of humans lack TLR5 function. We investigated whether the nonsense polymorphism influences obesity and susceptibility to type 2 diabetes (T2D). R392X was genotyped in two cohorts from Saudi Arabia, a region where obesity and type 2 diabetes (T2D) are highly prevalent. The nonsense allele was found to protect from obesity (pcombined = 0.0062; odds ratio, 0.51) and to associate with lower body mass index (BMI) (pcombined = 0.0061); this allele also correlated with a reduced production of proinflammatory cytokines. A significant interaction was noted between rs5744168 and sex in affecting BMI (pinteraction = 0.006), and stratification by gender revealed that the association is driven by females (pcombined = 0.0016 and 0.0006 for obesity and BMI, respectively). The nonsense polymorphism also associated with BMI in nonobese women. After correction for BMI, the 392X allele was found to represent a risk factor for T2D with a sex-specific effect (pinteraction = 0.023) mediated by females (p = 0.021; odds ratio, 2.60). Fasting plasma glucose levels in nondiabetic individuals were also higher in women carrying the nonsense allele (p = 0.012). Thus, in contrast to T5KO mice, loss of human TLR5 function protects from weight gain, but in analogy to the animal model, the nonsense allele predisposes to T2D. These effects are apparently sex-specific. Data in this study reinforce the hypothesis that metabolic diseases, including T2D, are associated with immune dysregulation.

An Evolutionary Analysis of Antigen Processing and Presentation across Different Timescales Reveals Pervasive Selection

The antigenic repertoire presented by MHC molecules is generated by the antigen processing and presentation (APP) pathway. We analyzed the evolutionary history of 45 genes involved in APP at the inter- and intra-species level. Results showed that 11 genes evolved adaptively in mammals. Several positively selected sites involve positions of fundamental importance to the protein function (e.g. the TAP1 peptide-binding domains, the sugar binding interface of langerin, and the CD1D trafficking signal region). In CYBB, all selected sites cluster in two loops protruding into the endosomal lumen; analysis of missense mutations responsible for chronic granulomatous disease (CGD) showed the action of different selective forces on the very same gene region, as most CGD substitutions involve aminoacid positions that are conserved in all mammals. As for ERAP2, different computational methods indicated that positive selection has driven the recurrent appearance of protein-destabilizing variants during mammalian evolution. Application of a population-genetics phylogenetics approach showed that purifying selection represented a major force acting on some APP components (e.g. immunoproteasome subunits and chaperones) and allowed identification of positive selection events in the human lineage. We also investigated the evolutionary history of APP genes in human populations by developing a new approach that uses several different tests to identify the selection target, and that integrates low-coverage whole-genome sequencing data with Sanger sequencing. This analysis revealed that 9 APP genes underwent local adaptation in human populations. Most positive selection targets are located within noncoding regions with regulatory function in myeloid cells or act as expression quantitative trait loci. Conversely, balancing selection targeted nonsynonymous variants in TAP1 and CD207 (langerin). Finally, we suggest that selected variants in PSMB10 and CD207 contribute to human phenotypes. Thus, we used evolutionary information to generate experimentally-testable hypotheses and to provide a list of sites to prioritize in follow-up analyses.

RIG-I-like receptors evolved adaptively in mammals, with parallel evolution at LGP2 and RIG-I.

RIG-I-like receptors (RLRs) are nucleic acid sensors that activate antiviral innate immune response. These molecules recognize diverse non-self RNA substrates and are antagonized by several viral inhibitors. We performed an evolutionary analysis of RLR genes (RIG-I, MDA5, and LGP2) in mammals. Results indicated that purifying selection had a dominant role in driving the evolution of RLRs. However, application of maximum-likelihood analyses identified several positions that evolved adaptively. Positively selected sites are located in all domains of MDA5 and RIG-I, whereas in LGP2 they are confined to the helicase domain. In both MDA5 and RIG-I, the linkers separating the caspase activation and recruitment domain and the helicase domain represented preferential targets of positive selection. Independent selective events in RIG-I and LGP2 targeted the corresponding site (Asp421 and Asp179, respectively) within a protruding α-helix that grips the V-shaped structure formed by the pincer. Most of the positively selected sites in MDA5 are in regions unique to this RLR, including a characteristic insertion within the helicase domain. Additional selected sites are located at the contact interface between MDA5 monomers, in spatial proximity to a positively selected human polymorphism (Arg843His) and immediately external to the parainfluenza virus 5 V protein binding region. Structural analyses suggested that the positively selected His834 residue is involved in parainfluenza virus 5 V protein binding. Data herein suggest that RLRs have been engaged in host-virus genetic conflict leading to diversifying selection and indicate parallel evolution at the same site in RIG-I and LGP2, a position likely to be of central importance in antiviral responses.

Genetic adaptation of the human circadian clock to day-length latitudinal variations and relevance for affective disorders

BackgroundThe temporal coordination of biological processes into daily cycles is a common feature of most living organisms. In humans, disruption of circadian rhythms is commonly observed in psychiatric diseases, including schizophrenia, bipolar disorder, depression and autism. Light therapy is the most effective treatment for seasonal affective disorder and circadian-related treatments sustain antidepressant response in bipolar disorder patients. Day/night cycles represent a major circadian synchronizing signal and vary widely with latitude.ResultsWe apply a geographically explicit model to show that out-of-Africa migration, which led humans to occupy a wide latitudinal area, affected the evolutionary history of circadian regulatory genes. The SNPs we identify using this model display consistent signals of natural selection using tests based on population genetic differentiation and haplotype homozygosity. Signals of natural selection driven by annual photoperiod variation are detected for schizophrenia, bipolar disorder, and restless leg syndrome risk variants, in line with the circadian component of these conditions.ConclusionsOur results suggest that human populations adapted to life at different latitudes by tuning their circadian clock systems. This process also involves risk variants for neuropsychiatric conditions, suggesting possible genetic modulators for chronotherapies and candidates for interaction analysis with photoperiod-related environmental variables, such as season of birth, country of residence, shift-work or lifestyle habits.

Ancient and Recent Selective Pressures Shaped Genetic Diversity at AIM2-Like Nucleic Acid Sensors

AIM2-like receptors (ALRs) are a family of nucleic acid sensors essential for innate immune responses against viruses and bacteria. We performed an evolutionary analysis of ALR genes (MNDA, PYHIN1, IFI16, and AIM2) by analyzing inter- and intraspecies diversity. Maximum-likelihood analyses indicated that IFI16 and AIM2 evolved adaptively in primates, with branch-specific selection at the catarrhini lineage for IFI16. Application of a population genetics–phylogenetics approach also allowed identification of positive selection events in the human lineage. Positive selection in primates targeted sites located at the DNA-binding interface in both IFI16 and AIM2. In IFI16, several sites positively selected in primates and in the human lineage were located in the PYD domain, which is involved in protein–protein interaction and is bound by a human cytomegalovirus immune evasion protein. Finally, positive selection was found to target nuclear localization signals in IFI16 and the spacer region separating the two HIN domains. Population genetic analysis in humans revealed that an IFI16 genic region has been a target of long-standing balancing selection, possibly acting on two nonsynonymous polymorphisms located in the spacer region. Data herein indicate that ALRs have been repeatedly targeted by natural selection. The balancing selection region in IFI16 carries a variant with opposite risk effect for distinct autoimmune diseases, suggesting antagonistic pleiotropy. We propose that the underlying scenario is the result of an ancestral and still ongoing host–pathogen arms race and that the maintenance of susceptibility alleles for autoimmune diseases at IFI16 represents an evolutionary trade-off.