Effects of delayed NSAID administration after experimental eccentric contraction injury – A cellular and proteomics study

Abstract

Background

Acute muscle injuries are exceedingly common and non-steroidal anti-inflammatory drugs (NSAIDs) are widely consumed to reduce the associated inflammation, swelling and pain that peak 1–2 days post-injury. While prophylactic use or early administration of NSAIDs has been shown to delay muscle regeneration and contribute to loss of muscle strength after healing, little is known about the effects of delayed NSAID use. Further, NSAID use following non-penetrating injury has been associated with increased risk and severity of infection, including that due to group A streptococcus, though the mechanisms remain to be elucidated. The present study investigated the effects of delayed NSAID administration on muscle repair and sought mechanisms supporting an injury/NSAID/infection axis.

Methods

A murine model of eccentric contraction (EC)-induced injury of the tibialis anterior muscle was used to profile the cellular and molecular changes induced by ketorolac tromethamine administered 47 hr post injury.

Results

NSAID administration inhibited several important muscle regeneration processes and down-regulated multiple cytoprotective proteins known to inhibit the intrinsic pathway of programmed cell death. These activities were associated with increased caspase activity in injured muscles but were independent of any NSAID effect on macrophage influx or phenotype switching.

Conclusions

These findings provide new molecular evidence supporting the notion that NSAIDs have a direct negative influence on muscle repair after acute strain injury in mice and thus add to renewed concern about the safety and benefits of NSAIDS in both children and adults, in those with progressive loss of muscle mass such as the elderly or patients with cancer or AIDS, and those at risk of secondary infection after trauma or surgery.

About the Authors:

Amy E. Bryant

* E-mail: [email protected], [email protected]

Affiliations U.S. Department of Veterans Affairs, Office of Research and Development, Boise, ID, United States of America, University of Washington School of Medicine, Seattle, WA, United States of America

ORCID http://orcid.org/0000-0002-4143-1718

Michael J. Aldape

Affiliations U.S. Department of Veterans Affairs, Office of Research and Development, Boise, ID, United States of America, Northwest Nazarene University, Nampa, ID, United States of America

Clifford R. Bayer

Affiliation: U.S. Department of Veterans Affairs, Office of Research and Development, Boise, ID, United States of America

Eva J. Katahira

Affiliation: U.S. Department of Veterans Affairs, Office of Research and Development, Boise, ID, United States of America

Laura Bond

Affiliation: Boise State University, Boise, ID, United States of America

Carrie D. Nicora

Affiliation: Pacific Northwest National Laboratory, Richland, WA, United States of America

Thomas L. Fillmore

Affiliation: Pacific Northwest National Laboratory, Richland, WA, United States of America

Therese R. W. Clauss

Affiliation: Pacific Northwest National Laboratory, Richland, WA, United States of America

Thomas O. Metz

Affiliation: Pacific Northwest National Laboratory, Richland, WA, United States of America

Bobbie-Jo Webb-Robertson

Affiliation: Pacific Northwest National Laboratory, Richland, WA, United States of America

Dennis L. Stevens

Affiliations U.S. Department of Veterans Affairs, Office of Research and Development, Boise, ID, United States of America, University of Washington School of Medicine, Seattle, WA, United States of America

Introduction

Non-steroidal anti-inflammatory drugs (NSAIDs) are cornerstones of pain management in homes and hospitals worldwide. These agents are consumed or prescribed for diverse conditions ranging from mild intermittent to chronic musculoskeletal pain, as well as for pain associated with cancer, AIDS, surgery and rehabilitation after placement of prosthetic devices. While NSAIDs have their place in these settings, they are not innocuous and can contribute to peptic ulcers, upper gastro-intestinal bleeding, renal disease and adverse cardiovascular events [1]. NSAID use has also been linked to poor wound healing and to increased risk of complications following surgery, including anastomotic failure following colorectal resection (reviewed in [2]).

NSAID use is widely advocated to manage the inflammation, pain and swelling associated with acute muscle strain injuries. Such injuries are exceedingly common among athletes, military personnel, and civilians of all ages. Within the first 24 to 48 hours after unaccustomed exercise or injury, a constellation of symptoms—referred to as delayed onset muscle soreness (DOMS)—develops that includes pain with muscle use, stiffness and swelling. Without treatment, these symptoms peak between 24 and 72 hours, and typically resolve within 5 to 7 days [3]. Yet, recent evidence suggests that the inflammatory response to injury is a necessary phase of soft tissue healing and its inhibition with NSAIDs can significantly delay muscle regeneration and decrease muscle strength after repair (reviewed in [4]). The deleterious NSAID effects on muscle repair have been attributed to their ability to block of cyclooxygenase-2 (COX-2)-derived prostaglandins which are known to stimulate muscle progenitor cell (satellite cell) responses to exercise [5–7]. These studies compared satellite cell activity in human volunteers or animal subjects when NSAIDs were administered prior to resistance training or to exercise-induced injury. In human studies, administration of non-selective NSAIDs (e.g., ibuprofen) blocked the exercise-induced increase in satellite cell number [7,8] and mixed muscle protein synthesis [9]. Similarly, animal studies demonstrate a negative effect of NSAIDs on satellite cell responses to exercise and on regeneration after injury [5,10]. Although NSAID use prior to exercise/injury is not a typical practice in the general population, these results clearly indicate that prophylactic use of NSAIDs, as is common among professional and elite athletes [4], exerts a negative effect on satellite cell activity.

In addition to their effects on muscle regeneration, recent clinical and epidemiological studies have suggested that NSAIDs may also predispose to severe bacterial infections including those due to Clostridium difficile [11,12], Streptococcus pneumoniae [13,14], Staphylococcus aureus [13,15], Gram negative organisms [15] and group A streptococcus (reviewed in [16]). Regarding the latter infection, a clinical entity termed “cryptogenic” or “spontaneous” group A streptococcal myonecrosis has been described in which patients initially present to their health care provider for increasing localized pain, typically at a site of non-penetrating soft tissue injury such as a muscle strain or bruise. Without an obvious portal of bacterial entry and in the absence of the classical cutaneous signs of developing infection, the correct diagnosis is missed or delayed. Instead, patients are prescribed NSAIDs and sent home, only to return 24–48 hrs later in multiple organ failure and intractable hypotension. Nearly 50% of streptococcal toxic shock syndrome (StrepTSS) patients with necrotizing soft tissue infection fall into this “cryptogenic” group [17]. Mortality is high (30–85%) and morbidity is extensive. By the time the systemic manifestations of StrepTSS are evident and the correct diagnosis of myonecrosis is made, emergent surgery—often including life-altering multiple limb amputations—is frequently required to ensure survival. Survivors are often left with profound cognitive impairments and endure prolonged hospitalization and rehabilitation at great emotional and financial expense.

Though a large body of clinical evidence supported a possible injury/NSAID/group A streptococcal axis, controversy remained. Opponents have argued that NSAID use merely masks the signs and symptoms of developing infection, thereby delaying diagnosis and treatment. Since no experimental studies had directly tested this possible association, we developed a murine model of exercise-induced muscle strain coupled with group A streptococcal bacteremia in which we investigated this notion [18]. We showed that NSAID administration 47 hrs post-injury significantly increased the trafficking of circulating group A streptococcus to the injured site [18]. In addition, using a murine model of established infection, we further demonstrated that three different non-selective NSAIDs each significantly increased disease severity and reduced the time to 100% mortality [19]. In total, this evidence suggested to us that NSAIDs directly contribute to poor outcomes in these infections.

As mentioned above, most studies examining the mechanisms by which NSAIDs affect muscle regeneration used models in which NSAIDs were given prior to resistance training or exercise-induced injury. Few studies investigated NSAID effects when given later–at times analogous to the peak of perceived muscle soreness in humans. Thus the present study utilized flow cytometric analysis of infiltrating inflammatory cells and an unbiased proteomics approach to investigate the cellular and protein changes in murine tibialis anterior (TA) muscle after eccentric contraction (EC)-induced muscle injury in the presence and absence of ketorolac tromethamine–a non-selective NSAID—given 47 hrs after muscle injury.

Findings document a predominantly inhibitory effect of ketorolac on muscle metabolism and biogenesis that was independent of any significant change in the numbers of infiltrating inflammatory cells or their functional phenotypes. Instead, results suggest that NSAID use at this time significantly down-regulates several anti-apoptosis proteins that inhibit the extrinsic pathway of programmed cell death. Such down-regulation of pro-survival proteins was accompanied by increased caspase activity in the muscle. By promoting a pro-apoptosis phenotype in regenerating muscles, NSAIDs could expand the nidus of injury, and thereby increase muscle susceptibility to post-injury infection. This activity could also contribute to the known ability of NSAIDs to delay muscle regeneration and reduce muscle strength post-healing. These findings warrant further studies in humans and if confirmed, this new knowledge may shift the current paradigm of pain management in multiple clinical settings.

Methods

Ethical approval

Studies involving animal subjects were approved by the Institutional Animal Care and Use Committee (IACUC) U.S. Department of Veterans Affairs Medical Center, Boise, ID and adhered to guidelines specified by the Veterans Health Administration (handbook 1200.07), the U.S. Public Health Service and the National Institutes of Health. Work was conducted in the Boise VA Medical Center’s Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC)-accredited Veterinary Medical Unit. Animals were housed in groups of 5–10 in standard large rodent cages and allowed food and drink ad libitum. All investigators understand this journal’s ethical principles regarding animal subjects in research and confirm that our work complies with these principles.

Animal model

This model has been previously described in detail [18]. In brief, adult Swiss Webster outbred mice (female, 20–25 g, Taconic Laboratories; N = 9/group) were anesthetized with 2% inhaled isoflurane in oxygen and positioned in the temperature-controlled exercise apparatus (modified model 360B, Aurora Scientific, Ontario, Canada) that permits control of both ankle dorsiflexion torque and muscle length. Anesthesia was maintained throughout the experimental procedure via delivery of 2% isoflurane via nose cone and breathing rates were carefully monitored. Adequacy of anesthesia was defined as no withdrawal response to toe-pinch.

The tibialis anterior (TA) muscle was stimulated electrically (10–12 volts, 100–150 Hz for 400 ms) by sterile 28-gauge needle electrodes placed subcutaneously ~2 mm apart near the right peroneal nerve just lateral to the midline and distal to the knee joint. Beginning at 200 ms after nerve stimulation and at maximal muscle contraction, the foot was rotated distally 78 degrees. At the end of the stimulation period, the foot was returned to the starting position. This exercise regimen was repeated every 30 seconds 50 times. Prior to, and immediately after exercise, the maximal mean isometric torque (MIT) generated by two isometric contractions (i.e., the foot held motionless) was recorded in each animal. Only animals displaying a post-EC reduction in MIT of >40% were included for study since we [18] and others [20] have previously shown that this level of functional impairment constitutes a moderate muscle injury. The contralateral leg of each animal was not exercised and served as control. After exercise, a single dose of buprenorphine (0.1 mg/kg in 0.5 mL sterile saline) was administered to minimize post-exercise distress and to replenish fluids lost during the 1 hr exercise regimen. This opioid analgesic agent is active at the central nervous system level; the single dose used is 20-fold less than that shown to have modest transient anti-inflammatory activity in the rat (2 mg/kg q 12 hr) [21].

At 47 hrs after EC, animals were treated with ketorolac tromethamine (trade name, Toradol; 7.5 mg/kg by IP injection) or its vehicle (saline). Ketorolac is a prescription non-selective NSAID commonly used in clinics and emergency room settings for short-term management of moderately severe, acute muscle pain in adults. It was chosen for these studies because 1) it is of the same chemical subclass as ibuprofen (propionic acid) with similar COX1/COX2 inhibitory-50 [IC₅₀] ratios (.35 and .5, respectively) [22] and can be given parenterally to experimental animals; 2) in our previously published experimental studies in mice, ketorolac significantly increased trafficking of circulating group A streptococcus to sites of EC-induced muscle injury [18] and 3) in established group A streptococcal myonecrosis in mice, ketorolac was equivalent to ibuprofen and indomethacin in the ability to significantly increase severity of infection and decrease the time to 100% mortality [19]. The time of NSAID administration was chosen for several reasons: 1) our data (below) and those from Lieber’s group [23,24] have shown that muscle force-generating capacity remains significantly reduced at 48 hrs post injury and that this dysfunction is associated with marked architectural and inflammatory changes and in myoregulatory gene expression [23,24]; 2) this time post-injury is related to the peak soreness in humans after unaccustomed exercise or muscle injury and when NSAID use is prevalent and 3) studies by us (presented here) and by Arnold et al [25] have shown that by this time after experimental strain injury, polymorphonuclear leukocyte (PMNL) influx has waned and infiltrated pro-inflammatory macrophages have begun conversion to an anti-inflammatory, pro-resolution phenotype; thus any NSAID effect that might be observed would not be attributable to its ability to limit the tissue leukocytic inflammatory response.

Seven hours after NSAID or vehicle administration (54 hr post-injury), mice were anesthetized as described above, heparinized blood specimens obtained by retro-orbital puncture and, while still under anesthesia, the animals were sacrificed by cervical dislocation. The TA muscles from the left and right legs were harvested and flash frozen for proteome analysis. The timing of muscle harvest was based on our previously published work demonstrating a significant NSAID-induced increase in trafficking of group A streptococcus to the TA muscle at this time post-injury [18]. Plasma samples were frozen at -70C. In separate studies, animals were exercised as described above and treated with NSAID or saline at 47 hrs post injury. At selected times thereafter as indicated in the figure legends, TA muscles were harvested post-mortem and were either flash frozen for quantitative RT-PCR analysis for genes of interest, fixed in neutral buffered formalin for routine histopathology or processed for leukocyte isolation or caspase activity assay as described below.

Proteomic analyses

The accurate mass and time (AMT) tag approach [26] was used to generate quantitative proteomics data for muscle tissues from the above described animals (see text in S1 Supporting Information: Additional Methodology for Proteomics Analysis for full experimental details). Comprehensive shotgun proteomics analyses were first performed to generate an AMT tag database of identified mouse peptides/proteins. This database was then used to identify injury- and NSAID-associated quantitative and qualitative differences in protein expression. Specifically, injured and non-injured TA muscles from each animal treated with or without NSAIDs at 47 h after injury (9 mice/group) as described above were harvested at 7 h after NSAID (or vehicle) administration, flash frozen and sent on dry ice to Pacific Northwest National Laboratory. Tissues from injured and non-injured TA muscles from all treatment groups were individually homogenized, digested using trypsin, and aliquots from each digest were pooled and subjected to offline fractionation [27]. Each fraction was then analyzed by nano capillary LC-MS/MS, and the identified peptides were used to populate the mouse muscle peptide/protein AMT tag database. Peptides were identified with a false discovery rate <0.1% and with a mass measurement error of +/- 4ppm.

The resultant comprehensive protein database was then used as a ‘look up table’ for the subsequent treatment group-specific quantitative proteomics analyses, which involved analysis of the individual tissue protein digests using the same LC-MS instrument platform, followed by data analysis using the AMT tag approach. Briefly, the observed masses and normalized elution times (NETs) of peptides detected in the individual samples were matched against the entries in the AMT tag database using the PRISM Data Analysis system [28], which is a series of software tools developed in-house (e.g. Decon2LS [29] and VIPER [30]) and freely available at omics.pnl.gov/software. Individual steps in this data processing approach are reviewed in [26]. The peptides identified from this matching process and the associated integrated LC-MS peak areas were retained as a matrix for subsequent data analysis.

LC-MS log transformed peak intensity (i.e., abundance) data was processed through appropriate quality control (QC) filters. One mouse sample did not meet adequate quality metrics associated with correlation and missing values and was removed from subsequent analyses [31]. Following QC, the dataset was appropriately normalized to reduce systematic effects of the analytical process, such as variation in total sample protein. The log transformed normalized peptide abundance dataset was appropriately filtered to remove peptides with inadequate data for statistical analysis. All proteins having 1 or more peptides passing the above-mentioned filters were assessed.

Each protein that met both the criteria of sample-level “Coverage” [Number of paired injured and control muscles ≥ 50% or (N_(injured) or N_(control)) ≥75%] and of sample-level “Difference” [|Fold Change| ≥1.25 or |N_(injured)—N_(control)| ≥ 50%) were subjected to statistical analysis to identify those proteins that were significantly affected by injury or by NSAID. Proteins were analyzed for quantitative effects associated with injury or with NSAID administration using a standard T-test and for qualitative effects using G-test methodology [32]. The T-test indicated proteins that are present in both treatment groups that have a quantitative change in expression, i.e., up- or down-expression, and G-test indicated proteins that are present in one group and absent from another for identification of qualitative changes. Protein quantitative measures were computed using BP-Quant [33] and standard averaging methods [34] and statistical analysis followed the same T-test and G-test strategies.

Functional analysis

Diseases and functions associated with injury +/- NSAID administration were identified in part using QIAGEN's Ingenuity Pathway Analysis (IPA, version 9.0, QIAGEN Redwood City, www.qiagen.com). Briefly, a dataset of significantly altered proteins containing the protein identifier, the P value, and the fold-change was uploaded into IPA and a core analysis was performed. For proteins that were significant by G-test (and for which a fold-change is not relevant), values of 2 or -2 were assigned to indicate those that were qualitatively increased or decreased with treatment, respectively, compared to non-treated controls.

In addition, initial protein functional categorization was performed by GeneOntology’s PANTHER classification system (www.pantherdb.org). Primary function was further refined based on descriptions provided by the NCBI Protein database and the available literature. Here, the functional categories were based on the updated Clusters of Orthologous Groups (COG) of proteins described by Tatusov et al [35] with some modifications. For multi-functional proteins, the primary category was assigned by considering a role in skeletal muscle structure, function and biogenesis.

Gene expression analyses

The TA muscles from both exercised and the contra-lateral, non-exercised legs of saline- or NSAID-treated animals (4-12/group) were harvested at 4 or 7 hrs post-treatment (51 and 54 hrs post-injury) and immediately flash-frozen in liquid nitrogen for subsequent analysis of gene expression by real time quantitative reverse transcriptase PCR (qRT-PCR). Frozen muscles, suspended in TriReagent and 4-bromoanisole (Molecular Research Center Inc, Cincinnati, OH), were disrupted in a TissueLyser II (Qiagen, Valencia, CA). Total RNA was extracted from the resulting homogenate using an RNeasy RNA isolation kit (Qiagen), then treated with TURBO DNA-free (Ambion, Grand Island, NY) to remove contaminating DNA. RNA quantity was determined by Qubit (Invitrogen, Grand Island, NY) analysis and quality was assessed by TapeStation (Agilent, Santa Clara, CA) analysis. cDNA was synthesized from 2 micrograms of isolated total RNA using MMLV reverse transcriptase (New England Biolabs, Ipswich, MA) then diluted 1:10. Following manufacturer’s protocol, 1ul of the diluted cDNA was added to 96 well reaction plates (MicroAmp™ Fast Optical plates; Applied Biosystems) containing 1X SYBR Green/Rox PCR master mix (Qiagen) and 0.4uM RT² qPCR Primer Assay mix for murine genes of interest (Qiagen, Valencia, CA). qRT-PCR was performed on either an ABI 7500 Fast PCR (Applied Biosystems) or Realplex 2 (Eppendorf) machine using the following parameters: 95°C for 10 min; 40 cycles of 95°C for 15 sec and 60°C for 1 min. A dissociation (melting) curve was run immediately after the PCR program to verify product integrity. The threshold cycle (Ct) for each transcript and change in threshold cycle (ΔCt) for the gene of interest/GAPDH pair of transcripts were monitored for each amplification reaction and the data was analyzed using the 2^-ΔΔCt method [36] to calculate fold change between the exercised and non-exercised muscles. Statistical differences were determined by inspection of 95% confidence intervals (CI) of the log₂-transformed fold-change values (ΔΔC_t; non-coverage of 0 by the CI indicates statistical difference between treatment and control at p < 0.05).

Leukocyte isolation

To assess the dynamics and phenotypic characteristics of infiltrating cellular response after injury, TA muscles from 3 animals/group were removed and digested following the methods of Arnold et al [25] with some modifications. Muscles were rinsed in DPBS then digested twice in DMEM containing 0.2% collagenase B and 0.2% trypsin at 37°C for 1 hr. Resultant material was filtered with a 70 μm cell strainer. Isolated cells were recovered by centrifugation and suspended in PBS pH 7.2, 0.5% BSA, and 2mM EDTA and passed through a 30 μm MACS pre-separation filter (Miltenyi Biotec, Auburn, CA). After treatment with an FcR blocking reagent (Miltenyi Biotec) to prevent nonspecific binding, CD11b+ cells were isolated by positive selection using magnetically labeled microbeads conjugated to anti-murine CD11b antibody (rat IgG2b; clone M1/70.15.11.5; Miltenyi Biotec) and selected by magnetic column according to the manufacturer’s instructions. Murine monocytes/macrophages (MO/MP) express abundant, high affinity surface CD11b. In contrast, resting granulocytes (including eosinophils) and NK cells express fewer, lower affinity CD11b molecules; upon activation, a functional up-regulation (i.e., increased affinity) as well as increased CD11b expression is observed. Total cell number was assessed by direct hemocytometer counting and viability was >98% by trypan blue exclusion.

A portion of the final cell suspension was fixed to a microscope slide by cytospin, stained with Wright’s stain, and a manual differential performed by a blinded animal pathologist at the Caine Veterinary Teaching Hospital, Caldwell, Idaho. To determine the phenotype of the isolated MO/MP, cells were incubated with either PE-Ly-6C (Miltenyi Biotec), FITC-F4/80 (eBioscience), or isotype control antibodies (eBioscience) at 4°C for 10 min. Cells were washed once and resuspended in 0.5 mLs PBS with 0.5% BSA and 2 mM EDTA. Analysis was immediately performed on an Epics XL flow cytometer with the monocyte gate drawn based on size and granularity. An average of 10,000 events was analyzed per sample. MO/MPs in the contra-lateral, non-exercised TA muscle were similarly isolated in parallel, however too few cells were recovered for subsequent analysis. Statistical differences were assessed using the Mann-Whitney test, based on an exact distribution [37,38], as an alternative to the parametric t test.

Caspase activity assay

TA muscles were harvested from injured and non-injured legs 7 hrs after NSAID or vehicle treatment (i.e., 54 hrs post-injury). Each TA muscle (4/group) was flash frozen and pulverized by Qiagen Tissuelyzer in 250 uL of extraction buffer (25 mM HEPES pH 7.5, 0.1% Triton X-100, 5 mM MgCl₂, 2 mM DTT, 74 uM antipain, 0.15 uM aprotinin, 1.3 mM EDTA, 20 mM leupeptin and 15 uM pepstatin). Digests were clarified by centrifugation (15,000xg, 20 min) and protein concentration determined by bicinchoninic acid assay (BCA) assay (Pierce). Combined activities of caspases 3 and 7 in the digests were measured by commercial ELISA (Promega) which utilizes a profluorescent DEVD peptide-rhodamine 110 substrate [(Z-DEVD)2-R110]. Data are reported as Relative Luminescence Units (RLU)/mg protein normalized to that of the non-injured control. Statistically significant differences were determined by ANOVA using a mixed model with random effect of animal within treatment. The level of significance was set at P < .05.

Results

Eccentric contraction model of muscle injury

To investigate the effects of injury with and without NSAID treatment on the regenerating skeletal muscle proteome and cellular inflammatory response, we utilized our established murine model of eccentric contraction (EC)-induced muscle injury [18]. This model was chosen since 1) EC injury is more relevant to the human condition than is freeze- or toxin-induced injury, 2) the TA muscle largely consists of Type II (fast twitch) fibers [39] thereby obviating the confounding influence of multiple fiber types in analysis; 3) fast-twitch muscles are considered more susceptible to EC-induced injury [40] and 4) because NSAID consumption for pain and swelling is highly associated with such injury [40].

Repeated bouts of EC exercise caused an ~50% loss in TA functional capacity (measured as a reduction in Mean Isometric Torque, MIT) which is indicative of moderate muscle injury [41]. This reduction in MIT was not different among animals subsequently randomized to “No NSAID” and “NSAID” groups (No NSAID: 57.7% ± 5.7% vs NSAID: 59.8% ± 5.0%). EC injury also resulted in a marked tissue inflammatory response at 24 hrs (Fig 1) and stimulated myogenic gene expression (Fig 2). Specifically, significant early and sustained increases were observed in myod1 (a transcription factor controlling satellite cell activation and proliferation [20]) and vim (a structural protein upregulated in proliferating myoblasts [42]), whereas ptgs2 (prostaglandin synthetase 2 [aka cyclooxygenase 2, COX-2], an early mediator of muscle repair [5]), was highest at 24 hr post-injury and declined thereafter (Fig 2). As noted by others [5], expression of ptgs1 (prostaglandin synthetase 1, [COX-1]) was only minimally increased with injury. By 48–54 hrs, myog (an inducer of myoblast maturation) was maximally upregulated suggesting that many proliferating myoblasts had become committed to the differentiation program.

[Figure omitted. See PDF.]

Fig 1. Experimental muscle injury disrupts myofiber architecture and causes marked influx of inflammatory cells.

Mice underwent the eccentric contraction (EC) exercise regimen as described in Methods. At 24 hr post-EC, animals were euthanized and both left and right tibialis anterior (TA) muscles were harvested post-mortem for routine histopathology. One representative animal of 3 is depicted. The exercised muscle, but not the control TA, displays destruction of normal muscle architecture and a marked inflammatory cell influx. These pathologies, and the corresponding reduction in muscle force, are accepted criteria that define muscle injury.

https://doi.org/10.1371/journal.pone.0172486.g001

[Figure omitted. See PDF.]

Fig 2. Experimental eccentric contraction injury stimulates physiological responses characteristic of muscle regeneration.

Mice (4-12/group) underwent the EC regimen described above. Non-injured and injured TA muscles were harvested at the indicated times post-injury, flash-frozen and used to evaluate the relative expression of genes involved in muscle regeneration by qRT-PCR as described in Methods. Data are given as the mean fold-change in gene expression relative to the non-injured muscle ± standard error of the mean (SEM). Calculation of the 95% confidence interval using the log₂-transformed fold-change values revealed that all genes were significantly increased at each time point relative to the non-injured control as indicated by the grouped asterisks, except the 24h ptgs1 expression. No corresponding proteomics analyses were conducted at 24 or 48 hrs; at 54 hrs post-injury the only corresponding protein that was significantly altered by injury was vimentin (Table in S1 Table: Proteins Significantly Altered by Injury; 1.745-fold increase; p = .0236).

https://doi.org/10.1371/journal.pone.0172486.g002

Thus, our model fulfills multiple criteria that define acute muscle injury–namely a loss in functional capacity, changes in muscle architecture, an early influx of inflammatory cells [20] and expression of key genes involved in muscle regeneration after injury [5,23,43]. This model is therefore appropriate for the study of the effects of delayed NSAID administration on muscle regeneration after acute strain injury.

Injury, NSAIDs and the cellular inflammatory response

Muscle regeneration depends on a highly synchronized influx of different leukocyte populations with shifting functional capacities (reviewed [44]). Specifically, after the initial neutrophil response has waned, infiltrated pro-inflammatory monocytes must switch into pro-resolution macrophages to support myogenesis [25]. To examine these dynamics in our model, total CD11b⁺ leukocytes were isolated over time from TA muscles by positive magnetic sorting. Too few cells were obtained from non-injured TA muscles to permit further analysis (not shown). In EC-injured TA muscles from untreated mice, the number of CD11b⁺ cells isolated at 24 and 48 hrs was not statistically different (Table 1); however the subpopulations were decidedly distinct. At 24 hr, neutrophils and pro-inflammatory M1 monocyte/macrophages (MO/MPs; Ly-6C⁺) predominated; few pro-resolution M2 MO/MPs (F4/80⁺) were detected. By 48 hrs, neutrophils were significantly reduced whereas the MO/MP population had increased and their phenotype was shifting from M1 to M2 as indicated by a >50% reduction in Ly-6C⁺ cells and an increase in F4/80 expression (Table 1).

[Figure omitted. See PDF.]

Table 1. Dynamics and Cellular Phenotypes of Infiltrating Leukocytes after Muscle Injury.

https://doi.org/10.1371/journal.pone.0172486.t001

To test whether delayed NSAID administration could alter the leukocytic inflammatory response in tissue, animals were treated 47 hrs after injury with either saline or NSAID; TA muscles were harvested 7 hrs later. At 54 hrs post-injury, vehicle-treated control animals showed a further reduction in the percent of M1 macrophages and a continued increase in M2s (Table 1) demonstrating that the initial pro-inflammatory response (PMNL, M1 macrophages) had peaked and the pro-resolution cellular response was well underway. Ketorolac administration at 47 hrs post-injury did not significantly alter the number of CD11b⁺ cells, the leukocyte differential, or the M1/M2 percentages in the tissues (Table 1).

Injury, NSAIDs and the skeletal muscle proteome

An unbiased proteomics approach was used to investigate the global effects of injury, with and without NSAID administration, on regenerating muscles. In non-NSAID treated mice, the AMT tag proteomic approach described in Methods identified 89 significant, differentially-expressed proteins in injured muscles at 54 hrs post-EC (Table in S1 Table: Proteins Significantly Altered by Injury); of these, 67 (75%) were increased compared to non-injured TA muscles and were predominantly associated with muscle repair processes. For example, calpain (CAN1), a calcium-activated protease associated with the degenerative phase of muscle repair [45], was significantly increased. Proteins essential for, or indicative of, early muscle regeneration were also increased including vimentin (VIME), muscle-type acylphosphatase-2 (ACYP2), and protein unc-45 homolog B (UN45B). Myoblast proliferation was also suggested by significant increases in proteins driving cell division [i.e., the 55 kDa regulatory subunit of serine/threonine protein phosphatase 2A (2ABD) and mitochondrial fission 1 protein (FIS1)]. Similarly, active repair of the neuromuscular junction and the sacrolemmal membrane appeared to be underway as evidenced by increased dihydropyrimidinase-related proteins 1 and 2 (DPYL1, 2) and dysferlin (DYSF), respectively. Revascularization was also promoted via a down regulation of collagen alpha-1(XV) chain (COFA1) whose C-terminal region contains the anti-angiogenic factor, endostatin. Endogenous defenses against stress-induced cytotoxicity appeared activated as indicated by increases in several chaperone proteins [clusterin (CLUS), DnaJ homolog subfamily A member 2 (DNJA2), heat-shock protein beta-2 (HSPB2)] and an inhibitor of ROS-induced apoptosis [glutaredoxin-related protein 5 (GLRX5) [46]]. Serum-associated proteins such as albumin (ALBU) and serotransferrin (TRFE) were also significantly increased indicating that edema remains a prominent feature at this point in the injury resolution process.

In contrast, NSAID administration significantly altered 277 proteins in the injured TA muscles (Fig 3A; also see Table in S2 Table: Proteins Significantly Altered by NSAID) compared to injured muscles from vehicle-treated animals. Of these, 71% (196/277) were decreased by NSAIDs including 14 metabolism and pro-biogenesis proteins that were found to be significantly increased by injury alone (Table 2 and S2 Table: Proteins Significantly Altered by NSAID). These findings support a direct NSAID-induced reversal of key repair responses to injury.

[Figure omitted. See PDF.]

Fig 3. AMT-tag proteomics identifies proteins from injured muscles that are significantly altered by NSAID treatment.

(A) Proteins that met both the criteria of “Coverage” (defined as: [Number of paired injured and control muscles ≥ 50% or (N_(injured) or N_(control)) ≥75%]) and of “Difference” (defined as: [|Fold Change| ≥1.25 or |N_(injured)—N_(control)| ≥ 50%)]) were subjected to statistical analysis to identify those that were significantly affected by NSAID treatment. Proteins were analyzed for quantitative effects using a standard T-test and for qualitative effects using G-test methodology [32]. In total, 277 differentially expressed proteins were found to be significantly altered by NSAID administration. (B) Functional categories of the 277 differentially expressed proteins (assigned as described in Methods; Table in S2 Table: Proteins Significantly Altered by NSAID) are depicted. PTM: post-translational modification; ECM: extracellular matrix.

https://doi.org/10.1371/journal.pone.0172486.g003

[Figure omitted. See PDF.]

Table 2. Common proteins significantly altered by injury and by NSAID administration.

https://doi.org/10.1371/journal.pone.0172486.t002

IPA core analysis of the 277 NSAID-regulated proteins returned “Cell Death and Survival” as the top molecular and cellular functional category (128 unique proteins; P < .001; Fig 4). The activation state of all subcategories involving cell death was predicted to be increased, whereas those involving cell survival and proliferation were predicted to be decreased (Fig 4). Interestingly, 19 proteins mapped to subcategories of “Necrosis of Muscle” and/or to “Cell Death of Muscle”. These subcategories had z-scores of 2.729 and 3.085, respectively, predicting an increased activation state (Table 3). In contrast, the subcategories of “Cell Survival” (43 proteins; z-score = -3.339) and “Cell Proliferation” (103 proteins; z-score = -2.721) were predicted to be decreased by NSAID administration (Table 3).

[Figure omitted. See PDF.]

Fig 4. IPA reveals “Cell Death and Survival” as the major biofunction affected by NSAID administration after muscle injury.

Heat maps of biological functions in which significantly altered proteins were associated with NSAID administration. Color scale is based on z-score with orange indicating predicted activation and blue indicating inhibition. Top: Global view of all returned biofunctions. Bottom: Enlarged view of “Cell Death and Survival” showing “Cell Death of Muscle Cells” (white star), “Necrosis of Muscle” (yellow star) and “Cell Survival” (blue star).

https://doi.org/10.1371/journal.pone.0172486.g004

[Figure omitted. See PDF.]

Table 3. IPA predicts an NSAID-induced Increase in Cell Death and a Reduction in Proliferation/Survival in Regenerating TA Muscles.

https://doi.org/10.1371/journal.pone.0172486.t003

Closer inspection by Gene Ontology analysis and available information in the NCBI protein database showed that the proteins affected by NSAID administration are broadly distributed and enriched in the nucleus, mitochondria, neuromuscular junction and plasma membrane. These proteins are involved in several key biological functions (Fig 3B; also see Table in S2 Table: Proteins Significantly Altered by NSAID) including muscle metabolism, biogenesis and function, and cellular processes including apoptosis as described below.

Metabolism and transport

This functional category contained the greatest number of NSAID-affected proteins (68/277; 24.5%). While not muscle-specific, these proteins are heavily involved in muscle energetics and 88% were down-regulated by NSAIDs including those involved in glycolysis, TCA cycle, electron transport and beta-oxidation. This finding is consistent with the known ability of NSAIDs to delay muscle regeneration after injury [6,8,10] and to suppress muscle metabolism after strength-training exercise (reviewed in [47]).

Muscle specific biogenesis/structure/function

Myogenesis: NSAID administration negatively impacted several proteins associated with post-injury myogenesis. Loss of desmin (DESM) is a hall-mark feature of muscle injury [48] and NSAID treatment significantly reduced this protein. Vimentin (VIME), an intermediate filament protein expressed in regenerating but not mature muscle fibers [49,50] was also down-regulated. Although vimentin is expressed by other cell types including macrophages[51], results from our cellular studies above suggest that the NSAID-induced reduction in vimentin was not due to reduced inflammatory cell influx. Proteins that regulate myoblast differentiation/maturation were also negatively affected by NSAID administration including ankrin repeat domain-containing protein-2 (ANKR2, a coordinator of proliferation and apoptosis during myoblast differentiation; also abbreviated ANKRD2) [52], cysteine and glycine-rich protein 3 (CSRP3, also known as muscle LIM protein [MLP]–a critical positive regulator of myogenic differentiation [53]) and nascent polypeptide-associated complex subunit alpha; muscle-specific form (NACAM—a muscle-specific transcription factor present in differentiated myotubes but not myoblasts). In contrast, injury alone increased ANKR2 and CSRP3 by 2.5 and 2.1-fold, respectively. Although these latter increases did not reach statistical significance (p = .07), they are consistent with studies by Barash et al whose genomic profiling of EC-injured mouse TA muscles showed significant increases in ANKRD2 and CSRP3/MLP gene expression that peaked at ~12 hrs post-injury and remained elevated at 48 hrs [23].

Sarcolemmal Repair: EC-induced injury is associated with disruptions in sarcolemmal integrity and repair of such lesions depends on the coordinated interaction of tripartite motif-containing protein 72 (TRI72), dysferlin (DYSF) [54,55] and microtubules [56]. Our data demonstrate that NSAID administration reduced TRI72, DYSF and multiple members of the beta tubulin family (TBB2A, 2B, 4A, 4B, 6). To our knowledge, these data provide the first evidence that NSAIDs may interfere with repair of sarcolemmal integrity after injury.

Cytoskeleton: NSAID administration also significantly impacted multiple myocytoskeletal proteins. Specifically, 6 isoforms of adult skeletal muscle myosin heavy chain (MYH), 1 perinatal myosin heavy chain and 1 myosin regulatory light chain were significantly increased by NSAID administration. Three of these (MYH4, MYH6, and MRLS) are characteristic of fast-twitch muscles, such as the TA muscle used here, whereas MYH1, 7 and 7B are largely associated with slow-twitch fibers. Actin bundling/capping proteins were also increased including F-actin-capping protein subunits alpha-1 and -2 (CAZA1, CAZA2), cofilin-1 (COF1), and two forms of transgelin (TAGL, TAGL2). NSAID administration also significantly increased 3 members of the actinin family (ACTN1, 2, 3), myozenin-3 (MYOZ3) and telethonin (TELT). ACTN2 and -3 are skeletal muscle-specific and, together with MYOZ3 and TELT, link Z-line proteins to the myofibrillar actin filaments of the sarcomere.

In contrast, other Z-line associated proteins that link the sarcomere to the costamere [57] were decreased by NSAID treatment including DESM and CSRP3/MLP (mentioned above), zinc finger domain-containing proteins known as four-and-a-half LIM domain proteins 1 and 3 (FHL1, 3), PDZ and LIM domain protein 3 (PDLI3), and xin actin-binding repeat-containing protein 2 (XIRP2).

Apoptosis

Several anti-apoptosis proteins were significantly down-regulated by NSAID administration, including two members of the 14-3-3 family (1433E, 1433Z; [58]), translationally-controlled tumor protein (TCTP) [59], the mitochondrial proteins 3-ketoacyl-CoA thiolase (THIM; aka acetyl-Coenzyme A acyltransferase; [60]) and thioredoxin (THIOM; [61]), mitochondrial fission process protein 1 (MTFP1; aka mitochondrial 18kDa protein; [62]), and methylthioribulose-1-phosphate dehydratase (MTNB)–an inhibitor of the caspase 9/cytochrome c/APAF-1-dependent cell death pathway [63]. Calreticulin (CALR) was also down-regulated; cell-surface CALR serves as a general recognition ligand on apoptotic cells that facilitates their recognition and removal during resolution of inflammation [64]. Although only represented by a single peptide in our dataset, the novel anti-apoptosis protein Bcl-2-associated athanogene-1 (BAG1) was highly down-regulated by NSAID administration (>13-fold reduction; p = 0.002).

The NSAID-induced decrease in anti-apoptosis proteins was associated with a significant increase in caspase 3/7 enzymatic activity in injured TA muscle homogenates from 4 of 4 ketorolac-treated animals compared to non-injured muscle (mean increase = 23.7%, p = .02, Fig 5). In contrast, caspase 3/7 activity in injured TA muscles from saline-treated animals was significantly reduced in all animals tested (mean decrease = 26.2%; p = .005, Fig 5).

[Figure omitted. See PDF.]

Fig 5. NSAID-induced decrease in anti-apoptosis proteins is associated with increased caspase activity.

Caspase 3/7activity in muscle homogenates from saline- and NSAID-treated animals (4/group) was measured at 54 hrs after injury (7 hrs post-treatment) by commercial ELISA. Data are given as the fold-change for each animal relative to its own non-injured TA muscle. Responses between saline versus NSAID-treated animals were statistically significant as determined by ANOVA using a mixed model with random effect of animal within treatment with the level of significance set at p < .05.

https://doi.org/10.1371/journal.pone.0172486.g005

Discussion

Antecedent non-penetrating soft tissue injury often precedes the development of severe bacterial infection, including the ~50% of patients with group A streptococcal necrotizing fasciitis/myonecrosis who lack an obvious portal of bacterial entry (reviewed in [65]). Although some epidemiologic factors that increase the risk of death have been defined, the risk factors for initiation of infection are less clear. Non-penetrating trauma and use of non-steroidal anti-inflammatory drugs (NSAIDs) have each been associated with infection onset and/or worse outcomes in humans and in animal models of this infection (reviewed in [16]). For instance, in 2007, a retrospective study found that patients with GAS necrotizing fasciitis (but not those with GAS cellulitis) were 6 times more likely than matched controls to have a recent history of blunt trauma [66]. Further, a 2008 prospective epidemiologic study showed that NSAID use was independently associated with increased risk for severe GAS infection [67]. Our own work has shown that NSAID administration promotes trafficking of circulating GAS to sites of injury [18], greatly accelerates progression of established infection and reduces antibiotic efficacy [19].

To investigate how NSAIDs might predispose injured muscle to infection, the present study exploited our established animal model of eccentric contraction (EC)-induced muscle injury. Unlike previous reports in which NSAIDs were administered prior to experimental injury or excessive exercise, the present studies were designed to more closely mimic the human scenario in which NSAIDs are typically used during the peak of perceived muscle soreness–typically 24–48 hrs post-injury. At this point post-injury, markers of inflammation and muscle dysfunction are objectively measurable. For example, Lieber’s group has shown in mouse and rat models that muscle force-generating capacity remains significantly reduced at 24–48 hrs following experimental strain-associated muscle injury [23,24] and that this dysfunction is associated with marked architectural and inflammatory changes and in myoregulatory gene expression [23,24]. Results of the genetic, histologic and functional studies presented here are consistent with these reports and provide evidence that active muscle regeneration is ongoing at the time chosen for NSAID administration (47 hr post-injury) in our model.

Also within this critical 24–48 hr timeframe after injury, Arnold and colleagues have shown that infiltrated pro-inflammatory M1 macrophages convert to an anti-inflammatory/pro-resolution M2 phenotype to sustain myogenic differentiation, myofiber growth and membrane repair [25]. Because such phenotype switching is critical for muscle regeneration after injury and because prostaglandin metabolites (e.g. 15d- PGJ₂) promote the M2 phenotype [68], we first investigated whether delayed NSAID use might directly affect the M1 to M2 transition. Ketorolac, given after the peak of M1 influx, had no significant effect on the numbers or functional phenotypes of macrophages in the injured TA muscles. These data suggest that the signaling events promoting macrophage phenotype switching either occurred prior to NSAID administration or are NSAID-insensitive. Further, these results suggest that any NSAID-induced deficit in functional muscle regeneration is not related to effects on pro-resolution M2 leukocyte dynamics at this NSAID administration time point.

Like macrophages, muscle fibers can alter their functional phenotype under certain conditions such as changes in nerve supply, loading/unloading and aging [69]. Our findings suggest that NSAID use should be included on this list. Specifically, we demonstrate a NSAID-associated increase in slow MYH isoforms in the fast-twitch TA muscle which could, in part, explain the purported loss of muscle strength associated with NSAID use. The mechanism responsible is not apparent here, though the observed reduction in the thyroid hormone transport protein, transthyretin (TTHY; 3.19-fold decrease) could play a role [70] Further, by increasing MYOZ3, NSAID administration could activate the calcineurin/NFAT signaling pathway that drives transcription of slow fiber type-specific genes even in fast-twitch muscles [69]. An NSAID-induced transition of muscle fiber type could have clinical import in children, in those lacking normal muscle regenerative capacity, in patients with progressive loss of muscle mass (elderly, AIDS), and those undergoing rehabilitation after trauma or surgery. More research is required to confirm this notion.

Muscle repair/regeneration also requires repair of the saracolemmal membrane, revascularization and re-innervation of the myofiber, and reconstitution of the extracellular matrix [71]. Our proteomics data suggest that at 48–54 hrs after EC injury, muscle repair is well underway and is dominated by proliferation of muscle cell precursors and upregulation of stress-response/damage control mechanisms. Our data further suggest a damage-control, pro-survival process is also engaged. Each of these restorative processes was deleteriously affected by NSAID administration.

The NSAID-induced down-regulation of multiple inhibitors of the mitochondrial-based intrinsic apoptosis pathway was particularly intriguing. The 14-3-3 family of proteins, via binding to pro-apoptotic moieties such as Bcl-2-associated death promoter (BAD) and other ligands, is critical for cell survival signaling [58]. TCTP prevents apoptosis by inserting into the mitochondrial membrane and blocking Bax dimerization [72,73] and may also destabilize the tumor suppressor (pro-apoptosis) protein p53 [74,75]. Translational regulation of TCTP occurs via the PI3K/Akt/mTORC1 pathway [76]. Zhang and colleagues have recently shown that mTOR is essential for myogenesis [77]. Thus, our data may supply a missing link in the pathway between injury and regeneration and suggest that an injury/TCTP/mTOR axis drives this process. MTNB (also known as Apaf-1-interacting protein, APIP [78]) was initially identified by Cho et al as an inhibitor of hypoxia-induced intrinsic apoptosis in skeletal muscle [63] where it competes with caspase 9 for binding to Apaf-1 of the active apoptosome. APIP/MTNB also inhibits cytochrome c-induced activation of caspase-9 [63,79]. These activities are independent of its enzymatic function in methionine salvage [80]. Lastly, BAG1 was significantly reduced by NSAIDs. BAG1 is a novel multifunctional protein that was first identified via its ability to functionally augment the anti-apoptotic protein, Bcl-2 [81]. Recently, Warren et al demonstrated that contraction-induced injury, but not freeze injury, induced high level gene expression of a related protein, BAG3 [82]. In our study, the functional consequence of NSAID-induced down-regulation of these anti-apoptosis proteins was increased caspase activity in the injured muscles.

Interestingly, studies have suggested that controlled apoptosis is required for myogenic differentiation [83,84] and Fernando et al have implicated activated caspase-3 in this process. In their studies, caspase-3 activity was maximally (but transiently) increased in cultured murine C2C12 myoblasts after 24 hrs of serum withdrawal-induced differentiation. Biochemical blockade or genetic ablation of caspase-3 dramatically reduced myotube formation. Other data from this study suggested that the differentiation-associated caspase-3 activity was not strictly apoptotic since traditional markers of apoptosis (i.e., PARP cleavage, Annexin V staining) were absent or unchanged in these cells. In contrast to their findings, our data demonstrate an overall down-regulation of caspase activity in muscles undergoing regeneration at 54 hrs after strain injury (Fig 5, left). Further, our immunohistologic studies and western blot analysis of injured vs non-injured TA muscles from untreated mice failed to demonstrate an injury-induced increase in activated caspase-3 at this time (not shown). These apparently dichotomous results might be explained by differences in the model systems used (i.e., in vitro cultured cells with serum starvation vs whole muscles undergoing regeneration after strain injury) or the single time of sampling in our in vivo model (i.e., 54 hr post injury). We also demonstrate here that NSAID administration significantly increased caspase activity in regenerating muscles (Fig 5, right). Here too, we found no evidence of caspase-3 activation in tissues or tissue homogenates from NSAID-treated mice (not shown). In the absence of detectable caspase-3 activity and because our enzymatic detection assay measures both active caspase-3 and -7, it suggests a possible role for caspase-7 in the NSAID-induced response. Further, it is also conceivable that some activity attributed by Fernando et al to caspase-3 may in fact be related to caspase-7 since 1) these enzymes share overlapping substrate specificities [85] and 2) because the pharmacologic inhibitor used by these investigators (Z-DEVD.fmk) also inhibits caspase-7 [86]. Caspase-7 is an endoplasmic reticulum (ER)-associated caspase and can be regulated by the ER stress response protein known as glucose-regulated protein 78 (GRP78) [87] (also shown here to be reduced by NSAID administration). Thus, we hypothesize that GRP78-based inhibition of caspase-7 activity is required to ensure muscle regeneration after strain-induced injury and that NSAIDs delay repair by relieving this inhibition. Elucidation of this potential mechanism this requires further study.

Conclusions

In summary, our data provide new evidence suggesting that NSAID use during the peak of post-injury pain and inflammation decreases muscle metabolism, prevents neuromuscular junction and sarcolemmal repair, stimulates muscle fiber type transition and promotes a pro-apoptosis phenotype. Such findings provide molecular evidence supporting the notion that NSAIDs have a negative influence on myogenesis, including reducing muscle strength post-healing. By promoting cell death, NSAIDs could expand the nidus of injury, thereby increasing muscle susceptibility to post-injury infection [18,19]. These findings support renewed concerns about the risks versus benefits of NSAID use, especially in those with increased susceptibility to contraction-induced injury, those with muscle wasting or poor regenerative capacity, and those at risk for life-threatening bacterial myonecrosis. Given the potential to change the current paradigm of pain management in these settings, studies to validate these results in humans are warranted and must include NSAIDs of different classes given at different times post-injury.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Proteins Significantly Altered by Injury.

https://doi.org/10.1371/journal.pone.0172486.s001

(XLSX)

S2 Table. Proteins Significantly Altered by NSAID.

https://doi.org/10.1371/journal.pone.0172486.s002

(XLSX)

S1 Supporting Information. Additional Methodology for Proteomics Analysis.

https://doi.org/10.1371/journal.pone.0172486.s003

(DOCX)

Acknowledgments

The authors wish to thank Kara Schmitt and Kristin Burnum-Johnson for expert technical assistance.

Author Contributions

1. Conceptualization: AEB TOM BJWR DLS.

2. Data curation: AEB MJA TOM BJWR.

3. Formal analysis: AEB TOM BJWR LB.

4. Funding acquisition: AEB MJA.

5. Investigation: AEB MJA CRB EJK CDN TLF TRWC TOM LB BJWR.

6. Methodology: AEB MJA CRB TOM BJWR.

7. Resources: AEB TOM BJWR LB DLS.

8. Supervision: AEB TOM BJWR.

9. Visualization: AEB.

10. Writing – original draft: AEB MJA TOM BJWR.

11. Writing – review & editing: AEB MJA CRB TOM BJWR LB DLS.

Citation: Bryant AE, Aldape MJ, Bayer CR, Katahira EJ, Bond L, Nicora CD, et al. (2017) Effects of delayed NSAID administration after experimental eccentric contraction injury – A cellular and proteomics study. PLoS ONE 12(2): e0172486. https://doi.org/10.1371/journal.pone.0172486

References

1. US Food and Drug Administration (2015) Safety Announcement: FDA strengthens warning that non-aspirin nonsteroidal anti-inflammatory drugs (NSAIDs) can cause heart attaches or strokes.

2. Hakkarainen TW, Steele SR, Bastaworous A, Dellinger EP, Farrokhi E, Farjah F et al. (2015) Nonsteroidal anti-inflammatory drugs and the risk for anastomotic failure: a report from Washington State's Surgical Care and Outcomes Assessment Program (SCOAP). JAMA Surg 150: 223–228. pmid:25607250

3. Lieber RL, Friden J (2002) Morphologic and mechanical basis of delayed-onset muscle soreness. J Am Acad Orthop Surg 10: 67–73. pmid:11809052

4. Mackey AL, Mikkelsen UR, Magnusson SP, Kjaer M (2012) Rehabilitation of muscle after injury—the role of anti-inflammatory drugs. Scand J Med Sci Sports 22: e8–14. pmid:22449131

5. Bondesen BA, Mills ST, Kegley KM, Pavlath GK (2004) The COX-2 pathway is essential during early stages of skeletal muscle regeneration. Am J Physiol Cell Physiol 287: C475–C483. pmid:15084473

6. Mendias CL, Tatsumi R, Allen RE (2004) Role of cyclooxygenase-1 and -2 in satellite cell proliferation, differentiation, and fusion. Muscle Nerve 30: 497–500. pmid:15372441

7. Mikkelsen UR, Langberg H, Helmark IC, Skovgaard D, Andersen LL, Kjaer M et al. (2009) Local NSAID infusion inhibits satellite cell proliferation in human skeletal muscle after eccentric exercise. J Appl Physiol (1985) 107: 1600–1611. pmid:19713429

8. Mackey AL, Kjaer M, Dandanell S, Mikkelsen KH, Holm L, Dossing S et al. (2007) The influence of anti-inflammatory medication on exercise-induced myogenic precursor cell responses in humans. J Appl Physiol 103: 425–431. pmid:17463304

9. Trappe TA, White F, Lambert CP, Cesar D, Hellerstein M, Evans WJ (2002) Effect of ibuprofen and acetaminophen on postexercise muscle protein synthesis. Am J Physiol Endocrinol Metab 282: E551–E556. pmid:11832356

10. Mishra DK, Friden J, Schmitz MC, Lieber RL (1995) Anti-inflammatory medication after muscle injury. A treatment resulting in short-term improvement but subsequent loss of muscle function. J Bone Joint Surg Am 77: 1510–1519. pmid:7593059

11. Permpalung N, Upala S, Sanguankeo A, Sornprom S (2016) Association between NSAIDs and Clostridium difficile-Associated Diarrhea: A Systematic Review and Meta-Analysis. Can J Gastroenterol Hepatol 2016: 7431838. pmid:27446866

12. Dial S, Delaney JA, Barkun AN, Suissa S (2005) Use of gastric acid-suppressive agents and the risk of community-acquired Clostridium difficile-associated disease. JAMA 294: 2989–2995. pmid:16414946

13. Leroy S, Marc E, Bavoux F, Treluyer JM, Gendrel D, Breart G et al. (2010) Hospitalization for severe bacterial infections in children after exposure to NSAIDs: a prospective adverse drug reaction reporting study. Clin Drug Investig 30: 179–185. pmid:20155990

14. Messika J, Sztrymf B, Bertrand F, Billard-Pomares T, Barnaud G, Branger C et al. (2014) Risks of nonsteroidal antiinflammatory drugs in undiagnosed intensive care unit pneumococcal pneumonia: younger and more severely affected patients. J Crit Care 29: 733–738. pmid:24997726

15. Souyri C, Olivier P, Grolleau S, Lapeyre-Mestre M (2008) Severe necrotizing soft-tissue infections and nonsteroidal anti-inflammatory drugs. Clin Exp Dermatol 33: 249–255. pmid:18261144

16. Bryant AE, Bayer CR, Aldape MJ, Stevens DL (2015) The roles of injury and nonsteroidal anti-inflammatory drugs in the development and outcomes of severe group A streptococcal soft tissue infections. Curr Opin Infect Dis 28: 231–239. pmid:25918957

17. Stevens DL, Tanner MH, Winship J, Swarts R, Reis KM, Schlievert PM et al. (1989) Reappearance of scarlet fever toxin A among streptococci in the Rocky Mountain West: Severe group A streptococcal infections associated with a toxic shock-like syndrome. N Eng J Med 321(1): 1–7.

18. Hamilton SM, Bayer CR, Stevens DL, Lieber RL, Bryant AE (2008) Muscle injury, vimentin expression, and nonsteroidal anti-inflammatory drugs predispose to cryptic group A streptococcal necrotizing infection. J Infect Dis 198: 1692–1698. pmid:18939933

19. Hamilton SM, Bayer CR, Stevens DL, Bryant AE (2014) Effects of selective and nonselective nonsteroidal anti-inflammatory drugs on antibiotic efficacy of experimental group A streptococcal myonecrosis. J Infect Dis 209: 1429–1435. pmid:24218498

20. Paulsen G, Mikkelsen UR, Raastad T, Peake JM (2012) Leucocytes, cytokines and satellite cells: what role do they play in muscle damage and regeneration following eccentric exercise? Exerc Immunol Rev 18: 42–97. pmid:22876722

21. Volker D, Bate M, Gentle R, Garg M (2000) Oral buprenorphine is anti-inflammatory and modulates the pathogenesis of streptococcal cell wall polymer-induced arthritis in the Lew/SSN rat. Lab Anim 34: 423–429. pmid:11072864

22. Brooks PM, Day RO (1991) Nonsteroidal antiinflammatory drugs—differences and similarities. N Engl J Med 324: 1716–1725. pmid:2034249

23. Barash IA, Mathew L, Ryan AF, Chen J, Lieber RL (2004) Rapid muscle-specific gene expression changes after a single bout of eccentric contractions in the mouse. Am J Physiol Cell Physiol 286: C355–C364. pmid:14561590

24. Peters D, Barash IA, Burdi M, Yuan PS, Mathew L, Friden J et al. (2003) Asynchronous functional, cellular and transcriptional changes after a bout of eccentric exercise in the rat. J Physiol 553: 947–957. pmid:14514871

25. Arnold L, Henry A, Poron F, Baba-Amer Y, van RN, Plonquet A et al. (2007) Inflammatory monocytes recruited after skeletal muscle injury switch into antiinflammatory macrophages to support myogenesis. J Exp Med 204: 1057–1069. pmid:17485518

26. Zimmer JS, Monroe ME, Qian WJ, Smith RD (2006) Advances in proteomics data analysis and display using an accurate mass and time tag approach. Mass Spectrom Rev 25: 450–482. pmid:16429408

27. Yang F, Shen Y, Camp DG, Smith RD (2012) High-pH reversed-phase chromatography with fraction concatenation for 2D proteomic analysis. Expert Rev Proteomics 9: 129–134. pmid:22462785

28. Kiebel GR, Auberry KJ, Jaitly N, Clark DA, Monroe ME, Peterson ES et al. (2006) PRISM: a data management system for high-throughput proteomics. Proteomics 6: 1783–1790. pmid:16470653

29. Jaitly N, Mayampurath A, Littlefield K, Adkins JN, Anderson GA, Smith RD (2009) Decon2LS: An open-source software package for automated processing and visualization of high resolution mass spectrometry data. BMC Bioinformatics 10: 87. pmid:19292916

30. Monroe ME, Tolic N, Jaitly N, Shaw JL, Adkins JN, Smith RD (2007) VIPER: an advanced software package to support high-throughput LC-MS peptide identification. Bioinformatics 23: 2021–2023. pmid:17545182

31. Schulz-Trieglaff O, Machtejevas E, Reinert K, Schluter H, Thiemann J, Unger K (2009) Statistical quality assessment and outlier detection for liquid chromatography-mass spectrometry experiments. BioData Min 2: 4. pmid:19351414

32. Webb-Robertson BJ, McCue LA, Waters KM, Matzke MM, Jacobs JM, Metz TO et al. (2010) Combined Statistical Analyses of Peptide Intensities and Peptide Occurrences Improves Identification of Significant Peptides from MS-Based Proteomics Data. J Proteome Res.

33. Webb-Robertson BJ, Matzke MM, Datta S, Payne SH, Kang J, Bramer LM et al. (2014) Bayesian proteoform modeling improves protein quantification of global proteomic measurements. Mol Cell Proteomics 13: 3639–3646. pmid:25433089

34. Polpitiya AD, Qian WJ, Jaitly N, Petyuk VA, Adkins JN, Camp DG et al. (2008) DAnTE: a statistical tool for quantitative analysis of -omics data. Bioinformatics 24: 1556–1558. pmid:18453552

35. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV et al. (2003) The COG database: an updated version includes eukaryotes. BMC Bioinformatics 4: 41. pmid:12969510

36. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25: 402–408. pmid:11846609

37. Sidney S (1956) Nonparametric Statistics for the Behavioral Sciences.

38. StatXact Software Version 8.0 (2007) [Computer Software].

39. Augusto V, Padovani CR, Campos GER (2004) Skeletal muscle fiber types in C67BL6J mice. Braz J Morphol Sci 21: 89–94.

40. Choi SJ (2014) Differential susceptibility on myosin heavy chain isoform following eccentric-induced muscle damage. J Exerc Rehabil 10: 344–348. pmid:25610817

41. Paulsen G, Mikkelsen UR, Raastad T, Peake JM (2012) Leucocytes, cytokines and satellite cells: what role do they play in muscle damage and regeneration following eccentric exercise? Exerc Immunol Rev 18: 42–97. pmid:22876722

42. Vater R, Cullen MJ, Harris JB (1994) The expression of vimentin in satellite cells of regenerating skeletal muscle in vivo. Histochem J 26: 916–928. pmid:7896567

43. Martineau LC, Gardiner PF (2001) Insight into skeletal muscle mechanotransduction: MAPK activation is quantitatively related to tension. J Appl Physiol 91: 693–702. pmid:11457783

44. Tidball JG (2005) Inflammatory processes in muscle injury and repair. Am J Physiol Regul Integr Comp Physiol 288: R345–R353. pmid:15637171

45. Belcastro AN, Shewchuk LD, Raj DA (1998) Exercise-induced muscle injury: a calpain hypothesis. Mol Cell Biochem 179: 135–145. pmid:9543356

46. Linares GR, Xing W, Govoni KE, Chen ST, Mohan S (2009) Glutaredoxin 5 regulates osteoblast apoptosis by protecting against oxidative stress. Bone 44: 795–804. pmid:19442627

47. Trappe TA, Liu SZ (2013) Effects of prostaglandins and COX-inhibiting drugs on skeletal muscle adaptations to exercise. J Appl Physiol (1985) 115: 909–919. pmid:23539318

48. Vaittinen S, Lukka R, Sahlgren C, Hurme T, Rantanen J, Lendahl U et al. (2001) The expression of intermediate filament protein nestin as related to vimentin and desmin in regenerating skeletal muscle. J Neuropathol Exp Neurol 60: 588–597. pmid:11398835

49. Sejersen T, Lendahl U (1993) Transient expression of the intermediate filament nestin during skeletal muscle development. J Cell Sci 106 (Pt 4): 1291–1300.

50. Bryant AE, Bayer CR, Huntington JD, Stevens DL (2006) Group A streptococcal myonecrosis: increased vimentin expression after skeletal-muscle injury mediates the binding of Streptococcus pyogenes. J Infect Dis 193: 1685–1692. pmid:16703512

51. Fuchs E, Weber K (1994) Intermediate filaments: structure, dynamics, function, and disease. Annu Rev Biochem 63: 345–382. pmid:7979242

52. Bean C, Verma NK, Yamamoto DL, Chemello F, Cenni V, Filomena MC et al. (2014) Ankrd2 is a modulator of NF-kappaB-mediated inflammatory responses during muscle differentiation. Cell Death Dis 5: e1002. pmid:24434510

53. Vafiadaki E, Arvanitis DA, Sanoudou D (2015) Muscle LIM Protein: Master regulator of cardiac and skeletal muscle functions. Gene 566: 1–7. pmid:25936993

54. Posey AD Jr., Demonbreun A, McNally EM (2011) Ferlin proteins in myoblast fusion and muscle growth. Curr Top Dev Biol 96: 203–230. pmid:21621072

55. Cai C, Masumiya H, Weisleder N, Matsuda N, Nishi M, Hwang M et al. (2009) MG53 nucleates assembly of cell membrane repair machinery. Nat Cell Biol 11: 56–64. pmid:19043407

56. McDade JR, Michele DE (2014) Membrane damage-induced vesicle-vesicle fusion of dysferlin-containing vesicles in muscle cells requires microtubules and kinesin. Hum Mol Genet 23: 1677–1686. pmid:24203699

57. Ervasti JM (2003) Costameres: the Achilles' heel of Herculean muscle. J Biol Chem 278: 13591–13594. pmid:12556452

58. Masters SC, Fu H (2001) 14-3-3 proteins mediate an essential anti-apoptotic signal. J Biol Chem 276: 45193–45200. pmid:11577088

59. Nagano-Ito M, Ichikawa S (2012) Biological effects of Mammalian translationally controlled tumor protein (TCTP) on cell death, proliferation, and tumorigenesis. Biochem Res Int 2012: 204960. pmid:22675633

60. Cao W, Liu N, Tang S, Bao L, Shen L, Yuan H et al. (2008) Acetyl-Coenzyme A acyltransferase 2 attenuates the apoptotic effects of BNIP3 in two human cell lines. Biochim Biophys Acta 1780: 873–880. pmid:18371312

61. Nonn L, Williams RR, Erickson RP, Powis G (2003) The absence of mitochondrial thioredoxin 2 causes massive apoptosis, exencephaly, and early embryonic lethality in homozygous mice. Mol Cell Biol 23: 916–922. pmid:12529397

62. Tondera D, Santel A, Schwarzer R, Dames S, Giese K, Klippel A et al. (2004) Knockdown of MTP18, a novel phosphatidylinositol 3-kinase-dependent protein, affects mitochondrial morphology and induces apoptosis. J Biol Chem 279: 31544–31555. pmid:15155745

63. Cho DH, Hong YM, Lee HJ, Woo HN, Pyo JO, Mak TW et al. (2004) Induced inhibition of ischemic/hypoxic injury by APIP, a novel Apaf-1-interacting protein. J Biol Chem 279: 39942–39950. pmid:15262985

64. Gardai SJ, McPhillips KA, Frasch SC, Janssen WJ, Starefeldt A, Murphy-Ullrich JE et al. (2005) Cell-surface calreticulin initiates clearance of viable or apoptotic cells through trans-activation of LRP on the phagocyte. Cell 123: 321–334. pmid:16239148

65. Bisno AL, Stevens DL (1996) Streptococcal infections in skin and soft tissues. N Engl J Med 334: 240–245. pmid:8532002

66. Nuwayhid ZB, Aronoff DM, Mulla ZD (2007) Blunt trauma as a risk factor for group A streptococcal necrotizing fasciitis. Ann Epidemiol 17: 878–881. pmid:17697787

67. Lamagni TL, Neal S, Keshishian C, Alhaddad N, George R, Duckworth G et al. (2008) Severe Streptococcus pyogenes Infections, United Kingdom, 2003–2004. Emerg Infect Dis 14: 201–209.

68. Alleva DG, Johnson EB, Lio FM, Boehme SA, Conlon PJ, Crowe PD (2002) Regulation of murine macrophage proinflammatory and anti-inflammatory cytokines by ligands for peroxisome proliferator-activated receptor-gamma: counter-regulatory activity by IFN-gamma. J Leukoc Biol 71: 677–685. pmid:11927655

69. Pette D, Staron RS (2001) Transitions of muscle fiber phenotypic profiles. Histochem Cell Biol 115: 359–372. pmid:11449884

70. Soukup T, Smerdu V (2014) Effect of altered innervation and thyroid hormones on myosin heavy chain expression and fiber type transitions: a mini-review. Histochem Cell Biol.

71. Charge SB, Rudnicki MA (2004) Cellular and molecular regulation of muscle regeneration. Physiol Rev 84: 209–238. pmid:14715915

72. Susini L, Besse S, Duflaut D, Lespagnol A, Beekman C, Fiucci G et al. (2008) TCTP protects from apoptotic cell death by antagonizing bax function. Cell Death Differ 15: 1211–1220. pmid:18274553

73. Susini L, Besse S, Duflaut D, Lespagnol A, Beekman C, Fiucci G et al. (2008) TCTP protects from apoptotic cell death by antagonizing bax function. Cell Death Differ 15: 1211–1220. pmid:18274553

74. Rho SB, Lee JH, Park MS, Byun HJ, Kang S, Seo SS et al. (2011) Anti-apoptotic protein TCTP controls the stability of the tumor suppressor p53. FEBS Lett 585: 29–35. pmid:21081126

75. Amson R, Pece S, Lespagnol A, Vyas R, Mazzarol G, Tosoni D et al. (2012) Reciprocal repression between P53 and TCTP. Nat Med 18: 91–99.

76. Bommer UA, Iadevaia V, Chen J, Knoch B, Engel M, Proud CG (2015) Growth-factor dependent expression of the translationally controlled tumour protein TCTP is regulated through the PI3-K/Akt/mTORC1 signalling pathway. Cell Signal 27: 1557–1568. pmid:25936523

77. Zhang P, Liang X, Shan T, Jiang Q, Deng C, Zheng R et al. (2015) mTOR is necessary for proper satellite cell activity and skeletal muscle regeneration. Biochem Biophys Res Commun 463: 102–108. pmid:25998386

78. Mary C, Duek P, Salleron L, Tienz P, Bumann D, Bairoch A et al. (2012) Functional identification of APIP as human mtnB, a key enzyme in the methionine salvage pathway. PLoS One 7: e52877. pmid:23285211

79. Cho DH, Lee HJ, Kim HJ, Hong SH, Pyo JO, Cho C et al. (2007) Suppression of hypoxic cell death by APIP-induced sustained activation of AKT and ERK1/2. Oncogene 26: 2809–2814. pmid:17086211

80. Kang W, Hong SH, Lee HM, Kim NY, Lim YC, Le le TM et al. (2014) Structural and biochemical basis for the inhibition of cell death by APIP, a methionine salvage enzyme. Proc Natl Acad Sci U S A 111: E54–E61. pmid:24367089

81. Takayama S, Sato T, Krajewski S, Kochel K, Irie S, Millan JA et al. (1995) Cloning and functional analysis of BAG-1: a novel Bcl-2-binding protein with anti-cell death activity. Cell 80: 279–284. pmid:7834747

82. Warren GL, Summan M, Gao X, Chapman R, Hulderman T, Simeonova PP (2007) Mechanisms of skeletal muscle injury and repair revealed by gene expression studies in mouse models. J Physiol 582: 825–841. pmid:17478534

83. Dee K, Freer M, Mei Y, Weyman CM (2002) Apoptosis coincident with the differentiation of skeletal myoblasts is delayed by caspase 3 inhibition and abrogated by MEK-independent constitutive Ras signaling. Cell Death Differ 9: 209–218. pmid:11840171

84. Fernando P, Kelly JF, Balazsi K, Slack RS, Megeney LA (2002) Caspase 3 activity is required for skeletal muscle differentiation. Proc Natl Acad Sci U S A 99: 11025–11030. pmid:12177420

85. Walsh JG, Cullen SP, Sheridan C, Luthi AU, Gerner C, Martin SJ (2008) Executioner caspase-3 and caspase-7 are functionally distinct proteases. Proc Natl Acad Sci U S A 105: 12815–12819. pmid:18723680

86. Margolin N, Raybuck SA, Wilson KP, Chen W, Fox T, Gu Y et al. (1997) Substrate and inhibitor specificity of interleukin-1 beta-converting enzyme and related caspases. J Biol Chem 272: 7223–7228. pmid:9054418

87. Reddy RK, Mao C, Baumeister P, Austin RC, Kaufman RJ, Lee AS (2003) Endoplasmic reticulum chaperone protein GRP78 protects cells from apoptosis induced by topoisomerase inhibitors: role of ATP binding site in suppression of caspase-7 activation. J Biol Chem 278: 20915–20924. pmid:12665508