Serum Metabolomic And Lipidomic Profiling Identifies Diagnostic Biomarkers For Seropositive And Seronegative Rheumatoid Arthritis Patients

doi:10.21203/rs.3.rs-816146/v1

Download PDF

Research Article

Serum Metabolomic And Lipidomic Profiling Identifies Diagnostic Biomarkers For Seropositive And Seronegative Rheumatoid Arthritis Patients

https://doi.org/10.21203/rs.3.rs-816146/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 01 Dec, 2021

Read the published version in Journal of Translational Medicine →

You are reading this latest preprint version

Background: Diagnosing seronegative rheumatoid arthritis (RA) can be challenging due to complex diagnostic criteria. We sought to discover diagnostic biomarkers for seronegative RA cases by studying metabolomic and lipidomic changes in RA patient serum.

Methods: We carried out metabolomic and lipidomic profiling of metabolites and lipids in serum of 225 RA patients and 100 normal controls. These samples were divided into a discovery set (n = 243) and a validation set (n = 82). A machine-learning-based multivariate classification model was constructed using distinctive metabolites and lipids signals.

Results: Twenty-six metabolites and lipids were identified from the discovery cohort to construct a RA diagnosis model. The model was subsequently tested on a validation set and achieved accuracy of 90.2%, with sensitivity of 89.7% and specificity of 90.6%. Both seropositive and seronegative patients were identified using this model. A co-occurrence network using serum omics profiles was built and parsed into six modules, showing significant association between the inflammation and immune activity markers and aberrant metabolism of energy metabolism, lipids metabolism and amino acid metabolism. Acyl carnitines (20:3), aspartyl-phenylalanine, pipecolic acid, phosphatidylethanolamine PE (18:1) and lysophosphatidylethanolamine LPE (20:3) were positively correlated with the RA disease activity, while histidine and phosphatidic acid PA (28:0) were negatively correlated with the RA disease activity.

Conclusions: A panel of 26 serum markers were selected from omics profiles to build a machine-learning-based prediction model that could aid in diagnosing seronegative RA patients. Potential markers were also identified in stratifying RA cases based on disease activity.

Molecular Biology

General Biochemistry

Rheumatoid arthritis

Seronegative

Metabolomic

Lipidomic

Rheumatoid arthritis (RA) is a complex chronic autoimmune disease with variable presenting symptoms and serum autoantibody test results. Currently diagnosing RA is primarily based on clinical symptoms and the presence of various serum autoantibodies including rheumatoid factor (RF) and anti-citrullinated protein antibody (ACPA) [1–3]. Although most patients with confirmed RA have an abnormal test for RF and/ or ACPA, about 15–20% of cases do not have the elevated levels of RF and ACPA[4]. Early recognition and assessment of the progression of RA are of paramount importance to initiate disease-modifying antirheumatic drugs (DMARDs). Therefore, developing reliable diagnostic biomarkers for evaluating seronegative RA is needed [1].

Chronic inflammatory state in RA causes metabolic changes that can serve as biomarkers for diagnosis, disease activity and treatment efficacy monitoring [5, 6]. The liquid chromatography-mass spectrometry (LC-MS)-based metabolomic and lipidomic approach is a powerful tool to discover potential biomarkers and metabolic pathway remodeling occurred in chronic diseases including RA [4, 7]. However, small sample size and inadequate number of analytes studied so far limit definitive conclusions drawing [4, 8–12]. The potential of using LC-MS-based metabolomic and lipidomic profiling for novel biomarker discovery to achieve more accurate diagnosis and disease activity monitoring remains largely unexplored.

In our study, we aimed to build a classification model that can be applied to both seropositive and seronegative patients by integrative LC-MS-based metabolomic and lipidomic profiling of serum samples from a large cohort of RA patients and normal control subjects.

Clinical samples collection

Serum samples were collected from 225 RA patients and 100 normal control subjects at the Affiliated Hospital of Nanjing University of Chinese Medicine. The clinical and demographic characteristics of the study were summarized in Table 1. Disease activity for RA patients was assessed at the time of serum collection through quantification of tender and swollen joints, erythrocyte sedimentation rate (ESR), C-reactive protein (CRP) levels and the visual analog scale-grading of the patients. Cytokine concentrations were determined by immunoassay. Subjects in any one or more of the following categories were excluded from our analysis: 1) the presence of type I or II diabetes; 2) active viral and/or bacterial infection; 3) the presence of osteoarthritis. The study was approved by the medical ethics committee of the Affiliated Hospital of Nanjing University of Chinese Medicine and followed the tenets of the Declaration of Helsinki (2018NL-106-02). Written informed consents were obtained from all study subjects. Patients were clinically diagnosed with RA according to the American College of Rheumatology and European League Against Rheumatism (EULAR) 2010 criteria [8,13]. Venous blood was collected in the morning before breakfast from all the participants, and then sera were separated and stored at −80 °C until use.

Sample preparation

The metabolites and lipids in serum were isolated and treated as previously reported [14] with slight modifications. Fifty microliters of thawed serum samples were precipitated by adding 200μL of cold acetonitrile. After centrifugation at 14000g for 10 min at 4°C, the supernatant was divided into four fractions: two for polar metabolites analysis by hydrophilic interaction liquid chromatography (HILIC)-LC-MS methods with positive ion mode electrospray ionization (ESI) and negative ion mode ESI, namely method 1 and method 2, respectively; two for lipids analysis by reverse phases (RP)-LC-MS methods with positive ion mode ESI and negative ion mode ESI, namely method 3 and method 4, respectively. After centrifugation, samples were dried and stored at −80°C until use. The quality control (QC) samples were prepared by mixing equal volumes of sera from RA patients and controls before sample preparation as they were aliquoted for analysis. These QC samples were utilized to estimate a “mean” profile representing all analytes encountered during the LC-MS analysis [15].

Metabolomic and lipidomic analysis

For the polar metabolite profiling, dried samples were reconstituted in 100μl of 80% methanol and analyzed by using a Dionex U3000 LC system coupled online to a Q Exactive Orbitrap MS instrument (Thermo Fisher Scientific, MA, USA) set at 35000 resolution (at m/z 200). For the lipids profiling, dried samples were reconstituted in 100μl of acetonitrile for instrument analysis. (HILIC)-LC-MS and (RP)-LC-MS were performed for metabolomic and lipidomic profiling, respectively. The details of the analytical experiments were described in Supplemental Materials. Data pre-processing was carried out as previously reported [16]. The XCMS package [17] was used for the extraction of peak abundances of metabolites and lipids. The background noises and contaminants in the peak tables were filtered by using CPVA with default parameters [18]. Quality assurance was achieved by the statTarget package [19] as follows. Briefly, peaks with more than 50% missing values were removed. The intensity of remaining peaks in samples was corrected according to the QC-RFSC algorithm. Principal component analysis (PCA)-based data quality evaluation was provided in Supplemental materials (Figure S1). Polar metabolites identification was according to the MS/MS spectra, the retention time of commercially available standard compounds, and the accurate mass of compounds from the Human Metabolome Database (www.hmdb.ca) and MassBank of North America (http://mona.fiehnlab.ucdavis.edu). Lipids were identified based on MS/MS match by using the LipidSearch software. Pathway analysis based on identified metabolites was carried out using Metaboanalyst 4.0 [20] according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (www.genome.jp/kegg/).

Data analyses

The Mann–Whitney–Wilcoxon test or Welch’s t-test, was performed to measure the significance of each peak, with results adjusted for multiple testing using false discovery rates (FDR) correction. The cross-validated partial least squares discriminant analysis (PLS-DA) and variable importance in the projection (VIP_plsda) were calculated using the statTarget package. Serum metabolites passing the threshold of VIP_plsda > 1, fold change > 1.2 or < 0.8, and adjusted p-value < 0.05 were regarded as potential markers between the RA and NCs groups. The alluvial plots and forest plot were performed using the ggalluvial package and forestplot package, respectively. The co-occurrence network was carried out and visualized in the igraph package (https://igraph.org). The ordinal regression was performed by using the ordinalgmifs package for fitting an ordinal response model [21]. Coupling the receiver operating characteristic curve (ROC) with its area under the curve (AUC), a widely used method to estimate the diagnostic potential of a classifier in clinical applications, was performed using the pROC package [22]. All packages were implemented using the freely available R language.

Clinical data and patient characteristics

The demographic information of the study participants was summarized in Table 1. Among 225 RA patients, the positive rate of RF and ACPA was 57.8% and 42.2%, respectively. These samples were randomly divided into two independent cohorts. The discovery set consisted of 172 RA patients and 71 normal controls (NCs), and the validation set consisted of 53 RA patients and 29 NCs (Figure 1A). The alluvial plot shown in the Figure 1B displayed the number of individuals across the disease status, gender, disease activity score (DAS28-CRP and DAS28-ESR), as well as the status of autoantibodies such as rheumatoid factor (RF), ACPA, anti-nuclear antibodies (ANA), and anti-keratin antibodies (AKA). The positive rates of RF, ACPA, ANA, AKA in the RA patients were 57.8%, 42.2%, 49.4% and 49.4% (Table 1 and Figure 1B).

Metabolomic and lipidomic profile (MLP)

In four datasets generated from four different LC-MS runs, there were 1697 and 1960 high-confident peaks identified for polar metabolites with method 1 and 2, and 4771 and 3973 high-confident peaks identified for lipids with method 3 and 4. From this 12401 high-confident peaks we identified 265 candidate disease-associated metabolites and lipids (Table S1), including 38 organic acids, 10 amines, 37 amino acids, 5 nucleotides, 4 bile acids, 26 acyl-carnitines (AcCa), 118 glycerophospholipids, 9 sphingolipids and 2 glycerolipids. These metabolites represent enriched metabolic pathways involving carnitine synthesis, oxidation of branched-chain fatty acids, biotin metabolism, malate-aspartate shuttle, citric acid cycle, urea cycle, phenylalanine and tyrosine metabolism, phospholipid biosynthesis and histidine metabolism (Figure S2).

Clustering of test samples using MLP and pathway enrichment analysis

To analyze whether any of the identified metabolites or lipids are associated with RA, we carried out the univariate and multivariate analysis of the identified metabolites and lipids. The well-established partial least squares-discriminant analysis (PLS-DA) demonstrated excellent separation of NCs from seropositive RA and seronegative RA based on MLP (Figure 2A). The permutation tests (n = 1000) were performed to validate the PLS-DA model we used (Figure S3). This separation clearly demonstrated the difference in serum metabolite and lipid levels that existed between RA patients and normal control subjects. Statistical analysis identified 68 significantly upregulated and 38 downregulated serum metabolites and lipids that were correlated with RA (Figure 2B). Pathway enrichment analysis (Figure 2C) highlighted the top enriched metabolic pathways in RA. Specifically, metabolic products associated with Warburg effect, pentose phosphate metabolism, glycolysis, and lipid metabolism products involving phospholipid, sphingolipid, oxidation of branched chain fatty acids and carnitine synthesis were upregulated in the RA group, while histidine metabolism was downregulated.

Constructing a multivariate classification model for RA

We constructed a multivariate classification model including 26 metabolites and lipids that could separate seropositive and seronegative RA cases from normal controls (Figure 2D-E). We evaluated three machine learning algorithms such as binary logistic regression, random forest and support vector machine on the above-mentioned metabolites and lipids for the classification of RA cases using the discovery set (n = 243). The binary logistic regression algorithm based model trained with leave-one-out cross-validation had an accuracy of 100% (Figure 2E, AUC = 1.00). This model was tested using the independent validation set (n = 82, Table S2), and showed AUC = 0.91 with test sensitivity of 89.7% and specificity of 90.6% (Figure 2E). We did an analysis of the 5 misclassified RA cases, and it is noteworthy that all of the 5 cases had borderline to positive RF antibody levels (Table S3).

Association of MLP and clinical parameters

To understand the association of MLP with clinical parameters of RA, we used a correlation-based network approach to explore the co-occurrence patterns between 28 clinical variables and serum concentrations of 265 identified metabolites and lipids in the RA patients, which could help to explore the correlation of metabolic products and clinical confounders. In the co-occurrence network (Figure 3), nodes represent clinical parameters or metabolites and lipids, and two nodes are connected if they are significantly correlated (adjusted p-value < 0.05 and r-value > 0.2). The resulting network consisted of 117 nodes and 274 edges. The modularity index was 0.467, which is above 0.4 suggested for a modularly structured network [23]. Overall, the entire network could be parsed into six modules. As shown in Figure 3, module 1 contained CRP, leukocyte, neutrophil and ACPA related to inflammatory activity and 17 metabolites and lipids. Modules 2 and 3 contained RF, IgG, IgM, IgA, globulin, ESR, ANA and AKA related to the immune response and 22 metabolites and lipids. The general clinical variables such as age, gender and BMI were clustered in modules 4-5 and were far away from modules 1-3. Therefore, the co-occurrence network analysis revealed that there was a clear association between the serum inflammation/immune markers and concentrations of metabolites and lipids in the RA patients.

Stratification of disease activity using MLP

To analyze the utility of using MLP for stratification of RA cases based on disease activity, twenty metabolites and lipids strongly correlated to the inflammatory/immune activity were selected from modules 1 and 2 in the co-occurrence network (Figure 3). Sixteen metabolites and lipids had significantly increased odd ratios for RA status, while 4 metabolites and lipids had significantly decreased odd ratios (Figure 4). Disease activity score DAS28-CRP is often utilized to evaluate the disease activity of RA patients [24,25]. Using DAS28-CRP to stratify RA patients, we further identify 7 metabolites and lipids that were strongly associated with disease activity categories by using the ordinal regression method (Figure 5). Among these, AcCa (20:3), aspartyl-phenylalanine (asp-phe), pipecolic acid, lysophosphatidylethanolamine (LPE 18:1) and LPE (20:3) appeared to be positively correlated with higher RA disease activity, while histidine and PA (28:0) were negatively correlated with RA disease activity (Figure 5A). AcCa (20:3), LPE (18:1) and LPE (20:3) were significantly increased in the remission-low risk (R-L) group compared to the NCs group (p < 0.05, AUC > 0.7). A small peptide asp-phe and lipid LPE (18:1) were significantly increased in the DAS28-CRP high disease activity (HIGH) group compared with the R-L group (p < 0.05, AUC > 0.7), indicating their potential role as biomarkers in RA disease activity stratification (Figure 5B).

RA is a highly heterogeneous disease with variable presenting symptoms, severity and response to treatment. Autoimmune damage may happen years before symptoms occur and clinical diagnosis is made. It is important for early and accurate diagnosis of RA and prompt initiation of effective treatment to prevent joint damage and functional loss [26]. Current diagnostic criteria are based on comprehensive evaluation of symptoms, serology status and acute phase reactant levels. Because a significant percentage of RA patients are negative for serologic markers, additional diagnostic methods are actively being developed to help increase diagnostic accuracy, particularly for seronegative RA patients. The chronic inflammation and joint destruction in RA patients may cause metabolic perturbations in the peripheral blood, providing opportunities to discover potential biomarkers to improve the clinical diagnosis of RA [4, 27].

Previous studies have reported metabolic changes in RA by using nuclear magnetic resonance (NMR) or RPLC-MS [4, 9, 10]. We developed a multi-platform strategy including two reverse phase and two HILIC chromatography hyphenated to high-resolution mass spectrometry (HRMS) in positive and negative ionization modes for the simultaneous quantitative analysis of metabolites and lipids. By means of the integrative metabolomic and lipidomic analysis, we were able to identify 26 serum metabolites and lipids that correctly classified RA patients from normal control subjects. A multivariate classification model was derived from the discovery set and was subsequently validated using an independent validation set. It is noteworthy that our model combined with RF test reached 100% validation accuracy, and all patients with seronegative RA were identified by using our model, thus showing the potential for providing a diagnostic tool in both seropositive and seronegative RA.

Similar findings pertaining to the differences in the levels of lysine, phenylalanine, proline, ornithine and salicylic acid between the RA patients and normal control subjects were observed previously [10, 28]. These metabolites have been associated with immune system activation [4, 29, 30]. On the other hand, we found serum level of histidine, which has anti-inflammatory effects and antioxidative stress effect via inhibiting peroxisome proliferator-activated receptor γ (PPARγ) -involved pathways [31], was significantly decreased in RA patients. The role of histidine in RA needs further studies since histidine supplementation did not show an advantage in a double-blind trial [32]. In our study histidine level was inversely correlated with the RA disease activity defined by DAS28-CRP. Histidine level was significantly decreased in patients with moderate and high disease activity, while no significant changes were seen in RA patients with remission and low disease activity (Fig. 5A).

Lipids are implicated in diverse biological functions, including being the major component of cell membranes and the regulation of cell migration and inflammation [7, 33]. Our study found that increased levels of lipid mediators such as sphingosine 1-phosphate (S1P) and 12-HETE and lysophospholipids such as LPE (18:1), LPE (18:2), lysophosphatidylcholine (LPC 18:1), LPC (18:4) and lysophosphatidylserine (LPS 18:0) in RA patients were positively correlated with CRP levels (Fig. 3), supporting their roles in inflammation. Other groups have reported that LPC was shown to induce cyclooxygenase-2 expression in the vascular endothelial cells, playing a possible pro-inflammatory role [34]. S1P is generated from sphingosine by activation of sphingosine kinase and has been implicated as a potential therapeutic target in RA [35]. Moreira et al. reported that 12-HETE, an arachidonic acid-derived metabolite, may play an important role in the chronic inflammatory process associated with RA by mediating proinflammatory actions [36].

Ongoing inflammation and immune activation are characterized by rising energy demand that is required for immune cell growth, proliferation and the production of proinflammatory molecules [37]. In our study, the increased levels of metabolites (e.g., alpha-ketoglutaric acid, AcCa (3:0), fructose-6-phosphate and glucose-6-phosphate) seen in the serum of RA patients indicate enhanced energy metabolism involving glycolysis, citric acid cycle and fatty acid oxidation (Fig. 4). Methotrexate (MTX), a first line treatment for RA, is an antimetabolite drug that targets folic acid metabolism. Many ongoing clinical trials are investigating potential drugs that target lipid/glucose metabolic pathways to improve inflammation control and disease outcome [5]. In addition, using metabolomic approach to delineate inflammation-associated metabolic status of RA patients may offer a promising method for disease activity and treatment effect monitoring.

Our study has both strengths and limitations. Its strengths include that it employs a comprehensive multi-platform approach for metabolomic and lipidomic analysis, large sample size and its potential to provide diagnostic value for both seropositive and seronegative RA patients. Its limitations include that it is still a discovery study that shows correlation of serum markers to RA, but does not directly address `the mechanisms of these markers in RA. For developing our method into clinical practice, our classification model still needs to be evaluated in other multinational/multiethnic cohorts.

In conclusion, we have built a serum metabolic/lipidomic- markers-based model that has potential diagnose value. Our study suggests the integrative comprehensive metabolomic and lipidomic profiling is a promising system biology approach for uncovering biomarkers useful for RA diagnosis and disease activity monitoring.

Ethics approval and consent to participate: The study was approved by the medical ethics committee of the Affiliated Hospital of Nanjing University of Chinese Medicine. Written informed consents were obtained from all study subjects.

Consent for publication: Written informed consent for publication of their clinical details and/or clinical images was obtained from the patient/parent/guardian/ relative of the patient.

Availability of data and materials: All data generated or analyzed during this study are included in this published article [and its supplementary information files].

Competing interests: The authors declare no conflict of interest.

Funding: This research received no external funding.

Author contribution: Hemi Luan: Conceptualization, Methodology, Investigation, Writing – Original draft, Writing – Review & editing. Wanjian Gu: Resources. Hua Li: Resources. Zi Wang: Investigation, Writing – Original draft. Lu Lu: Resources. Mengying Ke: Resources. Jiawei Lu: Resources. Wenjun Chen: Resources. Zhangzhang Lan: Writing – Review & editing. Yanlin Xiao: Writing – Review & editing. Jinyue Xu: Writing – Review & editing. Yi Zhang: Writing – Review & editing. Zongwei Cai: Conceptualization,. Shijia Liu: Resources, Conceptualization. Wenyong Zhang: Methodology, Writing – Review & editing.

Acknowledgements: The authors would like to acknowledge the financial support from the National Natural Science Foundation of China (ref. no.: 21904058, ref. no.: 21904057 and ref. no.: 81774096), the Open Projects of the Discipline of Chinese Medicine of Nanjing University of Chinese Medicine Supported by the Subject of Academic Priority Discipline of Jiangsu Higher Education Institutions (NO. ZYX03KF031), the Young Elite Scientists Sponsorship Program by CAST (ref. no.: QNRC2-B04).

Bugatti S, Manzo A, Montecucco C, Caporali R. The Clinical Value of Autoantibodies in Rheumatoid Arthritis. Front Med [Internet]. 2018;5:1–10. Available from: https://www.frontiersin.org/article/10.3389/fmed.2018.00339/full.
Saraux A, Berthelot JM, Gérard Chalès, Le Henaff C, Mary JY, Thorel JB, et al. Value of laboratory tests in early prediction of rheumatoid arthritis. Arthritis Care Res (Hoboken) [Internet]. 2002;47:155–65. Available from: http://doi.wiley.com/10.1002/art.10241.
Wasserman A. Rheumatoid Arthritis: Common Questions About Diagnosis and Management. Am Fam Physician [Internet]. 2018;97:455–62. Available from: http://www.ncbi.nlm.nih.gov/pubmed/29671563.
Souto-Carneiro M, Tóth L, Behnisch R, Urbach K, Klika KD, Carvalho RA, et al. Differences in the serum metabolome and lipidome identify potential biomarkers for seronegative rheumatoid arthritis versus psoriatic arthritis. Ann Rheum Dis [Internet]. 2020;79:499–506. Available from: https://ard.bmj.com/lookup/doi/10.1136/annrheumdis-2019-216374.
Chimenti MS, Triggianese P, Conigliaro P, Candi E, Melino G, Perricone R. The interplay between inflammation and metabolism in rheumatoid arthritis. Cell Death Dis [Internet]. Nature Publishing Group; 2015;6:e1887–e1887. Available from: http://dx.doi.org/10.1038/cddis.2015.246.
Kerekes G, Nurmohamed MT, González-Gay MA, Seres I, Paragh G, Kardos Z, et al. Rheumatoid arthritis and metabolic syndrome. Nat Rev Rheumatol [Internet]. Nature Publishing Group; 2014;10:691–6. Available from: http://dx.doi.org/10.1038/nrrheum.2014.121.
van Halm VP, Nielen MMJ, Nurmohamed MT, van Schaardenburg D, Reesink HW, Voskuyl AE, et al. Lipids and inflammation: serial measurements of the lipid profile of blood donors who later developed rheumatoid arthritis. Ann Rheum Dis [Internet]. 2006;66:184–8. Available from: https://ard.bmj.com/lookup/doi/10.1136/ard.2006.051672.
Carlson AK, Rawle RA, Wallace CW, Adams E, Greenwood MC, Bothner B, et al. Global metabolomic profiling of human synovial fluid for rheumatoid arthritis biomarkers. Clin Exp Rheumatol [Internet]. 2019;37:393–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/30620276.
Dubey D, Kumar S, Chaurasia S, Guleria A, Ahmed S, Singh R, et al. NMR-Based Serum Metabolomics Revealed Distinctive Metabolic Patterns in Reactive Arthritis Compared with Rheumatoid Arthritis. J Proteome Res [Internet]. 2018;18:acs.jproteome.8b00439. Available from: https://pubs.acs.org/doi/10.1021/acs.jproteome.8b00439.
Li J, Che N, Xu L, Zhang Q, Wang Q, Tan W, et al. LC-MS-based serum metabolomics reveals a distinctive signature in patients with rheumatoid arthritis. Clin Rheumatol [Internet]. 2018;37:1493–502. Available from: http://link.springer.com/10.1007/s10067-018-4021-6.
Sasaki C, Hiraishi T, Oku T, Okuma K, Suzumura K, Hashimoto M, et al Metabolomic approach to the exploration of biomarkers associated with disease activity in rheumatoid arthritis. Mudiam MKR, editor. PLoS One [Internet]. 2019;14:e0219400. Available from: https://dx.plos.org/10.1371/journal.pone.0219400.
Zabek A, Swierkot J, Malak A, Zawadzka I, Deja S, Bogunia-Kubik K, et al. Application of 1 H NMR-based serum metabolomic studies for monitoring female patients with rheumatoid arthritis. J Pharm Biomed Anal [Internet]. Elsevier B.V.; 2016;117:544–50. Available from: http://dx.doi.org/10.1016/j.jpba.2015.10.007.
Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO, et al. 2010 Rheumatoid arthritis classification criteria: An American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum [Internet]. 2010;62:2569–81. Available from: http://doi.wiley.com/10.1002/art.27584.
Luan H, Meng N, Liu P, Feng Q, Lin S, Fu J, et al. Pregnancy-Induced Metabolic Phenotype Variations in Maternal Plasma. J Proteome Res [Internet]. 2014;13:1527–36. Available from: https://pubs.acs.org/doi/10.1021/pr401068k.
Luan H, Liu L-F, Tang Z, Zhang M, Chua K-K, Song J-X, et al. Comprehensive urinary metabolomic profiling and identification of potential noninvasive marker for idiopathic Parkinson’s disease. Sci Rep [Internet]. Nature Publishing Group; 2015;5:13888. Available from: http://dx.doi.org/10.1038/srep13888.
Ji F, Sreenivasmurthy SG, Wei J, Shao X, Luan H, Zhu L, et al. Study of BDE-47 induced Parkinson’s disease-like metabolic changes in C57BL/6 mice by integrated metabolomic, lipidomic and proteomic analysis. J Hazard Mater [Internet]. Elsevier; 2019;378:120738. Available from: https://doi.org/10.1016/j.jhazmat.2019.06.015.
Mahieu NG, Genenbacher JL, Patti GJ. A roadmap for the XCMS family of software solutions in metabolomics. Curr Opin Chem Biol [Internet]. Elsevier Ltd; 2016;30:87–93. Available from: https://linkinghub.elsevier.com/retrieve/pii/S1367593115001374.
Luan H, Jiang X, Ji F, Lan Z, Cai Z, Zhang W. CPVA: A web-based metabolomic tool for chromatographic peak visualization and annotation. Bioinformatics. 2020.
Luan H, Ji F, Chen Y, Cai Z. statTarget: A streamlined tool for signal drift correction and interpretations of quantitative mass spectrometry-based omics data. Anal Chim Acta [Internet]. Elsevier B.V.; 2018;1036:66–72. Available from: https://doi.org/10.1016/j.aca.2018.08.002.
Chong J, Wishart DS, Xia J. Using. MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr Protoc Bioinforma [Internet]. 2019;68:1–128. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1002/cpbi.86.
Archer KJ, Hou J, Zhou Q, Ferber K, Layne JG, Gentry AE. ordinalgmifs: An R Package for Ordinal Regression in High-dimensional Data Settings. Cancer Inform [Internet]. 2014;13:CIN.S20806. Available from: http://journals.sagepub.com/doi/10.4137/CIN.S20806.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: an open-source package for R and S + to analyze and compare ROC curves. BMC Bioinformatics [Internet]. BioMed Central Ltd; 2011;12:77. Available from: http://www.biomedcentral.com/1471-2105/12/77.
MEHRLE D, STROSSER A, HARKIN A. Walk-modularity and community structure in networks. Netw Sci [Internet]. 2015;3:348–60. Available from: https://www.cambridge.org/core/product/identifier/S205012421500020X/type/journal_article.
Arts EEA, Fransen J, Den Broeder AA, van Riel PLCM, Popa CD. Low disease activity (DAS28 ≤ 3.2) reduces the risk of first cardiovascular event in rheumatoid arthritis: a time-dependent Cox regression analysis in a large cohort study. Ann Rheum Dis [Internet]. 2017;76:1693–9. Available from: https://ard.bmj.com/lookup/doi/10.1136/annrheumdis-2016-210997.
Wells G, Becker J-C, Teng J, Dougados M, Schiff M, Smolen J, et al. Validation of the 28-joint Disease Activity Score (DAS28) and European League Against Rheumatism response criteria based on C-reactive protein against disease progression in patients with rheumatoid arthritis, and comparison with the DAS28 based on erythr. Ann Rheum Dis [Internet]. 2009;68:954–60. Available from: https://ard.bmj.com/lookup/doi/10.1136/ard.2007.084459.
Smolen JS, Aletaha D, Barton A, Burmester GR, Emery P, Firestein GS, et al. Rheumatoid arthritis. Nat Rev Dis Prim [Internet]. Macmillan Publishers Limited; 2018;4:18001. Available from: http://dx.doi.org/10.1038/nrdp.2018.1.
Tomizawa T, Ito H, Murata K, Hashimoto M, Tanaka M, Murakami K, et al. Distinct biomarkers for different bones in osteoporosis with rheumatoid arthritis. Arthritis Res Ther [Internet]. Arthritis Research & Therapy; 2019;21:174. Available from: https://arthritis-research.biomedcentral.com/articles/10.1186/s13075-019-1956-1.
Kim S, Hwang J, Xuan J, Jung YH, Cha H-S, Kim KH. Global Metabolite Profiling of Synovial Fluid for the Specific Diagnosis of Rheumatoid Arthritis from Other Inflammatory Arthritis. Bahn Y-S, editor. PLoS One [Internet]. 2014;9:e97501. Available from: https://dx.plos.org/10.1371/journal.pone.0097501.
Li P, Yin Y-L, Li D, Woo Kim S, Wu G. Amino acids and immune function. Br J Nutr [Internet]. 2007;98:237–52. Available from: https://www.cambridge.org/core/product/identifier/S000711450769936X/type/journal_article.
Qi Y, Pi Z, Liu S, Song F, Lin N, Liu Z. A metabonomic study of adjuvant-induced arthritis in rats using ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry. Mol Biosyst [Internet]. Royal Society of Chemistry; 2014;10:2617. Available from: http://dx.doi.org/10.1039/c4mb00131a.
Sun X, Feng R, Li Y, Lin S, Zhang W, Li Y, et al. Histidine supplementation alleviates inflammation in the adipose tissue of high-fat diet-induced obese rats via the NF-κB- and PPARγ-involved pathways. Br J Nutr [Internet]. 2014;112:477–85. Available from: https://www.cambridge.org/core/product/identifier/S0007114514001056/type/journal_article.
Pinals RS, Harris ED, Burnett JB, Gerber DA. Treatment of rheumatoid arthritis with L-histidine: a randomized, placebo-controlled, double-blind trial. J Rheumatol [Internet]. 1977;4:414–419. Available from: http://europepmc.org/abstract/MED/342692.
Navarro-Millán I, Goyal P, Safford MM. Lipid screening and statins alongside disease-modifying anti-rheumatic drugs for patients with rheumatoid arthritis. Rheumatology [Internet]. 2019;58:933–4. Available from: https://academic.oup.com/rheumatology/article/58/6/933/5110150.
Brouwers H, von Hegedus J, Toes R, Kloppenburg M, Ioan-Facsinay A. Lipid mediators of inflammation in rheumatoid arthritis and osteoarthritis. Best Pract Res Clin Rheumatol [Internet]. 2015;29:741–55. Available from: https://linkinghub.elsevier.com/retrieve/pii/S152169421600005X.
Hu P, Chen Y, Cai P, Jiang L, Wu L. Sphingosine-1-phosphate: a potential therapeutic target for rheumatoid arthritis. Mol Biol Rep [Internet]. 2011;38:4225–30. Available from: http://link.springer.com/10.1007/s11033-010-0545-9.
Moreira V, Gutiérrez JM, Lomonte B, Vinolo MAR, Curi R, Lambeau G, et al. 12-HETE is a regulator of PGE2 production via COX-2 expression induced by a snake venom group IIA phospholipase A2 in isolated peritoneal macrophages. Chem Biol Interact [Internet]. Elsevier B.V.; 2020;317:108903. Available from: https://doi.org/10.1016/j.cbi.2019.108903.
Patel CH, Leone RD, Horton MR, Powell JD. Targeting metabolism to regulate immune responses in autoimmunity and cancer. Nat Rev Drug Discov [Internet]. Springer US; 2019;18:669–88. Available from: http://dx.doi.org/10.1038/s41573-019-0032-5.

Table 1. Clinical and demographic characteristics of the study participants.

ESR, Erythrocyte sedimentation rate; CRP, C-Reactive protein; IgA, Immunoglobulin A; IgG, Immunoglobulin G; IgM, Immunoglobulin M; ACPA, anti-citrullinated protein antibody; ANA, anti-nuclear antibodies; ANA, anti-keratin antibodies. DAS28-CRP, Disease activity score 28-joint count C reactive protein; DAS28-ESR, Disease activity score 28-joint count erythrocyte sedimentation rate. Statistical significance was determined using unpaired two-tailed Student’s t test.

AdditionalFile1.pdf

Download PDF

Journal Publication

published 01 Dec, 2021

Read the published version in Journal of Translational Medicine →

Reviews received at journal
13 Sep, 2021
Reviewers invited by journal
13 Sep, 2021
Editor assigned by journal
17 Aug, 2021
First submitted to journal
15 Aug, 2021

You are reading this latest preprint version

Serum Metabolomic And Lipidomic Profiling Identifies Diagnostic Biomarkers For Seropositive And Seronegative Rheumatoid Arthritis Patients

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

Discussion

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 1