Dataset Information

Plasma metabolomic profiling in patients with rheumatoid arthritis identifies biochemical features predictive of quantitative disease activity.

ABSTRACT:

Background

Rheumatoid arthritis (RA) is a chronic, autoimmune disorder characterized by joint inflammation and pain. In patients with RA, metabolomic approaches, i.e., high-throughput profiling of small-molecule metabolites, on plasma or serum has thus far enabled the discovery of biomarkers for clinical subgroups, risk factors, and predictors of treatment response. Despite these recent advancements, the identification of blood metabolites that reflect quantitative disease activity remains an important challenge in precision medicine for RA. Herein, we use global plasma metabolomic profiling analyses to detect metabolites associated with, and predictive of, quantitative disease activity in patients with RA.

Methods

Ultra-high-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was performed on a discovery cohort consisting of 128 plasma samples from 64 RA patients and on a validation cohort of 12 samples from 12 patients. The resulting metabolomic profiles were analyzed with two different strategies to find metabolites associated with RA disease activity defined by the Disease Activity Score-28 using C-reactive protein (DAS28-CRP). More specifically, mixed-effects regression models were used to identify metabolites differentially abundant between two disease activity groups ("lower", DAS28-CRP ≤ 3.2; and "higher", DAS28-CRP > 3.2) and to identify metabolites significantly associated with DAS28-CRP scores. A generalized linear model (GLM) was then constructed for estimating DAS28-CRP using plasma metabolite abundances. Finally, for associating metabolites with CRP (an indicator of inflammation), metabolites differentially abundant between two patient groups ("low-CRP", CRP ≤ 3.0 mg/L; "high-CRP", CRP > 3.0 mg/L) were investigated.

Results

We identified 33 metabolites differentially abundant between the lower and higher disease activity groups (P < 0.05). Additionally, we identified 51 metabolites associated with DAS28-CRP (P < 0.05). A GLM based upon these 51 metabolites resulted in higher prediction accuracy (mean absolute error [MAE] ± SD: 1.51 ± 1.77) compared to a GLM without feature selection (MAE ± SD: 2.02 ± 2.21). The predictive value of this feature set was further demonstrated on a validation cohort of twelve plasma samples, wherein we observed a stronger correlation between predicted and actual DAS28-CRP (with feature selection: Spearman's ρ = 0.69, 95% CI: [0.18, 0.90]; without feature selection: Spearman's ρ = 0.18, 95% CI: [-0.44, 0.68]). Lastly, among all identified metabolites, the abundances of eight were significantly associated with the CRP patient groups while controlling for potential confounders (P < 0.05).

Conclusions

We demonstrate for the first time the prediction of quantitative disease activity in RA using plasma metabolomes. The metabolites identified herein provide insight into circulating pro-/anti-inflammatory metabolic signatures that reflect disease activity and inflammatory status in RA patients.

SUBMITTER: Hur B

PROVIDER: S-EPMC8185925 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:IntroductionTo identify markers and mechanisms of resistance to adalimumab therapy, we studied global gene expression profiles in synovial tissue specimens obtained from severe rheumatoid arthritis (RA) patients before and after initiation of treatment.MethodsPaired synovial biopsies were obtained from the affected knee of 25 DMARD (disease-modifying antirheumatic drug)-resistant RA patients at baseline (T0) and 12 weeks (T12) after initiation of adalimumab therapy. DAS28-CRP (disease activity score using 28 joint counts-C-reactive protein) scores were computed at the same time points, and patients were categorized as good, moderate, or poor responders according to European League Against Rheumatism criteria. Global gene expression profiles were performed in a subset of patients by means of GeneChip Human Genome U133 Plus 2.0 Arrays, and confirmatory immunohistochemistry experiments were performed on the entire cohort.ResultsGene expression studies performed at baseline identified 439 genes associated with poor response to therapy. The majority (n = 411) of these genes were upregulated in poor responders and clustered into two specific pathways: cell division and regulation of immune responses (in particular, cytokines, chemokines, and their receptors). Immunohistochemistry experiments confirmed that high baseline synovial expression of interleukin-7 receptor alpha chain (IL-7R), chemokine (C-X-C motif) ligand 11 (CXCL11), IL-18, IL-18 receptor accessory (IL-18rap), and MKI67 is associated with poor response to adalimumab therapy. In vitro experiments indicated that genes overexpressed in poor responders could be induced in fibroblast-like synoviocytes (FLS) cultures by the addition of tumor necrosis factor-alpha (TNF-alpha) alone, IL-1beta alone, the combination of TNF-alpha and IL-17, and the combination of TNF-alpha and IL-1beta.ConclusionsGene expression studies of the RA synovium may be useful in the identification of early markers of response to TNF blockade. Genes significantly overexpressed at baseline in poor responders are induced by several cytokines in FLSs, thereby suggesting a role for these cytokines in the resistance to TNF blockade in RA.