P Dubois

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.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

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.

A meta-analysis of genome-wide association scans identifies IL18RAP, PTPN2, TAGAP, and PUS10 as shared risk loci for Crohn's disease and celiac disease

Crohns disease (CD) and celiac disease (CelD) are chronic intestinal inflammatory diseases, involving genetic and environmental factors in their pathogenesis. The two diseases can co-occur within families, and studies suggest that CelD patients have a higher risk to develop CD than the general population. These observations suggest that CD and CelD may share common genetic risk loci. Two such shared loci, IL18RAP and PTPN2, have already been identified independently in these two diseases. The aim of our study was to explicitly identify shared risk loci for these diseases by combining results from genome-wide association study (GWAS) datasets of CD and CelD. Specifically, GWAS results from CelD (768 cases, 1,422 controls) and CD (3,230 cases, 4,829 controls) were combined in a meta-analysis. Nine independent regions had nominal association p-value <1.0×10−5 in this meta-analysis and showed evidence of association to the individual diseases in the original scans (p-value <1×10−2 in CelD and <1×10−3 in CD). These include the two previously reported shared loci, IL18RAP and PTPN2, with p-values of 3.37×10−8 and 6.39×10−9, respectively, in the meta-analysis. The other seven had not been reported as shared loci and thus were tested in additional CelD (3,149 cases and 4,714 controls) and CD (1,835 cases and 1,669 controls) cohorts. Two of these loci, TAGAP and PUS10, showed significant evidence of replication (Bonferroni corrected p-values <0.0071) in the combined CelD and CD replication cohorts and were firmly established as shared risk loci of genome-wide significance, with overall combined p-values of 1.55×10−10 and 1.38×10−11 respectively. Through a meta-analysis of GWAS data from CD and CelD, we have identified four shared risk loci: PTPN2, IL18RAP, TAGAP, and PUS10. The combined analysis of the two datasets provided the power, lacking in the individual GWAS for single diseases, to detect shared loci with a relatively small effect.

Evolutionary and Functional Analysis of Celiac Risk Loci Reveals SH2B3 as a Protective Factor against Bacterial Infection

Celiac disease (CD) is an intolerance to dietary proteins of wheat, barley, and rye. CD may have substantial morbidity, yet it is quite common with a prevalence of 1%-2% in Western populations. It is not clear why the CD phenotype is so prevalent despite its negative effects on human health, especially because appropriate treatment in the form of a gluten-free diet has only been available since the 1950s, when dietary gluten was discovered to be the triggering factor. The high prevalence of CD might suggest that genes underlying this disease may have been favored by the process of natural selection. We assessed signatures of selection for ten confirmed CD-associated loci in several genome-wide data sets, comprising 8154 controls from four European populations and 195 individuals from a North African population, by studying haplotype lengths via the integrated haplotype score (iHS) method. Consistent signs of positive selection for CD-associated derived alleles were observed in three loci: IL12A, IL18RAP, and SH2B3. For the SH2B3 risk allele, we also show a difference in allele frequency distribution (Fst) between HapMap phase II populations. Functional investigation of the effect of the SH2B3 genotype in response to lipopolysaccharide and muramyl dipeptide revealed that carriers of the SH2B3 rs3184504*A risk allele showed stronger activation of the NOD2 recognition pathway. This suggests that SH2B3 plays a role in protection against bacteria infection, and it provides a possible explanation for the selective sweep on SH2B3, which occurred sometime between 1200 and 1700 years ago.

HLA-DQA1-HLA-DRB1 variants confer susceptibility to pancreatitis induced by thiopurine immunosuppressants

Pancreatitis occurs in approximately 4% of patients treated with the thiopurines azathioprine or mercaptopurine. Its development is unpredictable and almost always leads to drug withdrawal. We identified patients with inflammatory bowel disease (IBD) who had developed pancreatitis within 3 months of starting these drugs from 168 sites around the world. After detailed case adjudication, we performed a genome-wide association study on 172 cases and 2,035 controls with IBD. We identified strong evidence of association within the class II HLA region, with the most significant association identified at rs2647087 (odds ratio 2.59, 95% confidence interval 2.07–3.26, P = 2 × 10−16). We replicated these findings in an independent set of 78 cases and 472 controls with IBD matched for drug exposure. Fine mapping of the HLA region identified association with the HLA-DQA1*02:01–HLA-DRB1*07:01 haplotype. Patients heterozygous at rs2647087 have a 9% risk of developing pancreatitis after administration of a thiopurine, whereas homozygotes have a 17% risk.

Translational Mini-Review Series on the Immunogenetics of Gut Disease: Immunogenetics of coeliac disease

Recent advances in immunological and genetic research in coeliac disease provide new and complementary insights into the immune response driving this chronic intestinal inflammatory disorder. Both approaches confirm the central importance of T cell‐mediated immune responses to disease pathogenesis and have further begun to highlight other relevant components of the mucosal immune system, including innate immunity and the control of lymphocyte trafficking to the mucosa. In the last year, the first genome wide association study in celiac disease led to the identification of multiple new risk variants. These risk regions implicate genes involved in the immune system. Overlap with autoimmune diseases is striking with several of these regions being shown to confer susceptibility to other chronic immune‐mediated diseases, particularly type 1 diabetes.

Rare and functional SIAE variants are not associated with autoimmune disease risk in up to 66,924 individuals of European ancestry

Rare and functional SIAE variants are not associated with autoimmune disease risk in up to 66,924 individuals of European ancestry

Association study of the IL18RAP locus in three European populations with coeliac disease

Coeliac disease is caused by dietary gluten, triggering a chronic inflammation of the small intestine in genetically predisposed individuals. Recently, a risk locus on chromosome 2q11-q12, harbouring interleukin 18 receptor accessory protein (IL18RAP) and three other genes, was suggested for coeliac disease. IL18 has been shown to play an important role in T helper type 1 activity in coeliac disease, making this locus a highly interesting candidate. In this study, two previously indicated risk variants at the IL18RAP locus (rs13015714 and rs917997) were tested for genetic association in 1638 cases with coeliac disease and 1385 control individuals from the Finnish, Hungarian and Italian populations. The protein expression level of IL18RAP was also compared between risk allele carriers and non-carriers by Western blotting. Furthermore, immunohistochemical analysis was performed to study IL18RAP protein expression in small intestinal biopsies of untreated and treated coeliac patients and controls. We confirmed genetic association and dose effects of variants at the 2q12.1 locus with coeliac disease in the Hungarian population. The GA haplotype of the markers rs13015714 and rs917997 showed the strongest association (P = 0.0001, odds ratio = 1.475, 95% confidence interval 1.21-1.80). Two putative isoforms of IL18RAP were detected and the ratios and total levels of these isoforms may contribute to the aetiology of coeliac disease. Our study supports IL18RAP as a novel predisposing gene for coeliac disease and highlights the need for further functional studies on this relatively unknown gene in coeliac disease pathogenesis.

New susceptibility genes for ulcerative colitis

Genome-wide association studies in Crohns disease have convincingly identified multiple risk variants. Two new studies in ulcerative colitis, a related inflammatory bowel disease, now reveal a glimpse of the overlapping spectrum of risk variants, highlighting both the commonalities and the unique features of these disorders.