Brandon M. Roberts1
Sarah E. Deemer2
Daniel L. Smith Jr.3
James A. Mobley4
Nicolas Musi5,6
Eric P. Plaisance1*- 1Department of Human Studies, Division of Molecular and Translational Biomedicine, University of Alabama at Birmingham, Birmingham, AL, United States
- 2Department of Kinesiology, Health Promotion, and Recreation, University of North Texas, Denton, TX, United States
- 3Department of Nutrition Sciences, Division of Molecular and Translational Biomedicine, University of Alabama at Birmingham, Birmingham, AL, United States
- 4Department of Anesthesiology and Perioperative Medicine, Division of Molecular and Translational Biomedicine, University of Alabama at Birmingham, Birmingham, AL, United States
- 5Barshop Institute for Longevity and Aging Studies, University of Texas Health Science Center San Antonio, San Antonio, TX, United States
- 6San Antonio Geriatric Research, Education, and Clinical Center, San Antonio, TX, United States
Exogenous ketone ester supplementation provides a means to increase circulating ketone concentrations without the dietary challenges imposed by ketogenic diets. Our group has shown that oral R,S-1,3, butanediol diacetoacetate (BD-AcAc2) consumption results in body weight loss or maintenance with moderate increases in circulating ketones. We have previously shown a diet consisting of 25% BD-AcAc2 can maintain lean body mass (LBM) and induce fat mass (FM) loss in young, healthy male mice, but the underlying mechanisms are still unknown. Therefore, the purpose of this study was to determine if a diet consisting of 25% BD-AcAc2 (ketone ester, KE) would alter body composition, transcriptional regulation, the proteome, and the lipidome of skeletal muscle in aged mice. We hypothesized that the KE group would remain weight stable with improvements in body composition compared to controls, resulting in a healthy aging phenotype. Male C57BL/6J mice (n = 16) were purchased from Jackson Laboratories at 72 weeks of age. After 1 week of acclimation, mice were weighed and randomly assigned to one of two groups (n = 8 per group): control (CON) or KE. A significant group by time interaction was observed for body weight (P < 0.001), with KE fed mice weighing significantly less than CON. FM increased over time in the control group but was unchanged in the KE group. Furthermore, LBM was not different between CON and KE mice despite KE mice weighing less than CON mice. Transcriptional analysis of skeletal muscle identified 6 genes that were significantly higher and 21 genes that were significantly lower in the KE group compared to CON. Lipidomic analysis of skeletal muscle identified no differences between groups for any lipid species, except for fatty acyl chains in triacylglycerol which was 46% lower in the KE group. Proteomics analysis identified 44 proteins that were different between groups, of which 11 were lower and 33 were higher in the KE group compared to CON. In conclusion, 72-week-old male mice consuming the exogenous KE, BD-AcAc2, had lower age-related gains in body weight and FM compared to CON mice. Furthermore, transcriptional and proteomics data suggest a signature in skeletal muscle of KE-treated mice consistent with markers of improved skeletal muscle regeneration, improved electron transport chain utilization, and increased insulin sensitivity.
Introduction
Ketogenic diets (KDs) improve or maintain metabolic function and attenuate age-related increases in adiposity and bodyweight, presumably by increasing circulating ketone concentrations in both animal models and humans (1, 2). Clinical studies have demonstrated the beneficial effect of KDs as an efficient nutritional strategy in the treatment of obesity (3–5). This effect occurs not only in modifying body weight and body composition but also in modulating epigenetic markers, as well as circulating levels of myokines and circulating levels of cytokines and markers of oxidative stress and this effect appears to be induced by the action of ketone bodies synergistically with the weight loss induced by the VLCKD treatment (6–8). In mice, KDs have been shown to extend longevity and healthspan by decreasing oxidative and endoplasmic reticulum (ER) stress resulting in lower protein turnover in skeletal muscle, which may allow greater maintenance of muscle mass and function with both age and disease (9, 10). These findings have practical implications in humans, but by nature of the requirement for extensive carbohydrate restriction and high fat content (often > 75%), KDs present adherence challenges and may not be a long-term strategy for all humans.
Exogenous ketone ester (KE) supplementation provides a means to increase circulating ketone concentrations without the dietary challenges imposed by KDs. Ketone esters have been used in many forms such as R,S-1,3-butanediol diacetoacetate (BD-AcAc2, ketone diester), D-β-hydroxybutyrate-(R)-1,3 butanediol (ketone monoester), and more recently Bis-hexanoyl (R)-1,3-butanediol (BH-BD) (11–13). Studies in animal models show that oral BD-AcAc2 consumption results in body weight loss or maintenance with moderate increases in circulating ketones (0.5–1.0 mM) (14–16). Recently, we examined concentration-dependent effects of BD-AcAc2 on body weight, adiposity, energy intake, and energy expenditure in lean mice showing that on an ad libitum basis, mice consuming a 25% (by kcals) KE diet consumed the same amount of food as an ad libitum fed control group, but had a significant difference in body weight (BW) and fat mass (FM), and maintained lean body mass (LBM) (17). Furthermore, after adjustment for LBM and FM, there was no difference in resting energy expenditure (REE) compared to control (17). Our group has also shown that BD-AcAc2 induces loss of BW and FM in diet-induced obesity without changes in LBM or changes in skeletal muscle thermogenic activity (18). Furthermore, when housed in thermoneutral conditions, mice consuming BD-AcAc2 decreased body weight resulting in lower adiposity. The decrease in body weight observed in KE-fed mice transpired without an increase in REE or TEE (19). These findings suggest that a diet consisting of 25% BD-AcAc2 can maintain LBM and induce FM loss in young, healthy male mice, but the underlying mechanisms are still unknown.
Despite extensive evidence of reductions or maintenance of body weight and adiposity, little is known regarding the effects of BD-AcAc2 on skeletal muscle. Previously, we found that a KE diet does not produce transcriptional changes consistent with mitochondrial biogenesis or respiration (18). Therefore, the purpose of this study was to determine if a diet consisting of 25% BD-AcAc2 would alter body composition, transcriptional regulation, the proteome, and the lipidome of skeletal muscle in aged B6 mice. We hypothesized that the KE group would remain weight stable with improvements in body composition (higher LBM and lower FM) compared to controls resulting in a healthy aging phenotype.
Materials and methods
Study design and diets
Male C57BL/6J mice (n = 16) were purchased from Jackson Laboratories (Bar Harbor, ME, USA) at 72 weeks of age and acclimated for 7 days. At 73 weeks of age, mice were weighed and randomly assigned control (CON, n = 8) or KE (KE, n = 8) diet. Mice were fed ad libitum for 9 weeks. The KE diet (Dyets Inc., Bethlehem, PA, #104425) contained 25% kcals from BD-AcAc2, which replaced 25% of carbohydrate energy in the CON diet (Dyets Inc., #104419) with fat and protein content equal between the diets. Each diet was 3.7 kcal/g and details of the individual diets (CON and KE) are described in Deemer et al. (17). The UAB Institutional Animal Care and Use Committee (IACUC) approved the investigation.
Husbandry
Mice were single-housed and maintained on a standard 12:12 light-dark cycle. All animals had access to ad libitum water and food throughout the study. Food intake and body weight were recorded every day, and fresh food was given each day at the same time. Two body composition (fat and lean body mass) measurements were completed in the University of Alabama at Birmingham (UAB) Small Animal Phenotyping Core on day 33 and day 66 of the study by quantitative magnetic resonance (QMR; EchoMedical MRI, Houston, TX, USA). Due to COVID-19 shutdowns body composition was not obtained at baseline.
RNA isolation and nanostring analysis
Tissue samples were flash frozen in liquid nitrogen at the time of dissection and stored at −80°C. RNA from vastus lateralis was isolated according to manufacturer instructions (RNEasy Mini Kit; Qiagen; Germany). All procedures occurred at room temperature. Briefly, an approximately 30 mg sample was homogenized in Buffer RLT with metal beads using a bead mill (FisherbrandTM Bead Mill Homogenizer) for 2 min at 5 m/s. Lysate was centrifuged for 3 min and the supernatant was transferred to a new microcentrifuge tube. A 1:1 by volume of 70% ethanol was added to the supernatant and mixed by trituration. Approximately 700 μL of sample was transferred to the Qiagen spin column and isolation of mRNA was completed per manufacturer instructions. The A260/A280 ratio was measured to quantify mRNA purity (Thermo ScientificTM NanoDropTM Lite Spectrophotometer) and isolated mRNA was concentrated to 15 ng/μL prior to NanoString analysis. Concentrated mRNA was analyzed using a 770-plex Mouse PanCancer Pathways kit (nCounter; NanoString, Seattle, WA, United States) and raw image counts (RCC files) were obtained using a SPRINT Profiler (NanoString) following a 16-h hybridization assay at 65°C.
Data was analyzed by ROSALIND®,1 with a HyperScale architecture developed by ROSALIND, Inc., (San Diego, CA, USA). Read Distribution percentages, violin plots, identity heatmaps, and sample MDS plots were generated as part of the QC step. Normalization, fold changes, and p-values were calculated using criteria provided by Nanostring. Housekeeping probes to be used for normalization are selected based on the geNorm algorithm as implemented in the NormqPCR R library (20). Fold changes and p-values are calculated using the fast method as described in the nCounter® Advanced Analysis 2.0 User Manual. P-value adjustment is performed using the Benjamini-Hochberg method of estimating false discovery rates (FDR).
Proteomics
Proteomics analysis was carried out as previously referenced with minor differences (21). Proteomics samples were compared in the control (n = 7) and KE groups (n = 6). All protein extracts were obtained using T-PER™ Mammalian Protein Extraction Reagent (Thermo Fisher Scientific, Cat# 78510) containing Halt Protease Inhibitor Cocktail (Thermo Fisher Scientific, Cat# 78429) and Dounce-homogenized in lysis buffer, 20–30 strokes per sample, and centrifuged (∼12 Kg) for 10 min at 4oC to remove debris. Protein extracts were then quantified using a Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, Cat# PI23225). Forty (40) μg of protein was diluted in 35 μL of NuPAGE LDS sample buffer (1x final conc., Invitrogen, Cat# NP0007) for each sample. Proteins were then reduced with dithiothreitol (DTT) and denatured at 70°C for 10 min prior to loading onto Novex NuPAGE 10% Bis-Tris protein gels (Invitrogen, Cat.# NP0315BOX) and separated (35 min at 200 V). Gels were stained overnight with Novex Colloidal Blue Staining kit (Invitrogen, Cat# LC6025). Following de-staining, each lane was cut into 6-molecular weight (MW) fractions and equilibrated in 100 mM ammonium bicarbonate (AmBc). Each gel plug was then digested overnight with Trypsin Gold, Mass Spectrometry Grade (Promega, Cat# V5280) following the manufacturer’s instructions. Peptide extracts were reconstituted in 0.1% formic acid/ddH2O at 0.1 μg/μL. Peptide digests (8 μL each) were injected onto a 1,260 Infinity nHPLC stack (Agilent Technologies), and separated using a 75 micron I.D. x 15 cm pulled tip C-18 column (Jupiter C-18 300 Å, 5 micron, Phenomenex).
The XCalibur RAW files were collected in profile mode, centroided and converted to MzXML using ReAdW v. 3.5.1. Mgf files were then created using MzXML2Search (included in TPP v. 3.5) for all scans. The data was searched using SEQUEST (Thermo Fisher Scientific), which was set for three maximum missed cleavages, a precursor mass window of 20 ppm, trypsin digestion, variable modification C at 57.0293, and M @ 15.9949 as a base setting. Searches were performed with redundant sequences removed on the mus musculus specific subset of the UniProt100 database. The list of peptide IDs generated based on SEQUEST search results were filtered using Scaffold (Protein Sciences, Portland, OR, USA). The cut-off values included a minimum peptide length of > 5 AA’s, with no MH + 1 charge states, with peptide probabilities of > 80% C.I., and with the number of peptides per protein ≥ 2. The protein probabilities were set to a > 99.0% C.I. and an FDR < 1.0. Scaffold incorporates the two most common methods for statistical validation of large proteome datasets, the false discovery rate (FDR) and protein probability (22–24). Relative quantification across experiments was then performed via spectral counting (25, 26). Spectral count abundances were then normalized between samples (27).
For protein abundance ratios determined, we set a 1.5-fold change as the threshold for significance, determined empirically by analyzing the inner-quartile data from the control experiments using ln-ln plots, where the Pierson’s correlation coefficient (R) is 0.98, and > 95–99% of the normalized intensities fell between the set fold changes. In each case, both tests (t-test and fold change) were required to pass in order to be considered significant. All multivariate analyses, including 2D HCA HeatMaps and PCA plots were carried out using Qlucore Omics Explorer (Qlucore, Lund Sweden). Gene ontology assignments and pathway analysis were carried out using MetaCore (GeneGO Inc., St. Joseph, MI, USA). Interactions identified within MetaCore are manually correlated using full text articles. Detailed algorithms have been described previously (28, 29).
Lipidomics
Lipid species were analyzed using multidimensional mass spectrometry-based shotgun lipidomic analysis (30). In brief, muscle homogenates containing 0.6 mg of protein (quantified by Pierce BCA assay) was transferred to a disposable glass test tube. A premixture of lipid internal standards (IS) was added prior to conducting lipid extraction for quantification of the targeted lipid species. Lipid extraction was performed using a modified Bligh and Dyer procedure (31), and each lipid extract was reconstituted in chloroform:methanol (1:1, v:v) at a volume of 400 μL/mg protein.
For shotgun lipidomics, lipid extract was further diluted to a final concentration of ∼500 fmol total lipids per μL. Mass spectrometric analysis was performed on a triple quadrupole mass spectrometer (TSQ Altis, Thermo Fisher Scientific, San Jose, CA, USA) and a Q Exactive mass spectrometer (Thermo Scientific, San Jose, CA, USA), both of which were equipped with an automated nanospray device (TriVersa NanoMate, Advion Bioscience Ltd., Ithaca, NY, USA) as described (32). Identification and quantification of lipid species were performed using an automated software program (33, 34). Data processing (e.g., ion peak selection, baseline correction, data transfer, peak intensity comparison and quantitation) was performed as described (34). The result was normalized to protein content (nmol lipid/mg protein or pmol lipid/mg protein).
Statistical analysis of body weight, food intake, and body composition
Graphical and formal statistical tests performed using PROC UNIVARIATE revealed that all variables were normally distributed. Differences within and between groups were analyzed using a GROUP x TIME repeated measures analysis of variance (ANOVA) in SAS (version 9.4m7, Cary, NC, USA). For body weight and food intake, weekly mean values of mice within a group were compared over the 9 weeks. Body composition variables were compared between day 33 and day 66 between and within groups. Due to COVID-19 restrictions, a baseline body composition measurement was not performed. Tissue weights between groups were analyzed by one-way ANOVA. Tukey-Kramer post hoc testing was used to explore significant differences determined by RMANOVA. Significance was set a priori at P < 0.05. Data are expressed as means ± standard error of the mean (SEM).
Results
Body weight, muscle weight, body composition, and food intake
Groups were not significantly different in body weight at baseline. There was a significant group by time interaction in the overall analysis for body weight (P < 0.001), with KE fed mice weighing significantly less than CON at endpoint. However, these differences were not significant after post hoc adjustments for multiple testing. Interestingly, there was a within-group increase in body weight in the CON group in week 2 compared to week 8 (P < 0.01), week 9 (P < 0.001), and week 10 [(P < 0.001) Figure 1A]. There were no between or within group changes in energy intake (kcal/day) (Figure 1B). For body composition, there was an overall group by time interaction (P < 0.01), which was driven by a within-group increase in FM in the CON group (P < 0.01) while there was no significant within-group change in the KE group (Figures 1C,D). There were also no significant differences in LBM, muscle wet weight of vastus lateralis, gastrocnemius, soleus, and tibialis anterior between groups (data not shown).
Figure 1. (A) Weekly mean body weight of ketone ester (KE) mice was significantly less than control (CON) mice only during the last 3 weeks of the dietary intervention. (B) Weekly energy intake was similar between CON and KE groups. (C) Fat mass measured by quantitative magnetic resonance (QMR) was greater in CON compared to KE mice, and increased from day 33 to day 66 in CON mice only. (D) Lean mass measured by QMR was not different between CON and KE mice despite KE mice weighing less than CON mice.
Gene expression
A broad approach was employed to compare gene expression in vastus lateralis between groups using the ROSALIND gene array. This analysis was completed at the endpoint (week 83) of the study. Of the 773 genes analyzed, 6 were significantly higher in the KE group compared CON (Table 1). 23 genes were also significantly lower in the KE group compared to CON (Table 2).
Table 1. Genes that were significantly higher in ketone ester (KE) compared to control (CON).
Table 2. Genes that were significantly lower in ketone ester (KE) compared to control (CON).
Lipidomics
Lipidomics was performed on vastus lateralis in CON (n = 7) and KE groups (n = 6). There were no differences between groups for any lipid species, except for fatty acyl chains in triacylglycerol (nmol/mg protein) which was 46% lower in the KE group compared to CON (P = 0.05, Table 3).
Table 3. Skeletal muscle lipidomic analysis.
Proteomics
There were 1,020 protein IDs with > 99% C.I. and < 1% FDR from the initial analysis. Proteomics was performed on the gastrocnemius in CON (n = 7) and KE groups (n = 6). Many proteins identified were in low enough abundance and presented with a fair amount of heterogeneity “zero” values across each group and were therefore eliminated from the analysis. Of those proteins that were more homogeneous, 545 high confidence hits were observed in at least 3 or more specimens per group with a signal over zero. A non-parametric statistical analysis was then applied to these remaining proteins along with a fold change of > 1.5 or < −1.5, which led to 44 proteins that were statistically changed in abundance in the KE compared to CON. Of those, 11 were significantly higher in KE compared to CON (Table 4) and 33 were significantly lower (Table 5) in KE compared to CON.
Table 4. Proteomic differences in mouse muscle tissues: Ketone ester (KE) vs. control (CON) 11 proteins were significantly higher in KE compared to CON.
Table 5. Proteomic differences in mouse muscle tissues: Ketone ester (KE) vs. control diets 33 proteins were significantly lower in KE compared to CON.
Discussion
The purpose of this study was to determine the effects of a 25% KE diet on body weight, body composition, and the skeletal muscle transcriptional, proteomic, and lipidomic responses in aging 72-week-old C57BL/6J mice for nine weeks. We found mice consuming a KE diet remained weight stable, with no changes in FM or LBM, yet the control group increased body weight and fat mass over the course of the study. There were no differences in energy intake between the groups, which agrees with our previous findings in juvenile mice (19). Interestingly, the KE group displayed a unique skeletal muscle transcriptional and proteomic profile compared to the control group yet had a very similar lipidomic profile.
Transcriptional analysis indicated six genes were higher in the KE group compared to CON. Growth Arrest And DNA Damage Inducible Alpha (Gadd45a) expression was 1.85-fold higher in KE mice compared to CON. Previous longitudinal data indicates that Gadd45a is elevated in aged mouse skeletal muscle during the transition to sarcopenia, which involves progressive muscle atrophy (35). Others have reported that Gadd45a expression represents a protective negative feedback response to denervation, which delays the rate of atrophy and myofiber type transition, potentially preserving myofiber type during chronic denervation (36). We found that tenascin-C (Tnc) was upregulated 1.44-fold in the KE group, and previous findings have implicated Tnc in the formation, maturation, and stabilization of the neuromuscular junction (37). In line with these findings, RUNX Family Transcription Factor 1 (Runx1) was up-regulated 1.69-fold in the KE group. Runx1 is required to sustain skeletal muscle by preventing denervated myofibers from undergoing myofibrillar disorganization and autophagy (38). Runx1 ablation results in excessive autophagy during denervation which leads to severe atrophy, suggesting that these findings are protective with regards to age-related autophagy and atrophy of skeletal muscle (39). Other studies have found similar results in aged mice, showing Gadd45a and Runx1 are elevated with a KD, which has been shown to mitigate sarcopenia (10). Another transcript, calponin homology-associated smooth muscle protein (Chasm), which is required for tropomyosin binding, was 1.79-fold lower in KE compared to CON (40). Some have suggested the protein form of CHASM is a discrete marker of Type IIa muscle fibers, which supports the idea that the KE group may be undergoing a transition from Type II to Type I fibers to attenuate muscle denervation induced with aging (41, 42). Decreased calcineurin A (CnA) in our proteomic analysis (−2.82-fold compared to CON) also suggests muscle protection since CnA has been shown to indirectly protect muscle fibers from atrophy by raising the proportion of slow fibers in muscles (43). This occurs, in part, by maintaining peroxisome proliferator-activated receptor-gamma coactivator (PGC-1α) through activation of the myocyte enhancer factor-2 (MEF2) and Nuclear factor of activated T-cells (NFAT) transcription factors (43). Atrogin-1, a muscle-specific E3 ubiquitin ligase that is upregulated during atrophy, may contribute to this process because it initiates CnA degradation in cardiomyocytes (44, 45). Therefore, a reduction in CnA could influence the fiber switching that occurs in aging. Taken together, these finding suggest that skeletal muscle may be remodeling toward a slow fiber phenotype in the KE group, potentially protecting against denervation induced muscle fiber loss during sarcopenia. The implications of these findings are not entirely clear and require further exploration to determine if these responses to BD-AcAc2 mitigate age-related reductions in skeletal muscle mass and function.
Our results also highlighted the activation of atrophy pathways. For example, we found a small yet significant decrease in Calcium/Calmodulin Dependent Protein Kinase II Beta (Camk2b) mRNA expression in KE, which has been identified as a downstream target of p38 mitogen-activated protein kinase alpha (p38α MAPK) and positive regulator of muscle atrophy (46). Furthermore, pharmacological inhibition of CAMK2B activity suppresses denervation-induced muscle atrophy (47). SMAD Family Member 3 (Smad3) was also 1.25-fold lower in KE compared to control, which could translate to less atrophy since activation of Smad proteins inhibit the function of Akt and the expression of Atrogin-1 and Muscle-specific RING finger protein 1 (MuRF1) by Forkhead box O transcription factors (48). The COP9 singalosome (Cops4), which is similar in structure and function to that of the 19S regulatory particle of 26S proteasome, was 2.63-fold higher in KE and has been shown to interact with SCF-type E3 ubiquitin ligases and act as a positive regulator of E3 ubiquitin ligases. Taken together, these findings suggest that the KE may protect against muscle atrophy in aged mice even though skeletal muscle weight was not different between groups, nor were Atrogin-1 and MAFbx mRNA expression. Previous reports indicate that at least one form of KE can attenuate muscle wasting with some diseases in humans (49).
Several transcripts in the KE group were involved in myogenesis, muscle repair, and muscle stem cells. For example, Pim-1 oncogene protein (Pim1) was 1.56-fold lower in the current study and Pim1 is required for proper myogenesis (50). Furthermore, Prominin 1 (Prom1) expression was 1.37-fold higher in KE, which is required for stem cell maintenance and activation (51) along with a reduction in the transcripts for MDS1 And EVI1 Complex Locus Protein EV1 (Mecom) and Rac Family Small GTPase 3 (Rac3), both of which have been implicated in stem cell control. Our proteomic data support these findings, with Galectin-1 protein content lower in KE mice than CON. Galectin-1 regulates myotube growth in regenerating skeletal muscle (52) and treatment with galectin-1 improves myogenic potential in some mouse models (53). This pathway acts through the Annexin family, and we found that Annexin A1 and Annexin 2 mRNA expression were downregulated, in KE mice, which could partially explain our proteomic findings with galectin-1 (52). Annexin A1 affects myoblast fusion causing a slowdown in regeneration of injured muscle but does not cause muscle damage or decrease the ability of injured myofibers to repair (54), while a lack of Annexin 2 results in poor myofiber repair and progressive muscle weakening with age (55). Our proteomic data partially support these findings since Annexin A1 was lower in the KE group compared to CON. PDZ and LIM domain 5 (PDLIM5) plays a positive role in the proliferation and differentiation of skeletal muscle satellite cells, and our proteomic data show it was 1.5-fold lower in KE compared to control (56). The silencing of PDLIM5 increases the nuclear accumulation of differentiation inhibitor (Id2), which inhibits the proliferation and differentiation of myoblasts. Another protein, Fibrillin 1, was higher in KE mice and it influences the skeletal muscle stem cell microenvironment that interacts with bone (57). On the other hand, transgelin was lower in KE mice and previous studies have indicated overexpression or deficiency of transgelin results in significant changes in cell morphology, cytoskeleton, and functional capacity for migration that influence stem cell differentiation (58). Collectively, these data suggest that the KE may affect muscle stem cells or muscle regeneration capacity such that the KE diet is remodeling the muscle stem cell niche in a positive manner, but more research is needed to confirm this speculation.
The ability of the metabolism of ketone bodies to increase insulin sensitivity has been demonstrated previously in the working perfused rat heart (59). This may act through Casitas b-lineage (Cbl) lymphoma, which is involved in regulating the degradation of receptor tyrosine kinases. Genetic mutation of the Cbl gene in mice causes a lean phenotype and enhanced peripheral insulin action potentially due to elevated energy expenditure (60, 61). While our previous findings indicate no differences in energy expenditure in younger mice consuming a KE, we found Cbl mRNA lower in KE mice compared to control. Interleukin-1 receptor-associated kinase-like 2 (Irak2) mRNA expression was also lower in the KE compared to CON and others have shown that Irak2 reduces high-fat-diet-induced weight gain when inactivated (62, 63). Overall, these findings suggest that BD-AcAc2 could be influencing muscle metabolism to improve glucose tolerance and insulin sensitivity. Our findings that control mice gained fat while KE did not gain fat support this possibility and more in-depth analysis with other-omics is warranted.
The lipidomics analysis performed in the current study may also garner insight about the effects of BD-AcAc2 on insulin signaling/sensitivity. While BD-AcAc2 supplementation appears to have little effect on lipid species within skeletal muscle; the analysis did show that fatty acyl chains in triacylglycerol (TAG) were significantly lower (46%) in the KE group compared to CON. Given our findings on transcripts related to insulin resistance, and that intracellular TAG content is associated with insulin resistance in skeletal muscle, these findings support the overall concept that the KE may improve insulin sensitivity. Additional studies will be required to determine if these effects are direct or if these findings are due to decreased FM (64), or reductions in hepatic lipid content and inflammatory markers (65). Our proteomics data show Acyl-coenzyme A thioesterase 13 (Acot13), a member of the Acyl-CoA thioesterases family that hydrolyze fatty acyl-CoA esters, was lower in KE mice compared to CON (66, 67). Interestingly, mice with whole-body Acot13 deletion have improved glucose homeostasis on a high-fat diet (68). Therefore, lower Acot13 protein content in skeletal muscle of the KE group supports our transcriptomic and lipidomic data.
Our comprehensive characterization of skeletal muscle also suggests that at least some portions of autophagy signaling are affected by BD-AcAc2. For example, disheveled segment polarity protein 2 (Dvl2) and TNF receptor associated factor 7 (Traf7) were significantly lower in KE compared to CON. Dvl-2 functions as part of the upstream Wnt signal transduction pathway of β-catenin and GSK-3β, and can positively regulate the Wnt signal pathway (69). TRAF7 regulates skeletal muscle through its activity as a ubiquitin ligase and depletion of Traf7 accelerates myogenesis, in part through downregulation of nuclear factor-κB (NF-κB) activity (70, 71). The downregulation of Traf7, along with decreased Dvl2 in the KE group suggests that BD-AcAc2 may reduce cellular stress. Furthermore, profilin 1 protein was lower in the KE group. Increased expression of profilin-1 in vascular smooth muscle cells induces stress fiber formation and causes cardiac hypertrophy and fibrosis by modulating actin polymerization (72, 73). Histone deacetylase 6 (Hdac6) expression was also lower, and it controls fusion of autophagosomes to lysosomes (74). Evidence suggests that HDAC6 also plays a role in the eventual clearance of aggresomes, implying a functional connection between HDAC6 and autophagy (74–76). While we can only speculate about the possible role of BD-AcAc2 on cellular stress, others have found that a KD attenuates the increase in cellular stress due to aging (10).
Other proteins that were higher in the KE group compared to CON involved elements of the electron transport chain (ETC). For example, crystallin zeta (CryZ), a quinone oxidoreductase and threodoxi-like 1, which controls disulfide oxidoreductate activity was higher in KE compared to CON at the protein level. Another protein elevated in the KE group was voltage-dependent anion channel 1 (VDAC1) which controls the metabolic crosstalk between mitochondria and the cytosol, by regulating the influx and efflux of metabolites, cations, and nucleotides (77). Human VDAC1 deficiency compromises pyruvate oxidation and ATP production (78). Others have shown VDAC1 deficiency causes multiple defects in the mitochondria and ETC complexes in both oxidative and glycolytic striated muscle biopsies. This includes loss of ion homeostasis resulting in irreparable cell injury and concomitantly cell death, primarily through necrosis. Higher expression of electron transfer flavoprotein dehydrogenase (ETFDH) mRNA in KE mice suggests improvement of skeletal muscle mitochondria β-oxidation in the ETC (79) while an upregulation of Bcl-2-associated X protein (Bax) suggests that the KE may enhance regulation of mitochondrial morphology during apoptosis (80). Given findings that a basal amount of BAX is necessary to maintain energy production via aerobic respiration this may be important (81). Taken together, these data indicate that the KE is influencing the ETC, which is similar to what has been reported with a KD in aging mice (10).
Nutritional ketosis has been defined as a metabolic state where β-hydroxybutyrate (β-HB) blood concentrations rise above 0.5 millimole per liter (mM/L) (82, 83). In human clinical trials, circulating ketones such as β-HB reach levels ranging from 1 to 2 mM/L and up to 5 mM/L (6, 84, 85). Comparatively, we (and others) have shown the ability to induce similar levels of ketones in rodents (15, 17, 19, 86). Exogenous ketones offer a way to induce ketone production quickly and efficiently in humans without the need for caloric restriction or a KD and have been well-studied in exercise physiology (11, 12, 87–89). While more research is needed, KE have been suggested as a novel method to treat several diseases, including perturbations that occur in aging (90–94).
A strength of this research is that we took an -omics based approach to determine how a ketone ester influences skeletal muscle, providing multiple avenues for future research. A limitation of this approach is not identifying a specific molecular mechanism by which the KE is acting–although based on our findings it seems the BD-AcAc2 is influencing several pathways in skeletal muscle that may be responsible for the physiological and metabolic responses observed. Another limitation is the lack of baseline data on FM and LBM, which was due to COVID-19 shutdowns in the UAB Animal Resources Program. Nonetheless, mice were randomized so this should equally distribute confounding variables at baseline and our previous work in this rodent and dietary model indicates a very homogenous sample in these mice. In addition, body weights of the mice were similar among groups at the start of the study and in our experience with B6 mice, FM of these animals is highly correlated with overall body mass and should not represent a confounding variable for this study. A limitation that muddies the interpretation of our findings are the transcriptional and proteomic analysis being completed in different muscles due to insufficient tissue, so future experiments should include this outcome. The different tissue usage may explain the lack of correlation between our proteomic and mRNA data. Finally, the study is likely underpowered given the number of-omics approaches we utilized and given the lack of previous data on KE and muscle mass at the time of study initiation, we did not have sufficient data for a traditional power calculation.
In conclusion, our findings indicate that a diet containing 25% BD-AcAc2 in male aging B6 mice produces a unique skeletal muscle transcriptional, lipidomic and proteomic signature that influences signaling involved in muscle denervation as well as muscle atrophy and autophagy. Furthermore, KE may improve the satellite cell niche and β-oxidation in the ETC. Our lipidomic data and body composition data align with transcriptional findings that suggest that BD-AcAc2 may play a role in insulin sensitivity. Future studies are needed to examine whether insulin sensitivity is indeed improved in skeletal muscle or whether these findings are related to decreased adiposity and/or fat content in liver.
Data availability statement
The data presented in the study are deposited in the Open Science Framework repository (https://osf.io/e2yd5/).
Ethics statement
The animal study was reviewed and approved by UAB Institutional Animal Care and Use Committee (IACUC).
Author contributions
BR, SD, DS, and EP conceptualized the study. BR and SD conducted the experiments and collected and analyzed the data. NM contributed to the lipidomics data analysis and interpretation. JM contributed to the proteomics data analysis and interpretation. BR wrote the manuscript. All authors edited and approved the final version prior to submission for publication.
Funding
Research reported in this [publication/press release] was supported by the National Institute on Aging of the National Institutes of Health under award P30 AG050886, P30 AG003319, and P30 AG044271. NM also was supported by U01AR071130 and U01AR071150. Comprehensive Cancer Center Core Support Grant: Mass Spectrometry/Proteomics Shared Facility P30 CA013148. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
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.
Abbreviations
BD-AcAc2, R,S-1,3, butanediol diacetoacetate; KE, ketone ester; LBM, lean body mass; FM, fat mass; KD, ketogenic diet; Tnc, tenascin-C; ER, endoplasmic reticulum; Gadd45a, growth arrest and DNA damage inducible alpha; Runx1, family Transcription Factor 1; Chasm, calponin homology-associated smooth muscle protein; CnA, calcineurin A; PGC-1α, peroxisome proliferator-activated receptor-gamma coactivator; MEF2, myocyte enhancer factor-2; NFAT, nuclear factor of activated T-cells; Atrogin-1/MAFbx, Muscle Atrophy F-box gene; Pim1, pim-1 oncogene protein; Mecom, MDS1 And EVI1 Complex Locus Protein EV1; Camk2b, Calcium/Calmodulin Dependent Protein Kinase II Beta; P38α MAPK, p38 mitogen-activated protein kinases; Irak2, Interleukin-1 receptor-associated kinase-like 2; MuRF1, Muscle-specific RING finger protein 1; Cops9, COP9 singalosome; Dvl1, disheveled segment polarity protein 2 (Dvl2); Traf7, TNF receptor associated factor 7; Hdac6, histone deacetylase 6; CryZ, crystallin zeta; Bax, bcl-2-associated X protein; VDAC1, voltage-dependent anion channel 1.
Footnotes
References
1. Kennedy AR, Pissios P, Otu H, Roberson R, Xue B, Asakura K, et al. A high-fat, ketogenic diet induces a unique metabolic state in mice. Am J Physiol Endocrinol Metab. (2007) 292:E1724–39. doi: 10.1152/ajpendo.00717.2006
2. Murabito JM, Yuan R, Lunetta KL. The search for longevity and healthy aging genes: insights from epidemiological studies and samples of long-lived individuals. J Gerontol A Biol Sci Med Sci. (2012) 67:470–9. doi: 10.1093/gerona/gls089
3. Dowis K, Banga S. The potential health benefits of the ketogenic diet: a narrative review. Nutrients. (2021) 13:1654.
4. Castellana M, Conte E, Cignarelli A, Perrini S, Giustina A, Giovanella L, et al. Efficacy and safety of very low calorie ketogenic diet (VLCKD) in patients with overweight and obesity: a systematic review and meta-analysis. Rev Endocr Metab Disord. (2020) 21:5–16. doi: 10.1007/s11154-019-09514-y
5. Gomez-Arbelaez D, Bellido D, Castro AI, Ordoñez-Mayan L, Carreira J, Galban C, et al. Body composition changes after very-low-calorie ketogenic diet in obesity evaluated by 3 standardized methods. J Clin Endocrinol Metab. (2017) 102:488–98. doi: 10.1210/jc.2016-2385
6. Lorenzo PM, Sajoux I, Izquierdo AG, Gomez-Arbelaez D, Zulet MA, Abete I, et al. Immunomodulatory effect of a very-low-calorie ketogenic diet compared with bariatric surgery and a low-calorie diet in patients with excessive body weight. Clin Nutr. (2022) 41:1566–77. doi: 10.1016/j.clnu.2022.05.007
7. Sajoux I, Lorenzo PM, Gomez-Arbelaez D, Zulet MA, Abete I, Castro AI, et al. Effect of a very-low-calorie ketogenic diet on circulating myokine levels compared with the effect of bariatric surgery or a low-calorie diet in patients with obesity. Nutrients. (2019) 11:2368.
8. Crujeiras AB, Izquierdo AG, Primo D, Milagro FI, Sajoux I, Jácome A, et al. Epigenetic landscape in blood leukocytes following ketosis and weight loss induced by a very low calorie ketogenic diet (VLCKD) in patients with obesity. Clin Nutr. (2021) 40:3959–72. doi: 10.1016/j.clnu.2021.05.010
9. Roberts BM, Lavin KM, Many GM, Thalacker-Mercer A, Merritt EK, Bickel CS, et al. Human neuromuscular aging: sex differences revealed at the myocellular level. Exp Gerontol. (2018) 106:116–24. doi: 10.1016/j.exger.2018.02.023
10. Wallace MA, Aguirre NW, Marcotte GR, Marshall AG, Baehr LM, Hughes DC, et al. The ketogenic diet preserves skeletal muscle with aging in mice. Aging Cell. (2021) 20:e13322.
11. Margolis LM, O’Fallon KS. Utility of ketone supplementation to enhance physical performance: a systematic review. Adv Nutr. (2020) 11:412–9. doi: 10.1093/advances/nmz104
12. Clarke K, Tchabanenko K, Pawlosky R, Carter E, Todd King M, Musa-Veloso K, et al. Kinetics, safety and tolerability of (R)-3-hydroxybutyl (R)-3-hydroxybutyrate in healthy adult subjects. Regul Toxicol Pharmacol. (2012) 63:401–8. doi: 10.1016/j.yrtph.2012.04.008
13. Hashim SA, VanItallie TB. Ketone body therapy: from the ketogenic diet to the oral administration of ketone ester. J Lipid Res. (2014) 55:1818–26. doi: 10.1194/jlr.R046599
14. Kesl SL, Poff AM, Ward NP, Fiorelli TN, Ari C, Van Putten AJ, et al. Effects of exogenous ketone supplementation on blood ketone, glucose, triglyceride, and lipoprotein levels in Sprague-Dawley rats. Nutr Metab (Lond). (2016) 13:9. doi: 10.1186/s12986-016-0069-y
15. Poff AM, Ari C, Arnold P, Seyfried TN, D’Agostino DP. Ketone supplementation decreases tumor cell viability and prolongs survival of mice with metastatic cancer. Int J Cancer. (2014) 135:1711–20. doi: 10.1002/ijc.28809
16. Ari C, Murdun C, Goldhagen C, Koutnik AP, Bharwani SR, Diamond DM, et al. Exogenous ketone supplements improved motor performance in preclinical rodent models. Nutrients. (2020) 12:2459. doi: 10.3390/nu12082459
17. Deemer SE, Davis RAH, Gower BA, Koutnik AP, Poff AM, Dickinson SL, et al. Concentration-dependent effects of a dietary ketone ester on components of energy balance in mice. Front Nutr. (2019) 6:56. doi: 10.3389/fnut.2019.00056
18. Davis RAH, Deemer SE, Bergeron JM, Little JT, Warren JL, Fisher G, et al. Dietary R, S-1,3-butanediol diacetoacetate reduces body weight and adiposity in obese mice fed a high-fat diet. FASEB J. (2019) 33:2409–21. doi: 10.1096/fj.201800821RR
19. Deemer SE, Davis RAH, Roberts BM, Smith DL Jr., Koutnik AP, Poff AM, et al. Exogenous dietary ketone ester decreases body weight and adiposity in mice housed at thermoneutrality. Obesity (Silver Spring). (2020) 28:1447–55. doi: 10.1002/oby.22855
20. Perkins JR, Dawes J, McMahon S, Bennett D, Orengo C, Kohl M, et al. ReadqPCR and NormqPCR: R packages for the reading, quality checking and normalisation of RT-qPCR quantification cycle (Cq) data. BMC Genomics. (2012) 13:296. doi: 10.1186/1471-2164-13-296
21. Ludwig MR, Kojima K, Bowersock GJ, Chen D, Jhala NC, Buchsbaum DJ, et al. Surveying the serologic proteome in a tissue-specific kras(G12D) knockin mouse model of pancreatic cancer. Proteomics. (2016) 16:516–31. doi: 10.1002/pmic.201500133
22. Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. (2002) 74:5383–92. doi: 10.1021/ac025747h
23. Nesvizhskii AI, Keller A, Kolker E, Aebersold R. A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. (2003) 75:4646–58.
24. Weatherly DB, Atwood JAIII, Minning TA, Cavola C, Tarleton RL, Orlando R. A Heuristic method for assigning a false-discovery rate for protein identifications from Mascot database search results. Mol Cell Proteomics. (2005) 4:762–72. doi: 10.1074/mcp.M400215-MCP200
25. Old WM, Meyer-Arendt K, Aveline-Wolf L, Pierce KG, Mendoza A, Sevinsky JR, et al. Comparison of label-free methods for quantifying human proteins by shotgun proteomics. Mol Cell Proteomics. (2005) 4:1487–502.
26. Liu H, Sadygov RG, Yates JR III. A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem. (2004) 76:4193–201. doi: 10.1021/ac0498563
27. Beissbarth T, Hyde L, Smyth GK, Job C, Boon WM, Tan SS, et al. Statistical modeling of sequencing errors in SAGE libraries. Bioinformatics. (2004) 20(Suppl. 1):i31–9.
28. Bhatia VN, Perlman DH, Costello CE, McComb ME. Software tool for researching annotations of proteins: open-source protein annotation software with data visualization. Anal Chem. (2009) 81:9819–23.
29. Ekins S, Bugrim A, Brovold L, Kirillov E, Nikolsky Y, Rakhmatulin E, et al. Algorithms for network analysis in systems-ADME/Tox using the MetaCore and MetaDrug platforms. Xenobiotica. (2006) 36:877–901. doi: 10.1080/00498250600861660
30. Han X. Lipidomics: Comprehensive Mass Spectrometry of Lipids/Xianlin Han. Hoboken, NJ: John Wiley & Sons, Inc (2016).
31. Wang M, Han X. Multidimensional mass spectrometry-based shotgun lipidomics. Methods Mol Biol. (2014) 1198:203–20.
32. Han X, Yang K, Gross RW. Microfluidics-based electrospray ionization enhances the intrasource separation of lipid classes and extends identification of individual molecular species through multi-dimensional mass spectrometry: development of an automated high-throughput platform for shotgun lipidomics. Rapid Commun Mass Spectrom. (2008) 22:2115–24. doi: 10.1002/rcm.3595
33. Wang M, Wang C, Han RH, Han X. Novel advances in shotgun lipidomics for biology and medicine. Prog Lipid Res. (2016) 61:83–108. doi: 10.1016/j.plipres.2015.12.002
34. Yang K, Cheng H, Gross RW, Han X. Automated lipid identification and quantification by multidimensional mass spectrometry-based shotgun lipidomics. Anal Chem. (2009) 81:4356–68. doi: 10.1021/ac900241u
35. Barns M, Stewart B, Cernelc J, Yuminaga Y. Molecular analyses provide insight into mechanisms underlying sarcopenia and myofibre denervation in old skeletal muscles of mice. Int J Biochem Cell Biol. (2014) 53:174–85. doi: 10.1016/j.biocel.2014.04.025
36. Ehmsen JT, Kawaguchi R, Kaval D, Johnson AE, Nachun D, Coppola G, et al. GADD45A is a protective modifier of neurogenic skeletal muscle atrophy. JCI Insight. (2021) 6:e149381. doi: 10.1172/jci.insight.149381
37. Cifuentes-Diaz C, Faille L, Goudou D, Schachner M, Rieger F, Angaut-Petit D. Abnormal reinnervation of skeletal muscle in a tenascin-C-deficient mouse. J Neurosci Res. (2002) 67:93–9. doi: 10.1002/jnr.10109
38. Wang X, Blagden C, Fan J, Nowak SJ, Taniuchi I, Littman DR, et al. Runx1 prevents wasting, myofibrillar disorganization, and autophagy of skeletal muscle. Genes Dev. (2005) 19:1715–22. doi: 10.1101/gad.1318305
39. Grumati P, Bonaldo P. Autophagy in skeletal muscle homeostasis and in muscular dystrophies. Cells. (2012) 1:325–45.
40. Ulke-Lemee A, Ishida H, Borman MA, Valderrama A, Vogel HJ, MacDonald JA. Tropomyosin-binding properties of the CHASM protein are dependent upon its calponin homology domain. FEBS Lett. (2010) 584:3311–6. doi: 10.1016/j.febslet.2010.07.012
41. Rowan SL, Rygiel K, Purves-Smith FM, Solbak NM, Turnbull DM, Hepple RT. Denervation causes fiber atrophy and myosin heavy chain co-expression in senescent skeletal muscle. PLoS One. (2012) 7:e29082. doi: 10.1371/journal.pone.0029082
42. Aare S, Spendiff S, Vuda M, Elkrief D, Perez A, Wu Q, et al. Failed reinnervation in aging skeletal muscle. Skelet Muscle. (2016) 6:29.
43. Hudson MB, Price SR. Calcineurin: a poorly understood regulator of muscle mass. Int J Biochem Cell Biol. (2013) 45:2173–8.
44. Lin J, Wu H, Tarr PT, Zhang CY, Wu Z, Boss O, et al. Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature. (2002) 418:797–801. doi: 10.1038/nature00904
45. Li HH, Kedar V, Zhang C, McDonough H, Arya R, Wang DZ, et al. Atrogin-1/muscle atrophy F-box inhibits calcineurin-dependent cardiac hypertrophy by participating in an SCF ubiquitin ligase complex. J Clin Invest. (2004) 114:1058–71. doi: 10.1172/JCI22220
46. Clarke BA, Drujan D, Willis MS, Murphy LO, Corpina RA, Burova E, et al. The E3 ligase MuRF1 degrades myosin heavy chain protein in dexamethasone-treated skeletal muscle. Cell Metab. (2007) 6:376–85.
47. Yuasa K, Okubo K, Yoda M, Otsu K, Ishii Y, Nakamura M, et al. Targeted ablation of p38alpha MAPK suppresses denervation-induced muscle atrophy. Sci Rep. (2018) 8:9037. doi: 10.1038/s41598-018-26632-w
48. Egerman MA, Glass DJ. Signaling pathways controlling skeletal muscle mass. Crit Rev Biochem Mol Biol. (2014) 49:59–68.
49. Koutnik AP, Poff AM, Ward NP, DeBlasi JM, Soliven MA, Romero MA, et al. Ketone bodies attenuate wasting in models of atrophy. J Cachexia Sarcopenia Muscle. (2020) 11:973–96. doi: 10.1002/jcsm.12554
50. Liu Y, Shang Y, Yan Z, Li H, Wang Z, Liu Z, et al. Pim1 kinase positively regulates myoblast behaviors and skeletal muscle regeneration. Cell Death Dis. (2019) 10:773. doi: 10.1038/s41419-019-1993-3
51. Singer D, Thamm K, Zhuang H, Karbanová J, Gao Y, Walker JV, et al. Prominin-1 controls stem cell activation by orchestrating ciliary dynamics. EMBO J. (2019) 38:e99845. doi: 10.15252/embj.201899845
52. Kami K, Senba E. Galectin-1 is a novel factor that regulates myotube growth in regenerating skeletal muscles. Curr Drug Targets. (2005) 6:395–405. doi: 10.2174/1389450054021918
53. Vallecillo-Zuniga ML, Rathgeber M, Poulson D, Kartchner B, Luddington J, Gill H, et al. Treatment with galectin-1 improves myogenic potential and membrane repair in dysferlin-deficient models. PLoS One. (2020) 15:e0238441. doi: 10.1371/journal.pone.0238441
54. Leikina E, Defour A, Melikov K, Van der Meulen JH, Nagaraju K, Bhuvanendran S, et al. Annexin A1 Deficiency does not Affect Myofiber Repair but Delays Regeneration of Injured Muscles. Sci Rep. (2015) 5:18246.
55. Defour A, Medikayala S, Van der Meulen JH, Hogarth MW, Holdreith N, Malatras A, et al. Annexin A2 links poor myofiber repair with inflammation and adipogenic replacement of the injured muscle. Hum Mol Genet. (2017) 26:1979–91. doi: 10.1093/hmg/ddx065
56. He H, Yin H, Yu X, Zhang Y, Ma M, Li D, et al. PDLIM5 affects chicken skeletal muscle satellite cell proliferation and differentiation via the p38-MAPK pathway. Animals (Basel). (2021) 11:1016. doi: 10.3390/ani11041016
57. Nakatani M, Ito J, Koyama R, Iijima M, Yoshimoto N, Niimi T, et al. Scaffold protein enigma homolog 1 overcomes the repression of myogenesis activation by inhibitor of DNA binding 2. Biochem Biophys Res Commun. (2016) 474:413–20. doi: 10.1016/j.bbrc.2016.04.119
58. Elsafadi M, Manikandan M, Dawud RA, Alajez NM, Hamam R, Alfayez M, et al. Transgelin is a TGFbeta-inducible gene that regulates osteoblastic and adipogenic differentiation of human skeletal stem cells through actin cytoskeleston organization. Cell Death Dis. (2016) 7:e2321. doi: 10.1038/cddis.2016.196
60. Molero JC, Waring SG, Cooper A, Turner N, Laybutt R, Cooney GJ, et al. Casitas b-lineage lymphoma-deficient mice are protected against high-fat diet-induced obesity and insulin resistance. Diabetes. (2006) 55:708–15. doi: 10.2337/diabetes.55.03.06.db05-0312
61. Molero JC, Jensen TE, Withers PC, Couzens M, Herzog H, Thien CB, et al. c-Cbl-deficient mice have reduced adiposity, higher energy expenditure, and improved peripheral insulin action. J Clin Invest. (2004) 114:1326–33. doi: 10.1172/JCI21480
62. Sun XJ, Kim SP, Zhang D, Sun H, Cao Q, Lu X, et al. Deletion of interleukin 1 receptor-associated kinase 1 (Irak1) improves glucose tolerance primarily by increasing insulin sensitivity in skeletal muscle. J Biol Chem. (2017) 292:12339–50. doi: 10.1074/jbc.M117.779108
63. Zhou H, Wang H, Yu M, Schugar RC, Qian W, Tang F, et al. IL-1 induces mitochondrial translocation of IRAK2 to suppress oxidative metabolism in adipocytes. Nat Immunol. (2020) 21:1219–31. doi: 10.1038/s41590-020-0750-1
64. Watt MJ, Hoy AJ. Lipid metabolism in skeletal muscle: generation of adaptive and maladaptive intracellular signals for cellular function. Am J Physiol Endocrinol Metab. (2012) 302:E1315–28.
65. Moore MP, Cunningham RP, Davis RAH, Deemer SE, Roberts BM, Plaisance EP, et al. A dietary ketone ester mitigates histological outcomes of NAFLD and markers of fibrosis in high-fat diet fed mice. Am J Physiol Gastrointest Liver Physiol. (2021) 320:G564–72. doi: 10.1152/ajpgi.00259.2020
66. Bekeova C, Anderson-Pullinger L, Boye K, Boos F, Sharpadskaya Y, Herrmann JM, et al. Multiple mitochondrial thioesterases have distinct tissue and substrate specificity and CoA regulation, suggesting unique functional roles. J Biol Chem. (2019) 294:19034–47. doi: 10.1074/jbc.RA119.010901
67. Bakshi I, Brown SHJ, Brandon AE, Suryana E, Mitchell TW, Turner N, et al. Increasing Acyl CoA thioesterase activity alters phospholipid profile without effect on insulin action in skeletal muscle of rats. Sci Rep. (2018) 8:13967. doi: 10.1038/s41598-018-32354-w
68. Kang HW, Niepel MW, Han S, Kawano Y, Cohen DE, et al. Thioesterase superfamily member 2/acyl-CoA thioesterase 13 (Them2/Acot13) regulates hepatic lipid and glucose metabolism. FASEB J. (2012) 26:2209–21.
69. Smalley MJ, Sara E, Paterson H, Naylor S, Cook D, Jayatilake H, et al. Interaction of axin and Dvl-2 proteins regulates Dvl-2-stimulated TCF-dependent transcription. EMBO J. (1999) 18:2823–35. doi: 10.1093/emboj/18.10.2823
70. Tsikitis M, Acosta-Alvear D, Blais A, Campos EI, Lane WS, Sánchez I, et al. Traf7, a MyoD1 transcriptional target, regulates nuclear factor-kappaB activity during myogenesis. EMBO Rep. (2010) 11:969–76. doi: 10.1038/embor.2010.154
71. Suryadevara V, Willis MS. Walk the line: the role of ubiquitin in regulating transcription in myocytes. Physiology (Bethesda). (2019) 34:327–40. doi: 10.1152/physiol.00055.2018
72. Kooij V, Viswanathan MC, Lee DI, Rainer PP, Schmidt W, Kronert WA, et al. Profilin modulates sarcomeric organization and mediates cardiomyocyte hypertrophy. Cardiovasc Res. (2016) 110:238–48. doi: 10.1093/cvr/cvw050
73. Moustafa-Bayoumi M, Alhaj MA, El-Sayed O, Wisel S, Chotani MA, Abouelnaga ZA, et al. Vascular hypertrophy and hypertension caused by transgenic overexpression of profilin 1. J Biol Chem. (2007) 282:37632–9. doi: 10.1074/jbc.M703227200
74. Lee JY, Koga H, Kawaguchi Y, Tang W, Wong E, Gao YS, et al. HDAC6 controls autophagosome maturation essential for ubiquitin-selective quality-control autophagy. EMBO J. (2010) 29:969–80. doi: 10.1038/emboj.2009.405
75. Walsh ME, Van Remmen H. Emerging roles for histone deacetylases in age-related muscle atrophy. Nutr Healthy Aging. (2016) 4:17–30. doi: 10.3233/NHA-160005
76. Pandey UB, Nie Z, Batlevi Y, McCray BA, Ritson GP, Nedelsky NB, et al. HDAC6 rescues neurodegeneration and provides an essential link between autophagy and the UPS. Nature. (2007) 447:859–63. doi: 10.1038/nature05853
77. Camara AKS, Zhou Y, Wen PC, Tajkhorshid E, Kwok WM. Mitochondrial VDAC1: a key gatekeeper as potential therapeutic target. Front Physiol. (2017) 8:460. doi: 10.3389/fphys.2017.00460
78. Huizing M, Ruitenbeek W, Thinnes FP, DePinto V, Wendel U, Trijbels FJ, et al. Deficiency of the voltage-dependent anion channel: a novel cause of mitochondriopathy. Pediatr Res. (1996) 39:760–5. doi: 10.1203/00006450-199605000-00003
79. Missaglia S, Pegoraro V, Marozzo R, Tavian D, Angelini C. Correlation between ETFDH mutations and dysregulation of serum myomiRs in MADD patients. Eur J Transl Myol. (2020) 30:8880. doi: 10.4081/ejtm.2019.8880
80. Zhang Y, Iqbal S, O’Leary MF, Menzies KJ, Saleem A, Ding S, et al. Altered mitochondrial morphology and defective protein import reveal novel roles for Bax and/or Bak in skeletal muscle. Am J Physiol Cell Physiol. (2013) 305:C502–11. doi: 10.1152/ajpcell.00058.2013
81. Boohaker RJ, Zhang G, Carlson AL, Nemec KN, Khaled AR. BAX supports the mitochondrial network, promoting bioenergetics in nonapoptotic cells. Am J Physiol Cell Physiol. (2011) 300:C1466–78. doi: 10.1152/ajpcell.00325.2010
82. Poff AM, Koutnik AP, Egan B. Nutritional ketosis with ketogenic diets or exogenous ketones: features, convergence, and divergence. Curr Sports Med Rep. (2020) 19:251–9.
83. Newman JC, Verdin E. beta-hydroxybutyrate: a signaling metabolite. Annu Rev Nutr. (2017) 37:51–76.
84. Laeger T, Metges CC, Kuhla B. Role of beta-hydroxybutyric acid in the central regulation of energy balance. Appetite. (2010) 54:450–5.
85. Paoli A. Ketogenic diet for obesity: friend or foe? Int J Environ Res Public Health. (2014) 11:2092–107.
86. Koutnik AP, D’Agostino DP, Egan B. Anticatabolic effects of ketone bodies in skeletal muscle. Trends Endocrinol Metab. (2019) 30:227–9.
87. Evans M, Cogan KE, Egan B. Metabolism of ketone bodies during exercise and training: physiological basis for exogenous supplementation. J Physiol. (2017) 595:2857–71. doi: 10.1113/JP273185
88. Stubbs BJ, Cox PJ, Evans RD, Santer P, Miller JJ, Faull OK, et al. On the metabolism of exogenous ketones in humans. Front Physiol. (2017) 8:848. doi: 10.3389/fphys.2017.00848
89. Pinckaers PJ, Churchward-Venne TA, Bailey D, van Loon LJ. Ketone bodies and exercise performance: the next magic bullet or merely hype? Sports Med. (2017) 47:383–91.
90. Walsh JJ, Myette-Côté É, Neudorf H, Little JP. Potential therapeutic effects of exogenous ketone supplementation for type 2 diabetes: a review. Curr Pharm Des. (2020) 26:958–69. doi: 10.2174/1381612826666200203120540
91. Kovacs Z, Brunner B, Ari C. Beneficial effects of exogenous ketogenic supplements on aging processes and age-related neurodegenerative diseases. Nutrients. (2021) 13:2197. doi: 10.3390/nu13072197
92. Selvaraj S, Kelly DP, Margulies KB. Implications of altered ketone metabolism and therapeutic ketosis in heart failure. Circulation. (2020) 141:1800–12. doi: 10.1161/CIRCULATIONAHA.119.045033
93. Jensen NJ, Wodschow HZ, Nilsson M, Rungby J. Effects of ketone bodies on brain metabolism and function in neurodegenerative diseases. Int J Mol Sci. (2020) 21:8767.
Keywords: skeletal muscle, ketone ester, nutrition, sarcopenia, proteomics, lipidomics
Citation: Roberts BM, Deemer SE, Smith DL Jr, Mobley JA, Musi N and Plaisance EP (2022) Effects of an exogenous ketone ester using multi-omics in skeletal muscle of aging C57BL/6J male mice. Front. Nutr. 9:1041026. doi: 10.3389/fnut.2022.1041026
Received: 10 September 2022; Accepted: 25 October 2022;
Published: 15 November 2022.
Edited by:
Demin Cai, Yangzhou University, ChinaReviewed by:
Zsolt Kovacs, Eötvös Loránd University, HungaryAna B. Crujeiras, Complejo Hospitalario Universitario de Santiago, Spain
Copyright © 2022 Roberts, Deemer, Smith, Mobley, Musi and Plaisance. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eric P. Plaisance, cGxhaXNlcEB1YWIuZWR1