Autoimmune risk variants in ERAP2 are associated with gene-expression levels in thymus.

Genetic polymorphisms in the endoplasmic reticulum aminopeptidase (ERAP)1 and ERAP2 genes have been associated with several autoimmune diseases (AIDs) at a genome-wide significance level. In this study, we performed a cis expression quantitative trait locus (eQTL) screen to investigate whether seven fine-mapped AID single-nucleotide polymorphisms (SNPs) in the ERAP-region influence the gene-expression levels of ERAP1 and ERAP2 in thymus. After quality control, we identified six significant eQTLs. We further assessed the peak eQTL signals, and both genes showed highly significant and independent thymic eQTL signals (P=2.16 × 10-15 and P=8.22 × 10-23, respectively). Interestingly, the peak eQTL signal overlapped with the AID risk loci in ERAP2 (r2>0.94), but were distinct in ERAP1 (r2<0.4). Finally, among the SNPs showing the most significant eQTL associations with ERAP2 (P<3.4 × 10-20), six were located within transcription factor motifs in an enhancer region in thymus. Our study therefore reveals the fine-mapped AID risk variants that act as eQTLs with ERAP2 in thymus, and highlights the potential causal regulatory variants.

Coeliac disease-associated polymorphisms influence thymic gene expression.

Significant associations between coeliac disease (CD) and single nucleotide polymorphisms (SNPs) distributed over 40 genetic regions have been established. The majority of these SNPs are non-coding and 20 SNPs were, by expression quantitative trait loci (eQTL) analysis, found to harbour cis regulatory potential in peripheral blood mononuclear cells (PBMC). Almost all regions contain genes with an immunological relevant function, of which many act in the same biological pathways. One such pathway is T-cell development in the thymus, a pathway previously not explored in CD pathogenesis. The aim of our study was to explore the regulatory potential of the CD-associated SNPs (n=50) by eQTL analysis in thymic tissue from 42 subjects. In total, 43 nominal significant (P<0.05) eQTLs were found within 24 CD-associated chromosomal regions, corresponding to 27 expression-altering SNPs (eSNPs) and 40 probes (eProbes) that represents 39 unique genes (eGenes). Nine significant probe-SNP pairs (corresponding to 8 eSNPs and 7 eGenes) overlapped with previous findings in PBMC (rs12727642-PARK7, rs296547-DDX59, rs917997-IL18RAP, rs842647-AHSA2, rs13003464-AHSA2, rs6974491-ELMO1, rs2074404-NSF (two independent probes) and rs2298428-UBE2L3). When compared across more tissues, we found that 14 eQTLs could represent potentially novel thymus-specific eQTLs. This implies that CD risk polymorphisms could affect gene regulation in thymus.

Four novel coeliac disease regions replicated in an association study of a Swedish-Norwegian family cohort.

Recent genome-wide association studies have identified 1q31 (RGS1), 2q11-12 (IL18RAP), 3p21 (CCR1/CCR3/CCR2), 3q25-26 (IL12A/SCHIP1), 3q28 (LPP), 4q27 (IL2/IL21), 6q25 (TAGAP) and 12q24 (SH2B3) as susceptibility regions for coeliac disease (CD). We have earlier replicated association with the IL2/IL21 region. This study aimed at replicating the remaining regions in a family cohort using the transmission disequilibrium test, which is not prone to population stratification as a source of false-positive results. Nine single nucleotide polymorphisms (SNPs) within these regions were genotyped in 325 Swedish-Norwegian CD families. We found significant associations with the same alleles in the regions 1q31 (rs2816316; P(nc)=0.0060), 3p21 (rs6441961; P(nc)=0.0006), 3q25-26 (rs17810564; P(nc)=0.0316 and rs9811792; P(nc)=0.0434) and 3q28 (rs1464510; P(nc)=0.0037). Borderline, but non-significant, associations were found for rs917997 (IL18RAP), whereas no evidence for association could be obtained for rs13015714 (IL18RAP) or rs1738074 (TAGAP). The lack of replication of the latter SNPs could be because of limited power. rs3184504 (SH2B3) was not analysed because of assay failure. The most significantly associated region, 3p21 (CCR1/CCR3/CCR2), was further analysed by typing of 30 SNPs, with the aim of identifying the causal variant responsible for the initial association. Several SNPs showed association with CD, but none displayed associations stronger than rs6441961, nor did any of them add to the effect initially marked by rs6441961 in a conditional analysis. However, differential effects of rs6441961(*)C carrying haplotypes were indicated, and we thus cannot exclude the possibility that our inability to obtain evidence for multiple independent effects in the CCR1/CCR3/CCR2 gene region was related to a power issue.

Association study of IL2/IL21 and FcgRIIa: significant association with the IL2/IL21 region in Scandinavian coeliac disease families.

The first genome-wide association study performed in a UK coeliac disease (CD) case-control cohort revealed association with a linkage disequilibrium block containing the KIAA1109/Tenr/IL2/IL21 genes. Also recently, an association with a non-synonymous polymorphism in FcgammaRIIa (CD32a) was reported in CD with an unusually strong P-value. We aimed to replicate the reported associations with the single nucleotide polymorphisms rs13119723 A>G and rs6822844 G>T in the KIAA1109/Tenr/IL2/IL21 region and rs1801274 G>A in the FcgammaRIIa gene in a family sample consisting of 325 Swedish/Norwegian families using the robust transmission disequilibrium test. The family sample used in this study included 100 families with two or more children affected by CD and 225 families with one affected child. We could confirm significant association between the polymorphisms rs13119723 A>G and rs6822844 G>T located in the KIAA1109/Tenr/IL2/IL21 region and CD (P-value 0.001 and 0.002, respectively). However, we found no association with the FcgammaRIIa rs1801274 G>A polymorphism (P-value=0.3). In conclusion, our results support the KIAA1109/Tenr/IL2/IL21 region as a true CD susceptibility region.

HLA class II alleles in Norwegian patients with coexisting type 1 diabetes and celiac disease.

BACKGROUND: Type 1 diabetes (T1D) and celiac disease (CeD) are 2 distinct diseases, but there is an increased risk of developing CeD for T1D patients. Both diseases are associated with HLA-class II alleles, such as DQB1 *02:01 and DQB1 *03:02; however, their risk contribution vary between the diseases. MATERIALS AND METHODS: We genotyped HLA-DRB1 and - DQB1 in 215 patients with coexisting T1D and CeD identified from a T1D cohort, and compared them to patients with T1D (N = 487) and CeD (N = 327), as well as healthy controls (N = 368). RESULTS: The patients with coexisting T1D and CeD had an intermediate carrier frequency (72.8%) of the DRB1 *03:01- DQB1 *02:01- DQA1 *05:01 haplotype compared to T1D (64.1%) and CeD (88.7%) patients. The DRB1 *03:01- DQB1 *02:01- DQA1 *05:01/ DRB1 *04- DQB1 *03:02- DQA1 *03 haplotype combination, encoding DQ2.5 and DQ8 molecules, was equally frequent among patients with both T1D and CeD (52.6%) and T1D patients (46.8%) but significantly lower in CeD patients (9.5%). CONCLUSION: Overall, the patients with coexisting T1D and CeD had an HLA profile more similar to T1D patients than CeD patients.

Fine mapping study in Scandinavian families suggests association between coeliac disease and haplotypes in chromosome region 5q32.

The previous genome-wide scan in Scandinavian families supported earlier evidence for linkage of a region on chromosome 5 (5q31-33) to coeliac disease. This study deals with further genetic mapping of an 18 cM region, spanning from marker GAh18A (131.87 Mb) to D5S640 (149.96 Mb). Linkage and association analyses were performed in a two-step approach. First, seven microsatellites were added. Strong evidence for linkage was obtained with a Zlr score of 3.96, P(nc) = 4 x 10(-5) at marker D5S436. The strongest association was with a haplotype consisting of the markers D5S2033 and D5S2490 (P(nc) < 0.001). In the second step, we added 17 microsatellites and 69 single nucleotide polymorphisms (SNPs) to the analysis. These markers were located close to or within candidate genes across the region of approximately 7 Mb beneath the linkage peak marked by D5S2017 and D5S812. A substantial increase of the linkage signal with a maximum Zlr score of 4.6 at marker rs1972644 (P(nc) = 2 x 10(-6)) was obtained and several SNPs showed association. Seven SNPs that individually showed the strongest association were genotyped in a second independent family sample set (225 trios). In the trio family sample as well as in the multiplex family sample, the strongest association was found with SNPs within the region flanked by the associated microsatellites D5S2033 and D5S2490 at 5q32.

Association analysis of MYO9B gene polymorphisms and inflammatory bowel disease in a Norwegian cohort.

Association between single nucleotide polymorphisms (SNPs) within the MYO9B gene and celiac disease was recently reported. The role of MYO9B in celiac disease was suggested to relate to an epithelial barrier defect. The region to which MYO9B localize is also linked with inflammatory bowel disease (IBD). For these reasons, we hypothesize that MYO9B could also be a susceptibility gene in IBD. To address this, we performed an association study of a Norwegian IBD cohort (149 patients with Crohn's disease, 308 patients with ulcerative colitis and 562 healthy controls) using SNPs, which tagged the celiac disease associated MYO9B haplotype. No association between these SNPs and IBD was observed. Our results failed to support the notion that MYO9B is a susceptibility gene in IBD.

Association analysis of the 1858C>T polymorphism in the PTPN22 gene in juvenile idiopathic arthritis and other autoimmune diseases.

A functional single nucleotide polymorphism, 1858C>T, in the PTPN22 gene, encoding a tyrosine phosphatase, has been reported to be associated with type I diabetes and some other autoimmune diseases. To further investigate whether this polymorphism may be a general susceptibility factor for autoimmunity, we performed an association study in five different autoimmune diseases, three previously not tested. We found an association with juvenile idiopathic arthritis (OR=1.41; P=0.04), not previously reported, and a tendency for an association with coeliac disease (OR=1.35; P=0.08). In primary sclerosing cholangitis, no association was observed (OR=0.95; P=0.8). Furthermore, we confirmed the increased risk in rheumatoid arthritis (OR=1.58; P=0.001), but could not find support for an association with systemic lupus erythematosus (OR=0.94; P=0.8). Altogether, we have provided further evidence of an association between autoimmune diseases and the 1858C>T polymorphism in PTPN22.

Genetic analysis of the CD28/CTLA4/ICOS (CELIAC3) region in coeliac disease.

In order to extend our previous findings of genetic linkage to the CD28/CTLA4/ICOS region on chromosome 2q33 (CELIAC3) in coeliac disease (CD), we have investigated 22 genetic markers in 325 Norwegian/Swedish multiplex and simplex CD families. We found both linkage and association with several markers, primarily in the multiplex material. We observed strong linkage disequilibrium (LD) between SNPs (Single Nucleotide Polymorphisms) within an LD block delimited by MH30 and D2S72. A haplotype of this region marked by the alleles -1147*T: + 49*A:CT60*G:CT61*A was significantly associated with CD, suggesting that one or more polymorphisms of this haplotype, possibly -1147*T, are involved in CD susceptibility. The CT60 SNP, a polymorphism found to be most strongly associated with some other immune-mediated diseases, was not associated with CD, as this SNP was part of both associated and non-associated haplotypes. Moreover, our results suggest that CELIAC3 harbours several independent loci contributing to CD susceptibility.