Introduction: Discrimination of rheumatoid arthritis (RA) patients from patients with other inflammatory or degenerative joint diseases or healthy individuals purely on the basis of genes differentially expressed in high-throughput data has proven very difficult. Thus, the present study sought to achieve such discrimination by employing a novel unbiased approach using rule-based classifiers. Methods: Three multi-center genome-wide transcriptomic data sets (Affymetrix HG-U133 A/B) from a total of 79 individuals, including 20 healthy controls (control group - CG), as well as 26 osteoarthritis (OA) and 33 RA patients, were used to infer rule-based classifiers to discriminate the disease groups. The rules were ranked with respect to Kiendl’s statistical relevance index, and the resulting rule set was optimized by pruning. The rule sets were inferred separately from data of one of three centers and applied to the two remaining centers for validation. All rules from the optimized rule sets of all centers were used to analyze their biological relevance applying the software Pathway Studio. Results: The optimized rule sets for the three centers contained a total of 29, 20, and 8 rules (including 10, 8, and 4 rules for ‘RA’), respectively. The mean sensitivity for the prediction of RA based on six center-to-center tests was 96% (range 90% to 100%), that for OA 86% (range 40% to 100%). The mean specificity for RA prediction was 94% (range 80% to 100%), that for OA 96% (range 83.3% to 100%). The average overall accuracy of the three different rule-based classifiers was 91% (range 80% to 100%). Unbiased analyses by Pathway Studio of the gene sets obtained by discrimination of RA from OA and CG with rule-based classifiers resulted in the identification of the pathogenetically and/or therapeutically relevant interferon-gamma and GM-CSF pathways. Conclusion: First-time application of rule-based classifiers for the discrimination of RA resulted in high performance, with means for all assessment parameters close to or higher than 90%. In addition, this unbiased, new approach resulted in the identification not only of pathways known to be critical to RA, but also of novel molecules such as serine/threonine kinase 10.
Introduction: Rheumatoid arthritis (RA) is a chronic inflammatory and destructive joint disease characterized by overexpression of pro-inflammatory/pro-destructive genes and other activating genes (for example, proto-oncogenes) in the synovial membrane (SM). The gene expression in disease is often characterized by significant inter-individual variances via specific synchronization/desynchronization of gene expression. To elucidate the contribution of the variance to the pathogenesis of disease, expression variances were tested in SM samples of RA patients, osteoarthritis (OA) patients, and normal controls (NCs).Method: Analysis of gene expression in RA, OA, and NC samples was carried out using Affymetrix U133A/B oligonucleotide arrays, and the results were validated by real-time reverse transcription-polymerase chain reaction. For the comparison between RA and NC, 568 genes with significantly different variances in the two groups (P≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test) were selected. For the comparison between RA and OA, 333 genes were selected. By means of theKyoto Encyclopedia of Genes and Genomes, the pathways/complexes significantly affected by higher gene expression variances were identified in each group.Results: Ten pathways/complexes significantly affected by higher gene expression variances were identified in RA compared with NC, including cytokine–cytokine receptor interactions, the transforming growth factor-beta pathway, and anti-apoptosis. Compared with OA, three pathways with significantly higher variances were identified in RA (for example, B-cell receptor signaling and vascular endothelial growth factor signaling). Functionally, the majority of the identified pathways are involved in the regulation of inflammation, proliferation, cell survival, and angiogenesis.Conclusion: In RA, a number of disease-relevant or even disease-specific pathways/complexes are characterized by broad intra-group inter-individual expression variances. Thus, RA pathogenesis in different individuals may depend to a lesser extent on common alterations of the expression of specific key genes, and rather on individual-specific alterations of different genes resulting in common disturbances of key pathways.
Genome-wide gene expression was comparatively investigated in early-passage rheumatoid arthritis (RA) and osteoarthritis (OA) synovial fibroblasts (SFBs; n = 6 each) using oligonucleotide microarrays; mRNA/protein data were validated by quantitative PCR (qPCR) and western blotting and immunohistochemistry, respectively. Gene set enrichment analysis (GSEA) of the microarray data suggested constitutive upregulation of components of the transforming growth factor (TGF)-β pathway in RA SFBs, with 2 hits in the top 30 regulated pathways. The growth factor TGF-β1, its receptor TGFBR1, the TGF-β binding proteins LTBP1/2, the TGF-β-releasing thrombospondin 1 (THBS1), the negative effector SkiL, and the smad-associated molecule SARA were upregulated in RA SFBs compared to OA SFBs, whereas TGF-β2 was downregulated. Upregulation of TGF-β1 and THBS1 mRNA (both positively correlated with clinical markers of disease activity/severity) and downregulation of TGF-β2 mRNA in RA SFBs were confirmed by qPCR. TGFBR1 mRNA (only numerically upregulated in RA SFBs) and SkiL mRNA were not differentially expressed. At the protein level, TGF-β1 showed a slightly higher expression, and the signal-transducing TGFBR1 and the TGF-β-activating THBS1 a significantly higher expression in RA SFBs than in OA SFBs. Consistent with the upregulated TGF-β pathway in RA SFBs, stimulation with TGF-β1 resulted in a significantly enhanced expression of matrix-metalloproteinase (MMP)-11 mRNA and protein in RA SFBs, but not in OA SFBs. In conclusion, RA SFBs show broad, constitutive alterations of the TGF-β pathway. The abundance of TGF-β, in conjunction with an augmented mRNA and/or protein expression of TGF-β-releasing THBS1 and TGFBR1, suggests a pathogenetic role of TGF-β-induced effects on SFBs in RA, for example, the augmentation of MMP-mediated matrix degradation/remodeling.
MicroRNAs (miRNAs) are small non-coding RNA molecules acting as post-transcriptional regulators of gene expression. They are involved in many biological processes, and their dysregulation is implicated in various diseases, including multiple sclerosis (MS). Interferon-beta (IFN-beta) is widely used as a first-line immunomodulatory treatment of MS patients. Here, we present the first longitudinal study on the miRNA expression changes in response to IFN-beta therapy. Peripheral blood mononuclear cells (PBMC) were obtained before treatment initiation as well as after two days, four days, and one month, from patients with clinically isolated syndrome (CIS) and patients with relapsing-remitting MS (RRMS). We measured the expression of 651 mature miRNAs and about 19,000 mRNAs in parallel using real-time PCR arrays and Affymetrix microarrays. We observed that the up-regulation of IFN-beta-responsive genes is accompanied by a down-regulation of several miRNAs, including members of the mir-29 family. These differentially expressed miRNAs were found to be associated with apoptotic processes and IFN feedback loops. A network of miRNA-mRNA target interactions was constructed by integrating the information from different databases. Our results suggest that miRNA-mediated regulation plays an important role in the mechanisms of action of IFN-beta, not only in the treatment of MS but also in normal immune responses. miRNA expression levels in the blood may serve as a biomarker of the biological effects of IFN-beta therapy that may predict individual disease activity and progression.
Introduction: About 30% of rheumatoid arthritis patients fail to respond adequately to TNFα-blocking therapy. There is a medical and socioeconomic need to identify molecular markers for an early prediction of responders and nonresponders. Methods: RNA was extracted from peripheral blood mononuclear cells of 19 rheumatoid arthritis patients before the first application of the TNFα blocker etanercept as well as after 72 hours. Clinical response was assessed over 3 months using the 28-joint-count Disease Activity Score and X-ray scans. Supervised learning methods were applied to Affymetrix Human Genome U133 microarray data analysis to determine highly selective discriminatory gene pairs or triplets with prognostic relevance for the clinical outcome evinced by a decline of the 28-joint-count Disease Activity Score by 1.2. Results: Early downregulation of expression levels secondary to TNFα neutralization was associated with good clinical responses, as shown by a decline in overall disease activity 3 months after the start of treatment. Informative gene sets include genes (for example, NFKBIA, CCL4, IL8, IL1B, TNFAIP3, PDE4B, PPP1R15A and ADM) involved in different pathways and cellular processes such as TNFα signalling via NFκB, NFκB-independent signalling via cAMP, and the regulation of cellular and oxidative stress response. Pairs and triplets within these genes were found to have a high prognostic value, reflected by prediction accuracies of over 89% for seven selected gene pairs and of 95% for 10 specific gene triplets. Conclusion: Our data underline that early gene expression profiling is instrumental in identifying candidate biomarkers to predict therapeutic outcomes of anti-TNFα treatment regimes.
A subset of patients with multiple sclerosis (MS) shows an increased endogenous IFN-like activity before initiation of IFN-beta treatment. The molecular basis of this phenomenon and its relevance to predict individual therapy outcomes are not yet fully understood. We studied the expression patterns of these patients, the prognostic value of an elevated IFN-like activity, and the gene regulatory effects of exogenously administered IFN-beta.