慢性心功能衰竭患者尿液的代谢组学研究.docx

慢性心功能衰竭患者尿液的代谢组学研究 庞博+王斌+舒振波+姚小晓+张国栋+扈聪+吴绥生摘 要 采用基于液相色谱-质谱联用的方法对慢性心力衰竭（Chronic heart failure， CHF）患者和正常对照（Control）人群的尿液进行分析， 筛选慢性心力衰竭患者尿液中的差异代谢物， 研究其发病机制， 并为临床治疗提供科学依据选择15个慢性心力衰竭患者（年龄（62.273.14）岁）及15个正常人（年龄（65.414.63）岁）， 采用高分辨度快速液相色谱-四极杆-飞行时间串联质谱（RRLC-QTOF/MS）技术对尿液代谢物进行分析， 采用主成分分析（PCA）对两组代谢物进行分类， 并筛选潜在生物标记物； 运用偏最小二乘判别分析法（PLS-DA）建模， 考察生物标记物对疾病筛选的预测能力研究结果表明， CHF组和Control组尿液代谢物谱能得到很好的区分， 发现并鉴定了2种潜在生物标记物尿苷及丙氨酰色氨酸， 提示嘧啶代谢和色氨酸代谢可能在心力衰竭发生发展中有重要作用关键词 心力衰竭； 快速高分辨液相色谱-質谱； 尿苷； 丙氨酰色氨酸； 代谢组1 引 言心力衰竭（Heart failure， HF）是指由于心脏器质性或功能性疾病而导致心室充盈和射血能力的降低， 不能满足组织代谢需要， 从而引发的复杂临床综合征， 是冠心病、高血压、心肌梗死等多种心血管疾病的终末阶段。

据美国心脏学会统计， 在过去30年， 心力衰竭的发生率及住院率增加了约3倍， 在65岁以上人群中排名首位[1]心力衰竭的发生发展涉及代谢、神经激素、免疫系统等多个病理生理通路[2]生物标志物一般是指可供客观测评病理生理状态下某种特征性的生化指标， 通过测定生物标志物可以获知机体所处的生物学状态和疾病的进程， 并为诊断、预后及治疗提供依据[3]目前， 已发现许多心力衰竭相关的重要的生物标记物， 如脑钠肽（BNP）[4]、肿瘤因子TNF、白介素IL-1β、IL-6、IL-18、IL-33 及其受体sTNFR1 和sTNFR2等[5]， 但经验证并在临床上大范围使用的只有BNP， 虽然其它标记物浓度在心力衰竭患者中确实有一定程度的升高， 但因特异性和灵敏度均不理想， 尚未进入临床应用目前， 寻找心力衰竭发生发展过程中新的生物标志物， 发现新的治疗靶点， 从而帮助心力衰竭患者风险分层及改善预后是心力衰竭的重要研究方向之一代谢组学是考察生物体系受刺激或扰动后（例如某个基因变异或某一病理生理状态下）其代谢产物-内源性代谢物质种类、数量及其变化规律的科学[6～10]可以通过核磁共振、色谱-质谱联用技术等进行代谢物数据采集。

核磁共振技术对样本分析选择性高且无损伤性， 而液相色谱-质谱技术拥有更高的灵敏度和分辨率随着快速高分辨液相色谱（RRLC）的发展，其更高灵敏度、分辨率及分离度的特点使其更适合于代谢组学的研究[11]本研究采用高分辨快速液相色谱-四极杆-飞行时间串联质谱联用技术（Rapid resolution liquid chromatography-quadrupole time-of-flight mass spectrometry， RRLC-QTOF/MS）对慢性心力衰竭（Chronic heart failure， CHF）患者的尿液进行代谢组学分析， 研究心力衰竭的特征性代谢谱， 从中筛选代表性生物标志物， 以期能为心力衰竭提供新的诊断策略， 并对分析心力衰竭的预后及治疗提供新的依据， 为进一步研究心力衰竭的发生机制提供参考2 实验部分2.1 仪器与试剂Agilent1200快速分离液相色谱系统（美国Agilent公司）， 配备高压二元泵， 控温自动进样器； Agilent 6520 Q-TOF质谱仪（美国Agilent公司）， 配有ESI 离子源； Eppendorf 5804R 高速离心机（德国Eppendorf公司）； Milli-Q超纯水仪（美国Millipore公司）。

乙腈和甲酸（色谱纯， 美国TEDIA公司）； 其它试剂至少为分析纯； 实验用水为Milli-Q超纯水（18.2 MΩcm）2.2 样本来源慢性心力衰竭患者尿样和正常人尿样均由吉林大学附属第一医院提供研究对象选取我院心血管内科于2015年10～11月收治的CHF患者15例入选标准： （1）均符合2014年中国心力衰竭诊断和治疗指南中心力衰竭诊断标准[12]2）按美国纽约心脏病协会（New York Heart Disease Association， NYHA）心脏功能分级Ⅲ到Ⅳ级， 其中Ⅲ级组6例、Ⅳ级组9例排除标准： （1） 非心力衰竭主因入院者； （2）院内死亡者； （3）除冠心病、高心病、扩心病之外其它原因（血液系统， 慢性贫血性心脏病， 内分泌及代谢原因如甲亢性心脏病等）导致慢性心力衰竭者； （4）合并严重的需外科手术治疗的先天性心脏病者（如瓣膜病）； （5）肝肾功能衰竭者； （6）患有其它系统严重疾病如恶性肿瘤、出血性疾病等、严重的急性感染或代谢紊乱者； （7）由于意识障碍、智力低下或者言语含糊不能清晰、准确描述病史者选取同一时期在本院体检中心体检健康者15例为对照组， 入选标准： 既往身体健康， 无重大疾病史， 以及无近期感染史， 无吸毒及酗酒史。

样本的一般情况见表1此次研究已取得患者同意， 签订知情同意书， 经医院伦理委员会通过所有样本均由清晨空腹尿样， 采集后3000 r/min离心10 min， 上清液于80℃储存备用2.3 样本处理将样本从80℃取出， 室温下溶解， 采用离心沉降后四倍稀释法[13]， 在4℃以12000 r/min 离心5 min， 取50 μL上清液， 用水稀释至200 μL， 以0.45 μm滤膜过滤， 待LC-MS分析 2.4 分析条件2.4.1 色谱条件 Agilent Eclipse Plus C18色谱柱（150 mm 2.1 mm， 3.5 μm， 美国Agilent公司）； 柱温： 30℃； 流动相A： 超纯水（含0.1%（V/V）甲酸）； 流动相B： 乙腈； 梯度洗脱条件： 0～5 min， 10%～30% B； 5～10 min， 30%～80% B； 10～12 min， 80%～85% B； 12～15 min， 85%～95% B； 15～20 min， 95% B； 20～21 min， 95%～100% B； 流速： 0.4 mL/min2.4.2 质谱条件 分别采用电喷雾正离子模式（ESI+）和负离子模式（ESI-）对样品进行检测； 质量扫描范围m/z 100～1200， 干燥气流量（N2）为9 L/min， 干燥气温度为 250℃， 雾化电压为0.2758 MPa， 毛细管电压为3.5 kV， 碎裂电压为150 V， 锥孔电压为65 V。

在样品测试之前， 使用调谐液校正质量轴2.5 数据处理尿液样品用RRLC/MS进行检测， 得到样品的总离子流色谱图， 在Mass Hunter 软件的分子特征识别模式下， 进行峰校准、背景扣除、面积归一化和数据简化处理， 将结果转化为包含化合物的保留时间和质荷比信息的CEF格式文件将CEF格式文件导入Mass Profiler Professional （MPP， Agilent Technologies， USA）软件， 统计软件进行滤噪和归一化多变量数据进行主成分分析（Principal component analysis， PCA）， 对原始数据进行数据压缩， 将多维的数据压缩成几个主要的主成分（PCs）来描述数据内部的特征， 通过PCA 得分图（Score plot）获得样本分类信息， 通过PCA 载荷图（Loading plot）发现可作为生物标记物的化合物并予以鉴定运用偏最小二乘判别分析（Partial least square-discriminant analysis， PLS-DA）法建立預测模型， 选取基于PCA分析得出的特征项列表， 以原样本的2/3建立训练集， 剩余的1/3作为验证集， 考察生物标记物对慢性心力衰竭筛选的准确度。

鉴定潜在生物标时使用以下数据库[11]： HMDB （ http： //www. hmdb. ca/）； METLIN （ http： //metlin.scripps. edu/）3 结果与讨论3.1 慢性心力衰竭患者和正常人尿液LC-MS总离子流图尿液样品经RRLC-QTOF/MS分别在正负离子模式下进行检测分析图1为随机选取的一个正常人和一个慢性心力衰竭患者的尿液的总离子流图（Total ion current， TIC）， 可见在正负离子模式下， 二者的色谱峰存在明显差别， 推测二者尿液代谢产物中存在一些浓度差异的物质3.2 数据分析3.2.1 主成分分析（PCA） 采用主成分分析考察慢性心力衰竭患者和正常人尿样代谢物变化如图2所示， 在正离子模式和负离子模式下的PCA得分图能够很好的区分慢性心力衰竭患者和正常人尿样代谢物图谱图中每个点代表一个样本在二维平面上的投影， 每个样本的位置由其自身的代谢决定， 处于相同生理病理状态的样本通常具有相似的代谢物组成， 因此在得分图上也处于相似的位置， 彼此之间距离越远， 表示其生理病理状态相差越大3.2.2 差量变异研究 基于慢性心力衰竭患者和正常人的PCA分析结果， 确定两组之间存在代谢差异， 这些变量的PCA载荷图如图3所示， 图中每个点代表样本中检测到的代谢物相关信息， 距离原点越远表明该代谢物对两组间分类的贡献越大。

根据载荷图的结果， 选取其中距离远点比较远， 即对分组贡献较大的代谢物作为潜在生物标记物由于代谢组学考察的内源性代谢物的分子量在1000 Da范围内， 而质谱的质量精度在1000 Da 时误差为1010-6， 所以在使用数据库考察潜在生物标记物时应选择测量值与实际值分子量误差差距小于10106的化合物3.2.3 潜在标志物的鉴定 根据精确分子量和串联质谱结果以及与数据库或标准品的比较来进行潜在生物标记物的鉴定以负离子模式下的离子m/z 244.0674为例来说明生物标记物的鉴定过程负离子模式下该离子的提取离子色谱图及在保留时间1.602 min时的质谱图如图4A和图4B所示将精确分子量与数据库中化合物进行比对， 初步鉴定该化合物为尿苷如图4C所示， 串联质谱的结果也证明了这一点， 与标准品的串联质谱分析结果完全一致， 此化合物被鉴定为尿苷以同样方法鉴定其它化合物， 但是一些化合物未能鉴定出表2列出了潜在生物标记物的鉴定结果通过ANOVA分析得到其在健康对照组、CHF组之间的变化趋势3.2.4 尿样生物标记物的验证 采用偏最小二乘判别分析（PLS-DA）对3.2.2中发现的潜在生物标记物进行建模， 验证它们在慢性心力衰竭诊断上的预测能力。

通过计算得出标记物对慢性心力衰竭预测能力的正确率， 正负离子模式下分别为94.5%和93.1%， 具有较好的预测能力3.3 慢性心力衰竭潜在生物标志物分析本研究采用的RRLC-QTOF/MS方法， 具有灵敏度高、分辨率高及分离度好的特点， 近年来广泛用于基础和临床的代谢组学研究[15，16]通过此方法对慢性心力衰竭患者和健康人群尿液进行检测， 结合统计学处理发现并鉴定了7个潜在的生物标志物及其变化趋势， 研究结果表明， 尿苷和丙氨酰色氨酸的表达上调趋势可能与慢性心力衰竭的发生发展有关 1993年， Smoleński[17]和Yaoita等[18]提出缺血心肌嘧啶代谢异常， 并在心肌缺血大鼠模型发现尿苷的异常变化， 推测尿苷可能与缺血后血管累及的心肌细胞的生存有关2006年， Krylova等[19]在研究中发现， 尿苷和尿嘧啶可以通过。