Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease.

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

Host genetic control of natural killer cell diversity revealed in the Collaborative Cross.

Natural killer (NK) cells are innate effectors armed with cytotoxic and cytokine-secreting capacities whose spontaneous antitumor activity is key to numerous immunotherapeutic strategies. However, current mouse models fail to mirror the extensive immune system variation that exists in the human population which may impact on NK cell-based therapies. We performed a comprehensive profiling of NK cells in the Collaborative Cross (CC), a collection of novel recombinant inbred mouse strains whose genetic diversity matches that of humans, thereby providing a unique and highly diverse small animal model for the study of immune variation. We demonstrate that NK cells from CC strains displayed a breadth of phenotypic and functional variation reminiscent of that reported for humans with regards to cell numbers, key marker expression, and functional capacities. We took advantage of the vast genetic diversity of the CC and identified nine genomic loci through quantitative trait locus mapping driving these phenotypic variations. SNP haplotype patterns and variant effect analyses identified candidate genes associated with lung NK cell numbers, frequencies of CD94+ NK cells, and expression levels of NKp46. Thus, we demonstrate that the CC represents an outstanding resource to study NK cell diversity and its regulation by host genetics.

CancerImmunityQTL: a database to systematically evaluate the impact of genetic variants on immune infiltration in human cancer.

Tumor-infiltrating immune cells as integral component of the tumor microenvironment are associated with tumor progress, prognosis and responses to immunotherapy. Genetic variants have been demonstrated to impact tumor-infiltrating, underscoring the heritable character of immune landscape. Therefore, identification of immunity quantitative trait loci (immunQTLs), which evaluate the effect of genetic variants on immune cells infiltration, might present a critical step toward fully understanding the contribution of genetic variants in tumor development. Although emerging studies have demonstrated the determinants of germline variants on immune infiltration, no database has yet been developed to systematically analyze immunQTLs across multiple cancer types. Using genotype data from TCGA database and immune cell fractions estimated by CIBERSORT, we developed a computational pipeline to identify immunQTLs in 33 cancer types. A total of 913 immunQTLs across different cancer types were identified. Among them, 5 immunQTLs are associated with patient overall survival. Furthermore, by integrating immunQTLs with GWAS data, we identified 527 immunQTLs overlapping with known GWAS linkage disequilibrium regions. Finally, we constructed a user-friendly database, CancerImmunityQTL (http://www.cancerimmunityqtl-hust.com/) for users to browse, search and download data of interest. This database provides an informative resource to understand the germline determinants of immune infiltration in human cancer and benefit from personalized cancer immunotherapy.

Tensor decomposition of stimulated monocyte and macrophage gene expression profiles identifies neurodegenerative disease-specific trans-eQTLs.

Recent human genetic studies suggest that cells of the innate immune system have a primary role in the pathogenesis of neurodegenerative diseases. However, the results from these studies often do not elucidate how the genetic variants affect the biology of these cells to modulate disease risk. Here, we applied a tensor decomposition method to uncover disease associated gene networks linked to distal genetic variation in stimulated human monocyte and macrophage gene expression profiles. We report robust evidence that some disease associated genetic variants affect the expression of multiple genes in trans. These include a Parkinson's disease locus influencing the expression of genes mediated by a protease that controls lysosomal function, and Alzheimer's disease loci influencing the expression of genes involved in type 1 interferon signaling, myeloid phagocytosis, and complement cascade pathways. Overall, we uncover gene networks in induced innate immune cells linked to disease associated genetic variants, which may help elucidate the underlying biology of disease.

Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis.

Asthma is one of the most common chronic diseases affecting both children and adults. We report a genome-wide association meta-analysis of 69,189 cases and 702,199 controls from Iceland and UK biobank. We find 88 asthma risk variants at 56 loci, 19 previously unreported, and evaluate their effect on other asthma and allergic phenotypes. Of special interest are two low frequency variants associated with protection against asthma; a missense variant in TNFRSF8 and 3' UTR variant in TGFBR1. Functional studies show that the TNFRSF8 variant reduces TNFRSF8 expression both on cell surface and in soluble form, acting as loss of function. eQTL analysis suggests that the TGFBR1 variant acts through gain of function and together with an intronic variant in a downstream gene, SMAD3, points to defective TGFßR1 signaling as one of the biological perturbations increasing asthma risk. Our results increase the number of asthma variants and implicate genes with known role in T cell regulation, inflammation and airway remodeling in asthma pathogenesis.

An Allele-Specific Functional SNP Associated with Two Systemic Autoimmune Diseases Modulates IRF5 Expression by Long-Range Chromatin Loop Formation.

Both systemic lupus erythematosus (SLE) and systemic sclerosis (SSc) are autoimmune diseases sharing similar genetic backgrounds. Genome-wide association studies have constantly disclosed numerous genetic variants conferring to both disease risks at 7q32.1, but the functional mechanisms underlying them are still largely unknown. Through a series of bioinformatics and functional analyses, we prioritized a potential independent functional single-nucleotide polymorphism (rs13239597) within TNPO3 promoter region, residing in a putative enhancer element and validated that IRF5 is the distal target gene (∼118 kb) of rs13239597, which is a key regulator involved in pathogenic autoantibody dysregulation, increasing risk of both SLE and SSc. We experimentally validated the long-range chromatin interactions between rs13239597 and IRF5 using chromosome conformation capture assay. We further demonstrated that rs13239597-A acted as an allele-specific enhancer regulating IRF5 expression, independently of TNPO3 by using dual-luciferase reporter assays and CRISPR-Cas9. Particularly, the transcription factor EVI1 could preferentially bind to rs13239597-A allele and increase the enhancer activity to regulate IRF5 expression. Taken together, our results uncovered a mechanistic insight of a noncoding functional variant acting as an allele-specific distal enhancer to directly modulate IRF5 expression, which might obligate in understanding of complex genetic architectures of SLE and SSc pathogenesis.

Potentially functional genetic variants in the complement-related immunity gene-set are associated with non-small cell lung cancer survival.

The complement system plays an important role in the innate and adaptive immunity, complement components mediate tumor cytolysis of antibody-based immunotherapy, and complement activation in the tumor microenvironment may promote tumor progression or inhibition, depending on the mechanism of action. In the present study, we conducted a two-phase analysis of two independently published genome-wide association studies (GWASs) for associations between genetic variants in a complement-related immunity gene-set and overall survival of non-small cell lung cancer (NSCLC). The GWAS dataset from Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial was used as the discovery, and multivariate Cox proportional hazards regression with false-positive report probability for multiple test corrections were performed to evaluate associations between 14,699 single-nucleotide polymorphisms (SNPs) in 111 genes and survival of 1,185 NSCLC patients. The identified significant SNPs in a single-locus analysis were further validated with 984 NSCLC patients in the GWAS dataset from the Harvard Lung Cancer Susceptibility (HLCS) Study. The results showed that two independent, potentially functional SNPs in two genes (VWF rs73049469 and ITGB2 rs3788142) were significantly associated with NSCLC survival, with a combined hazards ratio (HR) of 1.22 [95% confidence interval (CI) = 1.07-1.40, P = 0.002] and 1.16 (1.07-1.27, 6.45 × 10-4 ), respectively. Finally, we performed expression quantitative trait loci (eQTL) analysis and found that survival-associated genotypes of VWF rs73049469 were also significantly associated with mRNA expression levels of the gene. These results indicated that genetic variants of the complement-related immunity genes might be predictors of NSCLC survival, particularly for the short-term survival, possibly by modulating the expression of genes involved in the host immunity.

Reduced monocyte and macrophage TNFSF15/TL1A expression is associated with susceptibility to inflammatory bowel disease.

Chronic inflammation in inflammatory bowel disease (IBD) results from a breakdown of intestinal immune homeostasis and compromise of the intestinal barrier. Genome-wide association studies have identified over 200 genetic loci associated with risk for IBD, but the functional mechanisms of most of these genetic variants remain unknown. Polymorphisms at the TNFSF15 locus, which encodes the TNF superfamily cytokine commonly known as TL1A, are associated with susceptibility to IBD in multiple ethnic groups. In a wide variety of murine models of inflammation including models of IBD, TNFSF15 promotes immunopathology by signaling through its receptor DR3. Such evidence has led to the hypothesis that expression of this lymphocyte costimulatory cytokine increases risk for IBD. In contrast, here we show that the IBD-risk haplotype at TNFSF15 is associated with decreased expression of the gene by peripheral blood monocytes in both healthy volunteers and IBD patients. This association persists under various stimulation conditions at both the RNA and protein levels and is maintained after macrophage differentiation. Utilizing a "recall-by-genotype" bioresource for allele-specific expression measurements in a functional fine-mapping assay, we localize the polymorphism controlling TNFSF15 expression to the regulatory region upstream of the gene. Through a T cell costimulation assay, we demonstrate that genetically regulated TNFSF15 has functional relevance. These findings indicate that genetically enhanced expression of TNFSF15 in specific cell types may confer protection against the development of IBD.

Mycobacterial infection induces a specific human innate immune response.

The innate immune system provides the first response to infection and is now recognized to be partially pathogen-specific. Mycobacterium tuberculosis (MTB) is able to subvert the innate immune response and survive inside macrophages. Curiously, only 5-10% of otherwise healthy individuals infected with MTB develop active tuberculosis (TB). We do not yet understand the genetic basis underlying this individual-specific susceptibility. Moreover, we still do not know which properties of the innate immune response are specific to MTB infection. To identify immune responses that are specific to MTB, we infected macrophages with eight different bacteria, including different MTB strains and related mycobacteria, and studied their transcriptional response. We identified a novel subset of genes whose regulation was affected specifically by infection with mycobacteria. This subset includes genes involved in phagosome maturation, superoxide production, response to vitamin D, macrophage chemotaxis, and sialic acid synthesis. We suggest that genetic variants that affect the function or regulation of these genes should be considered candidate loci for explaining TB susceptibility.

Dense genotyping of immune-related susceptibility loci reveals new insights into the genetics of psoriatic arthritis.

Psoriatic arthritis (PsA) is a chronic inflammatory arthritis associated with psoriasis and, despite the larger estimated heritability for PsA, the majority of genetic susceptibility loci identified to date are shared with psoriasis. Here, we present results from a case-control association study on 1,962 PsA patients and 8,923 controls using the Immunochip genotyping array. We identify eight loci passing genome-wide significance, secondary independent effects at three loci and a distinct PsA-specific variant at the IL23R locus. We report two novel loci and evidence of a novel PsA-specific association at chromosome 5q31. Imputation of classical HLA alleles, amino acids and SNPs across the MHC region highlights three independent associations to class I genes. Finally, we find an enrichment of associated variants to markers of open chromatin in CD8(+) memory primary T cells. This study identifies key insights into the genetics of PsA that could begin to explain fundamental differences between psoriasis and PsA.

Follicular lymphoma-protective HLA class II variants correlate with increased HLA-DQB1 protein expression.

Multiple follicular lymphoma (FL) susceptibility single-nucleotide polymorphisms in the human leukocyte antigen (HLA) class I and II regions have been identified, including rs6457327, rs3117222, rs2647012, rs10484561, rs9268853 and rs2621416. Here we validated previous expression quantitative trait loci results with real-time reverse transcription quantitative PCR and investigated protein expression in B-lymphoblastoid cell lines and primary dendritic cells using flow cytometry, cell-based enzyme-linked immunosorbent assay and western blotting. We confirmed that FL-protective rs2647012-linked variants, in high linkage disequilibrium with the extended haplotype DRB1*15:01-DQA1*01:02-DQB1*06:02, correlate with increased HLA-DQB1 expression. This association remained significant at the protein level and was reproducible across different cell types. We also found that differences in HLA-DQB1 expression were not related to changes in activation markers or class II, major histocompatibility complex, transactivator expression, suggesting the role of an alternative regulatory mechanism. However, functional analysis using RegulomeDB did not reveal any relevant regulatory candidates. Future studies should focus on the clinical relevance of increased HLA-DQB1 protein expression facilitating tumor cell removal through increased immune surveillance.

The genetics of autoimmune thyroiditis: the first decade.

Most of our current understanding of the genetic predisposition to autoimmune disease can be traced to experiments performed in the decade from 1971 to 1981. Chella David was a key contributor to this research. Many of these early steps came from studies of experimental autoimmune thyroiditis. This model has been especially valuable because essentially the same disease can occur spontaneously in selected strains of animals or can be induced by deliberate immunization. From a genetic point of view, the disease has been investigated in three different species: mice, rats and chickens. The same antigen, thyroglobulin, initiates the disease in all three species. Among the main discoveries were the relationship of autoimmune disease to the major histocompatibility complex (MHC), the interplay of different subregions within the MHC in promoting or retarding development of disease, the differing roles of MHC class II and MHC I class genes in induction and effector phases, respectively, and the cumulative effect of non-MHC genes, each of which represents a small addition to overall susceptibility. Other experiments revealed that genetic differences in thyroglobulin allotypes influence susceptibility to thyroiditis. Thyroid glands differed in different strains in vulnerability to passive transfer of antibody. The first evidence of modulatory genes on the sex-related X chromosome emerged. All of these genetic findings were concurrently translated to the human disease, Hashimoto's thyroiditis, where thyroglobulin is also the initiating antigen.

Severity of innate immune-mediated colitis is controlled by the cytokine deficiency-induced colitis susceptibility-1 (Cdcs1) locus.

Genetic modifier loci influence the phenotypic expression of many Mendelian traits; insight into disease pathogenesis gained from their identification in animal disease models may impact the treatment of human multigenic disorders. We previously described an innate immune-driven model of spontaneous ulcerative colitis in T-bet(-/-).Rag2(-/-) double-deficient mice that resembles human ulcerative colitis. On a BALB/c background, this disease is highly penetrant and results in the development of colorectal cancer. However, we observed that colitis in T-bet(-/-).Rag2(-/-) mice on a C57BL/6 background was significantly less severe. Quantitative trait locus analysis using an N2 backcross strategy revealed a single major quantitative trait locus on chromosome 3 that mapped to the Cdcs1 (cytokine deficiency-induced colitis susceptibility-1) locus previously identified in the Il10(-/-) and Gnai2(-/-) colitis models. Congenic introduction of the susceptible Cdcs1 interval from C3H/He into the C57BL/6 background restored colitis severity. Bone marrow reconstitution experiments further mapped the effect of host genetics on disease severity to the hematopoietic compartment. There were distinct differences in the expression of several Cdcs1 genes in bone marrow-derived dendritic cells from Cdcs1 congenic mice. We conclude that the Cdcs1 locus controls colitis severity in T-bet(-/-).Rag2(-/-) mice through innate immune cells.

Genetic control of IL-1 beta production and inflammatory response by the mouse Irm1 locus.

Genome-wide linkage analysis using single nucleotide polymorphism arrays was carried out in pedigrees of mice differing in the extent of acute inflammatory response (AIRmax or AIRmin). The AIR phenotype was determined by quantifying the number of infiltrating cells in the 24-h exudate induced by Biogel P-100 s.c. injection and by ex vivo IL-1beta production by leukocytes stimulated with LPS and ATP. We mapped the major inflammatory response modulator 1 locus on chromosome 7, at the 1-logarithm of odds (LOD) confidence interval from 116.75 to 139.75 Mb, linked to the number of infiltrating cells (LOD = 3.61) through the production of IL-1beta (LOD = 9.35). Of several interesting candidate genes mapping to the inflammatory response modulator 1 locus, 28 of these were differentially expressed in the bone marrow of AIRmax and AIRmin mice. These findings represent a step toward the identification of the genes underlying this complex phenotype.

TNF production in macrophages is genetically determined and regulates inflammatory disease in rats.

Dysregulation of TNF is an important pathophysiological phenotype for many diseases. Recently, certain genetically regulated loci have been identified to regulate several inflammatory diseases. We hypothesized that a region on rat chromosome 4 known to regulate experimental autoimmune encephalomyelitis, experimental arthritis and experimental autoimmune neuritis harbors a gene regulating central inflammatory molecules, such as TNF. We therefore mapped TNF production using linkage analysis in the 12th generation of an advanced intercross line between DA and PVG.AV1 rats, which differ in susceptibility to several inflammatory conditions. A single TNF-regulating quantitative trait locus with a logarithm of odds score of 6.2 was identified and its biological effect was confirmed in a congenic rat strain. The profound TNF regulation mapped in congenic strains to the macrophage population. Several TLR signaling cascades led to the same reduced proinflammatory phenotype in congenic macrophages, indicating control of a convergence point for innate inflammatory activity. The decreased TNF potential and reduced proinflammatory macrophage phenotype in congenic rats was also associated with reduced clinical severity in experimental autoimmune encephalomyelitis, pristane-induced arthritis and sepsis experimental models. Determination of genes and mechanisms involved in this genetically determined TNF regulation will be valuable in understanding disease pathogenesis and aid treatment development.

The major locus for mouse adenovirus susceptibility maps to genes of the hematopoietic cell surface-expressed LY6 family.

Susceptibility to mouse adenovirus type 1 is associated with the major quantitative trait locus Msq1. Msq1 was originally mapped to a 13-Mb region of mouse chromosome (Chr) 15 in crosses between SJL/J and BALB/cJ inbred mice. We have now narrowed Msq1 to a 0.75-Mb interval from 74.68 to 75.43 Mb, defined by two anonymous markers, rs8259436 and D15Spn14, using data from 1396 backcross mice. The critical interval includes 14 Ly6 or Ly6-related genes, including Ly6a (encoding Sca-1/TAP), Ly6e (Sca-2/Tsa1), Ly6g (Gr-1), and gpihbp1 (GPI-anchored high-density lipoprotein-binding protein 1), as well as the gene encoding an aldosterone synthase (Cyp11b2). The Ly6 family members are attractive candidates for virus susceptibility genes because their products are GPI-anchored membrane proteins expressed on lymphoid and myeloid cells, with proposed functions in cell adhesion and cell signaling. To determine interstrain variation in susceptibility and produce additional resources for cloning Msq1, we assayed the susceptibility phenotype of four previously untested inbred mouse strains. Susceptibility of strain 129S6/SvEvTac was subsequently localized to the Ly6 complex region, using polymorphic genetic markers on Chr 15 in a population of 271 (129S6/SvEvTac x BALB/cJ)F(1) x BALB/cJ backcross mice. We identified a major 129S6/SvEvTac susceptibility allele, Msq1(129S6), on Chr 15 in the same region as Msq1(SJL). The results indicate that a major host factor in mouse adenovirus type 1 susceptibility is likely to be a member of the Ly6 gene family.

Genetic control of severe egg-induced immunopathology and IL-17 production in murine schistosomiasis.

Infection with the trematode parasite Schistosoma mansoni results in a distinct heterogeneity of disease severity, both in humans and in an experimental mouse model. Severe disease is characterized by pronounced hepatic egg-induced granulomatous inflammation in a proinflammatory cytokine environment, whereas mild disease corresponds with reduced hepatic inflammation in a Th2 skewed cytokine environment. This marked heterogeneity indicates that genetic differences play a significant role in disease development, yet little is known about the genetic basis of dissimilar immunopathology. To investigate the role of genetic susceptibility in murine schistosomiasis, quantitative trait loci analysis was performed on F(2) progeny derived from SJL/J and C57BL/6 mice, which develop severe and mild pathology, respectively. In this study, we show that severe liver pathology in F(2) mice 7 wk after infection significantly correlated with an increase in the production of the proinflammatory cytokines IL-17, IFN-gamma, and TNF-alpha by schistosome egg Ag-stimulated mesenteric lymph node cells. Quantitative trait loci analysis identified several genetic intervals controlling immunopathology as well as IL-17 and IFN-gamma production. Egg granuloma size exhibited significant linkage to two loci, D4Mit203 and D17Mit82, both of which were inherited in a BL/6 dominant manner. Furthermore, a significant reduction of hepatic granulomatous inflammation and IL-17 production in interval-specific congenic mice demonstrated that the two identified genetic loci have a decisive effect on the development of immunopathology in murine schistosomiasis.

Distribution of BoLA-DRB3 allelic frequencies and identification of a new allele in the iranian cattle breed sistani (Bos indicus).

The distribution of the frequencies of BoLA-DRB3 gene alleles in the Iranian cattle breed Sistani was studied by the PCR-RFLP ("hemi-nested") assay using restriction endonucleases RsaI, HaeIII and BstYI. In the examined cattle breed (65 animals) 32 alleles have been identified one of which being described for the first time (6.15% frequency). The nucleotide sequence of the polymorphic region of exon 2 of this allele has been determined and submitted in the GeneBank database under accession number DQ486519. The submitted sequence has maximum homology (92%) with the previously described sequence DRB3-mRNA from Bos indicus (AccN X79346) and differs from it by 24 nucleotide substitutions which result in 16 amino acid substitutions. The peptide (on the basis of the reconstructed amino acid sequence) has 89% identity to the sequence encoded by the BIDRBF 188 locus (Bos indicus). The results obtained permit the sequence described by us to be considered as a new allele of the BoLA-DRB3 gene (DRB3.2**X). The total frequency of the main six alleles (DRB3.2*X, *10, *11, *20, *34 and *X) occurring with a frequency of over 5% is about 60% in Iranian Sistani cattle. Fifteen alleles have <1% frequency. The highest frequency was observed for DRB3.2*8 allele (21.54%) like in other previously described breeds of Bos indicus (up to 23.07%). The Iranian breed Sistani has a high level of similarity by the spectrum of BoLA-DRB3 alleles and their frequencies to other Bos indicus breeds and significantly differs by these criteria from the Bos taurus breeds. The Iranian Sistani herd under study includes alleles associated with to resistance to leukemia (DRB3.2*ll and *23) and to different forms of mastitis (DRB3.2*2, *7, *11, *23 and *24) although their frequencies are low (from 0.77 to 5.37%). On the whole, a high level of diversity of BoLA-DRB3 gene alleles and the availability of alleles associated with resistance to different diseases makes this breed of interest for breeding practice.