Imputing Amino Acid Polymorphisms in Human Leukocyte Antigens

DNA sequence variation within human leukocyte antigen (HLA) genes mediate susceptibility to a wide range of human diseases. The complex genetic structure of the major histocompatibility complex (MHC) makes it difficult, however, to collect genotyping data in large cohorts. Long-range linkage disequilibrium between HLA loci and SNP markers across the major histocompatibility complex (MHC) region offers an alternative approach through imputation to interrogate HLA variation in existing GWAS data sets. Here we describe a computational strategy, SNP2HLA, to impute classical alleles and amino acid polymorphisms at class I (HLA-A, -B, -C) and class II (-DPA1, -DPB1, -DQA1, -DQB1, and -DRB1) loci. To characterize performance of SNP2HLA, we constructed two European ancestry reference panels, one based on data collected in HapMap-CEPH pedigrees (90 individuals) and another based on data collected by the Type 1 Diabetes Genetics Consortium (T1DGC, 5,225 individuals). We imputed HLA alleles in an independent data set from the British 1958 Birth Cohort (N = 918) with gold standard four-digit HLA types and SNPs genotyped using the Affymetrix GeneChip 500 K and Illumina Immunochip microarrays. We demonstrate that the sample size of the reference panel, rather than SNP density of the genotyping platform, is critical to achieve high imputation accuracy. Using the larger T1DGC reference panel, the average accuracy at four-digit resolution is 94.7% using the low-density Affymetrix GeneChip 500 K, and 96.7% using the high-density Illumina Immunochip. For amino acid polymorphisms within HLA genes, we achieve 98.6% and 99.3% accuracy using the Affymetrix GeneChip 500 K and Illumina Immunochip, respectively. Finally, we demonstrate how imputation and association testing at amino acid resolution can facilitate fine-mapping of primary MHC association signals, giving a specific example from type 1 diabetes

Fine-mapping and functional studies highlight potential causal variants for rheumatoid arthritis and type 1 diabetes

© 2018, The Author(s), under exclusive licence to Springer Nature America, Inc. To define potentially causal variants for autoimmune disease, we fine-mapped1,276 rheumatoid arthritis (11,475 cases, 15,870 controls)3and type 1 diabetes loci (9,334 cases, 11,111 controls)4. After sequencing 799 1-kilobase regulatory (H3K4me3) regions within these loci in 568 individuals, we observed accurate imputation for 89% of common variants. We defined credible sets of ≤5 causal variants at 5 rheumatoid arthritis and 10 type 1 diabetes loci. We identified potentially causal missense variants at DNASE1L3, PTPN22, SH2B3, and TYK2, and noncoding variants at MEG3, CD28–CTLA4, and IL2RA. We also identified potential candidate causal variants at SIRPG and TNFAIP3. Using functional assays, we confirmed allele-specific protein binding and differential enhancer activity for three variants: the CD28–CTLA4 rs117701653 SNP, MEG3 rs34552516 indel, and TNFAIP3 rs35926684 indel

Interpreting an apoptotic corpse as anti-inflammatory involves a chloride sensing pathway

Apoptotic cell clearance (efferocytosis) elicits an anti-inflammatory response by phagocytes, but the mechanisms that underlie this response are still being defined. Here, we uncover a chloride-sensing signalling pathway that controls both the phagocyte 'appetite' and its anti-inflammatory response. Efferocytosis transcriptionally altered the genes that encode the solute carrier (SLC) proteins SLC12A2 and SLC12A4. Interfering with SLC12A2 expression or function resulted in a significant increase in apoptotic corpse uptake per phagocyte, whereas the loss of SLC12A4 inhibited corpse uptake. In SLC12A2-deficient phagocytes, the canonical anti-inflammatory program was replaced by pro-inflammatory and oxidative-stress-associated gene programs. This 'switch' to pro-inflammatory sensing of apoptotic cells resulted from the disruption of the chloride-sensing pathway (and not due to corpse overload or poor degradation), including the chloride-sensing kinases WNK1, OSR1 and SPAK-which function upstream of SLC12A2-had a similar effect on efferocytosis. Collectively, the WNK1-OSR1-SPAK-SLC12A2/SLC12A4 chloride-sensing pathway and chloride flux in phagocytes are key modifiers of the manner in which phagocytes interpret the engulfed apoptotic corpse

A method for gene-based pathway analysis using genomewide association study summary statistics reveals nine new type 1 diabetes associations.

Pathway analysis can complement point-wise single nucleotide polymorphism (SNP) analysis in exploring genomewide association study (GWAS) data to identify specific disease-associated genes that can be candidate causal genes. We propose a straightforward methodology that can be used for conducting a gene-based pathway analysis using summary GWAS statistics in combination with widely available reference genotype data. We used this method to perform a gene-based pathway analysis of a type 1 diabetes (T1D) meta-analysis GWAS (of 7,514 cases and 9,045 controls). An important feature of the conducted analysis is the removal of the major histocompatibility complex gene region, the major genetic risk factor for T1D. Thirty-one of the 1,583 (2%) tested pathways were identified to be enriched for association with T1D at a 5% false discovery rate. We analyzed these 31 pathways and their genes to identify SNPs in or near these pathway genes that showed potentially novel association with T1D and attempted to replicate the association of 22 SNPs in additional samples. Replication P-values were skewed (P=9.85×10-11) with 12 of the 22 SNPs showing P<0.05. Support, including replication evidence, was obtained for nine T1D associated variants in genes ITGB7 (rs11170466, P=7.86×10-9), NRP1 (rs722988, 4.88×10-8), BAD (rs694739, 2.37×10-7), CTSB (rs1296023, 2.79×10-7), FYN (rs11964650, P=5.60×10-7), UBE2G1 (rs9906760, 5.08×10-7), MAP3K14 (rs17759555, 9.67×10-7), ITGB1 (rs1557150, 1.93×10-6), and IL7R (rs1445898, 2.76×10-6). The proposed methodology can be applied to other GWAS datasets for which only summary level data are available.This is the final version. It was first published by Wiley at http://onlinelibrary.wiley.com/doi/10.1002/gepi.21853/abstract

Context-dependent compensation among phosphatidylserine-recognition receptors

Phagocytes express multiple phosphatidylserine (PtdSer) receptors that recognize apoptotic cells. It is unknown whether these receptors are interchangeable or if they play unique roles during cell clearance. Loss of the PtdSer receptor Mertk is associated with apoptotic corpse accumulation in the testes and degeneration of photoreceptors in the eye. Both phenotypes are linked to impaired phagocytosis by specialized phagocytes: Sertoli cells and the retinal pigmented epithelium (RPE). Here, we overexpressed the PtdSer receptor BAI1 in mice lacking MerTK (Mertk(-/-) Bai1(Tg)) to evaluate PtdSer receptor compensation in vivo. While Bai1 overexpression rescues clearance of apoptotic germ cells in the testes of Mertk(-/-) mice it fails to enhance RPE phagocytosis or prevent photoreceptor degeneration. To determine why MerTK is critical to RPE function, we examined visual cycle intermediates and performed unbiased RNAseq analysis of RPE from Mertk(+/+) and Mertk(-/-) mice. Prior to the onset of photoreceptor degeneration, Mertk(-/-) mice had less accumulation of retinyl esters and dysregulation of a striking array of genes, including genes related to phagocytosis, metabolism, and retinal disease in humans. Collectively, these experiments establish that not all phagocytic receptors are functionally equal, and that compensation among specific engulfment receptors is context and tissue dependent

A genome-wide assessment of the role of untagged copy number variants in type 1 diabetes.

Genome-wide association studies (GWAS) for type 1 diabetes (T1D) have successfully identified more than 40 independent T1D associated tagging single nucleotide polymorphisms (SNPs). However, owing to technical limitations of copy number variants (CNVs) genotyping assays, the assessment of the role of CNVs has been limited to the subset of these in high linkage disequilibrium with tag SNPs. The contribution of untagged CNVs, often multi-allelic and difficult to genotype using existing assays, to the heritability of T1D remains an open question. To investigate this issue, we designed a custom comparative genetic hybridization array (aCGH) specifically designed to assay untagged CNV loci identified from a variety of sources. To overcome the technical limitations of the case control design for this class of CNVs, we genotyped the Type 1 Diabetes Genetics Consortium (T1DGC) family resource (representing 3,903 transmissions from parents to affected offspring) and used an association testing strategy that does not necessitate obtaining discrete genotypes. Our design targeted 4,309 CNVs, of which 3,410 passed stringent quality control filters. As a positive control, the scan confirmed the known T1D association at the INS locus by direct typing of the 5' variable number of tandem repeat (VNTR) locus. Our results clarify the fact that the disease association is indistinguishable from the two main polymorphic allele classes of the INS VNTR, class I-and class III. We also identified novel technical artifacts resulting into spurious associations at the somatically rearranging loci, T cell receptor, TCRA/TCRD and TCRB, and Immunoglobulin heavy chain, IGH, loci on chromosomes 14q11.2, 7q34 and 14q32.33, respectively. However, our data did not identify novel T1D loci. Our results do not support a major role of untagged CNVs in T1D heritability

Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping.

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.We acknowledge use of DNA from The UK Blood Services collection of Common Controls (UKBS-CC collection), which is funded by the Wellcome Trust grant 076113/C/04/Z and by the USA National Institute for Health Research program grant to the National Health Service Blood and Transplant (RP-PG-0310-1002). We acknowledge the use of DNA from the British 1958 Birth Cohort collection, which is funded by the UK Medical Research Council grant G0000934 and the Wellcome Trust grant 068545/Z/02. This research utilized resources provided by the Type 1 Diabetes Genetics Consortium, a collaborative clinical study sponsored by the National Institute of Diabetes and Digestive and Kidney Diseases, the National Institute of Allergy and Infectious Diseases, the National Human Genome Research Institute, the National Institute of Child Health and Human Development and the JDRF and is supported by the USA National Institutes of Health grant U01-DK062418. The JDRF/Wellcome Trust Diabetes and Inflammation Laboratory is funded by the JDRF (9-2011-253), the Wellcome Trust (091157) and the National Institute for Health Research Cambridge Biomedical Centre. The research leading to these results has received funding from the European Union's 7th Framework Programme (FP7/2007-2013) under grant agreement no.241447 (NAIMIT). The Cambridge Institute for Medical Research (CIMR) is in receipt of a Wellcome Trust Strategic Award (100140). CW is supported by the Wellcome Trust (089989). We acknowledge the National Institute for Health Research Cambridge Biomedical Research Centre for funding.This is the final version of the article. It first appeared from PLOS via http://dx.doi.org/10.1371/journal.pgen.100527

Additive and interaction effects at three amino acid positions in HLA-DQ and HLA-DR molecules drive type 1 diabetes risk.

Variation in the human leukocyte antigen (HLA) genes accounts for one-half of the genetic risk in type 1 diabetes (T1D). Amino acid changes in the HLA-DR and HLA-DQ molecules mediate most of the risk, but extensive linkage disequilibrium complicates the localization of independent effects. Using 18,832 case-control samples, we localized the signal to 3 amino acid positions in HLA-DQ and HLA-DR. HLA-DQβ1 position 57 (previously known; P = 1 × 10(-1,355)) by itself explained 15.2% of the total phenotypic variance. Independent effects at HLA-DRβ1 positions 13 (P = 1 × 10(-721)) and 71 (P = 1 × 10(-95)) increased the proportion of variance explained to 26.9%. The three positions together explained 90% of the phenotypic variance in the HLA-DRB1-HLA-DQA1-HLA-DQB1 locus. Additionally, we observed significant interactions for 11 of 21 pairs of common HLA-DRB1-HLA-DQA1-HLA-DQB1 haplotypes (P = 1.6 × 10(-64)). HLA-DRβ1 positions 13 and 71 implicate the P4 pocket in the antigen-binding groove, thus pointing to another critical protein structure for T1D risk, in addition to the HLA-DQ P9 pocket.This research utilizes resources provided by the Type 1 Diabetes Genetics Consortium, a collaborative clinical study sponsored by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), National Institute of Allergy and Infectious Diseases (NIAID), National Human Genome Research Institute (NHGRI), National Institute of Child Health and Human Development (NICHD), and Juvenile Diabetes Research Foundation International (JDRF) and supported by U01 DK062418. This work is supported in part by funding from the National Institutes of Health (5R01AR062886-02 (PIdB), 1R01AR063759 (SR), 5U01GM092691-05 (SR), 1UH2AR067677-01 (SR), R01AR065183 (PIWdB)), a Doris Duke Clinical Scientist Development Award (SR), the Wellcome Trust (JAT) and the National Institute for Health Research (JAT and JMMH), and a Vernieuwingsimpuls VIDI Award (016.126.354) from the Netherlands Organization for Scientific Research (PIWdB). TLL was supported by the German Research Foundation (LE 2593/1-1 and LE 2593/2-1).This is the accepted manuscript. The final version is available at http://www.nature.com/ng/journal/v47/n8/full/ng.3353.html

Cytomegalovirus infection enhances the immune response to influenza

Cytomegalovirus (CMV) is a beta-herpes virus present in a latent form in most people worldwide. In immunosuppressed individuals, CMV can reactivate and cause serious clinical complications, but the effect of the latent state on healthy people remains elusive. We undertook a systems approach to understand the differences between seropositive and negative subjects and measured hundreds of immune system components from blood samples including cytokines and chemokines, immune cell phenotyping, gene expression, ex vivo cell responses to cytokine stimuli and the antibody response to seasonal influenza vaccination. As expected, we found decreased responses to vaccination and an overall down-regulation of immune components in aged individuals regardless of CMV serostatus. In contrast, CMV-infected young adults exhibited an overall up-regulation of immune components including enhanced antibody responses to influenza vaccination, increased CD8+ T cell sensitivity, and elevated levels of circulating IFN-γ compared to uninfected individuals. Experiments with young mice infected with murine CMV also showed significant protection from an influenza virus challenge compared with uninfected animals, although this effect declined with time. These data show that CMV and its murine equivalent can have a beneficial effect on the immune response of young, healthy individuals, which may explain the continued coexistence of CMV and mammals throughout their evolution