《GeroScience》：Quadriceps mitochondrial DNA quantity, quality, and gene expression after 2 years of calorie restriction: exploratory results from the CALERIE trial
编辑推荐：
为探讨热量限制(CR)对非肥胖人群肌肉线粒体DNA(mtDNA)完整性的影响，研究人员对CALERIE?试验的肌肉活检样本进行了为期两年的纵向分析，运用液滴数字PCR(ddPCR)技术检测mtDNA拷贝数与缺失突变频率，并关联了基因表达数据。研究发现了mtDNA拷贝数与性别、时间、基础代谢指标的相关性，并揭示了与线粒体功能相关的基因表达富集，为理解CR对线粒体基因组适应性和其作为代谢健康生物标志物的潜力提供了新见解。
衰老及其相关疾病是当今社会面临的主要健康挑战，严重影响生活质量并带来沉重的医疗负担。在此过程中，线粒体功能障碍被公认为是衰老的标志之一，并与多种年龄相关的代谢、退行性疾病乃至癌症的进展密切相关。线粒体不仅是细胞的“能量工厂”，还在细胞信号传导、分化和程序性死亡中扮演关键角色。而线粒体DNA(mtDNA)的完整性，作为线粒体质量的核心要素，其变化与衰老进程紧密相连。在人类骨骼肌中，mtDNA拷贝数随年龄增长而线性下降，而mtDNA缺失突变频率则呈指数级上升，这些损伤会破坏呼吸链功能，最终可能导致细胞死亡和肌肉衰减。那么，有没有一种干预措施能够“延缓”甚至“修复”这种与年龄相关的线粒体衰退呢？
热量限制(Caloric Restriction, CR)，即在保证营养均衡的前提下减少热量摄入，在多种模式生物中已被证明能够延长寿命并延缓年龄相关疾病的发生。其潜在的益处包括改善线粒体功能、促进线粒体生物合成、减少氧化应激等。然而，关于CR对人类骨骼肌mtDNA的影响，现有的研究结果却存在矛盾。有的短期(6个月)研究发现CR能显著增加肌肉mtDNA拷贝数，而另一项为期12个月的研究则发现拷贝数随时间下降，且CR组与随意进食(ad libitum, AL)对照组之间无显著差异。此外，这些研究大多未考察CR对mtDNA突变的影响。这些不确定性凸显了在更长期、更严谨的人体试验中，系统评估CR对线粒体基因组影响的必要性。
为此，一项发表于《GeroScience》杂志的最新研究，利用具有里程碑意义的人体长期热量限制试验——CALERIE?(Comprehensive Assessment of Long-term Effects of Reducing Intake of Energy)项目的生物样本，进行了深入的探索性分析。CALERIE?是一项为期两年的多中心随机对照试验，旨在比较对非肥胖成年人进行25%的热量限制与随意进食对照之间的差异。本研究的目标明确：评估利用最新技术分析线粒体DNA完整性的可行性，探索mtDNA指标与生理、转录组数据的关联，为未来假设提供信息，估算效应量为后续研究提供统计效力参考，并最终为庞大的CALERIE数据库增添宝贵的线粒体基因组数据。
主要技术方法
本研究对来自CALERIE试验参与者的股四头肌活检样本(基线n=93，12个月n=44，24个月n=31)进行了分析。关键技术包括：1) 使用液滴数字PCR(droplet digital PCR, ddPCR)对提取的总DNA样本进行绝对定量，精确检测每个二倍体细胞核对应的核DNA(nDNA)拷贝数、mtDNA拷贝数以及mtDNA缺失突变频率，该方法无需标准曲线，结果不依赖于PCR效率。2) 对来自相同肌肉活检的现有基因表达数据进行了关联分析。3) 在统计学上，采用线性混合效应模型(Linear Mixed-Effects Model, LME)评估mtDNA指标在组内和组间随时间的纵向变化，并在基线水平分析了mtDNA指标与生理指标、基因表达的相关性。4) 使用预排序基因集富集分析(PreRank Gene Set Enrichment Analysis, GSEAPreranked)对与mtDNA指标相关的基因列表进行通路和功能富集分析。
研究结果
1. 股四头肌mtDNA拷贝数的纵向分析
线性混合效应模型分析显示，时间(β = -27.309, p = 0.0001)和性别(β = 688.059, p = 0.0005)对mtDNA拷贝数有显著影响。具体而言，在研究的两年期间，mtDNA拷贝数总体呈下降趋势，且女性的拷贝数显著低于男性。然而，CR组与AL组之间(p = 0.2898)，以及年龄对拷贝数的影响(p = 0.4644)均未达到统计学显著性。
2. 股四头肌mtDNA缺失突变频率的纵向分析
对log10转换后的mtDNA缺失突变频率的分析显示，其随年龄增长而显著增加(p = 0.0041)。但时间、CR干预或性别对其均未产生统计学上的显著影响。
3. mtDNA拷贝数与突变频率预测肌肉相关生理指标
在基线数据中，mtDNA拷贝数与多项生理指标呈弱到中度正相关，包括去脂体重(r = 0.43)、自我报告的能量摄入(r = 0.34)、静息代谢率(r = 0.42)、总能量消耗(r = 0.31)以及最大摄氧量(V?O2max, r = 0.41)。mtDNA缺失突变频率则与V?O2max呈较弱的负相关(r = -0.23)。
4. mtDNA拷贝数与基线基因表达的关联
在调整了年龄、性别、种族和技术批次后，基线时mtDNA拷贝数与606个基因的表达相关。富集分析显示，上调的基因显著富集于跨膜电子传递、氧化磷酸化等线粒体功能相关通路，以及与神经退行性疾病相关的通路。而下调的基因则富集于核基因组维持、RNA剪接和核糖体相关功能。
5. 纵向研究中mtDNA拷贝数与基因表达的关联
在纵向分析中，与mtDNA拷贝数相关的基因数量减少至286个，其中仅6个与基线分析重叠。一个在基线和纵向分析中均相关的基因是MTRNR2L12，它属于人类线粒体衍生肽家族，具有多种保护作用。纵向分析中上调的通路转向细胞外基质、膜蛋白和免疫细胞激活，而下调的通路则转向线粒体和蛋白质翻译。
6. mtDNA缺失突变频率与基因表达的关联
基线时，mtDNA缺失突变频率与362个基因的表达相关，上调通路主要与炎症相关，下调通路则与蛋白质翻译和线粒体功能相关。纵向分析中，相关基因减少至212个，其中与基线重叠的仅2个。一个在纵向分析中相关的基因是IRGM，它参与自噬过程。与拷贝数分析类似，CR干预往往削弱了AL组中观察到的mtDNA突变频率与基因表达之间的变化关联。
研究结论与讨论
本研究揭示了关于人体股四头肌在CR干预期间线粒体基因组的几项新发现。首先，性别和时间显著影响mtDNA拷贝数，女性拷贝数更低，且拷贝数在两年期间呈下降趋势。其次，平均来看，CR组的mtDNA缺失突变频率较低，且在男女间存在差异，但由于样本量限制，这些差异未达统计学显著性。最后，也是最重要的发现是，mtDNA拷贝数和突变频率与大量基因的表达相关，其中一些基因富集于氧化磷酸化和线粒体功能通路，这提示了线粒体与细胞核基因组之间存在活跃的信号对话。
研究结果也验证了一些已知现象：即使在相对年轻的参与者中，年龄增长也与更高的mtDNA缺失突变频率相关，且mtDNA指标能够预测部分与骨骼肌生理和代谢相关的指标，进一步支持了其作为代谢健康生物标志物的潜力。
关于CR对mtDNA的影响，本研究在为期两年的干预中未发现CR对肌肉mtDNA拷贝数或缺失突变频率有统计学上的显著影响。这与之前一些短期、结果矛盾的研究形成了对比。差异可能源于CR持续时间、限制程度或mtDNA测量方法的不同。研究者指出，缺乏标准化的mtDNA测量方案是导致研究间难以直接比较的重要原因。
这项研究的深远意义在于，它为解决衰老生物学中的一个核心问题——如何量化“生物年龄”并评估干预措施效果——提供了新的思路。延长生存期并非临床可行的终点，因此需要可靠的生物标志物来预测健康寿命、筛选干预方案并评估其效果。线粒体生物能量学与寿命密切相关，是药物开发和治疗试验（如CR、二甲双胍、运动）的有吸引力的靶点。本研究表明，mtDNA的完整性指标，特别是与代谢指标和线粒体功能基因表达密切相关的拷贝数，有望发展成为评估代谢健康和衰老进程的生物标志物。
当然，研究也存在局限性。为期两年的CR干预在人类寿命尺度上仍属短期，可能不足以揭示CR对mtDNA突变累积的长期影响。此外，研究样本的多样性有限，参与者主要为健康、非肥胖的成年人，这限制了结论向患有慢性病、老年人或不同遗传背景人群的推广。未来需要更长期、参与者更多样化的CR研究，并改进活检技术以提高参与率和样本量，从而更全面地揭示CR、线粒体基因组与人类健康衰老之间复杂的相互作用网络。
总而言之，这项基于CALERIE珍贵样本的探索性研究，不仅为理解热量限制如何影响人类骨骼肌线粒体基因组适应性增添了关键数据，也突出了mtDNA完整性作为潜在生物标志物在推动健康衰老干预措施研发中的价值。它为未来的相关研究奠定了方法学基础并指明了方向，持续推动着人类对延缓衰老、提升生命质量的科学探索。
