Skip to main content

Altered skeletal muscle metabolic pathways, age, systemic inflammation, and low cardiorespiratory fitness associate with improvements in disease activity following high-intensity interval training in persons with rheumatoid arthritis



Exercise training, including high-intensity interval training (HIIT), improves rheumatoid arthritis (RA) inflammatory disease activity via unclear mechanisms. Because exercise requires skeletal muscle, skeletal muscle molecular pathways may contribute. The purpose of this study was to identify connections between skeletal muscle molecular pathways, RA disease activity, and RA disease activity improvements following HIIT.


RA disease activity assessments and vastus lateralis skeletal muscle biopsies were performed in two separate cohorts of persons with established, seropositive, and/or erosive RA. Body composition and objective physical activity assessments were also performed in both the cross-sectional cohort and the longitudinal group before and after 10 weeks of HIIT. Baseline clinical assessments and muscle RNA gene expression were correlated with RA disease activity score in 28 joints (DAS-28) and DAS-28 improvements following HIIT. Skeletal muscle gene expression changes with HIIT were evaluated using analysis of covariance and biological pathway analysis.


RA inflammatory disease activity was associated with greater amounts of intramuscular adiposity and less vigorous aerobic exercise (both p < 0.05). HIIT-induced disease activity improvements were greatest in those with an older age, elevated erythrocyte sedimentation rate, low cardiorespiratory fitness, and a skeletal muscle molecular profile indicative of altered metabolic pathways (p < 0.05 for all). Specifically, disease activity improvements were linked to baseline expression of RA skeletal muscle genes with cellular functions to (1) increase amino acid catabolism and interconversion (GLDC, BCKDHB, AASS, PYCR, RPL15), (2) increase glycolytic lactate production (AGL, PDK2, LDHB, HIF1A), and (3) reduce oxidative metabolism via altered beta-oxidation (PXMP2, ACSS2), TCA cycle flux (OGDH, SUCLA2, MDH1B), and electron transport chain complex I function (NDUFV3). The muscle mitochondrial glycine cleavage system (GCS) was identified as critically involved in RA disease activity improvements given upregulation of multiple GCS genes at baseline, while GLDC was significantly downregulated following HIIT.


In the absence of physical activity, RA inflammatory disease activity is associated with transcriptional remodeling of skeletal muscle metabolism. Following exercise training, the greatest improvements in disease activity occur in older, more inflamed, and less fit persons with RA. These exercise training-induced immunomodulatory changes may occur via reprogramming muscle bioenergetic and amino acid/protein homeostatic pathways.

Trial registration, NCT02528344. Registered on 19 August 2015.


Exercise training improves disease activity and systemic inflammation in rheumatoid arthritis (RA) [1,2,3,4]. In addition to traditional aerobic exercise programs, high-intensity interval training (HIIT)—where participants alternate bouts of near-maximal intensity aerobic exercise with bouts of lower intensity—significantly improves RA Disease Activity Score in 28 joints (DAS-28) and erythrocyte sedimentation rate (ESR) [5]. HIIT also significantly improves RA peripheral neutrophil chemotaxis and macrophage phagocytosis; findings suggest that exercise improves RA innate immune function and pathogen response [5]. In non-RA populations, chronic exercise training and lifelong physical activity improve regulatory T cell anti-inflammatory function, decrease CD4+ Th17 polarization, and ameliorate immunosenescence [6, 7]; however, the specific effects of exercise on improving RA chronic inflammation and adaptive immunity are unknown. Improved understanding of the mechanisms that connect exercise training with its immunomodulatory effects holds great promise to enhance future lifestyle and pharmacologic therapies for RA and other chronic inflammatory diseases.

One underexplored mediator of exercise-induced improvements in RA inflammation is skeletal muscle. Skeletal muscle is the predominant locomotive organ and is largely responsible for the beneficial health effects of exercise. Skeletal muscle, during acute exercise, communicates widely with other organ systems, including immune cells, via the secretion of proteins termed “myokines” [8]. Myokines have far-ranging benefits, including helping to regulate systemic energy metabolism [8,9,10,11,12]. For example, during exercise, interleukin (IL)-6 and IL-15 are secreted by muscle in a pulsatile fashion and, in contrast to sepsis or inflammation, independent from tumor necrosis factor (TNF)-α secretion. Muscle IL-6 secretion during exercise has direct anti-inflammatory effects by stimulating IL-10 production from CD4+ T cells and inhibiting systemic TNF-α secretion [13, 14]. While skeletal muscle contractile activity appears to improve immune function, chronic immune activation (i.e., RA) adversely impacts skeletal muscle.

In RA, chronic inflammation is strongly linked to impaired muscle function and “rheumatoid cachexia” or sarcopenic obesity (i.e., decreased muscle mass, increased fat mass, and intramuscular adiposity) [15]. Further, the RA skeletal muscle molecular profile is defined by altered inflammatory cytokines and impairments in satellite cell function and oxidative metabolism [15,16,17,18]. Of particular interest to this study, it is unclear how these RA muscle molecular pathways specifically (1) contribute to and are affected by inflammation and autoimmunity and (2) influence and predict immune system changes with exercise training. In a cross-sectional cohort of persons with established RA, we identified clinical and skeletal muscle molecular phenotypes associated with RA inflammatory disease activity. To better understand how exercise-trained skeletal muscle modulates inflammation, we then studied an independent RA cohort that completed a 10-week HIIT exercise program and identified baseline clinical and skeletal muscle molecular phenotypes associated with improvements in RA disease activity. We hypothesized that altered skeletal muscle immune and metabolic pathways are intricately linked to (1) RA inflammatory disease activity and (2) improvements in RA inflammation following exercise training.


Design and participants

For the cross-sectional RA cohort #1, persons with RA (n = 47) completed clinical assessments and underwent skeletal muscle biopsies in a cross-sectional design [16]. All RA participants met American College of Rheumatology (ACR) 1987 criteria [19]. Participants had a positive rheumatoid factor and/or anti-cyclic citrullinated peptide antibody or had imaging evidence of hand or foot erosions. They had no medication changes within 3 months of enrollment and were all using less than or equal to 5 mg of prednisone daily. RA participants were excluded with pregnancy, type 2 diabetes mellitus, or known cardiovascular disease.

For HIIT RA cohort #2, persons with RA (n = 12) completed clinical assessments and underwent skeletal muscle biopsies before and after a 10-week treadmill walking HIIT program [5, 17]. Similar to the cross-sectional RA cohort #1, all RA participants satisfied the ACR 1987 RA classification criteria [19], had no medication changes for 3 months prior to or during the intervention, and were taking less than or equal to prednisone 5 mg per day. RA medications used during the study were previously reported [17]. For the HIIT RA cohort #2, inclusion required participants to be sedentary, defined as exercising less than 2 days per week. Exclusions included diabetes mellitus, cardiovascular disease, and if they were unable to walk without assistance on a treadmill.

Exercise intervention

In the HIIT RA cohort #2, persons with RA completed supervised treadmill walking exercise sessions three times per week for 10 weeks as previously described [5, 17]. Exercise prescriptions were based on baseline cardiopulmonary exercise tests. During each treadmill walking exercise session, participants completed (1) a 5-min warm-up, (2) 10 alternating intervals (60–90 s each) of high-intensity (80–90% heart rate reserve) and low-intensity (50–60% heart rate reserve) exercise, and (3) a 5-min cool-down.

Clinical assessments

For RA cohorts #1 and #2, health history questionnaires were obtained and inflammatory disease activity was assessed by DAS-28 determined from patient-completed pain and overall health visual analog scales, physician assessed tender and swollen joint counts, and ESR [20]. For the cross-sectional RA cohort #1, physical activity was assessed via a wearable accelerometer for 7 days [16, 21]; body composition was assessed via abdominal and thigh CT scans [22]. For the HIIT RA cohort #2, cardiorespiratory fitness was assessed via graded treadmill exercise testing as maximal oxygen consumption during exercise (VO2 peak) [5, 16]; body composition was assessed via BodPod® (BodPod System; Life Measurement Corporation, Concord, CA, USA) [17].

Skeletal muscle gene expression profiling

Vastus lateralis needle biopsies were obtained from participants after an overnight fast by standard Bergstrom technique and stored at −80°C in the cross-sectional RA cohort #1 and before and after 10 weeks of HIIT in RA cohort #2 at between 24 and 48 h after the last exercise training bout [17]. Illumina Human HT-12v3 and HT-12v4 Expression BeadChips were used for muscle quantitative RNA analyses for a random subset of cross-sectional RA cohort #1 (n = 20) and the HIIT RA cohort #2 (n = 12), respectively. For each frozen sample, 20–30 mg of the muscle was homogenized; biotinylated total RNA was prepared using the Illumina TotalPrep RNA amplification kit (Life Technologies, Grand Island, NY, USA) and hybridized to the Human HT-12 BeadChips (Illumina, San Diego, CA, USA); and RNA profiling performed as previously described [16].

Statistical analyses

Differential Spearman correlations were performed using Partek Genomics Suite (Partek, Inc.; St. Louis, MO, USA) to assess associations of baseline muscle gene expression with (a) baseline DAS-28 in the cross-sectional RA cohort #1 and (b) HIIT-mediated change (i.e., improvement) in RA disease activity (ΔDAS-28) in the HIIT RA cohort #2. For the HIIT RA cohort #2, ΔDAS-28 was computed as pre-HIIT DAS-28 minus post-HIIT DAS-28 so that larger numbers and positive associations represent improved disease activity. To validate baseline DAS-28 correlations from the cross-sectional RA cohort #1, Spearman correlations were also performed between baseline muscle gene expression and baseline DAS-28 in the HIIT RA cohort #2. Genes with expression correlated with DAS-28 or ΔDAS28 (Spearman rho p < 0.05) were uploaded for Ingenuity Pathway Analyses (IPA,; positive and negative coefficients represented positive and inverse associations, respectively, between gene expression and DAS-28 (cross-sectional RA cohort #1) or ΔDAS-28 (HIIT RA cohort #2). IPA-defined canonical pathways—created based on right-handed Tukey’s T tests to analyze overrepresentation of genes and pathways based on molecular pathways in the current literature—were identified as pathways with the greatest proportional representation (gene set/total pathway genes) (Fisher Exact Test; p < 0.05). Of importance to the interpretation of our analyses, we purposely chose a less stringent significance cutoff (i.e., p < 0.05 as opposed to false discovery rate) when creating initial gene sets to minimize random false positive and highly related results with downstream pathway analyses. Top differentially associated gene lists (p < 0.001) from both RA cohorts were analyzed for functional classification using the PANTHER Gene List Analysis tool and for overrepresentation against a homo sapiens gene reference database using the PANTHER Overrepresentation Test ( with Reactome pathways ( annotation [23,24,25]. Changes in gene expression fold differences after HIIT in RA cohort #2 were assessed via analysis of covariance (ANCOVA) modeling in Partek Genomics Suite.


Clinical associations with RA disease activity

In the cross-sectional RA cohort #1, disease activity correlated with ESR (rho=0.80, p < 0.001), low thigh muscle density (i.e., greater intramuscular adipose tissue) [26] (rho=−0.32, p = 0.032), and fewer total daily minutes of high (rho=−0.36, p = 0.025) and very high intensity physical activity (rho=−0.37, p = 0.02) (Table 1).

Table 1 Rheumatoid arthritis clinical characteristics and disease activity relationships

Skeletal muscle gene expression associations with RA disease activity

Of 29,042 total mapped gene transcripts analyzed, expression of 25 genes was highly correlated (rho>0.68 or <−0.68; p < 0.001) with RA disease activity scores (DAS-28) in the cross-sectional RA cohort #1. The top highly correlated skeletal muscle gene was IL1RL2 (interleukin 1 receptor-like 2; rho=−0.84) (Supplementary Table 1). IL1RL2 expression in RA skeletal muscle was also significantly correlated with DAS-28 at baseline in the HIIT RA cohort #2 (Supplementary Table 1). The 25 highly correlated genes were categorized according to biological function (PANTHER gene ontology (GO)) (Fig. 1A) and pathways (Reactome) (Fig. 1B). Of the 25 highly associated genes, we identified Reactome pathways that were overrepresented. Without correction, 5 non-overlapping Reactome pathways were overrepresented (>2 genes/set) (p < 0.05): (1) Basigin interactions, (2) cell junction organization, (3) cell–cell communication, (4) interleukin-1 family signaling, and (5) ion channel transport. However, accounting for false discovery, no Reactome pathway was significantly overrepresented.

Fig. 1
figure 1

Top rheumatoid arthritis skeletal muscle gene associations with inflammatory disease activity. a Bar chart categorizing top rheumatoid arthritis (RA) skeletal muscle differentially expressed genes highly associated (Spearman’s rho p < 0.001) with RA disease activity score in 28 joints (DAS-28) by PANTHER gene ontology (GO) biological processes (bars represent total number genes per category; percent gene process hits/total listed next to bars). b Bar chart showing Reactome pathways overrepresented (>2 genes per pathway) within top RA skeletal muscle genes highly associated (p < 0.001) with disease activity. c Bar chart categorizing top baseline RA skeletal muscle differentially expressed genes highly associated (Spearman’s rho p < 0.001) with improvement in RA disease activity following high-intensity interval training (HIIT) by PANTHER GO biological processes (bars represent total number genes per category; percent gene process hits/total listed next to bars). d Bar chart showing Reactome pathways overrepresented (>2 genes per pathway) within top baseline RA skeletal muscle genes highly associated (p < 0.001) with improvement in disease activity following HIIT. *Represents pathway that reached significance at the calculated false discovery rate (p < 0.05)

To further explore relationships between the skeletal muscle transcriptome and RA systemic inflammatory disease, we uploaded skeletal muscle genes differentially associated with disease activity to biological pathway software (IPA) (Spearman’s correlation p < 0.05). These data included a total of 1417 molecules (684 negatively correlated, 733 positively correlated). The top canonical pathway was “DNA Methylation and Transcriptional Repression Signaling” (Table 2 and Supplementary Table 2), represented by 7 genes (CHD4, MTA2, SIN3A, DNMT3A, DNMT3B, H4C14, MBD2) of 33 total genes in that pathway. Of the 10 top canonical pathways in the cross-sectional RA cohort #1, 5 included at least one significant skeletal muscle gene relationship that was validated in the HIIT RA cohort #2 (Supplementary Table 2).

Table 2 Cross-sectional rheumatoid arthritis cohort #1: skeletal muscle canonical pathways associated with disease activity

Clinical associations with exercise training-related improvements in RA disease activity

In the HIIT RA cohort #2, HIIT-mediated improvements in DAS-28 correlated with older age (rho=0.67, p = 0.022), higher baseline ESR (rho=0.64, p = 0.035), and lower baseline cardiorespiratory fitness/VO2 peak (rho=−0.75, p = 0.001) (Table 1).

Baseline skeletal muscle gene expression associations with exercise training-related improvements in RA disease activity

Muscle expression for 43 genes of 27,911 total mapped gene transcripts analyzed in the HIIT RA cohort #2 was highly correlated with improvements in disease activity (rho>0.80 or <−0.80; p < 0.001). The strongest correlation was for KIF25 (kinesin family member 25; rho=0.915) (Supplementary Table 3). Within the 43 genes, all 4 non-overlapping, overrepresented Reactome pathways (without correction; p < 0.05) involved cellular substrate energy metabolism: (1) glycine degradation, (2) metabolism of amino acids and derivatives, (3) pyruvate metabolism, and (4) tricarboxylic acid (TCA) cycle and respiratory electron transport. The glycine degradation pathway (R-HSA-389661) was the only pathway significantly overrepresented at the false discovery rate (Fig. 1D). Of genes where baseline expression was associated with HIIT-mediated improvement in disease activity, KIF25 (kinesin family member 25) and MNT (MAX network transcriptional repressor) were significantly upregulated and GLDC (glycine decarboxylase) was significantly downregulated following HIIT (p < 0.05) (Supplementary Table 3).

We used IPA to better understand skeletal muscle pathways that may contribute to exercise-induced immune regulation. Using a p value cut off of 0.05 for Spearman’s correlations, 1527 molecules (1049 negatively correlated, 478 positively correlated) were uploaded. IPA identified 39 canonical pathways (p < 0.05). The top canonical pathway was “Proline Biosynthesis II” (Table 3 and Supplementary Table 4), represented by 4 genes (OAT, PYCR1, PYCR2, PYCR3) of 5 total genes in that pathway. The top canonical pathway not directly related to energy substrate metabolism was “nNOS Signaling in Skeletal Muscle Cells” (Table 3), which included 9 (CACNA1A, CACNA1C, CACNA1I, CACNA2D4, CACNB4, CACNG8, CHRNA1, RYR1, CACNB1) of 40 pathway genes. The only canonical pathway represented in both analyses for RA cohorts #1 and #2 was “Branched-chain α-keto acid Dehydrogenase Complex” (Tables 2 and 3).

Table 3 High-intensity interval training rheumatoid arthritis cohort #2: skeletal muscle canonical pathways associated with improvements in disease activity


Via analysis of two cohorts of established RA (one cross-sectional and one prior to and after exercise training), we identified a phenotype for an exercise training-induced anti-inflammatory response as being older, more inflamed, less aerobically fit, and with multiple alterations in skeletal muscle metabolic pathways. Based on our transcriptomic analyses of top highly associated genes and IPA canonical pathways, we have identified the apparent skeletal muscle cellular metabolic pathways that (1) are connected to RA inflammation and (2) are associated with exercise training modulation of systemic immune responses (Fig. 2). These skeletal muscle pathways are highlighted by major gene expression alterations in amino acid catabolism and the regulation of glycolysis and TCA cycle flux. Taken together, rewiring of protein homeostasis and oxidative metabolism in sedentary RA skeletal muscle may be connected to the perpetuation of systemic inflammation. This association of the RA muscle transcriptional profile with HIIT-mediated improvements in RA disease activity suggests that exercise modulation of inflammation occurs in concert with the reprogramming of skeletal muscle metabolism.

Fig. 2
figure 2

Is altered skeletal muscle metabolism linked to fueling chronic inflammation? Figure summarizes skeletal muscle cellular metabolic pathways that may contribute to rheumatoid arthritis (RA) chronic inflammation based on differentially expressed RA muscle genes that are highly associated with disease activity score in 28 joints (DAS-28) (*) or improvements in DAS-28 following high-intensity interval training (HIIT) (#). Genes included in the pathways diagram are either highly associated (Spearman’s rho p < 0.001) or part of significantly associated ingenuity canonical pathways (p < 0.05). Genes whose expression significantly (p < 0.05) changes following HIIT are listed ($). Based on direction of positive (↑) or negative (↓) associations, genes whose function is to promote/upregulate (blue) or inhibit/downregulate (red) a specific pathway are highlighted. Bolded arrows represent proposed upregulated pathways. Dashed arrows represent proposed downregulated pathways. Green end-dotted arrows represent hypothesized pathways in which altered skeletal muscle metabolism may contribute to chronic immune-activation and inflammation in RA

Amino acid homeostasis

In RA, amino acids are critical to sustain chronic inflammation. In RA synovium, chronic inflammation upregulates amino acid transporters, fueling maladaptive synovial proliferation and bone remodeling pathways [27]. In the circulation, activated RA immune cells have high energy demands that are met with glutamine, the majority of which is produced by skeletal muscle [28,29,30]. Further, inhibition of glutaminolysis is a promising therapeutic target for the management of Th17-driven autoimmune diseases [31, 32]. Our data demonstrate that muscle catabolism and interconversion of multiple amino acids—including glycine, lysine, glutamate/glutamine, and the branched-chain amino acids leucine, isoleucine, and valine—are linked to both RA inflammation and exercise training-induced reduction of inflammation (Fig. 2).

Chronic immune activity and exercise training-related modulation appear strongly related to the RA muscle mitochondrial glycine cleavage pathway. Baseline upregulation of multiple key glycine cleavage system genes (GLDC, GCSH, LIAS, AMT) associated with HIIT-related improvements in disease activity, while GLDC (glycine dehydrogenase) was significantly downregulated following HIIT. Glycine is a proteogenic amino acid and critical for multiple metabolic pathways in humans [33]. Catabolism of glycine via the glycine cleavage pathway to N5,N10-methylenetetrahydrofolate is utilized for the biosynthesis of purine nucleotides [34, 35], key regulators of T lymphocyte proliferation and survival [36, 37]. Thus, increased activity of the glycine cleavage system may be one way that impaired RA skeletal muscle metabolism contributes to and is impacted by chronic RA immune dysfunction.

Altered skeletal muscle metabolism of other amino acids also link to RA inflammation. Skeletal muscle branched-chain amino acid (BCAA) catabolism—via increased mitochondrial branched-chain α-ketoacid (BCKA) dehydrogenase (BCKDH) complex activity—was associated with both baseline RA disease activity and with exercise-mediated disease activity improvements. We theorize that increased RA muscle BCKDH activity might lead to glutamate/glutamine formation via multiple pathways. First, BCAA are metabolized to both BCKA and glutamate directly. Second, BCKA are converted to acetyl-CoA; subsequent α-ketoglutarate flux from the TCA cycle leads to further glutamate production via ornithine aminotransferase/OAT [25, 34]. Glutamate is then readily converted to glutamine with glutamine production augmented by lysine degradation via increased muscle AACC activity. Glutamine released into circulation can then sustain activation of RA immune cells [30, 34, 38] (Fig. 2).

Interestingly, opposing protein homeostatic pathways involved in both protein synthesis and muscle atrophy associated with subsequent improvement in RA inflammatory disease activity following HIIT. These pathways were highlighted by greater gene expression of RPL15 (a ribosomal protein critical for muscle protein synthesis) [34] and the nNOS calcium signaling canonical pathway important for driving muscle atrophy [39]. Taken together, these altered pathways suggest that RA skeletal muscle tissue is in a heightened state of protein turn-over; further, exercise training may work in part to regulate these protein catabolic and anabolic programs. Our findings of upregulated RA muscle protein turnover are in agreement with the prominent clinical phenotype of muscle loss due to “rheumatoid cachexia” or sarcopenic obesity [15]. Therefore, the increased ribosomal activity in RA muscle may represent a failed compensatory response to maintain muscle mass in the face of constant protein breakdown and atrophy signaling triggered by chronic inflammation. One unanswered question is whether this state of dysregulated RA muscle amino acid/protein homeostasis is merely a byproduct of inflammation or if muscle protein breakdown truly contributes to direct fueling of chronic immune activity [30].

Glycolysis regulation

In addition to amino acids, dysregulation of RA skeletal muscle glycolysis and lactate production appears to support chronic immune signaling. We previously showed that, as compared to healthy matched controls, RA skeletal muscle is characterized by increased glycolysis as evidenced by transcriptional downregulation of oxidative metabolism components and accumulation of pyruvate [16]. Here, in a separate RA cohort, additional pro-glycolytic programs were linked to exercise-mediated improvements in RA disease activity. These transcriptomic features include increased pre-training breakdown of glycogen to glucose, decreased conversion of pyruvate to acetyl CoA, and increased interconversion between pyruvate and lactate [34] (Fig. 2). It is unclear from our data alone if an increased reliance on muscle glycolysis and subsequent lactate production contributes to RA immune cytokine signaling. Of note, the release of lactate from tumor and inflamed tissue—including RA synovium—stimulates pro-inflammatory IL-17 secretion [40,41,42]. Further, our data show that increased RA muscle hypoxia-inducible factor 1 alpha (HIF1A) expression is also associated with exercise-related improvements in inflammation; increased lactate can trigger increased HIF1A [43], where HIF1A can then stimulate immune cell IL-1β and IL-17 production [44].

TCA cycle flux and oxidative metabolism

Increased RA muscle glycolysis appears to be connected to remodeling of oxidative metabolic pathways. Oxidative phosphorylation requires reducing agents (i.e., NADH and FADH2), typically generated through acetyl CoA and TCA cycle metabolism. Our findings suggest that RA muscle amino acid sources contribute more and glucose and fatty acids less than expected to acetyl CoA generation (Fig. 2). Remodeling of RA muscle oxidative metabolism was also evidenced by the downregulation of multiple TCA cycle enzymes (i.e., OGDH, SUCLA2, MDH1B). This apparent reduction in beta-oxidation is congruent with findings of increased lipid storage (i.e., intramuscular adiposity) in RA muscle. Intramuscular adiposity is linked to disability and an early aging phenotype in RA [45, 46]. Taken together, this suggests that exercise training is potentially important for modulating both RA fat metabolism and functional impairments.

Immune pathways

In addition to metabolic alterations, we hypothesized that upregulated RA skeletal muscle immune/inflammatory cytokine pathways would be associated with disease activity and improvements in disease activity following exercise training; however, our data did not consistently support this theory. In the cross-sectional RA cohort #1, IL1RL2 and TAB1 were the only skeletal muscle genes specifically involved in immune function or inflammation to correlate highly with RA disease activity (Supplementary Table 1). Skeletal muscle IL-15 production was the top canonical pathway associated with disease activity (12/128 genes; p = 0.12; data not shown). Interestingly, muscle IL-15 has important roles in regulating fat mass and systemic metabolism, as well as modulating lymphocyte development and acute inflammation [11]. In the HIIT RA cohort #2, only muscle expression of immune-related genes FCRL6, TNFRSF19, CMTM4, and NKG7 were highly correlated with improvements in RA disease activity following HIIT. In established RA, these data suggest that systemic inflammation impacts skeletal muscle metabolism to a greater extent than the direct effects of localized inflammation within skeletal muscle tissue.


Though we have identified multiple skeletal muscle pathways that are connected to RA inflammation and can potentially be modulated by exercise training, our findings should be reviewed in the context of a few key limitations. Primarily, the bulk of our analyses relied on correlations based on stored human tissue samples; and thus, true causative pathways or intervenable targets could not be interrogated. Further, in the HIIT RA cohort #2, we analyzed correlations between baseline/pre-HIIT factors—as opposed to changes in these factors—with changes in RA disease activity following exercise training. While this decision to focus on baseline factors may limit our understanding of co-occurring exercise training-induced muscle and immune modifications, our analyses better show how immune overactivity influences RA muscle at baseline and the potential for exercise training to rectify those interactions. The longitudinal HIIT RA cohort analyses are also limited by the small sample size of 12 RA participants and the lack of a non-RA control group completing the HIIT program; thus, it is unclear if the findings of the study are unique to RA or not. Finally, our HIIT cohort did not include a control RA group (i.e., not undergoing exercise training) which somewhat limits the interpretation of exercise training-specific effects on RA muscle and immune function. Thus, further detailed study of cross-talk between RA skeletal muscle and immune cells is necessary to dissect the intricacies of these pathways so they can eventually be exploited to improve patient care.


High-intensity interval training yielded the greatest improvements in inflammatory disease activity for persons with RA who were older, less aerobically fit, and with greater systemic inflammation. We identified alterations in multiple RA skeletal muscle metabolic pathways that associate strongly with disease activity as well as exercise training-related improvements in disease activity. Our data suggest that transcriptional remodeling of RA muscle amino acid, glucose, and fat metabolism is intricately linked to chronic inflammation. Exercise training may reprogram these muscle metabolic pathways to help maintain muscle-immune bioenergetic balance and modulate systemic immune responses. Ultimately, our findings provide further evidence for connections between impaired skeletal muscle and immune function in chronic inflammatory diseases such as RA and the ability of exercise training to rectify these underappreciated interorgan maladaptations.

Availability of data and materials

The raw data supporting the conclusions of this article will be made available by the authors upon request without undue reservation. Illumina array data will be deposited to the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO).


  1. Wadley AJ, van Zanten JJV, Stavropoulos-Kalinoglou A, Metsios GS, Smith JP, Kitas GD, et al. Three months of moderate-intensity exercise reduced plasma 3-nitrotyrosine in rheumatoid arthritis patients. Eur J Appl Physiol. 2014;114(7):1483–92.

    Article  CAS  PubMed  Google Scholar 

  2. Hakkinen A, Sokka T, Kotaniemi A, Hannonen P. A randomized two-year study of the effects of dynamic strength training on muscle strength, disease activity, functional capacity, and bone mineral density in early rheumatoid arthritis. Arthritis Rheum. 2001;44(3):515–22.<515::AID-ANR98>3.0.CO;2-5.

    Article  CAS  PubMed  Google Scholar 

  3. Stavropoulos-Kalinoglou A, Metsios GS, van Zanten JJV, Nightingale P, Kitas GD, Koutedakis Y. Individualised aerobic and resistance exercise training improves cardiorespiratory fitness and reduces cardiovascular risk in patients with rheumatoid arthritis. Ann Rheum Dis. 2013;72(11):1819–25.

    Article  CAS  PubMed  Google Scholar 

  4. van den Ende CH, Breedveld FC, le Cessie S, Dijkmans BA, de Mug AW, Hazes JM. Effect of intensive exercise on patients with active rheumatoid arthritis: a randomised clinical trial. Ann Rheum Dis. 2000;59(8):615–21.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Bartlett DB, Willis LH, Slentz CA, Hoselton A, Kelly L, Huebner JL, et al. Ten weeks of high-intensity interval walk training is associated with reduced disease activity and improved innate immune function in older adults with rheumatoid arthritis: a pilot study. Arthritis Res Ther. 2018;20(1):127.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Weinhold M, Shimabukuro-Vornhagen A, Franke A, Theurich S, Wahl P, Hallek M, et al. Physical exercise modulates the homeostasis of human regulatory T cells. J Allergy Clin Immunol. 2016;137(5):1607–10 e8.

    Article  PubMed  Google Scholar 

  7. Duggal NA, Pollock RD, Lazarus NR, Harridge S, Lord JM. Major features of immunesenescence, including reduced thymic output, are ameliorated by high levels of physical activity in adulthood. Aging Cell. 2018;17(2):e12750. Epub 2018 Mar 8.

  8. Severinsen MCK, Pedersen BK. Muscle-Organ crosstalk: the emerging roles of myokines. Endocr Rev. 2020;41(4):594-609.

  9. Pedersen BK, Febbraio MA. Muscle as an endocrine organ: focus on muscle-derived interleukin-6. Physiol Rev. 2008;88(4):1379–406.

    Article  CAS  PubMed  Google Scholar 

  10. Gudiksen A, Schwartz CL, Bertholdt L, Joensen E, Knudsen JG, Pilegaard H. Lack of skeletal muscle IL-6 affects pyruvate dehydrogenase activity at rest and during prolonged exercise. PLoS One. 2016;11(6):e0156460.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Leal LG, Lopes MA, Batista ML Jr. Physical exercise-induced myokines and muscle-adipose tissue crosstalk: a review of current knowledge and the implications for health and metabolic diseases. Front Physiol. 2018;9:1307.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Pal M, Febbraio MA, Whitham M. From cytokine to myokine: the emerging role of interleukin-6 in metabolic regulation. Immunol Cell Biol. 2014;92(4):331–9.

    Article  CAS  PubMed  Google Scholar 

  13. Jin JO, Han X, Yu Q. Interleukin-6 induces the generation of IL-10-producing Tr1 cells and suppresses autoimmune tissue inflammation. J Autoimmun. 2013;40:28–44.

    Article  CAS  PubMed  Google Scholar 

  14. Starkie R, Ostrowski SR, Jauffred S, Febbraio M, Pedersen BK. Exercise and IL-6 infusion inhibit endotoxin-induced TNF-alpha production in humans. FASEB J. 2003;17(8):884–6.

    Article  CAS  PubMed  Google Scholar 

  15. Andonian BJ, Huffman KM. Skeletal muscle disease in rheumatoid arthritis: the center of cardiometabolic comorbidities? Curr Opin Rheumatol. 2020;32(3):297–306.

    Article  PubMed  Google Scholar 

  16. Huffman KM, Jessee R, Andonian B, Davis BN, Narowski R, Huebner JL, et al. Molecular alterations in skeletal muscle in rheumatoid arthritis are related to disease activity, physical inactivity, and disability. Arthritis Res Ther. 2017;19(1):12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Andonian BJ, Bartlett DB, Huebner JL, Willis L, Hoselton A, Kraus VB, et al. Effect of high-intensity interval training on muscle remodeling in rheumatoid arthritis compared to prediabetes. Arthritis Res Ther. 2018;20(1):283.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Boutrup RJ, Farup J, Vissing K, Kjaer M, Mikkelsen UR. Skeletal muscle stem cell characteristics and myonuclei content in patients with rheumatoid arthritis: a cross-sectional study. Rheumatol Int. 2018;38(6):1031–41.

    Article  CAS  PubMed  Google Scholar 

  19. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.

    Article  CAS  PubMed  Google Scholar 

  20. Prevoo ML, van ‘t Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, van Riel PL. Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum. 1995;38(1):44–8.

    Article  CAS  PubMed  Google Scholar 

  21. Huffman KM, Pieper CF, Hall KS, St Clair EW, Kraus WE. Self-efficacy for exercise, more than disease-related factors, is associated with objectively assessed exercise time and sedentary behaviour in rheumatoid arthritis. Scand J Rheumatol. 2015;44(2):106–10.

    Article  CAS  PubMed  Google Scholar 

  22. Andonian BJ, Chou CH, Ilkayeva OR, Koves TR, Connelly MA, Kraus WE, et al. Plasma microRNAs in established rheumatoid arthritis relate to adiposity and altered plasma and skeletal muscle cytokine and metabolic profiles. Front Immunol. 2019;10:1475.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;47(D1):D419–D26.

    Article  CAS  PubMed  Google Scholar 

  24. Thomas PD, Kejariwal A, Guo N, Mi H, Campbell MJ, Muruganujan A, et al. Applications for protein sequence-function evolution data: mRNA/protein expression analysis and coding SNP scoring tools. Nucleic Acids Res. 2006;34(Web Server issue):W645–50.

    Article  Google Scholar 

  25. Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498–503.

    Article  CAS  PubMed  Google Scholar 

  26. AbouAssi H, Tune KN, Gilmore B, Bateman LA, McDaniel G, Muehlbauer M, et al. Adipose depots, not disease-related factors, account for skeletal muscle insulin sensitivity in established and treated rheumatoid arthritis. J Rheumatol. 2014;41(10):1974–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xu J, Jiang C, Cai Y, Guo Y, Wang X, Zhang J, et al. Intervening upregulated SLC7A5 could mitigate inflammatory mediator by mTOR-P70S6Ksignal in rheumatoid arthritis synoviocytes. Arthritis Res Ther. 2020;22(1):200.

  28. Newsholme EA, Parry-Billings M. Properties of glutamine release from muscle and its importance for the immune system. JPEN J Parenter Enteral Nutr. 1990;14(4 Suppl):63S–7S.

    Article  CAS  PubMed  Google Scholar 

  29. Labow BI, Souba WW, Abcouwer SF. Mechanisms governing the expression of the enzymes of glutamine metabolism--glutaminase and glutamine synthetase. J Nutr. 2001;131(9 Suppl):2467S–74S discussion 86S-7S.

    Article  CAS  Google Scholar 

  30. Kelly B, Pearce EL. Amino assets: how amino acids support immunity. Cell Metab. 2020;32(2):154–75.

    Article  CAS  PubMed  Google Scholar 

  31. Kono M, Yoshida N, Maeda K, Suarez-Fueyo A, Kyttaris VC, Tsokos GC. Glutaminase 1 inhibition reduces glycolysis and ameliorates lupus-like disease in MRL/lpr mice and experimental autoimmune encephalomyelitis. Arthritis Rheumatol. 2019;71(11):1869–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Palsson-McDermott EM, O'Neill LAJ. Targeting immunometabolism as an anti-inflammatory strategy. Cell Res. 2020;30(4):300–14.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Alves A, Bassot A, Bulteau AL, Pirola L, Morio B. Glycine metabolism and its alterations in obesity and metabolic diseases. Nutrients. 2019;11(6):1356.

  34. UniProt C. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47(D1):D506–D15.

    Article  Google Scholar 

  35. Kikuchi G, Motokawa Y, Yoshida T, Hiraga K. Glycine cleavage system: reaction mechanism, physiological significance, and hyperglycinemia. Proc Jpn Acad Ser B Phys Biol Sci. 2008;84(7):246–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Quemeneur L, Gerland LM, Flacher M, Ffrench M, Revillard JP, Genestier L. Differential control of cell cycle, proliferation, and survival of primary T lymphocytes by purine and pyrimidine nucleotides. J Immunol. 2003;170(10):4986–95.

    Article  CAS  PubMed  Google Scholar 

  37. Yin J, Ren W, Huang X, Deng J, Li T, Yin Y. Potential mechanisms connecting purine metabolism and cancer therapy. Front Immunol. 2018;9:1697.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Cruzat V, Macedo Rogero M, Noel Keane K, Curi R, Newsholme P. Glutamine: metabolism and immune function, supplementation and clinicaltranslation. Nutrients. 2018;10(11):1564.

  39. Percival JM. nNOS regulation of skeletal muscle fatigue and exercise performance. Biophys Rev. 2011;3(4):209–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Shime H, Yabu M, Akazawa T, Kodama K, Matsumoto M, Seya T, et al. Tumor-secreted lactic acid promotes IL-23/IL-17 proinflammatory pathway. J Immunol. 2008;180(11):7175–83.

    Article  CAS  PubMed  Google Scholar 

  41. Pucino V, Certo M, Bulusu V, Cucchi D, Goldmann K, Pontarini E, et al. Lactate buildup at the site of chronic inflammation promotes disease by inducing CD4(+) T cell metabolic rewiring. Cell Metab. 2019;30(6):1055–74 e8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Haas R, Smith J, Rocher-Ros V, Nadkarni S, Montero-Melendez T, D'Acquisto F, et al. Lactate regulates metabolic and pro-inflammatory circuits in control of T cell migration and effector functions. PLoS Biol. 2015;13(7):e1002202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. De Saedeleer CJ, Copetti T, Porporato PE, Verrax J, Feron O, Sonveaux P. Lactate activates HIF-1 in oxidative but not in Warburg-phenotype human tumor cells. PLoS One. 2012;7(10):e46571.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Talreja J, Talwar H, Bauerfeld C, Grossman LI, Zhang K, Tranchida P, et al. HIF-1alpha regulates IL-1beta and IL-17 in sarcoidosis. Elife. 2019;8.

  45. Baker JF, Mostoufi-Moab S, Long J, Zemel B, Ibrahim S, Taratuta E, et al. Intramuscular fat accumulation and associations with body composition, strength, and physical functioning in patients with rheumatoid arthritis. Arthritis Care Res (Hoboken). 2018;70(12):1727–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Khoja SS, Patterson CG, Goodpaster BH, Delitto A, Piva SR. Skeletal muscle fat in individuals with rheumatoid arthritis compared to healthy adults. Exp Gerontol. 2020;129:110768.

    Article  PubMed  Google Scholar 

Download references


Original figure art was created with We appreciate the support of the Duke University Division of Rheumatology and Immunology. We acknowledge Leslie Willis and the Duke Center for Living staff members for their help with participant recruitment, exercise training, and recording of data. We thank Milton Campbell for his assistance with tissue RNA preparation. We thank the Duke Center for Human Genetics for their assistance with molecular analyses. We acknowledge Deborah Muoio for her assistance with data analysis interpretation. We acknowledge and appreciate greatly all participants in the study.


This research was supported by the Duke Pepper Center REC Career Development Award (to BJA), the National Institute on Aging of the National Institutes of Health R03-AG067949 (to BJA), the National Institute of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of Health R21-AR076663 (to KMH), and an EU Marie Curie Outgoing Fellowship Award PIOF-GA-2013-629981 (to DBB). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Author information

Authors and Affiliations



BJA, AJ, WEK, DBB, and KMH designed the study and experimental approach. KMH, DBB, DMP, and MJH acquired data and completed laboratory studies. BJA, AJ, AK, DMP, and MJH conducted data analyses. BJA, AJ, and AK contributed to manuscript writing. The authors contributed to data interpretation and reviewed, edited, and approved the manuscript.

Corresponding author

Correspondence to Brian J. Andonian.

Ethics declarations

Ethics approval and consent to participate

All participants gave written informed consent. The study was approved by the Duke University Medical Center Institutional Review Board (IRB no. Pro00064057).

Consent for publication

Not applicable.

Competing interests

Author DMP is an employee at Phastar Inc. (Cambridge, MA, USA). All other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1.

Cross-sectional rheumatoid arthritis cohort #1 top skeletal muscle genes. Supplementary Table 2. Cross-sectional rheumatoid arthritis cohort #1: Skeletal muscle canonical pathways associated with disease activity. Supplementary Table 3. High-intensity interval training rheumatoid arthritis cohort #2 top skeletal muscle genes associated with improvements in disease activity. Supplementary Table 4. High-intensity interval training rheumatoid arthritis cohort #2: Skeletal muscle canonical pathways associated with improvements in disease activity.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Andonian, B.J., Johannemann, A., Hubal, M.J. et al. Altered skeletal muscle metabolic pathways, age, systemic inflammation, and low cardiorespiratory fitness associate with improvements in disease activity following high-intensity interval training in persons with rheumatoid arthritis. Arthritis Res Ther 23, 187 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: