Systematic genetic analysis of early-onset gout: ABCG2 is the only associated locus.

OBJECTIVE: The aim of this study was to examine whether serum urate-associated genetic variants are associated with early-onset gout. METHODS: Participants with gout in the Genetics of Gout in Aotearoa study with available genotyping were included (n = 1648). Early-onset gout was defined as the first presentation of gout <40 years of age. Single nucleotide polymorphisms (SNPs) for the 10 loci most strongly associated with serum urate were genotyped. Allelic association of the SNPs with early-onset gout was tested using logistic regression in an unadjusted model and in a model adjusted for sex, body mass index, tophus presence, flare frequency, serum creatinine and highest serum urate. The analysis was also done in two replication cohorts: Eurogout (n = 704) and Ardea (n = 755), and data were meta-analysed. RESULTS: In the Genetics of Gout in Aotearoa study, there were 638 (42.4%) participants with early-onset gout. The ABCG2 rs2231142 gout risk T-allele was present more frequently in participants with early-onset gout compared with the later-onset group. For the other SNPs tested, no differences in risk allele number were observed. In the allelic association analysis, the ABCG2 rs2231142 T-allele was associated with early-onset gout in unadjusted and adjusted models. Analysis of the replication cohorts confirmed the association of early-onset gout with the ABCG2 rs2231142 T-allele, but not with other serum urate-associated SNPs. In the meta-analysis, the odds ratio (95% CI) for early-onset gout for the ABCG2 rs2231142 T-allele was 1.60 (1.41, 1.83). CONCLUSION: In contrast to other serum urate-raising variants, the ABCG2 rs2231142 T-allele is strongly associated with early-onset gout.

Mast cells in early rheumatoid arthritis associate with disease severity and support B cell autoantibody production.

OBJECTIVES: Mast cells (MCs) are involved in the pathogenesis of rheumatoid arthritis (RA). However, their contribution remains controversial. To establish their role in RA, we analysed their presence in the synovium of treatment-naïve patients with early RA and their association and functional relationship with histological features of synovitis. METHODS: Synovial tissue was obtained by ultrasound-guided biopsy from treatment-naïve patients with early RA (n=99). Immune cells (CD3/CD20/CD138/CD68) and their relationship with CD117+MCs in synovial tissue were analysed by immunohistochemistry (IHC) and immunofluorescence (IF). The functional involvement of MCs in ectopic lymphoid structures (ELS) was investigated in vitro, by coculturing MCs with naïve B cells and anticitrullinated protein antibodies (ACPA)-producing B cell clones, and in vivo in interleukin-27 receptor alpha (IL27ra)-deficient and control mice during antigen-induced arthritis (AIA). RESULTS: High synovial MC counts are associated with local and systemic inflammation, autoantibody positivity and high disease activity. IHC/IF showed that MCs reside at the outer border of lymphoid aggregates. Furthermore, human MCs promote the activation and differentiation of naïve B cells and induce the production of ACPA, mainly via contact-dependent interactions. In AIA, synovial MC numbers increase in IL27ra deficient mice, in association with ELS and worse disease activity. CONCLUSIONS: Synovial MCs identify early RA patients with a severe clinical form of synovitis characterised by the presence of ELS.

Repeated FcεRI triggering reveals modified mast cell function related to chronic allergic responses in tissue.

BACKGROUND: Activation of mast cells through FcεRI plays an important role in acute allergic reactions. However, little is known about the function of mast cells in patients with chronic allergic inflammation or the effect of repeated FcεRI triggering occurring in such responses. OBJECTIVE: We aimed to identify changes in mast cell function after repeated FcεRI triggering and to correlate these changes to chronic allergic responses in tissue. METHODS: Human cord blood-derived mast cells were treated for 2 weeks with anti-IgE. The function of naive or treated mast cells was analyzed by means of RNA sequencing, quantitative RT-PCR, flow cytometry, and functional assays. Protein secretion was measured with ELISAs and multiplex assays. RESULTS: We observed several changes in mast cell function after repeated anti-IgE triggering. Although the acute response was dampened, we identified 289 genes significantly upregulated after repeated anti-IgE. Most of these genes (84%) were not upregulated after a single anti-IgE stimulus, indicating a significantly different response mode characterized by increased antigen presentation, response to bacteria, and chemotaxis. Changes in mast cell function were related to changes in expression of the transcription factors RXRA and BATF and others. Importantly, we found a substantial overlap between genes upregulated after repeated anti-IgE triggering and genes upregulated in tissue from patients with chronic allergy, in particular those of patients with chronic rhinosinusitis. CONCLUSION: Our study provides evidence for intrinsic modulation of mast cell function on repeated FcεRI-mediated activation. The overlap with gene expression in tissues is suggestive of a direct link between repeated IgE-mediated activation of mast cells and chronic allergy.

Rare, low-frequency, and common variants in the protein-coding sequence of biological candidate genes from GWASs contribute to risk of rheumatoid arthritis.

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p(enrichment) = 6.4 × 10(-4)). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10(-6)). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA.

Human genetics in rheumatoid arthritis guides a high-throughput drug screen of the CD40 signaling pathway.

Although genetic and non-genetic studies in mouse and human implicate the CD40 pathway in rheumatoid arthritis (RA), there are no approved drugs that inhibit CD40 signaling for clinical care in RA or any other disease. Here, we sought to understand the biological consequences of a CD40 risk variant in RA discovered by a previous genome-wide association study (GWAS) and to perform a high-throughput drug screen for modulators of CD40 signaling based on human genetic findings. First, we fine-map the CD40 risk locus in 7,222 seropositive RA patients and 15,870 controls, together with deep sequencing of CD40 coding exons in 500 RA cases and 650 controls, to identify a single SNP that explains the entire signal of association (rs4810485, P = 1.4×10(-9)). Second, we demonstrate that subjects homozygous for the RA risk allele have ∼33% more CD40 on the surface of primary human CD19+ B lymphocytes than subjects homozygous for the non-risk allele (P = 10(-9)), a finding corroborated by expression quantitative trait loci (eQTL) analysis in peripheral blood mononuclear cells from 1,469 healthy control individuals. Third, we use retroviral shRNA infection to perturb the amount of CD40 on the surface of a human B lymphocyte cell line (BL2) and observe a direct correlation between amount of CD40 protein and phosphorylation of RelA (p65), a subunit of the NF-κB transcription factor. Finally, we develop a high-throughput NF-κB luciferase reporter assay in BL2 cells activated with trimerized CD40 ligand (tCD40L) and conduct an HTS of 1,982 chemical compounds and FDA-approved drugs. After a series of counter-screens and testing in primary human CD19+ B cells, we identify 2 novel chemical inhibitors not previously implicated in inflammation or CD40-mediated NF-κB signaling. Our study demonstrates proof-of-concept that human genetics can be used to guide the development of phenotype-based, high-throughput small-molecule screens to identify potential novel therapies in complex traits such as RA.

Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12.

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power.

Genetic risk scores and number of autoantibodies in patients with rheumatoid arthritis.

OBJECTIVE: Certain HLA-DRB1 alleles and single-nucleotide polymorphisms (SNPs) are associated with rheumatoid arthritis (RA). Our objective was to examine the combined effect of these associated variants, calculated as a cumulative genetic risk score (GRS) on RA predisposition, as well as the number of autoantibodies (none, one or two present). METHOD: We calculated four GRSs in 4956 patients and 4983 controls from four European countries. All four scores contained data on 22 non-HLA-risk SNPs, and three scores also contained HLA-DRB1 genotypes but had different HLA typing resolution. Most patients had data on both rheumatoid factor (RF) and anti-citrullinated proteins antibodies (ACPA). The GRSs were standardised (std.GRS) to account for population heterogeneity. Discrimination between patients and controls was examined by receiveroperating characteristics curves, and the four std.GRSs were compared across subgroups according to autoantibody status. RESULTS: The std.GRS improved its discriminatory ability between patients and controls when HLA-DRB1 data of higher resolution were added to the combined score. Patients had higher mean std.GRS than controls (p=7.9×10(-156)), and this score was significantly higher in patients with autoantibodies (shown for both RF and ACPA). Mean std.GRS was also higher in those with two versus one autoantibody (p=3.7×10(-23)) but was similar in patients without autoantibodies and controls (p=0.12). CONCLUSIONS: The GRS was associated with the number of autoantibodies and to both RF and ACPA positivity. ACPA play a more important role than RF with regards to the genetic risk profile, but stratification of patients according to both RF and ACPA may optimise future genetic studies.

Immunogenetics of rheumatoid arthritis: Understanding functional implications.

The last decade has seen a dramatic technological revolution. The characterisation of the majority of the common variations in our genetic code in 2003 precipitated the discovery of the genetic risk factors predisposing to Rheumatoid Arthritis development and progression. Prior to 2007, only a handful of genetic risk factors had been identified, HLA, PTPN22 and CTLA4. Since then, over 100 genetic risk loci have been described, with the prediction that an ever-increasing number of risk alleles with consistently decreasing effect sizes will be discovered in the years to come. Each risk locus harbours multiple candidate genes and the proof of causality of each of these candidates is as yet unknown. An enrichment of these RA-associated genes is found in three pathways: T-cell receptor signalling, JAK-STAT signalling and the NF-κB signalling cascade, and currently drugs targeting these pathways are available for the treatment of RA. However, the role that RA-associated genes have in these pathways and how they contribute to disease is not always clear. Major efforts in understanding the contribution of genetic risk factors are currently under way with studies querying the role of genetic variation in gene expression of coding and non-coding genes, epigenetic marks and other regulatory mechanisms yielding ever more valuable insights into mechanisms of disease. Recent work has suggested a possible enrichment of non-coding RNAs as well as super-enhancers in RA genetic loci indicating possible new insights into disease mechanism. This review brings together these emerging genetic data with an emphasis on the immunogenetic links these findings have provided and what we expect the future will bring.

Genetic basis of autoantibody positive and negative rheumatoid arthritis risk in a multi-ethnic cohort derived from electronic health records.

Discovering and following up on genetic associations with complex phenotypes require large patient cohorts. This is particularly true for patient cohorts of diverse ancestry and clinically relevant subsets of disease. The ability to mine the electronic health records (EHRs) of patients followed as part of routine clinical care provides a potential opportunity to efficiently identify affected cases and unaffected controls for appropriate-sized genetic studies. Here, we demonstrate proof-of-concept that it is possible to use EHR data linked with biospecimens to establish a multi-ethnic case-control cohort for genetic research of a complex disease, rheumatoid arthritis (RA). In 1,515 EHR-derived RA cases and 1,480 controls matched for both genetic ancestry and disease-specific autoantibodies (anti-citrullinated protein antibodies [ACPA]), we demonstrate that the odds ratios and aggregate genetic risk score (GRS) of known RA risk alleles measured in individuals of European ancestry within our EHR cohort are nearly identical to those derived from a genome-wide association study (GWAS) of 5,539 autoantibody-positive RA cases and 20,169 controls. We extend this approach to other ethnic groups and identify a large overlap in the GRS among individuals of European, African, East Asian, and Hispanic ancestry. We also demonstrate that the distribution of a GRS based on 28 non-HLA risk alleles in ACPA+ cases partially overlaps with ACPA- subgroup of RA cases. Our study demonstrates that the genetic basis of rheumatoid arthritis risk is similar among cases of diverse ancestry divided into subsets based on ACPA status and emphasizes the utility of linking EHR clinical data with biospecimens for genetic studies.

Identification of a genetic variant for joint damage progression in autoantibody-positive rheumatoid arthritis.

BACKGROUND: Joint destruction is a hallmark of autoantibody-positive rheumatoid arthritis (RA), though the severity is highly variable between patients. The processes underlying these interindividual differences are incompletely understood. METHODS: We performed a genome-wide association study on the radiological progression rate in 384 autoantibody-positive patients with RA. In stage-II 1557 X-rays of 301 Dutch autoantibody-positive patients with RA were studied and in stage-III 861 X-rays of 742 North American autoantibody-positive patients with RA. Sperm-Associated Antigen 16 (SPAG16) expression in RA synovium and fibroblast-like synoviocytes (FLS) was examined using Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR) and immunohistochemistry. FLS secrete metalloproteinases that degrade cartilage and bone. SPAG16 genotypes were related to matrix metalloproteinase (MMP)-3 and MMP-1 expression by FLS in vitro and MMP-3 production ex vivo. RESULTS: A cluster of single nucleotide polymorphisms (SNPs) at 2q34, located at SPAG16, associated with the radiological progression rate; rs7607479 reached genome-wide significance. A protective role of rs7607479 was replicated in European and North American patients with RA. Per minor allele, patients had a 0.78-fold (95% CI 0.67 to 0.91) progression rate over 7 years. mRNA and protein expression of SPAG16 in RA synovium and FLS was verified. FLS carrying the minor allele secreted less MMP-3 (p=1.60×10(-2)). Furthermore, patients with RA carrying the minor allele had lower serum levels of MMP-3 (p=4.28×10(-2)). In a multivariate analysis on rs7607479 and MMP-3, only MMP-3 associated with progression (p=2.77×10(-4)), suggesting that the association between SPAG16-rs7607479 and joint damage is mediated via an effect on MMP-3 secretion. CONCLUSIONS: Genetic and functional analyses indicate that SPAG16 influences MMP-3 regulation and protects against joint destruction in autoantibody-positive RA. These findings could enhance risk stratification in autoantibody-positive RA.

Associations of autoantibodies, autoimmune risk alleles, and clinical diagnoses from the electronic medical records in rheumatoid arthritis cases and non-rheumatoid arthritis controls.

OBJECTIVE: The significance of non-rheumatoid arthritis (RA) autoantibodies in patients with RA is unclear. The aim of this study was to assess associations of autoantibodies with autoimmune risk alleles and with clinical diagnoses from the electronic medical records (EMRs) among RA cases and non-RA controls. METHODS: Data on 1,290 RA cases and 1,236 non-RA controls of European genetic ancestry were obtained from the EMRs of 2 large academic centers. The levels of anti-citrullinated protein antibodies (ACPAs), antinuclear antibodies (ANAs), anti-tissue transglutaminase antibodies (AGTAs), and anti-thyroid peroxidase (anti-TPO) antibodies were measured. All subjects were genotyped for autoimmune risk alleles, and the association between number of autoimmune risk alleles present and number of types of autoantibodies present was studied. A phenome-wide association study (PheWAS) was conducted to study potential associations between autoantibodies and clinical diagnoses among RA cases and non-RA controls. RESULTS: The mean ages were 60.7 years in RA cases and 64.6 years in non-RA controls. The proportion of female subjects was 79% in each group. The prevalence of ACPAs and ANAs was higher in RA cases compared to controls (each P < 0.0001); there were no differences in the prevalence of anti-TPO antibodies and AGTAs. Carriage of higher numbers of autoimmune risk alleles was associated with increasing numbers of autoantibody types in RA cases (P = 2.1 × 10(-5)) and non-RA controls (P = 5.0 × 10(-3)). From the PheWAS, the presence of ANAs was significantly associated with a diagnosis of Sjögren's/sicca syndrome in RA cases. CONCLUSION: The increased frequency of autoantibodies in RA cases and non-RA controls was associated with the number of autoimmune risk alleles carried by an individual. PheWAS of EMR data, with linkage to laboratory data obtained from blood samples, provide a novel method to test for the clinical significance of biomarkers in disease.

Meta-analysis of genome-wide association studies in celiac disease and rheumatoid arthritis identifies fourteen non-HLA shared loci.

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P<5 × 10(-8) in a combined analysis of all 50,266 samples, including four SNPs that have not been previously confirmed in either disease: rs10892279 near the DDX6 gene (P(combined) =  1.2 × 10(-12)), rs864537 near CD247 (P(combined) =  2.2 × 10(-11)), rs2298428 near UBE2L3 (P(combined) =  2.5 × 10(-10)), and rs11203203 near UBASH3A (P(combined) =  1.1 × 10(-8)). We also confirmed that 4 gene loci previously established in either CD or RA are associated with the other autoimmune disease at combined P<5 × 10(-8) (SH2B3, 8q24, STAT4, and TRAF1-C5). From the 14 shared gene loci, 7 SNPs showed a genome-wide significant effect on expression of one or more transcripts in the linkage disequilibrium (LD) block around the SNP. These associations implicate antigen presentation and T-cell activation as a shared mechanism of disease pathogenesis and underscore the utility of cross-disease meta-analysis for identification of genetic risk factors with pleiotropic effects between two clinically distinct diseases.

Genome-wide Association Study of Methotrexate-Induced Liver Injury in Rheumatoid Arthritis Patients.

Hepatotoxicity is a serious adverse drug reaction related to methotrexate (MTX). However, the cause of drug-induced liver injury (DILI) is still unclear and unpredictable. Genetic risk factors may predispose for MTX-DILI. Therefore, we conducted a nested case-control genome-wide association study to explore genetic risk factors associated with MTX-DILI. Seven international groups contributed blood samples and data of patients with rheumatoid arthritis who used MTX. MTX-DILI was defined as an alanine aminotransferase (ALT) level of at least three times the upper limit of normal (ULN), to increase contrast controls ALT levels did not raise above two times the ULN. Per study site, control subjects and patients with MTX-DILI (ratio 3:1) were matched for age, gender, and duration of MTX use. Patients were genotyped using Illumina GSA MD-24v1-0 and data were imputed using the 1000 Genomes reference panel. Single-nucleotide polymorphisms (SNPs) were analyzed using an additive genetic model, corrected for sex, country, and age. A P-value of ≤ 5 × 10-8 was considered significant, whereas a  P-value of ≤ 5 × 10-6 was considered suggestive. A total of 108 MTX-DILI cases and 311 controls were included for association analysis. None of the SNPs were significantly associated with MTX-DILI. However, we found seven suggestive genetic variants associated with MTX-DILI (P-values 7.43 × 10-8 to 4.86 × 10-6 ). Of those, five SNPs were in the intronic protein-coding regions of FTCDNL1, BCOR, FGF14, RBMS3, and PFDN4/DOK5. Investigation of candidates SPATA9 (rs72783407), PLCG2 (rs60427389), RAVER2 (rs72675408), JAK1 (rs72675451), PTPN2 (rs2476601), MTHFR C677T (rs1801133), and into the HLA region did not show significant findings. No genetic variants associated with MTX-DILI were found, whereas suggestive SNPs need further investigation.

A large-scale rheumatoid arthritis genetic study identifies association at chromosome 9q33.2.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease affecting both joints and extra-articular tissues. Although some genetic risk factors for RA are well-established, most notably HLA-DRB1 and PTPN22, these markers do not fully account for the observed heritability. To identify additional susceptibility loci, we carried out a multi-tiered, case-control association study, genotyping 25,966 putative functional SNPs in 475 white North American RA patients and 475 matched controls. Significant markers were genotyped in two additional, independent, white case-control sample sets (661 cases/1322 controls from North America and 596 cases/705 controls from The Netherlands) identifying a SNP, rs1953126, on chromosome 9q33.2 that was significantly associated with RA (OR(common) = 1.28, trend P(comb) = 1.45E-06). Through a comprehensive fine-scale-mapping SNP-selection procedure, 137 additional SNPs in a 668 kb region from MEGF9 to STOM on 9q33.2 were chosen for follow-up genotyping in a staged-approach. Significant single marker results (P(comb)<0.01) spanned a large 525 kb region from FBXW2 to GSN. However, a variety of analyses identified SNPs in a 70 kb region extending from the third intron of PHF19 across TRAF1 into the TRAF1-C5 intergenic region, but excluding the C5 coding region, as the most interesting (trend P(comb): 1.45E-06 --> 5.41E-09). The observed association patterns for these SNPs had heightened statistical significance and a higher degree of consistency across sample sets. In addition, the allele frequencies for these SNPs displayed reduced variability between control groups when compared to other SNPs. Lastly, in combination with the other two known genetic risk factors, HLA-DRB1 and PTPN22, the variants reported here generate more than a 45-fold RA-risk differential.

Association of the 6q23 region with the rate of joint destruction in rheumatoid arthritis.

BACKGROUND: Two novel genetic polymorphisms on chromosome 6q23 are associated with susceptibility to rheumatoid arthritis (RA). Both polymorphisms (rs6920220 and rs10499194) reside in a region close to the gene encoding tumour necrosis factor alpha-induced protein 3 (TNFAIP3). TNFAIP3 is a negative regulator of NF-kappaB and is involved in inhibiting TNF-receptor-mediated signalling effects. Interestingly, the initial associations were detected in patients with longstanding RA. However, no association was found for rs10499194 in a Swedish cohort with early arthritis. This might be caused by over-representation of patients with severe disease in cohorts with longstanding RA. OBJECTIVE: To analyse the effect of the 6q23 region on the rate of joint destruction. METHODS: Five single nucleotide polymorphisms in 6q23 were genotyped in 324 Dutch patients with early RA. Genotypes were correlated with progression of radiographic joint damage for a follow-up time of 5 years. RESULTS: Two polymorphisms (rs675520 and rs9376293) were associated with severity of radiographic joint damage in patients positive for anti-citrullinated protein/peptide antibodies (ACPA). Importantly, the effects were present after correction for confounding factors such as secular trends in treatment. CONCLUSIONS: These data associate the 6q23 region with the rate of joint destruction in ACPA+ RA.

The TRAF1-C5 region on chromosome 9q33 is associated with multiple autoimmune diseases.

OBJECTIVES: The TRAF1-C5 locus has recently been identified as a genetic risk factor for rheumatoid arthritis (RA). Since genetic risk factors tend to overlap with several autoimmune diseases, a study was undertaken to investigate whether this region is associated with type 1 diabetes (TID), celiac disease (CD), systemic sclerosis (SSc) and systemic lupus erythematosus (SLE). METHODS: The most consistently associated SNP, rs10818488, was genotyped in a total of 735 patients with T1D, 1049 with CD, 367 with SSc, 746 with SLE and 3494 ethnically- and geographically-matched healthy individuals. The replication sample set consisted of 99 patients with T1D, 272 with SLE and 482 healthy individuals from Crete. RESULTS: A significant association was detected between the rs10818488 A allele and T1D (OR 1.14, p=0.027) and SLE (OR 1.16, p=0.016), which was replicated in 99 patients with T1D, 272 with SLE and 482 controls from Crete (OR 1.64, p=0.002; OR 1.43, p=0.002, respectively). Joint analysis of all patients with T1D (N=961) and all patients with SLE (N=1018) compared with 3976 healthy individuals yielded an allelic common OR of 1.19 (p=0.002) and 1.22 (p=2.6 x 10(-4)), respectively. However, combining our dataset with the T1D sample set from the WTCCC resulted in a non-significant association (OR 1.06, p=0.087). In contrast, previously unpublished results from the SLEGEN study showed a significant association of the same allele (OR 1.19, p=0.0038) with an overall effect of 1.22 (p=1.02 x 10(-6)) in a total of 1577 patients with SLE and 4215 healthy individuals. CONCLUSION: A significant association was found for the TRAF1-C5 locus in SLE, implying that this region lies in a pathway relevant to multiple autoimmune diseases.

Long noncoding RNA H19X is a key mediator of TGF-ß-driven fibrosis.

TGF-ß is a master regulator of fibrosis, driving the differentiation of fibroblasts into apoptosis-resistant myofibroblasts and sustaining the production of extracellular matrix (ECM) components. Here, we identified the nuclear long noncoding RNA (lncRNA) H19X as a master regulator of TGF-ß-driven tissue fibrosis. H19X was consistently upregulated in a wide variety of human fibrotic tissues and diseases and was strongly induced by TGF-ß, particularly in fibroblasts and fibroblast-related cells. Functional experiments following H19X silencing revealed that H19X was an obligatory factor for TGF-ß-induced ECM synthesis as well as differentiation and survival of ECM-producing myofibroblasts. We showed that H19X regulates DDIT4L gene expression, specifically interacting with a region upstream of the DDIT4L gene and changing the chromatin accessibility of a DDIT4L enhancer. These events resulted in transcriptional repression of DDIT4L and, in turn, in increased collagen expression and fibrosis. Our results shed light on key effectors of TGF-ß-induced ECM remodeling and fibrosis.

OTUD6B-AS1 Might Be a Novel Regulator of Apoptosis in Systemic Sclerosis.

Antisense long non-coding RNAs (AS lncRNAs) have increasingly been recognized as important regulators of gene expression and they have been found to play key roles in several diseases. However, very little is known about the role of AS lncRNAs in fibrotic diseases such as systemic sclerosis (SSc). Our recent screening experiments by RNA sequencing showed that ovarian tumor domain containing 6B antisense RNA1 (OTUD6B-AS1) and its sense gene OTUD6B were significantly downregulated in SSc skin biopsies. Therefore, we aimed to identify key regulators of OTUD6B-AS1 and to analyze the functional relevance of OTUD6B-AS1 in SSc. OTUD6B-AS1 and OTUD6B expression in SSc and healthy control (HC) dermal fibroblasts (Fb) after stimulation with transforming growth factor-ß (TGFß), Interleukin (IL)-4, IL-13, and platelet-derived growth factor (PDGF) was analyzed by qPCR. To identify the functional role of OTUD6B-AS1, dermal Fb or human pulmonary artery smooth muscle cells (HPASMC) were transfected with a locked nucleic acid antisense oligonucleotide (ASO) targeting OTUD6B-AS1. Proliferation was measured by BrdU and real-time proliferation assay. Apoptosis was measured by Caspase 3/7 assay and Western blot for cleaved caspase 3. While no difference was recorded at the basal level between HC and SSc dermal Fb, the expression of OTUD6B-AS1 and OTUD6B was significantly downregulated in both SSc and HC dermal Fb after PDGF stimulation in a time-dependent manner. Only mild and inconsistent effects were observed with TGFß, IL-4, and IL-13. OTUD6B-AS1 knockdown in Fb and HPASMC did not affect extracellular matrix or pro-fibrotic/proinflammatory cytokine production. However, OTUD6B-AS1 knockdown significantly increased Cyclin D1 expression at the mRNA and protein level. Moreover, silencing of OTUD6B-AS1 significantly reduced proliferation and suppressed apoptosis in both dermal Fb and HPASMC. OTUD6B-AS1 knockdown did not affect OTUD6B expression at the mRNA level and protein level. Our data suggest that OTUD6B-AS1 regulates proliferation and apoptosis via cyclin D1 expression in a sense gene independent manner. This is the first report investigating the function of OTUD6B-AS1. Our data shed light on a novel apoptosis resistance mechanism in Fb and vascular smooth muscle cells that might be relevant for pathogenesis of SSc.

Antisense Long Non-Coding RNAs Are Deregulated in Skin Tissue of Patients with Systemic Sclerosis.

Systemic sclerosis is an autoimmune disease characterized by fibrosis of skin and multiple organs of which the pathogenesis is poorly understood. We studied differentially expressed coding and non-coding genes in relation to systemic sclerosis pathogenesis with a specific focus on antisense non-coding RNAs. Skin biopsy-derived RNAs from 14 early systemic sclerosis patients and six healthy individuals were sequenced with ion-torrent and analyzed using DEseq2. Overall, 4,901 genes with a fold change >1.5 and a false discovery rate <5% were detected in patients versus controls. Upregulated genes clustered in immunologic, cell adhesion, and keratin-related processes. Interestingly, 676 deregulated non-coding genes were detected, 257 of which were classified as antisense genes. Sense genes expressed opposite of these antisense genes were also deregulated in 42% of the observed sense-antisense gene pairs. The majority of the antisense genes had a similar effect sizes in an independent North American dataset with three genes (CTBP1-AS2, OTUD6B-AS1, and AGAP2-AS1) exceeding the study-wide Bonferroni-corrected P-value (PBonf < 0.0023, Pcombined = 1.1 × 10-9, 1.4 × 10-8, 1.7 × 10-6, respectively). In this study, we highlight that together with coding genes, (antisense) long non-coding RNAs are deregulated in skin tissue of systemic sclerosis patients suggesting a novel class of genes involved in pathogenesis of systemic sclerosis.