衰老同步削弱男性成年人骨骼肌蛋白质组可塑性与抗阻训练肥大反应

《Physiological Reports》：Aging concomitantly reduces skeletal muscle proteome plasticity and hypertrophic responses to resistance training in male adults

　　

抗阻训练（resistance training，RT）诱导的肌肉肥大随衰老而减弱，但这一减弱现象的蛋白质组学基础尚不清楚。研究人员因此探究骨骼肌蛋白质组可塑性的降低是否伴随与年龄相关的肥大反应迟钝。17名年轻（21.9±2.5岁）和15名老年（57.5±6.9岁）未经训练的男性完成了10–12周RT，并在训练前后采用数据非依赖性采集质谱法（data-independent acquisition mass spectrometry，DIA-MS）对股外侧肌（vastus lateralis，VL）可溶性蛋白质组进行 profiling（共2106种蛋白质）。尽管训练量相似（p=0.771），年轻人瘦体重增加更多（+2.73 vs. +1.18 kg，p=0.009），VL厚度增加也更多（+0.49 vs. +0.20 cm，p=0.002）。蛋白质组可塑性与此差异平行：RT改变了年轻人中135种蛋白质（FDR<0.05），而老年人中无一改变。这种减弱是幅度而非方向性的，因为年轻人中135种训练响应蛋白的81.5%在老年人中呈同向变化（ρ=0.711），且老年人中有75个基因集在RT后显著富集，以翻译机器（translational machinery）为首（NES最高达+2.85）。基线时，278种蛋白质在年龄组间存在差异（FDR<0.05），且年龄主导了共表达网络（Turquoise模块AUC=0.96）。总之，老年男性迟钝的肥大反应与RT诱导的蛋白质组可塑性丧失相伴随，反映的是压缩而非质变的RT适应。

论文解读文章

一、研究背景与问题

维持或增加骨骼肌质量对中老年人降低跌倒致命伤害风险、预防心血管代谢疾病及提升生活质量至关重要。骨骼肌质量通常在30–40岁达到峰值，此后进行性下降，60岁后流失加速。这种与年龄相关的肌肉流失（即肌少症，sarcopenia）影响高达45%的70岁以上成年人，每年造成约404亿美元的医疗费用负担。尽管抗阻训练（RT）是维持或增加肌肉质量与力量的可行策略，但老年人的肥大反应较年轻人减弱，表现为同等训练后瘦体重增幅减小，并伴随选择性II型纤维萎缩，后者构成年龄相关肌肉流失的主体。这种力量保持（部分通过增强神经适应实现）与肥大减弱（即合成代谢抵抗，anabolic resistance）之间的分离，构成以运动为基础的肌少症预防干预的根本障碍。

在分子层面，多种机制促成老年人减弱的肥大反应，其中许多反映更广泛的蛋白质稳态（proteostasis）丧失——这是生物衰老的主要标志之一。衰老可能钝化mTORC1通路对机械负荷和氨基酸可利用性的激活反应，降低肌肉蛋白质合成速率；自噬-溶酶体和泛素-蛋白酶体系统的年龄相关缺陷损害蛋白质质量控制并促进损伤蛋白积累；线粒体氧化能力的年龄相关缺陷可能限制肥大重塑所需的生物能量资源；卫星细胞功能障碍（以II型纤维特异性卫星细胞含量下降和负荷激活受损为特征）可能损害再生能力；累积蛋白质合成缺陷伴核糖体生物发生减弱也可能损害老年人对RT的反应。

研究人员近期提出假说：衰老与RT诱导的“蛋白质组可塑性”（proteome plasticity）降低相关，该术语用以表示肌肉蛋白质组对合成代谢刺激动态响应的能力，操作化为训练诱导蛋白质丰度变化的数量、幅度和功能一致性。多项证据支持这一概念：年轻与老年肌肉间的基线蛋白质组差异显著，涵盖线粒体代谢、分子伴侣和细胞骨架元件；老年人慢性转录组对RT的反应同样减弱，II型纤维中响应转录本少七倍。然而，直接比较年轻与老年人训练诱导蛋白质组重塑的研究仍然稀缺。在有限的跨年龄组RT蛋白质组学研究中，存在未对两个年龄组同时进行训练以直接比较、侧重肌肉功能而非大小变化、或采用第一代蛋白质组学技术导致动态范围检测受限等局限。本研究因此受到稀疏证据的驱动，旨在直接比较年轻与老年未经训练男性在RT前后的VL可溶性蛋白质组。

二、研究设计与方法

本研究为三项在奥本大学（Auburn University）进行的母研究（McIntosh等，2026；Ruple等，2021；Sexton等，2021）的VL肌肉活检标本的二次分析。共纳入17名年轻（21.9±2.5岁）和15名老年（57.5±6.9岁）未经训练男性，入组前≥6个月未进行结构化RT，BMI在18.5–35 kg/m2之间，无肌肉活检或抗阻运动禁忌症，未使用影响肌肉代谢的药物。所有参与者完成每周两次、持续10–12周的有监督渐进式RT，训练内容包括下肢（45°腿举、腿屈伸、俯卧腿弯举、深蹲或六角杠硬拉变式）和上肢（卧推变式、绳索下拉或划船）动作，每动作3–5组、每组6–15次，以主观疲劳等级或保留次数指导渐进超负荷。训练量负荷（组数×次数×负荷）全程记录并在年龄组间比较。

VL肌肉厚度采用B型超声（NextGen LOGIQe R8；GE Healthcare）测量，以训练前活检瘢痕为引导复现训练前后测量位置。身体成分采用双能X射线吸收法（DXA）（Lunar Prodigy；GE Corporation）评估，获得去骨瘦体重（LBM）和全身脂肪量。VL活检在基线和末次训练后≥72小时采用5号穿刺针负压技术从右侧VL肌腹中部获取，约50–100 mg组织，清除可见血液、脂肪和结缔组织后液氮速冻，?80°C保存。蛋白质组学分析中，采用含20 mM Tris–HCl（pH 7.5）、150 mM NaCl、1 mM Na2EDTA、1 mM EGTA、1% Triton、2.5 mM焦磷酸钠、1 mM β-甘油磷酸、1 mM Na3VO4、1 μg/mL亮肽酶的1X通用组织裂解缓冲液（Cell Signaling Technology；货号：9803）提取可溶性蛋白组分，500×g离心5分钟保留上清（即“裂解物”）。该裂解程序倾向于溶解膜、胞质、核及部分收缩和细胞外基质蛋白，但500×g离心保留了大部分肌原纤维蛋白（按干重计约占肌肉蛋白的80%）于沉淀中，因此分析组分相对于全肌肉 depleted 肌原纤维蛋白、过表达肌浆蛋白，故描述为可溶性蛋白质组而非肌浆富集组分。

蛋白裂解物（每样本50 μg）采用EasyPep Mini MS样品制备试剂盒（Thermo Fisher；货号：A4006）制备用于LC–MS/MS分析，包括还原、烷基化、Trypsin/Lys-C蛋白酶消化和肽段纯化。质谱分析采用外标校准的Thermo Orbitrap Exploris 240高分辨电喷雾串联质谱仪，配合Dionex UltiMate3000 RSLC Nano系统，135分钟线性梯度（3%–50% B），数据非依赖性采集（DIA）模式，10 m/z隔离窗口、1 m/z重叠、碰撞能量28、分辨率15,000，每样本三个技术重复。原始DIA文件采用DIA-NN（v 2.2.0）处理，所有原始文件单批次联合分析以实现保留时间对齐和跨样本匹配（MBR）定量，以DPHLv2谱图库指导肽段鉴定，采用经审查的智人UniProt参考蛋白质组（2020-04版本）补充iRT校准肽、反向诱饵和常见污染物进行蛋白质重注释，肽段和蛋白质水平FDR控制在1%。

质量控制流程包括：按Human Protein Atlas（25.0版本，2026年2月6日下载）注释为骨骼肌表达的蛋白质按UniProt登录号去重，保留在至少一个实验组中≥10个样本检出的蛋白质；去除高丰度血液蛋白和常见污染物；采用四方法共识（每样本缺失率、PCA马氏距离、中位强度MAD、样本间Pearson相关）识别样本离群值，移除两个单时间点样本并保留其未配对伙伴；丰度经proteoDA R包循环loess归一化。过滤后保留2106种蛋白质用于下游分析。缺失值采用missForest（ntree=100，maxiter=10）插补，该法从16种方法的基准测试中选出，在倍数变化保留、通路水平一致性和重构误差方面综合表现最佳。

差异表达蛋白（DEP）丰度采用limma R包经验贝叶斯调节评估，模型通过duplicateCorrelation函数估计受试者内相关以考虑配对PRE–POST设计。定义四个对比：年轻人训练效应（Young_Post?Young_Pre）、老年人训练效应（Old_Post?Old_Pre）、基线衰老效应（Old_Pre?Young_Pre）和年龄×训练交互（[Old_Post?Old_Pre]?[Young_Post?Young_Pre]）。显著性采用Π-score（整合统计显著性和倍数变化幅度的复合指标，Π=p^|log2FC|，Π<0.05为显著），同时以Benjamini–Hochberg校正FDR（<0.05）和名义p值（<0.05）作为补充阈值。蛋白质集富集分析（PSEA）采用两种主要方法：以limma调节t统计量排序的蛋白质列表对MSigDB Hallmark、KEGG Medicus、Reactome、GO:BP和GO Slim基因集（大小15–500；排除疾病、癌症和感染相关条目）进行fgseaMultilevel秩列表富集，每数据库分别进行Benjamini–Hochberg校正；以及对比先验Hallmark和GO Slim集合的对比水平NES值进行一致性和逆转分析。过表达分析（ORA）采用fgsea::fora对散点图象限划分的蛋白质进行，经BH校正和Jaccard去冗余（cutoff=0.5）。一致性分析计算训练效应logFC值间的Spearman秩相关，采用Rank-Rank Hypergeometric Overlap（RRHO2）进行无阈值一致性评估，全局蛋白质组组成通过主成分分析（PCA）可视化，年龄、训练时间点及其交互的贡献通过置换多元方差分析（PERMANOVA；adonis2，999次置换）量化。逆转分析将每蛋白质的衰老效应（Old_Pre?Young_Pre）与其在老年人中的训练响应（Old_Post?Old_Pre）比较，分类为逆转、加剧或可忽略（|log2FC|<0.2）；由于两对比共享Old_Pre且符号相反，其相关在引入任何生物学之前即为负，因此计算三种耦合感知统计量：受试者内Pre/Post标签交换零分布（1000次抽取）、老年人不相交半数的分裂半估计（500次抽取）、以及与不共享样本的年轻人训练响应的相关。加权基因共表达网络分析（WGCNA）构建2106种蛋白质在62个样本上的带符号网络，模块-性状关联采用线性混合模型（eigengene~group+(1|subject)）及Kenward-Roger自由度评估以考虑配对设计。所有统计分析和可视化在R（4.6.0版本）中完成，代码和已处理数据可在https://github.com/Dustyn-T-Lewis/YvO_2026获取。

三、研究结果

3.1 参与者特征与训练结果

训练量负荷在年龄组间相似（p=0.771）。两年龄组训练后DXA瘦体重均增加（年轻：+2.73 kg，老年：+1.18 kg；组间p=0.009），VL厚度均增加（年轻：+0.49 cm，老年：+0.20 cm）。年轻人VL厚度增幅显著大于老年人（Δ=0.284 cm，p=0.002），确认在可比训练刺激下VL存在差异性肥大反应。年轻人硬拉1RM增幅也更大（交互p<0.001；n=17年轻，9老年），II型纤维横截面积增幅更大（交互p=0.003；n=13年轻，10老年）。老年人基线脂肪量（30.2±6.7 vs. 19.4±9.0 kg，p=0.001）和脂肪-瘦体重比（0.52±0.11 vs. 0.34±0.15，p<0.001）均高于年轻人，两项指标训练后两组均无明显变化。

3.2 蛋白质组学概览

质量控制过滤后，2106种蛋白质在来自17名年轻和15名老年人的62份VL活检样本中被定量。主成分分析显示年龄组部分分离，但训练状态未分离。PERMANOVA检测到年龄（R2=0.089，p=0.001）和训练时间点（R2=0.029，p=0.001）的显著主效应，而年龄×时间交互不显著（R2=0.013，p=0.260）。各对比效应量分布差异显著：衰老对比中位|log2FC|=0.18，年轻人训练=0.13，老年人训练=0.09，显示从衰老特征到老年人训练响应的逐级压缩。核心结果是年龄组间在FDR水平上的差距：训练在年轻人中改变了135种蛋白质（FDR<0.05），在老年人中无一改变（最小校正p=0.156，由HPRT1和FSCN1共享）。相同不对称性在所有三个阈值下均存在：衰老对比在各阈值下产生最多DEP（Π<0.05时195种；FDR<0.05时278种；p<0.05时570种），其次为年轻人训练（99；135；448），老年人训练产生Π<0.05时18种和p<0.05时196种。交互对比识别出33种DEP（Π<0.05），确认老年人反应减弱在统计上稳健。衰老显著蛋白中最强 depleted 的是肌节组分（LMOD2，log2FC=?2.35；ACTN3，?2.12），最 elevated 的包括GPT2（+1.68，氨基酸分解代谢）和PLIN2（+1.12，脂滴储存）。UpSet交集分析显示294种独特Π-score DEP主要为对比特异性，衰老贡献最大独占集（174种蛋白质，59%），其次为年轻人训练（62种，21%）。蛋白质集富集分析（fGSEA）在四个对比中识别出916个显著富集的通路-对比组合（padj<0.05）。年轻人训练和交互对比产生最多通路水平信号（分别为324和282条通路），老年人训练产生4.3倍更少的富集通路（75 vs. 324），但75条仍具连贯性：52条正向富集，最强为翻译起始（NES=+2.85，padj=1.3×10?1?）、真核翻译延伸（+2.83）、eIF2AK4/GCN2氨基酸缺乏反应（+2.83）、无义介导衰变（+2.79）和rRNA加工（+2.77）。交互对比尽管仅33种个体DEP，却识别出282条通路的实质性通路水平信号，表明年龄依赖性训练差异表现为功能相关蛋白质的协调变化。

3.3 老年人训练反应方向保守但幅度压缩

将年轻人135种FDR<0.05的训练响应蛋白限制分析，110种（81.5%）在老年人中同向变化（Spearman ρ=0.711）。跨全部2094种在两年龄组均定量的蛋白质，同向比例为61.7%，ρ=0.288（p=3.4×10??1）。符号一致性随效应量单调上升，从最低|log2FC|五分位的52.5%升至最高五分位的74.0%（趋势p=1.7×10?12）。25种反向变化蛋白未形成独特生物学类别：11种为复合体I亚基或线粒体输入组分，无一在老年人中达到FDR<0.10，其在老年人中的中位|log2FC|为0.038（对比110种同向蛋白的0.151，p=3.6×10??），故反向符号被赋予一个与零不可区分的效应。110种同向蛋白富集于伴侣蛋白介导的微管蛋白折叠、细胞骨架依赖性胞内转运和胞质tRNA氨酰化（FDR<0.05时48个条目）。老年人效应量系统性压缩（中位|log2FC|：0.13 vs. 0.09），与完整但减弱的训练反应一致。通路水平一致性（Spearman ρ=0.571，p=1.65×10??）强于蛋白质水平一致性。五条通路在两年龄组均达显著，全部方向一致。最大年龄分化信号在氧化磷酸化（Hallmark；NES：年轻?2.08 vs. 老年?0.61）。63条精选Hallmark和GO Slim通路中24条在至少一个年龄组达显著。120种在任一训练或交互对比中达Π-score显著的蛋白质中，42种一致上调、38种一致下调、40种年龄组间反应不一致。一致下调蛋白富集于细胞外基质组织。旋转基因集检验（fry）确认年轻人在训练中上调的蛋白质（n=43 Π-DEPs）在老年人中一致响应（p=0.003）。

3.4 老年人训练对抗衰老特征

比较衰老对比（195种Π-score DEP）与老年人训练对比，209种Π-score DEP中157种（75.1%）在衰老与老年人训练对比间呈反号变化。通路水平环形插图显示显著镜像：随衰老下调的通路随训练上调，反之亦然。在衰老或老年人训练对比中达Π-score显著的209种蛋白质中，157种（75.1%）落在逆转象限。在278种衰老FDR<0.05蛋白质上，衰老和老年人训练系数相关r=?0.562，77.0%在反号象限；在组均值轴加top-n选择下观测值r=?0.488。针对通过交换每名老年人受试者内Pre和Post标签并重新选择蛋白质来保持耦合的零分布，该值落在零分布内（零均值=?0.267，95% CI ?0.715至0.241，p=0.25）。在老年人不相交半数中测量两轴得r=?0.171（95% CI ?0.548至0.344）。通路水平上，18条通路在任一或两对比中显著富集，胞质翻译显示最强逆转（衰老NES=?2.93，训练+2.77）。30种最高幅度衰老DEP中多数显示方向性逆转，富集于肌肉细胞发育和肌生成。由于两对比共享老年基线，此报告为描述性模式而非训练逆转衰老蛋白质组的证据。

3.5 加权基因共表达网络揭示年龄和训练相关蛋白模块

从2106种蛋白质在62个样本上构建带符号WGCNA网络，识别出9个共表达模块（70–340种蛋白质；1558种分配，548种未分配）。网络构建采用Pearson相关，软阈值幂12，依据两个一致标准：达到无标度拟合R2=0.877（对照0.87截断值，按Langfelder&Horvath对小型样本研究的指导从常规0.90放宽）和WGCNA FAQ对超过40个样本（此处62）带符号网络的推荐。模块-性状关联经线性混合模型（eigengene~group+(1|subject)）及Kenward-Roger自由度评估。年龄是主导组织原则，三个模块达BH校正显著：Turquoise（340种蛋白质；脂质分解代谢|TCA循环；r=+0.58，pBH=2.5×10??）、Black（107种；糖酵解|翻译因子；r=?0.60，pBH=2.3×10??）、Yellow（210种；肌肉收缩|肌节；r=?0.44，pBH=0.008）。Blue模块（272种；伴侣蛋白/蛋白质稳态）是唯一与年轻人训练显著相关的模块（r=+0.62，pBH=0.002）。无模块在老年人训练或交互对比中达显著。年龄升高Turquoise模块由有机酸分解代谢过程（padj=6.7×10?22）和脂肪酸分解代谢主导。年龄 depleted Yellow模块富集于横纹肌收缩（padj=2.0×10?2?）和肌生成。年龄 depleted Black模块注释最弱：无GO生物学过程条目达显著，最强为葡萄糖分解代谢（padj=0.083），但KEGG糖酵解（padj=1.2×10??）和Hallmark缺氧（padj=3.8×10?3）显著。训练响应Blue模块富集于蛋白质折叠过程（padj=3.4×10??）和Hsp90伴侣循环（padj=9.0×10??）。Brown富集于蛋白酶体降解（padj=1.1×10??）和糖异生（padj=1.1×10??）。Green由ATP合成偶联电子传递主导（padj=9.9×10???）。Magenta富集于翻译起始，Pink富集于ATP合酶和化学渗透偶联。八/九个模块在训练前后网络间强保守（Z>10），Magenta中度保守（Z=8.1）。

3.6 基线蛋白质组特征区分年龄组

使用模块特征基因对六个结局进行单变量模块-结局区分。最强信号为年龄区分：Black/M6（糖酵解|翻译因子）完全分离年轻与老年（AUC=1.00，置换p<0.001，q<0.010），其次Turquoise/M1（脂质分解代谢|TCA循环；AUC=0.96，q<0.010）、Blue/M2（蛋白质折叠|伴侣蛋白；AUC=0.90，q<0.010）和Yellow/M4（肌肉收缩|肌节；AUC=0.80，q=0.032）。受试者内训练前后区分较弱（合并AUC 0.47–0.73），无单一模块对ΔVL或ΔLBM响应者状态达BH显著。留一受试者外验证确认样本内AUC未受乐观偏差膨胀：固定模块投影和全WGCNA重拟合LOSO均复现Turquoise（0.96/0.96）、Blue（0.92/0.93）和Yellow（0.80/0.81）的折叠外AUC，Black两条路径均返回1.00但重拟合值基于六个折叠。受试者内变化（ΔME）不区分年龄（AUC=0.48，p=0.91），而Pre+Post+ΔME组合堆叠达AUC=0.96（p=0.017）。原始蛋白（k=10）年龄分类器达AUC=0.991（置换p<0.005），而仅表型区分不显著（AUC=0.78，p=0.17）。Brown和Black是留一受试者外网络重拟合下最不可复现模块，其余七模块（包括年龄区分Turquoise、Blue、Yellow、Red和Green）可靠复现。模块特征基因也追踪表型适应并具明显年龄分层：年轻人中氧化磷酸化特征基因变化追踪VL厚度变化（r=+0.67，原始p=0.006），基线肌肉收缩追踪基线瘦体重（r=?0.50，原始p=0.047，n=16）；老年人关联较弱或为零。模块-表型关联来自180检验筛选（9模块×2预测因子×5结局变化×2年龄层），基线横断面相关构成每性状家族10模块的独立家族，两家族均无BH校正存活，应视为探索性年龄依赖模式。年龄相关Turquoise模块（340种蛋白质）中，10个枢纽候选经模块成员度（kME>0.75）识别，顶部枢纽蛋白（ETFB、IDH2、ALDH5A1、ETFA、CRAT）以脂肪酸和有机酸分解代谢酶为主，强化该年龄相关模块的代谢身份。

四、讨论与结论

研究产生五项主要发现。第一，尽管两年龄组实施相似10–12周方案，老年人表现出显著减弱的蛋白质组重塑，训练响应蛋白少5.5倍且零FDR显著蛋白。第二，尽管压缩，年轻人识别出的训练响应蛋白81.5%在老年人中同向变化，表明老年人训练反应保守但幅度受限。第三，老年人训练诱导变化与衰老效应方向相反，在核糖体和翻译机器中最明显。第四，共表达网络分析识别出年龄分化生物学程序，氧化磷酸化成为对训练介导逆转最抵抗且在年龄组间最分化的通路。第五，基线蛋白质组特征区分年龄组，而无模块-表型关联在校正后达显著。综合而言，衰老限制蛋白质组重塑的幅度但不使潜在适应程序不可变，基线蛋白质组状态是否设定反应能力边界仍悬而未决，而训练可在这些边界内部分恢复个体蛋白质丰度。

讨论中，研究人员指出老年人训练反应缺陷幅度可观：Π-score阈值下年轻人产生99种DEP对比老年人18种（5.5倍减弱），FDR<0.05下绝对悬殊（135 vs. 0），且放宽至FDR<0.10不改变此结论（老年人训练对比中最小校正p=0.156，HPRT1）。老年人训练叙述因此基于通路水平证据，其中75个基因集达显著且核糖体和翻译机器协调上调（NES最高+2.85，padj=1.3×10?1?），同时无共表达模块在老年人训练中达显著。Π-score等复合指标在低功效对比中排序效应仍有用，但不替代多重性控制。所有显著性阈值的一致模式支持真实生物学减弱，与O'Leary等（2025）报道两年龄组均零FDR显著蛋白但检测到显著通路一致，也与研究人员既往蛋白质组学工作显示衰老主要影响非肌原纤维蛋白池且训练诱导变化边缘一致。尽管定量缺陷，训练反应未发生质变：年轻人135种FDR<0.05训练响应蛋白中蛋白质水平一致性为81.5%（ρ=0.711），全蛋白质组为61.7%（ρ=0.288，p=3.4×10??1），差距由效应过小致符号无信息性的蛋白质解释而非分化生物学。秩-秩图确认强一致性且可忽略不一致信号。这种“保守伴压缩”模式提示旨在放大现有反应而非重定向其轨迹的治疗策略可能对老年人长期肌肉发展有前景。通路水平一致性（ρ=0.571，p=1.65×10??）大幅超过蛋白质水平一致性，表明协调功能程序跨年龄保存优于个体蛋白质效应量。最大年龄分化信号在氧化磷酸化（Hallmark；NES：年轻?2.08 vs. 老年?0.61）。衰老和训练共享最小蛋白质水平重叠（85%的Π-score DEP为对比专属），但年龄组间通路一致性（ρ=0.571）近乎两倍于蛋白质一致性（ρ=0.288）。尽管年龄解释适度全局方差（PERMANOVA R2=8.9%），识别出195种DEP，表明衰老通过许多小而协调的变化重塑蛋白质组，最佳在通路和网络水平捕获。

关于训练部分逆转衰老蛋白质组，研究人员指出老年人训练诱导变化倾向于与衰老效应方向相反，但两对比不独立。209种在衰老和老年人训练对比中达Π<0.05的蛋白质中157种（75.1%）呈反号变化，前30种衰老影响DEP中22种（73%）显示方向性逆转。最强逆转在核糖体和翻译机器，胞质翻译显示完全方向反转（衰老NES=?2.93；训练NES=+2.77）。相比之下，氧化磷酸化显示最小逆转，表明翻译程序比线粒体代谢程序对运动介导恢复更响应。这些方向计数是描述性的。衰老对比和老年人训练对比均含Old_Pre且符号相反，故两者间负相关仅凭结构即可预期。一旦在零分布中保持该耦合，相关不再可区分（观测r=?0.488对比零均值?0.267，p=0.25），在不相交老年人半数中测量两轴的分裂半估计得r=?0.171且置信区间跨零。因此逆转模式被视为产生假说而非确立，将其与共享基线分离需在测量训练响应的受试者之外定义衰老特征。

关于共表达网络架构，WGCNA识别出组织1558/2106种蛋白质的九个模块（548种未分配），模块水平图景是每蛋白质和每通路所见压缩的再上一级聚合。年龄组织网络：Turquoise（脂质分解代谢|TCA循环；r=+0.58，pBH=2.5×10??）、Black（糖酵解|翻译因子；r=?0.60，pBH=2.3×10??）和Yellow（肌肉收缩|肌节；r=?0.44，pBH=0.008）在考虑配对设计的线性混合模型下均分离年龄组。合读描述一个转变而非三个独立发现：燃料氧化上调、翻译机器下调、收缩装置下调。训练仅达一个模块，Blue（蛋白质折叠|伴侣蛋白；r=+0.62，pBH=0.002），且仅在年轻人中；无模块在老年人训练或交互中达显著。模块水平富集显示相同两模式更粗粒度：训练富集在年龄组间分化，仅Blue（伴侣蛋白和蛋白质折叠）在年轻和老年人中均达显著（NES+3.95和+3.07），四模块仅年轻人显著、四模块仅老年人显著。对衰老对比模式反转：在衰老和老年人训练中均显著的五模块全部携带反号，Turquoise和Brown随年龄上升随训练下降，Blue、Magenta和Red反之。与蛋白质和通路水平一致，这些是描述性的，因两对比共享Old_Pre项且不独立。同一排序因此在检查的每一聚合水平均成立——大衰老信号、较小年轻人训练信号、以及数据无法解析的老年人训练信号——而架构本身保持不动：八/九个模块在训练前后网络间强保守（Z>10），Magenta中度保守（Z=8.1）。骨骼肌蛋白质组的共表达架构因此在干预期间基本稳定，即使个体蛋白质丰度变化。在此稳定框架内，Green（氧化磷酸化|复合体I）携带一个结局邻近信号：其特征基因变化追踪年轻人VL肥大（r=+0.67，原始p=0.006），未在该层内校正存活（pBH=0.058），故报告为探索性且不解读方向。可复现年龄信号由Turquoise和Yellow携带。Black返回最大年龄相关（r=?0.60，pBH=2.3×10??）但为网络中最不可复现模块，故解读为与翻译转变一致而非其证据，Red和Magenta同向但未达显著（pBH=0.11和0.13）；通路水平上六个最强年龄 depleted 基因集均为核糖体或翻译相关（NES?3.10至?2.94）。一致性是方向性而非组成性——Black贡献胞质翻译93个前缘蛋白中8个，对比Red 37个和Magenta 27个——但两水平均将翻译机器置于衰老信号中心。Brown和Black是留一受试者外重拟合下最不可复现模块，故其组成及从中读取的任何生物学应视为临时。

关于线粒体和代谢通路重塑，年轻人训练后显示强氧化磷酸化抑制（Hallmark OxPhos NES=?2.08），可能反映生物合成资源转向生长相关程序，而此反应在老年人中几乎完全缺失（NES=?0.61），产生所有通路中最强交互效应（NES=+2.02）。衰老本身与氧化磷酸化基因产物升高相关（NES=+2.04），可能代表对线粒体效率下降的补偿性上调。Robinson等（2017）报道RT增加老年但非年轻成人线粒体蛋白质合成，提示老年肌肉通过增强线粒体生物发生而非年轻人观察到的糖酵解转变来响应能量需求，尽管Jessen等（2026）最近显示RT诱导的OxPhos抑制主要发生在I型纤维，提示全肌肉信号反映纤维类型特异性重塑。这种年龄特异性未能重塑线粒体蛋白含量可能反映线粒体质量控制受损，线粒体功能障碍日益被认为是生物衰老的标志和潜在驱动因素。MoTrPAC线粒体图谱进一步证明耐力训练期间组织特异性线粒体多组学重塑，骨骼肌在RT背景下显示缺乏协调氧化蛋白上调。

关于基线蛋白质组特征区分年龄，模块特征基因将年龄和训练相关信号分解为不同生物学程序。单变量每模块ROC表明年龄分离主要由代谢连接模块（Turquoise（脂质分解代谢|TCA循环），AUC=0.96，q<0.010）和蛋白质稳态（Blue（蛋白质折叠|伴侣蛋白），AUC=0.90，q<0.010）携带，Yellow（肌肉收缩|肌节）提供较弱但显著信号（AUC=0.80，q=0.032），而受试者内训练时间区分在所有模块中较弱（合并Pre vs. Post AUC 0.47–0.73）。模块特征基因也追踪表型适应并具明显年龄分层：年轻人中氧化磷酸化特征基因变化追踪VL厚度变化（r=+0.67，原始p=0.006），基线肌肉收缩追踪基线瘦体重（r=?0.50，原始p=0.047，n=16），老年人关联较弱或为零。

关于机制和临床意义，多种独立分析方法——差异表达、一致性和网络分析——在共同生物学主题上的收敛提供了衰老如何影响蛋白质组的洞见。蛋白质水平显著减弱伴通路水平一致性保存，提示衰老在转录和转录后水平引入瓶颈，压缩个体蛋白质效应量同时保留功能程序协调。研究人员因此提出，最大蛋白质组重塑所需的合成代谢信号基础设施随衰老减少，但本研究设计无法将该约束定位至特定信号步骤。若保守伴压缩模式成立，将支持放大现有反应幅度而非重定向反应轨迹的治疗策略。核糖体生物发生通路成为特别有吸引力的干预靶点，因核糖体蛋白是老年人训练中最强恢复的蛋白质之一，且核糖体生物发生能力预测老年人肥大结局（Stec等，2016）；尽管翻译能力与效率的相对重要性可能随训练进展而转变，延长RT后更大肥大与更低核糖体相关基因表达相关（Brown等，2024）。虽然推测性，旨在增强核糖体生物发生和效率的药理学、营养学或其他干预策略——如亮氨酸补充或抗阻运动结合必需氨基酸供给——可能增强老年人蛋白质组重塑幅度。

关于局限性，第一，分析为二次分析，源自共享训练方案但不同营养补充范式（花生蛋白和甜菜根汁）的母RT研究，受试者水平随机截距解释受试者内相关。补充在此合并样本中完全嵌套于年龄组内，故补充项与年龄项混淆且未拟合。基线层差异因此无法与母研究间队列差异分离，而此处报告的训练响应似乎不依赖补充。第二，样本量（n=17年轻，n=15老年）为检测年轻人中等到大训练效应提供充分功效，但交互对比功效有限。老年人训练中完全缺失FDR显著发现可能既反映生物学减弱也反映功效限制。第三，老年人基线携带显著更多脂肪量和更高脂肪-瘦体重比，两者均未随训练差异变化。肥胖因此无法解释训练反应差异（受试者内测量），但在横断面衰老对比中无法与实足年龄分离且未调整。第四，衰老组分是横断面的，年龄组间比较可能反映队列差异而非生物学衰老本身，单一训练后时间点无法捕获蛋白质组重塑动力学。第五，所有参与者为男性，结果可能因骨骼肌生理和训练适应的已知性别差异而无法推广至女性。此外，所有活检取自VL，结果可能无法推广至不同纤维类型组成或功能需求的肌肉，特别是上肢或主要慢肌纤维肌肉如比目鱼肌。运动诱导肌肉损伤恢复的年龄相关损害也可能根据所用RT方案产生不同适应反应