Cancer and Associated Therapies Impact the Skeletal Muscle Proteome

Introduction: Both cancer and cancer associated therapies (CAT; including chemotherapy or concurrent chemoradiation) disrupt cellular metabolism throughout the body, including the regulation of skeletal muscle mass and function. Adjunct testosterone therapy during standard of care chemotherapy and chemoradiation modulates CAT-induced dysregulation of skeletal muscle metabolism and protects lean body mass during CAT. However, the extent to which the skeletal muscle proteome is altered under these therapeutic conditions is unknown.

Objective: We probed the skeletal muscle proteome of cancer patients as an ancillary analysis following a randomized, double-blind, placebo-controlled phase II trial investigating the effect of adjunct testosterone on body composition in men and women with advanced cancers undergoing CAT.

Methods: Men and women diagnosed with late stage (≥IIB) or recurrent head and neck or cervical cancer who were scheduled to receive standard of care CAT were administered an adjunct 7 weeks treatment of weekly intramuscular injections of either 100 mg testosterone (CAT+T, n = 7; 2M/5F) or placebo/saline (CAT+P, n = 6; 4M/2F). Biopsies were performed on the vastus lateralis before (PRE) and after (POST) the 7 weeks treatment. Extracted proteins were separated with 2-dimensional gel electrophoresis (2DE), and subjected to analyses of total protein abundance, phosphorylation and S-nitrosylation. Proteoforms showing significant 1.5 fold differences (t-test p ≤ 0.05) between PRE and POST timepoints were identified by mass spectroscopy (MS), and lists of altered proteins were subjected to Gene Set Enrichment Analysis (GSEA) to identify affected pathways.

Results: A total of 756 distinct protein spots were identified. Of those spots, 102 were found to be altered in terms of abundance, phosphorylation, or S-nitrosylation, and identified by mass spectroscopy analysis to represent 58 unique proteins. Among the biological processes and pathways identified, CAT+P predominantly impacted metabolic processes, cell assembly, oxygen transport, and apoptotic signaling, while CAT+T impacted transcription regulation, muscle differentiation, muscle development, and contraction.

Conclusion: Cancer and CAT significantly altered the skeletal muscle proteome in a manner suggestive of loss of structural integrity, reduced contractile function, and disrupted metabolism. Proteomic analysis suggests that the addition of adjunct testosterone minimized the structural and contractile influence of cancer and its associated therapies.

Introduction

Cancer progression is often associated with cachexia, a debilitating, and phenotypically recognizable loss in body mass and physical function. In comparison to several other common catabolic conditions including prolonged physical inactivity and sarcopenia of aging, cancer cachexia is a mechanically distinct multisystem wasting condition that extends beyond the typical losses in muscle mass and strength. Cancer cachexia is characterized by wasting of organs and tissues that includes atrophy of skeletal muscle mass as well as adipose tissue, whereas conditions such as sarcopenia of aging or sedentary lifestyles are typically accompanied by increased adiposity along with more gradual declines in skeletal mass (Ali and Garcia, 2014).

Distinctions between catabolic conditions have been observed at the molecular level, indicating that while various conditions can disrupt skeletal muscle proteostasis and trigger skeletal muscle atrophy, the catabolic triggers have distinct proteomic patterns (Hunt et al., 2021). Genomic and proteomic analyses of mouse muscle by Hunt et al. support the current understanding that more acute drivers of atrophy such as corticosteroids and cachexia are characterized by decreased protein synthesis and increased degradation while age related muscle loss involves reduced maintenance of protein quality but not necessarily increased protein degradation (Hunt et al., 2021). Characterizing and comparing the differences and similarities in the skeletal muscle proteome between distinct wasting conditions will help provide insight into the mechanisms that drive muscle atrophy and inform development of targeted interventions against cachexia (Peixoto da Silva et al., 2020). Given the frequency of cachexia in people with cancer it is important to understand the proteomic profile of skeletal muscle in this population, including how interventions affect the abundance and posttranslational modification of specific proteins as well as the pathways associated with these changes. Confounding factors (including but not limited to the type of cancer, recurrence or progression of the disease, type or length of standard of care cancer treatment, patient age, and sex of the patient) which affect the progression of cancer cachexia, are not always well-controlled and can vary immensely.

Key intervention strategies against losses of muscle mass and strength encountered as a consequence of sedentary lifestyles, and even sarcopenia, include modulation of diet and physical activity. While immensely important as part of the therapeutic approach during cancer treatment, these interventions alone have not been proven to be adequate in countering the progression of cachexia (Schütte et al., 2020; Kasvis and Kilgour, 2021; Kim et al., 2021; van der Meij et al., 2021). In addition, cancer associated therapies (CAT) including chemo and radiation, inevitably add to the complexity when exploring the mechanisms that drive cancer cachexia as these treatments are themselves physiologically disruptive to cellular processes (Damrauer et al., 2018; Pin et al., 2019). Thus, while these treatments are clinically beneficial for targeting cancerous cells, improving patient quality of life, and aiding survival, the treatment administered to the patient becomes part of the cellular milieu and may affect key biological processes.

Hormones, including androgens such as testosterone which are widely known for the anabolic effect on skeletal muscle, have proven clinical efficacy with treatment approaches demonstrated in male as well as female populations (Sheffield-Moore et al., 2006; Horstman et al., 2012; Dos Santos and Storer, 2022). However, testosterone treatment as an intervention against muscle wasting in cancer patients is novel. We recently investigated the effect of adjunct testosterone on body composition and quality of life in men and women undergoing CAT for late-stage squamous cell carcinoma of the cervix or head and neck (Wright et al., 2018). This National Cancer Institute (NCI) funded, randomized, double-blind, placebo-controlled phase II trial demonstrated that 7 weeks of testosterone supplementation ameliorated weight loss in late stage cancer patients by preferentially increasing lean body mass. The improvements in lean body mass were associated with clinically meaningful improvements in patient quality of life and physical performance, although adjuvant testosterone therapy did not increase overall patient survival. To better understand how the course of standard of care CAT with or without adjuvant testosterone treatment differentially alters the skeletal muscle proteome, we compared skeletal muscle biopsies from the patients before and after respective treatment.

Methods

Patients/Subjects

The study complied with the Declaration of Helsinki and was approved by The University of Texas Medical Branch (UTMB) Institutional Review Board (IRB). Written informed consent was obtained from all subjects prior to being studied at the ITS-CRC (Institute for Translational Sciences - Clinical Research Center) at UTMB. Details of the study design and primary outcomes from this intent-to-treat investigation have been published previously (Wright et al., 2018) and the data are made available at https://doi.org/10.18738/T8/RFRQE2. Only subjects with both pre and post treatment muscle tissue samples available for proteomics were included in these ancillary analyses. These analyses include data from 13 squamous cell carcinoma patients enrolled in a double-blinded testosterone vs. placebo treatment study funded by the NCI and was designed to not interfere with standard of care treatment. Volunteers included patients with advanced or recurrent squamous cell carcinoma of the cervix (stages IIB, IIIA, and IIIB) between the ages of 18 and 65, and patients with advanced (stage III or IV) or recurrent head and neck squamous cell carcinoma between the ages of 18 and 75 (Table 1). The CAT+T group includes three cervical cancer patients and four head and neck cancer patients. One cervical cancer patient also followed an amino acid supplementation regimen in addition to receiving CAT+T treatment and POST study data collection occurred 3 months after PRE study data collection due to interruptions and postponements in her treatment schedule. Because of the intent-to-treat nature of the investigation, this patient was not excluded from participation or analysis. Pre-post analyses for CAT+T include five female and two male patients (age 53.6 ± 5.1 years, BMI 21.7 ± 8.0 kg/m²). Tissues from four CAT+T subjects were available for S-nitrosylation analyses. The CAT+P group includes two cervical cancer patients and four head and neck cancer patients. Tissues from four CAT+P subjects were included in S-nitrosylation analyses. Pre-post analyses for CAT+P group include two female and four male patients (age 49.7 ± 10.3 years, BMI 25.9 ± 8.0 kg/m²).

TABLE 1

TABLE 1. Patient population. Baseline demographic information for cancer patients receiving either cancer associated therapy with adjunct testosterone (CAT+T) or placebo (CAT+P). Based on muscle biopsy tissue availability, samples were analyzed for protein abundance, phosphorylation, and S-nitrosylation. H&N, Head and Neck.

Testosterone Treatment

The testosterone treatment followed a common replacement schedule used to treat hypogonadal men. Patients received weekly intramuscular injections of either 100 mg testosterone enanthate or sterile saline (placebo) over a period of 7 weeks. Testosterone and placebo injections were given by a nurse using an opaque syringe to obscure visual differences between testosterone and placebo. Because of the late stages of cancer examined in this study, both men and women were treated with the same dose of testosterone as this dose and duration would not be expected to affect female secondary sexual characteristics.

Muscle Biopsy Procedure

Muscle tissue was collected before (PRE) and after (POST) completion of standard of care chemo/radiation therapy during which time the patients were randomized to also receive either placebo (CAT+P) or testosterone (CAT+T) treatment. Muscle biopsy procedures were performed as described elsewhere (5, 6). Briefly, a site was marked over the vastus lateralis and cleaned with Betadine. Lidocaine (1%) was administered to the skin and muscle. An approximately 5 mm incision was made through the skin and fascia and a 5 mm Bergström needle was advanced into the muscle. While suction was applied, 100–200 mg of skeletal muscle tissue was collected by opening and closing the cutting window of the biopsy needle. The incision was sutured and covered with Bacitracin and steri-strips. Ice was applied to the site and ibuprofen was provided to the subject to alleviate soreness.

Skeletal Muscle Proteomics

Abundance and Phosphorylation - Proteomic analyses were performed by the UTMB Biomolecular Resource Facility (BRF) in fractionated muscle extracts using a Biofluids Analytical Platform (BAP) as previously described (Sheffield-Moore et al., 2013; Dillon et al., 2019). These analyses were completed in one continuous effort once all the muscle samples had been collected. The BAP fractionation component combines Superdex S-75 size-exclusion chromatography (SEC) of biofluids with electronically triggered fraction collection to create protein and peptide pools for subsequent separation and analysis. Fractionated samples were subjected to 2D gel electrophoresis (2DE) and stained for phosphoproteins (ProQ Diamond dye) or total proteins (Sypro Ruby dye). The gels were imaged, and analyzed using SameSpots 2D Gel Image Analysis Software for Differential Expression (TotalLab, Newcastle upon Tyne, United Kingdom). Images were aligned using select reference spots and log-transformed spot intensities were quantitatively compared between timepoints within groups (PRE vs. POST). Spots with significant changes (uncorrected p ≤ 0.05 and absolute fold differences ≥1.50) between PRE and POST time points identified by SameSpots within either group (CAT+P or CAT+T) were identified by Matrix Assisted Laser Desorption Ionization Time of Flight Mass Spectrometry (MALDI TOF MS), and compared against the Uniprot protein database for identification.

Nitrosylation—Analyses of the S-nitrosoproteome were performed using a fluorescence saturation technique (SNOFlo) separately from the other skeletal muscle proteome analyses (abundances and phosphorylation) as described previously using a ratio of ratios method (Wiktorowicz et al., 2011; Sheffield-Moore et al., 2013; Wiktorowicz et al., 2017). Briefly, after quantification of cysteine content (cysteic acid), denatured sample aliquots were split and either labeled with Bodipy Fl-maleimide (BD) or treated with ascorbate to reverse S-nitrosylation in the samples. Both sets of samples were dialyzed against denaturing buffer followed by BD labeling of the ascorbate treated samples. After two-dimensional gel electrophoresis, spot quantification by fluorescence imaging, and image analysis, protein abundance was measured by calculating the ratio of spot volumes on PRE versus POST gels of the ascorbate-treated samples. The ratios for the samples not treated with ascorbate were calculated similarly. The levels of S-nitrosylation were then expressed as the ratio of the two ratios (ratio of ratios), with the non-ascorbate-treated ratios (representing the sum of S-nitrosylated and abundance differences) normalized against the ascorbate-treated ratios (abundance difference only). A negative ratio of ratios indicates increased cysteinyl S-nitrosylation, and a positive ratio of ratios decreased S-nitrosylation.

Data and Statistical Analyses

To identify potential biological processes and molecular pathways that were altered within each group, the sets of significantly altered protein IDs from each grouping were subjected to gene set enrichment analysis (GSEA) based on the 2021 Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using the online Enrichr platform (https://amp.pharm.mssm.edu/Enrichr/) (Chen et al., 2013; Kuleshov et al., 2016; Xie et al., 2021). Significantly altered biological pathways were identified by an adjusted p < 0.05 using the Fisher exact test, and Benjamini–Hochberg method for correction for multiple hypotheses testing. Data were visualized based on an Enrichr combined score calculated by multiplying the log of the p-value from the Fisher exact test by the z-score of the deviation from the expected rank. Separate analyses were performed to identify altered pathways based on significantly altered protein abundance, S-nitrosylation, and phosphorylation. The processes and pathways identified through Enricher were interpreted in context of the phenotypical changes observed in the placebo vs. testosterone treated groups as previously reported (Wright et al., 2018).

Results

Patients

Changes in patient body composition and strength were main outcomes of the original study and details have been reported previously (Wright et al., 2018). Briefly, total body weight decreased in the placebo treated group and remained unchanged during testosterone treatment. Lean mass decreased with placebo treatment while this increased with testosterone treatment. Fat mass decreased in both the placebo and testosterone groups. There were no significant changes in muscle strength or physical performance in either group. There were no differences in 1-year survival between the groups.

Proteomics

Of 756 distinct spots, 102 were determined to be significantly altered (>1.50 fold change AND p < 0.05) from pre to post treatment in abundance (Placebo, n = 7 spots; Testosterone, n = 8 spots), phosphorylation (Placebo, n = 6 spots; Testosterone, n = 7 spots), and/or S-nitrosylation (Placebo, n = 59 spots; Testosterone, n = 26 spots), and identified using MS. These significantly altered spots were identified by MS, and found to represent 58 unique proteins (Figure 1). To facilitate interpretation, the spots were categorized based on the predicted molecular weights (MW) published in the UNIPROT database (www.uniprot.org). Of the 102 spots, 41 spots reflected proteoforms with MWs that corresponded to full-size proteins (within 15% of predicted UNIPROT MW), 22 spots were assumed to reflect protein dimers and aggregates (exceeded the predicted MW by at least 15%), and 39 spots were assumed to be fragments (more than 15% below predicted MW).

FIGURE 1

FIGURE 1. Heat map. Heat map depicting the fold change in abundance, phosphorylation, and S-nitrosylation of skeletal muscle proteins in cancer patients following 7 weeks of treatment with testosterone (T) or placebo (P). Proteoforms marked with “X” changed significantly (p < 0.05) from pre to post treatment. Proteoforms are labeled (central column) as either intact (I), fragments (F), or aggregate (A) based on identification and predicted molecular weight as described in the text.

Protein Abundance Changes

In the CAT+P (placebo) group, cancer and related therapy was associated with significant changes (POST vs. PRE) in the abundance of seven proteoforms including four intact proteins, one aggregate, and two fragments (Figure 1, Supplementary Table S1). Significantly altered intact proteins included those involved with muscle contraction (MYL3, MYL1), cellular structure (DES), and calcium regulation (ANXA2). Of the intact proteoforms that changed, MYL3 and MYL1 decreased in abundance while DES and ANXA2 increased in abundance. Cancer combined with standard of care therapy plus testosterone treatment (CAT+T) was associated with changes in abundance of two intact proteoforms, one aggregate, and five fragments (Figure 1). Altered intact proteins included increased abundance of the structural protein ANKRD2 and decreased abundance of BCORP1.

Protein Phosphorylation Changes

The course of cancer plus standard of care treatment (CAT+P) led to significant changes in phosphorylation of two intact proteoforms, two fragments, and two aggregates (Figure 1, Supplementary Table S2). Increased phosphorylation occurred in intact proteoforms of the signaling protein YWHAG and the contractile myofilament protein ACTG2. In the CAT+T group, phosphorylation of six intact proteoforms, four fragments, and three aggregates changed significantly during the course of treatment. Altered intact proteoforms included increased phosphorylation of ANKRD2 (4 spots) and decreased phosphorylations in BCORP (2 spots). Phosphorylation changes in these two proteins coincided with similar changes in abundance in these groups.

Protein S-Nitrosylation Changes

In the CAT+P group, a significant change in S-nitrosylation was seen in 26 intact proteoforms, 11 aggregates, and 22 fragments (Figure 1, Supplementary Table S3). The intact proteins included metabolic proteins ALDOA (3 spots), ENO1, CKM (2 spots), AK1, and ACO2, transport protein CA3, transcriptional protein EEF1A2, structural proteins ANKRD2 (3 spots) and LDB3, contractile proteins TNNT3 (2 spots), MYL2, ACTA1, and ACTG2, signaling protein MYOZ1 (2 spots), calcium regulation proteins ANXA2 and BEST3, uncoupling protein SLC25A47, and protein with uncertain function BCORP1 (2 spots). In the CAT+T group, the proteins that were nitrosylated differently after treatment included 10 intact proteoforms, nine fragments, and seven aggregates. The intact proteins included metabolic proteins ALDOA, PYGM, and ENO1, structural proteins ANKRD2 and LDB3, contractile proteins MYL2 (2 spots) and TNNC2, transport protein HBA2 and degradation protein TRIM72. S-Nitrosylation of ALDOA, ENO1, ANKRD2, LDB3, and MYL2 changed significantly in both CAT+P and CAT+T. Overall, 36 proteoforms were upregulated in CAT+P and 23 proteoforms were downregulated over the course of treatment. In CAT+T, 16 proteoforms were downregulated and 11 were upregulated.

Gene Set Enrichment Analysis

Cancer and CAT with or without testosterone supplementation resulted in significant changes in the abundance and post-translational modifications (phosphorylation and S-nitrosylation) of numerous proteins within the skeletal muscle proteome. Because the net effects of individual protein modifications on either the up or down regulation of entire processes and pathways are complex or unknown, the significantly altered proteins (positive and negative fold changes in intact, fragmented, or aggregate proteoforms) were assessed together to identify affected gene sets (GO) and pathways (KEGG). Analyses of the proteomic responses in vastus lateralis from advanced cancer patients receiving CAT with or without adjuvant testosterone treatment revealed several muscle specific gene sets (GO, Figures 2A,B, Supplementary Figure S1) and pathways (KEGG, Figures 2A,C, Supplementary Figure S2) that were affected as demonstrated by changes in protein abundances, phosphorylation, and S-nitrosylation of key skeletal muscle regulators. Gene set (GO) analyses indicated that the course of CAT+P (control group) had protein abundance changes associated with muscle filament sliding, lipoprotein regulation, and cell membrane processes such as raft assembly and vesicle fusion. Muscle contraction was identified as an altered process regardless of treatment with testosterone or placebo. Additionally, treatment with testosterone (CAT+T) redirected cellular processes towards regulation of proteins involved in cell development, morphogenesis, and cell differentiation. Pathway analyses (KEGG) indicated that CAT+P altered abundances of proteins involved in myopathy related pathways while CAT+T affected cell metabolism pathways. Protein phosphorylation changes during CAT+P were largely directed at metabolic processes, oxygen transport, mitochondria apoptotic signaling, and skeletal muscle development. There were no shared biological processes identified between CAT+P and CAT+T based on phosphorylation changes. The effects of testosterone treatment on protein phosphorylation were in line with the changes in protein abundances. Testosterone treatment appeared to affect regulation of cell differentiation, morphogenesis, and assembly. Pathway analyses (KEGG) indicated that CAT+P marginally affected phosphorylation of proteins involved in butanoate and amino acid metabolism. Conversely, CAT+T modified smooth muscle contraction pathway proteins. The biological processes identified by changes in S-nitrosylation were fairly consistent between CAT+P and CAT+T suggesting that cancer and CAT may have been the predominant regulators of S-nitrosylation. Over the course of CAT, regardless of treatment with testosterone, S-nitrosylation changes appeared largely directed towards processes associated with regulation of muscle contraction and energy metabolism. Pathway analyses (KEGG) indicated that CAT+P affected S-nitrosylation of proteins involving numerous pathways related to myopathies and cell metabolism. CAT+T affected the same pathways identified in CAT+P but also affected proteins with roles in butanoate metabolism, adrenergic signaling, and insulin resistance pathways.

FIGURE 2

FIGURE 2. Gene set enrichment analysis (GSEA). (A) Volcano Plots. Volcano plots of protein spots identified by 2D electrophoresis represented as -log10 p-value vs. log2 of fold change from pre to post. Proteins were deemed significantly altered with a p < 0.05 and a log2 fold change >1.5 (increased) or < −1.5 (decreased). Analyses identified significant differences in abundance (red), phosphorylation (blue), and S-nitrosylation (purple) between pre and post treatment biopsies for a subset of proteins (significantly altered proteins are colored). (B) Gene Ontology (GO). Significantly altered proteins were used to identify affected gene sets using the 2021 GO Biological Process 2021 database. The top 10 pathways identified by Enrichr combined score (see text for description) are displayed. Pathways that were identified as significantly altered (adjusted p < 0.05) are denoted by a dark colored bar for abundance (red), phosphorylation (blue), and S-nitrosylation (purple). The number and identity of significantly altered proteins and total number of proteins in the identified pathway are denoted on each bar. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG). Significantly altered proteins were used to identify affected pathways using the 2021 KEGG pathway database. The top 10 pathways identified by Enrichr combined score (see text for description) are displayed. Pathways that were identified as significantly altered (adjusted p < 0.05) are denoted by a dark colored bar for abundance (red), phosphorylation (blue), and S-nitrosylation (purple). The number and identity of significantly altered proteins and total number of proteins in the identified pathway are denoted on each bar.

Discussion

Seven weeks of cancer and standard-of-care cancer treatment was associated with significantly altered skeletal muscle protein abundance, phosphorylation, and S-nitrosylation. Although adjunct testosterone treatment largely preserved skeletal muscle mass in these patients (Wright et al., 2018), proteostasis was disrupted both with and without adjunct testosterone. Supporting these observations, treatment with testosterone likely promoted several biological processes involved in muscle tissue development, including cell differentiation and morphogenesis. The most noticeable differences between placebo and testosterone treated muscle included protein abundance and phosphorylation. None of the same proteoforms changed in abundance when comparing these two groups. Cancer and CAT profoundly affected protein S-nitrosylation regardless of testosterone treatment. The pathways identified as being most altered tended to converge around myopathies, muscle contraction (including smooth and cardiac muscle), and substrate metabolism. The changes in the skeletal muscle proteomes support the phenotypic changes, such as preserved muscle mass in the CAT+T compared to the CAT+P group, that were previously reported in these patients (Wright et al., 2018).

Skeletal muscle atrophy, as seen in cancer cachexia, results from disrupted cellular proteostasis characterized by decreased protein synthesis, increased protein degradation, or some combination of both (Roy and Kumar, 2019). Testosterone and selective androgen receptor modulators (SARMS) have been suggested and studied as interventions against cancer induced wasting (Basaria et al., 2001; Bhasin et al., 2006; Narayanan et al., 2008; Li et al., 2021). Our group previously showed that the placebo treated group lost 3.3% lean body mass during the course of CAT while the testosterone treated group gained 3.2% (Wright et al., 2018). Testosterone treatment offered protection against the loss of lean mass involving a distinct set of genes compared to the CAT+P group. Cancer and associated therapies are associated with both skeletal muscle cachexia and concomitant cardiac atrophy and dysfunction (Kazemi-Bajestani et al., 2014; Murphy, 2016). Several proteomic shifts that occurred during CAT were alterations to pathways with known involvement in cardiomyopathy, and despite cardiac muscle being a distinct cell type from skeletal muscle, likely pointed towards similar myopathic processes in skeletal muscle of these patients (Rausch et al., 2021). Nevertheless, it should be noted that an additional ancillary analysis of this study demonstrated that testosterone treatment was associated with improvement in some measures of cardiac function including left ventricular ejection fraction (Scott et al., 2019).

The observed decrease in contractile proteins MYL3 and MYL1 in the CAT+P group are consistent with cachectic myopathy. However, the structural protein DES increased in CAT+P, which is counterintuitive to the catabolic phenotype of these patients. Expression of the intermediate filament DES has been shown to increase near carcinomatous glands and microvessels but it is difficult to explain the increase in expression within skeletal muscle observed here (Arentz et al., 2011; Mittal, 2020). It is also important to consider the potential contribution of standard of care (SOC) treatment on PRE to POST treatment changes. Muscle wasting can be synergistically driven by factors related to cancer-host interactions as well as toxicities contributed by CAT (Campelj et al., 2021). Chemotherapy can induce increased DES expression in rhabdomyosarcomas (Klunder et al., 2003) and it is expected that the SOC that was offered to the patients altered protein expression and post-translational modifications to skeletal muscle proteoforms over the course of treatment in this investigation.

The stretch-responsive ankyrin-repeat protein ANKRD2 (Kojic et al., 2011) increased along with lean body mass in response to CAT+T. ANKRD2 has been associated with increased hypertrophy (Kemp et al., 2000) which could be driven by a testosterone induced stimulation of the IGF-I signaling system (Urban et al., 1995). ANKRD2 is also a potent repressor of inflammatory responses through direct interaction with the NF-κB repressor subunit p50. ANKRD2 protein and phosphorylation levels are shown to modulate the balance between physiological and pathological inflammatory responses in muscle (Bean et al., 2014). The observation of testosterone induced upregulation of ANKRD2 also adds support to the anti-inflammatory properties of testosterone. There is evidence that the expression of the structural and mechanosensitive proteins DES and ANKRD1 (which is expressed in cardiac and smooth muscle, but not skeletal muscle) are inversely regulated in smooth muscle (Mohamed and Boriek, 2012). The mechanisms between DES and ANKRD2 expression in skeletal muscle are unclear but such an inverse relationship could further support the observed increases in DES in the placebo group vs. increased ANKRD2 in testosterone treated patients.

Protein functions can be mediated through reversible post-translational modifications such as phosphorylation or S-nitrosylation of specific protein sites. S-nitrosylation involves the post-translational modification of proteins through covalent attachment of a nitrogen monoxide group to the thiol side chain of a cysteine (S) amino acid within the peptide. Nitric oxide (NO) necessary for S-nitrosylation is produced by nitric oxide synthase (NOS) (Foster et al., 2009). S-nitrosylation plays an important regulatory role in human physiology and alterations have been implicated in pathologies including cardiovascular, pulmonary, musculoskeletal and neurological (dys)function, as well as in cancer. NO-mediated protein S-nitrosylation has been established as a major signaling mechanism for cancer initiation and development (Tan et al., 2019). NO appears to have varying effects in cancer either being inhibitory or stimulatory depending on the case (Xu et al., 2002; Hickok and Thomas, 2010). NOS activity has been detected in tumor cells, and is associated with tumor grade, proliferation rate, etc. High levels of NOS expression may be cytostatic or cytotoxic for tumor cells, while low levels of NOS activity can have the opposite effect, resulting in tumor growth. Among the three known types of NOS, inducible NOS (iNOS) activity is upregulated in response to inflammatory cytokines (Ramamoorthy et al., 2009) and has been demonstrated in skeletal muscle (Foster et al., 2009; Hall et al., 2011). It is expected that much of the S-nitrosylation observed in the present study is the result of iNOS due to the inflammatory response to cancer (Hall et al., 2011). These results also support an anti-inflammatory role of testosterone treatment as more than double the number of proteoform S-nitrosylation changes occurred in response to CAT+P over the course of the study compared to CAT+T. Increased nitrosative stress is thought to be advantageous for tumor biology and interventions that reduce protein S-nitrosylation are actively being investigated as cancer treatment strategies (Rizza and Filomeni, 2020). Testosterone has been shown to inhibit iNOS synthesis (Friedl et al., 2000; Chavoshan et al., 2012). Additionally, as smooth muscle contraction was identified as one of the possible pathways affected by testosterone treatment it is possible that changes in endothelial NOS (eNOS) activity contributed to group differences in the levels of protein S-nitrosylation given the inevitability of vascular tissue in skeletal muscle homogenate. In the present study, 1.56 times as many proteoforms increased in S-nitrosylation compared to decreasing in CAT+P (36/23) while 1.45 times as many proteoforms decreased in S-nitrosylation compared to increasing in CAT+T. (16/11). The precise mechanism through which testosterone treatment may be regulating S-nitrosylation events deserves further investigation and whether there is therapeutic benefit for testosterone mediated denitrosylation in cancer therapy remains to be determined.

Limitations

The primary objective of this randomized phase II study was to investigate whether testosterone administered during standard of care chemotherapy and/or radiation could help patients with squamous cell carcinoma maintain body weight, muscle size, and muscle strength during treatment. The patients were not classified on a cachectic scale, as described by Fearon et al. (Fearon et al., 2011). The timeframe of patients presenting to the clinic (and our study) was often very short and did not allow for proper classification while under observation, or the quantity of losses were not well known (even by the patient) before seeking treatment. Also, some patients were advanced or recurrent patients and weight loss had already occurred at an earlier stage. The authors also acknowledge the limited sample size and heterogeneity of the study population, including the range of ages, sex, cancer types, and treatments. This ancillary study was not powered to assess the contributions of age on treatment efficacy, although this is an acknowledged factor worth investigating (Horstman et al., 2015). As this was an intent-to-treat study within a population of squamous cell carcinoma patients, individual variations in treatments or daily activities of living were expected to occur. This included treatment interruptions and blood transfusions as directed by oncologists. Other factors besides CAT such as inactivity and poor nutrition could affect changes in protein expression during the observed treatment period. As previously reported, both placebo and testosterone treated patients remained consistently sedentary throughout treatment (vigorous activity below 0.5%) and were in general poor health and nutritional status (Wright et al., 2018). These confounding factors were similarly variable between treatment groups and were considered underlying conditions to CAT in context of the finding reported here.

As the analyses were targeted to skeletal muscle tissue, the authors decided against using murine specific background databases such as the murine based SysMyo Muscle Gene Set available in Enrichr or the Mus Musculus Gene Set available through GO Consortium that were muscle specific. We acknowledge that a consequence of comparing a human skeletal muscle proteome against more comprehensive databases such as KEGG and GO may be a biological bias, resulting in a ‘skeletal muscle signatures’ in the results (Timmons et al., 2015). However, the objective of this report was to compare the proteome in response to treatment without limiting the findings towards only previously identified skeletal muscle pathways. Despite the introduction of these muscle signatures, the results point towards differences in the signatures between treatment groups.

Conclusion

Seven weeks of CAT was associated with changes in the skeletal muscle proteome of patients with cervical or head and neck cancer. Proteostasis was disrupted in patients receiving both adjunct testosterone and placebo, collectively consistent with cachectic myopathy, including loss of structural integrity, decreased contractile function, and disrupted metabolism. The two treatment groups differed in specific proteomic alterations to abundance, phosphorylation, and S-nitrosylation resulting in differences in altered pathways. In conjunction with prior results showing that adjunct testosterone treatment preserved skeletal muscle mass of patients, the differential proteomic changes in these two groups suggests that testosterone at least partially directed the functional pathways towards muscle tissue development and minimized the negative structural and contractile impact of cancer and associated therapies. The shotgun proteomic analysis approach revealed how CAT with or without testosterone treatment differentially affected protein expression and post-translational modification within skeletal muscle of these cancer patients. Targeted studies are necessary to determine the specific interactions and mechanisms that contribute to alterations in the proteome and how they lead to muscle wasting and other phenotypic changes associated with cancer cachexia.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board, University of Texas Medical Branch. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Concept and Design: ELD, RJU, and MS-M. Study Execution: ELD, KMR, CPD, CRG, RJU, and MS-M. Sample Analyses: JEW and KVS. Data Analyses: ELD, TJW, ARF, ABP, and KVS. Writing and Editing: ELD, TJW, and MS-M. All authors contributed to manuscript revision, read, and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.879263/full#supplementary-material

Supplementary Figure S1 | GO Biological Prosesses. Proteins found to be altered in either abundance, phosphorylation, or S-nitrosylation were combined to identify affected gene sets using the 2021 GO Biological Process database.

Supplementary Figure S2 | KEGG Pathways. Proteins found to be altered in either abundance, phosphorylation, or S-nitrosylation were combined to identify affected pathways using the 2021 KEGG pathway database.

Supplementary Table S1 | Protein abundances (A) CAT+P (B) CAT+T.

Supplementary Table S2 | Protein Phosphorylation (A) CAT+P (B) CAT+T.

Supplementary Table S3 | Protein S-Nitrosylation (A) CAT+P (B) CAT+T.

References

Ali S., Garcia J. M. (2014). Sarcopenia, Cachexia and Aging: Diagnosis, Mechanisms and Therapeutic Options - a Mini-Review. Gerontology 60 (4), 294–305. doi:10.1159/000356760

PubMed Abstract | CrossRef Full Text | Google Scholar

Arentz G., Chataway T., Price T. J., Izwan Z., Hardi G., Cummins A. G., et al. (2011). Desmin Expression in Colorectal Cancer Stroma Correlates with Advanced Stage Disease and Marks Angiogenic Microvessels. Clin. Proteom 8 (1), 16. doi:10.1186/1559-0275-8-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Basaria S., Wahlstrom J. T., Dobs A. S. (2001). Clinical Review 138: Anabolic-Androgenic Steroid Therapy in the Treatment of Chronic Diseases. J. Clin. Endocrinol. Metab. 86 (11), 5108–5117. doi:10.1210/jcem.86.11.7983

PubMed Abstract | CrossRef Full Text | Google Scholar

Bean C., Verma N. K., Yamamoto D. L., Chemello F., Cenni V., Filomena M. C., et al. (2014). Ankrd2 Is a Modulator of NF-κB-Mediated Inflammatory Responses during Muscle Differentiation. Cell Death Dis. 5, e1002. doi:10.1038/cddis.2013.525

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhasin S., Calof O. M., Storer T. W., Lee M. L., Mazer N. A., Jasuja R., et al. (2006). Drug Insight: Testosterone and Selective Androgen Receptor Modulators as Anabolic Therapies for Chronic Illness and Aging. Nat. Rev. Endocrinol. 2 (3), 146–159. doi:10.1038/ncpendmet0120

CrossRef Full Text | Google Scholar

Campelj D. G., Goodman C. A., Rybalka E. (2021). Chemotherapy-Induced Myopathy: The Dark Side of the Cachexia Sphere. Cancers 13 (14), 3615. doi:10.3390/cancers13143615

PubMed Abstract | CrossRef Full Text | Google Scholar

Chavoshan B., Fournier M., Lewis M. I., Porszasz J., Storer T. W., Da X., et al. (2012). Testosterone and Resistance Training Effects on Muscle Nitric Oxide Synthase Isoforms in COPD Men. Respir. Med. 106 (2), 269–275. doi:10.1016/j.rmed.2011.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen E. Y., Tan C. M., Kou Y., Duan Q., Wang Z., Meirelles G. V., et al. (2013). Enrichr: Interactive and Collaborative HTML5 Gene List Enrichment Analysis Tool. BMC Bioinforma. 14, 128. doi:10.1186/1471-2105-14-128

PubMed Abstract | CrossRef Full Text | Google Scholar

Damrauer J. S., Stadler M. E., Acharyya S., Baldwin A. S., Couch M. E., Guttridge D. C. (2018). Chemotherapy-induced Muscle Wasting: Association with NF-κB and Cancer Cachexia. Eur. J. Transl. Myol. 28 (2), 7590. doi:10.4081/ejtm.2018.7590

PubMed Abstract | CrossRef Full Text | Google Scholar

Dillon E. L., Soman K. V., Wiktorowicz J. E., Sur R., Jupiter D., Danesi C. P., et al. (2019). Proteomic Investigation of Human Skeletal Muscle before and after 70 Days of Head Down Bed Rest with or without Exercise and Testosterone Countermeasures. PloS one 14 (6), e0217690. doi:10.1371/journal.pone.0217690

PubMed Abstract | CrossRef Full Text | Google Scholar

Dos Santos M. R., Storer T. W. (2022). Testosterone Treatment as a Function-Promoting Therapy in Sarcopenia Associated with Aging and Chronic Disease. Endocrinol. Metab. Clin. North Am. 51 (1), 187–204. doi:10.1016/j.ecl.2021.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Fearon K., Strasser F., Anker S. D., Bosaeus I., Bruera E., Fainsinger R. L., et al. (2011). Definition and Classification of Cancer Cachexia: an International Consensus. Lancet Oncol. 12, 489–495. doi:10.1016/s1470-2045(10)70218-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Foster M. W., Hess D. T., Stamler J. S. (2009). Protein S-Nitrosylation in Health and Disease: a Current Perspective. Trends Mol. Med. 15 (9), 391–404. doi:10.1016/j.molmed.2009.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedl R., Brunner M., Moeslinger T., Spieckermann P. G. (2000). Testosterone Inhibits Expression of Inducible Nitric Oxide Synthase in Murine Macrophages. Life Sci. 68 (4), 417–429. doi:10.1016/s0024-3205(00)00953-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall D. T., Ma J. F., Di Marco S., Gallouzi I.-E. (2011). Inducible Nitric Oxide Synthase (iNOS) in Muscle Wasting Syndrome, Sarcopenia, and Cachexia. Aging 3 (8), 702–715. doi:10.18632/aging.100358

PubMed Abstract | CrossRef Full Text | Google Scholar

Hickok J., Thomas D. (2010). Nitric Oxide and Cancer Therapy: the Emperor Has NO Clothes. Cpd 16 (4), 381–391. doi:10.2174/138161210790232149

PubMed Abstract | CrossRef Full Text | Google Scholar

Horstman A. M., Dillon E. L., Urban R. J., Sheffield-Moore M. (2012). The Role of Androgens and Estrogens on Healthy Aging and Longevity. Journals Gerontology Ser. A Biol. Sci. Med. Sci. 67 (11), 1140–1152. doi:10.1093/gerona/gls068

CrossRef Full Text | Google Scholar

Hunt L. C., Graca F. A., Pagala V., Wang Y.-D., Li Y., Yuan Z.-F., et al. (2021). Integrated Genomic and Proteomic Analyses Identify Stimulus-dependent Molecular Changes Associated with Distinct Modes of Skeletal Muscle Atrophy. Cell Rep. 37 (6), 109971. doi:10.1016/j.celrep.2021.109971

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasvis P., Kilgour R. D. (2021). Diet and Exercise Interventions in Patients with Pancreatic Cancer. Pancreas 50 (5), 657–666. doi:10.1097/mpa.0000000000001829

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazemi-Bajestani S. M. R., Becher H., Fassbender K., Chu Q., Baracos V. E. (2014). Concurrent Evolution of Cancer Cachexia and Heart Failure: Bilateral Effects Exist. J. Cachexia Sarcopenia Muscle 5 (2), 95–104. doi:10.1007/s13539-014-0137-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Kemp T. J., Sadusky T. J., Saltisi F., Carey N., Moss J., Yang S. Y., et al. (2000). Identification of Ankrd2, a Novel Skeletal Muscle Gene Coding for a Stretch-Responsive Ankyrin-Repeat Protein. Genomics 66 (3), 229–241. doi:10.1006/geno.2000.6213

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim A. J., Hong D. S., George G. C. (2021). Diet-related Interventions for Cancer-Associated Cachexia. J. Cancer Res. Clin. Oncol. 147 (5), 1443–1450. doi:10.1007/s00432-021-03592-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Klunder J. W., Komdeur R., van der Graaf W. T. A., de Bont E. J. S. M., Hoekstra H. J., van den Berg E., et al. (2003). Expression of Multidrug Resistance-Associated Proteins in Rhabdomyosarcomas before and after Chemotherapy: the Relationship between Lung Resistance-Related Protein (LRP) and Differentiation. Hum. Pathol. 34 (2), 150–155. doi:10.1053/hupa.2003.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Kojic S., Radojkovic D., Faulkner G. (2011). Muscle Ankyrin Repeat Proteins: Their Role in Striated Muscle Function in Health and Disease. Crit. Rev. Clin. Lab. Sci. 48 (5-6), 269–294. doi:10.3109/10408363.2011.643857

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuleshov M. V., Jones M. R., Rouillard A. D., Fernandez N. F., Duan Q., Wang Z., et al. (2016). Enrichr: a Comprehensive Gene Set Enrichment Analysis Web Server 2016 Update. Nucleic Acids Res. 44 (W1), W90–W97. doi:10.1093/nar/gkw377

PubMed Abstract | CrossRef Full Text | Google Scholar

Li Y., Jin H., Chen Y., Huang T., Mi Y., Zou Z. (2021). Cancer Cachexia: Molecular Mechanism and Pharmacological Management. Biochem. J. 478 (9), 1663–1688. doi:10.1042/bcj20201009

PubMed Abstract | CrossRef Full Text | Google Scholar

Meij B. S., Teleni L., McCarthy A. L., Isenring E. A. (2021). Cancer Cachexia: an Overview of Diagnostic Criteria and Therapeutic Approaches for the Accredited Practicing Dietitian. J. Hum. Nutr. Diet. 34 (1), 243–254. doi:10.1111/jhn.12811

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal B. (2020). Desmin Dysregulation in Gall Bladder Cancer. Indian J. Med. Res. 151, 273–274. doi:10.4103/ijmr.ijmr_1540_19

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohamed J. S., Boriek A. M. (2012). Loss of Desmin Triggers Mechanosensitivity and Up‐regulation of Ankrd1 Expression through Akt‐NF‐κB Signaling Pathway in Smooth Muscle Cells. FASEB J. 26 (2), 757–765. doi:10.1096/fj.10-160291

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy K. T. (2016). The Pathogenesis and Treatment of Cardiac Atrophy in Cancer Cachexia. Am. J. Physiology-Heart Circulatory PhysiologyHeart circulatory physiology 310 (4), H466–H477. doi:10.1152/ajpheart.00720.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Narayanan R., Coss C. C., Yepuru M., Kearbey J. D., Miller D. D., Dalton J. T. (2008). Steroidal Androgens and Nonsteroidal, Tissue-Selective Androgen Receptor Modulator, S-22, Regulate Androgen Receptor Function through Distinct Genomic and Nongenomic Signaling Pathways. Mol. Endocrinol. 22 (11), 2448–2465. doi:10.1210/me.2008-0160

PubMed Abstract | CrossRef Full Text | Google Scholar

Peixoto da Silva S., Santos J. M. O., Costa e Silva M. P., Gil da Costa R. M., Medeiros R. (2020). Cancer Cachexia and its Pathophysiology: Links with Sarcopenia, Anorexia and Asthenia. J. cachexia, sarcopenia muscle 11 (3), 619–635. doi:10.1002/jcsm.12528

PubMed Abstract | CrossRef Full Text | Google Scholar

Pin F., Barreto R., Couch M. E., Bonetto A., O'Connell T. M. (2019). Cachexia Induced by Cancer and Chemotherapy Yield Distinct Perturbations to Energy Metabolism. J. cachexia, sarcopenia muscle 10 (1), 140–154. doi:10.1002/jcsm.12360

CrossRef Full Text | Google Scholar

Ramamoorthy S., Donohue M., Buck M. (2009). Decreased Jun-D and Myogenin Expression in Muscle Wasting of Human Cachexia. Am. J. Physiology-Endocrinology MetabolismEndocrinology metabolism 297 (2), E392–E401. doi:10.1152/ajpendo.90529.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Rausch V., Sala V., Penna F., Porporato P. E., Ghigo A. (2021). Understanding the Common Mechanisms of Heart and Skeletal Muscle Wasting in Cancer Cachexia. Oncogenesis 10 (1), 1. doi:10.1038/s41389-020-00288-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizza S., Filomeni G. (2020). Exploiting S-Nitrosylation for Cancer Therapy: Facts and Perspectives. Biochem. J. 477 (19), 3649–3672. doi:10.1042/bcj20200064

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy A., Kumar A. (2019). ER Stress and Unfolded Protein Response in Cancer Cachexia. Cancers 11 (12), 1929. doi:10.3390/cancers11121929

PubMed Abstract | CrossRef Full Text | Google Scholar

Schütte K., Middelberg-Bisping K., Schulz C. (2020). ‘Nutrition and Gastroenterological Support in End of Life Care’, Best Practice & Research. Clin. Gastroenterol. 48-49, 101692.

Google Scholar

Scott J. M., Dillon E. L., Kinsky M., Chamberlain A., McCammon S., Jupiter D., et al. (2019). Effects of Adjunct Testosterone on Cardiac Morphology and Function in Advanced Cancers: an Ancillary Analysis of a Randomized Controlled Trial. BMC cancer 19 (1), 778. doi:10.1186/s12885-019-6006-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheffield-Moore M., Paddon-Jones D., Casperson S. L., Gilkison C., Volpi E., Wolf S. E., et al. (2006). Androgen Therapy Induces Muscle Protein Anabolism in Older Women. J. Clin. Endocrinol. metabolism 91 (10), 3844–3849. doi:10.1210/jc.2006-0588

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheffield-Moore M., Wiktorowicz J. E., Soman K. V., Danesi C. P., Kinsky M. P., Dillon E. L., et al. (2013). Sildenafil Increases Muscle Protein Synthesis and Reduces Muscle Fatigue. Clin. Transl. Sci. 6 (6), 463–468. doi:10.1111/cts.12121

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan C., Li Y., Huang X., Wei M., Huang Y., Tang Z., et al. (2019). Extensive Protein S-Nitrosylation Associated with Human Pancreatic Ductal Adenocarcinoma Pathogenesis. Cell Death Dis. 10 (12), 914. doi:10.1038/s41419-019-2144-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Urban R. J., Bodenburg Y. H., Gilkison C., Foxworth J., Coggan A. R., Wolfe R. R., et al. (1995). Testosterone Administration to Elderly Men Increases Skeletal Muscle Strength and Protein Synthesis. Am. J. Physiol. 269 (5 Pt 1), E820–E826. doi:10.1152/ajpendo.1995.269.5.E820

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiktorowicz J. E., Stafford S. J., Garg N. J. (2017). Protein Cysteinyl-S-Nitrosylation. Methods Enzym. 586, 1–14. doi:10.1016/bs.mie.2016.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiktorowicz J. E., Stafford S., Rea H., Urvil P., Soman K., Kurosky A., et al. (2011). Quantification of Cysteinyl S-Nitrosylation by Fluorescence in Unbiased Proteomic Studies. Biochemistry 50 (25), 5601–5614. doi:10.1021/bi200008b

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright T. J., Dillon E. L., Durham W. J., Chamberlain A., Randolph K. M., Danesi C., et al. (2018). A Randomized Trial of Adjunct Testosterone for Cancer-Related Muscle Loss in Men and Women. J. cachexia, sarcopenia muscle 9 (3), 482–496. doi:10.1002/jcsm.12295

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie Z., Bailey A., Kuleshov M. V., Clarke D. J. B., Evangelista J. E., Jenkins S. L., et al. (2021). Gene Set Knowledge Discovery with Enrichr. Curr. Protoc. 1 (3), e90. doi:10.1002/cpz1.90

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu W., Liu L. Z., Loizidou M., Ahmed M., Charles I. G. (2002). The Role of Nitric Oxide in Cancer. Cell Res. 12 (5-6), 311–320. doi:10.1038/sj.cr.7290133

PubMed Abstract | CrossRef Full Text | Google Scholar