Abstract
Autophagy is an essential component of proteostasis and a key pathway in aging. Identifying associations between autophagy gene expression patterns in skeletal muscle and physical performance outcomes would further our knowledge of mechanisms related with proteostasis and healthy aging. Muscle biopsies were obtained from participants in the Study of Muscle, Mobility and Aging (SOMMA). For 575 participants, RNA was sequenced and expression of 281 genes related to autophagy regulation, mitophagy and mTOR/upstream pathways were determined. Associations between gene expression and outcomes including mitochondrial respiration in muscle fiber bundles (MAX OXPHOS), physical performance (VO2 peak, 400m walking speed, and leg power), and thigh muscle volume were determined using negative binomial regression models. For autophagy, key transcriptional regulators including TFE3 and NFKB-related genes (RELA, RELB, NFKB1) were negatively associated with outcomes. On the contrary, regulators of oxidative metabolism that also promote overall autophagy, mitophagy and pexophagy (PPARGC1A, PPARA, EPAS1) were positively associated with multiple outcomes. In line with this, several mitophagy, fusion and fission related genes (NIPSNAP2, DNM1L, OPA1) were also positively associated with outcomes. For mTOR pathway and related genes, expression of WDR59 and WDR24, both subunits of GATOR2 complex (an indirect inhibitor of mTORC1) and PRKAG3, which is a regulatory subunit of AMPK, were negatively correlated with multiple outcomes. Our study identifies autophagy and selective autophagy such as mitophagy gene expression patterns in human skeletal muscle related to physical performance, muscle volume and mitochondrial function in older persons which may lead to target identification to preserve mobility and independence.
INTRODUCTION
Aging is defined by a gradual loss of physiological integrity, with accrual of cellular damage widely considered the general cause. Several molecular hallmarks of aging have been proposed that contribute to accumulation of cellular damage, including genomic instability, loss of proteostasis, macroautophagy dysregulation, DNA damage and mitochondrial dysfunction (Lopez-Otin et al. 2023). In skeletal muscle, aging results in a loss of muscle mass, strength and oxidative capacity that contributes to lower cardiorespiratory fitness, slower walking speed and ultimately mobility limitations. These changes in muscle are linked to impaired protein homeostasis or proteostasis (Powers et al. 2009) including the loss of skeletal muscle proteins (Rooyackers et al. 1996; Balagopal et al. 1997; Short et al. 2005; Drummond et al. 2009), and alterations in the balance between protein synthesis and degradation (Rooyackers et al. 1996; Rasmussen et al. 2006). The molecular basis for these changes, including changes in the expression of several key genes and proteins implicated in growth, atrophy, proteosome degradation and mitochondrial metabolism, is also documented (Sandri 2002; Giresi et al. 2005; Short et al. 2005; Masiero et al. 2009; Sandri 2010).
Autophagy is the cellular process of lysosomal degradation and recycling of cytoplasmic components to maintain cellular homeostasis. Autophagy clears and recycles damaged organelles and macromolecules and preserves DNA stability (Sandri 2010; Schneider & Cuervo 2014; Kaushik et al. 2021). Autophagy decreases with aging in several tissues, however its role in regulating muscle mass and function remains poorly understood, particularly in humans. On one hand, excessive activation of autophagy aggravates muscle wasting (Mammucari et al. 2007) by removing portions of cytoplasm, proteins, and organelles.
Conversely, it was demonstrated that inhibition of autophagy can result in muscle degeneration and weakness (Masiero et al. 2009). Blocking autophagy in muscle can also lead to impacts on innervation, mitochondrial function, and oxidative damage (Carnio et al. 2014; Franco-Romero & Sandri 2021).
Autophagy is highly coordinated and is positively regulated by the energy sensor, AMP-activated protein kinase (AMPK) pathway, and negatively regulated by the nutrient sensor, the mammalian target of rapamycin (mTOR) pathway (Jung et al. 2010; Egan et al. 2011). The transcription factors, TFEB and TFE3, are recognized for their role regulating a network of genes (the CLEAR network) involved in lysosomal biogenesis, autophagy and lysosomal exocytosis (Sardiello et al. 2009). The nuclear localization of these factors is modulated by phosphorylation via both AMPK and mTORC pathways. The study of autophagy has led to the identification of multiple autophagy sub-types that selectively degrade organelles including mitophagy, the process of mitochondrial degradation, and pexophagy, which is the targeted degradation of peroxisomes. PGC-1α is a regulator of mitochondrial biogenesis and has been shown to also activate mitochondria fusion-fission events (Martin et al. 2014), autophagy and in particular mitophagy (Vainshtein et al. 2015; Yeo et al. 2019). While the extensive study of the regulation of autophagic flux (the combined process of autophagosome formation and clearance) in model systems has revealed this regulation to be complex and multifactorial, the relevant pathways that regulate autophagy and mitophagy in human muscle aging have yet to be fully elucidated.
To further understand the gene expression profiles related to autophagy in human muscle and how they associate with muscle and physical function phenotypes, we performed sequencing of RNA obtained from muscle biopsies collected from participants in the Study of Muscle, Mobility and Aging (SOMMA). SOMMA is a prospective, longitudinal study of older people at risk of major mobility disability designed to understand the contributions of skeletal muscle mass and key properties of muscle tissue from biopsies to mobility. Decreases in autophagy occur with aging and a few small studies have correlated protein markers of autophagy with compromised muscle function (Zeng et al. 2018; Aas et al. 2019). In this investigation, we examined associations between expression levels of 281 genes involved in autophagy with muscle mitochondrial function, 400-m walking speed, VO2 peak, leg strength, and thigh muscle volume. Here, we took a targeted candidate gene approach and hypothesized that expression of autophagy genes in the muscle of older persons are generally associated with muscle mitochondrial function and tissues beyond muscle, thereby substantially contributing to overall fitness, walking speed, strength, and muscle mass.
RESULTS
Participant characteristics
A total of 879 participants provided consent and completed baseline measurements across both clinical sites (Figure 1). Of the 879 participants with complete baseline measures, 591 participants had RNA sequencing completed and 575 of these had high quality sequencing and complete covariate data. The characteristics of the study population with complete RNA sequencing and complete covariate data are presented in Table 1.
SOMMA participants selected for cross-sectional analyses.
Data shown as n (%), mean +/- Standard Deviation. P-values are presented for variables across tertiles of P1-Max OXPHOS. P-values for continuous variables from ANOVA for normally distributed data, a Kruskal-Wallis test for skewed data. P for linear trend across categories was calculated with linear regression models for those normally distributed variables, a Jonckheere-Terpstra test for skewed data. P-values for categorical data from a chi-square test for homogeneity. P for trend was calculated with the Jonckheere-Terpstra test.
RNA (human Ensembl genes (ENSG)) detection
The mean, median, and SD of the PCR duplicate percent per sample was 59%, 56% and 9%, respectively (Supplementary Table 1). After PCR duplicates were removed, the number of aligned reads per sample was high (mean=69,117,209, median = 71,313,059, SD = 14,444,848, range = 12,853,785-102,724,183).
Association of autophagy gene expression with multiple outcomes
We utilized a published curated list of autophagy genes (Bordi et al. 2021) to conduct a targeted analysis of the RNASeq dataset. A total of 281 genes related to autophagy regulation (68 genes), mitophagy (80 genes), and mTOR and upstream pathways (133 genes), were analyzed (Supplemental table 2). For each gene set, we examined the genes that were most significantly associated across multiple outcomes (Figure 2). All results report log base 2-fold changes reflecting the change in gene expression per one SD unit increase in each trait.
Heat map capturing top 20 statistically significant (p < 0.05 FDR adjusted) genes identified by our models: each color represents positive (red) or negative (blue) associations.
Analysis of the curated list of autophagy genes did not reveal any significant correlations between the core autophagy machinery (Beclin, Atg7, Atg5) and outcomes (Max OXPHOS, VO2 peak). For autophagy regulators, key transcriptional regulators including TFE3, TFEB, NFKB related genes (RELA, RELB, NFKB1), FOS, and FOXO1 were all significantly negatively associated with associated with 400m walk speed, VO2 peak, and Max OXPHOS. Alternatively, regulators of oxidative metabolism that also promote autophagy and mitophagy including PPARGC1A, PPARA, and EPAS1, a driver of pexophagy(Schonenberger et al. 2015), were all positively associated with VO2 peak, and Max OXPHOS. We also found that several mitochondria, fusion and fission related genes (NIPSNAP2, DNM1L, OPA1) were also positively associated with multiple outcomes including VO2 peak, Max OXPHOS and 400m walking speed. For mTOR pathway and related genes, expression of WDR59 and WDR24, both subunits of GATOR2 complex (an inhibitor of mTORC1) and PRKAG3, which is a regulatory subunit of AMPK, were negatively associated with multiple outcomes including 400m walk speed, VO2 peak, leg strength, leg power and thigh muscle volume. In contrast, MAPKAP1, which is a subunit of mTORC2 was positively associated with outcomes. A summary of all statistically significant associations for each gene set: autophagy regulation, mitophagy, and mTOR and upstream pathways with each trait is presented in Supplemental tables 3-8.
Association of mitochondrial respiration with gene expression
We next examined associations between gene expression for all three gene sets and specific outcomes. Maximal Complex I&II supported OXPHOS was positively associated with the expression of the deacetylases SIRT5, SIRT3, the regulators of mitochondrial dynamics MFN2, MUL1, and the component of the vacuolar proton pump ATPV1B1 (Figure 3). Conversely, OXPHOS was negatively associated with the transcription factors FOS, MYC and the regulator of protein translation EEF2K, among others. The full list of associations between autophagy genes and Max OXPHOS is presented in Supplementary Table 3.
Volcano plot capturing all statistically significant (p < 0.05 FDR adjusted) genes identified by our models: Each dot represents a gene; the dot color indicates significance level. Base model: gene expression∼ Max OXPHOS + age + gender + clinic site + race/ethnicity + height + weight + physical activities + multimorbidity count + sequencing batch.
Association of VO2 peak with gene expression
Autophagy genes that were positively associated with higher VO2 peak, included the autophagy regulators PPARGC1A, PPARA, the mitochondria dynamic regulators related to mitophagy DNM1L and OPA1, the mitophagy core gene NIPSNAP2, and the pexophagy-related hypoxia induced gene EPAS1 (Figure 4). Genes that were negatively associated with VO2 peak, included the autophagy regulators RELA, RELB, and the regulators of AMPK and mTOR, PRKAG3 and WDR24, respectively. The full list of associations between autophagy genes and Max OXPHOS is presented in Supplementary Table 4.
Volcano plot capturing all statistically significant (p < 0.05 FDR adjusted) genes identified by our models: Each dot represents a gene; the dot color indicates significance level. Base model: gene expression∼ VO2 Peak + age + gender + clinic site + race/ethnicity + height + weight + physical activities + multimorbidity count + sequencing batch.
Association of 400m walk speed with gene expression
Finally, we examined associations between autophagy-related gene expression and 400m walk speed. Genes that were positively associated with higher 400m walk speed, included mTOR and upstream pathway genes PRKAG2, PRKAB2, and TSC1 (Figure 5). Genes that were negatively associated with 400m walk speed, included positive (TFEC, TFE3) and negative (ATF5) autophagy regulators and the mTOR regulator SESN3. The full list of associations between autophagy genes and 400m walk speed is presented in Supplementary Table 5. The list of associations between autophagy genes and leg strength, leg power and thigh muscle volume are presented in Supplementary Tables 6-8.
Volcano plot capturing all statistically significant (p < 0.05 FDR adjusted) genes identified by our models: Each dot represents a gene; the dot color indicates significance level. Base model: gene expression∼ walk speed + age + gender + clinic site + race/ethnicity + height + weight + physical activities + multimorbidity count + sequencing batch.
DISCUSSION
Autophagy is a critical mechanism for maintaining cellular homeostasis, the dysregulation of which has been implicated as a key contributor to the biology of aging. Despite this, the relevance of skeletal muscle autophagy to measures of physical function, fitness and muscle mass and mitochondrial function in older adults has received relatively little attention. Here we leveraged transcriptomic data generated from muscle biopsy samples obtained from adults participating in the Study of Muscle, Mobility and Aging (SOMMA) to investigate the relationships between autophagy gene expression and multiple aging muscle and whole-body phenotypes (Cummings et al. 2023). We identified distinct patterns of association between expression of genes known to be key transcriptional regulators and signaling molecules of the autophagy program with several functional and physical domains. These data provide unique insight into the transcriptional regulation of autophagy in clinically relevant human aging phenotypes.
The principal finding was that several transcriptional regulators of the autophagy/lysosome pathway including TFE3, TFEB, FOS, and FOXO1 were all negatively associated with VO2 peak, Max OXPHOS, leg strength and thigh muscle volume. These data suggest that the downregulation of autophagy/lysosome transcriptional program is associated with better muscle and whole-body performance. However, because autophagy is part of the cellular response to several stressors including deprivation of nutrients or growth factors, and hypoxia and protects cellular homeostasis against oxidative stress and inflammation (Javali et al. 2023), lower expression of the autophagy/lysosomal transcriptional program could be a consequence of reduced overall cellular stress. Our results imply that older adults with higher fitness/respiration and walking speed may have lower levels of oxidative stress and inflammation and thus lower expression of drivers of the autophagy program. This finding is congruent with reports that older adults who engage in exercise and who are more fit have lower indices of inflammation including lower TNFα and IL-1β due to exercise-induced IL-6 release from muscle (Starkie et al. 2003; Pedersen & Fischer 2007). This data is also in-line with reports that higher levels of circulating inflammatory mediators are associated with poor mobility in older adults (Penninx et al. 2004). In further support of this possibility, we found modest or no correlation with expression of many of the core autophagy machinery, suggesting that it is mostly inducible autophagy (activated in response to stress) that is upregulated in adults with lower muscle fitness.
We found that expression of PPARA and PPARGC1A were positively associated with multiple outcomes including VO2 peak and Max OXPHOS. PPARGC1A, the gene that encodes PGC-1α, is most well-known for its role in mitochondrial biogenesis but has been reported to protect muscle mass in the context of muscle atrophy triggered by aging (Yang et al. 2020), chronic heart failure (Geng et al. 2011) and various other muscle wasting conditions (Sandri et al. 2006). More recently, it has been shown that PPARGC1A also plays a role in activating the autophagy-lysosome pathway (Takikita et al. 2010; Tsunemi et al. 2012) and specifically mitophagy in muscle (Vainshtein et al. 2015). In addition, PPARGC1A inhibits the expression of multiple pro-inflammatory mediators, including TNF-α and IL-6 (Handschin et al. 2007) and regulates antioxidant defense mechanisms (Kang & Li Ji 2012), both of which in turn can activate autophagy. In line with the role of PPARGC1A in mitophagy, we found that several mitochondria, fusion and fission related genes (DNM1L, OPA1) and the mitophagy “eat me” signal NIPSNAP2 were also positively associated with multiple outcomes, including VO2 peak, Max OXPHOS and 400m walking speed. In rodent models, deletion of the mitochondria fusion mediator OPA1 in skeletal muscle, caused mitochondrial dysfunction, oxidative stress and inflammation (Tezze et al. 2017; Rodriguez-Nuevo et al. 2018). DNM1L (or Drp-1) is a mediator of mitochondrial fission and muscle-specific Drp1 knockout in mice leads to a severe myopathic phenotype including muscle wasting and weakness (Favaro et al. 2019) and Drp1 knock down also decreases in mitochondrial respiration and increases markers of muscle denervation, fibrosis, and oxidative stress and also alters autophagy and mitophagy levels (Dulac et al. 2020). Taken together, these findings indicate that mitophagy and mitochondrial fusion and fission are essential for skeletal muscle health. However, studies in humans did not detect any impact of aging per se on Opa1 protein content in the skeletal muscle of sedentary young and old humans (Distefano et al. 2017), while physically active older adults had greater expression of Opa1 protein (Tezze et al. 2017). With respect to Drp1 expression, some studies report a decreased (Distefano et al. 2017) or an unchanged (Joseph et al. 2012) Drp1 content in skeletal muscle of older individuals, while physically active older adults had greater expression of Drp1 protein (Tezze et al. 2017). Based on these prior studies the consensus seems to be that mitochondrial dynamics do not play a role in muscle aging. However, unlike the previous study designs, our analysis focuses on associations between mitophagy gene expression and age-related phenotypes rather than aging per se.
mTORC1 is a master regulator of cell growth and metabolism through sensing and integrating different nutritional and environmental cues. Dysregulation and hyperactivation of mTOR contributes to several age-related diseases, such as cancer, neurodegenerative diseases, and type 2 diabetes mellitus. In skeletal muscle, mTORC1 signaling plays a critical role in protein synthesis and degradation, including regulating many steps of autophagy. The decline in skeletal muscle mass with aging is complex but anabolic resistance, or the limited ability to increase protein synthesis and inhibit autophagy in response to stimuli such as amino acids, is implicated as a contributing factor. In this study, we found that expression of mTOR pathway related genes WDR59 and WDR24 were both negatively correlated with 400m walk speed, VO2 peak, leg strength, leg power and thigh muscle volume. WDR59 and WDR24 are both subunits of GATOR2 complex which functions as a positive regulator of amino-acid-mediated mTORC1 activation (Bar-Peled et al. 2013) and WDR24 modulates mTOR through the previously described changes in lysosome cellular positioning and dynamics (Korolchuk et al. 2011), and thus indirectly modulates autophagic flux (Cai et al. 2016). This is in line with the observation that TFEB expression is also negatively associated with outcomes as mTORC1 regulates TFEB localization to the nucleus and activation by phosphorylating the transcription factor on several serine and threonine residues (Rabanal-Ruiz et al. 2017). Taken together, the associations of WDR59, WDR24 and TFEB with outcomes suggest that greater performance/fitness in older adults could be related to a blunted capacity of muscle mTORC1 to inhibit autophagy.
Our study has limitations. The analysis was performed on whole muscle tissue and the impact of mitochondrial content and muscle fiber type proportions on gene expression are not accounted for. For example, we observed that several genes whose expression correlated with Max OXPHOS are expressed exclusively within (or on) the mitochondria (SIRT5, SIRT3, MFN2, MUL1). That said, examining relationships between gene expression and relevant aging phenotypes in whole tissue initially is important to establish proof of concept/mechanism. A second limitation is that the study participants are of mostly White ancestry, thus potentially limiting generalizability of findings. Our data are cross-sectional and observational which limits our ability or prove causality. A strength of our study is that we paired muscle gene expression and mitochondrial function with several measures of fitness, strength, and muscle mass and accounted for potential confounding factors. Previous studies did not include older adults at risk of mobility disability and only analyzed cross-sectional associations between one or two properties and physical performance. Moreover, by focusing on a curated set of genes, rather than the broader gene ontology families, we were able to test specific hypotheses regarding the role of autophagy across mitochondrial function, fitness, mobility, strength, and muscle volume.
CONCLUSION
This study of autophagy genes in 575 SOMMA participants has revealed significant associations between autophagy regulation, mitophagy and mTOR pathways genes and a diverse range of clinically-relevant phenotypes that include walking speed, VO2 peak, maximal mitochondrial respiration, total thigh muscle volume, leg power and leg strength. Our results support and add to the evidence that transcriptional regulation of autophagy plays an important role in driving declines in various indices of muscle function that are key to mobility with aging. Additional studies are needed to further decipher the role of autophagy, including the contribution of post-transcriptional regulation, in skeletal muscle fitness during aging.
EXPERIMENTAL PROCEDURES
Study population
The Study of Muscle, Mobility and Aging (SOMMA) is a prospective cohort study of mobility in community-dwelling older adults. Participants for the current study were from the baseline cohort, enrolled between April 2019 and December 2021 (Cummings et al. 2023). SOMMA was conducted at 2 clinical sites: University of Pittsburgh (Pittsburgh, PA) and Wake Forest University School of Medicine (Winston-Salem, NC). Eligible participants were ≥70 years old at enrollment, had a body mass index (BMI) of 18–40 kg/m2, and were eligible for magnetic resonance (MR) imaging and a muscle tissue biopsy (Cummings et al. 2023). Individuals were further excluded if they had active cancer or were in the advanced stages of heart failure, renal failure on dialysis, dementia, or Parkinson’s disease. Participants must have been able to complete the 400-meter walk; those who appeared as they might not be able to complete the 400m walk at the in-person screening visit completed a short distance walk (4 meters) to ensure their walking speed as >=0.6m/s. The study protocol was approved by the Western Institutional Review Board Copernicus Group (WCG IRB; study number 20180764) and all participants provided written informed consent. In brief, baseline testing occurred across 3 separate days of clinic visits that were generally within 6-8 weeks of each other. The mean time between Day 1 and 3 was 42 days or ∼ 6 weeks. Day 1 included general clinic assessments (e.g., physical and cognitive tests; 5 hours), Day 2 included magnetic resonance imaging and Cardiopulmonary Exercise Testing (MR and CPET, 2–3 hours), and Day 3 included fasting specimen and tissue collection (2 hours). There were 879 participants who completed Day 1 of baseline testing and had at least one primary SOMMA measure: CPET, MR imaging, or muscle tissue biopsy.
Demographic, health, and functional measures
Cardiorespiratory fitness (VO2 peak)
Cardiorespiratory fitness was measured using gold standard VO2 peak (mL/min) from Cardiopulmonary Exercise Testing (CPET). A standardized CPET, using a modified Balke or manual protocol, was administered to participants to measure ventilatory gases, oxygen and carbon dioxide inhaled and exhaled during exercise (Balady et al. 2010). Two slow 5-minute walking tests were conducted before and after the maximal effort test to assess walking energetics at preferred walking speed and a slow fixed speed of 1.5 mph. Participants who were excluded from the maximal effort symptom-limited peak test had acute electrocardiogram (ECG) abnormalities, uncontrolled blood pressure or history of myocardial infarction, unstable angina or angioplasty in the preceding 6 months. Testing for VO2 peak began at the participant’s preferred walking speed with incremental rate (0.5 mph) and/or slope (2.5%) increased in 2-minute stages until respiratory exchange ratio, ratio between VCO2 and VO2, was ≥1.05 and self-reported Borg Rating of Perceived Exertion (Borg 1982) was ≥17. Blood pressure, pulse oximetry, and ECG were monitored throughout exercise. VO2 peak was determined in the BREEZESUITE software (MGC Diagnostics, St. Paul, MN) as the highest 30-second average of VO2 (L/min) achieved. The data were manually reviewed to ensure the correct VO2 peak was selected for each participant.
Other measures
Participants were asked to walk at their usual pace for 400 m from which walking speed (m/s) was calculated. Whole-body D3Cr muscle mass was measured in participants using a d3-creatine dilution protocol as previously described (Stimpson et al. 2012; Stimpson et al. 2013). Knee extensor leg power was assessed using a Keiser Air 420 exercise machine in the same leg as the muscle biopsy. Resistance to test power was based on determination of the 1 repetition maximum leg extensor strength. Weight was assessed by balance beam or digital scales and height by wall-mounted stadiometers. An approximately 6-minute-long MR scan was taken of the whole body to assess body composition including thigh muscle volume with image processing by AMRA Medical (Linge et al. 2018). The CHAMPS questionnaire (Stewart et al. 2001) was used to assess specific types and the context of physical activities. Participants were asked to self-report physician diagnosis of cancer (excluding nonmelanoma skin cancer), cardiac arrythmia, chronic kidney disease, chronic obstructive pulmonary disease, coronary artery disease, congestive heart failure, depression, diabetes, stroke, and aortic stenosis; from this list a multimorbidity count (0, 1, or 2+) was calculated.
Body Size
Weight was assessed by balance beam or digital scales and height by wall-mounted stadiometers.
Gene expression and mitochondrial respiration measurements
Skeletal muscle biopsy collection and processing
Percutaneous biopsies were collected from the middle region of the musculus vastus lateralis under local anesthesia using a Bergstrom canula with suction. The specimen was blotted dry of blood and interstitial fluid and dissected free of any connective tissue and intermuscular fat. Approximately 20 mg of the biopsy specimen was placed into ice-cold BIOPS media (10 mM Ca–EGTA buffer, 0.1 M free calcium, 20 mM imidazole, 20 mM taurine, 50 mM potassium 2-[N-morpholino]-ethanesulfonic acid, 0.5 mM dithiothreitol, 6.56 mM MgCl2, 5.77 mM ATP, and 15 mM phosphocreatine [PCr], pH 7.1) for respirometry, as previously described (Mau et al. 2023). Myofiber bundles of approximately 2–3 mg were teased apart using a pair of sharp tweezers and a small Petri dish containing ice-cold BIOPS media. After mechanical preparation, myofiber bundles were chemically permeabilized for 30 minutes with saponin (2 mL of 50 μg/mL saponin in ice-cold BIOPS solution) placed on ice and a rocker (25 rpm). Myofiber bundles were washed twice (10 min each) with ice-cold MiR05 media (0.5 mM ethylenediaminetetraacetic acid, 3 mM MgCl2·6H2O, 60 mM K-lactobionate, 20 mM taurine, 10 mM KH2PO4, 20 mM N-2-hydroxyethylpiperazine-NL-2-ethanesulfonic acid, 110 mM sucrose, and 1 g/L bovine serum albumin, pH 7.1) on an orbital shaker (25 rpm). The second wash in MiR05 contained blebbistatin (25 μM), a myosin II ATPase inhibitor, that was used to inhibit muscle contraction. Fiber bundle wet weight was determined immediately after permeabilization using an analytical balance (Mettler Toledo, Columbus, OH).
Mitochondrial respiration
Maximal complex I- and II-supported oxidative phosphorylation (P1-Max OXPHOS, also known as State 3 respiration) was measured in permeabilized muscle fiber bundles from biopsies as previously described (Mau et al. 2023).
RNA Library preparation and sequencing
Total RNA from frozen human skeletal muscle samples (∼5 to 30mg) was prepared using Trizol solution (Invitrogen) according to manufacturer’s direction in 2.0mL Eppendorf safe-lock tubes. Homogenization was performed using the Bullet Blender (NextAdvance, Raymertown NY USA) with an appropriate quantity of stainless-steel beads (autoclaved, 0.5∼2mm, NextAdvance, Raymertown NY USA) at 4°C on Setting 8 (Max is 12) in 30 second bouts. The homogenization step was repeated 5 times for a total of 3 minutes with at least 1 minute break between each bout. The removal of residual genomic DNA was performed by incubating the RNA sample with DNase (AM1907, Thermosci) plus RiboLock RNase inhibitor (EO0381, Thermisci) at 37°C for 25min in a heating block (400rpm). Cleanup of the RNA samples was done using the DNase inactivation reagent following instructions provided by the manufacturer (AM1907, Thermosci). The RNA concentration and integrity were determined by using ThermoSci Nanodrop and Agilent Tapestation.
To prepare RNAseq library, polyA mRNA was isolated from about 250ng total RNA by using NEBNext Poly(A) mRNA magnetic isolation module (E7490L, NEB) and mRNA library constructed by using NEBNext Ultra II directional RNA library Pre Kit for Illumina (E7760L, NEB). Equal molarity of RNAseq libraries were pooled and sequenced on Illumina NovaSeq (2X150bp) to reach 80M reads per sample.
Alignment and quality control
The reads from RNA-sequencing were aligned to the Genome Reference Consortium Human Build 38 (GRCh38) using the HISAT2 software(Kim et al. 2019). Duplicated aligned reads were further marked and excluded using the Picardtools software (http://broadinstitute.github.io/picard/). Expression count data were obtained using the HTseq software (Anders et al. 2015). Genes with a total count of ≤ 20 across all samples were filtered out to remove non-expressed genes. The quality of the RNASeq data was gauged by the alignment rate and the duplication rate.
Association analyses
Expression levels of 281 genes related to mTOR and upstream pathways (133 genes), autophagy regulation (68 genes), and selective autophagy of mitochondria (mitophagy) (80 genes) were analyzed (Bordi et al. 2021). Gene expression associations with traits were identified using negative binomial regression models as implemented by DESeq2 (Love et al. 2014) in R and adjusted for age, gender, clinic site, race/ethnicity, height, weight, hours per week in all exercise-related activities (CHAMPS), multimorbidity count category, and sequencing batch. DESeq2 uses a negative binomial generalized linear model for differential analysis and applies the Wald test for the significance of GLM coefficients. The Benjamini-Hochberg false discovery rate method was used for P-value adjustment. Genes were considered differentially expressed according to the significance criteria of FDR<0.05. In negative binomial models, traits were modeled using the number of standard deviations (SDs) from each trait’s mean. Consequently, the reported log base 2-fold changes reflect the change in gene expression per one SD unit increase in each trait. Volcano plots were created to visualize the differential expression of RNAs (ENSGs) associated with functional measures. Heat maps were created to summarize significant ENSG associations across all analyzed traits.
Statement relating to relevant ethics and integrity policies
The study protocol was approved by the Western Institutional Review Board Copernicus Group (WCG IRB; study number 20180764) and all participants provided written informed consent.
Conflict of Interest statement
S.R.C. is a consultant to Bioage Labs. P.M.Ca. is a consultant to and owns stock in MyoCorps. All other authors declare no conflict of interest. A.M.C is a consultant to Generian Therapeuthics and it is a consultant to and owns stock in Life Biosciences.
Funding statement
SOMMA is funded by the National Institute on Aging (NIA) grant number R01AG059416. Study infrastructure support was funded in part by NIA Claude D. Pepper Older American Independence Centers at University of Pittsburgh (P30AG024827) and Wake Forest University (P30AG021332) and the Clinical and Translational Science Institutes, funded by the National Center for Advancing Translational Science at Wake Forest University (UL10TR001420).
Authors’ contributions
Peggy M Cawthon, Paul M Coen, Steven R Cummings, Steven B Kritchevsky, Anne B Newman, Russell T Hepple, Karyn Esser and Gregory J Tranah designed the study. Paul M Coen and Gregory J Tranah drafted the manuscript. Haley N Barnes, Paul M Coen, Daniel S Evans, and Zhiguang Huo, conducted the genetic analyses. Paul M Coen, Steven R Cummings, Steven B Kritchevsky, Anne B Newman, Russell T Hepple, Karyn Esser, Olaya Santiago Fernandez and Ana Maria Cuervo and Gregory J Tranah contributed to the interpretation of the results. Peggy M Cawthon, Paul M Coen, Frederico GS Toledo, Steven B Kritchevsky, Anne B Newman, Steven R Cummings, Gregory J Tranah have critically revised the manuscript. All authors reviewed and approved the final version of the manuscript and Haley N Barnes and Zhiguang Huo had full access to the data in the study and accept responsibility to submit for publication.
Data availability statement
All SOMMA data are publicly available via a web portal. Updated datasets are released approximately every 6 months (https://www.sommastudy.com/for-investigators).
SUPPLEMENTARY INFORMATION
Listings of Ensembl gene (ENSG) associations identified by negative binomial regression for each trait (table 3: Max OXPHOS, table 4: VO2 peak, table 5: 400-meter Walk Speed, table 6: Leg Strength, table 7: Leg Power, table 8: Thigh Muscle Volume).
Acknowledgments
We acknowledge all SOMMA staff and investigators, and we thank all the SOMMA participants who enabled this research. For a full list of personnel who contributed to the SOMMA study, please see Cummings SR, Newman AB, Coen PM, Hepple RT, Collins R, Kennedy K, Danielson M, Peters K, Blackwell T, Johnson E, Mau T, Shankland EG, Lui LY, Patel S, Young D, Glynn NW, Strotmeyer ES, Esser KA, Marcinek DJ, Goodpaster BH, Kritchevsky S, Cawthon PM. The Study of Muscle, Mobility and Aging (SOMMA). A Unique Cohort Study about the Cellular Biology of Aging and Age-related Loss of Mobility. J Gerontol A Biol Sci Med Sci. 2023 Feb 9:glad052. doi: 10.1093/gerona/glad052. Epub ahead of print. PMID: 36754371.
Footnotes
Paul.coen{at}adventhealth.com, zhuo{at}ufl.edu, greg.tranah{at}ucsf.edu, Haley.Barnes{at}ucsf.edu, Peggy.Cawthon{at}ucsf.edu, rthepple{at}phhp.ufl.edu, toledofs{at}upmc.edu, Daniel.Evans{at}ucsf.edu, olaya.santiagofernandez{at}einsteinmed.edu, ana-maria.cuervo{at}einsteinmed.edu, skritche{at}wakehealth.edu, Steven.Cummings{at}ucsf.edu, anewmann{at}pitt.edu, kaesser{at}ufl.edu