Skip to main content
  • Research article
  • Open access
  • Published:

Identification of intra-group, inter-individual, and gene-specific variances in mRNA expression profiles in the rheumatoid arthritis synovial membrane

Abstract

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 the Kyoto 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.

Introduction

Human rheumatoid arthritis (RA) is characterized by chronic inflammation and destruction of multiple joints, perpetuated by an abnormally transformed and invasive synovial membrane (SM), forming the so-called pannus tissue [1]. Many activated cell types contribute to the development and progression of RA. Monocytes/macrophages, dendritic cells, T and B cells, endothelial cells, and synovial fibroblasts are major components of the pannus [2–8] and participate in maintaining joint inflammation, degradation of extracellular matrix (ECM) components, and invasion of cartilage and bone [2, 4] as well as fibrosis of the affected joints [9].

The extended analysis of gene expression profiles in RA SM during the last decades has revealed several relevant gene groups affecting development and progression of the disease. Central transcription factors involved as key players in RA pathogenesis are AP-1, NF-κB, Ets-1, and SMADs [10–12]. These factors show binding activity for their cognate recognition sites in the promoters of inflammation-related cytokines (for example, tumor necrosis factor-alpha [TNF-α], interleukin [IL]-1β, and IL-6 [3]) and matrix-degrading enzymes (for example, matrix metalloproteinase [MMP]-1 and MMP-3 [13, 14]). The latter contribute to tissue degradation by destruction of ECM components, including aggrecan or collagen type I-IV, X, and XI [15].

The analysis of those comprehensive expression data has become feasible due to the implementation of microarray-based methods [16]. Therefore, a variety of comparisons can be performed, including differences in gene expression among different groups and/or individuals. In contrast to conventional differential gene expression analyses, the determination of inter-individual gene expression variances, often affecting gene expression of members of the same patient/donor group, is generally not considered in rheumatology, although those variances are known to be a characteristic of many diseases. In trisomy 21, for instance, inter-individual expression variances affect a number of tightly regulated genes. In addition, the variances are independent of the respective level of gene expression, and although only a minority of genes are affected, these genes are thought to be involved in the symptoms of trisomy 21 with the highest phenotypical differences [17]. Significant inter-individual expression variances have also been reported to affect the expression of telomerase subunits in malignant glioma [18] as well as protein tyrosine kinases and phosphatases in human basophils in asthma and inflammatory allergy [19]. The latter implies that such alterations may also play an important role within inflammatory diseases, reflected in either synchronization (that is, a loss of inter-individual gene expression variances) or desynchronization (that is, increased inter-individual gene expression variances) of gene expression within a group of different individuals/patients.

In RA, differences in gene expression profiles for specific genes among two subgroups of RA patients have been reported, but within these subgroups, the differences are limited to distinct expression levels without significant intra-subgroup expression variances [12]. To the best of our knowledge, there are as yet no reports on broad intra-group inter-individual gene expression variations among RA patients.

Interestingly, although the majority of reports show expression variances in tissues from patients with different diseases, variances have also been reported in normal tissues (for example, the human retina [20] or human B-lymphoblastoid cells [21]). In contrast to expression variations in diseases, the variations in normal donors are generally limited to a small number of genes (for example, 2.6% in the human retina [20]). To analyze inter-individual mRNA expression variances in RA, the occurrence of gene-specific expression differences in the SM was analyzed using the Bonferroni/Holm corrected Brown-Forsythe version of the Levene test for variance analysis [22–24] on the basis of genome-wide mRNA expression data in RA (n = 12), osteoarthritis (OA) (n = 10), and normal control (NC) (n = 9) synovial tissue.

Materials and methods

Patients and tissue samples

SM samples were obtained within 10 minutes following tissue excision upon joint replacement/synovectomy from RA (n = 12) and OA (n = 10) patients at the Department of Orthopedics, University Hospital Jena, Waldkrankenhaus 'Rudolf Elle' (Eisenberg, Germany). Tissue samples from joint trauma surgery (n = 9) were used as NCs (Table 1). After removal, tissue samples were frozen and stored at -70°C. Informed patient consent was obtained and the study was approved by the Ethics Committee of University Hospital Jena (Jena, Germany). RA patients were classified according to the American College of Rheumatology criteria [25], OA patients according to the respective criteria for OA [26].

Table 1 Clinical characteristics of the patients at the time of synovectomy/sampling

Isolation of total RNA

Tissue homogenization, total RNA isolation, treatment with RNase-free DNase I (Qiagen, Hilden, Germany), and cDNA synthesis were performed as described previously [27].

Microarray data analysis

RNA probes were labeled according to the instructions of the supplier (Affymetrix, Santa Clara, CA, USA). Analysis of gene expression was carried out using U133A/B oligonucleotide arrays. Hybridization and washing procedures were performed according to the supplier's instructions and microarrays were analyzed by laser scanning (Hewlett-Packard Gene Scanner; Hewlett-Packard Company, Palo Alto, CA, USA). Background-corrected signal intensities were determined using the MAS 5.0 software (Affymetrix). Subsequently, signal intensities were normalized among arrays to facilitate comparisons between different patients. For this purpose, arrays were grouped according to patient/donor groups (RA, n = 12; OA, n = 10; and NC, n = 9). The arrays in each group were normalized using quantile normalization [28]. Original data from microarray analyses were deposited in the Gene Expression Omnibus of the National Center for Biotechnology Information (Bethesda, MD, USA) (accession number GSE12021 [29]).

Real-time reverse transcription-polymerase chain reaction

The data obtained by Affymetrix microarrays were validated for six selected genes (IL13, MAPK8, SMAD2, IL2RG, PLCB1, and ATF5) using real-time reverse transcription-polymerase chain reaction (RT-PCR). PCRs were performed as previously described using a Mastercycler® ep realplex (Eppendorf, Hamburg, Germany) and SYBR-green. To normalize the amount of cDNA in each sample, the expression of the housekeeping gene GAPDH (glyceraldehyde 3-phosphate dehydrogenase) was determined [27]. Product specificity was confirmed by (a) melting curve analysis, (b) agarose gel electrophoresis, and (c) cycle sequencing of the PCR products.

Statistical analysis of gene expression variance

This analysis did not concentrate on differently expressed genes, but on genes with different variances in the three patient groups [30]. The assumption of homogeneity of variance can be rejected by a variance analysis according to Levene [22]. The Brown-Forsythe version of this test was used [23]. For independent groups of data, the null hypothesis (that is, variances are equal) was tested.

To control the stability of the variance, the variance calculation was tested for 2, 3, 5, 7, and 10 samples per group. For fewer than 5 samples, the calculation did not reach stable results, but stable results were achieved for more than 5 patients. In addition, the results of the statistical tests were influenced by the number of samples in each group (that is, small groups did not reach statistical significance).

The P value can be obtained by calculating the value of the cumulative distribution function at the point F. This is equivalent to the integral of the probability density function of the normal distribution over the interval [0, F]. To prevent the accumulation of false-positives due to multiple comparisons, the very strict Bonferroni correction was used [31]. Alternatively, the less conservative Holm correction was applied for the correction of the data [24]. The application of the Holm correction yielded results comparable to those obtained by Bonferroni correction and pointed out only very few new genes.

The variance-fold is defined as the quotient of the variance of one group (for example, OA patients) and the variance of another group (for example, RA patients). If the variance in the second group is higher than 1, the result is the multiplicative inverse and the algebraic sign is inverted. This way, all groups can be compared:

V a r F o l d = { var x ≥ var y : var x / var y var x < var y : − 1 * ( var y / var x ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOvayLaemyyaeMaemOCaiNaemOrayKaem4Ba8MaemiBaWMaemizaqMaeyypa0ZaaiqaaeaafaqaaeGabaaabaGagiODayNaeiyyaeMaeiOCai3aaSbaaSqaaiabdIha4bqabaGccqGHLjYScyGG2bGDcqGGHbqycqGGYbGCdaWgaaWcbaGaemyEaKhabeaakiabcQda6iGbcAha2jabcggaHjabckhaYnaaBaaaleaacqWG4baEaeqaaOGaei4la8IagiODayNaeiyyaeMaeiOCai3aaSbaaSqaaiabdMha5bqabaaakeaacyGG2bGDcqGGHbqycqGGYbGCdaWgaaWcbaGaemiEaGhabeaakiabgYda8iGbcAha2jabcggaHjabckhaYnaaBaaaleaacqWG5bqEaeqaaOGaeiOoaOJaeyOeI0IaeGymaeJaeiOkaOIaeiikaGIagiODayNaeiyyaeMaeiOCai3aaSbaaSqaaiabdMha5bqabaGccqGGVaWlcyGG2bGDcqGGHbqycqGGYbGCdaWgaaWcbaGaemiEaGhabeaakiabcMcaPaaaaiaawUhaaaaa@7339@

The application of a variance filter before testing of the data (excluding variance-fold values between 2.5 and -2.5 from the analysis) yielded equivalent results compared with the initial data analysis including the a posteriori application of the Bonferroni or the Holm correction. Following Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis (see below), the same pathways/complexes were indicated and only the ranking of selected pathways/complexes was changed (for example, the ranking of cytokine–cytokine receptor interactions and the mitogen-activated protein kinase [MAPK] pathway were inverted).

Analysis of inter-individual gene expression variances

Relevant genes were selected using different criteria: (a) a significance level of P ≤ 0.05 (Bonferroni/Holm corrected Brown-Forsythe version of the Levene test) for variance-fold values and (b) a cutoff value for absolute variance-fold levels of greater than 2.5 for higher variances in RA, OA, and NC, respectively. Using these criteria, 568 genes were selected for the comparison between RA and NC (307 with higher variances in RA and 261 with higher variances in NC) while 542 genes were used for the comparison OA versus NC (314 with higher variances in OA and 228 with higher variances in NC). Finally, 333 genes were selected for the comparison between RA and OA (186 with higher variances in RA and 147 with higher variances in OA). All selected genes are presented in Supplementary Table 1 (sorted according to absolute variance-fold values). Inter-individual variances of gene expression among the different groups were analyzed using predefined pathways and functional categories annotated by KEGG [32].

Mapping of probesets onto gene names

Gene names used for KEGG inputs follow the nomenclature of the HUGO Genome Nomenclature Committee [33] and are mostly derived from the Affymetrix annotation feature 'Gene Symbol' for the respective probeset. If required, corresponding RefSeqs were manually inspected.

Statistical KEGG analysis

To ensure that only KEGG pathways with a significant enrichment of more variant genes were obtained for further analyses, the χ2 test statistic was used. Following the calculation of the expected frequency of affected genes in each pathway, the difference between the expected frequency and the absolute frequency was determined. All pathways with a difference of less than 2 were ignored. As a second criterion of the multilevel test, P values of less than or equal to 0.15 were considered statistically significant [34]. Pathways with insignificant P values were examined in detail and subdivided into two or more sub-pathways if possible. In some cases, P values for selected sub-pathways decreased considerably.

Results

Analysis of inter-individual gene expression variances in rheumatoid arthritis, osteoarthritis, and normal control synovial membrane

For the comparison of inter-individual gene expression variances between RA SM (n = 12) and NC SM (n = 9), 568 genes were used (307 with significantly higher variances in RA and 261 with significantly higher variances in NC; P ≤ 0.05, Bonferroni/Holm corrected Brown-Forsythe version of the Levene test), resulting in the identification of 129 affected KEGG pathways/complexes in total (Supplementary Table 1a; shown for IL13 and CXCL13 in Figure 1). These pathways include 10 pathways significantly affected by higher gene expression variances in RA and 6 pathways significantly affected by higher gene expression variances in NC (in both cases P ≤ 0.15, χ2 test).

Figure 1
figure 1

Gene-specific inter-individual gene expression variances. The graph shows the individual gene expression level of rheumatoid arthritis (RA) (n = 12) and osteoarthritis (OA) (n = 10) patients as well as normal control (NC) donors (n = 9) for IL13 and CXCL13 (cytokine–cytokine receptor interactions). The mean gene expression (blue line) and the intra-group inter-individual variances in RA and NC synovial membrane (red bar) are indicated, resulting in significantly enhanced variances among patients within the RA group (P < 0.001, Bonferroni/Holm corrected Brown-Forsythe version of the Levene test).

For the comparison of OA (n = 10) and NC (n = 9) SM, 542 genes were used (314 with significantly higher variances in OA and 228 with significantly higher variances in NC; Supplementary Table 1b). A total of 128 affected KEGG pathways/complexes were identified, including 7 pathways significantly affected by higher gene expression variances in OA and 4 pathways significantly affected by higher gene expression variances in NC.

The comparison of RA (n = 12) and OA (n = 10) SM was performed with 333 genes (186 with significantly higher variances in RA and 147 with significantly higher variances in OA; Supplementary Table 1c). This comparison culminated in the identification of 114 pathways, 3 of which were significantly affected by higher gene expression variances in RA and 4 of which were significantly affected by higher gene expression variances in OA.

Real-time reverse transcription-polymerase chain reaction validation

Validation of the microarray data by real-time RT-PCR was attempted in RA, OA, and NC samples for the genes IL13, MAPK8, SMAD2, IL2RG, PLCB1, and ATF5. In three cases (50%), the results of microarray analyses and real-time RT-PCR were equivalent for RA versus NC (MAPK8: variance-fold 9.8 versus 5.2; IL2RG: variance-fold 5.6 versus 8.9; ATF5: variance-fold 1.7 versus 2.3); in addition, two cases (33%) tended to result in comparable variance-fold values for microarray and real-time RT-PCR (IL13: variance-fold 12 versus 1.3; SMAD2: variance-fold 5 versus 1.1). In only one case (PLCB1; 17%), microarray analyses and real-time RT-PCR validation showed contradictory results (higher variance in NC versus higher variance in RA). For OA versus NC, comparable results were achieved (only IL2RG and ATF5 showed contradictory results).

KEGG pathways identified in the comparison between rheumatoid arthritis and normal control

Pathways significantly affected by inter-individual gene expression variances in rheumatoid arthritis

Ten pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in the comparison between RA and NC, 7 of which were specific for RA, that is, did not appear in the comparison between OA and NC (for example, cytokine–cytokine receptor interactions; Figure 2). The occurrence of gene expression variances in the complete MAPK, transforming growth factor-beta (TGF-β), and apoptosis pathways/complexes did not reach statistical significance. Interestingly, within these pathways, significantly affected sub-pathways/sub-complexes could be identified: the classical TGF-β sub-pathway (Figure 3), the classical and the c-jun kinase (JNK)/p38 MAPK sub-pathway(s) (Figure 4), and the sub-complex of anti-apoptosis (Figure 5). A complete list of significantly affected pathways/complexes is presented in Table 2.

Figure 2
figure 2

Inter-individual mRNA expression variances among cytokine–cytokine receptor interactions in rheumatoid arthritis (RA) compared with normal control (NC). The graph shows genes affected by significant intra-group inter-individual mRNA expression variances in RA compared with NC (P ≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test; labeled in red) among Kyoto Encyclopedia of Genes and Genomes (KEGG) cytokine–cytokine receptor interactions, including the respective sub-pathways (P ≤ 0.15, χ2 test; labeled in red). Cellular processes with potential influence on or relevance for RA pathogenesis (for example, inflammation, proliferation, and cell survival) are labeled in blue, and anti-inflammatory/anti-destructive processes are labeled in black.

Figure 3
figure 3

Inter-individual mRNA expression variances in the transforming growth factor-beta (TGF-β) signaling pathway in rheumatoid arthritis (RA) compared with normal control (NC). The graph shows genes affected by significant intra-group inter-individual mRNA expression variances in RA compared with NC (P ≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test; labeled in red) in the Kyoto Encyclopedia of Genes and Genomes (KEGG) TGF-β signaling pathway. Among the three TGF-β family sub-pathways, the classical TGF-β sub-pathway is significantly affected by gene expression variances (P ≤ 0.15, χ2 test; indicated in red). TGF-β-regulated cellular processes with potential influence on or relevance for RA pathogenesis (for example, angiogenesis and cell survival) are labeled in blue.

Figure 4
figure 4

Inter-individual mRNA expression variances in the mitogen-activated protein kinase (MAPK) signaling pathway in rheumatoid arthritis (RA) compared with normal control (NC). The graph shows genes affected by significant intra-group inter-individual mRNA expression variances in RA compared with NC (P ≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test; labeled in red) in the Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK signaling pathway. Among the three MAPK family sub-pathways, the classical and the c-jun kinase (JNK)/p38 MAPK sub-pathways were significantly affected by gene expression variances (P ≤ 0.15, χ2 test; indicated in red). MAPK-regulated cellular processes with potential influence on or relevance for RA pathogenesis (for example, proliferation, inflammation, and anti-apoptosis) are labeled in blue.

Figure 5
figure 5

Inter-individual mRNA expression variances in the complex of apoptosis in rheumatoid arthritis (RA) compared with normal control (NC). The graph shows genes affected by significant intra-group inter-individual mRNA expression variances in RA compared with NC (P ≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test; labeled in red) in the Kyoto Encyclopedia of Genes and Genomes (KEGG) complex of apoptosis. Among the three apoptosis sub-complexes, the survival factor-dependent sub-complex was significantly affected by gene expression variances (P ≤ 0.15, χ2 test; indicated in red). Cellular processes with potential influence on or relevance for RA pathogenesis (expression of survival genes and cell survival) are labeled in blue.

Table 2 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in rheumatoid arthritis (RA) compared with normal control (that is, higher variances in RA)

Pathways significantly affected by inter-individual gene expression variances in normal control

Six pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in NC compared with RA, including the cell cycle and the Wnt (wingless-type MMTV integration site family) signaling pathway. All pathways/complexes were specific for NC. A complete list of significantly affected pathways/complexes is presented in Table 3.

Table 3 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in normal control (NC) compared with rheumatoid arthritis (that is, higher variances in NC)

KEGG pathways identified in the comparison between osteoarthritis and normal control

Pathways significantly affected by inter-individual gene expression variances in osteoarthritis

Seven pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in OA compared with NC. Among these pathways/complexes, six were specific for OA, including the complexes of apoptosis. A complete list of significantly affected pathways/complexes is presented in Table 4.

Table 4 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in osteoarthritis (OA) compared with normal control (that is, higher variances in OA)

Pathways significantly affected by inter-individual gene expression variances in normal control

Four pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in NC compared with OA. Three of those were specific for NC, including the Toll-like receptor signaling pathway. A complete list of significantly affected pathways/complexes is presented in Table 5.

Table 5 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in normal control (NC) compared with osteoarthritis (that is, higher variances in NC)

KEGG pathways identified in the comparison between rheumatoid arthritis and osteoarthritis

Pathways significantly affected by inter-individual gene expression variances in rheumatoid arthritis

Three pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in RA compared with OA. All pathways/complexes were specific for RA, including the vascular endothelial growth factor (VEGF) and the B-cell receptor signaling pathways. A complete list of significantly affected pathways/complexes is presented in Table 6.

Table 6 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in rheumatoid arthritis (RA) compared with osteoarthritis (that is, higher variances in RA)

Pathways significantly affected by inter-individual gene expression variances in osteoarthritis

Four pathways/complexes significantly affected by inter-individual mRNA expression variances were identified in OA compared with RA (for example, the complex of oxidative phosphorylation). All of them were specific for OA. A complete list of significantly affected pathways/complexes is presented in Table 7.

Table 7 KEGG pathways/complexes significantly affected by intra-group inter-individual gene expression variance in osteoarthritis (OA) compared with rheumatoid arthritis (that is, higher variances in OA)

Discussion

The present microarray-based and real-time RT-PCR-validated, genome-wide mRNA expression analysis in RA, OA, and NC SM by KEGG mapping shows that gene-specific, significant, intra-group/inter-individual variances in gene expression profiles occur in RA. These variances affect a variety of genes involved in numerous pathways/complexes potentially relevant for RA pathogenesis. Since significant variance-fold values are observed for many genes with comparable mean expression levels among different patient/donor groups (data not shown), the manifestation of gene expression variances does not necessarily depend on the respective mean mRNA expression level.

To our knowledge, gene expression variances in RA samples have been reported only for distinct subgroup-specific differences in gene expression profiles of RA patients [12]. Consequently, the present data demonstrate for the first time broad intra-group/inter-individual gene expression variances in RA SM samples, previously observed in other severe diseases such as trisomy 21, malignant glioma, and inflammatory allergy [17–19]. It has been hypothesized that expression variances of regulatory key genes contribute to the individual phenotype of the given disease [17], whether independent of or depending on the expression level.

Currently, the causes for gene expression variances among RA patients are unknown. Possible external reasons may include the higher average age of the individuals in the RA group as well as medication influencing immunological processes and the expression of immunologically relevant genes (for example, methotrexate, prednisolone, sulfasalazine, and/or nonsteroidal anti-inflammatory drugs [35, 36]) or differences in nutrition, with general effects on individual gene expression [37]. The inflammatory status of the respective joint at the time of surgical intervention may also substantially influence gene expression in the RA SM [38]. However, an analysis of the differential gene expression shows that the present RA group is generally characterized by an expression profile highly compatible with previous gene expression studies [39], including the overexpression of several transcription factors (for example, FOS, FOSB, JUN, and STAT1 [10–12]), cytokines/chemokines (for example, IL2, IL4, CCL23, and CCL25 [40]), signal transduction molecules (for example, MAPK9, MAP3K2, PTPN7, and AKT2 [41, 42]), cell cycle regulators (for example, CDC12, CCNB2, and CCNE2 [43]), and heat shock proteins (DNAJ molecules; [44]; data not shown), indicating that the present RA cohort is representative for RA patients in general.

Regarding internal molecular changes in the individuals, a participation of mutations or single nucleotide polymorphisms in different genes is plausible, either directly [45, 46] or via mutated regulators (for example, transcription factors, mRNA stability modifiers, and so on [47]). This also includes broader genomic rearrangements (for example, chromosomal translocations or polysomies [48, 49]) as well as epigenomic modifications (for example, gene/promoter methylation [50]). In addition, the individual composition of cell types in the analyzed SM samples may influence the mRNA expression profile, depending on the inflammatory status and/or cell proliferation, potentially resulting in enhanced immigration/proliferation of T cells, B cells, or synovial fibroblasts [51].

In RA compared with NC, 10 KEGG pathways/complexes are specifically and significantly affected by gene expression variances. As expected, the importance of immunological processes for RA progression [8] is reflected in several pathways directly involved in such networks (Toll-like, T cell, and Fc ε receptor signaling [52–54]). In the SM, alterations in immunological pathways/complexes may contribute to the development of local (and systemic) inflammation, reflecting the highly inflamed status of the joint as one of the major characteristics of RA [2, 55].

RA-specific gene expression variances also occur in cytokine–cytokine receptor interactions. Within this complex, a striking involvement of sub-pathways can be observed, with relevance for chemotaxis (CXC family chemokines [56]), angiogenesis, proliferation, and cell survival (TGF-β family [57, 58]) as well as inflammation, joint destruction, and fibrosis (TNF family [59, 60] and IL2RG shared pathway [9, 61]; Figure 2). Sub-pathways influencing tissue protection (interferon family [62]) or anti-inflammation and anti-angiogenesis (IL13RA1 [interleukin-13 receptor alpha-1] shared pathway [63]) are scarcely affected. Therefore, a specific influence of gene expression variances on cytokine-mediated aspects of the RA can be assumed [64].

Although the following pathways/complexes are not significantly affected by gene expression variances in total, embedded sub-pathways include the majority of affected genes, thus reaching statistical significance. In the TGF-β pathway, only members of the classical TGF-β sub-pathway are significantly affected, thus potentially influencing angiogenesis [58], cell survival [65], and cell proliferation [66] amongst others (Figure 3). Indeed, this (sub-) pathway appears to occupy a central position for the RA pathogenesis, due to the integration of various RA-relevant cellular functions. This is further underlined by its prominent role within the framework of cytokine–cytokine receptor interactions (Figure 2) and its influence on pro-inflammatory/pro-destructive features, either independent of or via MAPK (Figures 3 and 4). Within the MAPK signaling pathway, the 'classical' and the JNK/p38 MAPK sub-pathways – regulating proliferation, anti-apoptosis, and inflammation – are significantly affected by gene expression variances (Figure 4). This may be an indication of a participation of variable gene expression in inflammatory processes via MAPK variants (especially via JNK/MAPK8 [67]) and proliferation of activated cells (for example, synovial fibroblasts and T cells) in RA [68, 69] and MAPK-mediated anti-apoptosis (Figure 4).

Regarding apoptosis, genes particularly involved in the regulation of cell survival and anti-apoptosis are significantly affected by expression variances (Figure 5) [70]. Interestingly, the respective genes in this particular pathway also show increased expression levels in RA SM (data not shown). Pro-apoptotic genes are not affected in this pathway, corresponding to the absence of gene expression variances within the complex of p53-induced apoptosis (data not shown).

Depending on the individual gene expression level in each patient, gene expression variances in regulatory pathways may lead to enhanced inflammation [53, 54], angiogenesis [71, 72], enhanced collagen synthesis and secretion [9], and/or a reduced rate of apoptosis [73], thus potentially contributing to hyperplasia of the SM [74], collagen-dependent fibrosis of the joints [64], and a prolonged life span of activated synovial cells in RA [73, 75].

Since RA and OA samples share many aspects of their respective mRNA expression profiles [76, 77], genes in a number of pathways show comparable variance-fold values in both RA and OA (for example, apoptosis; Tables 2 and 4), thus reflecting basic similarities of joint diseases. However, RA and OA SM samples can be clearly differentiated regarding gene expression variances in other pathways/complexes. In OA, the pathways/complexes affected by higher expression variances than in NC indicate an OA-specific desynchronization of metabolic processes (Table 7). In contrast, RA-specific pathways/complexes are involved in the regulation of VEGF-mediated angiogenesis [74, 75] and vascular permeability [78], as well as B cell-dependent auto-immunity and inflammation [79]. The latter represents the elevated activity status of B cells (including cytokine production and T-cell activation) and – in connection with the affection of the anti-apoptotic sub-pathway – the enhanced survival of self-reactive B cells [5, 6, 80]. This may result in a pronounced role of B cells for disease development in RA compared with OA, which is also reflected in the increasing impact of B cell-directed treatment in RA [81].

In summary, these pathways indicate limited but distinct molecular/cellular differences between RA and OA and demonstrate a major contribution of inflammation and angiogenesis in RA. It is reasonable to assume that the RA pathogenesis is influenced by broad alterations of gene expression in general. For years, only differential gene expression analyses have been performed, resulting in the identification of some key genes but leading to the disregard of several genes with a more limited influence on RA, whose collective influence may still be as large as that of the already-known key players. Therefore, besides ubiquitous elevated expression levels of exceptional pro-inflammatory/pro-destructive key regulators/mediators like TNF-α, IL-1β [82], or MMP-1 [83], elevated or reduced expression levels of many different genes in various pathways/complexes may also influence RA development and progression. In this process, the affection of pathologically relevant pathways with differentially expressed genes may be more important than the character of the respective genes, resulting in different gene expression profiles among individual RA patients as reflected in the gene expression variances of the present study. As a consequence, synchronized or desynchronized gene expression in RA potentially shifts cellular activity from the normal to an activated status.

Regarding diagnosis and therapy of RA, the present results indicate that a more individualized approach for different patients may represent the future of RA treatment. Thus, the determination of individual gene expression patterns may facilitate the selection of the best medication or, more ambitiously, may allow directed modulation of (individually) selected pathways/complexes instead of broad suppression of inflammation by anti-inflammatory/anti-rheumatic drugs [84]. In addition, the present study helped to identify the TGF-β pathway as an accessory key player in RA, due to its central position within the regulatory networks. This suggestion is strongly supported by an emerging number of publications reporting a decisive impact of TGF-β on RA development/progression [57, 58, 85, 86]. The affected pathways (and the respective genes) reported here may provide the basis for further analyses of the RA pathogenesis and the differences between RA and OA on a cellular and molecular level.

Conclusion

In RA, a number of disease-relevant or even disease-specific KEGG pathways/complexes (for example, TGF-β signaling and anti-apoptosis) are characterized by broad intra-group inter-individual expression variances. This indicates that 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. Numerous affected pathways, including TGF-β signaling in a central position, are involved in inflammation, angiogenesis, proliferation, and cell survival, thus potentially influencing characteristic features of RA pathology.

Abbreviations

ECM:

extracellular matrix

IL:

interleukin

IL2RG:

interleukin 2 receptor gamma

JNK:

c-jun kinase

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MAPK:

mitogen-activated protein kinase

MMP:

matrix metalloproteinase

NC:

normal control

OA:

osteoarthritis

PCR:

polymerase chain reaction

RA:

rheumatoid arthritis

RT-PCR:

reverse transcription-polymerase chain reaction

SM:

synovial membrane

TGF-β:

transforming growth factor-beta

TNF:

tumor necrosis factor

VEGF:

vascular endothelial growth factor.

References

  1. Grassi W, De Angelis R, Lamanna G, Cervini C: The clinical features of rheumatoid arthritis. Eur J Radiol. 1998, 27 (Suppl 1): S18-S24. 10.1016/S0720-048X(98)00038-2.

    Article  PubMed  Google Scholar 

  2. Kinne RW, Palombo-Kinne E, Emmrich F: Activation of synovial fibroblasts in rheumatoid arthritis. Ann Rheum Dis. 1995, 54: 501-504. 10.1136/ard.54.6.501-b.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Abeles AM, Pillinger MH: The role of the synovial fibroblast in rheumatoid arthritis: cartilage destruction and the regulation of matrix metalloproteinases. Bull NYU Hosp Jt Dis. 2006, 64: 20-24.

    PubMed  Google Scholar 

  4. Karouzakis E, Neidhart M, Gay RE, Gay S: Molecular and cellular basis of rheumatoid joint destruction. Immunol Lett. 2006, 106: 8-13. 10.1016/j.imlet.2006.04.011.

    Article  CAS  PubMed  Google Scholar 

  5. Weyand CM, Seyler TM, Goronzy JJ: B cells in rheumatoid synovitis. Arthritis Res Ther. 2005, 7 (Suppl 3): S9-12. 10.1186/ar1737.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Keystone E: B cell targeted therapies. Arthritis Res Ther. 2005, 7 (Suppl 3): S13-S18. 10.1186/ar1738.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Firestein GS: Evolving concepts of rheumatoid arthritis. Nature. 2003, 423: 356-361. 10.1038/nature01661.

    Article  CAS  PubMed  Google Scholar 

  8. Firestein GS: Immunologic mechanisms in the pathogenesis of rheumatoid arthritis. J Clin Rheumatol. 2005, 11: S39-S44. 10.1097/01.rhu.0000166673.34461.33.

    Article  PubMed  Google Scholar 

  9. Postlethwaite AE, Holness MA, Katai H, Raghow R: Human fibroblasts synthesize elevated levels of extracellular matrix proteins in response to interleukin 4. J Clin Invest. 1992, 90: 1479-1485. 10.1172/JCI116015.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Firestein GS, Manning AM: Signal transduction and transcription factors in rheumatic disease. Arthritis Rheum. 1999, 42: 609-621. 10.1002/1529-0131(199904)42:4<609::AID-ANR3>3.0.CO;2-I.

    Article  CAS  PubMed  Google Scholar 

  11. Han Z, Boyle DL, Manning AM, Firestein GS: AP-1 and NF-kappaB regulation in rheumatoid arthritis and murine collagen-induced arthritis. Autoimmunity. 1998, 28: 197-208. 10.3109/08916939808995367.

    Article  CAS  PubMed  Google Scholar 

  12. Pouw Kraan van der TC, van Gaalen FA, Kasperkovitz PV, Verbeet NL, Smeets TJ, Kraan MC, Fero M, Tak PP, Huizinga TW, Pieterman E, Breedveld FC, Alizadeh AA, Verweij CL: 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. 10.1002/art.11096.

    Article  PubMed  Google Scholar 

  13. Beer S, Oleszewski M, Gutwein P, Geiger C, Altevogt P: Metalloproteinase-mediated release of the ectodomain of L1 adhesion molecule. J Cell Sci. 1999, 112 (Pt 16): 2667-2675.

    CAS  PubMed  Google Scholar 

  14. McCachren SS: Expression of metalloproteinases and metalloproteinase inhibitor in human arthritic synovium. Arthritis Rheum. 1991, 34: 1085-1093. 10.1002/art.1780340904.

    Article  CAS  PubMed  Google Scholar 

  15. Wu JJ, Lark MW, Chun LE, Eyre DR: Sites of stromelysin cleavage in collagen types II, IX, X, and XI of cartilage. J Biol Chem. 1991, 266: 5625-5628.

    CAS  PubMed  Google Scholar 

  16. Hardiman G: Microarrays Technologies 2006: an overview. Pharmacogenomics. 2006, 7: 1153-1158. 10.2217/14622416.7.8.1153.

    Article  PubMed  Google Scholar 

  17. Sultan M, Piccini I, Balzereit D, Herwig R, Saran NG, Lehrach H, Reeves RH, Yaspo ML: Gene expression variation in Down's syndrome mice allows prioritization of candidate genes. Genome Biol. 2007, 8: R91-10.1186/gb-2007-8-5-r91.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Shervington A, Patel R, Lu C, Cruickshanks N, Lea R, Roberts G, Dawson T, Shervington L: Telomerase subunits expression variation between biopsy samples and cell lines derived from malignant glioma. Brain Res. 2007, 1134: 45-52. 10.1016/j.brainres.2006.11.093.

    Article  CAS  PubMed  Google Scholar 

  19. MacGlashan DW: Relationship between spleen tyrosine kinase and phosphatidylinositol 5' phosphatase expression and secretion from human basophils in the general population. J Allergy Clin Immunol. 2007, 119: 626-633. 10.1016/j.jaci.2006.09.040.

    Article  CAS  PubMed  Google Scholar 

  20. Chowers I, Liu D, Farkas RH, Gunatilaka TL, Hackam AS, Bernstein SL, Campochiaro PA, Parmigiani G, Zack DJ: Gene expression variation in the adult human retina. Hum Mol Genet. 2003, 12: 2881-2893. 10.1093/hmg/ddg326.

    Article  CAS  PubMed  Google Scholar 

  21. Storey JD, Madeoy J, Strout JL, Wurfel M, Ronald J, Akey JM: Gene-expression variation within and among human populations. Am J Hum Genet. 2007, 80: 502-509. 10.1086/512017.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Levene H: Robust tests for equality of variances. Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling. Edited by: Olkin I, Ghurye SG, Hoeffding W, Madow WG, Mann HB. 1960, Palo Alto, CA: Stanford University Press, 278-292.

    Google Scholar 

  23. Brown MB, Forsythe AB: Robust tests for equality of variances. J Amer Statist Assoc. 1974, 69: 364-367. 10.2307/2285659.

    Article  Google Scholar 

  24. Holm S: A simple sequentially rejective multiple test procedure. Scand J Statist. 1979, 6: 65-70.

    Google Scholar 

  25. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.

    Article  CAS  PubMed  Google Scholar 

  26. Altman R, Asch E, Bloch D, Bole G, Borenstein D, Brandt K, Christy W, Cooke TD, Greenwald R, Hochberg M, Howell D, Kaplan D, Koopman W, Longley S, Mankin H, McShane DJ, Medsger T, Meenan R, Mikkelsen W, Moskowitz R, Murphy W, Rothschild B, Segal M, Sokoloff L, Wolfe F: Development of criteria for the classification and reporting of osteoarthritis. Classification of osteoarthritis of the knee. Diagnostic and Therapeutic Criteria Committee of the American Rheumatism Association. Arthritis Rheum. 1986, 29: 1039-1049. 10.1002/art.1780290816.

    Article  CAS  PubMed  Google Scholar 

  27. Huber R, Kunisch E, Gluck B, Egerer R, Sickinger S, Kinne RW: Comparison of conventional and real-time RT-PCR for the quantitation of jun protooncogene mRNA and analysis of junB mRNA expression in synovial membranes and isolated synovial fibroblasts from rheumatoid arthritis patients. Z Rheumatol. 2003, 62: 378-389. 10.1007/s00393-003-0472-4.

    Article  CAS  PubMed  Google Scholar 

  28. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.

    Article  CAS  PubMed  Google Scholar 

  29. Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo]

  30. Fisher RA: The correlation between relatives on the supposition of mendelian inheritance. Trans Roy Soc Edinb. 1918, 52: 399-433.

    Article  Google Scholar 

  31. Abdi H: Bonferroni and Sidak corrections for multiple comparisons. Encyclopedia of Measurement and Statistics. Edited by: Salkind NJ. 2007, Thousand Oaks, CA: SAGE Publications

    Google Scholar 

  32. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27: 29-34. 10.1093/nar/27.1.29.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. HUGO Gene Nomenclature Committee. [http://www.genenames.org]

  34. Edwards BJ, Haynes C, Levenstien MA, Finch SJ, Gordon D: Power and sample size calculations in the presence of phenotype errors for case/control genetic association studies. BMC Genet. 2005, 6: 18-10.1186/1471-2156-6-18.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Haupl T, Yahyawi M, Lubke C, Ringe J, Rohrlach T, Burmester GR, Sittinger M, Kaps C: Gene expression profiling of rheumatoid arthritis synovial cells treated with antirheumatic drugs. J Biomol Screen. 2007, 12: 328-340. 10.1177/1087057107299261.

    Article  PubMed  Google Scholar 

  36. Ospelt C, Gay S: Antirheumatic drugs and gene signatures. Curr Opin Investig Drugs. 2007, 8: 385-389.

    CAS  PubMed  Google Scholar 

  37. Paoloni-Giacobino A, Grimble R, Pichard C: Genetics and nutrition. Clin Nutr. 2003, 22: 429-435. 10.1016/S0261-5614(03)00064-5.

    Article  CAS  PubMed  Google Scholar 

  38. Hahn G, Stuhlmuller B, Hain N, Kalden JR, Pfizenmaier K, Burmester GR: Modulation of monocyte activation in patients with rheumatoid arthritis by leukapheresis therapy. J Clin Invest. 1993, 91: 862-870. 10.1172/JCI116307.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Batliwalla FM, Baechler EC, Xiao X, Li W, Balasubramanian S, Khalili H, Damle A, Ortmann WA, Perrone A, Kantor AB, Gulko PS, Kern M, Furie R, Behrens TW, Gregersen PK: Peripheral blood gene expression profiling in rheumatoid arthritis. Genes Immun. 2005, 6: 388-397. 10.1038/sj.gene.6364209.

    Article  CAS  PubMed  Google Scholar 

  40. Sweeney SE, Firestein GS: Rheumatoid arthritis: regulation of synovial inflammation. Int J Biochem Cell Biol. 2004, 36: 372-378. 10.1016/S1357-2725(03)00259-0.

    Article  CAS  PubMed  Google Scholar 

  41. Hammaker DR, Boyle DL, Chabaud-Riou M, Firestein GS: Regulation of c-Jun N-terminal kinase by MEKK-2 and mitogen-activated protein kinase kinase kinases in rheumatoid arthritis. J Immunol. 2004, 172: 1612-1618.

    Article  CAS  PubMed  Google Scholar 

  42. Liagre B, Vergne-Salle P, Leger DY, Beneytout JL: Inhibition of human rheumatoid arthritis synovial cell survival by hecogenin and tigogenin is associated with increased apoptosis, p38 mitogen-activated protein kinase activity and upregulation of cyclooxygenase-2. Int J Mol Med. 2007, 20: 451-460.

    CAS  PubMed  Google Scholar 

  43. Taranto E, Leech M: Expression and function of cell cycle proteins in rheumatoid arthritis synovial tissue. Histol Histopathol. 2006, 21: 205-211.

    CAS  PubMed  Google Scholar 

  44. Kurzik-Dumke U, Schick C, Rzepka R, Melchers I: Overexpression of human homologs of the bacterial DnaJ chaperone in the synovial tissue of patients with rheumatoid arthritis. Arthritis Rheum. 1999, 42: 210-220. 10.1002/1529-0131(199902)42:2<210::AID-ANR2>3.0.CO;2-U.

    Article  CAS  PubMed  Google Scholar 

  45. Kim SY, Han SW, Kim GW, Lee JM, Kang YM: TGF-beta1 polymorphism determines the progression of joint damage in rheumatoid arthritis. Scand J Rheumatol. 2004, 33: 389-394. 10.1080/03009740410010344.

    Article  CAS  PubMed  Google Scholar 

  46. Han SW, Kim GW, Seo JS, Kim SJ, Sa KH, Park JY, Lee J, Kim SY, Goronzy JJ, Weyand CM, Kang YM: VEGF gene polymorphisms and susceptibility to rheumatoid arthritis. Rheumatology (Oxford). 2004, 43: 1173-1177. 10.1093/rheumatology/keh281.

    Article  CAS  Google Scholar 

  47. Martinez A, Valdivia A, Pascual-Salcedo D, Balsa A, Fernandez-Gutierrez B, De la CE, Urcelay E: Role of SLC22A4, SLC22A5, and RUNX1 genes in rheumatoid arthritis. J Rheumatol. 2006, 33: 842-846.

    CAS  PubMed  Google Scholar 

  48. Kinne RW, Liehr T, Beensen V, Kunisch E, Zimmermann T, Holland H, Pfeiffer R, Stahl HD, Lungershausen W, Hein G, Roth A, Emmrich F, Claussen U, Froster UG: Mosaic chromosomal aberrations in synovial fibroblasts of patients with rheumatoid arthritis, osteoarthritis, and other inflammatory joint diseases. Arthritis Res. 2001, 3: 319-330. 10.1186/ar322.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Kinne RW, Kunisch E, Beensen V, Zimmermann T, Emmrich F, Petrow P, Lungershausen W, Hein G, Braun RK, Foerster M, Kroegel C, Winter R, Liesaus E, Fuhrmann RA, Roth A, Claussen U, Liehr T: Synovial fibroblasts and synovial macrophages from patients with rheumatoid arthritis and other inflammatory joint diseases show chromosomal aberrations. Genes Chromosomes Cancer. 2003, 38: 53-67. 10.1002/gcc.10242.

    Article  PubMed  Google Scholar 

  50. Shin HJ, Park HY, Jeong SJ, Park HW, Kim YK, Cho SH, Kim YY, Cho ML, Kim HY, Min KU, Lee CW: STAT4 expression in human T cells is regulated by DNA methylation but not by promoter polymorphism. J Immunol. 2005, 175: 7143-7150.

    Article  CAS  PubMed  Google Scholar 

  51. Hoffmann M, Pohlers D, Koczan D, Thiesen HJ, Wolfl S, Kinne RW: Robust computational reconstitution – a new method for the comparative analysis of gene expression in tissues and isolated cell fractions. BMC Bioinformatics. 2006, 7: 369-10.1186/1471-2105-7-369.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Andreakos E, Sacre S, Foxwell BM, Feldmann M: The toll-like receptor-nuclear factor kappaB pathway in rheumatoid arthritis. Front Biosci. 2005, 10: 2478-2488. 10.2741/1712.

    Article  CAS  PubMed  Google Scholar 

  53. Zhang Z, Gorman C, Clark JM, Cope AP: Rheumatoid arthritis: a disease of chronic, low-amplitude signals transduced through T cell antigen receptors?. Wien Med Wochenschr. 2006, 156: 2-10. 10.1007/s10354-005-0246-5.

    Article  PubMed  Google Scholar 

  54. Takai T: Fc receptors and their role in immune regulation and autoimmunity. J Clin Immunol. 2005, 25: 1-18. 10.1007/s10875-005-0353-8.

    Article  CAS  PubMed  Google Scholar 

  55. Huber LC, Distler O, Tarner I, Gay RE, Gay S, Pap T: Synovial fibroblasts: key players in rheumatoid arthritis. Rheumatology (Oxford). 2006, 45: 669-675. 10.1093/rheumatology/kel065.

    Article  CAS  Google Scholar 

  56. Pierer M, Rethage J, Seibl R, Lauener R, Brentano F, Wagner U, Hantzschel H, Michel BA, Gay RE, Gay S, Kyburz D: Chemokine secretion of rheumatoid arthritis synovial fibroblasts stimulated by Toll-like receptor 2 ligands. J Immunol. 2004, 172: 1256-1265.

    Article  CAS  PubMed  Google Scholar 

  57. Pohlers D, Beyer A, Koczan D, Wilhelm T, Thiesen HJ, Kinne RW: Constitutive upregulation of the TGF-b pathway in rheumatoid arthritis synovial fibroblasts. Arthritis Res Ther. 2007, 9: R59-10.1186/ar2217.

    Article  PubMed Central  PubMed  Google Scholar 

  58. Maruotti N, Cantatore FP, Crivellato E, Vacca A, Ribatti D: Angiogenesis in rheumatoid arthritis. Histol Histopathol. 2006, 21: 557-566.

    CAS  PubMed  Google Scholar 

  59. Kollias G, Douni E, Kassiotis G, Kontoyiannis D: The function of tumour necrosis factor and receptors in models of multi-organ inflammation, rheumatoid arthritis, multiple sclerosis and inflammatory bowel disease. Ann Rheum Dis. 1999, 58 (Suppl 1): I32-I39.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Rannou F, Francois M, Corvol MT, Berenbaum F: Cartilage breakdown in rheumatoid arthritis. Joint Bone Spine. 2006, 73: 29-36. 10.1016/j.jbspin.2004.12.013.

    Article  CAS  PubMed  Google Scholar 

  61. Young DA, Hegen M, Ma HL, Whitters MJ, Albert LM, Lowe L, Senices M, Wu PW, Sibley B, Leathurby Y, Brown TP, Nickerson-Nutter C, Keith JC, Collins M: Blockade of the interleukin-21/interleukin-21 receptor pathway ameliorates disease in animal models of rheumatoid arthritis. Arthritis Rheum. 2007, 56: 1152-1163. 10.1002/art.22452.

    Article  CAS  PubMed  Google Scholar 

  62. Tak PP: IFN-beta in rheumatoid arthritis. Front Biosci. 2004, 9: 3242-3247. 10.2741/1475.

    Article  CAS  PubMed  Google Scholar 

  63. Haas CS, Amin MA, Ruth JH, Allen BL, Ahmed S, Pakozdi A, Woods JM, Shahrara S, Koch AE: In vivo inhibition of angiogenesis by interleukin-13 gene therapy in a rat model of rheumatoid arthritis. Arthritis Rheum. 2007, 56: 2535-2548. 10.1002/art.22823.

    Article  CAS  PubMed  Google Scholar 

  64. Szekanecz Z, Koch AE: Update on synovitis. Curr Rheumatol Rep. 2001, 3: 53-63. 10.1007/s11926-001-0051-0.

    Article  CAS  PubMed  Google Scholar 

  65. Kawakami A, Urayama S, Yamasaki S, Hida A, Miyashita T, Kamachi M, Nakashima K, Tanaka F, Ida H, Kawabe Y, Aoyagi T, Furuichi I, Migita K, Origuchi T, Eguchi K: Anti-apoptogenic function of TGFbeta1 for human synovial cells: TGFbeta1 protects cultured synovial cells from mitochondrial perturbation induced by several apoptogenic stimuli. Ann Rheum Dis. 2004, 63: 95-97. 10.1136/ard.2003.014159.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Bira Y, Tani K, Nishioka Y, Miyata J, Sato K, Hayashi A, Nakaya Y, Sone S: Transforming growth factor beta stimulates rheumatoid synovial fibroblasts via the type II receptor. Mod Rheumatol. 2005, 15: 108-113. 10.1007/s10165-004-0378-2.

    Article  CAS  PubMed  Google Scholar 

  67. Han Z, Boyle DL, Aupperle KR, Bennett B, Manning AM, Firestein GS: Jun N-terminal kinase in rheumatoid arthritis. J Pharmacol Exp Ther. 1999, 291: 124-130.

    CAS  PubMed  Google Scholar 

  68. Ospelt C, Neidhart M, Gay RE, Gay S: Synovial activation in rheumatoid arthritis. Front Biosci. 2004, 9: 2323-2334. 10.2741/1399.

    Article  CAS  PubMed  Google Scholar 

  69. Forre O, Waalen K, Natvig JB, Kjeldsen-Kragh J: Evidence for activation of rheumatoid synovial T lymphocytes – development of rheumatoid T cell clones. Scand J Rheumatol Suppl. 1988, 76: 153-160. 10.3109/03009748809102965.

    Article  CAS  PubMed  Google Scholar 

  70. Vermeulen K, Berneman ZN, Van Bockstaele DR: Cell cycle and apoptosis. Cell Prolif. 2003, 36: 165-175. 10.1046/j.1365-2184.2003.00267.x.

    Article  CAS  PubMed  Google Scholar 

  71. Hayes AJ: Angioneogenesis in rheumatoid arthritis. Lancet. 1999, 354: 423-424. 10.1016/S0140-6736(05)75838-4.

    Article  CAS  PubMed  Google Scholar 

  72. Hirohata S, Sakakibara J: Angioneogenesis as a possible elusive triggering factor in rheumatoid arthritis. Lancet. 1999, 353: 1331-10.1016/S0140-6736(98)05912-1.

    Article  CAS  PubMed  Google Scholar 

  73. Gaur U, Aggarwal BB: Regulation of proliferation, survival and apoptosis by members of the TNF superfamily. Biochem Pharmacol. 2003, 66: 1403-1408. 10.1016/S0006-2952(03)00490-8.

    Article  CAS  PubMed  Google Scholar 

  74. Malemud CJ: Growth hormone, VEGF and FGF: involvement in rheumatoid arthritis. Clin Chim Acta. 2007, 375: 10-19. 10.1016/j.cca.2006.06.033.

    Article  CAS  PubMed  Google Scholar 

  75. Byrne AM, Bouchier-Hayes DJ, Harmey JH: Angiogenic and cell survival functions of vascular endothelial growth factor (VEGF). J Cell Mol Med. 2005, 9: 777-794. 10.1111/j.1582-4934.2005.tb00379.x.

    Article  CAS  PubMed  Google Scholar 

  76. Firestein GS, Alvaro-Gracia JM, Maki R: Quantitative analysis of cytokine gene expression in rheumatoid arthritis. J Immunol. 1990, 144: 3347-3353.

    CAS  PubMed  Google Scholar 

  77. Nakamura Y, Nawata M, Wakitani S: Expression profiles and functional analyses of Wnt-related genes in human joint disorders. Am J Pathol. 2005, 167: 97-105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  78. Middleton J, Americh L, Gayon R, Julien D, Aguilar L, Amalric F, Girard JP: Endothelial cell phenotypes in the rheumatoid synovium: activated, angiogenic, apoptotic and leaky. Arthritis Res Ther. 2004, 6: 60-72. 10.1186/ar1156.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  79. Amu S, Stromberg K, Bokarewa M, Tarkowski A, Brisslert M: CD25-expressing B-lymphocytes in rheumatic diseases. Scand J Immunol. 2007, 65: 182-191. 10.1111/j.1365-3083.2006.01889.x.

    Article  CAS  PubMed  Google Scholar 

  80. Szodoray P, Alex P, Frank MB, Turner M, Turner S, Knowlton N, Cadwell C, Dozmorov I, Tang Y, Wilson PC, Jonsson R, Centola M: A genome-scale assessment of peripheral blood B-cell molecular homeostasis in patients with rheumatoid arthritis. Rheumatology (Oxford). 2006, 45: 1466-1476. 10.1093/rheumatology/kel095.

    Article  CAS  Google Scholar 

  81. Anolik JH, Ravikumar R, Barnard J, Owen T, Almudevar A, Milner EC, Miller CH, Dutcher PO, Hadley JA, Sanz I: Cutting edge: anti-tumor necrosis factor therapy in rheumatoid arthritis inhibits memory B lymphocytes via effects on lymphoid germinal centers and follicular dendritic cell networks. J Immunol. 2008, 180: 688-692.

    Article  CAS  PubMed  Google Scholar 

  82. Dayer JM: Interleukin 1 or tumor necrosis factor-alpha: which is the real target in rheumatoid arthritis?. J Rheumatol Suppl. 2002, 65: 10-15.

    CAS  PubMed  Google Scholar 

  83. Pardo A, Selman M: MMP-1: the elder of the family. Int J Biochem Cell Biol. 2005, 37: 283-288. 10.1016/j.biocel.2004.06.017.

    Article  CAS  PubMed  Google Scholar 

  84. Smolen J, Aletaha D: The burden of rheumatoid arthritis and access to treatment: a medical overview. Eur J Health Econ. 2008, 8 (Suppl 2): S39-S47. 10.1007/s10198-007-0087-9.

    Article  PubMed  Google Scholar 

  85. Hammaker DR, Boyle DL, Inoue T, Firestein GS: Regulation of the JNK pathway by TGF-beta activated kinase 1 in rheumatoid arthritis synoviocytes. Arthritis Res Ther. 2007, 9: R57-10.1186/ar2215.

    Article  PubMed Central  PubMed  Google Scholar 

  86. Szekanecz Z, Haines GK, Harlow LA, Shah MR, Fong TW, Fu R, Lin SJ, Rayan G, Koch AE: Increased synovial expression of transforming growth factor (TGF)-beta receptor endoglin and TGF-beta 1 in rheumatoid arthritis: possible interactions in the pathogenesis of the disease. Clin Immunol Immunopathol. 1995, 76: 187-194. 10.1006/clin.1995.1114.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Ernesta Palombo-Kinne for critical reading of the manuscript and Bärbel Ukena, Ulrike Körner, and Ildiko Toth for excellent technical assistance. We are grateful to Andreas Roth, Rando Winter, Renée Fuhrmann, and Rudolf-Albrecht Venbrocks (Department of Orthopedics, University Hospital Jena, Eisenberg, Germany) as well as Wolfgang Lungershausen (Department of Traumatology, University Hospital Jena, Jena, Germany) for providing patient/donor material. The study was supported by the German Federal Ministry of Education and Research (BMBF) (grant FKZ 010405 to RWK), the Interdisciplinary Center for Clinical Research (IZKF) Jena (grant FKZ 0313652A to RG), and the Jena Centre for Bioinformatics (grant FKZ 0313652B and grant 01GS0413, NGFN-2 to RWK). RH was supported by a grant from the German National Academic Foundation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Raimund W Kinne.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RH performed the KEGG analyses, contributed to the real-time RT-PCR analyses, and participated in the writing of the manuscript. CH analyzed the microarray data, performed the bioinformatic analyses, and participated in the writing of the manuscript. RH and CH contributed equally to this work. UG participated in the data analyses. DP performed the real-time RT-PCR analyses. DK performed the Affymetrix microarray experiments. RG participated in the design and coordination of the study, including supervision of the bioinformatic analyses. RWK contributed to the design and coordination of the study and participated in the writing of the manuscript. All authors read and approved the final version of the manuscript.

Electronic supplementary material

13075_2007_2334_MOESM1_ESM.doc

Additional file 1: 'Supplementary Table 1A: Genes affected by intra-group, inter-individual mRNA expression variances (RA compared to NC)', 'Supplementary Table 1B: Genes affected by intra-group, inter-individual mRNA expression variances (OA compared to NC)', 'Supplementary Table 1C: Genes affected by intra-group, inter-individual mRNA expression variances (RA compared to OA)'. For KEGG analyses, relevant genes were selected according to (i) a significance level of p ≤ 0.05 (Bonferroni/Holm corrected Brown-Forsythe version of the Levene test) for variance-fold values and (ii) a cutoff value for absolute variance-fold levels of > 2.5 for higher variances in RA, OA, and NC, respectively. (A) 568 genes were selected for the comparison between RA and NC (307 with higher variances in RA, 261 with higher variances in NC), (B) 542 genes were used for the comparison OA versus NC (314 with higher variances in OA, 228 with higher variances in NC), and (C) 333 genes were selected for the comparison between RA and OA (186 with higher variances in RA, 147 with higher variances in OA). All genes are sorted according to absolute variance-fold values. (DOC 2 MB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Huber, R., Hummert, C., Gausmann, U. et al. Identification of intra-group, inter-individual, and gene-specific variances in mRNA expression profiles in the rheumatoid arthritis synovial membrane. Arthritis Res Ther 10, R98 (2008). https://doi.org/10.1186/ar2485

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/ar2485

Keywords