TRAF1/C5 but not PTPRC variants are potential predictors of rheumatoid arthritis response to anti-tumor necrosis factor therapy

Identifying genetic variants associated with chromatin looping and genome function

Here we present a comprehensive HiChIP dataset on naïve CD4 T cells (nCD4) from 30 donors and identify QTLs that associate with genotype-dependent and/or allele-specific variation of HiChIP contacts defining loops between active regulatory regions (iQTLs). We observe a substantial overlap between iQTLs and previously defined eQTLs and histone QTLs, and an enrichment for fine-mapped QTLs and GWAS variants. Furthermore, we describe a distinct subset of nCD4 iQTLs, for which the significant variation of chromatin contacts in nCD4 are translated into significant eQTL trends in CD4 T cell memory subsets. Finally, we define connectivity-QTLs as iQTLs that are significantly associated with concordant genotype-dependent changes in chromatin contacts over a broad genomic region (e.g., GWAS SNP in the RNASET2 locus). Our results demonstrate the importance of chromatin contacts as a complementary modality for QTL mapping and their power in identifying previously uncharacterized QTLs linked to cell-specific gene expression and connectivity.

Epigenetic targets of Janus kinase inhibitors are linked to genetic risks of rheumatoid arthritis

BACKGROUND

Current strategies that target cytokines (e.g., tumor necrosis factor (TNF)-α), or signaling molecules (e.g., Janus kinase (JAK)) have advanced the management for allergies and autoimmune diseases. Nevertheless, the molecular mechanism that underpins its clinical efficacy have largely remained elusive, especially in the local tissue environment. Here, we aimed to identify the genetic, epigenetic, and immunological targets of JAK inhibitors (JAKis), focusing on their effects on synovial fibroblasts (SFs), the major local effectors associated with destructive joint inflammation in rheumatoid arthritis (RA).

METHODS

SFs were activated by cytokines related to inflammation in RA, and were treated with three types of JAKis or a TNF-α inhibitor (TNFi). Dynamic changes in transcriptome and chromatin accessibility were profiled across samples to identify drug targets. Furthermore, the putative targets were validated using luciferase assays and clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing.

RESULTS

We found that both JAKis and the TNFi targeted the inflammatory module including IL6. Conversely, specific gene signatures that were preferentially inhibited by either of the drug classes were identified. Strikingly, RA risk enhancers for CD40 and TRAF1 were distinctively regulated by JAKis and the TNFi. We performed luciferase assays and CRISPR-based genome editing, and successfully fine-mapped the single causal variants in these loci, rs6074022-CD40 and rs7021049-TRAF1.

CONCLUSIONS

JAKis and the TNFi had a direct impact on different RA risk enhancers, and we identified nucleotide-resolution targets for both drugs. Distinctive targets of clinically effective drugs could be useful for tailoring the application of these drugs and future design of more efficient treatment strategies.

Gene Ontology Analysis Highlights Biological Processes Influencing Non-Response to Anti-TNF Therapy in Rheumatoid Arthritis

Anti-TNF therapy has significantly improved disease control in rheumatoid arthritis, but a fraction of rheumatoid arthritis patients do not respond to anti-TNF therapy or lose response over time. Moreover, the mechanisms underlying non-response to anti-TNF therapy remain largely unknown. To date, many single biomarkers of response to anti-TNF therapy have been published but they have not yet been analyzed as a system of interacting nodes. The aim of our study is to systematically elucidate the biological processes underlying non-response to anti-TNF therapy in rheumatoid arthritis using the gene ontologies of previously published predictive biomarkers. Gene networks were constructed based on published biomarkers and then enriched gene ontology terms were elucidated in subgroups using gene ontology software tools. Our results highlight the novel role of proteasome-mediated protein catabolic processes (p = 2.91 × 10⁻¹⁵) and plasma lipoproteins (p = 4.55 × 10⁻¹¹) in anti-TNF therapy response. The results of our gene ontology analysis help elucidate the biological processes underlying non-response to anti-TNF therapy in rheumatoid arthritis and encourage further study of the highlighted processes.

Validation of GWAS-Identified Variants for Anti-TNF Drug Response in Rheumatoid Arthritis: A Meta-Analysis of Two Large Cohorts

We aimed to validate the association of 28 GWAS-identified genetic variants for response to TNF inhibitors (TNFi) in a discovery cohort of 1361 rheumatoid arthritis (RA) patients monitored in routine care and ascertained through the REPAIR consortium and DANBIO registry. We genotyped selected markers and evaluated their association with response to TNFi after 6 months of treatment according to the change in disease activity score 28 (ΔDAS28). Next, we confirmed the most interesting results through meta-analysis of our data with those from the DREAM cohort that included 706 RA patients treated with TNFi. The meta-analysis of the discovery cohort and DREAM registry including 2067 RA patients revealed an overall association of the LINC02549_rs7767069 SNP with a lower improvement in DAS28 that remained significant after correction for multiple testing (per-allele OR_Meta=0.83, P_Meta=0.000077; P_Het=0.61). In addition, we found that each copy of the LRRC55_rs717117G allele was significantly associated with lower improvement in DAS28 in rheumatoid factor (RF)-positive patients (per-allele OR_Meta=0.67, P=0.00058; P_Het=0.06) whereas an opposite but not significant effect was detected in RF-negative subjects (per-allele OR_Meta=1.38, P=0.10; P_Het=0.45; P_Interaction=0.00028). Interestingly, although the identified associations did not survive multiple testing correction, the meta-analysis also showed overall and RF-specific associations for the MAFB_rs6071980 and CNTN5_rs1813443 SNPs with decreased changes in DAS28 (per-allele OR_{Meta_rs6071980} = 0.85, P=0.0059; P_Het=0.63 and OR_{Meta_rs1813443_RF+}=0.81, P=0.0059; P_Het=0.69 and OR_{Meta_rs1813443_RF-}=1.00, P=0.99; P_Het=0.12; P_Interaction=0.032). Mechanistically, we found that subjects carrying the LINC02549_rs7767069T allele had significantly increased numbers of CD45RO+CD45RA+ T cells (P=0.000025) whereas carriers of the LINC02549_rs7767069T/T genotype showed significantly increased levels of soluble scavengers CD5 and CD6 in serum (P=0.00037 and P=0.00041). In addition, carriers of the LRRC55_rs717117G allele showed decreased production of IL6 after stimulation of PBMCs with B burgdorferi and E coli bacteria (P=0.00046 and P=0.00044), which suggested a reduced IL6-mediated anti-inflammatory effect of this marker to worsen the response to TNFi. In conclusion, this study confirmed the influence of the LINC02549 and LRRC55 loci to determine the response to TNFi in RA patients and suggested a weak effect of the MAFB and CNTN5 loci that need to be further investigated.

Development and application of a PCR-HRM molecular diagnostic method of SNPs linked with TNF inhibitor efficacy

Background Clinical evidence indicates that genetic variations may interfere with the mechanism of drug action. Recently, it has been reported that the single nucleotide polymorphisms (SNPs) of STAT4, PTPN2, PSORS1C1 and TRAF3IP2RA genes are associated with the clinical efficacy of tumor necrosis factor (TNF) inhibitors in the treatment of rheumatoid arthritis (RA) patients. Therefore, the detection of the SNPs linked with TNF inhibitor efficacy may provide an important basis for the treatment of RA. This study intended to establish molecular diagnostic methods for genotyping the linked SNPs based on high resolution melting (HRM) curve analysis. Methods The polymerase chain reaction-HRM (PCR-HRM) curve analysis detecting systems were established by designing the primers of the four SNPs, rs7574865G>T, rs7234029A>G, rs2233945C>A and rs33980500C>T, and the performance and clinical applicability of which were evaluated by using the Sanger sequencing method and genotyping test for 208 clinical samples. Results The self-developed molecular diagnostic methods of PCR-HRM were confirmed to be able to correctly genotype the four SNPs, the sensitivity and specificity of which were 100% in this study. The repeatability and reproducibility tests showed that there is little variable in intra-assay and inter-assay (the coefficient of variation ranged from 0.01% to 0.07%). The slight changes of DNA template and primer concentrations, PCR cycle number and reaction system volume had no significant effect on the genotyping performance of the method. The PCR-HRM assays were also applied to other PCR thermocyclers with HRM function and use different saturation fluorescent dyes. Conclusions The PCR-HRM genotyping method established in this study can be applied to the routine molecular diagnosis of rs7574865, rs7234029, rs2233945 and rs33980500.

Association of TRAF1/C5 Locus Polymorphisms with Epilepsy and Clinical Traits in Mexican Patients with Neurocysticercosis

Neurocysticercosis is caused by the establishment of Taenia solium cysts in the central nervous system. Murine cysticercosis by Taenia crassiceps is a useful model of cysticercosis in which the complement component 5 (C5) has been linked to infection resistance/permissiveness. This work aimed to study the possible relevance for human neurocysticercosis of single nucleotide polymorphisms (SNPs) in the C5-TRAF1 region (rs17611 C/T, rs992670 G/A, rs25681 G/A, rs10818488 A/G, and rs3761847 G/A) in a Mexican population and associated with clinical and radiological traits related to neurocysticercosis severity (cell count in the cerebrospinal fluid [CSF cellularity], parasite location and parasite load in the brain, parasite degenerating stage, and epilepsy). The AG genotype of the rs3761847 SNP showed a tendency to associate with multiple brain parasites, while the CT and GG genotypes of the rs17611 and rs3761847 SNPs, respectively, showed a tendency to associate with low CSF cellularity. The rs3761847 SNP was associated with epilepsy under a dominant model, whereas rs10818488 was associated with CSF cellularity and parasite load under dominant and recessive models, respectively. For haplotypes, C5- and the TRAF1-associated SNPs were, respectively, in strong linkage disequilibrium with each other; thus, these haplotypes were studied independently. For C5 SNPs, carrying the CAA haplotype increases the risk of showing high CSF cellularity 3-fold and the risk of having extraparenchymal parasites 4-fold, two conditions that are related to severe disease. For TRAF1 SNPs, the GA and AG haplotypes were associated with CSF cellularity, and the AG haplotype was associated with epilepsy. Overall, these findings support the clear participation of C5 and TRAF1 in the risk of developing severe neurocysticercosis in the Mexican population.

TRAF1 Signaling in Human Health and Disease

Tumor necrosis factor receptor (TNFR) associated factor 1 (TRAF1) is a signaling adaptor first identified as part of the TNFR2 signaling complex. TRAF1 plays a key role in pro-survival signaling downstream of TNFR superfamily members such as TNFR2, LMP1, 4-1BB, and CD40. Recent studies have uncovered another role for TRAF1, independent of its role in TNFR superfamily signaling, in negatively regulating Toll-like receptor and Nod-like receptor signaling, through sequestering the linear ubiquitin assembly complex, LUBAC. TRAF1 has diverse roles in human disease. TRAF1 is overexpressed in many B cell related cancers and single nucleotide polymorphisms (SNPs) in TRAF1 have been linked to non-Hodgkin's lymphoma. Genome wide association studies have identified an association between SNPs in the 5' untranslated region of the TRAF1 gene with increased incidence and severity of rheumatoid arthritis and other rheumatic diseases. The loss of TRAF1 from chronically stimulated CD8 T cells results in desensitization of the 4-1BB signaling pathway, thereby contributing to T cell exhaustion during chronic infection. These apparently opposing roles of TRAF1 as both a positive and negative regulator of immune signaling have led to some confusion in the literature. Here we review the role of TRAF1 as a positive and negative regulator in different signaling pathways. Then we discuss the role of TRAF1 in human disease, attempting to reconcile seemingly contradictory roles based on current knowledge of TRAF1 signaling and biology. We also discuss avenues for future research to further clarify the impact of TRAF1 in human disease.

TRAF molecules in inflammation and inflammatory diseases

Purpose of Review

This review presents an overview of the current knowledge of TRAF molecules in inflammation with an emphasis on available human evidence and direct in vivo evidence of mouse models that demonstrate the contribution of TRAF molecules in the pathogenesis of inflammatory diseases.

Recent Findings

The tumor necrosis factor receptor (TNF-R)-associated factor (TRAF) family of cytoplasmic proteins was initially identified as signaling adaptors that bind directly to the intracellular domains of receptors of the TNF-R superfamily. It is now appreciated that TRAF molecules are widely employed in signaling by a variety of adaptive and innate immune receptors as well as cytokine receptors. TRAF-dependent signaling pathways typically lead to the activation of nuclear factor-κBs (NF-κBs), mitogen-activated protein kinases (MAPKs), or interferon-regulatory factors (IRFs). Most of these signaling pathways have been linked to inflammation, and therefore TRAF molecules were expected to regulate inflammation and inflammatory responses since their discovery in 1990s. However, direct in vivo evidence of TRAFs in inflammation and especially in inflammatory diseases had been lacking for many years, partly due to the difficulty imposed by early lethality of TRAF2^-/-, TRAF3^-/-, and TRAF6^-/- mice. With the creation of conditional knockout and lineage-specific transgenic mice of different TRAF molecules, our understanding about TRAFs in inflammation and inflammatory responses has rapidly advanced during the past decade.

Summary

Increasing evidence indicates that TRAF molecules are versatile and indispensable regulators of inflammation and inflammatory responses and that aberrant expression or function of TRAFs contributes to the pathogenesis of inflammatory diseases.

Systematic review and meta-analysis: pharmacogenetics of anti-TNF treatment response in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic inflammatory disease that affects ~1% of the Caucasian population. Over the last decades, the availability of biological drugs targeting the proinflammatory cytokine tumour necrosis factor α, anti-TNF drugs, has improved the treatment of patients with RA. However, one-third of the patients do not respond to the treatment. We wanted to evaluate the status of pharmacogenomics of anti-TNF treatment. We performed a PubMed literature search and all studies reporting original data on associations between genetic variants and anti-TNF treatment response in RA patients were included and results evaluated by meta-analysis. In total, 25 single nucleotide polymorphisms were found to be associated with anti-TNF treatment response in RA (19 from genome-wide association studies and 6 from the meta-analyses), and these map to genes involved in T cell function, NFκB and TNF signalling pathways (including CTCN5, TEC, PTPRC, FCGR2A, NFKBIB, FCGR2A, IRAK3). Explorative prediction analyses found that biomarkers for clinical treatment selection are not yet available.

Genetics of rheumatoid arthritis susceptibility, severity, and treatment response

A decade after the first genome-wide association study in rheumatoid arthritis (RA), a plethora of genetic association studies have been published on RA and its clinical or serological subtypes. We review the major milestones in the study of the genetic architecture of RA susceptibility, severity, and response to treatment. We set the scientific context necessary for non-geneticists to understand the potential clinical applications of human genetics and its significance for a stratified approach to the management of RA in the future.

Polymorphisms in STAT4, PTPN2, PSORS1C1 and TRAF3IP2 Genes Are Associated with the Response to TNF Inhibitors in Patients with Rheumatoid Arthritis

Objective

Rheumatoid Arthritis (RA) is a progressive autoimmune disease characterized by chronic joint inflammation and structural damage. Remission or at least low disease activity (LDA) represent potentially desirable goals of RA treatment. Single nucleotide polymorphisms (SNPs) in several genes might be useful for prediction of response to therapy. We aimed at exploring 4 SNPs in candidate genes (STAT4, PTPN2, PSORS1C1 and TRAF3IP2) in order to investigate their potential role in the response to therapy with tumor necrosis factor inhibitors (TNF-i) in RA patients.

Methods

In 171 RA patients we investigated the following SNPs: rs7574865 (STAT4), rs2233945 (PSORS1C1), rs7234029 (PTPN2) and rs33980500 (TRAF3IP2). Remission, LDA, and EULAR response were registered at 6 months and 2 years after initiation of first line TNF-i [Adalimumab (ADA) and Etanercept (ETN)].

Results

STAT4 variant allele was associated with the absence of a good/moderate EULAR response at 2 years of treatment in the whole RA group and in ETN treated patients. The PTPN2 SNP was associated with no good/moderate EULAR response at 6 months in ADA treated patients. Patients carrying PSORS1C1 variant allele did not reach LDA at 6 months in both the whole RA group and ETN treated patients. TRAF3IP2 variant allele was associated with the lack of LDA and remission achievement at 6 months in all RA cohort while an association with no EULAR response at 2 years of treatment occurred only in ETN treated patients.

Conclusions

For the first time, we reported that SNPs in STAT4, PTPN2, PSORS1C1, and TRAF3IP2 are associated with response to TNF-i treatment in RA patients; however, these findings should be validated in a larger population.

A longitudinal genome-wide association study of anti-tumor necrosis factor response among Japanese patients with rheumatoid arthritis

Background

Studies of Caucasian patients with rheumatoid arthritis (RA) to identify genetic biomarkers of anti-tumor necrosis factor (TNF) response have used response at a single time point as the phenotype with which single nucleotide polymorphism (SNP) associations have been tested. The findings have been inconsistent across studies. Among Japanese patients, only a few SNPs have been investigated. We report here the first genome-wide association study (GWAS) to identify genetic biomarkers of anti-TNF response among Japanese RA patients, using response at 2 time-points for a more reliable clinical phenotype over time.

Methods

Disease Activity Scores based on 28 joint counts (DAS28) were assessed at baseline (before initial therapy), and after 3 and 6 months in 487 Japanese RA patients starting anti-TNF therapy for the first time or switching to a new anti-TNF agent. A genome-wide panel of SNPs was genotyped and additional SNPs were imputed. Using change in DAS28 scores from baseline at both 3 (ΔDAS-3) and 6 months (ΔDAS-6) as the response phenotype, a longitudinal genome-wide association analysis was conducted using generalized estimating equations (GEE) models, adjusting for baseline DAS28, treatment duration, type of anti-TNF agent and concomitant methotrexate. Cross-sectional analyses were performed using multivariate linear regression models, with response from a single time point (ΔDAS-3 or ΔDAS-6) as phenotype; all other variables were the same as in the GEE models.

Results

In the GEE models, borderline significant association was observed at 3 chromosomal regions (6q15: rs284515, p = 6.6x10⁻⁷; 6q27: rs75908454, p = 6.3x10⁻⁷ and 10q25.3: rs1679568, p = 8.1x10⁻⁷), extending to numerous SNPs in linkage disequilibrium (LD) across each region. Potential candidate genes in these regions include MAP3K7, BACH2 (6q15), GFRA1 (10q25.3), and WDR27 (6q27). The association at GFRA1 replicates a previous finding from a Caucasian dataset. In the cross-sectional analyses, ΔDAS-6 was significantly associated with the 6q15 locus (rs284511, p = 2.5x10⁻⁸). No other significant or borderline significant associations were identified.

Conclusion

Three genomic regions demonstrated significant or borderline significant associations with anti-TNF response in our dataset of Japanese RA patients, including a locus previously associated among Caucasians. Using repeated measures of response as phenotype enhanced the power to detect these associations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-016-0920-6) contains supplementary material, which is available to authorized users.