Next Article in Journal
Long COVID-19 Pulmonary Sequelae and Management Considerations
Next Article in Special Issue
Serum Interleukin (IL)-23 and IL-17 Profile in Inflammatory Bowel Disease (IBD) Patients Could Differentiate between Severe and Non-Severe Disease
Previous Article in Journal
Elevated De Ritis Ratio as a Predictor for Acute Kidney Injury after Radical Retropubic Prostatectomy
Previous Article in Special Issue
IL17F: A Possible Risk Marker for Spondyloarthritis in HLA-B*27 Negative Brazilian Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Biomarker Candidates for Skeletal Muscle Loss in the Collagen-Induced Arthritis (CIA) Model

by
Paulo V. G. Alabarse
1,*,
Jordana M. S. Silva
2,
Rafaela C. E. Santo
2,
Marianne S. Oliveira
2,
Andrelise S. Almeida
2,
Mayara S. de Oliveira
2,
Mônica L. Immig
2,
Eduarda C. Freitas
2,
Vivian O. N. Teixeira
2,
Camilla L. Bathurst
3,
Claiton V. Brenol
2,
Lidiane I. Filippin
4,
Stephen P. Young
3,
Priscila S. Lora
5 and
Ricardo M. Xavier
2
1
VA San Diego Healthcare System, Veteran Medical Research Foundation, San Diego, CA 92161, USA
2
Hospital de Clínicas de Porto Alegre, Serviço de Reumatologia, Laboratório de Doenças Autoimunes, Porto Alegre 90035-903, Brazil
3
Queen Elizabeth Hospital and National Institute for Health Research, Birmingham Biomedical Research Centre, University of Birmingham, Birmingham B15 2TT, UK
4
Universidade La Salle, Canoas 92010-000, Brazil
5
Universidade do Vale do Rio dos Sinos, São Leopoldo 93022-750, Brazil
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(9), 837; https://doi.org/10.3390/jpm11090837
Submission received: 20 July 2021 / Revised: 10 August 2021 / Accepted: 16 August 2021 / Published: 26 August 2021
(This article belongs to the Special Issue Precision Medicine for Inflammatory and Autoimmune Diseases)

Abstract

:
There is no consensus for diagnosis or treatment of RA muscle loss. We aimed to investigate metabolites in arthritic mice urine as biomarkers of muscle loss. DBA1/J mice comprised collagen-induced arthritis (CIA) and control (CO) groups. Urine samples were collected at 0, 18, 35, 45, 55, and 65 days of disease and subjected to nuclear magnetic resonance spectroscopy. Metabolites were identified using Chenomx and Birmingham Metabolite libraries. The statistical model used principal component analysis, partial least-squares discriminant analysis, and partial least-squares regression analysis. Linear regression and Fisher’s exact test via the MetaboAnalyst website were performed (VIP-score). Nearly 100 identified metabolites had CIA vs. CO and disease time-dependent differences (p < 0.05). Twenty-eight metabolites were muscle-associated: carnosine (VIPs 2.8 × 102) and succinyl acetone (VIPs 1.0 × 10) showed high importance in CIA vs. CO models at day 65; CIA pair analysis showed histidine (VIPs 1.2 × 102) days 55 vs. 65, histamine (VIPs 1.1 × 102) days 55 vs. 65, and L-methionine (VIPs 1.1 × 102) days 0 vs. 18. Carnosine was fatigue- (0.039) related, creatine was food intake- (−0.177) and body weight- (−0.039) related, and both metabolites were clinical score- (0.093; 0.050) and paw edema- (0.125; 0.026) related. Therefore, muscle metabolic alterations were detected in arthritic mice urine, enabling further validation in RA patient’s urine, targeting prognosis, diagnosis, and monitoring of RA-mediated muscle loss.

Graphical Abstract

1. Introduction

Rheumatoid arthritis (RA) is an autoimmune inflammatory disease characterized by symmetric polyarthritis and systemic involvement [1]. It affects about 1% of people under 35 years old and more than 2% of adults over 60 years in the United States, and its prevalence has been reported worldwide [2]. While the joints are the main target of the disease, there are many extra-articular manifestations, such as body composition changes, which are strongly associated with long-term disability and premature mortality [3]. In fact, significant loss of muscle mass is widely described in the literature as one of the major metabolism-related changes seen in RA [4].
As a common and important complication of RA, the loss of muscle mass has been associated with the intensity of inflammation and the severity of disease [5]. Low muscle quantity, quality, and strength are the main features of sarcopenia, a muscle disease (muscle failure) associated with several other conditions [6]. In RA, the presence of sarcopenia has been reported in around 28–37% of patients [7,8,9]. Moreover, the RA muscle loss may occur with or without loss of fat mass, resulting in limited or absent changes in body mass index (BMI), characterizing rheumatoid cachexia (RC) [10]. The RC definition contrasts with classic cachexia, which is characterized by severe weight loss, mainly loss of muscle mass, and increased protein catabolism due to an underlying illness [11]. In a systematic review with meta-analysis, the prevalence of RC was 15–32%, while classical cachexia was absent [12].
The precise assessment of skeletal muscle status and its changes over time is the main problem to diagnose and treat body composition changes. Currently, there are several ways to estimate muscle loss, as can be listed: anthropometrics (e.g., BMI), bioelectrical impedance analysis, imaging techniques (computer tomography, nuclear resonance imaging, dual X-ray absorptiometry, ultrasound), biochemical analysis (assessment of total potassium), urinary creatinine amount, and different tests to quantify muscle function. However, they either lack diagnostic value (sensitivity, specificity) or are high in cost [13,14]. In this context, the search for novel biomarkers related to muscle loss is important to provide the attending physician with better ways to predict development, stage, and progression of muscle involvement during the routine follow-up of RA patients.
Metabolomics is a systems biology technique that uses a ‘top-down’ approach in which data is gathered at the systemic level [15]. The metabolome is defined as the complete set of metabolites from 100 to 1000 Da present in each biological system and currently represents a strong new tool for research [16,17]. Metabolic footprinting comprises the analysis of the extracellular metabolites produced by the organism and can be used to discriminate differences between health and diseased states and to investigate promising biomarkers for diseases [18].
As a comprehensive and sensible technique, the metabolomics analysis is reliable and reproducible and has been used in RA in a comprehensive way, e.g., for disease characterization and treatment response prediction [15,16,18,19]. Considering that the collagen-induced arthritis (CIA) model shares many similarities with RA, the purpose of this study was to analyze the urinary metabolomic profile associated with loss of muscle mass in CIA animals over time.

2. Materials and Methods

This study demonstrates the in vivo data of disease induction, animal follow-up and urine collection, as well as metabolomic analysis. Detailed in vivo experimentation methods were published and are described in the Appendix A in the Journal of Cachexia, Sarcopenia and Muscle Wasting. As all CIA mice presented muscle loss following the disease induction, in the present study the urine metabolomic profile of these mice was analyzed [20].
Animals: Male DBA/1J mice from 8 to 12 weeks of age were used. The mice were reared alone at 20 °C, with 12-h light–dark cycle, and free access to food and water. The animals were randomly divided into two experimental groups: (i) healthy animals (CO, n = 11) and (ii) collagen-induced arthritic animals (CIA, n = 13). Animals were followed up for 65 days and all measurements were done prior to the arthritis induction and thereafter at days 18, 25, 35, 45, 55, and 65. All experiments were performed in accordance with the Guiding Principles for Research Involving Animals and in accordance with the Ethics Committee of Research and Postgraduate Group of the Clinical Hospital of Porto Alegre (Hospital de Clínicas de Porto Alegre—HCPA; number. 14-0297).
Disease induction: Arthritis was induced with bovine type II collagen (CII, Chondrex, Inc., Woodinville, USA; 2 mg/mL) dissolved in 0.1 M acetic acid at 4 °C for 12 h, and Complete Freund’s Adjuvant (CFA; Sigma, St. Louis, MO, USA; 2 mg/mL) containing inactivated Mycobacterium tuberculosis. Fifty microliters of emulsion (CII + CFA) were intradermally injected at the base of the tail to induce arthritis; it was set as the day zero in this experiment. Eighteen days after the first injection, the animals received a reinforcement of CII emulsified with incomplete Freund’s adjuvant (IFA—without Mycobacterium tuberculosis) in another site of the tail (booster injection). During the procedures, mice were anesthetized with isoflurane 10% (Abbott) and 90% of oxygen [21]. Healthy control mice were manipulated and anesthetized; however, no injection was made. Animals were euthanized 65 days after the first injection.
Animal models follow up: The following experiments were performed and the results, which had already been published by the group, were used for correlation with the metabolites in the present study. At days 18, 25, 35, 45, 55, and 65, clinical severity score, hind paw edema, body weight, food intake, free exploratory locomotion, grip strength, and endurance exercise performance of the mice were evaluated. At the end of the experimental period, animals were euthanized, and the tibialis anterior and gastrocnemius muscles were dissected and weighed to confirm cachexia development. The tibiotarsal joint was collected to confirm the development of arthritis by histological analysis with Hematoxylin-Eosin staining [20].
Urine collection: Animals were confined individually in metabolic cages (Tecniplast S.p.A., Buguggiate, Italy) for 6 h to collect urine before the arthritis induction and thereafter at days 18, 25, 35, 45, 55, and 65. Urine was filtered in 0.22-µm filters and stored in −80 °C until further nuclear magnetic resonance spectroscopy (NMR) analysis.
NMR sample preparation and analysis: Urine samples were thawed on ice and centrifuged for 5 min, at 13,000× g at 4 °C. Supernatants were collected and 4× NMR buffer (final concentration: 100 mM of phosphate, 2 mM of difluorotrimethylsilanylphosphonic acid (DFTMP, Manchester Organics, Manchester, UK), 10% D2O, 0.1% azide, 0.5 mM of 4,4-dimethyl-4-silapentane-1-sulfonic acid) was added and mixed. Then samples were centrifuged for 5 min, at 13,000× g at 4 °C, and supernatants were transferred into glass champagne vials, capped and frozen at −80 °C until the NMR analysis. Samples were thawed and transferred to 1.7 mm NMR tubes (Bruker Biospin, Coventry, UK) using an Anachem Autosampler. One-dimensional 1H spectra were acquired at 300 K using a standard 1D-1H-Nuclear Overhauser Effect spectroscopy (NOESY) pulse sequence with water saturation using pre-sat in a Bruker AVANCE II 600 MHz NMR spectrometer (Bruker Corp., Billerica, USA) equipped with a 1.7-mm cryoprobe. Spectral width was set to 12 ppm, and the scans were repeated 128 times. Sample’s series were loaded into 96-tube racks and held at 6 °C in the SampleJet sample handling device until processed. Two-dimensional 1H J-resolved (JRES) spectra were also acquired to aid metabolite identification.
NMR spectra metabolic identification and pathway analysis: Chenomx software (Chenomx Inc., Edmonton, Canada) and the associated metabolites library were used to identify the compounds in 1D NMR spectra. The 2D-Jres NMR spectra were submitted to Birmingham Metabolome Library (BML) to identify and quantify their intensities. MetaboAnalyst website (http://www.metaboanalyst.ca/, 1 July 2021) [22] was used to elucidate which pathways the most relevant metabolites were related to [16,19,20].
Metabolomic network: Cytoscape analysis (the Cytoscape 3.3; cytoscape.org) was used to visualize and to analyze the biological networks between metabolic profiles and clinical parameters of the animal model follow-up. Data were imported in Metascape, a plugin of Cytoscape, and subjected to pathway analysis. All the metabolic pathways generated were further subjected to enrichment and topological analysis. Only the pathways having significant p-value were selected and further subjected to identification of biological processes associated. Compounds, reactions, enzymes, genes, pathways, and the relationships among them provided the initial framework for the metabolomic data analysis [23].
Statistics: Sample size was based on the previous research of our group in which the main outcome was muscle atrophy accessed by myofiber area in CIA model [24]. Considering an alpha of 5% and the statistical power of 90% to detect differences in tibialis anterior myofiber area, a sample size of seven animals per group was required. Due to the possibility of death of about 10% of the animals during the disease induction procedure, or during the disease course, the final sample size was 11 mice for the CO group and 13 mice for the CIA group. 1D NMR spectra were binned, normalized, glog-transformed, and subjected to principal component analysis (PCA) and partial least-squares discriminant analysis (PLSDA) to group data, followed by a venetian blind test to obtain specificity and sensibility of the model. PCA and PLSDA were also used to analyze BML results. Partial least-squares regression analysis (PLSR) was performed using the BML metabolic profile with each clinical data, as previously described [18]. Briefly, PLSR is a regression method that identifies which metabolites can predict a continuous variable (in vivo analysis was considered a continuous variable). This analysis yields R2, a measure of the cross-validated goodness-of-fit of the linear regression, and each metabolite contribution to the model, while permutation testing (multiple analyses using random data subsets) was used to assess the significance of this prediction. Pathway analysis was performed applying the Fisher’s exact test by the MetaboAnalyst website and to reach pathway impact [25,26,27]. Statistical significance was always set for a p-value under 0.05. The PLSDA model of the variable of importance in the projection score critical limit was set for a value higher than 1 × 10°, which means a p-value under 0.05. Metabolite list picking was manually done through KEGG (www.genome.jp/kegg, 1 July 2021) [28], PUBCHEM (pubchem.ncbi.nlm.nih.gov) [29], and HMDB (www.hmdb.ca, 1 July 2021) [30] websites in which it was checked if there was a direct relationship with muscle metabolism.

3. Results

In vivo experimental: CIA animals had significantly higher arthritis scores and hind paw edema volumes than CO. In addition, from histopathology, the disease was confirmed in all CIA animals. The CIA group showed no weight loss or decreased appetite. However, the normalized weights (i.e., divided by the body weight of the animal) for visceral and brown fat were reduced in CIA animals (54% and 39%, respectively) and ankle joint normalized weights and spleen normalized weights were increased in the CIA group when compared to CO (18% and 40%, respectively). CIA animals showed lower free exploratory locomotion, endurance exercise performance total time, and grip strength compared to the CO group. The dissected tibialis anterior and gastrocnemius muscles weighed less in the CIA than in the CO group (25% and 24%, respectively) and sarcoplasmic ratios were also smaller in CIA when compared with CO (23% and 22% less sarcoplasmic ratio, respectively). Thus, the myofiber diameter reduction was 45% in TA and 41% in GA [20]. After the confirmation of muscle loss in CIA mice, the non-targeted metabolome was acquired to identify all metabolites present in the urine of these animals during the development of arthritis [20]. To search for potential biomarkers of muscle loss in cachectic mice, urine samples from these same cohort of animals were collected for metabolic profiling analysis at seven time points (0, 18, 25, 35, 45, 55 and 65 days) during the 65 days of disease. Profile comparisons between CIA and CO mice, and between each time point within the CIA group, were performed aiming to characterize metabolic profiles of the diseased animals and their changes over time.
PCA and PLSDA models: Firstly, we developed an unsupervised PCA model using all the NMR processed spectra and it was followed by PLSDA supervised analysis to group them into clusters for group homogeneity (Figure 1A,B). The main component and latent variable of each sample was plotted with mean ± standard deviation of the mean (Figure 1C,D).
To reduce complexity, we performed a PCA model followed by PLSDA at each time point comparing CIA and CO (paired analysis; Figure 2 and Figure 3, PCA and PLSDA, respectively). In addition, to elucidate the disease course within CIA animals, we conducted paired analyses between the different time points (Figure 4 and Figure 5) in the CIA group. As can be seen, these PLSDA models could segregate CO from CIA groups as well as CIA groups in a time-dependent manner. Lastly, all the previously mentioned models from PLSDA were submitted to a venetian blind test to provide sensitivity and specificity values to be considered, as shown in Table 1. As observed by PCA, method samples have been homogeneously grouped with almost no samples outside the confidence interval (CI; Figure 2 and Figure 4). Based on the venetian blind test done after PLSDA, we decided to maintain all samples that could be considered outliers because those samples added more specificity and sensibility to the model than when removed (Figure 3 and Figure 5). Taken together, the built models were robust enough to provide a long list of statistical significance (VIP score higher than 1 × 10°) metabolites that were analyzed for muscle loss process.
Metabolic profiling: Considering the whole set of urine samples, almost 400 metabolites were identified, and time-dependent significant differences between CIA vs. CO, and CIA vs. CIA were found (Figure 6). From the list of all metabolites with statistically significant differences (data not shown), we speculate that these differences were related to the arthritis disease and to the processes involved in its development, which may have effects in all organs and tissues, in addition to the muscle loss process.
From all metabolites, 28 metabolites were associated with pathways related to muscle tissues, including muscle catabolic and anabolic processes based on information from KEGG, PUBCHEM, and HMDB websites (Table 2). The following metabolites showed high importance in the models: carnosine (VIP score 2.8 × 102) from CIA vs. CO at 65 days; histidine (VIP score 1.2 × 102) from CIA pair analysis days 55 vs. 65; histamine (VIP score 1.1 × 102) from CIA pair analysis days 55 vs. 65; L-methionine (VIP score 1.1 × 102) from CIA pair analysis days 0 vs. 18; and succinyl acetone (VIP score 1.0 × 102) from CIA vs. CO at 65 days. As expected, given their importance in muscle protein structure, there were also several amino acids with high statistical differences between the compared groups, i.e., high variable importance in the projection score.
Biomarker time-dependent identification:
After assessing the statistically significant differences in metabolites levels, these were then linked to pathways known to occur within muscle tissues. These metabolites were displayed in a timeline, according to the analysis results of different time points (Figure 7a). Before 45 days of established muscle loss presence, we identified the 3-methylhistidine, 4-aminobutyric acid, acetylcholine, arginine, aspartate, carnosine, creatine, creatinine, glutamine, histamine, histidine, isoleucine, leucine, l-methionine, lysine, myo-inositol, N,N-dimethylglycine, N-acetyl alanine, N-acetylmethionine, pantothenate, phenylalanine, phosphocholine, phosphocreatine, pyridoxine, sarcosine, succinyl acetone, and thiamine. After 45 days, from the previous list, aspartate, phenylalanine, and thiamine were excluded for the purpose of follow up. On the other hand, urocanate was included. The metabolites elected also are shown in a heatmap (Figure 7b) as well as a connection network with the weight each connection has (Figure 7c).
Several metabolites were related to the same 18 pathways that occur within muscle tissues, and 11 pathways were related to amino acid metabolism (Figure 8). Furthermore, the four most frequent pathways were: histidine metabolism (e.g., metabolites: 3-methylhistidine, carnosine, histamine, histidine, and urocanate); arginine and proline metabolism (arginine, 4-aminobutyrate, creatine, and phosphocreatine); glycine, serine, and threonine metabolism (N,N-dimethylglycine, sarcosine, phosphocholine); and creatine phosphate metabolism (creatine, creatinine, phosphocreatine). Some of these pathways also presented statistical differences between several time points, as shown in Table 3.
Briefly, histidine metabolism was the most frequent pathway altered during disease development compared to controls. Concerning metabolite frequency within the same metabolic pathway, histidine had the highest impact of 0.60 comparing the days 55 and 65 in CIA. Arginine and proline metabolisms also appeared frequently, with an impact of up to 0.28; valine, leucine, and isoleucine metabolisms appeared in almost every time point comparing CIA and CO, and had an impact of up to 0.33.
Regression analysis: Finally, we performed a regression analysis between the clinical data and metabolites from urine samples, and identified them by BML (Table 4). Although the regressions models generated with these data presented statistically significant differences, the highest R2 found in the models were weak, i.e., from 0.04 to 0.3. PLSR for grip strength showed that with 24 metabolites (including N-dimethylglycine, thiamine, and arginine) the model reached a maximum contribution R2 of 0.30 (p < 0.01); endurance exercise performance showed that with 50 metabolites (including N-acetyl alanine, sarcosine, N-acetylmethionine, phosphocholine, carnosine, and glutamine) the model reached a maximum contribution R2 of 0.11 (p < 0.17); free locomotion required 40 metabolites (including thiamine, N-acetyl alanine, sarcosine, arginine, and carnosine) to reach a maximum contribution R2 of 0.04 (p < 0.03); clinical score required 26 metabolites (including creatinine, methionine, carnosine, leucine, and glutamine) to reach a maximum contribution R2 of 0.05 (p < 0.03); hind paw edema showed that with 12 metabolites (including carnosine, leucine, creatinine, and 4-aminobutyrate) the model reached a maximum contribution R2 of 0.084 (p < 0.08); food intake required 15 metabolites (including N-acetyl alanine, creatinine, and 3-methylhistidine) to reach a maximum contribution R2 of 0.30 (p < 0.08); and body weight required 10 metabolites (including creatinine) to reach maximum contribution R2 of 0.06 (p < 0.01). Table 4 shows the R2 contribution value from the PLSR model for the statistically significant correlations between the 22 muscle-related metabolites and the clinical data. The metabolites R2 not related to muscle metabolism are not shown. Additionally, the 3-methylhistidine, 4-aminobutyric acid, acetylcholine, arginine, aspartate, carnosine, creatine, creatinine, glutamine, histamine, histidine, isoleucine, leucine, l-methionine, lysine, myo-inositol, N,N-dimethylglycine, N-acetyl alanine, N-acetylmethionine, pantothenate, phenylalanine, phosphocholine, phosphocreatine, pyridoxine, sarcosine, succinyl acetone, and thiamine were associated with muscle loss processes (Figure 8).

4. Discussion

Several metabolites and pathways related to muscle catabolism and anabolism were altered in urine of CIA mice compared to controls and relatively to the stage of disease development. Namely, the metabolites 3-methylhistamine, 4-aminobutyric acid, acetylcholine, arginine, carnosine, creatine, creatinine, glutamine, histamine, histidine, isoleucine, leucine, l-methionine, lysine, myo-inositol, N,N-dimethylglycine, N-acetyl alanine, N-acetylmethionine, pantothenate, phosphocholine, phosphocreatine, pyridoxine, sarcosine, succinyl acetone, and urocanate, which are involved in important signaling pathways for muscle tissues, showed more impact following the metabolomic analysis. Although muscle loss is an important feature in RA and in experimental arthritis, the number of studies is limited, and the correlation of muscle impairment in human and animal models is unclear [31]. Thus, the altered profiles found in this study may be valuable for the clinical assessment of RA muscle loss.
As observed in the Venn diagram (Figure 6), the analysis brought a list of 54 unique identified metabolites by the BML, 265 uniquely identified by ChenomX-NMR, and 154 common metabolites found in both automatic (BML) and manual (ChenomX—NMR) assessments. With all these metabolites, we could access an initial list of nearly 100 metabolites after the first steps of statistical analysis (PCA, PLSDA, and Venetian blind test). Afterwards, using the pathway analysis, and confirming the metabolites list in the KEGG, PUBCHEM, and HMDB databases, and performing the PLS-regression, we finally could reduce this list down to 27 metabolites previously mentioned. All these steps favor building a robust statistical set to reach the most muscle-loss linked list of metabolites in this research.
As previously demonstrated by our group, DBA/1J mice with CIA present approximately 30% decrease in myofiber cross-sectional area after 45 days of disease induction and are considered cachectic at the sixty-fifth day [20,24]. Continuing the investigation about arthritis muscle loss, we performed metabolomics analysis to explore potential biomarkers in the urine of these cachectic mice. The use of animal models has the advantage of circumventing the great heterogeneity implicit in human studies, increasing the power to prospect useful biomarkers. Previous metabolomics studies in RA have addressed profiles for disease activity, drug exposure, and prognosis in early disease, demonstrating the feasibility of this approach [15,17,18,19,32]. A systematic review compiled metabolomics studies, which used serum, urine, and synovial fluid for the analysis, and designated the following metabolites as important in RA: glucose, lactic acid, citric acid, leucine, methionine, isoleucine, valine, phenylalanine, threonine, serine, proline, glutamate, histidine, alanine, cholesterol, glycerol, and ribose [33]. Furthermore, studies reported that muscle features can also be measured using a metabolomics approach. Reduced amine metabolites have been found in plasma of CIA mice indicating that disordered amine response may be linked to the muscle wasting and to the increased resting energy expenditure mediated by RA [34,35]. The analysis of urine collected from CIA rats treated with tetrandrine described 23 potential biomarkers associated with CIA, which were mainly linked to metabolism of energy, amino acid, lipid, and gut microbe [36]. These findings in RA are similar to ours, demonstrating that metabolomics may be an interesting approach not only for disease prognosis and monitoring, but also to assess the muscle loss related to disease pathogenesis.
The histidine metabolism was the pathway most frequently altered in CIA mice. Furthermore, all metabolites related to this pathway (3-methylhistidine, carnosine, histamine, histidine, and urocanate) had major statistical significance in the PLSDA model analysis between CIA and CO groups. Histidine was described as one of the amino acids’ end products, highly excreted by RA patients, and as a discriminatory metabolite for RA treatment prediction [19]. Compared to controls, histidine serum levels were found lower in RA patients at baseline, and significantly higher 3 months after the treatment with TNF-α inhibitors [37]. Additionally, both RA and experimental arthritis showed better outcomes following the supplementation with carnosine [38,39]. Concerning the regression analysis, carnosine was related to fatigue, free locomotion, clinical score, and hind paw edema best models, while 3-methylhistidine was related to food intake. Accordingly, both metabolite 3-methyl-histidine and carnosine were identified as markers of muscle loss, and oral administration of histidine was able to improve grip strength and walking speed in chronic kidney disease patients [40,41]. In addition, lower levels of carnosine dipeptidase 1were associated with cancer cachexia, compared with weight stable patients [42]. Therefore, as histidine pathway has a strong relationship with muscle deficits, the metabolites of this pathway may potentially be validated as biomarkers for muscle loss in RA.
Arginine and proline metabolisms were also affected by CIA, and specifically the metabolites arginine, 4-aminobutyrate, creatine, and phosphocreatine had high statistical significance in PLSDA models between CIA and CO groups. Previously, metabolomics profiles of RA synovial fluid showed that arginine and proline pathways are downregulated with the disease, while serum levels of L-proline were increased in RA patients compared to controls [43,44]. Additionally, the metabolite 4-aminobutyrate, also known as GABA, was suggested to have influence in RA due to a regulatory role on inflammation [45]. The influence of the semi-essential amino acid, arginine, on muscle tissues, is related to the production of creatine, urea, and nitric oxide (NO), and to the synthesis of new proteins [46]. Through its vasodilatory function, NO may play a role in nutrient delivery to the muscle, as well as in healing or fibrosis fate [46,47,48,49,50]. Regarding muscle function, arginine was related to grip strength and free locomotion in the best regression models. Accordingly, it has been demonstrated that serum arginine was suitable to differentiate RA patients into classes of physical disability [51]. Thus, we can suggest a role of arginine and proline metabolisms in RA muscle loss process.
As creatine phosphate metabolism is deeply important for muscle function it was not surprising that metabolites related to this pathway were observed several times in our data. Creatinine levels are strongly associated with muscle energetic metabolism and wasting and are increased in individuals with muscle loss [52]. Creatinine is the degradation product of creatine when it is properly phosphorylated when used in the muscle energy processes. Although creatine supplementation in RA patients increases muscle mass, it does not change strength or physical function [53]. This metabolite has appeared in the best regression models for clinical score, hind paw edema, food intake, and body weight. As formerly mentioned, this pathway was used to monitor muscle mass change, being a strong candidate for validation as a biomarker also in RA.
Glycine, serine, and threonine compose metabolic pathways related to amino acids, and the metabolites N,N-dimethylglycine, sarcosine, and phosphocholine were frequent in PLSDA models. N,N-dimethylglycine has been positively correlated with fat-free mass in chronic obstructive pulmonary disease patients [40]. Otherwise, increased phosphocholine levels may have indirect effects on muscle loss, since murine cancer cells from tumors, which are capable of inducing cachexia, presented increased phosphocholine levels compared to tumors unable to induce cachexia [54]. Regarding muscle function, N,N-dimethylglycine contributed to the best model of grip strength, which agrees with the potential role of glycine related to fatigue previously described [55]. In RA, serum glycine was suitable to differentiate RA patients into classes of physical disability [51]. Therefore, directly, or indirectly, this pathway also characterizes a good biomarker candidate for muscle loss in RA.
In addition to pathways discussed before, the influence of other amino acids within muscle metabolism have been reported. Glutamine, in which muscle levels are known to be an important supply of this amino acid to the body, appeared in both endurance exercise performance and clinical score best regression models [56]. The availability of glutamine is reduced by conditions such as inactivity, and its low muscle concentration is associated with decreased protein synthesis in acute disease states [56,57]. Regarding methionine, which appeared in the best regression model for clinical score, it is the primary amino acid required to initiate protein synthesis, as well as acute phase proteins synthesis [58]. Leucine contributed to the best PLSR model for edema; this amino acid is suggested as a stimulator of protein synthesis, with potential to reduce muscle loss in mice with adenocarcinoma after supplementation [59,60]. Notwithstanding, metabolites not related to the amino acid metabolism, such as B vitamin complex components, also have anabolic roles in skeletal muscle tissue [61]. The best regression models for grip strength and free locomotion were related to thiamine.
Some statistically significant metabolites found in PLSDA in our model were not extensively discussed because of the paucity of literature regarding muscle metabolism process and/or related to chronic inflammatory process. These metabolites are listed here as follow: urocanate (histidine metabolism), succinyl acetone (tyrosine metabolite), N-acetylmethionine (methionine derivative), N-acetyl alanine (substrate from a variety of cellular reactions), acetylcholine (a neurotransmitter at neuromuscular junction), sarcosine (a derivative of glycine), and aspartate (a non-essential proteinogenic amino acid).
The present study has some limitations: lack of kidney function data or proteinuria to confirm any interference to the metabolic profile. Another limitation is that BML is an automatic detection system that cannot be controlled, resulting in a lack of detection of some metabolites, and leading to a less precise analysis; however, performing the NMR analysis functioned as a complementary analysis, and confirmed the automated BML findings. Despite this, our findings should not be neglected.

5. Conclusions

In conclusion, several of the metabolites found in the urine of CIA mice with muscle loss were related to muscle metabolism. Obtaining these results in an animal model highlights the possibility of relevant findings following the analysis of urine metabolites in RA patients. Of note, our group is already performing research that aims to investigate the urine metabolic profiles of RA subjects, and relate them to RC and sarcopenia phenotypes, in addition to analyzing the association of these conditions with the disease activity in patients (data not yet published). The addressed issue is important, and if confirmed in human studies, may be important for early individuation of subjects at risk for muscle wasting in RA. The identification and quantification of metabolites in RA patients would be useful to predict and to monitor the muscle involvement of the disease, since the alterations in pathways may differ among patients with diverse genetic and phenotypic backgrounds. Additionally, muscle loss is not limited to RA, but occurs in other several chronic inflammatory diseases, and in the aging process as sarcopenia. Thus, the metabolites could be validated as biomarkers for these diseases as well.

Author Contributions

P.V.G.A.: planning, cachexia model follow-up and metabolomics experiments, statistics, scientific discussion, and writing; J.M.S.S.: planning, cachexia model follow-up, statistics, scientific discussion, and writing; R.C.E.S.: planning, cachexia model follow-up, and scientific discussion; M.S.O.: metabolomics pathways analysis, scientific discussion, and writing; A.S.A.: cachexia model follow-up and writing; M.S.d.O.: cachexia model follow-up and writing; M.L.I.: cachexia model follow-up; E.C.F.: planning, cachexia model follow-up, and scientific discussion; V.O.N.T.: planning, statistics, scientific discussion, and writing; C.L.B.: metabolomics, and statistics; C.V.B.: paper corrections; L.I.F.: original idea of cachexia model; S.P.Y.: supervision of metabolomics steps, planning, scientific discussion, and paper corrections; P.S.L.; planning, cachexia model follow-up, statistics, scientific discussion, and writing; R.M.X.: supervision of cachexia model experiments, planning, scientific discussion, and paper corrections. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by grants from Coordination for the Improvement of Higher Education Personnel (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—CAPES), Foundation of Research Support of the State of Rio Grande do Sul (Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul—FAPERGS), Funding for Research and Events of Clinical Hospital of Porto Alegre (Fundo de Incentivo à Pesquisa e Eventos do Hospital de Clínicas de Porto Alegre—FIPE-HCPA), and National Council for Scientific and Technological Development (Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq).

Institutional Review Board Statement

This study was performed in accordance with the Guiding Principles for Research Involving Animals and in accordance with the Ethics Committee of Research and Postgraduate Group of the Clinical Hospital of Porto Alegre (Hospital de Clínicas de Porto Alegre—HCPA; number: 14-0297).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the animal facility professionals for their help during animal experimentation—Unidade de Experimentação Animal do Hospital de Clínicas de Porto Alegre—and the University of Birmingham for supporting NMR procedures.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Clinical severity score and measurement of edema: Arthritis severity was clinically determined for each paw, three times a week, according to the a scale of 0 to 4 (0, no evidence of erythema and swelling; 1, erythema and mild swelling confined to the tarsals or metatarsals; 2, erythema and moderate swelling of tarsal and the metatarsal or tarsal and ankle joints; 3, erythema and severe swelling extending from the ankle to metatarsal joints; and 4, erythema and severs welling encompassing the ankle, foot and digits, or ankylosis of the limb). The highest sum score that a mouse could reach was 16. Hind paw edema volume was measured using a plethysmometer (Insight Ltd.a., Ribeirão Preto, Brazil). Briefly, it is a small cylinder filled with a buffer connected to a device capable to measure the total fluid volume. We immersed the hind paw of the animal inside the cylinder, then added the total volume added was measured, and the difference between the final volume minus the initial volume resulted in pawn total volume.
Animal weight and food intake: Animals were weighed for total body mass three times a week starting before the disease induction. Animal food was weighed three times a week starting before the disease induction. The remaining food was subtracted from the offered food to follow up the food intake of each animal. Mean food intake was calculated through the sum of the total feed period between each interval (0, 18, 25,35, 45, 55, and 65 days after the first injection) and divided by the total sampling time of each period of time.
Free exploratory locomotion: Animals were put in an acrylic cage with motion sensors (Monitor de Atividades IR; Insight Equipaments Ltd.a., Ribeirão Preto, Brazil). After 30 s of adaptation, all free exploratory locomotion was recorded for 5 min. The total distance covered by each animal was sent to a computer for posterior analysis.
Grip strength: Animals were tested for maximum grip strength with a test adapted from Deacon. Briefly, first we used mesh with proper loads, each one amounting to 5, 20, 35, 50, 65, 80, and 95 g. Each mouse was held by the first third of the tail and suspended until it grasped the lighter weight with all paws. The animal had to hold the load for at least 3 recorded seconds. If the animal succeeded, it rested for 30 s before trying the next weight. If the animal failed three times with a 10-s rest between each attempt, the longest time it was able to hold the weight was recorded. The following equation was used: Fmax = P3seg + (5 × t < 3 seg), where Fmax is the maximum calculated grip strength, P3seg is the heaviest load the animal held for 3 s, and t < 3 seg is the longest time the animal held the heaviest load. The final result was expressed in grams (g).
Endurance exercise performance: The endurance exercise performance test was conducted at the beginning of the dark cycle of the mice, when animals are active. The endurance exercise performance test consists in inducing mice into exhaustion. The mice were separately set to run on treadmill developed by Hospital de Clínicas de Porto Alegre engineers. We followed the previous protocol. Briefly, animals were left on the treadmill for 5 min to adapt prior to the test. After that, the test started at a speed of 8.5 m·min−1 for 9 min, for the animals to warm up, and then every 3 min, the speed was increased in a constant ratio of 2.5 m·min−1. Animals reach exhaustion when they are not able to continue running up to maintain the speed (in this case, when the mice stayed for 10 s at the first third of the rear part of the equipment). The exhausted animals were put back into their cages, and the following measurements were recorded: speed at fatigue, total experiment time, and total covered distance. Only the total experiment time was used in this research.
Organs and tissue dissection: Sixty-five days after the first injection, the tibialis anterior and gastrocnemius muscles were dissected immediately after euthanasia, then weighed, and collected to measure muscle myofiber diameter by histological analysis with hematoxylin–eosin (HE) staining. The tibiotarsal joint was collected to confirm the development of arthritis by histological analysis with HE staining. The following organs were dissected: heart, liver, spleen, quadriceps, soleus muscle, and visceral and brown fat. All dissected organs and tissues not used for histological analyses were weighed in precision balance and frozen under 80 °C.
Histological analysis: The tibiotarsal joint, tibialis anterior, and gastrocnemius muscle of the DBA/1J animals were dissected and immersed in 10% buffered formalin for fixation for up to 3 days. Next, the tibiotarsal joints were decalcified in 10% nitric acid for 24 h. All these tissues were dehydrated and embedded in paraffin blocks. Slices 6-μm thick were arranged on microscope slides. Tibialis anterior and gastrocnemius muscle slices stained with HE were used for myofiber diameter measurement. One transverse section of each muscle was stained with HE and analyzed under an optic microscope (×400). Two straight lines crossing at a right angle at the fiber center were drawn in each myofiber. The mean of these diameters (in micrometers) was used to calculate the transverse section mean, based on circle area. For measuring the myofiber diameter of the whole muscle, we took 10 pictures of each section, and 20 fibers were measured from each picture using the Image-Pro Express software (version5.1.0.12, Media Cybernetics, Rockville, MD, USA). We used a histological score system to evaluate individual joints and assess arthritis severity. For synovial inflammation, five high-power magnification fields were scored for the percentage of infiltrating mononuclear cells as follows: 0, absent; 1, mild (1–10%); 2, moderate (11–50%); 3, severe (51–100%); for synovial hyperplasia: 0, absent; 1, mild (5–10 layers); 2,moderate (11–50 layers); 3, severe (>20 layers); for extension of pannus formation based on the reader’s impression:0, absent; 1, mild; 2, moderate; 3, severe; for synovial fibrosis: 0, absent; 1, mild (1–10%); 2, moderate (11–50%); 3, severe (51–100%); for cartilage erosion, that is, the percentage of the cartilage surface that was eroded: 0, absent; 1, mild (1–10%); 2, moderate (11–50%); 3, severe (51–100%); and for bone erosion: 0, none; 1, minor erosion(s)observed only at high-power magnification fields; 2, moderate erosion(s) observed at low magnification; and 3, severe transcortical erosion(s).

References

  1. Smolen, J.S.; Aletaha, D.; Barton, A.; Burmester, G.R.; Emery, P.; Firestein, G.S.; Kavanaugh, A.; McInnes, I.B.; Solomon, D.H.; Strand, V.; et al. Rheumatoid arthritis. Nat. Rev. Dis. Primers 2018, 4, 18001. [Google Scholar] [CrossRef]
  2. Navarro-Millan, I.; Singh, J.A.; Curtis, J.R. Systematic Review of Tocilizumab for Rheumatoid Arthritis: A New Biologic Agent Targeting the Interleukin-6 Receptor. Clin. Ther. 2012, 34, 788–802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Giles, J.T. Extra-articular Manifestations and Comorbidity in Rheumatoid Arthritis: Potential Impact of Pre-Rheumatoid Arthritis Prevention. Clin. Ther. 2019, 41, 1246–1255. [Google Scholar] [CrossRef] [PubMed]
  4. Masuko, K. Rheumatoid cachexia revisited: A metabolic co-morbidity in rheumatoid arthritis. Front. Nutr. 2014, 1, 20. [Google Scholar] [CrossRef] [Green Version]
  5. Uutela, T.I.; Kautiainen, H.J.; Häkkinen, A.H. Decreasing muscle performance associated with increasing disease activity in patients with rheumatoid arthritis. PLoS ONE 2018, 13, e0194917. [Google Scholar] [CrossRef]
  6. Cruz-Jentoft, A.J.; Bahat, G.; Bauer, J.; Boirie, Y.; Bruyère, O.; Cederholm, T.; Cederholm, T.; Cooper, C.; Landi, F.; Rolland, Y.; et al. Sarcopenia: Revised European consensus on definition and diagnosis. Age Ageing 2019, 48, 16–31. [Google Scholar] [CrossRef] [Green Version]
  7. Mochizuki, T.; Yano, K.; Ikari, K.; Okazaki, K. Sarcopenia-associated factors in Japanese patients with rheumatoid arthritis: A cross-sectional study. Geriatr. Gerontol. Int. 2019, 19, 907–912. [Google Scholar] [CrossRef] [PubMed]
  8. Torii, M.; Hashimoto, M.; Hanai, A.; Fujii, T.; Furu, M.; Ito, H.; Uozumi, R.; Hamaguchi, M.; Terao, C.; Yamamoto, W.; et al. Prevalence and factors associated with sarcopenia in patients with rheumatoid arthritis. Mod. Rheumatol. 2019, 29, 589–595. [Google Scholar] [CrossRef] [PubMed]
  9. Yamada, Y.; Tada, M.; Mandai, K.; Hidaka, N.; Inui, K.; Nakamura, H. Glucocorticoid use is an independent risk factor for developing sarcopenia in patients with rheumatoid arthritis: From the CHIKARA study. Clin. Rheumatol. 2020, 39, 1757–1764. [Google Scholar] [CrossRef]
  10. Roubenoff, R. Rheumatoid cachexia: A complication of rheumatoid arthritis moves into the 21st century. Arthritis Res. Ther. 2009, 11, 108. [Google Scholar] [CrossRef] [Green Version]
  11. Evans, W.J.; Morley, J.E.; Argilés, J.; Bales, C.; Baracos, V.; Guttridge, D.; Jatoi, A.; Kalantar-Zadeh, K.; Lochs, H.; Mantovani, G.; et al. Cachexia: A new definition. Clin. Nutr. 2008, 27, 793–799. [Google Scholar] [CrossRef]
  12. Santo, R.C.; Fernandes, K.Z.; Lora, P.S.; Filippin, L.I.; Xavier, R.M. Prevalence of rheumatoid cachexia in rheumatoid arthritis: A systematic review and meta-analysis. J. Cachexia Sarcopenia Muscle 2018, 9, 816–825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kuriyan, R. Body composition techniques. Indian J. Med. Res. 2018, 148, 648–658. [Google Scholar] [CrossRef]
  14. Ali, S.; Garcia, J.M. Sarcopenia, cachexia and aging: Diagnosis, mechanisms and therapeutic options—A mini-review. Gerontology 2014, 60, 294–305. [Google Scholar] [CrossRef] [Green Version]
  15. Fitzpatrick, M.; Young, S.P. Metabolomics—a novel window into inflammatory disease. Swiss Med. Wkly. 2013, 143, w13743. [Google Scholar] [CrossRef] [PubMed]
  16. Semerano, L.; Romeo, P.-H.; Boissier, M.-C. Metabolomics for rheumatic diseases: Has the time come? Ann. Rheum. Dis. 2015, 74, 1325–1326. [Google Scholar] [CrossRef] [Green Version]
  17. Priori, R.; Scrivo, R.; Brandt, J.; Valerio, M.; Casadei, L.; Valesini, G.; Manetti, C. Metabolomics in rheumatic diseases: The potential of an emerging methodology for improved patient diagnosis, prognosis, and treatment efficacy. Autoimmun. Rev. 2013, 12, 1022–1030. [Google Scholar] [CrossRef]
  18. Young, S.P.; Kapoor, S.R.; Viant, M.R.; Byrne, J.J.; Filer, A.; Buckley, C.D.; Kitas, G.D.; Raza, K. The Impact of Inflammation on Metabolomic Profiles in Patients with Arthritis. Arthritis Rheum. 2013, 65, 2015–2023. [Google Scholar] [CrossRef] [Green Version]
  19. Kapoor, S.R.; Filer, A.; Fitzpatrick, M.A.; Fisher, B.A.; Taylor, P.C.; Buckley, C.D.; McInnes, I.B.; Raza, K.; Young, S.P. Metabolic Profiling Predicts Response to Anti-Tumor Necrosis Factor alpha Therapy in Patients With Rheumatoid Arthritis. Arthritis Rheum. 2013, 65, 1448–1456. [Google Scholar] [CrossRef]
  20. Alabarse, P.V.; Lora, P.; Silva, J.M.; Santo, R.C.; Freitas, E.C.; de Oliveira, M.S.; Almeida, A.S.; Immig, M.; Teixeira, V.O.; Filippin, L.I.; et al. Collagen-induced arthritis as an animal model of rheumatoid cachexia. J. Cachexia Sarcopenia Muscle 2018, 9, 603–612. [Google Scholar] [CrossRef]
  21. Oliveira, P.G.; Grespan, R.; Pinto, L.G.; Meurer, L.; Brenol, J.C.; Roesler, R.; Schwartsmann, G.; Cunha, F.Q.; Xavier, R.M. Protective effect of RC-3095, an antagonist of the gastrin-releasing peptide receptor, in experimental arthritis. Arthritis Rheum. 2011, 63, 2956–2965. [Google Scholar] [CrossRef] [PubMed]
  22. MetaboAnalyst. Available online: http://www.metaboanalyst.ca/ (accessed on 1 July 2021).
  23. Gao, J.; Tarcea, V.G.; Karnovsky, A.; Mirel, B.R.; Weymouth, T.E.; Beecher, C.W.; Cavalcoli, J.D.; Athey, B.D.; Omenn, G.S.; Burant, C.F.; et al. Metscape: A Cytoscape plug-in for visualizing and interpreting metabolomic data in the context of human metabolic networks. Bioinformatics 2010, 26, 971–973. [Google Scholar] [CrossRef] [Green Version]
  24. Filippin, L.I.; Teixeira, V.N.; Viacava, P.R.; Lora, P.S.; Xavier, L.L.; Xavier, R.M. Temporal development of muscle atrophy in murine model of arthritis is related to disease severity. J. Cachexia Sarcopenia Muscle 2013, 4, 231–238. [Google Scholar] [CrossRef]
  25. Xia, J.; Psychogios, N.; Young, N.; Wishart, D.S. MetaboAnalyst: A web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009, 37, W652–W660. [Google Scholar] [CrossRef] [Green Version]
  26. Xia, J.; Mandal, R.; Sinelnikov, I.V.; Broadhurst, D.; Wishart, D.S. MetaboAnalyst 2.0—A comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012, 40, W127–W133. [Google Scholar] [CrossRef] [Green Version]
  27. Xia, J.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. MetaboAnalyst 3.0-making metabolomics more meaningful. Nucleic Acids Res. 2015, 43, W251–W257. [Google Scholar] [CrossRef] [Green Version]
  28. KEGG. Available online: www.genome.jp/kegg (accessed on 1 July 2021).
  29. PUBCHEM. Available online: pubchem.ncbi.nlm.nih.gov (accessed on 1 July 2021).
  30. HMDB. Available online: www.hmdb.ca (accessed on 1 July 2021).
  31. Romanick, M.; Thompson, L.V.; Brown-Borg, H.M. Murine models of atrophy, cachexia, and sarcopenia in skeletal muscle. Biochim. Et Biophys. Acta-Mol. Basis Dis. 2013, 1832, 1410–1420. [Google Scholar] [CrossRef] [Green Version]
  32. Sasaki, C.; Hiraishi, T.; Oku, T.; Okuma, K.; Suzumura, K.; Hashimoto, M.; Ito, H.; Aramori, I.; Hirayama, Y. Metabolomic approach to the exploration of biomarkers associated with disease activity in rheumatoid arthritis. PLoS ONE 2019, 14, e0219400. [Google Scholar] [CrossRef] [Green Version]
  33. Li, C.; Chen, B.; Fang, Z.; Leng, Y.-F.; Wang, D.-W.; Chen, F.-Q.; Xu, X.; Sun, Z.-L. Metabolomics in the development and progression of rheumatoid arthritis: A systematic review. Jt. Bone Spine 2020, 87, 425–430. [Google Scholar] [CrossRef]
  34. He, M.; Harms, A.C.; Van Wijk, E.; Wang, M.; Berger, R.; Koval, S.; Hankemeier, T.; Van Der Greef, J. Role of amino acids in rheumatoid arthritis studied by metabolomics. Int. J. Rheum. Dis. 2019, 22, 38–46. [Google Scholar] [CrossRef] [PubMed]
  35. He, M.; van Wijk, E.; van Wietmarschen, H.; Wang, M.; Sun, M.; Koval, S.; van Wijk, R.; Hankemeier, T.; van der Greef, J. Spontaneous ultra-weak photon emission in correlation to inflammatory metabolism and oxidative stress in a mouse model of collagen-induced arthritis. J. Photochem. Photobiol. B: Biol. 2017, 168, 98–106. [Google Scholar] [CrossRef]
  36. Yun, X.; Dong, S.; Hu, Q.; Dai, Y.; Xia, Y. 1H NMR-based metabolomics approach to investigate the urine samples of collagen-induced arthritis rats and the intervention of tetrandrine. J. Pharm. Biomed. Anal. 2018, 154, 302–311. [Google Scholar] [CrossRef]
  37. Zabek, A.; Swierkot, J.; Malak, A.; Zawadzka, I.; Deja, S.; Bogunia-Kubik, K.; Mlynarz, P. Application of 1 H NMR-based serum metabolomic studies for monitoring female patients with rheumatoid arthritis. J. Pharm. Biomed. Anal. 2016, 117, 544–550. [Google Scholar] [CrossRef] [Green Version]
  38. Bauerová, K.; Poništ, S.; Mihalová, D.; Dráfi, F.; Kuncírová, V. Utilization of adjuvant arthritis model for evaluation of new approaches in rheumatoid arthritis therapy focused on regulation of immune processes and oxidative stress. Interdiscip. Toxicol. 2011, 4, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Drafi, F.; Bauerova, K.; Kuncirova, V.; Ponist, S.; Mihalova, D.; Fedorova, T.; Harmatha, J.; Nosál, R. Pharmacological influence on processes of adjuvant arthritis: Effect of the combination of an antioxidant active substance with methotrexate. Interdiscip. Toxicol. 2012, 5, 84–91. [Google Scholar] [CrossRef] [PubMed]
  40. Ubhi, B.K.; Riley, J.H.; Shaw, P.A.; Lomas, D.A.; Tal-Singer, R.; MacNee, W.; Griffin, J.L.; Connor, S.C. Metabolic profiling detects biomarkers of protein degradation in COPD patients. Eur. Respir. J. 2011, 40, 345–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Watanabe, M.; Suliman, M.E.; Qureshi, A.R.; Garcia-Lopez, E.; Barany, P.; Heimbürger, O.; Stenvinkel, P.; Lindholm, B. Consequences of low plasma histidine in chronic kidney disease patients: Associations with inflammation, oxidative stress, and mortality. Am. J. Clin. Nutr. 2008, 87, 1860–1866. [Google Scholar] [CrossRef]
  42. Arner, P.; Henjes, F.; Schwenk, J.M.; Darmanis, S.; Dahlman, I.; Iresjö, B.-M.; Naredi, P.; Agustsson, T.; Lundholm, K.; Nilsson, P.; et al. Circulating Carnosine Dipeptidase 1 Associates with Weight Loss and Poor Prognosis in Gastrointestinal Cancer. PLoS ONE 2015, 10, e0123566. [Google Scholar] [CrossRef]
  43. Carlson, A.K.; Rawle, R.A.; Wallace, C.W.; Adams, E.; Greenwood, M.C.; Bothner, B.; June, R.K. Global metabolomic profiling of human synovial fluid for rheumatoid arthritis biomarkers. Clin. Exp. Rheumatol. 2019, 37, 393–399. [Google Scholar]
  44. Li, J.; Che, N.; Xu, L.; Zhang, Q.; Wang, Q.; Tan, W.; Zhang, M. LC-MS-based serum metabolomics reveals a distinctive signature in patients with rheumatoid arthritis. Clin. Rheumatol. 2018, 37, 1493–1502. [Google Scholar] [CrossRef]
  45. Kelley, J.M.; Hughes, L.B.; Bridges, S.L., Jr. Does gamma-aminobutyric acid (GABA) influence the development of chronic inflammation in rheumatoid arthritis? J. Neuroinflamm. 2008, 5, 1. [Google Scholar] [CrossRef] [Green Version]
  46. Ham, D.J.; Caldow, M.K.; Lynch, G.; Koopman, R. Arginine protects muscle cells from wasting in vitro in an mTORC1-dependent and NO-independent manner. Amino Acids 2014, 46, 2643–2652. [Google Scholar] [CrossRef]
  47. Filippin, L.I.; Moreira, A.J.; Marroni, N.P.; Xavier, R.M. Nitric oxide and repair of skeletal muscle injury. Nitric Oxide Biol. Chem. 2009, 21, 157–163. [Google Scholar] [CrossRef]
  48. Filippin, L.I.; Cuevas, M.J.; Lima, E.; Marroni, N.P.; Gonzalez-Gallego, J.; Xavier, R.M. Nitric oxide regulates the repair of injured skeletal muscle. Nitric Oxide 2011, 24, 43–49. [Google Scholar] [CrossRef] [Green Version]
  49. Filippin, L.I.; Cuevas, M.J.; Lima, E.; Marroni, N.P.; González-Gallego, J.; Xavier, R.M. The role of nitric oxide during healing of trauma to the skeletal muscle. Inflamm. Res. 2010, 60, 347–356. [Google Scholar] [CrossRef]
  50. Merry, T.L.; Steinberg, G.R.; Lynch, G.S.; McConell, G.K. Skeletal muscle glucose uptake during contraction is regulated by nitric oxide and ROS independently of AMPK. Am. J. Physiol. Endocrinol. Metab. 2010, 298, E577–E585. [Google Scholar] [CrossRef]
  51. Urbaniak, B.; Plewa, S.; Klupczynska, A.; Sikorska, D.; Samborski, W.; Kokot, Z.J. Serum free amino acid levels in rheumatoid arthritis according to therapy and physical disability. Cytokine 2018, 113, 332–339. [Google Scholar] [CrossRef]
  52. Patel, S.S.; Molnar, M.Z.; Tayek, J.A.; Ix, J.H.; Noori, N.; Benner, D.; Heymsfield, S.; Kopple, J.D.; Kovesdy, C.P.; Kalantar-Zadeh, K. Serum creatinine as a marker of muscle mass in chronic kidney disease: Results of a cross-sectional study and review of literature. J. Cachexia Sarcopenia Muscle 2012, 4, 19–29. [Google Scholar] [CrossRef] [PubMed]
  53. Wilkinson, T.J.; Lemmey, A.B.; Jones, J.G.; Sheikh, F.; Ahmad, Y.A.; Chitale, S.; Maddison, P.J.; O’Brien, T.D. Can Creatine Supplementation Improve Body Composition and Objective Physical Function in Rheumatoid Arthritis Patients? A Randomized Controlled Trial. Arthritis Care Res. 2016, 68, 729–737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Penet, M.-F.; Gadiya, M.M.; Krishnamachary, B.; Nimmagadda, S.; Pomper, M.G.; Artemov, D.; Bhujwalla, Z.M. Metabolic Signatures Imaged in Cancer-Induced Cachexia. Cancer Res. 2011, 71, 6948–6956. [Google Scholar] [CrossRef] [Green Version]
  55. Ma, H.; Liu, X.; Wu, Y.; Zhang, N. The Intervention Effects of Acupuncture on Fatigue Induced by Exhaustive Physical Exercises: A Metabolomics Investigation. Evid. Based Complement. Altern. Med. 2015, 2015, 508302. [Google Scholar] [CrossRef]
  56. Biolo, G.; Cederholm, T.; Muscaritoli, M. Muscle contractile and metabolic dysfunction is a common feature of sarcopenia of aging and chronic diseases: From sarcopenic obesity to cachexia. Clin. Nutr. 2014, 33, 737–748. [Google Scholar] [CrossRef]
  57. Rennie, M.J.; Ahmed, A.; Khogali, S.E.; Low, S.Y.; Hundal, H.S.; Taylor, P.M. Glutamine metabolism and transport in skeletal muscle and heart and their clinical relevance. J. Nutr. 1996, 126 (Suppl. 4), 1142S–1149S. [Google Scholar] [CrossRef] [PubMed]
  58. Smith, H.J.; Greenberg, N.A.; Tisdale, M.J. Effect of eicosapentaenoic acid, protein and amino acids on protein synthesis and degradation in skeletal muscle of cachectic mice. Br. J. Cancer 2004, 91, 408–412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Chevalier, S.; Winter, A. Do patients with advanced cancer have any potential for protein anabolism in response to amino acid therapy? Curr. Opin. Clin. Nutr. Metab. Care 2014, 17, 213–218. [Google Scholar] [CrossRef]
  60. Van Norren, K.; Peters, S.J.; Van Helvoort, A.; Kegler, D.; Argiles, J.M.; Luiking, Y.C.; Laviano, A.; Van Bergenhenegouwen, J.; Deutz, N.; Haagsman, H.P.; et al. Dose-dependent effects of leucine supplementation on preservation of muscle mass in cancer cachectic mice. Oncol. Rep. 2011, 26, 247–254. [Google Scholar] [CrossRef]
  61. Liu, M.; Alimov, A.; Wang, H.; Frank, J.; Katz, W.; Xu, M.; Ke, Z.-J.; Luo, J. Thiamine deficiency induces anorexia by inhibiting hypothalamic AMPK. Neuroscience 2014, 267, 102–113. [Google Scholar] [CrossRef] [Green Version]
Figure 1. PCA (a) and PLSDA (b) analysis of CIA and CO groups with all-time points. (c) and (d) represent the plotted data in mean ± standard deviation of the mean from PCA and PLSDA analysis, respectively. 0: all animals before the beginning of the experiment (blue triangle); CIA18: CIA group 18 days after disease induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the experiment beginning of the (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO 55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). PC: principal component; LV: latent variable.
Figure 1. PCA (a) and PLSDA (b) analysis of CIA and CO groups with all-time points. (c) and (d) represent the plotted data in mean ± standard deviation of the mean from PCA and PLSDA analysis, respectively. 0: all animals before the beginning of the experiment (blue triangle); CIA18: CIA group 18 days after disease induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the experiment beginning of the (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO 55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). PC: principal component; LV: latent variable.
Jpm 11 00837 g001
Figure 2. PCA analysis at each time point between CIA and CO (pair analysis). Graphics from (af) show comparison among 18, 25, 35, 45, 55, and 65 days, respectively. CIA18: CIA group 18 days after induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the beginning of the experiment (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). PC: principal component.
Figure 2. PCA analysis at each time point between CIA and CO (pair analysis). Graphics from (af) show comparison among 18, 25, 35, 45, 55, and 65 days, respectively. CIA18: CIA group 18 days after induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the beginning of the experiment (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). PC: principal component.
Jpm 11 00837 g002
Figure 3. PLSDA analysis at each time point between CIA and CO (pair analysis). Graphics from (af) show comparison among 18, 25, 35, 45, 55, and 65 days, respectively. CIA18: CIA group 18 days after disease induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the beginning of the experiment (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). LV: latent variable.
Figure 3. PLSDA analysis at each time point between CIA and CO (pair analysis). Graphics from (af) show comparison among 18, 25, 35, 45, 55, and 65 days, respectively. CIA18: CIA group 18 days after disease induction (pink square); CIA25: CIA group 25 days after disease induction (light pink triangle); CIA35: CIA group 35 days after disease induction (dark blue inverted triangle); CIA45: CIA group 45 days after disease induction (red star); CIA55: CIA group 55 days after disease induction (light blue circle); CIA65: CIA group 65 days after disease induction (light blue diamond); CO18: CO group 18 days after the beginning of the experiment (light pink square); CO25: CO group 25 days after the beginning of the experiment (green triangle); CO35: CO group 35 days after the beginning of the experiment (orange inverted triangle); CO45: CO group 45 days after the beginning of the experiment (blue star); CO55: CO group 55 days after the beginning of the experiment (purple circle); CO65: CO group 65 days after the beginning of the experiment (yellow diamond). LV: latent variable.
Jpm 11 00837 g003
Figure 4. PCA pair analysis of CIA groups from one-time point to the next one. Graphics from (af) show comparison among 0 vs. CIA 18, CIA 18 vs. CIA 25, CIA 25 vs. CIA 35, CIA 35 vs. CIA 45, CIA 45 vs. CIA 55, and CIA 55 vs. CIA 65 days, respectively. 0: all animals before the experiment beginning (blue triangle); CIA18: CIA group 18 days after disease induction (light blue inverted triangle); CIA25: CIA group 25 days after disease induction (red star); CIA35: CIA group 35 days after disease induction (green circle); CIA45: CIA group 45 days after disease induction (purple diamond); CIA55: CIA group 55 days after disease induction (pink square); CIA65: CIA group 65 days after disease induction (light pink triangle). PC: principal component.
Figure 4. PCA pair analysis of CIA groups from one-time point to the next one. Graphics from (af) show comparison among 0 vs. CIA 18, CIA 18 vs. CIA 25, CIA 25 vs. CIA 35, CIA 35 vs. CIA 45, CIA 45 vs. CIA 55, and CIA 55 vs. CIA 65 days, respectively. 0: all animals before the experiment beginning (blue triangle); CIA18: CIA group 18 days after disease induction (light blue inverted triangle); CIA25: CIA group 25 days after disease induction (red star); CIA35: CIA group 35 days after disease induction (green circle); CIA45: CIA group 45 days after disease induction (purple diamond); CIA55: CIA group 55 days after disease induction (pink square); CIA65: CIA group 65 days after disease induction (light pink triangle). PC: principal component.
Jpm 11 00837 g004
Figure 5. PLSDA pair analysis of CIA groups from one-time point to the next one. Graphics from (af) show comparison among 0 vs. CIA 18, CIA 18 vs. CIA 25, CIA 25 vs. CIA 35, CIA 35 vs. CIA 45, CIA 45 vs. CIA 55, and CIA 55 vs. CIA 65 days, respectively. 0: all animals before the experiment beginning (blue triangle); CIA18: CIA group 18 days after disease induction (light blue inverted triangle); CIA25: CIA group 25 days after disease induction (red star); CIA35: CIA group 35 days after disease induction (green circle); CIA45: CIA group 45 days after disease induction (purple diamond); CIA55: CIA group 55 days after disease induction (pink square); CIA65: CIA group 65 days after disease induction (light pink triangle). LV: latent variable.
Figure 5. PLSDA pair analysis of CIA groups from one-time point to the next one. Graphics from (af) show comparison among 0 vs. CIA 18, CIA 18 vs. CIA 25, CIA 25 vs. CIA 35, CIA 35 vs. CIA 45, CIA 45 vs. CIA 55, and CIA 55 vs. CIA 65 days, respectively. 0: all animals before the experiment beginning (blue triangle); CIA18: CIA group 18 days after disease induction (light blue inverted triangle); CIA25: CIA group 25 days after disease induction (red star); CIA35: CIA group 35 days after disease induction (green circle); CIA45: CIA group 45 days after disease induction (purple diamond); CIA55: CIA group 55 days after disease induction (pink square); CIA65: CIA group 65 days after disease induction (light pink triangle). LV: latent variable.
Jpm 11 00837 g005
Figure 6. Venn diagram showing source of identified metabolites from BML (red) and ChenomX (blue).
Figure 6. Venn diagram showing source of identified metabolites from BML (red) and ChenomX (blue).
Jpm 11 00837 g006
Figure 7. (A) The timeline according to the analysis results of different time points. (B) Heatmap of the selected metabolite from BML related to muscle loss. (C) Network of the selected metabolite (blue diamonds) with intensity of connections (red: higher intensity; blue: lower intensity) from BML. Pathway analysis.
Figure 7. (A) The timeline according to the analysis results of different time points. (B) Heatmap of the selected metabolite from BML related to muscle loss. (C) Network of the selected metabolite (blue diamonds) with intensity of connections (red: higher intensity; blue: lower intensity) from BML. Pathway analysis.
Jpm 11 00837 g007aJpm 11 00837 g007b
Figure 8. Pathways that are related to amino acid metabolism.
Figure 8. Pathways that are related to amino acid metabolism.
Jpm 11 00837 g008
Table 1. Specificity and sensibility of NMR data after PLSDA model built obtained by the venetian blind test.
Table 1. Specificity and sensibility of NMR data after PLSDA model built obtained by the venetian blind test.
GroupCalibrationCross Validated
CO vs. CIA Same Time PointSensitivitySpecificitySensitivitySpecificity
CO 18 and CIA 181.001.000.920.67
CO 25 and CIA 250.900.830.900.67
CO 35 and CIA 351.000.900.750.70
CO 45 and CIA 450.820.780.730.56
CO 55 and CIA 551.001.000.550.33
CO 65 and CIA 65110.4440.636
Pair analysis
0 and CIA 18110.7920.727
CIA 18 and CIA 250.7270.90.6360.8
CIA 25 and CIA 35110.90.917
CIA 35 and CIA 45110.6670.727
CIA 45 and CIA 550.7270.7270.2730.091
CIA 55 and CIA 65110.5450.25
ALL GROUPS
T 00.880.800.790.77
CIA 180.830.660.750.66
CIA 250.700.570.600.62
CIA 350.830.620.830.64
CIA 451.000.560.820.55
CIA 550.9090.4390.6360.485
CIA 650.6670.4030.6670.425
CO 180.4440.7840.3330.776
CO 250.8330.7880.6670.781
CO 350.90.7520.80.789
CO 450.6670.6640.5560.657
CO 550.7780.3810.6670.373
CO 650.9090.5230.6360.53
CIA ONLY
CIA 180.920.930.670.89
CIA 250.600.730.500.75
CIA 350.920.760.750.70
CIA 450.910.520.730.56
CIA 550.820.540.640.52
CIA 650.6670.550.560.54
Table 2. Metabolites associated with muscle catabolic and anabolic pathways.
Table 2. Metabolites associated with muscle catabolic and anabolic pathways.
MetaboliteStatistics Comparison Origin (NMR or BML) aVariable Importance in Projection Score
3-MethylhistidineCIA pair analysis 45 × 55 (NMR)1.0 × 101
CIA vs. CO at 18, 25, 45, 55, and 65 (BML)6.6 × 10°, 4.8 × 10°, 2.7 × 10°, 1.9 × 10°, 2.0 × 10°
CIA pair analysis 25 × 35, 35 × 45 (BML)2.0 × 10°, 1.3 × 10°
4-AminobutyrateCIA vs. CO at 18, 25, 35, 45, and 55 days (BML)1.3 × 10°, 3.6 × 101, 1.7 × 101, 1.4 × 10°, 1.5 × 10°
CIA pair analysis 18 × 25, 25 × 35, 35 × 45 (BML)1.8 × 10°, 3.5 × 10°, 4.9 × 10°
AcetylcholineCIA vs. CO at 45, and 65 days (BML)1.6 × 10°, 1.4 × 10°
CIA pair analysis 35 × 45, and 45 × 55 (NMR)4.9 × 101, 3.8 × 101
ArginineCIA vs. CO at 18, 35, and 65 days (BML)1.9 × 10°, 2.4 × 10°, 1.8 × 10°
AspartateCIA vs. CO at 18 days (BML)2.0 × 10°
CarnosineCIA vs. CO at 65 days (NMR)2.8 × 102
CIA pair analysis 25 × 35 (NMR)3.6 × 10°
CIA vs. CO at 18, 25, 35, 45, 55, and 65 days (BML)2.1 × 10°, 1.7 × 101, 4.6 × 101, 1.2 × 101, 2.9 × 101, 5.0 × 10°
CIA pair analysis 25 × 35, 35 × 45, 45 × 55, 55 × 65 (BML)5.9 × 10°, 1.6 × 101, 3.1 × 10°, 4.4 × 101
CreatineCIA pair analysis 35 × 45, 45 × 55 (NMR)2.3 × 10°, 8.4 × 101
CIA vs. CO at 18, 35, 55, and 65 days (BML)1.8 × 10°, 1.4 × 10°, 3.3 × 10°, 1.4 × 10°
CreatinineCIA vs. CO at 45 days (NMR)3.4 × 101
CIA pair analysis 35 × 45 (NMR)3.1 × 101
CIA vs. CO at 25, 35, 45, and 65 days (BML)3.3 × 10°, 9.2 × 10°, 1.6 × 10°, 4.0 × 10°,
CIA pair analysis 25 × 353.7 × 10°
GlutamineCIA vs. CO at 18, 25, 45, 55, and 65 (BML)3.1 × 10°, 1.8 × 10°, 3.2 × 10°, 3.4 × 10°, 2.4 × 10°
CIA pair analysis 25 × 35, 45 × 55 (BML)5.7 × 10°, 1.2 × 10°
HistamineCIA vs. CO at 35, 45, 55, and 65 days (NMR)2.8 × 101, 3.7 × 101, 6.2 × 101, 5.2 × 101
CIA pair analysis 0 × 18, 18 × 25, 25 × 35, 35 × 45, 45 × 55, 55 × 65 (NMR)8.7 × 10°, 8.5 × 10°, 4.9 × 10°, 1.3 × 101, 1.9 × 101, 1.1 × 102
CIA vs. CO at 35, 45, and 65 days (BML)3.9 × 10°, 3.7 × 10°, 4.9 × 10°
CIA pair analysis 25 × 35 (BML)4.2 × 10°
HistidineCIA vs. CO at 18 days (NMR)2.6 × 101
CIA pair analysis 55 × 65 (NMR)1.2 × 102
CIA vs. CO at 35, 55, and 65 days (BML)2.5 × 10°, 1.2 × 10°, 2.5 × 10°
IsoleucineCIA vs. CO at 45 days (NMR)2.1 × 101
CIA vs. CO at 18, 25, 35, 45, and 65 days (BML)5.3 × 10°, 1.1 × 10°, 1.5 × 10°, 1.7 × 10°, 1.7 × 10°,
LeucineCIA vs. CO at 18, 25, 35, 45, 55, and 65 days (BML)1.0 × 10°, 2.7 × 10°, 3.5 × 10°, 1.1 × 10°, 1.5 × 10°, 2.2 × 10°
CIA pair analysis 25 × 35, 35 × 45 (BML)1.8 × 10°, 1.6 × 10°
L-MethionineCIA vs. CO at 18 days (NMR)5.4 × 101
CIA pair analysis 0 × 18, 18 × 25 (NMR)1.1 × 102, 2.1 × 10°
CIA vs. CO at 25, and 55 days (BML)1.2 × 10°, 5.1 × 10°
LysineCIA pair analysis 0 × 18, 25 × 35, 45 × 55, 55 × 65 (NMR)2.8 × 101, 5.3 × 10°, 3.1 × 10°, 6.1 × 10°, 2.5 × 10°
CIA vs. CO at 18, and 45 days (BML)5.4 × 10°, 3.2 × 10°
myo-InositolCIA vs. CO at 18, 35, 45, 55, and 65 days (BML)1.7 × 10°, 4.9 × 10°, 1.1 × 10°, 1.6 × 10°, 1.2 × 10°
CIA pair analysis 25 × 35 (BML)2.1 × 101
N,N-DimethylglycineCIA vs. CO at 18, 25, 55, and 65 days (BML)3.5 × 10°, 2.4 × 10°, 2.6 × 10°, 1.3 × 10°
CIA pair analysis 25 × 35 (BML)4.1 × 10°
N-AcetylalanineCIA vs. CO at 18, 25, 45, and 55 days (BML)1.5 × 10°, 2.8 × 10°, 1.7 × 10°, 5.4 × 10°
N-AcetylmethionineCIA vs. CO at 18, 25, and 65 days (BML)3.1 × 10°, 8.7 × 10°, 3.4 × 10°
PantothenateCIA vs. CO at 35 days (NMR)2.1 × 101
CIA pair analysis 35 × 45, 45 × 55, 55 × 65 (NMR)1.4 × 101, 3.3 × 101, 8.8 × 10°
CIA vs. CO at 18, 25, 35, and 55 days (BML)2.7 × 10°, 1.8 × 10°, 1.0 × 10°, 1.8 × 10°
PhenylalanineCIA vs. CO at 18 days (BML)1.8 × 10°
PhosphocholineCIA vs. CO at 25, 35, and 55 days (BML)4.3 × 10°, 1.6 × 10°, 1.7 × 10°
PhosphocreatineCIA vs. CO at 45 days (NMR)3.4 × 101
CIA pair analysis 35 × 45, 45 × 55 (NMR)3.1 × 101, 8.4 × 101
CIA vs. CO 25, 35, 55, and 65 days (BML)1.9 × 10°, 2.9 × 10°, 6.2 × 10°, 3.3 × 10°
PyridoxineCIA vs. CO at 35, and 55 days (NMR)4.6 × 101, 5.2 × 101
SarcosineCIA vs. CO at 25, 35, 45, and 65 days (BML)1.0 × 10°, 3.3 × 10°, 2.3 × 10°, 1.9 × 10°
SuccinylacetoneCIA vs. CO at 18, 25, 55, and 65 days (NMR)7.0 × 101, 4.7 × 101, 6.8 × 101, 3.7 × 101, 1.0 × 102
CIA pair analysis 55 × 65 (NMR)5.2 × 101
ThiamineCIA pair analysis 25 × 35 (BML)2.8 × 10°
UrocanateCIA pair analysis 55 × 65 (NMR)5.0 × 10°
a, NMR: 1D nuclear magnetic resonance; BML: 2D J-res nuclear magnetic resonance identified by the Birmingham Metabolome Library.
Table 3. Metabolic pathway analysis of metabolite ranking data from PLSDA models comparing CIA vs. CO at each time-point and CIA pair analysis at each time. Data expressed as “% of metabolic pathway impact (p-value)”.
Table 3. Metabolic pathway analysis of metabolite ranking data from PLSDA models comparing CIA vs. CO at each time-point and CIA pair analysis at each time. Data expressed as “% of metabolic pathway impact (p-value)”.
BMLAlanine, Aspartate, and Glutamate MetabolismArginine and Proline MetabolismAscorbate and Aldarate MetabolismB Vitamin ComplexButanoate MetabolismGlycerophospholipid MetabolismGlycine, Serine, and Threonine MetabolismHistidine MetabolismLysine MetabolismNitrogen MetabolismPantothenate and CoA BiosynthesisValine, Leucine, and Isoleucine Degradation
CIA pair analysis 0 × 18 22 (0.002)0 (0.066)
CIA pair analysis 18 × 25 14 (0.040) 3 (0.326) 22 (0.027)
CIA pair analysis 25 × 3526 (0.070)28 (p < 0.001)0 (0.154) 3 (0.337) 3 (0.109)22 (0.029) 0 (0.154) 33 (0.185)
CIA pair analysis 35 × 45 22 (0.007) 3 (0.337)11 (0.103)2 (0.109)0 (0.244) 33 (0.185
CIA pair analysis 45 × 5515 (0.074)21 (0.001) 11 (0.110)0 (0.453)22 (0.032)0 (0.074)0 (0.159)
CIA pair analysis 55 ×65 0 (0.002)
NMR
CIA pair analysis 0 × 18 22 (0.166)0 (0.030)
CIA pair analysis 18 × 25 22 (0.130)
CIA pair analysis 25 × 35 0 (0.284)22 (0.010)0 (0.023) 2 (0.336)
CIA pair analysis 35 × 45 1 (0.487) 0 (0.364)0 (0.374)22 (0.202) 2 (0.202)2 (0.437)
CIA pair analysis 45 × 55 1 (0.503) 0 (0.378) 22 (0.210)0 (0.001) 2 (0.021)
CIA pair analysis 55 × 65 61 (0.001)0 (0.001) 2 (0.184)
CIA vs. CO at 18 days6 (0.113) 0 (0.104) 0 (0.144)24 (0.072)0 (0.020)0 (0.044)
CIA vs. CO at 25 days0 (0.113) 0 (0.104) 0 (0.009) 0 (0.053)
CIA vs. CO at 35 days 0 (0.010) 0 (0.199)22 (0.101) 2 (0.101)0 (0.075)
CIA vs. CO at 45 days0 (0.227)0 (0.378) 0 (0.021) 0 (0.284)22 (0.148)0 (0.042) 0 (0.111)
CIA vs. CO at 55 days0 (0.157) 8 (0.062)0 (0.145)4 (0.193)0 (0.199)22 (0.101) 0 (0.075)
CIA vs. CO at 65 days 0 (0.144)22 (0.002) 2 (0.174)
Table 4. Regression analysis between clinical data and metabolic profile found in urine and identified by BML. Metabolites shown are solely the ones related to muscle loss process. Value expresses the metabolite contribution to the model.
Table 4. Regression analysis between clinical data and metabolic profile found in urine and identified by BML. Metabolites shown are solely the ones related to muscle loss process. Value expresses the metabolite contribution to the model.
StrengthFatigueLocomotionScoreEdemaIntakeWeight
3-Methylhistidine0.044−0.033−0.020−0.020−0.0510.130−0.048
4-Aminobutyate−0.073−0.071−0.0160.0560.0970.0220.001
Acetylcholine0.0630.0150.015−0.024−0.0360.068−0.018
Arginine0.083−0.0040.031−0.093−0.0350.0040.029
Carnosine−0.0390.0390.0310.0930.125−0.051−0.026
Creatine−0.063−0.014−0.0510.0500.026−0.177−0.039
Creatinine−0.056−0.061−0.0620.1530.1070.1420.058
Glutamine−0.0190.0240.0030.0580.0660.0660.024
Histamine−0.029−0.035−0.0070.0010.036−0.0040.019
Histidine−0.048−0.0400.0220.0330.0590.036−0.050
Isoleucine−0.046−0.028−0.0180.0520.0490.000−0.030
L-Methionine−0.0820.003−0.0060.1090.048−0.0200.024
Leucine−0.101−0.069−0.0280.0830.1210.041−0.022
Lysine0.020−0.006−0.020−0.033−0.070−0.180−0.080
myo-Inositol−0.011−0.026−0.0440.0390.0080.050−0.012
N-Acetylalanine0.0280.0930.0630.031−0.0080.2140.032
N-Acetylmethionine0.0140.0560.021−0.060−0.0010.0350.014
N,N-Dimethylglycine0.107−0.022−0.014−0.044−0.070−0.006−0.020
Phosphocholine0.0160.0410.0050.022−0.0250.0220.000
Phosphocreatine0.0330.023−0.0100.0370.0600.000−0.007
Sarcosine0.0370.0660.053−0.017−0.0230.0830.028
Thiamine0.0840.0120.126−0.047−0.1560.0250.000
Values express the metabolite contribution to the model.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alabarse, P.V.G.; Silva, J.M.S.; Santo, R.C.E.; Oliveira, M.S.; Almeida, A.S.; de Oliveira, M.S.; Immig, M.L.; Freitas, E.C.; Teixeira, V.O.N.; Bathurst, C.L.; et al. Metabolomic Biomarker Candidates for Skeletal Muscle Loss in the Collagen-Induced Arthritis (CIA) Model. J. Pers. Med. 2021, 11, 837. https://doi.org/10.3390/jpm11090837

AMA Style

Alabarse PVG, Silva JMS, Santo RCE, Oliveira MS, Almeida AS, de Oliveira MS, Immig ML, Freitas EC, Teixeira VON, Bathurst CL, et al. Metabolomic Biomarker Candidates for Skeletal Muscle Loss in the Collagen-Induced Arthritis (CIA) Model. Journal of Personalized Medicine. 2021; 11(9):837. https://doi.org/10.3390/jpm11090837

Chicago/Turabian Style

Alabarse, Paulo V. G., Jordana M. S. Silva, Rafaela C. E. Santo, Marianne S. Oliveira, Andrelise S. Almeida, Mayara S. de Oliveira, Mônica L. Immig, Eduarda C. Freitas, Vivian O. N. Teixeira, Camilla L. Bathurst, and et al. 2021. "Metabolomic Biomarker Candidates for Skeletal Muscle Loss in the Collagen-Induced Arthritis (CIA) Model" Journal of Personalized Medicine 11, no. 9: 837. https://doi.org/10.3390/jpm11090837

APA Style

Alabarse, P. V. G., Silva, J. M. S., Santo, R. C. E., Oliveira, M. S., Almeida, A. S., de Oliveira, M. S., Immig, M. L., Freitas, E. C., Teixeira, V. O. N., Bathurst, C. L., Brenol, C. V., Filippin, L. I., Young, S. P., Lora, P. S., & Xavier, R. M. (2021). Metabolomic Biomarker Candidates for Skeletal Muscle Loss in the Collagen-Induced Arthritis (CIA) Model. Journal of Personalized Medicine, 11(9), 837. https://doi.org/10.3390/jpm11090837

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop