Skip to main content

Molecular alterations in skeletal muscle in rheumatoid arthritis are related to disease activity, physical inactivity, and disability

Abstract

Background

To identify molecular alterations in skeletal muscle in rheumatoid arthritis (RA) that may contribute to ongoing disability in RA.

Methods

Persons with seropositive or erosive RA (n = 51) and control subjects matched for age, gender, race, body mass index (BMI), and physical activity (n = 51) underwent assessment of disease activity, disability, pain, physical activity and thigh muscle biopsies. Muscle tissue was used for measurement of pro-inflammatory markers, transcriptomics, and comprehensive profiling of metabolic intermediates. Groups were compared using mixed models. Bivariate associations were assessed with Spearman correlation.

Results

Compared to controls, patients with RA had 75% greater muscle concentrations of IL-6 protein (p = 0.006). In patients with RA, muscle concentrations of inflammatory markers were positively associated (p < 0.05 for all) with disease activity (IL-1β, IL-8), disability (IL-1β, IL-6), pain (IL-1β, TNF-α, toll-like receptor (TLR)-4), and physical inactivity (IL-1β, IL-6). Muscle cytokines were not related to corresponding systemic cytokines. Prominent among the gene sets differentially expressed in muscles in RA versus controls were those involved in skeletal muscle repair processes and glycolytic metabolism. Metabolic profiling revealed 46% higher concentrations of pyruvate in muscle in RA (p < 0.05), and strong positive correlation between levels of amino acids involved in fibrosis (arginine, ornithine, proline, and glycine) and disability (p < 0.05).

Conclusion

RA is accompanied by broad-ranging molecular alterations in skeletal muscle. Analysis of inflammatory markers, gene expression, and metabolic intermediates linked disease-related disruptions in muscle inflammatory signaling, remodeling, and metabolic programming to physical inactivity and disability. Thus, skeletal muscle dysfunction might contribute to a viscous cycle of RA disease activity, physical inactivity, and disability.

Background

Despite a vast array of pharmacologic agents available to treat rheumatoid arthritis (RA), management is often complicated by insufficient treatment response, drug toxicity and contraindications, poor access to care and/or medications, and/or damage that predates medical intervention. These barriers lead to or are accompanied by systemic manifestations, disease-associated co-morbidities, chronic pain, physical inactivity, dysmobility, and poor physical function. Thus, further advances in RA care require identification of factors contributing to persistent deficiencies in quality of life and physical function, despite access to excellent anti-rheumatic medications.

Importantly, inactivity and muscle wasting are two important contributors to RA-related morbidity and mortality. Approximately half of patients with RA do not perform even a single bout of weekly physical exercise [1]. The sedentary lifestyle common to patients with RA gives rise to physical deconditioning and muscle atrophy, both of which are associated with osteoporosis, impaired immune function, glucose intolerance, insulin resistance, loss of independence, and increased mortality [2].

In addition to physical inactivity, other factors that likewise promote muscle loss and disability in patients with RA include inadequate protein ingestion, glucocorticoid treatment, and pro-inflammatory cytokines, all resulting in reduced myocyte protein synthesis and increased protein degradation [2, 3]. Inflammation can impact normal muscle turnover and responses to injury, both of which require an exquisitely coordinated remodeling process involving activation, proliferation and differentiation of muscle stem cells—also known as satellite cells. These processes are mediated largely by signals from intramuscular immune cells: neutrophils, regulatory T cells, pro-inflammatory M1 macrophages, and anti-inflammatory M2 macrophages.

The established roles of inflammation in both skeletal muscle remodeling and RA pathophysiology raise obvious questions regarding the potential interplay between muscle dysfunction and RA morbidity. Whereas the link between pro-inflammatory cytokines and muscle dysfunction has been investigated intensely in the context of diseases such as diabetes and cancer cachexia, this topic has remained surprisingly unexplored in RA. In the current study we sought to identify molecular perturbations in muscle specimens from individuals with RA, and to test the hypothesis that skeletal muscle inflammatory markers and derangements in tissue remodeling might contribute to metabolic decline and disability in these patients. Herein, we report that disease-activity-related muscle inflammatory markers are related to physical inactivity, and moreover, that disrupted skeletal muscle repair processes are associated with greater disability. These findings support a model in which skeletal muscle deterioration contributes to a vicious cycle of disease activity, muscle inflammatory signaling and disrupted remodeling, physical inactivity, and disability in patients with RA.

Methods

Design and participants

This was a cross-sectional investigation of individuals with RA and matched controls collected from the Durham, NC area. The RA group met the following criteria: (1) RA diagnosis meeting American College of Rheumatology 1987 criteria [4]; (2) seropositive disease (positive rheumatoid factor or anti-cyclic citrullinated peptide) or evidence of erosions on hand or foot imaging; (3) no medication changes within the three months prior to study enrollment; and (4) daily prednisone use ≤5 mg. Healthy participants without a diagnosis of RA, without joint pain, and without joint swelling lasting more than a week were matched to individual participants with RA by gender, race, age within 3 years, and body mass index (BMI) within 3 kg/m2. Exclusions included current pregnancy, type 2 diabetes mellitus, and known coronary artery disease. Further specific details on questionnaires and measurement protocols have previously been described [5]. This study was in compliance with the Helsinki Declaration and was approved by the Duke University Institutional Review Board.

Assessments of both groups included questionnaires, physical exams for disease status, fasting blood collection, intravenous glucose tolerance tests for insulin sensitivity, 7 days of accelerometer-measured physical activity, computed tomography (CT) imaging of abdomen and thigh, and vastus lateralis muscle biopsies [5]. Disability (health assessment questionnaire-disability index (HAQ-DI) and co-morbidities (co-morbidity index) were assessed by previously published questionnaires [6, 7]. Disease activity assessed by the disease activity score in 28 joints (DAS-28) was determined from a patient-completed visual analog scale, physician-determined numbers of tender and swollen joints, and erythrocyte sedimentation rate [8]. Plasma concentrations of inflammatory markers and cytokines were determined by immunoassay [5] and nuclear magnetic resonance (NMR) spectroscopy (GlycA) [9]. Insulin sensitivity was determined using Bergman’s minimal model [10] and concentrations of glucose and insulin (glucose: Beckman-CoulterDXC600; insulin: electrochemiluminscent assay from Meso Scale Discovery) at each of 29 time points during the intravenous glucose tolerance test.

Physical activity was measured with 7 days of accelerometry. After completing assessments, accelerometers (RT3, Stayhealthy, Inc., Monrovia, CA, USA) were provided to participants. Participants also received a pre-addressed and postage-applied box for return and directions for wearing on the waist above the right knee during waking hours for 7 days. Accelerometer data were evaluated for validity and non-wear time, and categorized into metabolic equivalents (METs) as previously described [11]. After data cleaning, valid data were available for 41 persons with RA and 31 controls. Time spent exercising was defined as the sum of time spent performing activity at METs equal to or greater than 3. CT scan analyses were performed using OsiriX (Pixmeo) to determine adipose and muscle tissue size and muscle tissue density (greater tissue density is indicative of less inter-muscular adipose tissue) [5]. Standard Bergstrom needle muscle biopsies were performed on the vastus lateralis in the fasting state; participants consumed only water during the 12 hours overnight prior to the biopsy [12]. Tissue was flash frozen in liquid nitrogen and stored at −80 ° C.

Skeletal muscle inflammatory marker measurements

Flash frozen muscle samples (5–10 mg) were homogenized in a buffer consisting of 1% Nonidet-P40, 1 mM EDTA, 150 mM NaCl, and 20 mM Tris-Cl for ELISA-based measures of muscle (m) interleukin (IL)-1β, mIL-6, mIL-8, m-tumor necrosis factor (TNF)-α (MSD 4-plex; K15053D-1) and m-toll like receptor (TLR)-4 (Abnova; KA1238). Assays were performed according to the manufacturers’ directions except for the addition of a 30-minute, room temperature, blocking step with 5% BSA followed by three PBS-T washes. Concentrations were normalized to starting masses. Spike-and-recovery assays for all analytes achieved 80–100% recovery confirming lack of assay interference by muscle homogenates. For each cytokine, the mean intra-assay and inter-assay coefficients of variation were: mIL-1β 8.5%, 13.2%; mIL-6 3.5%, 1.5%; mIL-8 4.0%, 4.0%; mTNF-α 8.4%, 10.4%; and mTLR-4 1.7% (only one plate was used for analyses).

Gene expression analyses

Muscle samples were selected for gene expression analyses in an effort to span the range of RA disease activity seen in the larger sample; these corresponded to the following DAS-28 categories: remission (n = 7), low (n = 4), moderate (n = 6), and high activity (n = 3). For each RA muscle sample, the corresponding sample from a control matched by age, gender, and BMI was included.

For RNA preparation, flash frozen muscle samples (20–30 mg) were homogenized in 1 mL TRIzol® (Thermo Fisher Scientific, Inc, Waltham, MA, USA). Biotinylated total RNA was generated using the Illumina TotalPrep RNA amplification kit (Life Technologies, Grand Island, NY, USA); 200 nanograms of RNA were used for the kit. The quality of the RNA was determined using the Bioanalyzer RNA Nano chip assay (Agilent, Santa Clara, CA, USA). Quantification of the RNA was determined using the Quant-iT RiboGreen RNA Assay Kit. The Human HT-12v3 Expression BeadChip (Illumina, San Diego, CA, USA) was used for quantitative whole genome RNA profiling. Biotinylated RNA (750 ng) was hybridized to the BeadChip and washed; Cy3-SA was then introduced to the hybridized samples and the BeadChips scanned on the Illumina iScan system according to manufacturer’s protocol. Quality control was performed using the Illumina GenomeStudio tools.

Gene expression fold-differences between groups were compared in Partek Genomics Suite (Partek, Inc.; St. Louis, MO, USA). For pathway analyses, differentially expressed genes (p < 0.02) were evaluated using the Ingenuity Pathway Analysis software (IPA, www.ingenuity.com). IPA identified the canonical pathways containing the greatest number of significant, differentially expressed genes in the dataset. IPA also generated novel networks of related genes and molecules based on the relationships present in the current literature.

Skeletal muscle metabolic intermediate measurements

Metabolites were measured in muscle from all participants (n = 102). Flash frozen muscle biopsies weighing approximately 25 mg were diluted 20 times (wt:vol) in ice-cold 50% acetonitrile containing 0.3% formate and homogenized for 120 sec in a TissueLyser II (Qiagen) at 30 Hz. Amino acids, organic acids, and acylcarnitines were analyzed using stable isotope dilution techniques in the Duke Molecular Physiology Metabolomics Core. Amino acids and acylcarnitine measurements were made by flow injection tandem mass spectrometry (MS) as previously described [13, 14]. The data were acquired using a Micromass Quattro Micro liquid chromatography (LC)-MS system running MassLynx 4.0 software (Waters Corporation, Milford, MA, USA). Organic acids were quantified using methods described previously [15] employing Trace Ultra GC coupled to ISQ MS operating under Xcalibur 2.2 (Thermo Fisher Scientific, Austin, TX, USA). All data are expressed as picomoles/mg tissue.

Statistical analyses

Accounting for the repeated measures in matched participants, patients with RA and controls were compared using mixed models. Muscle inflammatory molecules and metabolic intermediates were logarithmically transformed prior to group comparisons. Bivariate associations were assessed with Spearman correlation. Gene expression fold-changes were compared in Partek using analysis of variance (ANOVA). All other statistical analyses were performed using SAS 9.4 (SAS, Cary, NC). All data are available from the corresponding author upon reasonable request.

Results

Clinical measures and skeletal muscle inflammatory markers

As shown in Table 1, persons with RA were well-matched to controls by age, gender, and BMI. Patients with RA were recruited based on the inclusion criteria described and without respect to physical activity levels, body mass or body composition; similarly controls were included upon matching a patient with RA by age, gender, and BMI. Despite this, patients with RA and controls were similar with respect to physical activity levels, abdominal and thigh adipose depot size, muscle area, and muscle density [5, 11]. In those with RA, there was more co-morbidity, disability, and systemic inflammation; specifically, greater serum concentrations of high sensitivity C-reactive protein (hs-CRP), IL-6, and TNF-α (p < 0.05 for all) [5]. When skeletal muscle inflammatory markers were compared, there was approximately two times greater concentrations of the muscle cytokines, mIL-6 (p = 0.006) and mIL-8 in RA (p = 0.059) (Table 1).

Table 1 Participant characteristics

Akin to disease activity, RA muscle inflammatory markers exhibited variation across a broad range (Table 1). Muscle inflammatory marker concentrations were positively associated with disease activity (mIL-1β, mIL-8), disability (mIL-1β, mIL-6), and pain (mIL-1β, mTNF-α, mTLR-4) (p < 0.05 for all) (Table 2). Muscle cytokines, mIL-1β and mIL-8, were negatively correlated with use of biological agents; mTNF-α was negatively correlated with use of non-biological disease-modifying therapy (p < 0.05 for all) (Table 2). Importantly, there were no correlation between muscle inflammatory marker concentrations and prednisone treatment.

Table 2 Skeletal muscle inflammatory marker correlations in patients with rheumatoid arthritis

In addition to disease-related factors, muscle cytokine concentrations (mIL-1β, mIL-6, and mIL-8) were negatively associated with exercise minutes (p < 0.05 for all) (Table 2). Higher mIL-1β and mIL-6 concentrations were associated with less total physical activity (total METs; p < 0.05 for both) (Table 2). Although, both mIL-6 and mIL-8 were positively correlated with the systemic inflammatory marker, GlycA (p < 0.05 for both) (Table 2), muscle inflammatory marker concentrations were not related to insulin sensitivity or systemic cytokine concentrations.

Skeletal muscle gene expression

To better understand the etiology of RA muscle inflammatory markers, we compared RA (n = 20) and control (n = 20) skeletal muscle gene expression: 1939 genes were significantly upregulated or downregulated in RA samples (p < 0.05); 445 genes were identified when using a more stringent definition of significance (p < 0.02).

To identify other relationships between differentially expressed RA muscle genes, pathway analyses were performed using IPA, which has thousands of canonical pathways onto which our experimental gene expression differences were overlaid. Of those canonical pathways, IPA identified several pathways impacted by differential gene expression in muscle in RA (p < 0.05) (Table 3). Except for glycolysis and methionine degradation, these canonical pathways were identified because of reduced RA muscle gene expression for nuclear factor (NF)-kβ2, both nuclear factor of activated T cells (NFAT)5 and NFATC4, or all three. Also, none of the canonical pathways was predicted to be activated or inhibited by gene expression differences in muscle in RA (Z-scores < |2|) (Table 3) [16].

Table 3 Canonical pathways implicated in gene expression in muscle in rheumatoid arthritis

In addition to canonical networks, pathway analyses generate novel networks connecting differentially regulated molecules based on published associations. The IPA-generated novel network with the highest connection score depicted significant differences in expression of genes associated with connective tissue, dental, and dermatological diseases (Fig. 1; Table 4). The prominent pathway connections in muscle in RA were centered on regulation of the NF-kB complex, specifically NF-kB2. These were in the setting of differential regulation of genes in muscle repair and glycolysis.

Fig. 1
figure1

Novel network identified by muscle gene expression in rheumatoid arthritis (RA): gene expression was determined in muscle from 20 persons with RA and 20 controls matched by age, gender, and body mass index. The network shows connections between genes with differential expression in RA relative to control muscle. Genes in red were upregulated and genes in green were downregulated in muscle in RA

Table 4 Novel network genes

To augment traditional pathway analyses, we evaluated the 20 genes with the largest muscle expression differences in RA and control samples (Table 5) and examined gene members of well-established skeletal muscle anabolic, catabolic, and inflammatory pathways (Table 6). The top 20 upregulated and downregulated genes by fold difference were associated with muscle remodeling, satellite cell maturation, exercise intolerance, and/or energy metabolism; for these genes, the range of differences in expression was 20–50% (Table 5). Except for NF-kB2, there was no differential expression of canonical genes involved in skeletal muscle anabolic, catabolic, or inflammatory pathways (Table 6).

Table 5 Genes with the greatest differences in expression between patients with rheumatoid arthritis (RA) and controls
Table 6 Genes involved in skeletal muscle anabolic, catabolic, and inflammatory pathways

Skeletal muscle metabolic intermediates

When concentrations of skeletal muscle metabolic intermediates were compared between RA (n = 51) and controls (n = 51), muscle pyruvate concentrations were 46% greater in muscle in RA than in controls (p < 0.001) (Table 7). There were no significant differences in the concentrations of muscle amino acids, other organic acids, or acylcarnitines in RA compared to controls (Table 7). However, several muscle amino acids and acylcarnitines were related to RA disease activity and disability. For instance, greater concentrations of glycine, serine, aspartate/asparagine, and ornithine and lower muscle concentrations of alanine and fumarate were related to greater disease activity (p < 0.05) (Table 8). Greater muscle concentrations of glycine, proline, ornithine, arginine, and aspartate/asparagine were related to greater disability (p < 0.05) (Fig. 2); in contrast, lower concentrations of several long-chain unsaturated acylcarnitines were related to greater disease activity and disability (p < 0.05) (Table 8).

Table 7 Skeletal muscle metabolic intermediate concentrations
Table 8 Relationships between rheumatoid arthritis clinical features and muscle metabolic intermediates
Fig. 2
figure2

Schematic depiction of muscle injury repair showing potential impact of cytokine, gene expression, and amino acids on satellite cell activation, macrophage function, and fibrosis in muscle from patients with rheumatoid arthritis (RA). Boxes show gene IDs for genes differentially regulated in patients with RA and in controls. See Table 3 for gene descriptions

Discussion

Here, we report that in RA, skeletal muscle exhibits molecular alterations in inflammatory markers, transcriptional profiles, and metabolic signatures. Both at protein and transcriptional levels, muscle had a pro-inflammatory phenotype in RA. Additionally, differential gene expression in muscle in RA was indicative of dysregulation of muscle repair, promotion of glycolysis, and poor mitochondrial function. Upregulated glycolysis and mitochondrial inefficiency were supported by greater muscle concentrations of the glycolytic end-product pyruvate in RA. Further, disease activity and disability were related to lesser concentrations of long-chain acylcarnitines and greater concentrations of amino acid precursors for muscle fibrosis. Taken together, these alterations in proteins, gene expression, and metabolic intermediates were indicative of muscle in RA in a state of chronically activated, yet dysregulated remodeling with increased glycolysis, mitochondrial inefficiency, and fibrotic material (Fig. 2).

This represents the first report of significant markers of inflammation in muscle in RA. The clinical importance of these molecules is demonstrated by the significant association of several muscle cytokines with RA disease activity, disability, pain, and physical inactivity. The IPA-generated novel network centered on downregulation of NF-kB2, a protein that promotes non-canonical NF-kB signaling and opposes inflammatory signaling [17]. Downregulation of NF-kB2 would be predicted to favor coordinated upregulation of pro-inflammatory NF-kB signaling in muscle in RA. We were unable to determine if the muscle cytokines and pro-inflammatory transcripts in RA were derived from myocytes, inflammatory cells, or other cellular sources. Interestingly, muscle cytokine concentrations did not reflect those measured in circulation, suggesting these disease-associated inflammatory markers stem from local rather than systemic events.

Based on the strong relationships between muscle inflammatory markers and disability, pain and physical inactivity, we suspected that increased intramuscular cytokines may be indicative of a disrupted muscle remodeling process. In fact, muscle gene expression alterations in RA were consistent with promotion of satellite cell differentiation and upregulation of several facets of the normally well-coordinated process of muscle remodeling. For instance, muscle in RA was characterized by downregulation of ZFP36, the gene that encodes tristetraprolin (TTP), which reduces inflammation by destabilizing pro-inflammatory cytokine transcripts [18, 19] and prevents satellite cell activation by destabilizing myogenic regulatory factor, MyoD, mRNA [20]. Thus, the reduction in ZFP36 expression in muscle in RA would be expected to promote pro-inflammatory cytokine production and satellite cell activation.

Other gene expression changes also suggest both chronic activation and temporal dysregulation of muscle remodeling. For instance, downregulation of ZBTB16 would promote inflammation and proliferation of autoreactive T cells [21]. In contrast to the reduced ZFP36 expected to promote satellite cell activation, the increased CITED2 would be expected to reduce satellite cell activation [22, 23]. Increased expression of ABRA, RCAN1, VGLL2, MYF6 and decreased expression of RRAD would promote differentiation of satellite cells [2430]. More descriptions of differentially expressed genes are provided in Additional file 1.

Gene expression alterations indicative of glycolysis promotion and poor mitochondrial function were supported by greater muscle concentrations of the glycolytic end-product pyruvate in RA. Further, disease activity and disability were related to lower concentrations of fatty-acid-derived long-chain acylcarnitines. One plausible explanation for this relationship is that fewer long-chain acylcarnitines indicate less oxidative metabolism and fewer mitochondria, consistent with a glycolytic phenotype. RA disease activity and disability were also related to higher concentrations of amino acid precursors for muscle fibrosis. M2-type macrophages contain arginase, which metabolizes arginine to ornithine [31]. Ornithine is converted to proline, which provides a substrate for resident fibroblasts to generate collagen. In addition to proline, collagen formation also requires glycine; glycine and proline each account for a third of the collagen amino acids. While collagen is critical for extracellular matrix production, in the setting of a chronically activated remodeling process, excess collagen production leads to fibrosis [31]. Thus, the relationships between these amino acids and disease activity and disability may indicate a fibrotic process in muscle associated with active disease that contributes to RA-associated disability.

There were several limitations to this study. RA medication regimens were not uniform among participants, and effects of these medications on skeletal muscle are unclear. Twenty-five percent of patients with RA used prednisone at low doses, which is not expected to have significant myopathic effects; despite this, they had significant alterations in muscle inflammatory markers and systemic inflammation relative to controls. Without histopathologic assessment or single cell isolations, we were unable to determine the cellular source of muscle cytokines, transcripts, or metabolites. Our findings indicate that either RA medication regimens or the RA disease process itself alters skeletal muscle inflammatory molecules, transcriptional profiles, and metabolic pathways.

Conclusions

Taken together, these alterations in pro-inflammatory cytokines, gene expression, and metabolic intermediates are indicative of RA muscle in a state of chronically activated, yet dysregulated remodeling, with increased glycolysis, mitochondrial inefficiency, and fibrosis. It is very likely these changes contribute to the ongoing issues of exercise intolerance and disability in persons with RA. Future work should be directed at understanding whether these deficits may be mitigated by combining pharmacologic treatment with physical activity, to reduce inflammatory signaling and/or fibrosis while promoting skeletal muscle efficiency. Therefore, to improve the lives of patients with RA, future work should be directed toward understanding the role of skeletal muscle in RA, and interactions between treatment regimens, physical activity, and influences of skeletal muscle on the clinical status in RA.

Abbreviations

BMI:

body mass index

BSA:

bovine serum albumin

DAS-28:

disease activity score

DMARD:

disease modifying anti-rheumatic drugs

ELISA:

enzyme-linked immunosorbent assay

HAQ-DI:

health assessment questionnaire-disability index

HOMA:

homeostasis model assessment

hsCRP:

high sensitivity C-reactive protein

Hu:

Houndsfield units

IL:

interleukin

IPA:

Ingenuity Pathway Analysis

LC:

liquid chromatography

m:

muscle

MET:

metabolic equivalent

MS:

mass spectrometry

NaN:

not a number

NF-kβ:

nuclear factor-kβ

NMR:

nuclear magnetic resonance

PBS:

phosphate-buffered saline

RA:

rheumatoid arthritis

TLR:

toll-like receptor

TNF-α:

tumor necrosis factor-alpha

TTP:

tristetraprolin

VAS:

visual analog scale

References

  1. 1.

    Lee J, Dunlop D, Ehrlich-Jones L, Semanik P, Song J, Manheim L, Chang RW. Public health impact of risk factors for physical inactivity in adults with rheumatoid arthritis. Arthritis Care Res (Hoboken). 2012;64:488–93.

    Article  Google Scholar 

  2. 2.

    Kotler DP. Cachexia. Ann Intern Med. 2000;133:622–34.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Glass DJ. Skeletal muscle hypertrophy and atrophy signaling pathways. Int J Biochem Cell Biol. 2005;37:1974–84.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31:315–24.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    AbouAssi H, Tune KN, Gilmore B, Bateman LA, McDaniel G, Muehlbauer M, Huebner JL, Hoenig HM, Kraus VB, St Clair EW, et al. Adipose depots, not disease-related factors, account for skeletal muscle insulin sensitivity in established and treated rheumatoid arthritis. J Rheumatol. 2014;41:1974–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Rigler SK, Studenski S, Wallace D, Reker DM, Duncan PW. Co-morbidity adjustment for functional outcomes in community-dwelling older adults. Clin Rehabil. 2002;16:420–8.

    Article  PubMed  Google Scholar 

  7. 7.

    Bruce B, Fries JF. The Stanford health assessment questionnaire: dimensions and practical applications. Health Qual Life Outcomes. 2003;1:20.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    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:44–8.

  9. 9.

    Bartlett DB, Connelly MA, AbouAssi H, Bateman LA, Tune KN, Huebner JL, Kraus VB, Winegar DA, Otvos JD, Kraus WE, Huffman KM. A novel inflammatory biomarker, GlycA, associates with disease activity in rheumatoid arthritis and cardio-metabolic risk in BMI-matched controls. Arthritis Res Ther. 2016;18:86.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Bergman RN, Ider YZ, Bowden CR, Cobelli C. Quantitative estimation of insulin sensitivity. Am J Physiol. 1979;236:E667–77.

    CAS  PubMed  Google Scholar 

  11. 11.

    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:106–10.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Bergstrom J. Percutaneous needle biopsy of skeletal muscle in physiological and clinical research. Scand J Clin Lab Invest. 1975;35:609–16.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    An J, Muoio DM, Shiota M, Fujimoto Y, Cline GW, Shulman GI, Koves TR, Stevens R, Millington D, Newgard CB. Hepatic expression of malonyl-CoA decarboxylase reverses muscle, liver and whole-animal insulin resistance. Nat Med. 2004;10:268–74.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Wu JY, Kao HJ, Li SC, Stevens R, Hillman S, Millington D, Chen YT. ENU mutagenesis identifies mice with mitochondrial branched-chain aminotransferase deficiency resembling human maple syrup urine disease. J Clin Invest. 2004;113:434–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Jensen MV, Joseph JW, Ilkayeva O, Burgess S, Lu D, Ronnebaum SM, Odegaard M, Becker TC, Sherry AD, Newgard CB. Compensatory responses to pyruvate carboxylase suppression in islet beta-cells. Preservation of glucose-stimulated insulin secretion. J Biol Chem. 2006;281:22342–51.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Kramer A, Green J, Pollard Jr J, Tugendreich S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 2014;30:523–30. 510.1093/bioinformatics/btt1703. Epub 2013 Dec 1013.

  17. 17.

    Cubillos-Zapata C, Hernandez-Jimenez E, Toledano V, Esteban-Burgos L, Fernandez-Ruiz I, Gomez-Pina V, Del Fresno C, Siliceo M, Prieto-Chinchina P, Perez de Diego R, et al. NFkappaB2/p100 is a key factor for endotoxin tolerance in human monocytes: a demonstration using primary human monocytes from patients with sepsis. J Immunol. 2014;193:4195–202.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Geyer BC, Ben Ari S, Barbash S, Kilbourne J, Mor TS, Soreq H. Nicotinic stimulation induces Tristetraprolin over-production and attenuates inflammation in muscle. Biochim Biophys Acta. 2012;1823:368–78.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Carballo E, Lai WS, Blackshear PJ. Feedback inhibition of macrophage tumor necrosis factor-alpha production by tristetraprolin. Science. 1998;281:1001–5.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Hausburg MA, Doles JD, Clement SL, Cadwallader AB, Hall MN, Blackshear PJ, Lykke-Andersen J, Olwin BB. Post-transcriptional regulation of satellite cell quiescence by TTP-mediated mRNA decay. Elife. 2015;4:e03390. doi:10.7554/eLife.03390.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Liu TM, Lee EH, Lim B, Shyh-Chang N. Concise review: balancing stem cell self-renewal and differentiation with PLZF. Stem Cells. 2016;34:277–87.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kranc KR, Schepers H, Rodrigues NP, Bamforth S, Villadsen E, Ferry H, Bouriez-Jones T, Sigvardsson M, Bhattacharya S, Jacobsen SE, Enver T. Cited2 is an essential regulator of adult hematopoietic stem cells. Cell Stem Cell. 2009;5:659–65.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Bhattacharya S, Michels CL, Leung MK, Arany ZP, Kung AL, Livingston DM. Functional role of p35srj, a novel p300/CBP binding protein, during transactivation by HIF-1. Genes Dev. 1999;13:64–75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Lamon S, Wallace MA, Russell AP. The STARS signaling pathway: a key regulator of skeletal muscle function. Pflugers Arch. 2014;466:1659–71.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Rybkin II OM, Copeland V, Czubryt MP, Shelton JM, van Rooij E, Richardson JA, Hill JA, De Windt LJ, Bassel-Duby R, et al. Calcineurin is necessary for the maintenance but not embryonic development of slow muscle fibers. Mol Cell Biol. 2005;25:6629–38.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Maeda T, Chapman DL, Stewart AF. Mammalian vestigial-like 2, a cofactor of TEF-1 and MEF2 transcription factors that promotes skeletal muscle differentiation. J Biol Chem. 2002;277:48889–98. Epub 42002 Oct 48889.

  27. 27.

    Bonnet A, Dai F, Brand-Saberi B, Duprez D. Vestigial-like 2 acts downstream of MyoD activation and is associated with skeletal muscle differentiation in chick myogenesis. Mech Dev. 2010;127:120–36.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Yang Y, Creer A, Jemiolo B, Trappe S. Time course of myogenic and metabolic gene expression in response to acute exercise in human skeletal muscle. 1985) 2005; 98:1745-1752.

  29. 29.

    Ma Z, Sun X, Xu D, Xiong Y, Zuo B. MicroRNA, miR-374b, directly targets Myf6 and negatively regulates C2C12 myoblasts differentiation. Biochem Biophys Res Commun. 2015;467:670–5.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Hawke TJ, Kanatous SB, Martin CM, Goetsch SC, Garry DJ. Rad is temporally regulated within myogenic progenitor cells during skeletal muscle regeneration. Am J Physiol Cell Physiol. 2006;290:C379–87.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Wehling-Henricks M, Jordan MC, Gotoh T, Grody WW, Roos KP, Tidball JG. Arginine metabolism by macrophages promotes cardiac and muscle fibrosis in mdx muscular dystrophy. PLoS One. 2010;5:e10763.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Spang N, Feldmann A, Huesmann H, Bekbulat F, Schmitt V, Hiebel C, Koziollek-Drechsler I, Clement AM, Moosmann B, Jung J, et al. RAB3GAP1 and RAB3GAP2 modulate basal and rapamycin-induced autophagy. Autophagy. 2014;10:2297–309.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Su J, Ekman C, Oskolkov N, Lahti L, Strom K, Brazma A, Groop L, Rung J, Hansson O. A novel atlas of gene expression in human skeletal muscle reveals molecular changes associated with aging. Skelet Muscle. 2015;5:35.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Pandey SN, Cabotage J, Shi R, Dixit M, Sutherland M, Liu J, Muger S, Harper SQ, Nagaraju K, Chen YW. Conditional over-expression of PITX1 causes skeletal muscle dystrophy in mice. Biol Open. 2012;1:629–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Merico D, Roifman M, Braunschweig U, Yuen RK, Alexandrova R, Bates A, Reid B, Nalpathamkalam T, Wang Z, Thiruvahindrapuram B, et al. Compound heterozygous mutations in the noncoding RNU4ATAC cause Roifman Syndrome by disrupting minor intron splicing. Nat Commun. 2015;6:8718.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Nagy R, Wang H, Albrecht B, Wieczorek D, Gillessen-Kaesbach G, Haan E, Meinecke P, de la Chapelle A, Westman JA. Microcephalic osteodysplastic primordial dwarfism type I with biallelic mutations in the RNU4ATAC gene. Clin Genet. 2012;82:140–6.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Tobimatsu K, Noguchi T, Hosooka T, Sakai M, Inagaki K, Matsuki Y, Hiramatsu R, Kasuga M. Overexpression of the transcriptional coregulator Cited2 protects against glucocorticoid-induced atrophy of C2C12 myotubes. Biochem Biophys Res Commun. 2009;378:399–403. 310.1016/j.bbrc.2008.1011.1062.

  38. 38.

    Hountondji C, Bulygin K, Crechet JB, Woisard A, Tuffery P, Nakayama J, Frolova L, Nierhaus KH, Karpova G, Baouz S. The CCA-end of P-tRNA contacts both the human RPL36AL and the A-site bound translation termination factor eRF1 at the peptidyl transferase center of the human 80s ribosome. Open Biochem J. 2014;8:52–67.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Tillmann H, Eschrich K. Isolation and characterization of an allelic cDNA for human muscle fructose-1,6-bisphosphatase. Gene. 1998;212:295–304.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Rakus D, Maciaszczyk E, Wawrzycka D, Ulaszewski S, Eschrich K, Dzugaj A. The origin of the high sensitivity of muscle fructose 1,6-bisphosphatase towards AMP. FEBS Lett. 2005;579:5577–81.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Pirog M, Gizak A, Rakus D. Changes in quaternary structure of muscle fructose-1,6-bisphosphatase regulate affinity of the enzyme to mitochondria. Int J Biochem Cell Biol. 2014;48:55–9.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Herrer I, Rosello-Lleti E, Rivera M, Molina-Navarro MM, Tarazon E, Ortega A, Martinez-Dolz L, Trivino JC, Lago F, Gonzalez-Juanatey JR, et al. RNA-sequencing analysis reveals new alterations in cardiomyocyte cytoskeletal genes in patients with heart failure. Lab Invest. 2014;94:645–53.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Qiao S, Dennis M, Song X, Vadysirisack DD, Salunke D, Nash Z, Yang Z, Liesa M, Yoshioka J, Matsuzawa S, et al. A REDD1/TXNIP pro-oxidant complex regulates ATG4B activity to control stress-induced autophagy and sustain exercise capacity. Nat Commun. 2015;6:7014.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Tsukahara M, Suemori H, Noguchi S, Ji ZS, Tsunoo H. Novel nucleolar protein, midnolin, is expressed in the mesencephalon during mouse development. Gene. 2000;254:45–55.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Hofmeister-Brix A, Kollmann K, Langer S, Schultz J, Lenzen S, Baltrusch S. Identification of the ubiquitin-like domain of midnolin as a new glucokinase interaction partner. J Biol Chem. 2013;288:35824–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Stuart CA, Yin D, Howell ME, Dykes RJ, Laffan JJ, Ferrando AA. Hexose transporter mRNAs for GLUT4, GLUT5, and GLUT12 predominate in human muscle. Am J Physiol Endocrinol Metab. 2006;291:E1067–73.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Karunakaran D, Thrush AB, Nguyen MA, Richards L, Geoffrion M, Singaravelu R, Ramphos E, Shangari P, Ouimet M, Pezacki JP, et al. Macrophage mitochondrial energy status regulates cholesterol efflux and is enhanced by anti-miR33 in atherosclerosis. Circ Res. 2015;117:266–78.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Anunciado-Koza RP, Zhang J, Ukropec J, Bajpeyi S, Koza RA, Rogers RC, Cefalu WT, Mynatt RL, Kozak LP. Inactivation of the mitochondrial carrier SLC25A25 (ATP-Mg2+/Pi transporter) reduces physical endurance and metabolic efficiency in mice. J Biol Chem. 2011;286:11659–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Plaisier CL, Bennett BJ, He A, Guan B, Lusis AJ, Reue K, Vergnes L. Zbtb16 has a role in brown adipocyte bioenergetics. Nutr Diabetes. 2012;2:e46.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Gordon EA, Whisenant TC, Zeller M, Kaake RM, Gordon WM, Krotee P, Patel V, Huang L, Baldi P, Bardwell L. Combining docking site and phosphosite predictions to find new substrates: identification of smoothelin-like-2 (SMTNL2) as a c-Jun N-terminal kinase (JNK) substrate. Cell Signal. 2013;25:2518–29.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by NIH/NIAMS K23AR054904; a pilot project part of NIH/NIA P30AG028716; an Early Career Development Award from the Central Society of Clinical Research; and an ACR-REF/ASP Junior Career Development Award in Geriatric Medicine funded via Atlantic Philanthropies, ACR-REF, John A. Hartford Foundation and ASP. LipoScience, Inc. (now LabCorp) provided GlycA determinations at no cost.

Availability of data and materials

All data are available from the corresponding author upon reasonable request.

Authors’ contributions

KMH, RJ, BA, BND, and MJH contributed to the data analysis and data interpretation, and wrote the manuscript. RN, JLH, JM, BFG, KNT, and MC participated in acquisition of data and laboratory studies, and reviewed/edited the manuscript. VBK, TRK, DMM, and WEK participated in conceptual design and data interpretation, and reviewed/edited the manuscript. All authors read and approved the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was in compliance with the Helsinki Declaration and was approved (Protocol 7701) by the Duke University Institutional Review Board. We obtained written informed consent from all study participants prior to all study activities.

Author information

Affiliations

Authors

Corresponding author

Correspondence to Kim M. Huffman.

Additional file

Additional file 1:

Supplemental gene information. Additional detail on genes that were most differentially expressed in muscle from patients with rheumatoid arthritis and from control participants. Genes are described and categorized by the proposed function of their respective gene products (DOCX 88 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Huffman, K.M., Jessee, R., Andonian, B. et al. Molecular alterations in skeletal muscle in rheumatoid arthritis are related to disease activity, physical inactivity, and disability. Arthritis Res Ther 19, 12 (2017). https://doi.org/10.1186/s13075-016-1215-7

Download citation

Keywords

  • Gene expression
  • Metabolomics
  • Satellite cells
  • Fibrosis
  • Inflammation