Interaction between rheumatoid arthritis and pregnancy: correlation of molecular data with clinical disease activity measures
1Department of Rheumatology and Clinical Immunology, Charité-University Medicine, Berlin, Germany, 2Department of Rheumatology and Clinical Immunology/Allergology, University Hospital of Bern, Bern, Switzerland and 3German Arthritis Research Center (DRFZ), Berlin, Germany.
Correspondence to: T. Häupl, Department of Rheumatology and Clinical Immunology, Charité CCM, Tucholskystr. 2, 10117 Berlin, Germany. E-mail: thomas.haeupl{at}charite.de
| Abstract |
|---|
|
|
|---|
Objective. The factors that induce remission of RA during pregnancy and the relapse occurring after delivery remain an enigma. In a previous study, we investigated gene-expression profiles of peripheral blood mononuclear cells (PBMC) in patients with RA and healthy women in late pregnancy and postpartum. Profiles of samples from both groups were similar in late pregnancy with elevated monocyte and decreased lymphocyte signatures. Postpartum, in RA PBMC the high level of monocyte transcripts persisted. Further increase was observed in adhesion, migration and signalling processes related to monocytes but also in lymphocytes despite similar clinical activity due to intensified drug treatment. This prompted us to investigate correlations between clinical parameters of disease activity and gene profiles.
Methods. Transcriptome data were correlated with RADAI, CRP, monocyte and lymphocyte counts. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotations, monocytes and lymphocytes signatures were used as reference information.
Results. Comparative analysis of PBMC expression profiles from RA patients during and after pregnancy with RADAI and CRP revealed a correlation of these disease activity parameters predominantly with monocyte transcripts. Genes related to cellular programs of adhesion, migration and response to infections were upregulated. Comparing clinically active and not-active RA patients postpartum revealed a cluster of 19 genes that could also identify active disease during pregnancy.
Conclusion. The data suggest that an increase of the RADAI and an elevation of CRP is a consequence of molecular activation of monocytes. Furthermore, they indicate that molecular activation of T lymphocytes may remain clinically unrecognized postpartum. It is conceivable that a set of 19 genes may qualify as molecular disease activity marker.
KEY WORDS: Rheumatoid arthritis, Pregnancy, Gene expression, Correlation analysis, Molecular pathways, Biomarkers, Monocyte activation, Postpartum flare
| Introduction |
|---|
|
|
|---|
The remission of rheumatoid arthritis (RA) during pregnancy occurring in the majority of patients remains an enigma. Improvement of signs and symptoms during pregnancy is followed by a relapse after delivery [1].
Recently we have studied gene expression of peripheral blood mononuclear cells (PBMC) in the third trimester and 24 weeks after delivery in six patients with RA and in eight age-matched healthy women (Häupl et al., submitted). Gene-expression profiles performed by affymetrix analysis revealed an elevated monocyte and a decreased lymphocyte transcriptome during pregnancy in healthy women. Postpartum this proportion was inverted. This may reflect a suppression of the adaptive immune response as part of the tolerance mechanisms towards the semiallogenic fetus. The complementary increase of phagocytes during pregnancy also found by Crocker et al. [2] suggests compensatory mechanisms provided by innate immune functions. In RA patients, the PBMC transcriptome during pregnancy was comparable to the one in healthy pregnant women, but it became activated postpartum with a rise in the level of monocyte and a concomitant increase in lymphocyte transcripts (Häupl et al., submitted).
Traditionally, disease activity in RA is defined by clinical symptoms. The only objective variable of the commonly used Disease Activity Score (DAS) is the erythrocyte sedimentation rate or the acute phase reactant CRP. Recent data clearly showed that behind improvement of clinical symptoms subclinical disease activity—for example, progressive joint destruction—may persist [3, 4]. This constellation prompted a search for correlations between clinical parameters of RA and expressed genes and the comparison of the setting during and after pregnancy.
| Patients and methods |
|---|
|
|
|---|
Patients and controls
Six patients with RA (fulfilling the ACR criteria [5]) and eight age-matched healthy women were studied at gestational week 32–34 and 24 weeks postpartum. The study was performed at the Department of Rheumatology and Clinical Immunology/Allergology of the University Hospital of Bern after approval by the institutional review board of Bern and informed patient consent was obtained. Mean age of pregnant patients was 31 (range 21–38), and of pregnant healthy women 33 (21–40). Disease activity by the RA Disease Activity Index (RADAI) was 1.7 (range 0–8.6). Tender and swollen joint count were 3.8 (range 0–10) and 2.7 (range 0–6), respectively. Mean level of CRP was 20.1 mg/l (range 9–52), monocyte counts 0.63 cells/pl (range 0.41–0.88) and lymphocyte counts 1.91 cells/pl (range 1.28–3.31). During pregnancy the following medications were allowed: NSAIDs until week 32 of gestation, prednisone maximum 10 mg/day throughout pregnancy, antimalarials and sulfasalazine. Flares postpartum were treated according to standard protocols.
RNA isolation and Affymetrix analysis
Gene-expression profiles were determined in PBMC derived from the patients and the healthy women (ND) during the third trimester of pregnancy and 24 weeks after delivery and processed as described (Häupl et al., submitted). Reference data from different highly purified cell types were used to estimate cell-type association of each individual gene as described elsewhere (Häupl et al., submitted). In brief, after extraction by Qiagen RNeasy Mini Kit according to the manufacturer's protocol (Qiagen, Hilden, Germany), 3 µg of total RNA were used for in vitro transcription and biotin-labeling of cRNA (ENZO Biochem, New York, NY, USA). Total 50 µg/ml cRNA was hybridized on HG-U133A GeneChips (Affymetrix, Santa Clara, CA, USA) for 16 h at 45°C, subsequently washed, stained (fluidic station) and scanned (Hewlett Packard Genearray Scanner, Affymetrix).
Statistical analysis
Spearman rank correlation was performed for all RA patients during and after pregnancy between clinical parameters and expression levels for all genes on the microarray. Cell-type association scores for individual genes were calculated as previously defined (Häupl et al., submitted), and compared for the top 100 genes with best correlation to RADAI, CRP, monocytes and lymphocytes. Functional association with the 34 KEGG pathways on cell communication, immune system, infectious diseases, signal transduction, signaling molecules and interaction as well as cell growth and death was estimated by the percentage of genes that revealed a correlation coefficient of R
0.5. For hierarchical clustering, the Genesis software tool was applied [6].
| Results and discussion |
|---|
|
|
|---|
Correlation of molecular data with clinical criteria of disease activity
Disease activity of RA was absent or very low in the third trimester in all but one patient. It increased in most patients after delivery and required immunosuppressive treatment. Spearman rank correlation was performed for all RA patients during and after pregnancy comparing clinical parameters (RADAI, CRP, lymphocyte and monocyte counts) and expression levels for all genes on the microarray. Each of the top 100 genes identified for correlation either with RADAI, CRP, monocyte or lymphocyte blood count revealed coefficients between R = 0.99 and 0.63. Scores for association of these genes with the major blood cell types were calculated based on comparative analysis of reference signatures for each cell type as described elsewhere (Häupl et al., submitted). RADAI- and CRP-related gene sets both displayed a dominant role of monocyte transcripts (mean score ± S.E.M; RADAI: 0.16 ± 0.02; CRP: 0.16 ± 0.02) compared to CD4 T-cell (RADAI: –0.11 ± 0.02; CRP: –0.10 ± 0.02) or CD19 B-cell transcripts (RADAI: –0.08 ± 0.02; CRP: –0.07 ± 0.02).
The analysis of the data demonstrated for the first time a correlation between RADAI as well as CRP with monocyte but not with lymphocyte profiles. Of note, the correlation could be found despite treatment of patients with different immunosuppressive drugs. One might argue that this correlation could be expected as drugs targeting primarily monocytic functions (e.g. TNF blocking agents [7], tocilizumab [8]) show a decrease in CRP levels and clinical signs of inflammation much faster and more potently than drugs targeting lymphocytes or lymphocyte functions (rituximab, cyclosporin, abatacept). Despite differences in kinetics lymphocyte inhibition is also effective to repress disease progression in RA [9–11]. With our current findings of a dominant monocyte profile in active disease but also a broad panel of expression changes related to lymphocytes during postpartum disease reactivation (Häupl et al., submitted), one is tempted to speculate that a defined molecular profile would perform better in the prediction of therapeutic efficacy of a given drug than clinical symptoms.
Correlation with genes of molecular pathways
For functional interpretation, 1822 genes of six relevant pathway groups with 34 pathways defined in the KEGG database were investigated. The portion of genes per pathway was determined that correlated with RADAI, CRP or the calculated PBMC lymphocyte or monocyte fraction by a coefficient of R
0.5. Based on this portion, pathways were ranked all together or for each of the six groups separately (Table 1).
|
Genes correlating with the RADAI were associated with adhesion, leucocyte endothelial migration, infection-related processes including TLR signaling and cholera-associated responses as well as signal transduction in the Notch-, phosphatidylinositol- and mTOR-signalling pathway. Correlation with CRP revealed genes also associated with cell adhesion and communication, responses in infections (cholera, H. pylori) and Notch signalling as well as genes associated with the hematopoietic lineage, the complement and coagulation cascade, B-cell receptor and ErbB signalling.
The monocyte fraction correlated with genes predominantly associated with infection-related processes, with migration, ECM interaction, hematopoietic lineage, antigen processing, Calcium and Notch signalling. The lymphocyte fraction correlated with genes of the T-cell and B-cell receptor signalling pathway, with cell growth and death-related gene activity and with the phosphatidylinositol-, mTOR-, ErbB- and VEGF-signalling pathways.
Taken together, correlation with functional processes revealed differences between the four parameters. While RADAI, CRP and monocyte fraction-related processes focused on adhesion, migration, response to infections and Notch-signalling, lymphocyte fractions were predominantly associated with typical lymphoid receptor signalling pathways, different pathways of the signalling cascade and cell growth.
As lymphocytes are suppressed during pregnancy (Häupl et al., submitted) [12] and increase postpartum, one can speculate that lymphocytes are unleashed after pregnancy, proliferate and differentiate, in order to upregulate adaptive immune functions. In light of the discrepancy that RA patients benefit from pregnancy with its dominance of phagocytes but flare postpartum when lymphocyte function is recurring [12] by further increase of the monocyte dominance, a role of lymphocytes in this process seems obvious. While expression of lymphocyte marker genes were not or even negatively correlated with markers of disease activity, phosphatidylinositol- and mTOR-signalling were correlating with both, PBMC lymphocyte fractions and disease activity and not with the PBMC monocyte fraction. This suggests that both signalling processes reflect lymphocyte activation. Targeting lymphocytes in RA via the mTOR pathway has been discussed earlier [13]. Genes of the phosphatidylinositol pathway are involved in many different cellular processes. In our study, especially expression of phosphoinositide 3-kinases was correlating with the lymphocyte fraction, kinases that have been implicated in inflammatory processes of RA including T cells but also other cell types [14].
Genes differentially expressed in active disease
In a third step of analysis, we addressed the question whether a discrete panel of genes might serve as molecular disease activity marker set. A pattern consisting of 19 genes was identified that reflected disease activity in this cohort of patients. It was derived by comparative analysis of RA patients postpartum with high and low disease activity as defined by increased levels of CRP and/or RADAI. With this gene selection we identified also one patient with active disease throughout pregnancy and excluded all other patients as well as all normal controls during and after pregnancy (Fig. 1). No relevant overlap existed with genes identified by others except from members of the HLA class II and the S100 family [15–18]. The increase of S100 family members in our study was concordant with results reported by others [15, 18]. In contrast, HLA expression was found decreased in RA compared to healthy controls by Bovin et al. [15] and van der Pouw Kraan et al. [17] using PBMC and PAX gene, respectively. However, analyses by Lequerré et al. [16] support our findings of increased expression of HLA class II which they found elevated in non-responders towards infliximab.
|
| Conclusion |
|---|
|
|
|---|
Taken together this study provides new insights into the correlation between gene-expression profiles of PBMC and disease activity in RA. In our previous analysis (Häupl et al., submitted) we observed an interesting and unexpected finding. There was concordance between molecular and clinical disease activity of RA in the third trimester. However, in spite of low disease activity under therapy with immunosuppressive drugs postpartum, we found a consistent difference of gene expression between all RA patients and all healthy women after delivery. In our correlation analysis we could identify particular sets of marker genes that can indicate clinical disease activity, not only postpartum but also during pregnancy. Nevertheless, prospective studies will have to prove whether these sets of genes identified as markers of active RA qualify for quantification of disease activity and eventually may help to predict disease outcome.
| Acknowledgements |
|---|
|
|
|---|
Funding: The work was supported by the German Science Foundation (DFG grant HA-2267/2), the German Federal Ministry of Education and Research (BMBF) through the National Genome Research Network (Infection & Inflammation Network SIPAGE, grant 01GS0413), by the Swiss National Science Foundation (P.V. and M.O.), by the Kurt and Senta Herrmann Foundation, Liechtenstein (P.V. and M.O.) and the EU project AutoCure (grant 018661).
Disclosure statement: The authors have declared no conflicts of interest.
| References |
|---|
|
|
|---|
- Østensen M, Fuhrer L, Mathieu R, Seitz M, Villiger PM. A prospective study of pregnant patients with rheumatoid arthritis and ankylosing spondylitis using validated clinical instruments. Ann Rheum Dis (2004) 63:1212–7.
[Abstract/Free Full Text] - Crocker IP, Baker PN, Fletcher J. Neutrophil function in pregnancy and rheumatoid arthritis. Ann Rheum Dis (2000) 59:555–64.
[Abstract/Free Full Text] - Wais T, Fierz W, Stoll T, Villiger PM. Subclinical disease activity in systemic lupus erythematosus: immunoinflammatory markers do not normalize in clinical remission. J Rheumatol (2003) 30:2133–9.
[Abstract/Free Full Text] - Baraliakos X, Listing J, Brandt J, et al. Radiographic progression in patients with ankylosing spondylitis after 4 yrs of treatment with the anti-TNF-alpha antibody infliximab. Rheumatology (2007) 46:1450–3.
[Abstract/Free Full Text] - Arnett FC, Edworthy SM, Bloch DA, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum (1988) 31:315–24.[Web of Science][Medline]
- Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics (2002) 18:207–8.
[Abstract/Free Full Text] - Elliott MJ, Maini RN, Feldmann M, et al. Treatment of rheumatoid arthritis with chimeric monoclonal antibodies to tumor necrosis factor alpha. Arthritis Rheum (1993) 36:1681–90.[Web of Science][Medline]
- Maini RN, Taylor PC, Szechinski J, et al. Double-blind randomized controlled clinical trial of the interleukin-6 receptor antagonist, tocilizumab, in European patients with rheumatoid arthritis who had an incomplete response to methotrexate. Arthritis Rheum (2006) 54:2817–29.[CrossRef][Web of Science][Medline]
- Edwards JC, Szczepanski L, Szechinski J, et al. Efficacy of B-cell-targeted therapy with rituximab in patients with rheumatoid arthritis. N Engl J Med (2004) 350:2572–81.
[Abstract/Free Full Text] - Genovese MC, Becker JC, Schiff M, et al. Abatacept for rheumatoid arthritis refractory to tumor necrosis factor alpha inhibition. N Engl J Med (2005) 353:1114–23.
[Abstract/Free Full Text] - Tugwell P, Pincus T, Yocum D, et al. Combination therapy with cyclosporine and methotrexate in severe rheumatoid arthritis. The Methotrexate-Cyclosporine combination study group. N Engl J Med (1995) 333:137–41.
[Abstract/Free Full Text] - Tallon DF, Corcoran DJ, O'Dwyer EM, Greally JF. Circulating lymphocyte subpopulations in pregnancy: a longitudinal study. J Immunol (1984) 132:1784–7.[Abstract]
- Temsirolimus: CCI 779, CCI-779, cell cycle inhibitor-779. Drugs R D 2004. 5:363–7.
- Rommel C, Camps M, Ji H. PI3K delta and PI3K gamma: partners in crime in inflammation in rheumatoid arthritis and beyond? Nat Rev Immunol (2007) 7:191–201.[CrossRef][Web of Science][Medline]
- Bovin LF, Rieneck K, Workman C, et al. Blood cell gene expression profiling in rheumatoid arthritis. Discriminative genes and effect of rheumatoid factor. Immunol Lett (2004) 93:217–26.[CrossRef][Web of Science][Medline]
- Lequerré T, Gauthier-Jauneau AC, Bansard C, et al. Gene profiling in white blood cells predicts infliximab responsiveness in rheumatoid arthritis. Arthritis Res Ther (2006) 8:R105.[CrossRef][Medline]
- van der Pouw Kraan TC, van Baarsen LG, Wijbrandts CA, et al. Expression of a pathogen-response program in peripheral blood cells defines a subgroup of Rheumatoid Arthritis patients. Genes Immun (2007) (in press).
- Parker A, Izmailova ES, Narang J, et al. Peripheral blood expression of nuclear factor-kappab-regulated genes is associated with rheumatoid arthritis disease activity and responds differentially to anti-tumor necrosis factor-alpha versus methotrexate. J Rheumatol (2007) 34:1817–22.
[Abstract/Free Full Text]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
M. Cutolo, M. Matucci-Cerinic, M. Lockshin, and M. Ostensen Introduction: new trends in pregnancy and rheumatic diseases Rheumatology, June 1, 2008; 47(suppl_3): iii1 - iii1. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

