SNP Variants in Major Histocompatibility Complex Are Associated with Sarcoidosis Susceptibility-A Joint Analysis in Four European Populations.

Sarcoidosis is a multiorgan inflammatory disorder with heritability estimates up to 66%. Previous studies have shown the major histocompatibility complex (MHC) region to be associated with sarcoidosis, suggesting a functional role for antigen-presenting molecules and immune mediators in the disease pathogenesis. To detect variants predisposing to sarcoidosis and to identify genetic differences between patient subgroups, we studied four genes in the MHC Class III region (LTA, TNF, AGER, BTNL2) and HLA-DRA with tag-SNPs and their relation to HLA-DRB1 alleles. We present results from a joint analysis of four study populations (Finnish, Swedish, Dutch, and Czech). Patients with sarcoidosis (n = 805) were further subdivided based on the disease activity and the presence of Löfgren's syndrome. In a joint analysis, seven SNPs were associated with non-Löfgren sarcoidosis (NL; the strongest association with rs3177928, P = 1.79E-07, OR = 1.9) and eight with Löfgren's syndrome [Löfgren syndrome (LS); the strongest association with rs3129843, P = 3.44E-12, OR = 3.4] when compared with healthy controls (n = 870). Five SNPs were associated with sarcoidosis disease course (the strongest association with rs3177928, P = 0.003, OR = 1.9). The high linkage disequilibrium (LD) between SNPs and an HLA-DRB1 challenged the result interpretation. When the SNPs and HLA-DRB1 alleles were analyzed together, independent association was observed for four SNPs in the HLA-DRA/BTNL2 region: rs3135365 (NL; P = 0.015), rs3177928 (NL; P < 0.001), rs6937545 (LS; P = 0.012), and rs5007259 (disease activity; P = 0.002). These SNPs act as expression quantitative trait loci (eQTL) for HLA-DRB1 and/or HLA-DRB5. In conclusion, we found novel SNPs in BTNL2 and HLA-DRA regions associating with sarcoidosis. Our finding further establishes that polymorphisms in the HLA-DRA and BTNL2 have a role in sarcoidosis susceptibility. This multi-population study demonstrates that at least a part of these associations are HLA-DRB1 independent (e.g., not due to LD) and shared across ancestral origins. The variants that were independent of HLA-DRB1 associations acted as eQTL for HLA-DRB1 and/or -DRB5, suggesting a role in regulating gene expression.

T-cell receptor-HLA-DRB1 associations suggest specific antigens in pulmonary sarcoidosis.

In pulmonary sarcoidosis, CD4(+) T-cells expressing T-cell receptor Vα2.3 accumulate in the lungs of HLA-DRB1*03(+) patients. To investigate T-cell receptor-HLA-DRB1*03 interactions underlying recognition of hitherto unknown antigens, we performed detailed analyses of T-cell receptor expression on bronchoalveolar lavage fluid CD4(+) T-cells from sarcoidosis patients.Pulmonary sarcoidosis patients (n=43) underwent bronchoscopy with bronchoalveolar lavage. T-cell receptor α and ß chains of CD4(+) T-cells were analysed by flow cytometry, DNA-sequenced, and three-dimensional molecular models of T-cell receptor-HLA-DRB1*03 complexes generated.Simultaneous expression of Vα2.3 with the Vß22 chain was identified in the lungs of all HLA-DRB1*03(+) patients. Accumulated Vα2.3/Vß22-expressing T-cells were highly clonal, with identical or near-identical Vα2.3 chain sequences and inter-patient similarities in Vß22 chain amino acid distribution. Molecular modelling revealed specific T-cell receptor-HLA-DRB1*03-peptide interactions, with a previously identified, sarcoidosis-associated vimentin peptide, (Vim)429-443 DSLPLVDTHSKRTLL, matching both the HLA peptide-binding cleft and distinct T-cell receptor features perfectly.We demonstrate, for the first time, the accumulation of large clonal populations of specific Vα2.3/Vß22 T-cell receptor-expressing CD4(+) T-cells in the lungs of HLA-DRB1*03(+) sarcoidosis patients. Several distinct contact points between Vα2.3/Vß22 receptors and HLA-DRB1*03 molecules suggest presentation of prototypic vimentin-derived peptides.

High-Density Genetic Mapping Identifies New Susceptibility Variants in Sarcoidosis Phenotypes and Shows Genomic-driven Phenotypic Differences.

RATIONALE: Sarcoidosis is a multisystem disease of unknown cause. Löfgren's syndrome (LS) is a characteristic subgroup of sarcoidosis that is associated with a good prognosis in sarcoidosis. However, little is known about its genetic architecture or its broader phenotype, non-LS sarcoidosis. OBJECTIVES: To address the genetic architecture of sarcoidosis phenotypes, LS and non-LS. METHODS: An association study in a white Swedish cohort of 384 LS, 664 non-LS, and 2,086 control subjects, totaling 3,134 subjects using a fine-mapping genotyping platform was conducted. Replication was performed in four independent cohorts, three of white European descent (Germany, n = 4,975; the Netherlands, n = 613; and Czech Republic, n = 521), and one of black African descent (United States, n = 1,657), totaling 7,766 subjects. MEASUREMENTS AND MAIN RESULTS: A total of 727 LS-associated variants expanding throughout the extended major histocompatibility complex (MHC) region and 68 non-LS-associated variants located in the MHC class II region were identified and confirmed. A shared overlap between LS and non-LS defined by 17 variants located in the MHC class II region was found. Outside the MHC region, two LS-associated loci, in ADCY3 and between CSMD1 and MCPH1, were observed and replicated. CONCLUSIONS: Comprehensive and integrative analyses of genetics, transcription, and pathway modeling on LS and non-LS indicates that these sarcoidosis phenotypes have different genetic susceptibility, genomic distributions, and cellular activities, suggesting distinct molecular mechanisms in pathways related to immune response with a common region.

Identification of Immune-Relevant Factors Conferring Sarcoidosis Genetic Risk.

RATIONALE: Genetic variation plays a significant role in the etiology of sarcoidosis. However, only a small fraction of its heritability has been explained so far. OBJECTIVES: To define further genetic risk loci for sarcoidosis, we used the Immunochip for a candidate gene association study of immune-associated loci. METHODS: Altogether the study population comprised over 19,000 individuals. In a two-stage design, 1,726 German sarcoidosis cases and 5,482 control subjects were genotyped for 128,705 single-nucleotide polymorphisms using the Illumina Immunochip for the screening step. The remaining 3,955 cases, 7,514 control subjects, and 684 parents of affected offspring were used for validation and replication of 44 candidate and two established risk single-nucleotide polymorphisms. MEASUREMENTS AND MAIN RESULTS: Four novel susceptibility loci were identified with genome-wide significance in the European case-control populations, located on chromosomes 12q24.12 (rs653178; ATXN2/SH2B3), 5q33.3 (rs4921492; IL12B), 4q24 (rs223498; MANBA/NFKB1), and 2q33.2 (rs6748088; FAM117B). We further defined three independent association signals in the HLA region with genome-wide significance, peaking in the BTNL2 promoter region (rs5007259), at HLA-B (rs4143332/HLA-B*0801) and at HLA-DPB1 (rs9277542), and found another novel independent signal near IL23R (rs12069782) on chromosome 1p31.3. CONCLUSIONS: Functional predictions and protein network analyses suggest a prominent role of the drug-targetable IL23/Th17 signaling pathway in the genetic etiology of sarcoidosis. Our findings reveal a substantial genetic overlap of sarcoidosis with diverse immune-mediated inflammatory disorders, which could be of relevance for the clinical application of modern therapeutics.

Gene-gene interaction and RNA splicing profiles of MAP2K4 gene in rheumatoid arthritis.

We performed gene-gene interaction analysis, with HLA-DRB1 shared epitope (SE) alleles for 195 SNPs within immunologically important MAP2K, MAP3K and MAP4K gene families, in 2010 rheumatoid arthritis (RA) patients and 2280 healthy controls. We found a significant statistical interaction for rs10468473 with SE alleles in autoantibody-positive RA. Individuals heterozygous for rs10468473 demonstrated higher expression of total MAP2K4 mRNA in blood, compared to A-allele homozygous. We discovered a novel, putatively translated, "cassette exon" RNA splice form of MAP2K4, differentially expressed in peripheral blood mononuclear cells from 88 RA cases and controls. Within the group of RA patients, we observed a correlation of MAP2K4 isoform expression with carried SE alleles, autoantibody, and rheumatoid factor profiles. TNF-dependent modulation of isoform expression pattern was detected in the Jurkat cell line. Our data suggest a genetic interaction between MAP2K4 and HLA-DRB1, and the importance of rs10468473 and MAP2K4 splice variants in the development of autoantibody-positive RA.

Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis.

Epigenetic mechanisms integrate genetic and environmental causes of disease, but comprehensive genome-wide analyses of epigenetic modifications have not yet demonstrated robust association with common diseases. Using Illumina HumanMethylation450 arrays on 354 anti-citrullinated protein antibody-associated rheumatoid arthritis cases and 337 controls, we identified two clusters within the major histocompatibility complex (MHC) region whose differential methylation potentially mediates genetic risk for rheumatoid arthritis. To reduce confounding factors that have hampered previous epigenome-wide studies, we corrected for cellular heterogeneity by estimating and adjusting for cell-type proportions in our blood-derived DNA samples and used mediation analysis to filter out associations likely to be a consequence of disease. Four CpGs also showed an association between genotype and variance of methylation. The associations for both clusters replicated at least one CpG (P < 0.01), with the rest showing suggestive association, in monocyte cell fractions in an independent cohort of 12 cases and 12 controls. Thus, DNA methylation is a potential mediator of genetic risk.

Genome-wide association analysis reveals 12q13.3-q14.1 as new risk locus for sarcoidosis.

Sarcoidosis is a systemic inflammatory disease of unknown aetiology, influenced by genetic and environmental factors. However, the loci so far identified for sarcoidosis explain only a part of its assumed heritability. To identify further susceptibility loci, we performed a genome-wide association analysis using the Affymetrix 6.0 Human GeneChip followed by validation and replication stages. After quality control, 637 cases, 1233 controls and 677 619 single-nucleotide polymorphisms (SNPs) were available for an initial screening. 99 SNPs were selected for validation in an independent study panel (1664 patients, 2932 controls). SNP rs1050045 was significantly associated with sarcoidosis (corrected p=0.0215) in the validation panel and yielded a p-value of 9.22 × 10(-8) (OR 1.24) in the meta-analysis of the screening and validation stage. A meta-analysis of three populations from Germany, the Czech Republic and Sweden confirmed this finding (p = 0.024; OR 1.14). Fine-mapping and mRNA expression studies pointed to osteosarcoma amplified 9 (OS9) as the most likely candidate for the underlying risk factor. The OS9 protein plays an important role in endoplasmic reticulum-associated protein degradation and acts during Toll-like receptor induced activation of myeloid cells. Expression analyses of OS9 mRNA provide evidence for a functional mechanism underlying the detected association signal.

A novel sarcoidosis risk locus for Europeans on chromosome 11q13.1.

RATIONALE: Sarcoidosis is a complex inflammatory disease with a heterogeneous clinical picture. Among others, an acute and chronic clinical course can be distinguished, for which specific genetic risk factors are known. OBJECTIVES: To identify additional risk loci for sarcoidosis and its acute and chronic subforms, we analyzed imputed data from a genome-wide association scan for these phenotypes. METHODS: After quality control, the genome-wide association scan comprised nearly 1.3 million imputed single-nucleotide polymorphisms based on an Affymetrix 6.0 Gene Chip dataset of 564 German sarcoidosis cases, including 176 acute and 354 chronic cases and 1,575 control subjects. MEASUREMENTS AND MAIN RESULTS: We identified chromosome 11q13.1 (rs479777) as a novel locus influencing susceptibility to sarcoidosis with genome-wide significance. The marker was significantly associated in three distinct German case-control populations and in an additional German family sample with odds ratios ranging from 0.67 to 0.77. This finding was further replicated in two independent European case-control populations from the Czech Republic (odds ratio, 0.75) and from Sweden (odds ratio, 0.79). In a meta-analysis of the included European case-control samples the marker yielded a P value of 2.68 × 10(-18). The locus was previously reported to be associated with Crohn disease, psoriasis, alopecia areata, and leprosy. For sarcoidosis, fine-mapping and expression analysis suggest KCNK4, PRDX5, PCLB3, and most promising CCDC88B as candidates for the underlying risk gene in the associated region. CONCLUSIONS: This study provides striking evidence for association of chromosome 11q13.1 with sarcoidosis in Europeans, and thus identified a further genetic risk locus shared by sarcoidosis, Crohn disease and psoriasis.

Interaction analysis between HLA-DRB1 shared epitope alleles and MHC class II transactivator CIITA gene with regard to risk of rheumatoid arthritis.

HLA-DRB1 shared epitope (SE) alleles are the strongest genetic determinants for autoantibody positive rheumatoid arthritis (RA). One of the key regulators in expression of HLA class II receptors is MHC class II transactivator (CIITA). A variant of the CIITA gene has been found to associate with inflammatory diseases.We wanted to explore whether the risk variant rs3087456 in the CIITA gene interacts with the HLA-DRB1 SE alleles regarding the risk of developing RA. We tested this hypothesis in a case-control study with 11767 individuals from four European Caucasian populations (6649 RA cases and 5118 controls).We found no significant additive interaction for risk alleles among Swedish Caucasians with RA (n = 3869, attributable proportion due to interaction (AP) = 0.2, 95%CI: -0.2-0.5) or when stratifying for anti-citrullinated protein antibodies (ACPA) presence (ACPA positive disease: n = 2945, AP = 0.3, 95%CI: -0.05-0.6, ACPA negative: n = 2268, AP = -0.2, 95%CI: -1.0-0.6). We further found no significant interaction between the main subgroups of SE alleles (DRB1*01, DRB1*04 or DRB1*10) and CIITA. Similar analysis of three independent RA cohorts from British, Dutch and Norwegian populations also indicated an absence of significant interaction between genetic variants in CIITA and SE alleles with regard to RA risk.Our data suggest that risk from the CIITA locus is independent of the major risk for RA from HLA-DRB1 SE alleles, given that no significant interaction between rs3087456 and SE alleles was observed. Since a biological link between products of these genes is evident, the genetic contribution from CIITA and class II antigens in the autoimmune process may involve additional unidentified factors.

The balance of expression of PTPN22 splice forms is significantly different in rheumatoid arthritis patients compared with controls.

BACKGROUND: The R620W variant in protein tyrosine phosphatase non-receptor 22 (PTPN22) is associated with rheumatoid arthritis (RA). The PTPN22 gene has alternatively spliced transcripts and at least two of the splice forms have been confirmed to encode different PTPN22 (LYP) proteins, but detailed information regarding expression of these is lacking, especially with regard to autoimmune diseases. METHODS: We have investigated the mRNA expression of known PTPN22 splice forms with TaqMan real-time PCR in relation to ZNF592 as an endogenous reference in peripheral blood cells from three independent cohorts with RA patients (n = 139) and controls (n = 111) of Caucasian origin. Polymorphisms in the PTPN22 locus (25 SNPs) and phenotypic data (gender, disease activity, ACPA and RF status) were used for analysis. Additionally, we addressed possible effects of methotrexate treatment on PTPN22 expression. RESULTS: We found consistent differences in the expression of the PTPN22 splice forms in unstimulated peripheral blood mononuclear cells between RA patients and normal controls. This difference was more pronounced when comparing the ratio of splice forms and was not affected by methotrexate treatment. CONCLUSIONS: Our data show that RA patients and healthy controls have a shift in balance of expression of splice forms derived from the PTPN22 gene. This balance seems not to be caused by treatment and may be of importance during immune response due to great structural differences in the encoded PTPN22 proteins.