Jennie H. M. Yang

Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes

The Wellcome Trust Case Control Consortium (WTCCC) primary genome-wide association (GWA) scan on seven diseases, including the multifactorial autoimmune disease type 1 diabetes (T1D), shows associations at P < 5 × 10−7 between T1D and six chromosome regions: 12q24, 12q13, 16p13, 18p11, 12p13 and 4q27. Here, we attempted to validate these and six other top findings in 4,000 individuals with T1D, 5,000 controls and 2,997 family trios independent of the WTCCC study. We confirmed unequivocally the associations of 12q24, 12q13, 16p13 and 18p11 (Pfollow-up ≤ 1.35 × 10−9; Poverall ≤ 1.15 × 10−14), leaving eight regions with small effects or false-positive associations. We also obtained evidence for chromosome 18q22 (Poverall = 1.38 × 10−8) from a GWA study of nonsynonymous SNPs. Several regions, including 18q22 and 18p11, showed association with autoimmune thyroid disease. This study increases the number of T1D loci with compelling evidence from six to at least ten.

Shared and Distinct Genetic Variants in Type 1 Diabetes and Celiac Disease

BACKGROUND Two inflammatory disorders, type 1 diabetes and celiac disease, cosegregate in populations, suggesting a common genetic origin. Since both diseases are associated with the HLA class II genes on chromosome 6p21, we tested whether non-HLA loci are shared. METHODS We evaluated the association between type 1 diabetes and eight loci related to the risk of celiac disease by genotyping and statistical analyses of DNA samples from 8064 patients with type 1 diabetes, 9339 control subjects, and 2828 families providing 3064 parent-child trios (consisting of an affected child and both biologic parents). We also investigated 18 loci associated with type 1 diabetes in 2560 patients with celiac disease and 9339 control subjects. RESULTS Three celiac disease loci--RGS1 on chromosome 1q31, IL18RAP on chromosome 2q12, and TAGAP on chromosome 6q25--were associated with type 1 diabetes (P<1.00x10(-4)). The 32-bp insertion-deletion variant on chromosome 3p21 was newly identified as a type 1 diabetes locus (P=1.81x10(-8)) and was also associated with celiac disease, along with PTPN2 on chromosome 18p11 and CTLA4 on chromosome 2q33, bringing the total number of loci with evidence of a shared association to seven, including SH2B3 on chromosome 12q24. The effects of the IL18RAP and TAGAP alleles confer protection in type 1 diabetes and susceptibility in celiac disease. Loci with distinct effects in the two diseases included INS on chromosome 11p15, IL2RA on chromosome 10p15, and PTPN22 on chromosome 1p13 in type 1 diabetes and IL12A on 3q25 and LPP on 3q28 in celiac disease. CONCLUSIONS A genetic susceptibility to both type 1 diabetes and celiac disease shares common alleles. These data suggest that common biologic mechanisms, such as autoimmunity-related tissue damage and intolerance to dietary antigens, may be etiologic features of both diseases.

Cell-specific protein phenotypes for the autoimmune locus IL2RA using a genotype-selectable human bioresource

Genome-wide association studies (GWAS) have identified over 300 regions associated with more than 70 common diseases. However, identifying causal genes within an associated region remains a major challenge. One approach to resolving causal genes is through the dissection of gene-phenotype correlations. Here we use polychromatic flow cytometry to show that differences in surface expression of the human interleukin-2 (IL-2) receptor alpha (IL2RA, or CD25) protein are restricted to particular immune cell types and correlate with several haplotypes in the IL2RA region that have previously been associated with two autoimmune diseases, type 1 diabetes (T1D) and multiple sclerosis. We confirm our strongest gene-phenotype correlation at the RNA level by allele-specific expression (ASE). We also define key parameters for the design and implementation of post-GWAS gene-phenotype investigations and demonstrate the usefulness of a large bioresource of genotype-selectable normal donors from whom fresh, primary cells can be analyzed.

Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing

Many disease-associated variants identified by genome-wide association (GWA) studies are expected to regulate gene expression. Allele-specific expression (ASE) quantifies transcription from both haplotypes using individuals heterozygous at tested SNPs. We performed deep human transcriptome-wide resequencing (RNA-seq) for ASE analysis and expression quantitative trait locus discovery. We resequenced double poly(A)-selected RNA from primary CD4+ T cells (n = 4 individuals, both activated and untreated conditions) and developed tools for paired-end RNA-seq alignment and ASE analysis. We generated an average of 20 million uniquely mapping 45 base reads per sample. We obtained sufficient read depth to test 1371 unique transcripts for ASE. Multiple biases inflate the false discovery rate which we estimate to be ∼50% for random SNPs. However, after controlling for these biases and considering the subset of SNPs that pass HapMap QC, 4.6% of heterozygous SNP-sample pairs show evidence of imbalance (P < 0.001). We validated four findings by both bacterial cloning and Sanger sequencing assays. We also found convincing evidence for allelic imbalance at multiple reporter exonic SNPs in CD6 for two samples heterozygous at the multiple sclerosis-associated variant rs17824933, linking GWA findings with variation in gene expression. Finally, we show in CD4+ T cells from a further individual that high-throughput sequencing of genomic DNA and RNA-seq following enrichment for targeted gene sequences by sequence capture methods offers an unbiased means to increase the read depth for transcripts of interest, and therefore a method to investigate the regulatory role of many disease-associated genetic variants.

Association of the vitamin D metabolism gene CYP27B1 with type 1 diabetes

OBJECTIVE—Epidemiological studies have linked vitamin D deficiency with the susceptibility to type 1 diabetes. Higher levels of the active metabolite 1α,25-dihydroxyvitamin D (1α,25(OH)2D) could protect from immune destruction of the pancreatic β-cells. 1α,25(OH)2D is derived from its precursor 25-hydroxyvitamin D by the enzyme 1α-hydroxylase encoded by the CYP27B1 gene and is inactivated by 24-hydroxylase encoded by the CYP24A1 gene. Our aim was to study the association between the CYP27B1 and CYP24A1 gene polymorphisms and type 1 diabetes. RESEARCH DESIGN AND METHODS—We studied 7,854 patients with type 1 diabetes, 8,758 control subjects from the U.K., and 2,774 affected families. We studied four CYP27B1 variants, including common polymorphisms −1260C>A (rs10877012) and +2838T>C (rs4646536) and 16 tag polymorphisms in the CYP24A1 gene. RESULTS—We found evidence of association with type 1 diabetes for CYP27B1 −1260 and +2838 polymorphisms, which are in perfect linkage disequilibrium. The common C allele of CYP27B1 −1260 was associated with an increased disease risk in the case-control analysis (odds ratio for the C/C genotype 1.22, P = 9.6 × 10−4) and in the fully independent collection of families (relative risk for the C/C genotype 1.33, P = 3.9 × 10−3). The combined P value for an association with type 1 diabetes was 3.8 × 10−6. For the CYP24A1 gene, we found no evidence of association with type 1 diabetes (multilocus test, P = 0.23). CONCLUSIONS—The present data provide evidence that common inherited variation in the vitamin D metabolism affects susceptibility to type 1 diabetes.

Type 1 Diabetes-Associated IL2RA Variation Lowers IL-2 Signaling and Contributes to Diminished CD4+CD25+ Regulatory T Cell Function

Numerous reports have demonstrated that CD4+CD25+ regulatory T cells (Tregs) from individuals with a range of human autoimmune diseases, including type 1 diabetes, are deficient in their ability to control autologous proinflammatory responses when compared with nondiseased, control individuals. Treg dysfunction could be a primary, causal event or may result from perturbations in the immune system during disease development. Polymorphisms in genes associated with Treg function, such as IL2RA, confer a higher risk of autoimmune disease. Although this suggests a primary role for defective Tregs in autoimmunity, a link between IL2RA gene polymorphisms and Treg function has not been examined. We addressed this by examining the impact of an IL2RA haplotype associated with type 1 diabetes on Treg fitness and suppressive function. Studies were conducted using healthy human subjects to avoid any confounding effects of disease. We demonstrated that the presence of an autoimmune disease-associated IL2RA haplotype correlates with diminished IL-2 responsiveness in Ag-experienced CD4+ T cells, as measured by phosphorylation of STAT5a, and is associated with lower levels of FOXP3 expression by Tregs and a reduction in their ability to suppress proliferation of autologous effector T cells. These data offer a rationale that contributes to the molecular and cellular mechanisms through which polymorphisms in the IL-2RA gene affect immune regulation, and consequently upon susceptibility to autoimmune and inflammatory diseases.

Clinical-Grade Multipotent Adult Progenitor Cells Durably Control Pathogenic T Cell Responses in Human Models of Transplantation and Autoimmunity

A major goal of immunotherapy remains the control of pathogenic T cell responses that drive autoimmunity and allograft rejection. Adherent progenitor cells, including mesenchymal stromal cells (MSCs) and multipotent adult progenitor cells (MAPCs), represent attractive immunomodulatory cell therapy candidates currently active in clinical trials. MAPCs can be distinguished from MSCs on the basis of cellular phenotype, size, transcriptional profile, and expansion capacity. However, despite their ongoing evaluation in autoimmune and allogeneic solid organ transplantation settings, data supporting the immune regulatory potential of clinical-grade MAPCs are limited. In this study, we used allogeneic islet transplantation as a model indication to assess the ability of clinical-grade MAPCs to control T cell responses that drive immunopathology in human autoimmune disease and allograft rejection. MAPCs suppressed T cell proliferation and Th1 and Th17 cytokine production while increasing secretion of IL-10 and were able to suppress effector functions of bona fide autoreactive T cells from individuals with type 1 diabetes mellitus, including killing of human islets. Furthermore, MAPCs favored the proliferation of regulatory T cells during homeostatic expansion driven by γ-chain cytokines and exerted a durable, yet reversible, control of T cell function. MAPC suppression required licensing and proceeded via IDO-mediated tryptophan catabolism. Therefore, the common immune modulatory characteristics of clinical-grade MAPCs shown in this study suggest that they can be regarded as an alternative source of adult progenitor cells with similar clinical usefulness to MSCs. Taken collectively, these findings may guide the successful deployment of both MSCs and MAPCs for the amelioration of human autoimmunity and allograft rejection.

Statistical colocalization of monocyte gene expression and genetic risk variants for type 1 diabetes

One mechanism by which disease-associated DNA variation can alter disease risk is altering gene expression. However, linkage disequilibrium (LD) between variants, mostly single-nucleotide polymorphisms (SNPs), means it is not sufficient to show that a particular variant associates with both disease and expression, as there could be two distinct causal variants in LD. Here, we describe a formal statistical test of colocalization and apply it to type 1 diabetes (T1D)-associated regions identified mostly through genome-wide association studies and expression quantitative trait loci (eQTLs) discovered in a recently determined large monocyte expression data set from the Gutenberg Health Study (1370 individuals), with confirmation sought in an additional data set from the Cardiogenics Transcriptome Study (558 individuals). We excluded 39 out of 60 overlapping eQTLs in 49 T1D regions from possible colocalization and identified 21 coincident eQTLs, representing 21 genes in 14 distinct T1D regions. Our results reflect the importance of monocyte (and their derivatives, macrophage and dendritic cell) gene expression in human T1D and support the candidacy of several genes as causal factors in autoimmune pancreatic beta-cell destruction, including AFF3, CD226, CLECL1, DEXI, FKRP, PRKD2, RNLS, SMARCE1 and SUOX, in addition to the recently described GPR183 (EBI2) gene.

Confirmation of novel type 1 diabetes risk loci in families

Aims/hypothesisOver 50 regions of the genome have been associated with type 1 diabetes risk, mainly using large case/control collections. In a recent genome-wide association (GWA) study, 18 novel susceptibility loci were identified and replicated, including replication evidence from 2,319 families. Here, we, the Type 1 Diabetes Genetics Consortium (T1DGC), aimed to exclude the possibility that any of the 18 loci were false-positives due to population stratification by significantly increasing the statistical power of our family study.MethodsWe genotyped the most disease-predicting single-nucleotide polymorphisms at the 18 susceptibility loci in 3,108 families and used existing genotype data for 2,319 families from the original study, providing 7,013 parent–child trios for analysis. We tested for association using the transmission disequilibrium test.ResultsSeventeen of the 18 susceptibility loci reached nominal levels of significance (p < 0.05) in the expanded family collection, with 14q24.1 just falling short (p = 0.055). When we allowed for multiple testing, ten of the 17 nominally significant loci reached the required level of significance (p < 2.8 × 10−3). All susceptibility loci had consistent direction of effects with the original study.Conclusions/interpretationThe results for the novel GWA study-identified loci are genuine and not due to population stratification. The next step, namely correlation of the most disease-associated genotypes with phenotypes, such as RNA and protein expression analyses for the candidate genes within or near each of the susceptibility regions, can now proceed.

Regulatory T Cell Responses in Participants with Type 1 Diabetes after a Single Dose of Interleukin-2: A Non-Randomised, Open Label, Adaptive Dose-Finding Trial

Background Interleukin-2 (IL-2) has an essential role in the expansion and function of CD4+ regulatory T cells (Tregs). Tregs reduce tissue damage by limiting the immune response following infection and regulate autoreactive CD4+ effector T cells (Teffs) to prevent autoimmune diseases, such as type 1 diabetes (T1D). Genetic susceptibility to T1D causes alterations in the IL-2 pathway, a finding that supports Tregs as a cellular therapeutic target. Aldesleukin (Proleukin; recombinant human IL-2), which is administered at high doses to activate the immune system in cancer immunotherapy, is now being repositioned to treat inflammatory and autoimmune disorders at lower doses by targeting Tregs. Methods and Findings To define the aldesleukin dose response for Tregs and to find doses that increase Tregs physiologically for treatment of T1D, a statistical and systematic approach was taken by analysing the pharmacokinetics and pharmacodynamics of single doses of subcutaneous aldesleukin in the Adaptive Study of IL-2 Dose on Regulatory T Cells in Type 1 Diabetes (DILT1D), a single centre, non-randomised, open label, adaptive dose-finding trial with 40 adult participants with recently diagnosed T1D. The primary endpoint was the maximum percentage increase in Tregs (defined as CD3+CD4+CD25highCD127low) from the baseline frequency in each participant measured over the 7 d following treatment. There was an initial learning phase with five pairs of participants, each pair receiving one of five pre-assigned single doses from 0.04 × 106 to 1.5 × 106 IU/m2, in order to model the dose-response curve. Results from each participant were then incorporated into interim statistical modelling to target the two doses most likely to induce 10% and 20% increases in Treg frequencies. Primary analysis of the evaluable population (n = 39) found that the optimal doses of aldesleukin to induce 10% and 20% increases in Tregs were 0.101 × 106 IU/m2 (standard error [SE] = 0.078, 95% CI = −0.052, 0.254) and 0.497 × 106 IU/m2 (SE = 0.092, 95% CI = 0.316, 0.678), respectively. On analysis of secondary outcomes, using a highly sensitive IL-2 assay, the observed plasma concentrations of the drug at 90 min exceeded the hypothetical Treg-specific therapeutic window determined in vitro (0.015–0.24 IU/ml), even at the lowest doses (0.040 × 106 and 0.045 × 106 IU/m2) administered. A rapid decrease in Treg frequency in the circulation was observed at 90 min and at day 1, which was dose dependent (mean decrease 11.6%, SE = 2.3%, range 10.0%–48.2%, n = 37), rebounding at day 2 and increasing to frequencies above baseline over 7 d. Teffs, natural killer cells, and eosinophils also responded, with their frequencies rapidly and dose-dependently decreased in the blood, then returning to, or exceeding, pretreatment levels. Furthermore, there was a dose-dependent down modulation of one of the two signalling subunits of the IL-2 receptor, the β chain (CD122) (mean decrease = 58.0%, SE = 2.8%, range 9.8%–85.5%, n = 33), on Tregs and a reduction in their sensitivity to aldesleukin at 90 min and day 1 and 2 post-treatment. Due to blood volume requirements as well as ethical and practical considerations, the study was limited to adults and to analysis of peripheral blood only. Conclusions The DILT1D trial results, most notably the early altered trafficking and desensitisation of Tregs induced by a single ultra-low dose of aldesleukin that resolves within 2–3 d, inform the design of the next trial to determine a repeat dosing regimen aimed at establishing a steady-state Treg frequency increase of 20%–50%, with the eventual goal of preventing T1D. Trial Registration ISRCTN Registry ISRCTN27852285; ClinicalTrials.gov NCT01827735