结核分枝杆菌潜伏性感染诊断的潜在标志物蛋白研究

张海晴,张雪迪,刘成永,周冬青,魏素梅,贾 彤, 张礼茂,郭 岩,李 莉,黄海滨

(江苏省徐州市传染病医院 结核二科,江苏 徐州,221004)

结核病是由结核分枝杆菌感染所致，有潜伏期及活动期2种疾病状态[1-2],潜伏性感染患者早期无明显症状[3],因此早期发现和治疗潜伏性感染患者对控制结核病的传播意义重大。目前，检测结核潜伏性感染的方法仍不完善，结核菌素试验不能明确区分卡介苗接种者、非结核分枝杆菌感染者、结核潜伏性感染和活动性结核病者,γ-干扰素 (IFN-γ)释放试验有局限性，亦不能鉴别潜伏性感染及活动性结核病[4]。因此，寻求结核潜伏性感染更为有效的诊断方法至关重要。潜伏性感染及活动性结核病患者中结核分枝杆菌蛋白表达差异，可作为潜伏性感染诊断的新思路。本研究通过同位素标记相对和绝对定量(iTRAQ)技术对结核分枝杆菌的休眠株及标准株进行蛋白质组差异表达筛选，并结合生物信息学技术进行分析，获得参与潜伏性感染菌株中直接与休眠相关的能量代谢通路中的蛋白P9WPC9,从而为筛选出可用于结核潜伏性感染诊断的潜在标志物提供理论依据。

1 材料与方法

1.1 材料与试剂

选取结核分枝杆菌标准株H37Rv和由标准株诱导的潜伏性感染菌株模型(休眠株)各三皿，分别编号为sta1、sta2、sta3;dor1、dor2、dor3,数量为1×107;其中sta3和dor3因样本蛋白降解，所以未能检测出相关数据。试剂为蛋白酶抑制剂(protease inhibitor,美国 Calbiochem公司)、胰酶(trypsinPromega,德国Merk公司)、乙腈(acetonitrile、超纯水，德国 Fisher Chemical公司)、三氟乙酸(trifluoroacetic acid)、二硫苏糖醇(dithiothreitol)、三乙基碳酸氢铵(TEAB)、尿素、碘代乙酰胺(iodoacetamide)(三氟乙酸、二硫苏糖醇、TEAB、尿素、碘代乙酰胺均购自美国 Sigma-Aldrich Sigma公司)。

1.2 实验方法

(1)菌株造模。① 标准菌种培养:将结核分枝杆菌标准株H37Rv接种于提前配置好的中性罗氏培养基(珠海贝索生物技术有限公司)上，培养3～4周后将菌株转种于含0.2%甘油、10% ADC的7H9液体培养基(BI)上，并置于摇床中震荡培养，温度为37 ℃,160转/min,12～15 h待达到对数生长期即可(吸光度值约为0.4)。② 休眠菌种模型构建:取已经达到对数生长期的标准菌株进行研磨，研磨均匀后，按照菌株与培养基1∶100进行配液，加入亚甲蓝，终浓度为1.5 μg/mL。将此混合液制成菌悬液，吸取10 mL于15 mL试管中，上层气体与下层液体的体积比为1∶2,封口胶密封后置于37 ℃恒温箱静止培养。培养12～15 d后，溶液中的亚甲蓝完全处于无色状态时，说明菌株已进入完全缺氧休眠状态。

(2)iTRAQ技术。① 菌株样本蛋白抽提:向不同菌株样本中分别加入4倍体积(尿素、1%蛋白酶抑制剂)裂解缓冲液，然后进行超声裂解。将裂解后的混合液至于4 ℃离心机中离心10 min,转速为12 000 g。去除细胞碎片后，将上清液转移至新的离心管中，利用BCA试剂盒进行蛋白浓度测定。② 样品及实验预判。③ 蛋白质经过还原及巯基封闭后，使用胰蛋白酶进行酶切。④ 肽段混合物分别使用4种不同的iTRAQ试剂进行标记，使用Strata X C18(Phenomenex)将胰酶酶解后肽段进行盐份祛除，待进行真空干燥及冷冻处理后，按照0.5 M TEAB将所获得的肽段全部进行溶解处理。根据iTRAQ操作说明书步骤，将肽段进行不同同位素的标记，样本分离信息及相关标记名称见表1。⑤ 等量混合各种iTRAQ试剂标记的肽段。⑥ 混合样品的HPLC分级处理:选取的色谱柱为Agilent 300 Extend C18,长度250 mm,内径4.6 mm,粒径5 μm。在高pH值状态下，对已混合处理好的肽段进行反向高效液相色谱法(HPLC)分级。⑦ 相关质谱数据分析:使用EASY-nLC 1200超高效液相系统对处理好的肽段进行分离，所有已处理好的肽段均使用液相色谱流动相A进行溶解，并将其注入NSI离子源中集中进行电离操作，然后进入Q Exactive HF-X质谱进行分析。⑧ 生物信息学分析:使用Maxquant (v1.5.2.8)对二级质谱所获得的数据进行检索分析，相关参数设置的数据库为Mycobacterium_tuberculosis_Uniprot (6 529条序列)，数据中的假阳性率(FDR)可通过添加反库计算方式进行优化，在数据库中添加常见污染库，以消除样本鉴定结果中污染蛋白的影响,iTRAQ-4plex参数为1%,蛋白鉴定参数为1%,PSM鉴定的FDR参数为1%,相关参数设置见表2。

表1 样本分离信息及相关标记名称

表2 质谱数据分析相关参数设置

2 结 果

2.1 蛋白鉴定

iTRAQ定量质谱分析共得到365 454张二级谱图。蛋白理论数据进行搜库及对比后共得到可利用有效谱图数量为5 624个，利用率为1.5%。通过谱图解析鉴定出518条肽段，其中514条肽段为特异性肽段，共鉴定出238个蛋白，其中187个可定量。

2.2 蛋白差异表达分析

通过iTRAQ质谱定量分析检测出肽段对应在每个样本中的定量值，每个蛋白会对应多个肽段，计算蛋白在每个样本中的定量值时，采用蛋白所对应特异性肽段的定量值中值作为该蛋白的定量值。对于每一次重复实验,2个不同样本之间蛋白定量比值作为差异表达量(Ratio)。本项目分别进行了2次不同重复实验，可取2次重复Ratio值的平均值作为最终Ratio值。同时，可以取2次不同重复实验的Ratio值变异系数CV值作为最终CV值。当CV值<0.1时，差异表达量变化超过1.2作为显著上调的变化阈值，小于1/1.2作为显著下调的变化阈值。差异表达蛋白统计信息见表3。

表3 差异表达蛋白信息

2.3 基因的本体论(GO)分析

使用GO分析对差异表达的蛋白进行统计分析，通过数据库对比，找出标准株与休眠株中不同差异蛋白在不同信号通路中的功能，根据参与条款能量代谢途径的差异蛋白进行分析。差异表达蛋白在GO二级分类中统计分布见图1。

2.4 同源蛋白簇分析(COG/KOG)功能分类

将所获得的质谱数据导入到数据库进行比对分析，通过数据库比对，找出标准株与休眠株中不同差异蛋白在不同信号通路中的功能，将差异蛋白进行COG/KOG功能分类统计(见图2)，发现不同生理状态下的结核杆菌菌株间，差异性表达的蛋白质同样与能量产生和转换、蛋白质翻译等关系。

2.5 差异蛋白互作网络分析

对结核分枝杆菌标准株H37Rv和诱导的潜伏性感染菌株样本中所获得表达差异的蛋白进行统一编码，将数据导入到STRING(v.10.5)蛋白网络互作数据库进行差异蛋白分析，按照参数confidence score >0.7 (high confidence)提取数据，从而分析标准株和潜伏性感染菌株中差异表达蛋白的相互作用网络。差异蛋白互作网络可通过R package“networkD3”工具构建展示，差异表达蛋白使用圆圈表示，图中标记的红颜色表示上调蛋白，蓝颜色表示下调蛋白。圆圈大小表示2种样本中差异表达蛋白及所对应的相互作用蛋白个数(圆圈越大说明与其互作的蛋白越多，圆圈越小说明与其互作的蛋白越少)。图中编号为C的部分为能量代谢相关通路蛋白，可通过差异蛋白的功能富集分析，清晰找出哪些差异蛋白参与了代谢调控潜伏性感染菌株能量代谢相关通路调节，从而进一步确定关键标志物蛋白。取排名前50个且互作关系最紧密的蛋白，最终通过数据库对比分析绘制出差异表达蛋白互作网络图，见图3。

2.6 具有鉴别诊断价值的蛋白质及互作网络

蛋白质组学差异表达结果以及生物信息学分析后得到详细蛋白及相关通路信息，见表4。筛选出的P9WPC9(基因名称为ClpC1)是潜伏期肺结核具有潜在诊断价值的蛋白质。休眠菌株P9WPC9表达水平是标准结核杆菌菌株的2.2倍，其网络互作情况见图4。

3 讨 论

近年来，虽然医疗水平已有大幅提高，但结核病的诊断及治疗仍是严重的公共卫生问题[6-7]。耐药菌株的不断出现[8]以及糖尿病、免疫缺陷性疾病等合并结核感染更增加了结核病的防治难度。传统结核菌素试验无法明确区分卡介苗接种者、非结核分枝杆菌感染者、结核潜伏性感染者及结核病患者[9];IFN-γ释放试验亦无法鉴别结核潜伏性感染及活动性结核病[10],若患者伴有其他免疫缺陷性疾病[11-12],则可会出现假阴性结果。因此寻求更加准确、有效的结核潜伏性感染诊断标志物至关重要。

从蛋白质组学分析结果中可知，所测差异蛋白主要与细胞能量代谢相关。进一步将差异蛋白进行COG/KOG功能分类统计，发现不同生理状态下结核杆菌菌株间，休眠株中P9WPC9表达上调2.2倍,P9WPC9为ATP依赖性Clp蛋白酶ATP结合亚基(基因名称为ClpC1),可以将蛋白酶引导至特定的底物处，亦可在没有ClpP的情况下执行伴侣功能，在ClpP2存在的情况下降解抗σ-E因子RseA，从而调节能量代谢途径，其调控主要与能量产生和转换、蛋白质翻译等相关。

因此,P9WPC9有望成为结核潜伏性感染的潜在诊断性标志物。下一步可应用western-blot对P9WPC9蛋白质水平进行验证，应用实时荧光定量PCR(qRT-PCR)等技术观察ClpC1 mRNA水平的改变。研究潜伏性感染及活动性结核病患者P9WPC9表达的改变，从而进一步探讨P9WPC9在结核潜伏性感染诊断中的临床应用价值。