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.

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.

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.

Brief Report: The Role of Rare Protein-Coding Variants in Anti-Tumor Necrosis Factor Treatment Response in Rheumatoid Arthritis

OBJECTIVE

In many rheumatoid arthritis (RA) patients, disease is controlled with anti-tumor necrosis factor (anti-TNF) biologic therapies. However, in a significant number of patients, the disease fails to respond to anti-TNF therapy. We undertook the present study to examine the hypothesis that rare and low-frequency genetic variants might influence response to anti-TNF treatment.

METHODS

We sequenced the coding region of 750 genes in 1,094 RA patients of European ancestry who were treated with anti-TNF. After quality control, 690 genes were included in the analysis. We applied single-variant association and gene-based association tests to identify variants associated with anti-TNF treatment response. In addition, given the key mechanistic role of TNF, we performed gene set analyses of 27 TNF pathway genes.

RESULTS

We identified 14,420 functional variants, of which 6,934 were predicted as nonsynonymous 2,136 of which were further predicted to be "damaging." Despite the fact that the study was well powered, no single variant or gene showed study-wide significant association with change in the outcome measures disease activity or European League Against Rheumatism response. Intriguingly, we observed 3 genes, of 27 with nominal signals of association (P < 0.05), that were involved in the TNF signaling pathway. However, when we performed a rigorous gene set enrichment analysis based on association P value ranking, we observed no evidence of enrichment of association at genes involved in the TNF pathway (P_enrichment  = 0.15, based on phenotype permutations).

CONCLUSION

Our findings suggest that rare and low-frequency protein-coding variants in TNF signaling pathway genes or other genes do not contribute substantially to anti-TNF treatment response in patients with RA.

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.

TYK2 protein-coding variants protect against rheumatoid arthritis and autoimmunity, with no evidence of major pleiotropic effects on non-autoimmune complex traits

Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10^-21), A928V (rs35018800, OR = 0.53, P = 1.2x10^-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10^-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, P_omnibus = 6x10^-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; P_omnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.

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.

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⁻⁹). 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⁻⁹), 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.

A current challenge in human genetics is to follow-up “hits” from genome-wide association studies (GWAS) to guide drug discovery for complex traits. Previously, we identified a common variant in the CD40 locus as associated with risk of rheumatoid arthritis (RA). Here, we fine-map the CD40 signal of association through a combination of dense genotyping and exonic sequencing in large patient collections. Further, we demonstrate that the RA risk allele is a gain-of-function allele that increases the amount of CD40 on the surface of primary human B lymphocyte cells from healthy control individuals. Based on these observations, we develop a high-throughput assay to recapitulate the biology of the RA risk allele in a system suitable for a small molecule drug screen. After a series of primary screens and counter screens, we identify small molecules that inhibit CD40-mediated NF-kB signaling in human B cells. While this is only the first step towards a more comprehensive effort to identify CD40-specific inhibitors that may be used to treat RA, our study demonstrates a successful strategy to progress from a GWAS to a drug screen for complex traits such as 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.

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.