A novel deep proteomic approach in human skeletal muscle unveils distinct molecular signatures affected by aging and resistance training

The skeletal muscle proteome alterations to aging and resistance training have been reported in prior studies. However, conventional proteomics in skeletal muscle typically yields wide protein abundance ranges that mask the detection of lowly expressed proteins. Thus, we adopted a novel deep proteomics approach whereby myofibril (MyoF) and non-MyoF fractions were separately subjected to protein corona nanoparticle complex formation prior to digestion and Liquid Chromatography Mass Spectrometry (LC-MS). Specifically, we investigated MyoF and non-MyoF proteomic profiles of the vastus lateralis muscle of younger (Y, 22±2 years old; n=5) and middle-aged participants (MA, 56±8 years old; n=6). Additionally, MA muscle was analyzed following eight weeks of resistance training (RT, 2d/week). Across all participants, the number of non-MyoF proteins detected averaged to be 5,645±266 (range: 4,888–5,987) and the number of MyoF proteins detected averaged to be 2,611±326 (range: 1,944–3,101). Differences in the non-MyoF (8.4%) and MyoF (2.5%) proteomes were evident between age cohorts, and most differentially expressed non-MyoF proteins (447/543) were more enriched in MA versus Y. Biological processes in the non-MyoF fraction were predicted to be operative in MA versus Y including increased cellular stress, mRNA splicing, translation elongation, and ubiquitin-mediated proteolysis. RT in MA participants only altered ~0.3% of MyoF and ~1.0% of non-MyoF proteomes. In summary, aging and RT predominantly affect non-contractile proteins in skeletal muscle. Additionally, marginal proteome adaptations with RT suggest more rigorous training may stimulate more robust effects or that RT, regardless of age, subtly alters basal state skeletal muscle protein abundances.


INTRODUCTION
Aging adversely affects skeletal muscle physiology as evidenced by a reduction in muscle stem (or satellite) cell content, a loss of myofibrillar protein, and a loss in motor units and myofibers [1][2][3][4].Muscle aging is also associated with impairments in mitochondrial function, dysfunctional redox balance, and altered proteostasis [5][6][7].The culmination of these events likely contributes to a loss in muscle mass, which according to a recent AGING review, is accelerated in all the body's musculature past the age of 50 years old [8].
Resistance training can reverse certain aspects of skeletal muscle aging.For instance, weeks to months of resistance training in older participants has been shown to increase tissue-level and myofiber hypertrophy and muscle satellite cell content [9].Resistance training also catalyzes skeletal muscle mitochondrial biogenesis and remodeling in older participants [10][11][12][13][14], and weeks to months of resistance training alters nuclear and mitochondrial DNA methylation patterns in older participants which may lead to "rejuvenating" effects on global mRNA expression patterns [11,15,16].
Proteomic investigations intend to examine the entire detectable protein expression signature of a given tissue under various experimental conditions.While other -omics-based approaches exist (i.e., genomics, epigenomics, transcriptomics, and metabolomics), it has been posited that proteomic signatures likely best translate to cellular and tissue phenotypes [17].Past proteomic investigations have provided tremendous insight as to how myofiber type, aging, and exercise training affect the skeletal muscle molecular milieu [1,14,[18][19][20][21][22].Notwithstanding, skeletal muscle-based proteomics poses technical challenges.For instance, skeletal muscle tissue processing with general lysis buffers results in the clearance of insoluble (e.g., contractile) proteins [23], and if standard bottom-up proteomics is employed thereafter, the relative abundances of these proteins will ultimately be inaccurate.Even if care is taken in isolating the poorly soluble contractile and soluble non-contractile protein fractions, another pitfall lies in lowly-abundant proteins being masked by highly abundant proteins in each fraction [1].Single fiber isolation techniques have increased the depth of proteins detected [21,24].However, certain disadvantages with this method exist including the burdensome process of tissue digestion and fiber dissection, the need for pooling myofibers to obtain adequate protein for proteomics, and the inability to detect proteins enriched in the extracellular matrix or stromal cells.
A novel deep proteomics approach in human plasma was recently published whereby unique nano-bio interaction properties of multiple magnetic nanoparticles (NPs) was leveraged for automated protein separation (referred to as the Proteograph assay; Seer, Inc. Redwood, CA, USA) [25].Downstream digestion followed by liquid chromatography coupled to mass spectrometry (LC-MS) analyses enabled the identification of over 2,000 plasma proteins and this provided approximately a 10fold increase in depth compared to prior studies that utilized other strategies to deplete plasma of highly abundant proteins [26,27].However, this approach has not been performed in human skeletal muscle.Thus, we sought to leverage this technology, along with our prior method of muscle tissue fractionation [23], to examine the proteomic signatures of the myofibrillar (MyoF) and non-myofibrillar (non-MyoF) fractions from the vastus lateralis (VL) muscle of a subset of younger (Y, 22±2 years old, n=5, all of whom were females) and middleaged (MA, 56±8 years old, n=6, 3 men and 3 women) participants.We also sought to determine how eight weeks of unilateral knee extensor resistance training affected the MyoF and non-MyoF proteomic signatures in the MA cohort.Given some of our past work in this area [1], we hypothesized that more non-MyoF proteins would be altered by aging when comparing MA and Y participants.We also hypothesized that resistance training in MA participants would affect more non-MyoF versus MyoF proteins.However, we did not adopt an a priori hypothesis regarding which proteins or biological processes would be affected between comparisons given the novelty of interrogating skeletal muscle using the Proteograph assay.

Characteristics of proteins identified in the MyoF and non-MyoF fractions
A total of 6,445 non-MyoF proteins and 4,421 MyoF proteins were identified in at least one participant (Figure 1A).Across all participants, the number of non-MyoF proteins detected averaged to be 5,645 ± 266 (range: 4,888-5,987) and the number of MyoF proteins detected averaged to be 2,611 ± 326 (range: 1,944-3,101).Notably, there was a numerically lower number of MyoF proteins detected in MA (preintervention) versus Y participants that approached statistical significance (p=0.066), and a non-significant greater number of non-MyoF proteins detected in MA (pre-intervention) versus Y participants (p=0.112; Figure 1C).

MYH isoform peptide identification information
Myosin heavy chain isoforms have been intensely studied in human skeletal muscle for fiber typing purposes and prominent isoforms include the slowtwitch type I isoform (encoded by the MYH7 gene) as well as the fast-twitch IIA (encoded by the MYH2 gene) and IIX (encoded by the MYH1 gene) isoforms [28].However, other MYH isoforms were highly enriched in the MyoF fraction according to data presented in Figure 1C.Because of this, we opted to provide the peptide sequences used for detecting some of these isoforms in Table 1 below.

Alternative proteins isoforms in the MyoF and non-MyoF fractions of MA and Y
The enhanced depth of detection provided by proteomics revealed the presence of numerous isoforms in both protein fractions; specifically, there were 175 isoforms for 82 MyoF proteins and 375 isoforms for 173 non-MyoF proteins.The MyoF proteins with the most isoforms included titin (TTN), myosin-binding protein C (MYBPC), and MICOS complex subunit MIC60 (IMMT); each of these targets had four isoforms detected.Given the vast research interest in titin [29,30], associated isoform data are plotted in Figure 3A; notably no significant aging or training effects were noted (p>0.01 for all comparisons).The non-MyoF proteins with the most isoforms included Gelsolin (GSN, four isoforms), IMMT (4 isoforms), and Reticulon-4 (4 isoforms); again, no significant age effects were noted for these targets (data not plotted).
We plotted significantly different alternative protein isoform abundances between Y and MA in both protein fractions given that bioinformatics on the non-MyoF fraction indicated "regulation of RNA splicing" (GO:0043484) was predicted to be upregulated in the older cohort (Figure 3B, 3C).MA (pre-intervention) Legend: these data contain the peptide sequences used for alignment to identify the several highly abundant myosin heavy chain isoforms in the MyoF fraction.
versus Y comparisons indicated that only three alternative MyoF protein isoforms were different between age groups (all higher in Y, p<0.01).However, 14 alternative non-MyoF protein isoforms were significantly different between age groups (11 higher in MA, p<0.01), and while not depicted in Figure 3C, 32 additional alternative non-MyoF protein isoforms were numerically different between age groups (25 higher in MA, p<0.05).Further, when examining the abundances of small nuclear ribonucleoproteins belonging to    spliceosome complexes in the non-MyoF fraction, three reached the p<0.01 significance threshold as being more enriched in MA (pre-intervention) versus Y participants (SNRP200, SNRPE, SNRPF), and several others were numerically greater in MA participants (SNRP40/70/A/A1/B/C/D2/D3, p<0.05; Figure 3D).Notably, training did not alter the expression of any SNRP in Figure 3D (p>0.05 for all), and most of these proteins were not detected in the MyoF fraction.

MyoF and non-MyoF protein differences prior to and following resistance training in MA participants
In MA, knee extensor resistance training significantly altered 13 MyoF proteins (11 upregulated and two downregulated, p<0.01; Figure 4A), and 64 non-MyoF proteins (56 upregulated and eight downregulated, p<0.01; Figure 4B, 4C).These alterations represented ~0.3% of the MyoF proteome and ~1.0% of the non-MyoF proteome.Bioinformatics within each fraction were attempted, albeit no pathways were predicted to be significantly affected.

Proteolysis targets manually interrogated in both fractions
Based on bioinformatics indicating that proteostasis was predicted to be altered with aging (Table 3), we manually interrogated proteolysis-related protein targets (i.e., calpain-1/2, and the summed spectra of 26S proteasome subunits).Figure 5

DISCUSSION
Using a novel analytical approach, we examined the deep proteomic signatures of the MyoF and non-MyoF fractions in younger adults as well as middleaged participants before and after eight weeks of knee extensor resistance training.More non-MyoF proteins differed between age cohorts compared to MyoF proteins (8.4% versus 2.5% of the respective protein pools).More non-MyoF proteins (447/543) were also more highly abundant in MA versus Y and bioinformatics predicted that several biological processes were more operative in the older participants.A greater abundance in alternative variants, proteins associated with spliceosomes, and proteolysis-related proteins were also evident in the non-MyoF fraction of MA versus Y, and these observations corroborated certain bioinformatics findings.Although resistance training in MA non-significantly increased VL crosssectional area (+6.5%) and significantly increased knee extensor strength (+8.7%), training marginally affected the MyoF and non-MyoF proteomes and no biological processes were predicted to be affected in either fraction.These findings will be expanded upon in the paragraphs below.
As stated above, several studies have performed proteomic analyses on skeletal muscle to compare molecular signatures that exist between younger and older adults or to examine how resistance training affects this aspect of the muscle-molecular milieu [1, 14,[18][19][20]22].The novelty of the current study was the proteomics approach utilized and the knowledge gained relative to these prior investigations.
Fractionation of muscle into solubilized MyoF and non-MyoF homogenates enabled the detection of unique proteins in each fraction, which has only been attempted in one other study to our knowledge [1].In this prior study, we performed bottom-up LC-MS-based proteomics on each fraction from younger resistancetrained, younger untrained, and older untrained men (n=6 per group).We identified a total of 810 proteins in both fractions that were expressed in at least one participant.In the current study we detected a total of 10,866 proteins in both fractions.While most of identified proteins were present in both fractions, we were able to identify 2,217 unique non-MyoF proteins and 193 unique MyoF proteins.This robust increase in detection depth (~13.4-fold) is insightful for numerous reasons.First, it was revealed that metabolic enzymes constituted the top class of proteins in both fractions as well as proteins that overlapped in both fractions.Hence, although myosin heavy chain isoforms, troponin, titin, and actin were the most enriched in the MyoF fraction, these data counter the notion that the MyoF fraction contains mainly contractile proteins.Several nuclear proteins were also identified in the MyoF fraction (e.g., histones and other chromatinbinding proteins) indicating that our MyoF isolation method likely pellets nuclei.Finally, we were able to identify numerous proteins that are not commonly reported in previous skeletal muscle proteomic studies (see Table 1 and Figure 2 for example).To this end, several MYH isoforms beyond the three common 7/2/1 isoforms were highly enriched in the MyoF fraction, and this may be due to the persistence of non-conventional or developmental isoforms in certain Notes: Certain p-values not meeting the significance criteria for these data were presented due to the visual differences observed between cohorts.AGING regions of adult myofibers or in transitioning fibers as discussed by Schiaffino et al. [31].Both fractions contained most of the large (~40) and small (~30) ribosomal subunit proteins, mitochondrial oxidative phosphorylation enzymes, and mitochondrial ribosomal proteins (>70), all of which likely represents the presence of sarcoplasmic and intermyofibrillar mitochondria.Both fractions also contained numerous alternative isoforms for several proteins, various transcription factors (e.g., MEF2D, F-box proteins, SMAD1/2/3, NFAT isoforms, and several others), DNA and RNA polymerase subunits, various growth factors and their receptors (e.g., EGFR, VEGFA, FGF2/13, PDGFRa/b, TGFB1/2, and others), dozens of eukaryotic initiation/ elongation factors, nearly 50 solute carrier family member proteins (i.e., nutrient and metabolite transporters), cytoplasmic and mitochondrial aminoacyl tRNA ligases, and hundreds of signal transduction proteins of interest to skeletal muscle biologists (e.g., mTOR, RPTOR, p70s6k, AMPK subunits, cyclindependent kinases and inhibitors, and others).We believe this enhanced level of detection was due to muscle fractionation, and more importantly, the utilization of NPs prior to mass spectrometry.Indeed, this same contention has been posited by others using this technology to increase detection depth of circulating proteins in human plasma [25,32].
Notable MyoF and non-MyoF proteome signatures between age groups were also evident.For instance, aging seemingly affects the non-MyoF protein signature more so than the MyoF fraction.This finding agrees with our past proteomic study where we reported that 37 non-MyoF proteins (versus only 18 MyoF proteins) were differentially expressed between college-aged and older men (mean age 62 years old) [1].However, the increased detection depth in the current study indicated that 112 of the 4,421 identified MyoF proteins met the p<0.01 significance threshold between age cohorts.Moreover, all but one of these proteins (TMPO, Lamina-associated polypeptide 2 isoforms beta/gamma) were greater in Y versus MA indicating either a loss or decreased expression of ~2.5% proteins belonging to the MyoF fraction.Although no biological processes were predicted to be affected between age cohorts based on this list of 111 proteins, several of these targets were notable.For example, three heat shock proteins (HSPB1, HSP90AB1, HSPD1) were more lowly abundant in MA participants and this agrees in principle with past literature indicating the expression of heat shock proteins in skeletal muscle is dysregulated with aging [33,34].Several MyoF mitochondrial proteins were also lower in MA participants (DECR1, NDUFS8, ETFDH, GPD2, STOML2, ALDH2, COQ5, IARS2, ALDH1B1, MRPL32).This also agrees with past literature indicating either a decrease in mitochondrial volume density or decreased mitochondrial function with skeletal muscle aging [35], and more specifically agrees with a study by Callahan et al. [36] who reported that older participants presented reductions in the size of intermyofibrillar mitochondria.
Strikingly, 543 proteins (~8.4%) of the non-MyoF proteome were different between MA and Y participants, and unlike the trends observed in the MyoF fraction, most of these proteins (447) were significantly enriched in the MA cohort.These figures agree with a proteomic investigation by Robinson et al. [14] who reported that more muscle proteins (220/347) were higher in untrained older versus younger individuals prior to a period of exercise training.These robust differences between age cohorts in the current study also revealed that several biological processes were predicted to be upregulated in MA participants.Some of these processes either contradict each other or agree with past literature reporting similar aging phenotypes.
Regarding the former, while more proteins associated with translation elongation were more abundant in MA participants (which would potentially promote muscle anabolism), proteins associated with ubiquitin-mediated proteolysis were also more enriched.One interpretation of these data could be that aging increases skeletal muscle protein turnover.However, this is likely not the case given that a variety of studies ranging from human tracer studies to nematode models have indicated that protein turnover in response to feeding or in a basal state is impaired with aging [37,38].This aspect of our data also agrees with a report by Ubaida-Mohien et al. [22] who showed that 31% of proteins related to proteostasis were altered with age (24 underrepresented and 50 overrepresented, p<0.05) in healthy older versus younger adults.Hence, we posit that proteins associated with these processes may have been more abundant in MA participants in a compensatory attempt to counter age-related declines in muscle protein turnover.It is also interesting that proteasome subunits and calpains were more enriched in the MyoF fraction of Y versus MA, whereas these same proteins were more enriched in the non-MyoF MA versus Y.While speculative, an enrichment of proteolytic proteins in the MyoF fraction might play a role in functional proteostasis, while an enrichment of these proteins in the non-MyoF fraction might be indicative of a gradual dysregulation in proteostasis.
Non-MyoF proteins associated with mRNA export and splicing were also elevated in MA participants, and this agrees with other reports.For instance, a recent review by Park et al. [39] cites a variety of cell culture evidence to suggest that the nuclear pore complex is disrupted with aging and that this leads to a dysregulation in mRNA export.A rodent study by Mobley et al. [40] AGING also suggests that mRNA levels linearly decrease in skeletal muscle with increasing age.Hence, again, a higher abundance in non-MyoF proteins associated with mRNA export could also be a compensatory response to offset these age-associated effects burdened by myonuclei.The greater abundance of proteins associated with mRNA splicing is striking and agrees in principle with a report by Rodriguez et al. [41] who showed that the skeletal muscle of aged mice possessed ~4 times more RNA splice variants than younger counterparts.Our data also agree with the abovementioned proteomics report by Ubaida-Mohien et al. [22] who showed that proteins related to alternative splicing were more abundant in healthy older (versus younger) adults.We also performed a follow-up analysis showing that several small nuclear ribonucleoproteins (SNRPs, or snRNPs) that make up spliceosomes and alternative non-MyoF protein isoforms (indicative of increased spliceosome activity) were elevated in MA versus Y participants.This is particularly insightful given that dysfunctional spliceosome activity and the aberrant RNA and protein expression of splice variants have been linked to age-associated maladies such as cellular senescence [42,43].Hence, these independent reports of ageassociated increases in splice RNA and protein variants, along with the current data suggesting that the relative abundances of spliceosomes are greater in older participants, warrant future research elucidating the causes and consequences of this phenomenon.
A final noteworthy topic was the observation of marginal MyoF and non-MyoF alterations with eight weeks of knee extensor training in MA participants.Although this may have been due to the modest training regimen that only lasted eight weeks in duration, limited proteome plasticity with aging and/or the limited ability of resistance training to alter the muscle proteome could also be plausible explanations.Support for both phenomena come from Robinson et al. [14] who employed deep proteomics to report that ~200 muscle proteins were altered in older participants after 12 weeks of resistance training (p<0.05), and this was less than the ~300 proteins that were altered in the younger participants.Hence, an aging effect was noted.However, assuming the authors identified >3,000 muscle proteins, which was not reported to our knowledge, this represents less than 10% of the detectable proteome being altered with resistance training.Deane et al. [13] used a different proteomics approach to examine the non-MyoF proteome adaptations in older participants following 20 weeks of resistance training in younger and older adults.Although their depth of detection was limited to ~160 proteins, resistance training only increased five non-MyoF proteins in older participants when a p<0.05 significance threshold was employed (i.e., ~3% of the detectable proteome).Interestingly, this effect was not confined to older participants given that the younger participants in their study only presented an elevation in four non-MyoF proteins with training.Our laboratory also used proteomics to examine non-MyoF protein adaptations in college-aged men following ten weeks of resistance training [1].Only 13 proteins were shown to be altered with training (12 up, one down, p<0.05) and this represented ~3.4% of the detectable non-MyoF proteome.Hence, these two latter studies do not support the aging hypothesis and, instead, provide evidence of limited muscle proteome plasticity with resistance training.Despite marginal non-MyoF proteome alterations in MA participants with training, there were interesting targets that were altered.For instance, the UBR7 E3 ligase was upregulated, and recent evidence suggests that an E3 ligase in this same protein family (UBR5) is required for load-induced skeletal muscle hypertrophy [44].Additionally, the knockdown of another member of this family (UBR4) promotes hypertrophy in Drosophila and mice [45].SRC was upregulated and this non-receptor tyrosine kinase has been implicated in interacting with vitamin D to promote anabolic signaling in skeletal muscle [46].HDAC4 was upregulated and this adaptation could be operative in muscle-metabolic adaptations and ultrastructural remodeling to resistance training given that HDACs have been implicated in controlling the expression of various metabolic and contractile protein genes [47,48].Two protein phosphatases were also upregulated (PPP1R12A and PIP4P2), and both have been shown to be involved in aspects of insulin and growth factor signaling [49,50].
There are notable limitations to the current study.First, the MA and Y participant pools were small, and Y participants were all women.Furthermore, the lack of training data in Y participants to examine differential age-associated training effects is a limitation.We lacked remaining skeletal muscle to perform downstream analyses (e.g., examining RNA splice variants and/or proteasome activity assays) which may have provided additional insight.Moreover, having additional age groups (e.g., older participants) would have provided more depth to the current study.Despite these limitations, the primary objective of this publication was to feature our novel proteomic approach as we believe that this will add tremendous insight into the field of skeletal muscle biology.Likewise, the data from this study can be used to generate hypotheses for other age-related or resistance training proteomic or targeted protein approaches moving forward.
In conclusion, we provide preliminary evidence to support that muscle aging predominantly affects the non-MyoF protein pool and that this is associated with AGING biological processes which may act to counteract dysfunctional cellular homeostasis.We also provide preliminary evidence of limited MyoF and non-MyoF proteome plasticity to shorter-term resistance training in middle-aged participants, and this agrees with prior proteomic investigations.Finally, and most importantly, we believe that the utilization of skeletal muscle tissue fractionation protocols and NP-based protein corona formation prior to downstream proteomics has the potential to add incredible insight in identifying novel protein targets affected by exercise training, aging, and various disease states.

Ethical approval and study design
Muscle specimens were obtained from two studies whereby approval was obtained from the Auburn University Institutional Review Board.The first protocol in untrained MA participants (approved protocol #21-461 MR 2110) involved investigating the effects of a dietary supplement (312 mg of combined Wasabia japonica extract, theacrine, and copper (I) niacin chelate) versus a placebo on potential blood marker responses over an eight-week period.A unilateral leg resistance training (two days/week) protocol was implemented to perform non-supplementation secondary analyses as presented herein.The six MA participants included in the current study were in the placebo group; thus, no confounding effects of dietary supplementation were expected.Additionally, both the pre-intervention and post-intervention biopsies were collected the same time of day following at least a four-hour fast, and the post-intervention muscle biopsy was collected 72 hours following the last training bout.Y participant muscle tissue was banked from a prior study examining how ten weeks of daily peanut protein supplementation affected resistance training outcomes in untrained individuals (approved protocol #19-249 MR 1907) [51].Notably, muscle tissue from these participants was collected in the basal state prior to the intervention.Hence, again, there were no potential confounding effects of supplementation.Study procedures for both projects were in accordance with the most recent revisions of the Declaration of Helsinki except for the MA study not being pre-registered as a clinical trial.

Knee extensor resistance training
The resistance training intervention consisted of supervised unilateral leg extensions (two days/week for eight weeks), and the intervention was preceded and followed by strength and VL muscle assessments (described in later paragraphs).All MA participants trained their right legs whereby each training session consisted of five sets of 12 repetitions.The beginning training load was established at ~40% of the participants' three-repetition maximum (3RM).After each set, participants verbally articulated their perceived repetitions in reserve (RIR) [52], and training load was adjusted accordingly.RIR values of 0-2 after a set resulted in no training load change in each session.RIR values of 3-5 for consecutive sets resulted in the training load being increased by 5-10%.For RIR values ≥6 after one set, the training load was increased by 10-20%.If the weight could not be performed with full range of motion, or the participant could not complete 12 repetitions for a given set, the training load was decreased accordingly.

Strength testing
The first and last workout of the eight-week training paradigm consisted of maximal leg extensorflexion torque assessments using isokinetic dynamometry (Biodex System 4; Biodex Medical Systems, Inc., Shirley, NY, USA) and 3RM leg extensor strength testing.Prior to dynamometer testing, the participant's lateral epicondyle was aligned with the axis of the dynamometer's lever arm, and the hip was positioned at 90°.The participant's shoulders, hips, and leg were strapped and secured for isolation during testing.Following three warm-up trials at a submaximal effort, participants completed five maximal voluntary isokinetic knee extension and flexion actions at 60 degrees/second.Participants were provided verbal encouragement during each contraction.The isokinetic contraction resulting in the greatest peak torque value was used for analyses.Approximately five minutes following isokinetic dynamometry testing, participants performed 3RM strength testing using a free-weight apparatus.Prior to testing, participants were given a warm-up load and instructed to complete 10 repetitions.After participants recorded their RIR for the warmup set, the weight was adjusted accordingly for another warm-up set of five repetitions.RIR was recorded again to determine the participants starting load for a 3RM attempt.The load was incrementally increased 5-10% per 3RM attempt until 3RM testing concluded, indicated by failure of full range of motion on any of the repetitions, or if RIR recorded was 0. Participants were allowed a full three minutes of recovery between attempts.The isokinetic dynamometry and 3RM testing described was similar for both the first and final workout.

Urine specific gravity testing for hydration
Participants performed a testing battery prior to the start of training (PRE) and 3-5 days following the last resistance training workout (POST).Participants arrived AGING for testing at a minimum of 4 hours fasted and well hydrated.Upon arrival participants submitted a urine sample (~5 mL) for urine specific gravity assessment (USG).Measurements were performed using a handheld refractometer (ATAGO; Bellevue, WA, USA).USG levels in all participants were ≤ 1.020, indicating sufficient hydration [53].

Body composition testing
Body composition was assessed using multi-frequency bioelectrical impedance analysis (InBody 520, Biospace, Inc., Seoul, Korea).From the scan, body fat percentage was recorded.Previously determined test-retest reliability yielded an intraclass correlation coefficient (ICC3,1) of 0.99, standard error of the measurement (SEM) of 0.87%, and minimal difference (MD) of 1.71% for body fat percentage.

Ultrasonography assessment for muscle morphology
A detailed description of VL assessments using ultrasonography has been published previously by our laboratory [54,55].Briefly, real-time B-mode ultrasonography (NextGen LOGIQe R8, GE Healthcare; Chicago, IL, USA) using a multifrequency lineararray transducer (L4-12T, 4-12 MHz, GE Healthcare) was used to capture VL muscle cross-sectional area (mCSA).Prior to scans, the mid-thigh location was determined by measuring the total distance from the mid-inguinal crease in a straight line to the proximal patella, with the knee and hip flexed at 90°, a mark was made using a permanent marker at 50% of the total length.From that location, a permanent marker was used transversely to mark the mid-belly of the VL.This marking is where all pre-intervention ultrasound images were taken as well as the muscle biopsy (described below).All post-intervention images were taken at the pre-intervention biopsy scar to ensure location consistency between scans.During mCSA scans, a flexible, semirigid pad was placed around the thigh and secured with an adjustable strap to allow the probe to move in the transverse plane.Using the panoramic function of the device (LogicView, GE Healthcare), images were captured starting at the lateral aspect of the VL and moving medially until rectus femoris was visualized, crossing the marked location.All ultrasound settings were held constant across participants and laboratory visits (frequency: 10 MHz, gain: 50 dB, dynamic range: 75), and scan depth was noted and held constant across time points per participant.Images were downloaded and analyzed offline using ImageJ software (National Institutes of Health, Bethesda, MD, USA).All ultrasound images were captured and analyzed by the same investigators at each timepoint.Previously determined test-retest reliability on 10 participants measured twice within 24 hours (where BAR captured images and JSG analyzed images) yielded an intraclass correlation of 0.99 and standard error of measurement of 0.60 cm 2 .

Collection of muscle tissue
Muscle biopsies from all participants were obtained from the mid-belly of the right VL, and sampling time of day was standardized for MA participants at pre and post resistance training intervention.Lidocaine (1%, 1.0 mL) was injected subcutaneously above the skeletal muscle fascia at the previously marked location.After five minutes of allowing the anesthetic to take effect, a small pilot incision was made using a sterile Surgical Blade No. 11 (AD Surgical; Sunnyvale, CA, USA), and the 5-gauge biopsy needle was inserted into the pilot incision ~1 cm below the fascia.Approximately 30-50 mg of skeletal muscle was removed using a double chop method and applied suction.Following biopsies, tissue was rapidly teased of blood and connective tissue, placed in pre-labeled foils, flash frozen in liquid nitrogen, and subsequently stored at −80° C until processing described below.

MyoF and non-MyoF protein fractionation
The MyoF and non-MyoF protein fractions were isolated per methods published by our laboratory and others [23,56].On the day of homogenization, muscle tissue was powdered on a liquid nitrogencooled ceramic mortar and pestle.Approximately 30 mg of tissue was homogenized using tight-fitting pestles in 500 µL of 25 mM Tris, pH 7.2, 0.5% Triton X-100, with added protease inhibitors (Promega, cat# G6521; Madison, WI, USA).Samples were centrifuged at 1,500 g for 10 minutes at 4° C, supernatants (non-MyoF fraction) were transferred to new 1.7 mL tubes, and tubes were stored at -80ºC until shipment on dry ice to Seer, Inc. Remaining MyoF pellets were kept on ice and thoroughly aspirated with micro-pipet tips to remove residual supernatant.Thereafter, 300 μL of solubilization buffer was added which contained 20 mM Tris-HCl, pH 7.2, 100 mM KCl, 20% glycerol, 1 mM DTT, 50 mM spermidine with added protease inhibitors (Promega, cat# G6521).Samples were then homogenized using tight-fitting pestles and stored at -80ºC until shipment on dry ice to Seer, Inc.

Proteograph assay
Proteomics analysis was performed at Seer, Inc.
(Redwood City, CA, USA).For each sample, 250 μL of received sample was subjected to the Seer Proteograph Assay protocol.After loading samples onto the SP100 Automation Instrument, protein corona formation and processing was initiated to generate desalted purified AGING peptides for protein identification using Reversed Phase (RP) LC-MS.To form the protein corona, Seer's proprietary NPs were mixed with the samples and incubated at 37° C for 1 hour.Unbound proteins were removed prior to downstream wash, reduction, alkylation, and protein digestion steps which were performed according to Seer's Proteograph Assay protocol [25].

LC-MS configuration
Peptides obtained from each of the five NP mixtures were separately reconstituted according in a solution of 0.1% formic acid and 3% acetonitrile [57] spiked with 5 fmol μL PepCalMix from SCIEX (Framingham, MA, USA).Reconstitution volumes varied by NP types to allow for constant peptide quantity for MS injection between samples regardless of starting volume (240 ng: NP1, 400 ng: NP2, 360 ng: NP3, 120 ng: NP4, and 320 ng: NP5). 4 µL of each sample were analyzed with a Ultimate3000 RLSCnano LC system coupled with a Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher Scientific; Waltham, MA, USA).Peptides were loaded on an Acclaim PepMap 100 C18 (0.3 mm ID × 5 mm) trap column and then separated on a 50 cm μPAC analytical column (PharmaFluidics, Zwijnaarde, Belgium) at a flow rate of 1 μL/min using a gradient of 5-25% solvent B (100% ACN) mixed into solvent A (100% water) over 26 minutes.The mass spectrometer was operated in Data Independent Acquisition (DIA) mode using 10 m/z isolation windows from 380-1200 m/z and 3-second cycle time.MS1 scans were acquired at 60k resolution and MS2 at 30k resolution.

Data processing
DIA LC-MS data were processed using Proteograph Analysis Suite (PAS) v2.1 (Seer, Inc) using the DIA-NN search engine (version 1.8.1) in library-free mode searching MS/MS spectra against an in silico predicted library based on Uniprot's Homo Sapiens reference database (UP000005640_9606, download December 9, 2022).Library-free search parameters included trypsin digestion allowing for one missed cleavage, N-terminal methionine excision, fixed modification of cysteine carbamidomethylation, peptide length of 7-30 amino acids, precursor range of 300-1800 m/z, and fragment ion range of 200-1800 m/z.Heuristic protein inference was enabled, MS1 and MS2 mass accuracy was set to 10 ppm.Precursor FDR was set to 0.01, and PG q-value was set to 0.01.Quantification was performed on summed abundances of all unique peptides considering only precursors passing the q-value cutoff.PAS summarizes all NP values for a single protein into a single quantitative value.Specifically, a single protein may have been measured up to five times, once for each nanoparticle.To derive the single measurement value, PAS uses a maximum representation approach, whereby the single quantification value for a particular peptide or protein group represents the quantitation value of the NP which most frequently has measured any given proteins across all samples.
The relative abundances of protein targets were obtained by normalizing raw spectra values for each identified protein to total spectra within-subject.After normalization, undetected protein abundance values were set at zero.Protein values are presented as spectra-normalized values in all figures and results.MyoF and non-MyoF proteome data can be found in Supplementary Tables 1, 2.

Statistics and bioinformatics
Data processing and statistical analysis were performed

A total of 4 ,
228 proteins overlapped in both fractions yielding 2,217 unique non-MyoF proteins, 193 unique MyoF proteins, and 6,638 unique proteins identified.Using the PANTHER Classification System classifications, the top five protein classes of MyoF proteins, non-MyoF proteins, and proteins in both fractions are presented in Figure 1B.The top 15 enriched MyoF and non-MyoF proteins in MA (preintervention) and Y are presented in Figure 1D, 1E.None of the 15 MyoF or non-MyoF proteins met the p<0.01 significance criteria between age cohorts.

Figure 3 .
Figure 3. MyoF and non-MyoF alternative protein isoform differences detected with proteomics.Legend: Data presented for Y and MA (pre-and post-intervention) include the identified titin isoforms in the MyoF fraction (A), alternative MyoF protein isoforms affected by aging (B), and alternative non-MyoF protein isoforms affected by aging (C), and small nuclear ribonucleoproteins that make up spliceosomes between cohorts (D).Data are presented as mean ± standard deviations for individual protein spectra values (normalized to total run spectra values) and y-axes were scaled as log10 for improved visualization.Symbols: #, indicates lower in MA versus Y at one or both time points (p<0.01);*, indicates greater in MA versus Y at one or both time points (p<0.01);Φ, indicates greater in MA versus Y at one or both time points for panel d only (p<0.05).Notes: (ND), indicates that the isoform number was not provided from the Uniprot's Homo Sapiens reference database (UP000005640_9606).

Figure 4 .
Figure 4. MyoF and non-MyoF proteins altered with resistance training in MA participants.Legend: Data presented for MA prior to and following eight weeks of knee extensor training include proteins in the MyoF fraction (11 up-regulated, 2 down-regulated; (A), the top 15 up-regulated proteins in the non-MyoF fraction (B), and all 8 down-regulated proteins in the non-MyoF fraction (C).Data are presented as mean ± standard deviations for individual protein spectra values (normalized to total run spectra values), and y-axes were scaled as log10 for improved visualization.Notes: No biological processes were predicted to be affected with training based on these alterations.

Figure 5 .
Figure 5. MyoF and non-MyoF proteasome and calpain proteins.Legend: Data presented for Y and MA prior to and following eight weeks of knee extensor training include proteasome subunits and calpains 1/2 in the MyoF fraction (A-C) and non-MyoF fraction (D-F).Data are presented as mean ± standard deviations for individual protein spectra values (normalized to total run spectra values).Symbols: #, indicates lower in MA versus Y at one or both time points (p<0.01);*, indicates greater in MA versus Y at one or both time points (p<0.01).Notes: Certain p-values not meeting the significance criteria for these data were presented due to the visual differences observed between cohorts.

Table 2 . Top 15 of 111 MyoF proteins greater in Y versus MA participants.
Legend: Data are presented as mean ± standard deviations for individual protein spectra values (normalized to total run spectra values) in younger (Y) and pre-intervention middle-aged (MA) participants.Protein targets are sorted from most to least abundant in the Y cohort.

Table 3 . Top 30 of 543 non-MyoF proteins different between Y versus MA participants.
Legend: Data are presented as mean ± standard deviations for individual protein spectra values (normalized to total run spectra values) in younger (Y) and pre-intervention middle-aged (MA) participants.Protein targets more abundant in Y versus MA are sorted from most to least abundant in the Y cohort.Protein targets more abundant in MA versus Y are sorted from most to least abundant in the MA cohort.

Table 4 . Biological processes affected based on non-MyoF proteins different between Y versus MA participants.
Legend: Biological pathways presented to be up-and down-regulated between age groups based on differential non-MyoF protein expression differences between age cohorts (i.e., Table3data).
). AGING protein.In the non-MyoF fraction, the summed spectra of the 42 detected 26S proteasome subunits was not significantly different between Y versus MA participants at either time point.However, both calpains were higher in the non-MyoF fraction of MA at both time points versus Y participants (CAPN1 was significant in MA pre-intervention versus Y, p<0.01), and training did not significantly affect either protein.
[58,59]icrosoft Excel for Microsoft 365 (Redmond, WA, USA) and GraphPad Prism version 9.2.0 (San Diego, CA, USA).Independent samples t-tests were used for Y versus MA (pre-intervention) to determine age effects, and dependent samples t-tests were used to determine training effects in MA.All data in tables and figures are presented as mean ± standard deviation (SD) values.Training phenotypes were considered significantly different at p<0.05, although approaching values (i.e., p<0.100) were discussed as "numerical" changes due to limited n-sizes.Conversely, significant aging and training effects for protein targets were established as p<0.01 for enhanced stringency given the high number of identified proteins, although again approaching values (i.e., p<0.05) were discussed in certain circumstances due to limited n-sizes.Bioinformatics was performed using PANTHER v17.0[58,59].First, protein lists from each fraction were characterized using the functional classification tool.Next, overrepresentation tests of PANTHER GO-Slim biological processes were performed between Y and MA participants and in MA participants from pre-to-post training.Parameters for statistical overrepresentation tests included the following: i) entered proteins had to meet the aforementioned p<0.01 significance threshold, ii) protein lists were entered separately based on being up-or downregulated to generate a list of biological processes that were predicted to be directionally affected, and iii) Fisher tests with Bonferroni adjusted p<0.05 values were used as significance thresholds.variousaspects of proteomics analysis.C. B. Mobley, T. N. Ziegenfuss, A. D. Fruge, and A. N. Kavazis were intellectually involved in various stages of the project.All co-authors critically contributed to the writing of this manuscript, and all co-authors approved the final submitted version.