小鼠衰老的代谢足迹
,1 ,2 ,三 ,1 ,1 ,1 ,4 ,4 ,1和a、,1
Riekelt H.Houtkooper公司
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
卡门·阿格曼
2荷兰阿姆斯特丹AZ 1105号学术医学中心实验室医学生物化学
Sander M.Houten先生
三荷兰阿姆斯特丹亚利桑那州1105号AMC实验室遗传代谢疾病
卡莱斯·坎特
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
艾伦·H·杰宁加
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
佩内洛普·A·安德烈
1瑞士洛桑大学综合技术与系统生理学实验室(EPFL),1015 Lausanne
查尔斯·托马斯
4瑞士洛桑1015号EPFL表型基因组学中心
拉斐尔·多恩伦
4瑞士洛桑1015号EPFL表型基因组学中心
克里斯蒂娜·肖恩詹斯
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
约翰·奥韦克斯
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
1瑞士洛桑大学洛桑分校综合与系统生理学实验室(EPFL),1015洛桑
2荷兰阿姆斯特丹AZ 1105号学术医学中心实验室医学生物化学
三荷兰阿姆斯特丹亚利桑那州1105号AMC实验室遗传代谢疾病
4瑞士洛桑1015号EPFL表型基因组学中心
2011年9月12日收到;2011年10月14日接受。
版权©2011,麦克米伦出版有限公司。保留所有权利 - 补充资料
补充信息补充信息
GUID:4FFBDA13-50C2-4500-816B-F1D83CB8CC33
补充信息表S1
指南:44DF2B61-74B2-43A6-A105-D99026C4191E
补充信息表S2
GUID:D77EEA6B-807B-42ED-AE65-6890EDBB607B
摘要
衰老的特点是细胞功能普遍下降,最终会影响全身的体内平衡。虽然DNA损伤和氧化应激都会导致衰老,但代谢功能障碍是衰老的共同标志,至少在无脊椎动物中是如此。由于缺乏对健康衰老哺乳动物代谢变化的全面综述,我们在这里比较了幼年和2岁小鼠的代谢参数。我们将体内表型与基因表达、生化分析和代谢组学进行了系统整合,从而确定了衰老的显著代谢足迹。在肝脏和肌肉中受影响的途径中,我们发现葡萄糖和脂肪酸代谢以及氧化还原稳态。这些变化转化为长链酰基肉碱减少,游离脂肪酸水平增加,老龄小鼠血浆中各种氨基酸显著减少。因此,这些代谢物可作为衰老和健康寿命的生物标志物。
衰老(最终死亡)是生命中不可避免的一部分,伴随着各种生理和心理疾病,包括常见的代谢、炎症、心血管和神经退行性疾病,这些疾病会缩短健康寿命。这类疾病的例子包括帕金森氏病和阿尔茨海默氏病,以及从心脏代谢的角度来看,肥胖、2型糖尿病和动脉粥样硬化。尽管这些临床表现众所周知,但导致衰老过程的复杂生物分子网络才刚刚开始被发现1有人建议几种常见的过程会导致或至少有助于衰老,包括DNA损伤、活性氧(ROS)积累和一般代谢功能障碍。到目前为止,这些大多被视为独立事件,但证据表明其中一些通路是相互关联的,正如最近DNA损伤和代谢控制之间的联系所强调的那样2,三.
无论其机制如何,衰老相关疾病的一个共同特征是代谢系统的参与,尤其是线粒体的参与1,4,5,6与衰老有关的最具特征的代谢途径是胰岛素/IGF1信号途径7。两个秀丽线虫和D.黑腹果蝇胰岛素受体突变体显示寿命延长,其作用是通过FOXO1转录因子、热休克因子HSF1和SKN1介导的5). 胰岛素/IGF1途径在哺乳动物寿命调节中的参与尚存在争议,但似乎有可能7第二条长寿途径以哺乳动物雷帕霉素(mTOR)为靶点,该靶点将胰岛素信号与其他营养物质(最显著的是氨基酸)的传感结合在一起,从而调节蛋白质翻译和自噬8,9与此一致,mTOR被证明参与了与老化相关的酮体生成减少10通过服用雷帕霉素抑制mTOR可延长小鼠寿命11与这些营养过剩系统相对的是营养限制途径,如sirtuin和AMP-activated protein kinase(AMPK)途径12,13被能量应激激活,例如在长时间禁食或运动后,它们抑制能量需求过程,有利于能量生产14sirtuins和AMPK都被认为是低等生物体寿命的积极调节因子(在1)但对于胰岛素/IGF1途径,需要证实其与哺乳动物寿命的相关性。在秀丽线虫,AMPK被证明直接磷酸化CREB调节的转录辅激活子-1,然后它就无法辅激活CREB转录因子以缩短寿命15这表明哺乳动物中可能存在一种衰老途径。此外,在鼠标中,AMPK已与SIRT1连接16,17和PGC-1α信号17,18巩固该系统与线粒体代谢的联系,尤其是线粒体呼吸。然而,由于抑制和刺激线粒体代谢似乎都能延长寿命,线粒体功能如何影响衰老仍存在争议(在1,5). 最有可能的是,生理环境和操作本身的性质可能决定衰老反应。这种表型似乎至少部分由线粒体未折叠蛋白反应介导,并以细胞非自主方式作用,涉及一种尚未确定的“有丝分裂因子”19.
一些研究报告了衰老的分子效应,例如在生物化学20,基因组21,或代谢组学水平22然而,这些研究的局限性在于这些不同数据集的多级集成问题。因此,我们的目的是获得一个无偏见的图谱,并了解衰老期间的代谢下降,这容易导致哺乳动物发生与衰老相关的疾病。我们对野生型C57BL/6J小鼠(代谢研究中最常用的菌株)衰老的临床、分子和代谢组学方面进行了表征。这使我们能够通过识别代谢物集和分子变化来定义衰老的一致表型,这些代谢物和分子变化可以用作健康跨度生物标记物,从而构成衰老的详细代谢足迹。
结果
代谢特征和身体性能
为了描述哺乳动物衰老过程中的代谢变化,我们使用了3个月龄(“年轻”)和22个月龄的C57BL/6J小鼠。我们通过广泛的临床表型研究验证了我们的衰老模型。正如预期的那样,老龄小鼠出现了生理衰退的迹象。老年小鼠明显比年轻小鼠重()主要是因为脂肪量增加(). 此外,老年小鼠出现离散性葡萄糖不耐受,这可能反映了胰岛素抵抗的发展(补充图S1在线)。在间接量热法中,当校正了瘦体重后,老年小鼠的能量消耗有减少的趋势()或与多元回归相比(第页=0.128),部分原因可能是体力活动减少(补充图S1在线)。老年小鼠的呼吸交换比(RER)低于年轻小鼠,表明老龄小鼠对脂肪有一定比例的底物偏好(). 当接受体育锻炼时,最大VO2显著减少()
衰老小鼠的临床和分子表型。(A) 体重。(B) 身体成分。(C) 间接量热法显示耗氧量(VO2)以及24小时记录期间的呼吸交换比(RER)。(D) 基本和最大VO2在代谢跑步机上。(E) 20小时自愿运动记录。(F) 年轻和老年小鼠在腓肠肌中的基因表达(G)柠檬酸合成酶(CS)的酶活性,以及腓肠肌氧化磷酸化复合物I和复合物IV。(H) 腓肠肌相关代谢信号通路的Western blot分析。pAMPKα代表AMPK的α亚单位的磷酸化/活化;免疫沉淀IRS1检测其对丝氨酸残基(p-Ser)的磷酸化作用;p-S6K1反映S6K1磷酸化/活化。Tubulin用作加载控制。数值表示为平均值±SEM;n=6–10(,n=20)。*第页≤0.05; **第页≤0.01; ***第页≤0.001.
在一次家庭式监测中证实,在间接热量测定期间,老年小鼠的体力活动减少,在此期间,小鼠可以自由使用轮子(). 老老鼠在最初的发现阶段和后来的夜间阶段跑得更少。在方形开放场地试验中也观察到了老年小鼠的这种低动力表型,该试验表明老年小鼠不如年轻小鼠敏捷(补充图S1在线)。猫道分析(用于确定小鼠的步态障碍)表明,在迈出下一步之前,单爪站立的时间会随着年龄的增长而显著增加(补充图S1在线)。
衰老组织中的分子和生化变化
衰老小鼠组织的分子特征证实了代谢-尤其是线粒体-下降的观点,这被认为是衰老过程的一部分基础。在腓肠肌、棕色脂肪组织(BAT)和肝脏中,我们观察到编码氧化磷酸化相关蛋白的各种转录物的表达下降(Ndufa2、Ndufb3、Cyt c、Cox IV、Atp5g1、Ucp3),脂肪酸氧化(FAO)(Cpt1b,阿卡德姆,阿卡德)线粒体生物发生(Tfam公司)、氧化应激防御(草皮2)和凋亡(Bax公司) (,补充图S1在线)。进一步证实了我们的研究主要捕获了与衰老相关的变化这一事实,即基因表达数据与之前发表的对C57BL/6N腓肠肌基因表达的调查完全一致,该调查比较了5个月龄和25个月龄小鼠21(补充图S1在线)。有趣的是,在这个数据集中,许多参与氧化磷酸化的基因的表达减少(例如。阿特帕2,Cox4i2公司)、粮农组织(例如。Cpt1b型,学院)和氧化应激防御(草皮2)通过热量限制(CR)恢复,这是迄今为止已知的唯一一种延长寿命的生理干预(补充图S1在线)。相反,各种抗氧化剂的表达(像素3,Gpx4型,Gpx6型)在老年小鼠中增加,在CR小鼠中恢复正常(补充图S1在线)。
虽然编码线粒体蛋白的核基因的表达普遍减少,但有些基因没有改变(例如。阿卡德尔和附件5g1腓肠肌,以及Mfn2公司和Cpt1a公司肝脏),BAT中几个基因的表达随着年龄的增长而增加(Tfam公司,Mfn2公司,Pgc-1α,考克斯IV) (,补充图S1在线)。我们假设这是由于老年小鼠激活了代偿机制来维持适当的线粒体功能。酶活性分析为维持线粒体容量可能的代偿机制提供了另一条证据。在肌肉和肝脏中,老年小鼠的柠檬酸合成酶(CS)活性增加,表明线粒体活性或数量增加(,补充图S1在线)。相反,复合物I和IV的相对活性在两个年龄组之间是相等的(,补充图S1在线)。肌肉和肝脏中CS活性的增加并非源于线粒体的生物生成,因为这些组织中线粒体DNA(mtDNA)的定量显示没有变化(补充图S1在线)。然而,在BAT中观察到线粒体DNA增加(补充图S1在线)在其他基因表达增加后Tfam公司和Pgc-1α(补充图S1在线)。线粒体功能的恶化通常反映在活性氧水平的增加。因此,我们通过4-羟基壬醛(4-HNE)(脂质过氧化的标志物)测量了活性氧损伤。正如预期的那样,老年小鼠肌肉和肝脏中的4-HNE水平增加(补充图S1在线)。
衰老中的信号通路
为了补充我们的基因表达和酶活性数据,我们还评估了代谢信号通路,重点是与衰老相关的既定代谢调节通路。在老年肝脏和肌肉中,胰岛素受体底物1(IRS-1)丝氨酸磷酸化增加,这是胰岛素抵抗的标志,尽管在这两个组织中增加的程度不同(,补充图S1在线)。在肝脏和肌肉中,p70 S6激酶磷酸化显著降低(,补充图S1在线),反映mTOR活性降低,可能是胰岛素信号受损的结果,如IRS-1丝氨酸磷酸化升高和Akt磷酸化降低所示。除了这些营养过剩途径外,我们还分析了营养限制传感器AMPK的表达和活性,发现其磷酸化(指示能量应激)在老年小鼠的肝脏和肌肉中没有改变(,补充图S1在线)。
血液中的代谢紊乱
为了进一步表征衰老时的代谢变化,我们对血浆进行了临床化学分析。血浆标志物发生了各种变化,如游离脂肪酸增加和甘油三酯水平降低(). 为了扩大我们对衰老小鼠血液代谢物的表征,并确定可能的新型衰老生物标记物,我们采取了一种针对性较低的方法,并在年轻和老年小鼠的血液中生成了广泛的代谢组学图谱。测量了75种以上的分析物,包括红细胞脂肪酸和血浆酰肉碱、氨基酸和酮体的特征,并对这些数据进行了聚类,以确定在衰老过程中反应相似的代谢物组(). 利用代谢物的这一图谱,根据年龄将小鼠分离出来,几个有趣的代谢物簇变得明显。其中一个簇含有长链和中链酰基肉碱。这些代谢产物是FAO的中间产物,在老年小鼠中减少(). 红细胞脂肪酸谱反映了脂肪酸稳态的长期变化,通常用于诊断多不饱和脂肪酸(PUFA)缺乏症。对于饱和脂肪酸,C16略有下降,而细长的C18和C20有所增加。此外,PUFA生物合成的各种中间体和产物如C20:4ω6增加,而底物C18:2ω6和C18:3ω3减少。PUFAs和血浆白蛋白的总体水平没有降低,这两种物质都与ROS的保护有关().
血浆生物化学和靶向代谢组学。(A) 血浆临床生物化学的双向分层聚类。(B) 目标代谢组学数据的双向分层聚类。FA描述的代谢物是红细胞脂肪酸,而AC表示血浆中的酰肉碱。(C) 选定代谢物的血浆水平。数值表示为方框图和胡须图;n=10;*第页≤0.05; **第页≤0.01; ***第页≤0.001.
除了脂肪代谢外,氨基酸代谢在衰老过程中也会发生显著变化。之前的研究表明,补充支链氨基酸(BCAAs)23或饮食中蛋氨酸含量的改变24,25影响寿命。我们对血浆氨基酸的表征表明,随着年龄的增长,脯氨酸、丙氨酸、丝氨酸、酪氨酸和蛋氨酸的水平明显降低,但BCAA的缬氨酸、亮氨酸和异亮氨酸水平没有降低(). 类似地,先前已经表明mTOR途径抑制生酮,并且生酮能力随着年龄的增长而降低10。与生酮减少一致,我们的分析发现,反映总酮体水平的β-羟基丁酸(BHBA)水平往往随着年龄的增长而下降().
新的衰老生物标志物
年轻和年老小鼠的血浆代谢物水平的明显差异使我们测试了这些是否可以可靠地区分两组,这是开发生物标记物的先决条件。
选择样本分类的相关变量是高维数据集(如代谢组学和基因表达研究)的一项常见任务,在这些数据集中,人们试图确定仍能实现良好预测性能的尽可能小的代谢物或基因集(即生物标记物)。随机森林(RF)是这样一种学习算法,它可以通过估计每个变量对类的重要性来实现这一点。RF分析使用随机的小鼠训练集和随机选择的少量代谢物生成许多分类树。被排除在树外的老鼠被用来估计分类错误,提供内部交叉验证。因此,“平均下降准确度”表示某种代谢物对两个测试组的分离有多大贡献,而总体“预测准确度”表示一组代谢物组合可以分离组的准确度。血液靶向代谢组学数据的RF分析定义了一组10种代谢物,这些代谢物构成了衰老的最佳预测因子,具有100%的预测准确性(). 特别是,红细胞亚油酸(C18:2ω6)降低和血浆C16-酰基肉碱降低是强有力的预测因素().
靶向代谢组学和微阵列的代谢产物和途径富集。(A) 为了将衰老状态与血液代谢物表达数据联系起来,我们使用了随机森林预测因子。根据代谢物对衰老状态的预测重要性,使用随机森林重要性度量对代谢物进行排序。我们发现,10种最重要的代谢物具有100%的明显预测准确性。(B) 监督层次聚类图,用于可视化RF预测因子的表达(z-score标准化)。RF图和监督层次聚类图表明,15种最重要的代谢物根据其老化状态对样品进行分层。FA描述的代谢物是红细胞脂肪酸,而AC表示血浆中的酰肉碱。色码描述了相应代谢物所属的途径。
在我们的血液代谢物分析和微阵列数据中出现这些明确的代谢变化后,我们还使用全球代谢组学绘制了年轻和老年小鼠肝脏和股四头肌的代谢变化图。我们在肝脏和肌肉中分别检测到473和263种代谢物,其中94种和65种在各组之间发生了显著变化(参见补充表S1和S2在线查看相对代谢物水平)。我们再次进行RF分析,对这些代谢变化进行分类(——). 在肝脏中,代谢物的变化可以93%准确预测年龄组()前25种代谢物的层次聚类证实了这两个年龄组的分离(). RF分析旨在从高维代谢组学数据中确定单个代谢物的相对贡献,随后我们进行了代谢物集富集分析(MSEA)26确定哪些途径受到影响。在肝脏中,预测年龄的最佳代谢物是葡萄糖代谢、磷脂代谢和氧化还原稳态(). 此外,我们进行了代谢组视图分析,以证明代谢产物“中枢”节点是否受到影响,这可能比“死端”代谢产物节点具有更多的生物学后果,从而证实主要是碳水化合物和甘油脂质代谢以及氧化还原平衡受到干扰(). 例如,在氧化还原稳态方面,老年小鼠的各种代谢物发生了变化,尤其是代谢物中心还原型谷胱甘肽(补充图S2在线)。
通过全球代谢组学在肝脏中进行老化生物标志物识别和途径分析。为了将衰老状态与肝代谢产物表达数据联系起来,我们使用了随机森林预测因子(A)。代谢产物根据年龄对组分离的重要性进行排名。我们发现,25个也是最重要的问题导致93%的明显预测准确率。彩色符号用于指示代谢物发挥作用的途径。监督层次聚类图(B),用于可视化RF预测因子的表达(z评分标准化)。RF图和监督层次聚类图表明,25种最重要的代谢物根据其老化状态对样品进行分层。(C) 代谢物集富集分析(MSEA)的汇总图(左面板),其中代谢物集合根据Holm p值进行排序,阴影线显示Holm p的截止值。代谢组视图(右侧面板)反映了根据中间中心性度量,x轴增加的代谢途径影响,这反映了代谢途径中的关键节点随着年龄的增长而发生了显著变化。
通过肌肉中的全球代谢组学进行衰老生物标志物识别和途径分析。为了将衰老状态与肌肉代谢产物表达数据联系起来,我们使用了随机森林预测因子(A)。代谢物根据其对按年龄分组的重要性进行排名。我们发现,这25个也是最重要的问题导致了100%的明显预测准确率。彩色符号用于指示代谢物发挥作用的途径。监督层次聚类图(B),用于可视化RF预测因子的表达(z评分标准化)。RF图和监督层次聚类图表明,25种最重要的代谢物根据其老化状态对样品进行分层。(C) 代谢产物集富集分析(MSEA)(左面板)和代谢组视图(右面板)的汇总图,如.
在肌肉中,前25位主要是脂肪酸代谢物()前20名已经显示出100%的预测准确率和明确的分组分离(). MSEA和代谢组观点证实了脂肪酸代谢的参与,但也强调了葡萄糖代谢的强大贡献,多不饱和脂肪酸的代谢亚麻酸和亚油酸是肌肉年龄的最佳预测因子(,补充图S3在线)。这些代谢物是多不饱和脂肪酸代谢的中枢,因此对预测能力有很大影响(补充图S3在线)。
讨论
老龄化及其相关疾病是社会的沉重负担。目前预期寿命的增加不仅影响到我们的社会系统,而且还伴随着常见慢性病的出现,包括神经系统、免疫系统和心血管系统的疾病,这些疾病往往达到流行病的程度。在这方面,重要的是要了解自然衰老过程,并阐明生活方式和/或药物干预可能产生的影响。近年来,人们提出了许多新的治疗方案来预防与衰老相关的疾病。尽管其中一些药物干预被证明可以延长寿命,即使是在哺乳动物中11,27由于尚未确定老化生物标记物,因此没有很好的数据显示健康寿命的改善。
在这里,我们描述了年轻和老年小鼠的临床以及分子、生化和代谢组的变化,从而定义了衰老的一致表型。除了血液和血浆中的代谢组学外,我们还采用了一种无偏见的方法来确定肌肉和肝脏中的代谢变化。我们使用C57BL/6J小鼠,因为这些小鼠特别适合研究代谢疾病28,29我们意识到,我们的研究可能会因伴随衰老的身体成分变化而产生偏差,但这一现象是所有物种衰老所固有的。虽然可以认为,没有脂肪量或葡萄糖稳态变化的衰老模型是最佳的,但这种模型并不存在,因为生理衰退是衰老的标志。
我们制作的自然衰老小鼠变化纲要强调了各种各样的分子和代谢变化。我们的数据表明线粒体功能受损的严重参与。我们不仅观察到与线粒体功能有关的各种转录物的表达显著下降,例如氧化磷酸化、FAO、凋亡,而且我们还测量了酶活性的变化,如4-HNE水平升高所证明的那样,可能最终导致氧化损伤的积累。这与我们在肝脏和肌肉代谢组学分析中确定氧化还原稳态及其中枢谷胱甘肽为受影响的途径相一致,表明活性氧物种在衰老表型的启动或进展中发挥作用。有趣的是,不同组织中基因表达和氧化磷酸化活性的降低似乎被线粒体生物发生和/或活性的补偿性增加所抵消,因为肌肉中的CS和BAT中的mtDNA丰度在老年小鼠中增加。
从信号传导的角度来看,胰岛素和mTOR途径受到了影响。IRS-1丝氨酸磷酸化增加与胰岛素抵抗有关30,这也是我们研究中观察到的衰老迹象(补充图S1在线)。IRS-1丝氨酸磷酸化阻碍其与胰岛素受体的适当结合,因此胰岛素促进IRS-1酪氨酸磷酸化并激活Akt和PI3-激酶及其下游通路的能力31,32,33与此一致,老年小鼠组织中S6K1磷酸化降低。S6K1磷酸化依赖于mTOR,mTOR由营养素和PI3-激酶协同激活34因此,IRS-1的丝氨酸磷酸化为mTOR信号传导减少提供了一个合理的解释。几种激酶可以诱导IRS-1丝氨酸磷酸化,包括JNK35、IKK36以及新型PKC36,37随着年龄增长,JNK和IKK的激活与慢性氧化和炎症应激有关。新的PKC,如PKC-θ,被二酰基甘油激活,当血清游离脂肪酸增加和线粒体功能缺陷时,二酰基丙三醇会积累,这两种情况都会随着年龄增长而再次出现。与之前的观察结果一致38,39,反映AMPK活化的基础AMPK磷酸化没有随着年龄的增长而改变,这表明在衰老组织中观察到的线粒体功能缺陷不会影响基础状态下ATP的生成。
我们研究的主要目标之一是识别新的衰老生物标志物。因此,我们进行了广泛的血液和血浆代谢组学研究,主要关注营养途径。在这种方法中,我们发现氨基酸和脂肪酸代谢物水平的显著变化可以准确预测衰老。除了血液代谢组学之外,我们还采用了一种无偏见的方法来确定衰老组织中的代谢变化。在肌肉和肝脏中,我们确定了在衰老时发生变化的代谢组,包括脂肪酸和葡萄糖代谢。
我们发现葡萄糖和糖酵解和糖原代谢的中间产物,如麦芽糖和麦芽四糖,是肝脏和肌肉衰老的生物标志物(——,). 在肝脏和肌肉中,糖原中间产物的积累表明老年小鼠的糖原代谢发生了改变,而乳酸增加和糖酵解中间产物减少则表明厌氧糖酵分解升高。在肌肉中,葡萄糖、葡萄糖-6-磷酸和麦芽糖水平的升高也表明葡萄糖和糖原代谢发生了变化(). 事实上,肌肉麦芽糖水平的增加与糖原分解的增加相一致。重要的是,这些代谢物数据与临床上的糖耐量受损和信号转导增加的IRS-1丝氨酸磷酸化数据相匹配。最近对糖尿病患者血液样本的代谢组学调查也发现麦芽糖水平升高40,再次支持我们发现的一些代谢变化与阿吉格相关的胰岛素抵抗有关。
衰老相关的糖和脂肪代谢变化模型。(A) 综述肌肉和肝脏中葡萄糖代谢的生化途径。年轻人和老年人之间变化的代谢物用粗体表示;肌肉和肝脏的代表性代谢物数据分别显示在左边缘和右边缘。总之,糖原中间产物麦芽糖和麦芽四糖的积累表明老年小鼠的肝糖原分解增加。另一方面,乳酸的增加和糖酵解中间产物的减少表明糖酵分解作用增强。肌肉中葡萄糖和葡萄糖-磷酸盐的升高,再加上乳酸水平的不变,表明肌肉中的葡萄糖处置减少,无论是通过分解还是糖原合成。事实上,肌肉麦芽糖水平的增加与糖原分解的增加是一致的。(B) 脂肪代谢的示意图,显示有代表性的变化。老年小鼠的血浆甘油三酯降低,而游离脂肪酸增加,表明组织中甘油三酸酯沉积增加,脂肪酸释放增强和/或分解减少。脂肪酸在肌肉中的积累、酰基肉碱向血浆中泄漏的减少(虚线)和酮生成减少支持了FAO的减少。肌肉中Cpt1b表达降低可能是导致FAO降低的原因之一,这也解释了血浆中乙酰肉碱水平降低的原因。老年小鼠的RER较低,表明其相对依赖脂肪氧化,这表明葡萄糖氧化也受到阻碍,脂肪代谢仍然是ATP生成的主要来源。数值表示为平均值±SEM或盒须图。*第页≤0.05**第页≤0.01; ***第页≤0.001.
在老年小鼠中明显受到影响的第二条途径是脂肪酸代谢(). 我们不仅观察到FAO相关基因的基因表达减少(Cpt1b型肌肉,阿卡德姆在肝脏,阿卡德尔BAT中),但来自该途径的代谢物也受到干扰。我们发现老年小鼠血浆游离脂肪酸增加,而血浆中链和长链酰肉碱降低。相反,众所周知,在禁食期间,血浆中的游离脂肪酸以及中长链酰基肉碱都会增加,这反映了FAO流量的增加41因此,我们的数据与线粒体FAO途径的进入和流量减少相吻合42,43,导致FA水平增加,同时减少了酰基肉碱“泄漏”回血浆室(中的虚线). 值得注意的是,在老年人中,几种脂肪酸在血浆中也更丰富22,尽管没有关于酰基肉碱的信息。我们的假设也与肌肉中脂肪酸的积累和Cpt1b型肌肉中的表达,表明线粒体膜上的肉碱酯化减少,从而减少了FAO进入线粒体的次数。这似乎与老年小鼠的RER降低这一事实相矛盾,这表明能量生成对脂肪的依赖性增加。然而,间接量热法测量全身代谢,反映所有代谢途径(如葡萄糖和脂肪氧化)和所有器官(如肝脏、肌肉等)中导致能量消耗的相对变化的组合。虽然我们不能排除FAO中组织特异性差异,但我们的数据表明脂肪代谢仍然是ATP生成的主要来源。
我们确定了各种组织中衰老的分子和生化特征,具有很高的预测准确性。我们的各种方法之间对这些途径的交叉验证——例如血液和组织代谢组学中的FAO代谢物,以及基因表达中的FAA基因——加强了它们作为生物标记物的潜在价值。然而,对于生物标记物的鉴定,我们建议分析血浆中的代谢物水平,因为它很容易获得,而无需牺牲动物或对人类具有太大的侵袭性。如上所述,血浆游离脂肪酸和酰基肉碱水平的变化对年轻和老年小鼠的分离具有很高的预测价值,并可能证明是衰老和健康寿命的一个强有力的候选生物标志物。值得注意的是,多不饱和脂肪酸的生物合成途径在老年肌肉中受到强烈影响,但在红细胞膜中也发现了类似的变化,再次证明了血液和血浆作为生物标记物来源的价值。此外,所选血浆氨基酸的测量增加了预测值。重要的是,最近证实了氨基酸作为生物标记物的相关性,因为包括BCAA在内的几种氨基酸的组合与胰岛素抵抗有关44并预测未来糖尿病45。我们在衰老的肝脏和肌肉中发现了BCAA代谢产物的一些变化,但在血浆中没有,这使得它们不太适合作为衰老的生物标志物。
尽管从临床角度来看可能不太实用,但一些额外的途径被证明对监测衰老过程和通过预防和治疗策略来延缓衰老和/或衰老相关病理是有用的,因为它们被证明同样可以预测衰老。这些包括增加组织氧化损伤和多不饱和脂肪酸中间产物的积累,特别是在肌肉中。此外,一组选定基因的mRNA表达(考克斯IV和Tfam公司)在所有受试组织中发生变化,可以预测衰老。与这些在大多数分析组织中发生变化的生物标记物相比,一些变化表明特定组织中存在老化,例如BAT中线粒体基因表达和mtDNA含量的增加,可能使该组织更好地适应老化的代谢挑战。
总之,我们的数据提供了一个控制良好的哺乳动物系统中代谢衰老的足迹。我们通过综合方法确定的途径也开辟了一条新的研究途径,旨在确定这些分子和代谢组系统对衰老过程的个人贡献。虽然其中一些通路已经被研究与衰老有关,但还需要进一步研究以充分确定其相关性,也包括在人类中。未来的研究还必须解决生理挑战,如禁食/喂养、运动或老年小鼠的高脂肪喂养对这些生物标记物的影响。此外,我们希望这些老化生物标志物的识别可以刺激合理设计策略,以防止或扭转某些受影响网络的衰退,从而延长健康寿命并防止加速老化。
方法
老鼠
从Janvier(法国圣贝切芬)购买13周龄(“年轻”)或93周龄((“老年”)的C57BL/6J小鼠。小鼠被分组饲养,并接受标准食物(2018号,含18%蛋白质、50%碳水化合物和6.0%脂肪;Harlan Laboratories,Madison,WI,USA)。在体内表型分析后,小鼠在通宵禁食后被处死,年龄达到24或103周。取肝素化血浆,将组织冷冻在液氮中进行生化和分子分析。
体内表型分析
大多数临床试验都是根据优吗啡项目中建立和验证的严格标准操作规程进行的28,46身体成分由Echo-MRI(Echo Medical Systems,Houston,TX,USA)和耗氧量(VO)测定2)使用综合实验动物监测系统(CLAMS)通过间接量热法监测呼吸交换比(RER)、活性水平和活性(美国俄亥俄州哥伦布市哥伦布仪器公司)47隔夜禁食后腹腔注射2 g/kg葡萄糖,通过测量血糖和胰岛素来分析葡萄糖耐量48.最大呼吸潜力(VO2max)是在代谢跑步机(美国俄亥俄州哥伦布市哥伦布仪器公司)中采用递增速度方案测定的。在跑步过程中,VO2和VCO2已测量。当小鼠精疲力竭、RER超过一个时间间隔大于1或VO时,实验停止2到达一个清晰的高原。将一个记录跑步的轮子放在家里的笼子里24小时来记录志愿运动。在此期间,老鼠被单独关在笼子里。使用CatWalk系统进行步态分析(Noldus Technology,Wageningen,荷兰)49。将小鼠放置在一个方形的开阔场地中30分钟,并使用视频跟踪记录它们的运动(荷兰瓦格宁根Noldus Technology公司的Ethovision XT)50.
组织学和生物化学分析
线粒体复合体I和IV的活性(美国犹吉尼亚州线粒体研究所)、柠檬酸合成酶(西格玛)根据制造商的协议进行测定。使用基于ELISA的方法(Cell Biolabs,San Diego,CA,USA)在PBS组织裂解液中测定4-HNE。如前所述进行蛋白质印迹2其他抗体包括Tubulin(Santa Cruz Biotechnology)、p-Ser(Sigma)、p-ACC(Millipore)、IRS-1、p-AMPKα、AMPKβ、S6K1和p-S6K1(Cell Signaling Technology)51.
基因表达与线粒体DNA丰度
使用Trizol从组织中提取RNA。使用QuantiTect RT试剂盒(Qiagen)合成cDNA,然后在罗氏Lightcycler 480上进行qPCR,并使用罗氏SYBR绿色母料。qPCR引物在扩展实验程序中进行了描述。线粒体DNA丰度(mtDNA)如前所述进行量化47简而言之,使用标准苯酚提取法提取DNA,并使用扩展实验程序中的引物将所得基因组DNA(gDNA)用于qPCR。来自5个月(幼)、25个月(大)或CR饮食的C57BL/6N小鼠腓肠肌的微阵列数据21,检索自公共基因老化Nexus数据库(http://gan.usc.edu/public/index.jsp)52.
血浆生物化学和靶向代谢组学
根据优吗啡SOP进行基本血浆临床化学分析28,46用串联质谱法(MS)测定血浆中的酰肉碱和氨基酸n个)和LC-MSn个分别地53,54.红细胞中的脂肪酸直接酯交换,并通过气相色谱和火焰离子化检测进行分析55使用既定程序在高氯酸脱蛋白血浆中测定β-羟基丁酸56
全球代谢组学
Metabolon(美国北卡罗来纳州达勒姆)根据公布的方法对年轻和老年小鼠(每组分别为8只和7只)的股四头肌和肝脏进行代谢组学研究57简而言之,样品制备是使用一系列专有的有机和水萃取液进行的,以去除蛋白质部分,同时最大限度地回收小分子。将提取的样品分成等份,在GC/MS和LC-MS/MS平台上进行分析。对于LC-MS/MS,样品分为两等分样品,在正(酸性溶剂)或负(碱性溶剂)电离模式下进行分析。在5%苯基GC柱中对双三甲基硅酰三氟乙酰胺衍生物样品进行GC-MS。
生物信息学和统计学
图中的数据以平均值±SEM或箱须图的形式表示,表示样本最小值、下四分位数、中位数、上四分位数和样本最大值,异常值用小圆圈表示。进行学生T检验或韦尔奇T检验(肝脏和肌肉代谢组学数据),对年轻和老年小鼠进行统计比较。A类第页-0.05或更低的值被认为是显著的。
使用Genepattern中实现的软件,使用完全连锁和Pearson秩相关距离对归一化代谢物进行无监督的层次聚类(http://www.broadinstitute.org/cancer/software/genepattern/)58,59通过从每个值中减去每个代谢物的平均表达值,然后将结果值除以标准偏差来计算z评分。热图中的颜色反映了代谢物的相对丰度水平,红色高于平均代谢物丰度值,蓝色低于平均代谢品丰度值。代谢物和动物的排序使用距离函数1-相关性在层次聚类中确定。
基于web的代谢组学数据处理工具MetaboAnalyst60用于组织代谢物数据分析。请参见http://www.metaboanalyst.ca详细方法。简言之,对于RF分析,代谢物数据进行对数转换(log 2),每个节点尝试的预测因子数量设置为每个组织变量总数的平方根,待生长的树木数量设置为10.000。使用随机森林包执行RF,并显示由此产生的混淆矩阵、袋外(OOB)误差和变量重要性(由数据置换时精度的平均降低确定)。对于MSEA,代谢物数据根据HMDB绘制,并选择“代谢物途径相关代谢物集”库(目前有88个条目)进行富集分析,这是使用包装全局测试进行的。代谢组视图显示了路径拓扑分析,这是一种在确定哪些路径更有意义时考虑路径结构的方法。为该分析选择的参数包括全局测试和节点中心性度量“中间中心性”,以估计节点重要性。该度量考虑了全局网络结构,使得出现在两个密集集群之间的节点将具有高介数中心性,即使其度中心性不高。根据中心性度量计算的节点重要性值通过路径的重要性之和进一步归一化,以在不同路径之间进行比较。因此,每条路径的总/最大重要性为1;每个代谢节点的重要性度量实际上是相对于总路径重要性的百分比,而路径影响是匹配代谢节点的累积百分比。
作者贡献
RHH和JA在CA、SMH、CC和KS的投入下构思和设计了该项目。RH、CA、SMH、CC、EHJ和PAA进行了实验。CT和RD计划并监督小鼠实验。RHH、CA、SMH和JA撰写了这份手稿,所有作者都提供了稿件。
致谢
我们感谢Auwerx实验室和小鼠表型单位的团队成员提供的技术支持和有益的讨论。RHH得到了荷兰科学研究组织Rubicon奖学金的支持。JA是雀巢能量代谢主席。作者的实验室工作得到了洛桑理工学院、欧盟创意项目(ERC-2008-AdG-)、Velux基金会、瑞士国家科学基金会(31003A-124713和31003A-125487)和荷兰科学研究组织(NWO VIDI No.016.086.336 to SMH)的资助。
工具书类
- Houtkooper R.H.、Williams R.W.和Auwerx J。长寿的代谢网络.单元格 142, 9–14 (2010).[PMC免费文章][公共医学][谷歌学者]
- 白P。等。.PARP-1抑制通过SIRT1激活增加线粒体代谢.单元格元 13, 461–468 (2011).[PMC免费文章][公共医学][谷歌学者]
- 萨欣E。等。.端粒功能障碍导致代谢和线粒体损伤.自然 470, 359–365 (2011).[PMC免费文章][公共医学][谷歌学者]
- Fontana L.、Partridge L.和Longo V.D。延长健康寿命——从酵母到人类.科学类 328, 321–326 (2010).[PMC免费文章][公共医学][谷歌学者]
- 凯尼恩·C·J。衰老的遗传学.自然 464, 504–512 (2010). [公共医学][谷歌学者]
- Lin M.T.和Beal M.F。神经退行性疾病中的线粒体功能障碍和氧化应激.自然 443, 787–795 (2006). [公共医学][谷歌学者]
- 罗素·S·J·和卡恩·C·R·。衰老的内分泌调节.Nat Rev Mol细胞生物学 8, 681–691 (2007). [公共医学][谷歌学者]
- Wullschleger S.、Loewith R.和Hall M.N。生长和代谢中的TOR信号.单元格 124, 471–484 (2006). [公共医学][谷歌学者]
- Zoncu R.、Efeyan A.和Sabatini D.M。mTOR:从生长信号整合到癌症、糖尿病和衰老.Nat Rev Mol细胞生物学 12, 21–35 (2011).[PMC免费文章][公共医学][谷歌学者]
- Sengupta S.、Peterson T.R.、Laplante M.、Oh S.和Sabatini D.M。mTORC1控制禁食诱导的酮生成及其随年龄的调节.自然 468, 1100–1104 (2010). [公共医学][谷歌学者]
- 哈里森·D·E。等。.在晚年喂食雷帕霉素可延长基因异质小鼠的寿命.自然 460,392-395(2009年)。[PMC免费文章][公共医学][谷歌学者]
- Haigis M.C.和Sinclair D.A。哺乳动物sirtuins:生物学见解和疾病相关性.病理学年度回顾 5, 253–295 (2010).[PMC免费文章][公共医学][谷歌学者]
- 担保人L。线粒体——衰老、热量限制和sirtuins之间的联系? 单元格 132, 171–176 (2008).[PMC免费文章][公共医学][谷歌学者]
- Canto C.和Auwerx J。PGC-1alpha、SIRT1和AMPK,一个控制能源消耗的能源传感网络.当前Opin Lipidol 20, 98–105 (2009).[PMC免费文章][公共医学][谷歌学者]
- 迈尔·W。等。.AMPK和钙调神经磷酸酶诱导的寿命延长由CRTC-1和CREB介导.自然 470, 404–408 (2011).[PMC免费文章][公共医学][谷歌学者]
- 富尔科·M。等。.葡萄糖限制通过AMPK介导的Nampt调节激活SIRT1抑制骨骼肌成肌细胞分化.开发人员单元格 14, 661–673 (2008).[PMC免费文章][公共医学][谷歌学者]
- 第C章。等。.AMPK通过调节NAD+代谢和SIRT1活性调节能量消耗.自然 458, 1056–1060 (2009).[PMC免费文章][公共医学][谷歌学者]
- 第C章。等。.AMPK和SIRT1在骨骼肌禁食和运动代谢适应中的相互依赖性.单元格元 11, 213–219 (2010).[PMC免费文章][公共医学][谷歌学者]
- Durieux J.、Wolff S.和Dillin A。电子传递链介导寿命的细胞非自主性.单元格 144, 79–91 (2011).[PMC免费文章][公共医学][谷歌学者]
- Lemieux H.、Vazquez E.J.、Fujioka H.和Hoppel C.L。大鼠心脏通透性纤维线粒体功能下降与衰老表型相关.老年医学杂志 65, 1157–1164 (2010).[PMC免费文章][公共医学][谷歌学者]
- 爱德华兹M.G。等。.衰老的基因表达谱揭示了p53介导的转录程序的激活.BMC基因组学 8,80(2007年)。[PMC免费文章][公共医学][谷歌学者]
- 劳顿K.A。等。.成人血浆代谢组分析.药物基因组学 9, 383–397 (2008). [公共医学][谷歌学者]
- D’Antona G。等。.补充支链氨基酸可促进中年小鼠的存活并支持心肌和骨骼肌线粒体的生物生成.单元格元 12,362–372(2010年)。[公共医学][谷歌学者]
- Grandison R.C.、Piper M.D.和Partridge L。氨基酸失衡解释了果蝇通过饮食限制延长寿命的原因.自然 462, 1061–1064 (2009).[PMC免费文章][公共医学][谷歌学者]
- 米勒R.A。等。.缺乏蛋氨酸的饮食延长小鼠寿命,减缓免疫和晶状体老化,改变葡萄糖、T4、IGF-I和胰岛素水平,增加肝细胞MIF水平和抗应激能力.老化细胞 4, 119–125 (2005).[PMC免费文章][公共医学][谷歌学者]
- Xia J.&Wishart D.S。MSEA:一种基于网络的工具,用于识别定量代谢组数据中具有生物学意义的模式.核酸研究 38第71–77页(2010年)。[PMC免费文章][公共医学][谷歌学者]
- 皮尔逊·K·J。等。.白藜芦醇可延缓与年龄相关的恶化,并在不延长寿命的情况下模拟饮食限制的转录方面.单元格元 8, 157–168 (2008).[PMC免费文章][公共医学][谷歌学者]
- 香槟M.F。等。.遗传背景决定小鼠的代谢表型.哺乳动物基因组 19, 318–331 (2008). [公共医学][谷歌学者]
- 投标人S.B。等。.饮食和遗传背景对固醇调节元件结合蛋白-1c、硬脂酰辅酶A去饱和酶1和代谢综合征发生的影响.糖尿病 54, 1314–1323 (2005). [公共医学][谷歌学者]
- Morino K.、Petersen K.F.和Shulman G.I。人胰岛素抵抗的分子机制及其与线粒体功能障碍的潜在联系.糖尿病 55补充2,S9–S15(2006)。[PMC免费文章][公共医学][谷歌学者]
- 帕兹·K。等。.胰岛素抵抗的分子基础。IRS-1和IRS-2的丝氨酸/苏氨酸磷酸化升高抑制了它们与胰岛素受体膜旁区域的结合,并损害了它们接受胰岛素诱导的酪氨酸磷酸化的能力.生物化学杂志 272, 29911–29918 (1997). [公共医学][谷歌学者]
- 阿吉雷五世。等。.胰岛素受体底物-1中Ser307的磷酸化阻断与胰岛素受体的相互作用并抑制胰岛素作用.生物化学杂志 277,1531-1537(2002年)。[公共医学][谷歌学者]
- Moeschel K.公司。等。.蛋白激酶C-zeta诱导胰岛素受体底物-1(IRS-1)中Ser318的磷酸化减弱与胰岛素受体的相互作用和IRS-1的酪氨酸磷酸化.生物化学杂志 279, 25157–25163 (2004). [公共医学][谷歌学者]
- Pearce L.R.、Komander D.和Alessi D.R。AGC蛋白激酶的基本原理.Nat Rev Mol细胞生物学 11, 9–22 (2010). [公共医学][谷歌学者]
- 鲁伊·L。等。.胰岛素/IGF-1和TNF-α通过不同途径刺激抑制性Ser307的IRS-1磷酸化.临床研究杂志 107, 181–189 (2001).[PMC免费文章][公共医学][谷歌学者]
- 高Z。等。.抑制剂κB激酶复合物对胰岛素受体底物1的丝氨酸磷酸化.生物化学杂志 277, 48115–48121 (2002). [公共医学][谷歌学者]
- 李毅(Li Y.)。等。.蛋白激酶C Theta通过磷酸化Ser(1101)的IRS1抑制胰岛素信号.生物化学杂志 279, 45304–45307 (2004). [公共医学][谷歌学者]
- 雷兹尼克·R·M。等。.Aging相关的AMP活化蛋白激酶活性降低和线粒体生物发生.单元格元 5, 151–156 (2007).[PMC免费文章][公共医学][谷歌学者]
- Gonzalez A.A.、Kumar R.、Mulligan J.D.、Davis A.J.和Saupe K.W。衰老对小鼠心肌和骨骼肌AMPK活性的影响:基础活性、变构激活和对体内低氧血症的反应.美国生理学杂志Regul Integr Comp Physiol 287,R1270–1275(2004)。[公共医学][谷歌学者]
- 苏赫雷·K。等。.糖尿病的代谢足迹:流行病学背景下的多平台代谢组学研究.公共科学图书馆 5,e13953(2010)。[PMC免费文章][公共医学][谷歌学者]
- Costa C.C.、de Almeida I.T.、Jakobs C.、Poll-The B.T.和Duran M。禁食和葵花籽油激发试验诱导儿童血浆乙酰肉碱水平的动态变化.儿科研究 46, 440–444 (1999). [公共医学][谷歌学者]
- Tucker M.Z.和Turcotte L.P。基础条件下灌注老龄大鼠肌肉脂肪酸氧化损伤.美国生理内分泌代谢杂志 282,E1102–1109(2002)。[公共医学][谷歌学者]
- Hyyti O.M.、Ledee D.、Ning X.H.、Ge M.和Portman M.A。衰老损害小鼠心肌脂肪酸和酮类氧化并改变心脏功能和代谢对胰岛素的反应.美国生理学杂志心脏循环生理学 299,H868–875(2010)。[PMC免费文章][公共医学][谷歌学者]
- 纽加德·C.B。等。.一种分支链氨基酸相关代谢特征,可区分肥胖和瘦人并导致胰岛素抵抗.单元格元 9, 311–326 (2009).[PMC免费文章][公共医学][谷歌学者]
- 王天杰。等。.代谢谱与糖尿病风险.自然·医学 17,448–453(2011年)。[PMC免费文章][公共医学][谷歌学者]
- Champy M.F.公司。等。.小鼠功能基因组学要求小鼠操作和居住条件标准化.哺乳动物基因组 15, 768–783 (2004). [公共医学][谷歌学者]
- 拉古格M。等。.白藜芦醇通过激活SIRT1和PGC-1α改善线粒体功能并预防代谢性疾病.单元格 127, 1109–1122 (2006). [公共医学][谷歌学者]
- Heikkinen S.、Argmann C.A.、Champy M.F.和Auwerx J。葡萄糖稳态的评估.电流质子分子生物学 第章 29第29B 23单元(2007年)。[公共医学][谷歌学者]
- Angeby-Moller K.、Berge O.G.和Hamers F.P。使用CatWalk方法评估角叉菜胶诱导的踝关节单关节炎行走大鼠的负重和疼痛行为:吗啡和罗非昔布的影响.神经科学方法杂志 174, 1–9 (2008). [公共医学][谷歌学者]
- 曼迪略S。等。.小鼠行为表型的可靠性、稳健性和再现性:一项跨实验室研究.基因组学杂志 34, 243–255 (2008).[PMC免费文章][公共医学][谷歌学者]
- 嗯S.H。等。.不含S6K1可以防止年龄和饮食诱导的肥胖,同时提高胰岛素敏感性.自然 431, 200–205 (2004). [公共医学][谷歌学者]
- 平移F。等。.基因老化Nexus:用于老化微阵列数据的网络数据库和数据挖掘平台.核酸研究 35,D756–759(2007)。[PMC免费文章][公共医学][谷歌学者]
- Chace D.H.、Kalas T.A.和Naylor E.W。串联质谱法在新生儿干血样多分析物筛选中的应用.临床化学 49, 1797–1817 (2003). [公共医学][谷歌学者]
- Piraud M。等。.未激活氨基酸的ESI-MS/MS分析:一种诊断氨基酸代谢遗传性疾病的新工具。正负电离模式下79个生物相关分子的碎片化研究.快速公共质谱 17, 1297–1311 (2003). [公共医学][谷歌学者]
- Dacremont G.和Vincent G。红细胞和成纤维细胞中血浆白蛋白和多不饱和脂肪酸(PUFA)的测定.J继承元疾病 18补遗1,84-89(1995年)。[公共医学][谷歌学者]
- Bergmeyer H.U.、Bergmeyer J.和Grassl M。酶分析方法(VCH Verlagsgesellschaft mbH,1986)。[谷歌学者]
- 加尔·W·E。等。.α-羟基丁酸是非糖尿病人群胰岛素抵抗和葡萄糖不耐受的早期生物标志物.公共科学图书馆 5,e10883(2010)。[PMC免费文章][公共医学][谷歌学者]
- Reich M。等。.基因模式2.0.自然基因 38, 500–501 (2006). [公共医学][谷歌学者]
- de Hoon M.J.、Imoto S.、Nolan J.和Miyano S。开源集群软件.生物信息学 20, 1453–1454 (2004). [公共医学][谷歌学者]
- Xia J.、Psychogios N.、Young N.和Wishart D.S。MetaboAnalyst:用于代谢组学数据分析和解释的web服务器.核酸研究 37,W652-660(2009年)。[PMC免费文章][公共医学][谷歌学者]