Prior to extension, Transcriptomes of fibroblast-like Synoviocytes from extended and Polyarticular juvenile idiopathic arthritis are indistinguishable
© The Author(s). 2018
Received: 21 September 2017
Accepted: 19 December 2017
Published: 8 January 2018
Our intent was to identify differences between the transcriptome of fibroblast-like synoviocytes (FLS) in oligoarticular juvenile idiopathic arthritis (JIA) before extension when compared to persistent subtype of JIA, when the two are clinically indistinguishable. Additionally, we sought to determine if differences between the transcriptomes of FLS from extended-to-be and polyarticular course JIA could be detected. Our hypothesis was that intrinsic differences in the transcriptome of the FLS from extended-to-be JIA would distinguish them from persistent oligoarticular JIA, before the course is clinically apparent.
Global gene expression was defined in cultured FLS from 6 controls, 12 JIA with persistent course, 7 JIA prior to extension (extended-to-be), 4 JIA with extended course and 6 polyarticular onset, using Affymetrix Human GeneChips 133plus2.0.
Bioconductor Linear Models for Microarray Analysis revealed 22 probesets with differential expression between persistent and extended-to-be FLS at 15% FDR, however only 2 probesets distinguished extended-to-be from extended and none distinguished extended-to-be and polyarticular at 15% FDR. Differences in extended and polyarticular gene expression profiles were not detected. Confirmation of select genes was done on the RNA level by RT-qPCR and on the protein level in synovial fluid by ELISA.
The transcriptome of FLS from extended-to-be juvenile idiopathic arthritis is distinct from persistent course before a clinical distinction can be made. Additionally, the transcriptome of extended-to-be and polyarticular course, including those who have already extended, are indistinguishable. These gene expression data suggest that FLS already reflect a polyarticular behavior early in disease course, suggesting that extended-to-be may be “latent polyarticular” at onset. These differences can be used to develop early biomarkers of disease course, allowing for better-informed treatment decisions.
Juvenile Idiopathic Arthritis (JIA), the most common rheumatic disease of childhood , carries risk of permanent disability. Oligoarticular subtype is the most benign, except when it evolves into polyarticular course, termed extended oligoarticular disease. Persistent oligoarthritis affects up to 4 joints, while extended oligoarthritis affects a cumulative total of 5 or more joints after the first 6 months. Between 21 and 50% of patients show extension to polyarticular course, of whom only 13–23% achieve remission . Functional outcome for children with extended oligoarticular JIA can be poor and the key to achieving improvement may be in prediction of extension so that early treatment can be initiated .
In rheumatoid arthritis (RA), fibroblast-like synoviocytes (FLS) are cytokine producing key effector cells that both directly destroy cartilage and promote perpetuation of inflammation . These cells have a central role in the pathogenesis of rheumatoid synovitis and may provide insight into differences in the pathogenic processes that allow extension to polyarticular course in JIA. We hypothesize that FLS contain central signals involved in the “switch” from oligoarticular to polyarticular course.
This is the first study to look at the comparative transcriptomes of JIA FLS of different disease subtypes. Transcriptome analysis revealed important differences in gene expression of FLS from oligoarticular JIA prior to extension (extended-to-be) when compared to that of FLS from persistent oligoarticular subtype. Moreover, no differences were detected in the transcriptomes of FLS from polyarticular courses (extended and polyarticular subtypes) and FLS of extended-to-be JIA, early in disease course. Differences in JIA FLS transcriptomes between persistent and extended-to-be challenge the classification criteria and can provide clues for candidate biomarker development to predict disease course in JIA.
Patients and samples
As part of an ongoing Institutional Review Board (IRB) -approved protocol, remnant synovial fluids/tissues were obtained from arthrocenteses and synovectomies from five groups:  controls (C),  persistent oligoarticular (PR),  oligoarticular prior to extension (Ext-to-be),  oligoarticular with extended course (E), and  polyarticular onset (Poly). Medical record review was conducted to select appropriate samples from our synovial fluid and tissue repository. All patients met the ILAR classification criteria for JIA , including persistent and extended oligoarticular and polyarticular subtypes, specifically excluding enthesitis-related and psoriatic subtypes. Clinical information on individual patients is in supplemental table (Additional file 1 Clinical Table). There were 6 C, 12 PR, 7 Ext-to-be, 4 E and 6 Poly included in the study. All control samples were obtained from orthopedic procedures. Patients with persistent oligoarticular course were followed for between 2 and 18 years from diagnosis without extension. The majority of patients were antinuclear antibody (ANA) positive. All patients who had documented testing for HLA-B27 allele and cyclic citrullinated peptide (CCP) antibody were negative. Of those tested, one patient with extended oligoarticular course was rheumatoid factor (RF) positive.
Synovial fluid/tissue was plated in 6-well plates with 15% fetal bovine serum/Dulbecco’s modified Eagle’s medium to establish primary cultures of FLS, passaged 3–6 times, then harvested at confluence. By the third passage, FLS are the predominant cell type in culture . As indicated in the supplemental table, four of the control FLS were grown from synovial tissue and the remainder of the FLS were grown from synovial fluid.
RNA was extracted from cultured FLS using TRIzol, purified using Qiagen RNeasy Mini Prep kit, amplified, fluorescence labeled and hybridized to GeneChip Human Genome U133 Plus 2.0 Arrays per Affymetrix protocol. Arrays were stained using GeneChip Hybridization, Wash, and Stain kit (Affymetrix), and scanned and analyzed using Affymetrix command and expression console software. Log2 expression values were calculated for the 54,675 probesets with GeneChip Robust Multiarray Analysis (GC-RMA), a method of microarray normalization and summarization, using the gcRMA Bioconductor package .
Gene expression analysis
Data was filtered for log2 expression >4 in all hybridizations in at least one group in each pairwise comparison, then for |1.5|-fold change. Probesets with a statistically significant difference in expression in pairwise comparisons were identified using Bioconductor package Linear Models for Microarray Analysis (LIMMA) with Benjamini and Hochberg’s method to control the false-discovery rate (FDR) , with significance of 15%. Hierarchical clustering was performed using Euclidean correlation and average linkage in MeV (MultiExperiment Viewer) .
Quantitative reverse transcription-polymerase chain reaction (qRT-PCR)
Agilent Bioanalzer and Nanodrop 2000 were used for quantification of RNA. High-Capacity cDNA RT Kit from LifeTechnologies was used (manufacturer’s procedures) with ABI Veriti thermal cycler. After reverse transcription, 20ul nuclease-free water was added (total volume 40ul). cDNA was assayed on 7900HT Real-Time PCR system with TaqMan Gene Expression Assays. TaqMan IDs for each gene are as follows: CLDN11 Hs00194440_m1, ANKRD44 Hs00403517_m1, ICAM2 Hs00609563_m1, KLHL13 Hs01006506_m1, LY6K Hs03988347_m1, KIF11 Hs00189698_m1, MAMLD1 Hs00193976_m1, RAB27B Hs00188156_m1, RGS2 Hs01009070_g1, ZNF204P Hs00394812_m1. Samples (in triplicate) on 384 well plate were cycled on 7900HT Real-Time PCR system following manufacturer’s default conditions using SDS software v2.4 relative quantification run. Data was analyzed using the RQ (relative quantification) Manager to calculate relative quantification using the delta-delta-Ct (ddCt) method .
Enzyme-linked immunosorbent assay (ELISA)
Using synovial fluid samples from 12 PR and 7 Ext-to-be, CD14 ELISA was completed using Quantikine Colormetric ELISA kits from R&D Systems (manufacturer’s protocol). For the other proteins, we performed standard sandwich ELISAs optimized by our lab . Reagents: MBP (Novus NB110-79873B, NBP2–46632), KLHL13 (SCBT 138378, bs-7752R), HspBAP1 (Novus NBP1–92014, LifeSpan LS-C371907–50), ANKRD44 (Novus NBP1–80887, LifeSpan LS-C251405–200). Averages were obtained using Excel. Standard errors and p-values were calculated using SigmaStat.
The transcriptome of FLS from extended-to-be JIA was distinct from the transcriptome of FLS of persistent oligoarticular JIA. We sought to identify potential candidate biomarkers for polyarticular course through differences in FLS gene expression. Gene expression profiling was done of cultured FLS from 6 controls (C), 12 persistent oligoarticular (PR), 7 extended-to-be, obtained prior to extension (Ext-to-be), 4 extended, obtained after disease had extended (E) and 6 polyarticular, with ≥5 joints involved in first 6 months of disease (Poly) samples. A large number of differences between controls and the JIA groups were identified. (Additional file 2 FDR_results_table).
Differentially Expressed Genes Between Persistent and Extended-to-be Oligoarticular JIA at 15% False Discovery Rate
Zinc finger protein 204, pseudogene
Regulator of G-protein signaling 2
Kelch-like family member 13
G protein-couple receptor associated sorting protein 1
Serine peptidase inhibitor, Kazal type 13 (putative)
Mastermind-like domain containing 1
Ankyrin repeat domain 44
Intercellular adhesion molecule 2
ATP binding cassette subfamily C member 6
Kinesin family member 11
RAB27B, member RAS oncogene family
Topoisomerase (DNA) II alpha
Lymphocyte antigen 6 complex, locus K
GEN1 Holliday junction 5′ flap endonuclease
Six probesets were more highly expressed in the persistent and 10 probesets more highly expressed in the extended-to-be FLS. RNA level confirmation by qRT-PCR was done for a representative group of 10 genes from the 16 differentially expressed genes (Fig. 1c). The real-time PCR confirmation was done on RNA from the same samples as used for the microarray analysis. The PCR results were consistent with the microarray finding for 8 of the 10 tested genes, with the exceptions of MAMLD1 and RAB27B.
No differences were detected between the transcriptome of FLS of extended-to-be oligoarticular JIA and that of FLS from polyarticular course JIA
Once we demonstrated that gene expression of FLS from Ext-to-be was distinct from that of FLS from persistent course, we sought to determine if FLS from Ext-to-be were distinguishable from those of polyarticular course. In a pairwise comparison of extended-to-be and extended samples, there were only 2 significant differentially expressed probesets by 15% FDR. Pairwise comparison of polyarticular versus either Ext-to-be or E revealed no differentially expressed probesets (lowest adjusted p-values 0.42 and 0.33, respectively).
Expression pattern of differentially expressed genes in FLS from polyarticular and extended JIA
Since there were no detectible differences between gene expression of Ext-to-be and polyarticular course, we turned our attention to the expression patterns of the 16 differentially expressed genes between PR and Ext-to-be. Patterns of expression of these genes in the E and Poly samples demonstrated close correlation with the pattern seen in the Ext-to-be when compared with PR (Fig. 1b).
Secreted proteins as candidate synovial fluid biomarkers for prediction of course in JIA
Our hypothesis was that intrinsic differences in the transcriptome of the FLS from extended-to-be JIA early in the course would distinguish them from persistent oligoarticular JIA, before the course is clinically apparent.
We performed gene expression profiling using FLS from controls and patients with differing disease courses, specifically persistent oligoarticular, extended-to-be and polyarticular course (extended and polyarticular onset). The extended-to-be samples were all taken from the very first sample available, which preceded extension in all patients. For the 16 differentially expressed genes, Ext-to-be already display expression patterns similar to the FLS of patients with polyarticular course. This pattern of gene expression is distinct from oligoarticular who don’t extend (persistent), demonstrating that there are biologic differences prior to detectible clinical differences between PR and Ext-to-be. Taken together, these data indicate a specific pattern of gene expression in FLS of polyarticular course (Poly and E) that is displayed by FLS of Ext-to-be even before that course is clinically apparent.
Ext-to-be FLS demonstrated increased expression of HLX, a homeobox transcription factor, which is a positive regulator of Th1 cell differentiation and a negative regulator of Th2 cell differentiation. Conversion of Th17 cells to non-classic Th1 phenotype has been demonstrated in the inflamed joint in JIA  and FLS may influence this altered T-cell repertoire. ICAM-2 (CD102) is an intercellular adhesion molecule that binds LFA-1 protein. It mediates adhesion and has a role in antigen-specific immune response, NK-cell mediated clearance and lymphocyte recirculation. Singh et al. have proposed ICAM-2 as a potential therapeutic target to inhibit FLS activation in RA as they demonstrated that T-cell induced activation of Akt in FLS is mediated by ICAM-2 . Unexpectedly, CD14 was increased in Ext-to-be FLS on the RNA level, which translated into increased protein in the synovial fluid. A novel procollagen + CD45 + CD14 + IL17RA + CD34- fibrocyte-like cell (FLC) has been reported in the synovial fluid of JIA and the authors postulated an FLC-CD8 T cell role in perpetuation of inflammation in JIA . In RA, increased soluble CD14 (sCD14) has been described in plasma and synovial fluid, with thoughts that sCD14 was produced by RA synovial macrophages through cleavage of membranous CD14. These authors postulated that RA synovial fibroblasts were sensitive to LPS in the presence of sCD14 and LPS-binding protein . One could argue that increased monocytes in Ext-to-be synovial fluid could increase CD14 level, however, the proportion of monocytes/macrophages in the synovial fluid mononuclear cells was not significantly different between PR and Ext-to-be , pointing towards the FLS as important contributors to differences in CD14 levels in the synovial fluid.
There are small numbers of samples and we may be underpowered to detect differences between certain groups, however, we did identify differences between the transcriptomes of PR and Ext-to-be FLS. This is the first study looking at direct patient samples of JIA FLS in this way. Even with small sample numbers in the comparisons, there were large numbers of differentially expressed genes between controls and each of the other groups. We were very stringent in our selection of samples to ensure that samples used met the ILAR criteria and extended-to-be were obtained prior to extension. In 8 of 10 genes tested by qRT-PCR, we were able to confirm the microarray findings. We removed XIST to remove some bias due to sex differences. Samples PR11 and PR12 cluster more closely with the Ext-to-be samples. PR11 was from a patient on naproxen who was followed for 9 years without extension. PR12 was from a patient who remained persistent in course (7 years follow up), however, was already on methotrexate. Ultimately, this patient required advanced therapy with anti-TNF agent for refractory uveitis, raising the question of altered course of arthritis through biologic therapy.
Although significant work has been done in FLS from RA, there is minimal comparable work done in JIA in general and none in the FLS from JIA. This is the first study to compare global gene expression of FLS from persistent, extended-to-be and polyarticular course JIA. Through this transcriptome level comparison, this study demonstrates that prior to extension, the FLS of JIA show a gene expression pattern that is more closely aligned to that of polyarticular course JIA than to that of persistent oligoarticular JIA. Using this information, we can develop prognostic synovial fluid biomarkers for predicting extension of oligoarticular JIA to polyarticular course, thereby allowing for more informed treatment decisions earlier in the course, perhaps improving outcomes for children with arthritis.
This correlation of early gene expression in the extended-to-be FLS with that of already extended and polyarticular FLS highlights the importance of FLS in reflecting ultimate course of disease. Rather than being persistent JIA that later evolves into an extended pattern, these data are the first pieces of evidence to suggest that extended-to-be patients may have “latent polyarticular” arthritis and perhaps should be treated more aggressively earlier in disease course. The FLS may hold the key for determining whether there is a pre-programmed prognosis or whether extended-to-be JIA already has smoldering polyarticular involvement that is below the limit of clinical detection. Our future goal is to develop a rapid throughput synovial fluid biomarker panel for predicting course in oligoarticular JIA. Demonstrating the translation of increased RNA expression into increased protein levels allows us to move forward with screening of potential synovial fluid biomarkers.
We would like to thank Edward Behrens MD for helpful discussions, Nemours/AI DuPont Hospital for Children, Division of Orthopedics, and the Nemours Biobank for acquisition of samples. We are grateful for the technical support of the COBRE-funded Biomolecular Core Laboratory, which is supported by grants P30GM114736 (COBRE) and P20GM103446 (INBRE) at Nemours Biomedical Research.
There was no financial support from commercial sources. This work was supported by the National Institutes of Health, National Institute of Arthritis, Musculoskeletal and Skin Research [grant number K23AR066724]; National Institutes of Health [grant number COBRE 8P20-GM-103464, ARRA P20-RR-020173], Arthritis Foundation Innovative Research Grant, The Nancy Taylor Foundation for Chronic Diseases, Delaware Community Foundation and The Open Net Foundation.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
All authors read and approved the final manuscript. ACB: conception and design of study, analysis and interpretation of data, drafting and revising manuscript.
MMS: acquisition, analysis and interpretation of data. SMM: analysis and interpretation of data. KES: interpretation of data, revising critically. CDR: conception and design of study, interpretation of data, revising critically.
Ethics approval and consent to participate
This study was approved by the Nemours Institutional Review Boards. Samples were selected from our repository, which is formed under a separate Nemours Institutional Review Boards protocol, with consents and assents.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lawrence RC, Helmick CG, Arnett FC, Deyo RA, Felson DT, Giannini EH, et al. Estimates of the prevalence of arthritis and selected musculoskeletal disorders in the United States. Arthritis Rheum. 1998;41(5):778–99. Epub 1998/05/20View ArticlePubMedGoogle Scholar
- Guillaume S, Prieur AM, Coste J, Job-Deslandre C. Long-term outcome and prognosis in oligoarticular-onset juvenile idiopathic arthritis. Arthritis Rheum. 2000;43(8):1858–65. Epub 2000/08/16View ArticlePubMedGoogle Scholar
- Foster HE, Marshall N, Myers A, Dunkley P, Griffiths ID. Outcome in adults with juvenile idiopathic arthritis: a quality of life study. Arthritis Rheum. 2003;48(3):767–75. Epub 2003/03/13View ArticlePubMedGoogle Scholar
- Huber LC, Distler O, Tarner I, Gay RE, Gay S, Pap T. Synovial fibroblasts: key players in rheumatoid arthritis. Rheumatology (Oxford). 2006;45(6):669–75. Epub 2006/03/29View ArticleGoogle Scholar
- Petty RE, Southwood TR, Manners P, Baum J, Glass DN, Goldenberg J, et al. International league of associations for rheumatology classification of juvenile idiopathic arthritis: second revision, Edmonton, 2001. J Rheumatol. 2004;31(2):390–2. Epub 2004/02/05PubMedGoogle Scholar
- Stebulis JA, Rossetti RG, Atez FJ, Zurier RB. Fibroblast-like synovial cells derived from synovial fluid. J Rheumatol. 2005;32(2):301–6. Epub 2005/02/05PubMedGoogle Scholar
- Wu Z, Irizarry RA, Gentleman R, Martinez-Murillo F, Spencer FA. Model-based background adjustment for oligonucleotide expression arrays. J Am Stat Assoc. 2004;99(468):909–17.View ArticleGoogle Scholar
- Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic acids research. 2015;43(7):e47. Epub 2015/01/22View ArticlePubMedPubMed CentralGoogle Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. BioTechniques. 2003;34(2):374–8. Epub 2003/03/05PubMedGoogle Scholar
- Zhang JD RM, Biczok R. ddCt Method for qRT-PCR Data Analysis. http://wwwbioconductororg/packages/37/bioc/vignettes/ddCt/inst/doc/rtPCRpdf. 2017.Google Scholar
- Cox KL, Devanarayan V, Kriauciunas A, Manetta J, Montrose C, Sittampalam S. Immunoassay Methods. In: Sittampalam GS, Coussens NP, Brimacombe K, Grossman A, Arkin M, Auld D, et al., editors. Assay Guidance Manual. Bethesda (MD)2004.Google Scholar
- Nistala K, Adams S, Cambrook H, Ursu S, Olivito B, de Jager W, et al. Th17 plasticity in human autoimmune arthritis is driven by the inflammatory environment. Proc Natl Acad Sci U S A. 2010;107(33):14751–6. Epub 2010/08/04View ArticlePubMedPubMed CentralGoogle Scholar
- Singh K, Colmegna I, He X, Weyand CM, Goronzy JJ. Synoviocyte stimulation by the LFA-1-intercellular adhesion molecule-2-Ezrin-Akt pathway in rheumatoid arthritis. J Immunol. 2008;180(3):1971–8. Epub 2008/01/23View ArticlePubMedGoogle Scholar
- Ferguson IYH, Griffin P, Michel J, Kietz D, de Vallejo A. A novel subset of fibrocyte-like cells in synovial fluid communicate with CD8 T cells and shape the local inflammatory milieu in juvenile idiopathic arthritis (JIA). J Immunol. 2015;194(1 Supplement):121.13.Google Scholar
- Yu S, Nakashima N, Xu BH, Matsuda T, Izumihara A, Sunahara N, et al. Pathological significance of elevated soluble CD14 production in rheumatoid arthritis: in the presence of soluble CD14, lipopolysaccharides at low concentrations activate RA synovial fibroblasts. Rheumatol Int. 1998;17(6):237–43. Epub 1998/05/21View ArticlePubMedGoogle Scholar
- Hunter PJ, Nistala K, Jina N, Eddaoudi A, Thomson W, Hubank M, et al. Biologic predictors of extension of oligoarticular juvenile idiopathic arthritis as determined from synovial fluid cellular composition and gene expression. Arthritis Rheum. 2010;62(3):896–907. Epub 2010/02/04View ArticlePubMedPubMed CentralGoogle Scholar