Balancing selection is common in the extended MHC region but most alleles with opposite risk profile for autoimmune diseases are neutrally evolving

Emerging Roles of MHC Class I Region-Encoded E3 Ubiquitin Ligases in Innate Immunity

The major histocompatibility complex (MHC) class I (MHC-I) region contains a multitude of genes relevant to immune response. Multiple E3 ubiquitin ligase genes, including tripartite motif 10 (TRIM10), TRIM15, TRIM26, TRIM27, TRIM31, TRIM38, TRIM39, TRIM40, and RING finger protein 39 (RNF39), are organized in a tight cluster, and an additional two TRIM genes (namely TRIM38 and TRIM27) telomeric of the cluster within the MHC-I region. The E3 ubiquitin ligases encoded by these genes possess important roles in controlling the intensity of innate immune responses. In this review, we discuss the E3 ubiquitin ligases encoded within the MHC-I region, highlight their regulatory roles in innate immunity, and outline their potential functions in infection, inflammatory and autoimmune diseases.

Long Noncoding RNAs in Host-Pathogen Interactions

Long noncoding RNAs (lncRNAs) are key molecules that regulate gene expression in a variety of organisms. LncRNAs can drive different transcriptional and post-transcriptional events that impact cellular functions. Recent studies have identified many lncRNAs associated with immune cell development and activation; however, an understanding of their functional role in host immunity to infection is just emerging. Here, we provide a detailed and updated review of the functional roles of lncRNAs in regulating mammalian immune responses during host-pathogen interactions, because these functions may be either beneficial or detrimental to the host. With increased mechanistic insight into the roles of lncRNAs, it may be possible to design and/or improve lncRNA-based therapies to treat a variety of infectious and inflammatory diseases.

Influence of ERAP1 and ERAP2 gene polymorphisms on disease susceptibility in different populations

The endoplasmic reticulum aminopeptidases (ERAPs), ERAP1 and ERAP2, makes a role in shaping the HLA class I peptidome by trimming peptides to the optimal size in MHC-class I-mediated antigen presentation and educating the immune system to differentiate between self-derived and foreign antigens. Association studies have shown that genetic variations in ERAP1 and ERAP2 genes increase susceptibility to autoimmune diseases, infectious diseases, and cancers. Both ERAP1 and ERAP2 genes exhibit diverse polymorphisms in different populations, which may influence their susceptibly to the aforementioned diseases. In this article, we review the distribution of ERAP1 and ERAP2 gene polymorphisms in various populations; discuss the risk or protective influence of these gene polymorphisms in autoimmune diseases, infectious diseases, and cancers; and highlight how ERAP genetic variations can influence disease associations.

Population Genetics and Natural Selection in Rheumatic Disease

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Natural selection is an important influence on human genetic variation. Because immune and inflammatory function genes are enriched for signals of positive selection, the prevalence of rheumatic disease-risk alleles seen in different populations is partially the result of differing selective pressures (eg, due to pathogens). This review summarizes the genetic regions associated with susceptibility to different rheumatic diseases and concomitant evidence for natural selection, including known agents of selection exerting selective pressure in these regions.

Natural Selection in the Great Apes

Natural selection is crucial for the adaptation of populations to their environments. Here, we present the first global study of natural selection in the Hominidae (humans and great apes) based on genome-wide information from population samples representing all extant species (including most subspecies). Combining several neutrality tests we create a multi-species map of signatures of natural selection covering all major types of natural selection. We find that the estimated efficiency of both purifying and positive selection varies between species and is significantly correlated with their long-term effective population size. Thus, even the modest differences in population size among the closely related Hominidae lineages have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments. Most signatures of balancing and positive selection are species-specific, with signatures of balancing selection more often being shared among species. We also identify loci with evidence of positive selection across several lineages. Notably, we detect signatures of positive selection in several genes related to brain function, anatomy, diet and immune processes. Our results contribute to a better understanding of human evolution by putting the evidence of natural selection in humans within its larger evolutionary context. The global map of natural selection in our closest living relatives is available as an interactive browser at http://tinyurl.com/nf8qmzh.

An Ancient Fecundability-Associated Polymorphism Switches a Repressor into an Enhancer of Endometrial TAP2 Expression

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, is heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits has been challenging. Here, we explore the functional significance and evolutionary history of a T/C polymorphism of SNP rs2071473, which we have previously shown is an eQTL for TAP2 and significantly associated with fecundability (time to pregnancy). We replicated the association between the rs2071473 genotype and TAP2 expression by using GTEx data and demonstrated that TAP2 is expressed by decidual stromal cells at the maternal-fetal interface. Next, we showed that rs2071473 is located within a progesterone-responsive cis-regulatory element that functions as a repressor with the T allele and an enhancer with the C allele. Remarkably, we found that this polymorphism arose before the divergence of modern and archaic humans, segregates at intermediate to high frequencies across human populations, and has genetic signatures of long-term balancing selection. This variant has also previously been identified in genome-wide association studies of immune-related disease, suggesting that both alleles are maintained as a result of antagonistic pleiotropy.

Ancient pathogen-driven adaptation triggers increased susceptibility to non-celiac wheat sensitivity in present-day European populations

BACKGROUND

Non-celiac wheat sensitivity is an emerging wheat-related syndrome showing peak prevalence in Western populations. Recent studies hypothesize that new gliadin alleles introduced in the human diet by replacement of ancient wheat with modern varieties can prompt immune responses mediated by the CXCR3-chemokine axis potentially underlying such pathogenic inflammation. This cultural shift may also explain disease epidemiology, having turned European-specific adaptive alleles previously targeted by natural selection into disadvantageous ones.

METHODS

To explore this evolutionary scenario, we performed ultra-deep sequencing of genes pivotal in the CXCR3-inflammatory pathway on individuals diagnosed for non-celiac wheat sensitivity and we applied anthropological evolutionary genetics methods to sequence data from worldwide populations to investigate the genetic legacy of natural selection on these loci.

RESULTS

Our results indicate that balancing selection has maintained two divergent CXCL10/CXCL11 haplotypes in Europeans, one responsible for boosting inflammatory reactions and another for encoding moderate chemokine expression.

CONCLUSIONS

This led to considerably higher occurrence of the former haplotype in Western people than in Africans and East Asians, suggesting that they might be more prone to side effects related to the consumption of modern wheat varieties. Accordingly, this study contributed to shed new light on some of the mechanisms potentially involved in the disease etiology and on the evolutionary bases of its present-day epidemiological patterns. Moreover, overrepresentation of disease homozygotes for the dis-adaptive haplotype plausibly accounts for their even more enhanced CXCR3-axis expression and for their further increase in disease risk, representing a promising finding to be validated by larger follow-up studies.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

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.

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.

Antigen-presenting cells digest intracellular and extracellular proteins and display the resulting antigenic repertoire on cell surface molecules for recognition by T cells. This process initiates cell-mediated immune responses and is essential to detect infections. The antigenic repertoire is generated by the antigen processing and presentation pathway. Because several pathogens evade immune recognition by hampering this process, genes involved in antigen processing and presentation may represent common natural selection targets. Thus, we analyzed the evolutionary history of these genes during mammalian evolution and in the more recent history of human populations. Evolutionary analyses in mammals indicated that positive selection targeted a very high proportion of genes (24%), and revealed that many selected sites affect positions of fundamental importance to the protein function. In humans, we found different signatures of natural selection acting both on regions that are expected to regulate gene expression levels or timing and on coding variants; two human selected polymorphisms may modulate the susceptibility to Crohn's disease and to HIV-1 infection. Therefore, we provide a comprehensive evolutionary analysis of antigen processing and we show that evolutionary studies can provide useful information concerning the location and nature of functional variants, ultimately helping to clarify phenotypic differences between and within species.

Candidate gene approach for parasite resistance in sheep--variation in immune pathway genes and association with fecal egg count

Sheep chromosome 3 (Oar3) has the largest number of QTLs reported to be significantly associated with resistance to gastro-intestinal nematodes. This study aimed to identify single nucleotide polymorphisms (SNPs) within candidate genes located in sheep chromosome 3 as well as genes involved in major immune pathways. A total of 41 SNPs were identified across 38 candidate genes in a panel of unrelated sheep and genotyped in 713 animals belonging to 22 breeds across Asia, Europe and South America. The variations and evolution of immune pathway genes were assessed in sheep populations across these macro-environmental regions that significantly differ in the diversity and load of pathogens. The mean minor allele frequency (MAF) did not vary between Asian and European sheep reflecting the absence of ascertainment bias. Phylogenetic analysis revealed two major clusters with most of South Asian, South East Asian and South West Asian breeds clustering together while European and South American sheep breeds clustered together distinctly. Analysis of molecular variance revealed strong phylogeographic structure at loci located in immune pathway genes, unlike microsatellite and genome wide SNP markers. To understand the influence of natural selection processes, SNP loci located in chromosome 3 were utilized to reconstruct haplotypes, the diversity of which showed significant deviations from selective neutrality. Reduced Median network of reconstructed haplotypes showed balancing selection in force at these loci. Preliminary association of SNP genotypes with phenotypes recorded 42 days post challenge revealed significant differences (P<0.05) in fecal egg count, body weight change and packed cell volume at two, four and six SNP loci respectively. In conclusion, the present study reports strong phylogeographic structure and balancing selection operating at SNP loci located within immune pathway genes. Further, SNP loci identified in the study were found to have potential for future large scale association studies in naturally exposed sheep populations.

A novel single nucleotide polymorphism in the protein tyrosine phosphatase N22 gene (PTPN22) is associated with Type 1 diabetes in a Chinese population

AIMS

To examine single nucleotide polymorphisms in the protein tyrosine phosphatase N22 gene (PTPN22) and to study their association with Type 1 diabetes in a Chinese cohort.

METHODS

Three hundred and sixty-four young patients with Type 1 diabetes and 719 healthy children were included in this case-controlled study. The genotypes of rs1217385, rs2488457 (-1123C>G), rs1217414, rs1217419, rs3765598 and rs2476601 (1858C>T) in the PTPN22 gene were determined using the SNaPshot method. Alleles, genotypes and haplotype frequencies were compared between patients with Type 1 diabetes and healthy control subjects. The association between single nucleotide polymorphisms and clinical traits/autoantibody status was also analysed.

RESULTS

The single nucleotide polymorphism, rs1217419, located in the second intron of the PTPN22 gene was associated with Type 1 diabetes (odds ratio 1.5, 95% CI 1.14-1.97, P = 0.003). An additional single nucleotide polymorphism, rs1217385, was also associated with Type 1 diabetes; however, the association was secondary to that of rs1217419. The previously reported single nucleotide polymorphism that is associated with Type 1 diabetes (-1123G>C) had only marginal association with Type 1 diabetes in our study. A marginal association was also identified between -1123G>C and glutamic acid decarboxylase autoantibody positivity in patients with Type 1 diabetes. There was no association between the single nucleotide polymorphism 1858C>T and Type 1 diabetes in our studied cohort.

CONCLUSIONS

Our study confirmed that PTPN22 is a gene that contributes to Type 1 diabetes susceptibility. The primary association occurs with single nucleotide polymorphism rs1217419 and there is clear heterogeneity of the association between PTPTN22 polymorphisms and Type 1 diabetes in a Chinese population compared with other populations.

Genes associated with SLE are targets of recent positive selection

The reasons for the ethnic disparities in the prevalence of systemic lupus erythematosus (SLE) and the relative high frequency of SLE risk alleles in the population are not fully understood. Population genetic factors such as natural selection alter allele frequencies over generations and may help explain the persistence of such common risk variants in the population and the differential risk of SLE. In order to better understand the genetic basis of SLE that might be due to natural selection, a total of 74 genomic regions with compelling evidence for association with SLE were tested for evidence of recent positive selection in the HapMap and HGDP populations, using population differentiation, allele frequency, and haplotype-based tests. Consistent signs of positive selection across different studies and statistical methods were observed at several SLE-associated loci, including PTPN22, TNFSF4, TET3-DGUOK, TNIP1, UHRF1BP1, BLK, and ITGAM genes. This study is the first to evaluate and report that several SLE-associated regions show signs of positive natural selection. These results provide corroborating evidence in support of recent positive selection as one mechanism underlying the elevated population frequency of SLE risk loci and supports future research that integrates signals of natural selection to help identify functional SLE risk alleles.

Signatures of natural selection on genetic variants affecting complex human traits

It has recently been hypothesized that polygenic adaptation, resulting in modest allele frequency changes at many loci, could be a major mechanism behind the adaptation of complex phenotypes in human populations. Here we leverage the large number of variants that have been identified through genome-wide association (GWA) studies to comprehensively study signatures of natural selection on genetic variants associated with complex traits. Using population differentiation based methods, such as FST and phylogenetic branch length analyses, we systematically examined nearly 1300 SNPs associated with 38 complex phenotypes. Instead of detecting selection signatures at individual variants, we aimed to identify combined evidence of natural selection by aggregating signals across many trait associated SNPs. Our results have revealed some general features of polygenic selection on complex traits associated variants. First, natural selection acting on standing variants associated with complex traits is a common phenomenon. Second, characteristics of selection for different polygenic traits vary both temporarily and geographically. Third, some studied traits (e.g. height and urate level) could have been the primary targets of selection, as indicated by the significant correlation between the effect sizes and the estimated strength of selection in the trait associated variants; however, for most traits, the allele frequency changes in trait associated variants might have been driven by the selection on other correlated phenotypes. Fourth, the changes in allele frequencies as a result of selection can be highly stochastic, such that, polygenic adaptation may accelerate differentiation in allele frequencies among populations, but generally does not produce predictable directional changes. Fifth, multiple mechanisms (pleiotropy, hitchhiking, etc) may act together to govern the changes in allele frequencies of genetic variants associated with complex traits.

The NeST long ncRNA controls microbial susceptibility and epigenetic activation of the interferon-γ locus

Long noncoding RNAs (lncRNAs) are increasingly appreciated as regulators of cell-specific gene expression. Here, an enhancer-like lncRNA termed NeST (nettoie Salmonella pas Theiler's [cleanup Salmonella not Theiler's]) is shown to be causal for all phenotypes conferred by murine viral susceptibility locus Tmevp3. This locus was defined by crosses between SJL/J and B10.S mice and contains several candidate genes, including NeST. The SJL/J-derived locus confers higher lncRNA expression, increased interferon-γ (IFN-γ) abundance in activated CD8(+) T cells, increased Theiler's virus persistence, and decreased Salmonella enterica pathogenesis. Transgenic expression of NeST lncRNA alone was sufficient to confer all phenotypes of the SJL/J locus. NeST RNA was found to bind WDR5, a component of the histone H3 lysine 4 methyltransferase complex, and to alter histone 3 methylation at the IFN-γ locus. Thus, this lncRNA regulates epigenetic marking of IFN-γ-encoding chromatin, expression of IFN-γ, and susceptibility to a viral and a bacterial pathogen.

Dissecting structural and nucleotide genome-wide variation in inbred Iberian pigs

BACKGROUND

In contrast to international pig breeds, the Iberian breed has not been admixed with Asian germplasm. This makes it an important model to study both domestication and relevance of Asian genes in the pig. Besides, Iberian pigs exhibit high meat quality as well as appetite and propensity to obesity. Here we provide a genome wide analysis of nucleotide and structural diversity in a reduced representation library from a pool (n=9 sows) and shotgun genomic sequence from a single sow of the highly inbred Guadyerbas strain. In the pool, we applied newly developed tools to account for the peculiarities of these data.

RESULTS

A total of 254,106 SNPs in the pool (79.6 Mb covered) and 643,783 in the Guadyerbas sow (1.47 Gb covered) were called. The nucleotide diversity (1.31x10-3 per bp in autosomes) is very similar to that reported in wild boar. A much lower than expected diversity in the X chromosome was confirmed (1.79x10-4 per bp in the individual and 5.83x10-4 per bp in the pool). A strong (0.70) correlation between recombination and variability was observed, but not with gene density or GC content. Multicopy regions affected about 4% of annotated pig genes in their entirety, and 2% of the genes partially. Genes within the lowest variability windows comprised interferon genes and, in chromosome X, genes involved in behavior like HTR2C or MCEP2. A modified Hudson-Kreitman-Aguadé test for pools also indicated an accelerated evolution in genes involved in behavior, as well as in spermatogenesis and in lipid metabolism.

CONCLUSIONS

This work illustrates the strength of current sequencing technologies to picture a comprehensive landscape of variability in livestock species, and to pinpoint regions containing genes potentially under selection. Among those genes, we report genes involved in behavior, including feeding behavior, and lipid metabolism. The pig X chromosome is an outlier in terms of nucleotide diversity, which suggests selective constraints. Our data further confirm the importance of structural variation in the species, including Iberian pigs, and allowed us to identify new paralogs for known gene families.

Cytokine gene polymorphisms and human autoimmune disease in the era of genome-wide association studies

Cytokine (receptor) genes have traditionally attracted great interest as plausible genetic risk factors for autoimmune disease. Since 2007, the implementation of genome-wide association studies has facilitated the robust identification of allelic variants in more than 35 cytokine loci as susceptibility factors for a wide variety of over 15 autoimmune disorders. In this review, we catalog the gene loci of interleukin, chemokine, and tumor necrosis factor receptor superfamily and ligands that have emerged as autoimmune risk factors. We examine recent progress made in the clarification of the functional mechanisms by which polymorphisms in the genes coding for interleukin-2 receptor alpha (IL2RA), IL7R, and IL23R may alter risk for autoimmune disease, and discuss opposite autoimmune risk alleles found, among others, at the IL10 locus.