Top

02-02-2017 | Rheumatoid arthritis | Article

Discovery of new candidate genes for rheumatoid arthritis through integration of genetic association data with expression pathway analysis

Journal: Arthritis Research & Therapy

Authors: Klementy Shchetynsky, Lina-Marcella Diaz-Gallo, Lasse Folkersen, Aase Haj Hensvold, Anca Irinel Catrina, Louise Berg, Lars Klareskog, Leonid Padyukov

Publisher: BioMed Central

Abstract

Background

Here we integrate verified signals from previous genetic association studies with gene expression and pathway analysis for discovery of new candidate genes and signaling networks, relevant for rheumatoid arthritis (RA).

Method

RNA-sequencing-(RNA-seq)-based expression analysis of 377 genes from previously verified RA-associated loci was performed in blood cells from 5 newly diagnosed, non-treated patients with RA, 7 patients with treated RA and 12 healthy controls. Differentially expressed genes sharing a similar expression pattern in treated and untreated RA sub-groups were selected for pathway analysis. A set of “connector” genes derived from pathway analysis was tested for differential expression in the initial discovery cohort and validated in blood cells from 73 patients with RA and in 35 healthy controls.

Results

There were 11 qualifying genes selected for pathway analysis and these were grouped into two evidence-based functional networks, containing 29 and 27 additional connector molecules. The expression of genes, corresponding to connector molecules was then tested in the initial RNA-seq data. Differences in the expression of ERBB2, TP53 and THOP1 were similar in both treated and non-treated patients with RA and an additional nine genes were differentially expressed in at least one group of patients compared to healthy controls. The ERBB2, TP53. THOP1 expression profile was successfully replicated in RNA-seq data from peripheral blood mononuclear cells from healthy controls and non-treated patients with RA, in an independent collection of samples.

Conclusion

Integration of RNA-seq data with findings from association studies, and consequent pathway analysis implicate new candidate genes, ERBB2, TP53 and THOP1 in the pathogenesis of RA.

de Klerk E, t Hoen PA. Alternative mRNA transcription, processing, and translation: insights from RNA sequencing. Trends Genet. 2015;31:128–39.CrossRefPubMed

Finotello F, Di Camillo B. Measuring differential gene expression with RNA-seq: challenges and strategies for data analysis. Brief Funct Genomics. 2015;14:130–42.CrossRefPubMed

Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee HS, Jia X, et al. Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nat Genet. 2012;44:291–6.CrossRefPubMedPubMedCentral

Padyukov L, Silva C, Stolt P, Alfredsson L, Klareskog L. A gene-environment interaction between smoking and shared epitope genes in HLA-DR provides a high risk of seropositive rheumatoid arthritis. Arthritis Rheum. 2004;50:3085–92.CrossRefPubMed

Klareskog L, Stolt P, Lundberg K, Kallberg H, Bengtsson C, Grunewald J, et al. A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination. Arthritis Rheum. 2006;54:38–46.CrossRefPubMed

Tedeschi SK, Bermas B, Costenbader KH. Sexual disparities in the incidence and course of SLE and RA. Clin Immunol. 2013;149:211–8.CrossRefPubMed

Barragan-Martinez C, Amaya-Amaya J, Pineda-Tamayo R, Mantilla RD, Castellanos-dela Hoz J, Bernal-Macias S, et al. Gender differences in Latin-American patients with rheumatoid arthritis. Gend Med. 2012;9:490–510. e5.CrossRefPubMed

Freudenberg J, Lee HS, Han BG, Shin HD, Kang YM, Sung YK, et al. Genome-wide association study of rheumatoid arthritis in Koreans: population-specific loci as well as overlap with European susceptibility loci. Arthritis Rheum. 2011;63:884–93.CrossRefPubMed

Okada Y, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, et al. Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population. Nat Genet. 2012;44:511–6.CrossRefPubMed

10.

Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, et al. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat Genet. 2010;42:508–14.CrossRefPubMedPubMedCentral

11.

Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. 2012;44:1336–40.CrossRefPubMedPubMedCentral

12.

Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014;506:376–81.CrossRefPubMed

13.

Sahlen P, Abdullayev I, Ramskold D, Matskova L, Rilakovic N, Lotstedt B, et al. Genome-wide mapping of promoter-anchored interactions with close to single-enhancer resolution. Genome Biol. 2015;16:156.CrossRefPubMedPubMedCentral

14.

Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.CrossRefPubMedPubMedCentral

15.

NCBI Gene Expression Omnibus, project number: GSE90081. http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE90081. Accessed 24 Nov 2016.

16.

Folkersen L, Brynedal B, Diaz-Gallo LM, Ramskold D, Shchetynsky K, Westerlind H, et al. Integration of known DNA, RNA and protein biomarkers provides prediction of anti-TNF response in rheumatoid arthritis: results from the COMBINE study. Mol Med. 2016;22:322–8.

17.

Chen J, Weiss WA. Alternative splicing in cancer: implications for biology and therapy. Oncogene. 2015;34:1–14.CrossRefPubMed

18.

Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.CrossRefPubMedPubMedCentral

19.

Yosef BYH. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Statist Soc Ser. 1995;B 57:289–300.

20.

Abe T, Takeuchi T. Rheumatoid arthritis and tumor necrosis factor alpha. Autoimmunity. 2001;34:291–303.CrossRefPubMed

21.

Raanani P, Ben-Bassat I. Immune-mediated complications during interferon therapy in hematological patients. Acta Haematol. 2002;107:133–44.CrossRefPubMed

22.

Moots RJ, Naisbett-Groet B. The efficacy of biologic agents in patients with rheumatoid arthritis and an inadequate response to tumour necrosis factor inhibitors: a systematic review. Rheumatology. 2012;51:2252–61.CrossRefPubMed

23.

Li Q, Hu J, Ding J, Zheng G. Fisher’s method of combining dependent statistics using generalizations of the gamma distribution with applications to genetic pleiotropic associations. Biostatistics. 2014;15:284–95.CrossRefPubMed

24.

Infante-Rivard C, Mirea L, Bull SB. Combining case-control and case-trio data from the same population in genetic association analyses: overview of approaches and illustration with a candidate gene study. Am J Epidemiol. 2009;170:657–64.CrossRefPubMed

25.

Hess A, Iyer H. Fisher’s combined p-value for detecting differentially expressed genes using Affymetrix expression arrays. BMC Genomics. 2007;8:96.CrossRefPubMedPubMedCentral

26.

Tervaniemi MH, Katayama S, Skoog T, Siitonen HA, Vuola J, Nuutila K, et al. NOD-like receptor signaling and inflammasome-related pathways are highlighted in psoriatic epidermis. Sci Rep. 2016;6:22745.CrossRefPubMedPubMedCentral

27.

Skopelja S, Hamilton BJ, Jones JD, Yang ML, Mamula M, Ashare A, et al. The role for neutrophil extracellular traps in cystic fibrosis autoimmunity. JCI Insight. 2016;1:e88912.CrossRefPubMedPubMedCentral

28.

Cheung AK, Ko JM, Lung HL, Chan KW, Stanbridge EJ, Zabarovsky E, et al. Cysteine-rich intestinal protein 2 (CRIP2) acts as a repressor of NF-kappaB-mediated proangiogenic cytokine transcription to suppress tumorigenesis and angiogenesis. Proc Natl Acad Sci U S A. 2011;108:8390–5.CrossRefPubMedPubMedCentral

29.

Kennedy JM, Fodil N, Torre S, Bongfen SE, Olivier JF, Leung V, et al. CCDC88B is a novel regulator of maturation and effector functions of T cells during pathological inflammation. J Exp Med. 2014;211:2519–35.CrossRefPubMedPubMedCentral

30.

Sun Y, Cheung HS. p53, proto-oncogene and rheumatoid arthritis. Semin Arthritis Rheum. 2002;31:299–310.CrossRefPubMed

31.

Maas K, Westfall M, Pietenpol J, Olsen NJ, Aune T. Reduced p53 in peripheral blood mononuclear cells from patients with rheumatoid arthritis is associated with loss of radiation-induced apoptosis. Arthritis Rheum. 2005;52:1047–57.CrossRefPubMed

32.

Kitamura T, Sekimata M, Kikuchi S, Homma Y. Involvement of poly(ADP-ribose) polymerase 1 in ERBB2 expression in rheumatoid synovial cells. Am J Physiol Cell Physiol. 2005;289:C82–8.CrossRefPubMed

33.

Yamaoka T, Yan F, Cao H, Hobbs SS, Dise RS, Tong W, et al. Transactivation of EGF receptor and ErbB2 protects intestinal epithelial cells from TNF-induced apoptosis. Proc Natl Acad Sci U S A. 2008;105:11772–7.CrossRefPubMedPubMedCentral

34.

Omenn GS, Guan Y, Menon R. A new class of protein cancer biomarker candidates: differentially expressed splice variants of ERBB2 (HER2/neu) and ERBB1 (EGFR) in breast cancer cell lines. J Proteomics. 2014;107:103–12.CrossRefPubMed

35.

Saric T, Beninga J, Graef CI, Akopian TN, Rock KL, Goldberg AL. Major histocompatibility complex class I-presented antigenic peptides are degraded in cytosolic extracts primarily by thimet oligopeptidase. J Biol Chem. 2001;276:36474–81.CrossRefPubMed