Rachele Cagliani

Parasites represent a major selective force for interleukin genes and shape the genetic predisposition to autoimmune conditions

Many human genes have adapted to the constant threat of exposure to infectious agents; according to the “hygiene hypothesis,” lack of exposure to parasites in modern settings results in immune imbalances, augmenting susceptibility to the development of autoimmune and allergic conditions. Here, by estimating the number of pathogen species/genera in a specific geographic location (pathogen richness) for 52 human populations and analyzing 91 interleukin (IL)/IL receptor genes (IL genes), we show that helminths have been a major selective force on a subset of these genes. A population genetics analysis revealed that five IL genes, including IL7R and IL18RAP, have been a target of balancing selection, a selection process that maintains genetic variability within a population. Previous identification of polymorphisms in some of these loci, and their association with autoimmune conditions, prompted us to investigate the relationship between adaptation and disease. By searching for variants in IL genes identified in genome-wide association studies, we verified that six risk alleles for inflammatory bowel (IBD) or celiac disease are significantly correlated with micropathogen richness. These data support the hygiene hypothesis for IBD and provide a large set of putative targets for susceptibility to helminth infections.

Widespread balancing selection and pathogen-driven selection at blood group antigen genes

Historically, allelic variations in blood group antigen (BGA) genes have been regarded as possible susceptibility factors for infectious diseases. Since host-pathogen interactions are major determinants in evolution, BGAs can be thought of as selection targets. In order to verify this hypothesis, we obtained an estimate of pathogen richness for geographic locations corresponding to 52 populations distributed worldwide; after correction for multiple tests and for variables different from selective forces, significant correlations with pathogen richness were obtained for multiple variants at 11 BGA loci out of 26. In line with this finding, we demonstrate that three BGA genes, namely CD55, CD151, and SLC14A1, have been subjected to balancing selection, a process, rare outside MHC genes, which maintains variability at a locus. Moreover, we identified a gene region immediately upstream the transcription start site of FUT2 which has undergone non-neutral evolution independently from the coding region. Finally, in the case of BSG, we describe the presence of a highly divergent haplotype clade and the possible reasons for its maintenance, including frequency-dependent balancing selection, are discussed. These data indicate that BGAs have been playing a central role in the host-pathogen arms race during human evolutionary history and no other gene category shows similar levels of widespread selection, with the only exception of loci involved in antigen recognition.

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.

Genetic diversity at endoplasmic reticulum aminopeptidases is maintained by balancing selection and is associated with natural resistance to HIV-1 infection

Human ERAP1 and ERAP2 encode two endoplasmic reticulum aminopeptidases. These enzymes trim peptides to optimal size for loading onto major histocompatibility complex class I molecules and shape the antigenic repertoire presented to CD8(+) T cells. Therefore, ERAP1 and ERAP2 may be considered potential selection targets and modulators of infection susceptibility. We resequenced two genic regions in ERAP1 and ERAP2 in three HapMap populations. In both cases, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, and reduced population genetic differentiation. The genealogy of ERAP1 and ERAP2 haplotypes was split into two major branches with deep coalescence times. These features suggest that long-standing balancing selection has acted on these genes. Analysis of the Lys528Arg (rs30187 in ERAP1) and Asn392Lys (rs2549782 in ERAP2) variants in an Italian population of HIV-1-exposed seronegative (ESN) individuals and a larger number of Italian controls indicated that rs2549782 significantly deviates from Hardy-Weinberg equilibrium (HWE) in ESN but not in controls. Technical errors were excluded and a goodness-of-fit test indicated that a recessive model with only genetic effects adequately explains HWE deviation. The genotype distribution of rs2549782 is significantly different in the two cohorts (P = 0.004), mainly as the result of an over-representation of Lys/Lys genotypes in the ESN sample (P-value for a recessive model: 0.00097). Our data suggest that genetic diversity in ERAP1 and ERAP2 has been maintained by balancing selection and that variants in ERAP2 confer resistance to HIV-1 infection possibly via the presentation of a distinctive peptide repertoire to CD8(+) T cells.

The landscape of human genes involved in the immune response to parasitic worms

BackgroundMore than 2 billion individuals worldwide suffer from helminth infections. The highest parasite burdens occur in children and helminth infection during pregnancy is a risk factor for preterm delivery and reduced birth weight. Therefore, helminth infections can be regarded as a strong selective pressure.ResultsHere we propose that candidate susceptibility genes for parasitic worm infections can be identified by searching for SNPs that display a strong correlation with the diversity of helminth species/genera transmitted in different geographic areas. By a genome-wide search we identified 3478 variants that correlate with helminth diversity. These SNPs map to 810 distinct human genes including loci involved in regulatory T cell function and in macrophage activation, as well as leukocyte integrins and co-inhibitory molecules. Analysis of functional relationships among these genes identified complex interaction networks centred around Th2 cytokines. Finally, several genes carrying candidate targets for helminth-driven selective pressure also harbour susceptibility alleles for asthma/allergy or are involved in airway hyper-responsiveness, therefore expanding the known parallelism between these conditions and parasitic infections.ConclusionsOur data provide a landscape of human genes that modulate susceptibility to helminths and indicate parasitic worms as one of the major selective forces in humans.

The signature of long-standing balancing selection at the human defensin β-1 promoter

BackgroundDefensins, small endogenous peptides with antimicrobial activity, are pivotal components of the innate immune response. A large cluster of defensin genes is located on human chromosome 8p; among them the beta defensin 1 (DEFB1) promoterhas been extensively studied since discovery that specific polymorphisms and haplotypes associate with asthma and atopy, susceptibility to severe sepsis, as well as HIV and Candida infection predisposition.ResultsHere, we characterize the sequence variation and haplotype structure of the DEFB1 promoter region in six human populations. In all of them, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, reduced population differentiation and a genealogy with common haplotypes separated by deep branches. Indeed, a significant departure from the expectation of evolutionary neutrality was observed in all populations and the possibility that this is due to demographic history alone was ruled out. Also, we verified that the selection signature is restricted to the promoter region and not due to a linked balanced polymorphism. A phylogeny-based estimation indicated that the two major haplotype clades separated around 4.5 million years ago, approximately the time when the human and chimpanzee lineages split.ConclusionAltogether, these features represent strong molecular signatures of long-term balancing selection, a process that is thought to be extremely rare outside major histocompatibility complex genes. Our data indicate that the DEFB1 promoter region carries functional variants and support previous hypotheses whereby alleles predisposing to atopic disorders are widespread in modern societies because they conferred resistance to pathogens in ancient settings.

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.

Both selective and neutral processes drive GC content evolution in the human genome

BackgroundMammalian genomes consist of regions differing in GC content, referred to as isochores or GC-content domains. The scientific debate is still open as to whether such compositional heterogeneity is a selected or neutral trait.ResultsHere we analyze SNP allele frequencies, retrotransposon insertion polymorphisms (RIPs), as well as fixed substitutions accumulated in the human lineage since its divergence from chimpanzee to indicate that biased gene conversion (BGC) has been playing a role in within-genome GC content variation. Yet, a distinct contribution to GC content evolution is accounted for by a selective process. Accordingly, we searched for independent evidences that GC content distribution does not conform to neutral expectations. Indeed, after correcting for possible biases, we show that intron GC content and size display isochore-specific correlations.ConclusionWe consider that the more parsimonious explanation for our results is that GC content is subjected to the action of both weak selection and BGC in the human genome with features such as nucleosome positioning or chromatin conformation possibly representing the final target of selective processes. This view might reconcile previous contrasting findings and add some theoretical background to recent evidences suggesting that GC content domains display different behaviors with respect to highly regulated biological processes such as developmentally-stage related gene expression and programmed replication timing during neural stem cell differentiation.

Genome-Wide Identification of Susceptibility Alleles for Viral Infections through a Population Genetics Approach

Viruses have exerted a constant and potent selective pressure on human genes throughout evolution. We utilized the marks left by selection on allele frequency to identify viral infection-associated allelic variants. Virus diversity (the number of different viruses in a geographic region) was used to measure virus-driven selective pressure. Results showed an excess of variants correlated with virus diversity in genes involved in immune response and in the biosynthesis of glycan structures functioning as viral receptors; a significantly higher than expected number of variants was also seen in genes encoding proteins that directly interact with viral components. Genome-wide analyses identified 441 variants significantly associated with virus-diversity; these are more frequently located within gene regions than expected, and they map to 139 human genes. Analysis of functional relationships among genes subjected to virus-driven selective pressure identified a complex network enriched in viral products-interacting proteins. The novel approach to the study of infectious disease epidemiology presented herein may represent an alternative to classic genome-wide association studies and provides a large set of candidate susceptibility variants for viral infections.