The gene expression profile in the synovium as a predictor of the clinical response to infliximab treatment in rheumatoid arthritis

Johan Lindberg, Carla A Wijbrandts, Lisa G van Baarsen, Gustavo Nader, Lars Klareskog, Anca Catrina, Rogier Thurlings, Margriet Vervoordeldonk, Joakim Lundeberg, Paul P Tak, Johan Lindberg, Carla A Wijbrandts, Lisa G van Baarsen, Gustavo Nader, Lars Klareskog, Anca Catrina, Rogier Thurlings, Margriet Vervoordeldonk, Joakim Lundeberg, Paul P Tak

Abstract

Background: Although the use of TNF inhibitors has fundamentally changed the way rheumatoid arthritis (RA) is treated, not all patients respond well. It is desirable to facilitate the identification of responding and non-responding patients prior to treatment, not only to avoid unnecessary treatment but also for financial reasons. In this work we have investigated the transcriptional profile of synovial tissue sampled from RA patients before anti-TNF treatment with the aim to identify biomarkers predictive of response.

Methodology/principal findings: Synovial tissue samples were obtained by arthroscopy from 62 RA patients before the initiation of infliximab treatment. RNA was extracted and gene expression profiling was performed using an in-house spotted long oligonucleotide array covering 17972 unique genes. Tissue sections were also analyzed by immunohistochemistry to evaluate cell infiltrates. Response to infliximab treatment was assessed according to the EULAR response criteria. The presence of lymphocyte aggregates dominated the expression profiles and a significant overrepresentation of lymphocyte aggregates in good responding patients confounded the analyses. A statistical model was set up to control for the effect of aggregates, but no differences could be identified between responders and non-responders. Subsequently, the patients were split into lymphocyte aggregate positive- and negative patients. No statistically significant differences could be identified except for 38 transcripts associated with differences between good- and non-responders in aggregate positive patients. A profile was identified in these genes that indicated a higher level of metabolism in good responding patients, which indirectly can be connected to increased inflammation.

Conclusions/significance: It is pivotal to account for the presence of lymphoid aggregates when studying gene expression patterns in rheumatoid synovial tissue. In spite of our original hypothesis, the data do not support the notion that microarray analysis of whole synovial biopsy specimens can be used in the context of personalized medicine to identify non-responders to anti-TNF therapy before the initiation of treatment.

Conflict of interest statement

Competing Interests: P.P. Tak has served as a consultant for Abbott, Amgen, Centocor, Schering-Plough, UCB, and Wyeth. Patients were treated with infliximab, a commercial antibody marketed by Schering-Plough. There are no patents, marketed products, or products in development related to this research. Also, the use of infliximab does not alter the adherence of this work to all the PLoS ONE policies on sharing data and materials.

Figures

Figure 1. Lymphocyte aggregates.
Figure 1. Lymphocyte aggregates.
The presence of lymphocyte aggregates was assessed on anti-CD3-stained sections. Aggregates were counted and graded by size. Aggregate size was assessed by counting the number of cells in a radius starting from an estimated center of the aggregate. Aggregate size was then classified as grade 1 (1–5 cells in the radius [panel B]), grade 2 (5–10 cells in the radius [panel C]), or grade 3 (>10 cells in the radius [panel D]). Tissue sections with no lymphocyte aggregates were graded as 0 (panel A).
Figure 2. Hierarchical cluster of patients using…
Figure 2. Hierarchical cluster of patients using features with high variation irrespective of response group.
Patients were clustered in a hierarchical dendrogram using features with a log2-ratio (ratio = sample intensity/reference intensity) interquartile range >1 (1321 features in total). Abbreviations in the dendrogram: PXX = Patient number; G,  = Good responder; M = Moderate responder; N = Non responder; L+/− = presence of large lymphocyte aggregates; S+/−presence of small aggregates. Colors in the dendrograms: Red = patients positive for either small or large lymphocyte aggregates; Blue = patients negative for small and large lymphocyte aggregates; Black = patients with no lymphocyte aggregate assessment.

References

    1. Tak PP, Bresnihan B. The pathogenesis and prevention of joint damage in rheumatoid arthritis: advances from synovial biopsy and tissue analysis. Arthritis Rheum. 2000;43:2619–2633.
    1. van der Helm-van Mil AH, Huizinga TW, de Vries RR, Toes RE. Emerging patterns of risk factor make-up enable subclassification of rheumatoid arthritis. Arthritis Rheum. 2007;56:1728–1735.
    1. van der Pouw Kraan TC, van Gaalen FA, Kasperkovitz PV, Verbeet NL, Smeets TJ, et al. Rheumatoid arthritis is a heterogeneous disease: evidence for differences in the activation of the STAT-1 pathway between rheumatoid tissues. Arthritis Rheum. 2003;48:2132–2145.
    1. Tracey D, Klareskog L, Sasso EH, Salfeld JG, Tak PP. Tumor necrosis factor antagonist mechanisms of action: a comprehensive review. Pharmacol Ther. 2008;117:244–279.
    1. Listing J, Strangfeld A, Kary S, Rau R, von Hinueber U, et al. Infections in patients with rheumatoid arthritis treated with biologic agents. Arthritis Rheum. 2005;52:3403–3412.
    1. Braun-Moscovici Y, Markovits D, Zinder O, Schapira D, Rozin A, et al. Anti-cyclic citrullinated protein antibodies as a predictor of response to anti-tumor necrosis factor-alpha therapy in patients with rheumatoid arthritis. J Rheumatol. 2006;33:497–500.
    1. Hyrich KL, Watson KD, Silman AJ, Symmons DP. Predictors of response to anti-TNF-alpha therapy among patients with rheumatoid arthritis: results from the British Society for Rheumatology Biologics Register. Rheumatology (Oxford) 2006;45:1558–1565.
    1. Lequerre T, Jouen F, Brazier M, Clayssens S, Klemmer N, et al. Autoantibodies, metalloproteinases and bone markers in rheumatoid arthritis patients are unable to predict their responses to infliximab. Rheumatology (Oxford) 2007;46:446–453.
    1. Ulfgren AK, Andersson U, Engstrom M, Klareskog L, Maini RN, et al. Systemic anti-tumor necrosis factor alpha therapy in rheumatoid arthritis down-regulates synovial tumor necrosis factor alpha synthesis. Arthritis Rheum. 2000;43:2391–2396.
    1. Wijbrandts CA, Dijkgraaf MG, Kraan MC, Vinkenoog M, Smeets TJ, et al. The clinical response to infliximab in rheumatoid arthritis is in part dependent on pretreatment tumour necrosis factor alpha expression in the synovium. Ann Rheum Dis. 2008;67:1139–1144.
    1. Marotte H, Maslinski W, Miossec P. Circulating tumour necrosis factor-alpha bioactivity in rheumatoid arthritis patients treated with infliximab: link to clinical response. Arthritis Res Ther. 2005;7:R149–155.
    1. Wolbink GJ, Voskuyl AE, Lems WF, de Groot E, Nurmohamed MT, et al. Relationship between serum trough infliximab levels, pretreatment C reactive protein levels, and clinical response to infliximab treatment in patients with rheumatoid arthritis. Ann Rheum Dis. 2005;64:704–707.
    1. Koczan D, Drynda S, Hecker M, Drynda A, Guthke R, et al. Molecular discrimination of responders and nonresponders to anti-TNFalpha therapy in rheumatoid arthritis by etanercept. Arthritis Res Ther. 2008;10:R50.
    1. Lequerre T, Gauthier-Jauneau AC, Bansard C, Derambure C, Hiron M, et al. Gene profiling in white blood cells predicts infliximab responsiveness in rheumatoid arthritis. Arthritis Res Ther. 2006;8:R105.
    1. Sekiguchi N, Kawauchi S, Furuya T, Inaba N, Matsuda K, et al. Messenger ribonucleic acid expression profile in peripheral blood cells from RA patients following treatment with an anti-TNF-alpha monoclonal antibody, infliximab. Rheumatology (Oxford) 2008;47:780–788.
    1. Michiels S, Koscielny S, Hill C. Prediction of cancer outcome with microarrays: a multiple random validation strategy. Lancet. 2005;365:488–492.
    1. Lindberg J, af Klint E, Catrina AI, Nilsson P, Klareskog L, et al. Effect of infliximab on mRNA expression profiles in synovial tissue of rheumatoid arthritis patients. Arthritis Res Ther. 2006;8:R179.
    1. van der Pouw Kraan TC, Wijbrandts CA, van Baarsen LG, Rustenburg F, Baggen JM, et al. Responsiveness to anti-tumour necrosis factor alpha therapy is related to pre-treatment tissue inflammation levels in rheumatoid arthritis patients. Ann Rheum Dis. 2008;67:563–566.
    1. Badot V, Galant C, Nzeusseu Toukap A, Theate I, Maudoux AL, et al. Gene expression profiling in the synovium identifies a predictive signature of absence of response to adalimumab therapy in rheumatoid arthritis. Arthritis Res Ther. 2009;11:R57.
    1. Felson DT, Anderson JJ, Boers M, Bombardier C, Furst D, et al. American College of Rheumatology. Preliminary definition of improvement in rheumatoid arthritis. Arthritis Rheum. 1995;38:727–735.
    1. van Gestel AM, Prevoo ML, van 't Hof MA, van Rijswijk MH, van de Putte LB, et al. Development and validation of the European League Against Rheumatism response criteria for rheumatoid arthritis. Comparison with the preliminary American College of Rheumatology and the World Health Organization/International League Against Rheumatism Criteria. Arthritis Rheum. 1996;39:34–40.
    1. Kraan MC, Reece RJ, Smeets TJ, Veale DJ, Emery P, et al. Comparison of synovial tissues from the knee joints and the small joints of rheumatoid arthritis patients: Implications for pathogenesis and evaluation of treatment. Arthritis Rheum. 2002;46:2034–2038.
    1. Thurlings RM, Wijbrandts CA, Mebius RE, Cantaert T, Dinant HJ, et al. Synovial lymphoid neogenesis does not define a specific clinical rheumatoid arthritis phenotype. Arthritis Rheum. 2008;58:1582–1589.
    1. Lindberg J, Gry M, Klevebring D, Wirta W. 2005. KTH microarray core facility webpage.
    1. Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res. 2005;33:D54–58.
    1. Lindberg J, af Klint E, Ulfgren AK, Stark A, Andersson T, et al. Variability in synovial inflammation in rheumatoid arthritis investigated by microarray technology. Arthritis Res Ther. 2006;8:R47.
    1. Yang YH, Speed T. Design issues for cDNA microarray experiments. Nat Rev Genet. 2002;3:579–588.
    1. Team RDC. R: A language and environment for statistical computing. Vienna, Austria: 2008.
    1. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
    1. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–5121.
    1. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
    1. Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–14868.
    1. Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–258.
    1. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–29.
    1. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–484.
    1. Kim SY, Volsky DJ. PAGE: parametric analysis of gene set enrichment. BMC Bioinformatics. 2005;6:144.
    1. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–15550.
    1. Reiner A, Yekutieli D, Benjamini Y. Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003;19:368–375.
    1. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, et al. NCBI GEO: mining tens of millions of expression profiles–database and tools update. Nucleic Acids Res. 2007;35:D760–765.
    1. Klaasen R, Thurlings RM, Wijbrandts CA, van Kuijk AWR, Baeten D, et al. The Relationship between Synovial Lymphocyte Aggregates and the Clinical Response to Infliximab in Rheumatoid Arthritis: a Prospective Study. 2009. Submitted for publication.
    1. van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415:530–536.
    1. Cantaert T, Kolln J, Timmer T, van der Pouw Kraan TC, Vandooren B, et al. B lymphocyte autoimmunity in rheumatoid synovitis is independent of ectopic lymphoid neogenesis. J Immunol. 2008;181:785–794.
    1. van Vollenhoven RF, Klareskog L. Clinical responses to tumor necrosis factor alpha antagonists do not show a bimodal distribution: data from the Stockholm tumor necrosis factor alpha followup registry. Arthritis Rheum. 2003;48:1500–1503.
    1. Smolen JS, Han C, Bala M, Maini RN, Kalden JR, et al. Evidence of radiographic benefit of treatment with infliximab plus methotrexate in rheumatoid arthritis patients who had no clinical improvement: a detailed subanalysis of data from the anti-tumor necrosis factor trial in rheumatoid arthritis with concomitant therapy study. Arthritis Rheum. 2005;52:1020–1030.
    1. Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S. Stochastic mRNA synthesis in mammalian cells. PLoS Biol. 2006;4:e309.
    1. Shankavaram UT, Reinhold WC, Nishizuka S, Major S, Morita D, et al. Transcript and protein expression profiles of the NCI-60 cancer cell panel: an integromic microarray study. Mol Cancer Ther. 2007;6:820–832.

Source: PubMed

3
Suscribir