一种急性胰腺炎预后标志物、急性胰腺炎预后预测模型及其
应用
技术领域
本发明涉及生物技术领域,涉及基因诊断技术领域,涉及一种疾病预后标志物、疾病预后预测模型及其应用,具体涉及一种急性胰腺炎预后标志物、急性胰腺炎预后预测模型及其在预测持续性器官衰竭状态中的应用。
背景技术
急性胰腺炎(AP)是常见的胃肠道疾病,在全世界范围内发病率呈逐年递增的趋势,给当前的医疗保健系统带来了沉重的负担。大多数AP患者的临床病程较轻,可以治愈且没有后遗症。尽管如此,10%~20%的AP患者发生有持续性器官衰竭(POF),被认为是不良预后的主要决定因素,死亡率高达30%。早期预测AP患者在预后过程中是否发生POF具有非常重要的临床意义,有利于医生在早期即可预测急性胰腺炎的转归情况,显著改善AP患者的预后。
目前,部分临床实践指南中建议使用评分系统预测AP发病的严重程度,例如Ranson评分系统、急性生理和慢性健康检测(APACHE)系统II、以及AP严重程度的床旁指数(BISAP)。但是,这些评分系统存在许多不足,包括准确性差、效率较低、耗时较长、操作复杂等。尽管人们也在关注血尿素氮(BUN)和肌酐作为参数在预测急性胰腺炎严重程度和预后中的作用,但是目前仍然缺乏有效的早期预测分子标志物。
微小RNA(hsa-miRNA)是一类内源性非编码RNA,可以在转录后调控基因的表达。研究发现,hsa-miR-216a在AP患者的血浆中显著升高,并且可以作为标志物诊断具有高度特异性的AP。但是,目前尚未有研究探索循环hsa-miRNA在预测AP患者是否发生POF方面的潜能。外泌体是直径为30~150nm的细胞外囊泡,体液中的外泌体是疾病生物标志物的稳定来源。
CN101827941A公开了用于区分胰腺癌与正常胰腺功能和/或慢性胰腺炎的方法,所述方法包括检测hsa-miR基因产物在组织样品中的表达水平,所述hsa-miR基因包括hsa-miR-21、hsa-miR-221、hsa-miR-222、hsa-miR-181a、hsa-miR-181b、hsa-miR-181d和hsa-miR-155及其组合,但是上述hsa-miR基因不是AP预后标志物,不能用于预测AP患者在预后过程中是否发生POF。
CN109777874A公开了一种适用于胰腺导管腺癌诊断及预后判断的血浆外泌体hsa-miRNA标志物及应用,对胰腺导管腺癌具有较高诊断及鉴别诊断价值的生物标志物包括血浆外泌体hsa-miRNA-95-3p、hsa-miRNA-26b-5p、hsa-miRNA-335-5p;通过血浆外泌体hsa-miRNA标志物和诊断试剂盒的研制和应用,为胰腺导管腺癌做出更精确的诊断,尤其是辅助临床医生将胰腺导管腺癌与慢性胰腺炎进行鉴别诊断,同时对患者的预后进行评估,为临床医生准确掌握患者病情,提高临床治疗效果奠定基础。但是上述hsa-miRNA同样不是AP预后标志物,不能用于预测AP患者在预后过程中是否发生POF。
因此,目前的发生POF的AP疾病尚无典型的诊断标志物,提供针对发生POF的AP疾病的特异性标志物,用于临床AP患者的预后评价和早期筛查,有利于为患者提供更加精确的预诊和诊断信息,具有重要意义。
发明内容
针对现有技术的不足和实际需求,本发明提供了一种急性胰腺炎预后标志物、急性胰腺炎预后预测模型及其应用,通过分析发生或未发生持续性器官衰竭的急性胰腺炎患者中血清外泌体hsa-miRNA的表达水平,获得的外泌体hsa-miRNA标志物与急性胰腺炎预后显著相关,由此构建的胰腺炎预后预测模型可以在早期预测急性胰腺炎患者是否发生持续性器官衰竭,有利于指导医生制定适宜的救助策略,改善急性胰腺炎患者的预后。
为达此目的,本发明采用以下技术方案:
第一方面,本发明提供了一种急性胰腺炎预后标志物,所述标志物包括hsa-miR-4265、hsa-miR-1208、hsa-miR-516a或hsa-miR-3127中的任意一种或至少两种的组合。本发明中,hsa-miR-4265、hsa-miR-1208、hsa-miR-516a和hsa-miR-3127与发生持续性器官衰竭(POF)的急性胰腺炎具有显著的相关性,根据上述标志物的相对表达量有助于构建急性胰腺炎预后预测模型,个体化预测急性胰腺炎患者是否发生持续性器官衰竭,为急性胰腺炎患者提供准确、有用的预后信息。
本发明中,hsa-miR-4265相较于hsa-miR-1208、hsa-miR-516a或hsa-miR-3127对急性胰腺炎发生持续性器官衰竭的预测能力更强,与急性胰腺炎不良预后显著相关,预测能力较现有的多因素评分系统APACHE II、BISAP和生化检验指标血尿素氮(BUN)、肌酐(Cr)得到显著提高。
优选地,所述标志物包括外泌体hsa-miR-3127、hsa-miR-4265、hsa-miR-1208或hsa-miR-516a中的任意一种或至少两种的组合。
根据本发明,外泌体中富含hsa-miRNA,外泌体hsa-miRNA与非外泌体hsa-miRNA相比,灵敏度更高,在临床诊断中具有更高的预测价值。
优选地,所述标志物包括血清外泌体hsa-miR-3127、hsa-miR-4265、hsa-miR-1208或hsa-miR-516a中的任意一种或至少两种的组合。
第二方面,本发明提供了一种急性胰腺炎预后预测模型,所述预测模型采用第一方面所述的标志物进行预后预测。
优选地,所述预测模型采用第一方面所述的标志物进行急性胰腺炎持续性器官衰竭状态的预测。
优选地,所述预测模型的输入变量为第一方面所述的标志物的相对表达量。
本发明中,第一方面所述的标志物的相对表达量采用定量PCR和/或高通量测序方法进行检测。
优选地,所述预测模型的输入变量为hsa-miR-4265、hsa-miR-1208和hsa-miR-3127的相对表达量。
优选地,所述预测模型的方程式为:
Cmi=6.043+0.244×hsa-miR-3127+0.511×hsa-miR-4265+0.309×hsa-miR-1208。
本发明中,以hsa-miR-4265、hsa-miR-1208和hsa-miR-3127的相对表达量作为输入变量构建的预测模型Cmi,与发生持续性器官衰竭的急性胰腺炎具有显著的相关性,对急性胰腺炎中持续性器官衰竭状态具有预测能力,与急性胰腺炎不良预后显著相关,预测能力较现有的多因素评分系统APACHE II、BISAP和生化检验指标血尿素氮(BUN)、肌酐(Cr)得到显著提高。
第三方面,本发明提供了一种第二方面所述的预测模型的构建方法,所述构建方法包括以下步骤:
(1)在已知持续性器官衰竭状态的急性胰腺炎样本中,进行hsa-miRNA表达谱分析,得到与持续性器官衰竭状态显著相关的hsa-miRNA;
(2)对所述与持续性器官衰竭状态显著相关的hsa-miRNA进行二项逻辑回归建模,构建得到所述预测模型。
优选地,步骤(1)所述样本包括血液、尿液、腹水、唾液或组织中的任意一种或至少两种的组合,优选为血液,进一步优选为血清外泌体。
优选地,步骤(1)所述与持续性器官衰竭状态显著相关的hsa-miRNA包括hsa-miR-4265、hsa-miR-1208、hsa-miR-516a、hsa-miR-3127、hsa-miR-631、hsa-miR-3180、hsa-miR-5194、hsa-miR-4294、hsa-miR-4513、hsa-miR-4685、hsa-miR-4776和hsa-miR-6083。
第四方面,本发明提供了一种预测急性胰腺炎预后的系统,所述系统采用第一方面所述的标志物进行预后预测。
优选地,所述系统采用第一方面所述的标志物进行急性胰腺炎持续性器官衰竭状态的预测。
优选地,所述系统的输入变量为第一方面所述的标志物的相对表达量。
优选地,所述系统的输入变量为hsa-miR-4265、hsa-miR-1208和hsa-miR-3127的相对表达量。
优选地,所述系统的方程式为:
Cmi=6.043+0.244×hsa-miR-3127+0.511×hsa-miR-4265+0.309×hsa-miR-1208。
第五方面,本发明提供了一种预测急性胰腺炎预后的装置,所述装置包括检测单元和分析单元。
优选地,所述检测单元检测样本中第一方面所述的标志物的相对表达量。
优选地,所述检测方法为定量PCR。
优选地,所述分析单元将第一方面所述的标志物的相对表达量作为输入变量,输入第二方面所述的预测模型和/或第四方面所述的系统,分析持续性器官衰竭状态。
第六方面,本方面提供了一种第一方面所述的标志物、第二方面所述的预测模型、第四方面所述的系统或第五方面所述的装置在制备急性胰腺炎预后诊断试剂和/或急性胰腺炎预后治疗药物中的应用。
与现有技术相比,本发明具有如下有益效果:
(1)本发明的标志物hsa-miR-3127、hsa-miR-1208、hsa-miR-516a和hsa-miR-4265与急性胰腺炎预后具有显著的相关性,是判断急性胰腺炎是否发生持续性器官衰竭最为灵敏和特异的、先于临床病理学特征的标志物;
(2)本发明以hsa-miR-4265、hsa-miR-1208和hsa-miR-3127的相对表达量作为输入变量构建的预测模型Cmi,与发生持续性器官衰竭的急性胰腺炎具有显著的相关性,可以准确预测胰腺炎预后是否发生持续性器官衰竭(POF),有助于个体化预测胰腺炎患者的预后情况,对临床医生进行急性胰腺炎患者管理具有指导作用。
附图说明
图1(A)为外泌体的电子显微镜图,图1(B)为外泌体的纳米粒子跟踪分析图;
图2(A)为hsa-miR-4265在训练集中的箱型图,图2(B)为hsa-miR-1208在训练集中的箱型图,图2(C)为hsa-miR-516a在训练集中的箱型图,图2(D)为Cmi模型在训练集中的箱型图;
图3(A)为hsa-miR-4265在第一验证集中的箱型图,图3(B)为hsa-miR-1208在第一验证集中的箱型图,图3(C)为hsa-miR-516a在第一验证集中的箱型图,图3(D)为Cmi模型在第一验证集中的箱型图;
图4(A)为hsa-miR-4265在第二验证集中的箱型图,图4(B)为hsa-miR-1208在第二验证集中的箱型图,图4(C)为hsa-miR-516a在第二验证集中的箱型图,图4(D)为Cmi模型在第二验证集中的箱型图;
图5(A)为hsa-miR-4265在合并的V1+V2中的箱型图,图5(B)为hsa-miR-1208在合并的V1+V2中的箱型图,图5(C)为hsa-miR-516a在合并的V1+V2中的箱型图,图5(D)为Cmi模型在合并的V1+V2中的箱型图;
图6(A)为hsa-miRNA-4265与AP患者的不良预后显著相关,图6(B)为Cmi模型与AP患者的不良预后显著相关;
图7(A)为训练集中hsa-miRNA-4265和Cmi模型与现有方法的POF预测能力比较,图7(B)为第一验证集中hsa-miRNA-4265和Cmi模型与现有方法的POF预测能力比较,图7(C)为第二验证集中hsa-miRNA-4265和Cmi模型与现有方法的POF预测能力比较,图7(D)为合并的V1+V2中hsa-miRNA-4265和Cmi模型与现有方法的POF预测能力比较。
具体实施方式
为进一步阐述本发明所采取的技术手段及其效果,以下结合实施例和附图对本发明作进一步地说明。可以理解的是,此处所描述的具体实施方式仅仅用于解释本发明,而非对本发明的限定。
实施例中未注明具体技术或条件者,按照本领域内的文献所描述的技术或条件,或者按照产品说明书进行。所用试剂或仪器未注明生产厂商者,均为可通过正规渠道商购获得的常规产品。
实施例1外泌体的提取和鉴定
本实施例从急性胰腺炎患者的血清样本中分离、筛选得到血清外泌体,电子显微镜结果如图1(A)所示,外泌体为存在双层膜结构的圆形颗粒,纳米粒子跟踪分析实验结果如图1(B)所示,大多数粒子的直径在30~150nm范围内,符合外泌体鉴定标准,说明从血清样本中成功提取得到血清外泌体。
实施例2样本的hsa-miRNA表达谱
本发明的研究样本包括790例急性胰腺炎(AP)患者的血清样本,分为训练集和两个验证集,其中,发生持续性器官衰竭(POF)的AP患者与非POF的AP患者相比,预后差,采用APACHE或BISAP系统进行评分后获得的分数高,血尿素氮(BUN)和肌酐参数高。
在发现阶段,采用微阵列对5例发生POF的AP样本和5例非POF的AP样本进行分析,鉴定得到90种差异表达的外泌体hsa-miRNA;对这90种外泌体hsa-miRNA表达谱进行聚类分析,可以实现上述10例样本的绝对分离,说明血清外泌体hsa-miRNA对于发生POF的AP具有潜在的预测能力。为了最大程度减少假阳性结果,从90种外泌体hsa-miRNA中排除了其中的55种。
实施例3预测模型的建立
在包含41例POF的AP样本和186例非POF的AP样本的训练集中,采用定量PCR对剩余的35种外泌体hsa-miRNA进行验证,根据检出率小于75%,排除了其中的17种hsa-miRNA;
在剩余的18种hsa-miRNA中,有12种hsa-miRNA(hsa-miR-4265、1208、516a、3127、631、3180、5194、4294、4513、4685、4776和6083)与POF状态显著相关(p<0.001);对这12种hsa-miRNA进行二项逻辑回归建模,构建得到基于3种hsa-miRNA(hsa-miR-4265、1208、3127)的预测AP患者是否发生POF的数学模型,并确定了模型的方程式:
Cmi=6.043+0.244×hsa-miR-3127+0.511×hsa-miR-4265+0.309×hsa-miR-1208
本实施例同时选择了hsa-miR-4265、1208和516a作为单独的候选标志物进行POF预测分析。
hsa-miR-4265、1208、516a和Cmi模型在训练集中的箱型图如图2(A)、图2(B)、图2(C)和图2(D)所示,可以看出,hsa-miR-4265、1208、516a和Cmi与POF具有显著的相关性,AUROC分别为:hsa-miR-4265:0.89(95%CI,0.84-0.95,P<0.0001),hsa-miR-1208:0.86(95%CI,0.79-0.93,P<0.0001),hsa-miR-516a:0.85(95%CI,0.78-0.91,P<0.0001),Cmi:0.90(95%CI,0.85-0.96,P<0.0001)。
实施例4验证hsa-miRNA-4265和Cmi模型对POF状态的预测价值
在两个独立的验证集中,hsa-miR-4265、3127、1208和516a在发生POF的AP样本中持续升高。
在第一验证集V1中,如图3(A)、图3(B)、图3(C)和图3(D)所示,hsa-miR-4265的AUC为0.91(95%CI,0.86-0.95),hsa-miR-1208的AUC为0.85(95%CI,0.78-0.91),hsa-miR-516a的AUC为0.78(95%CI,0.70-0.86),Cmi模型的AUC为0.89(95%CI,0.84-0.94);
在第二验证集V2中,如图4(A)、图4(B)、图4(C)和图4(D)所示,hsa-miR-4265的AUC为0.89(95%CI,0.83-0.94),hsa-miR-1208的AUC为0.80(95%CI,0.71-0.89),hsa-miR-516a的AUC为0.76(95%CI,0.67-0.85),Cmi模型的AUC为0.86(95%CI,0.80-0.92);
在合并的V1+V2中,如图5(A)、图5(B)、图5(C)和图5(D)所示,hsa-miR-4265的AUC为0.89(95%CI,0.86-0.93),hsa-miR-1208的AUC为0.83(95%CI,0.77-0.88),hsa-miR-516a的AUC为0.77(95%CI,0.71-0.83),Cmi模型的AUC为0.88(95%CI,0.84-0.91);
可以看出,在验证集中,hsa-miR-4265、1208、516a和Cmi与POF同样具有显著的相关性,与hsa-miR-1208和516a相比,hsa-miR-4265和Cmi的预测能力更强,可能作为预测POF的标志物。
实施例5hsa-miRNA-4265和Cmi模型与临床结果的相关性
鉴于POF是AP严重程度的关键决定因素,与AP不良预后密切相关,本实施例检测了hsa-miRNA-4265和Cmi模型与临床结果的相关性。
如图6(A)和图6(B)所示,hsa-miRNA-4265和Cmi模型与AP患者的不良预后显著相关。
实施例6比较hsa-miRNA-4265和Cmi模型与现有方法的POF预测能力
本实施例比较了hsa-miRNA-4265和Cmi模型与多因素评分系统APACHE II、BISAP和生化检验指标血尿素氮(BUN)、肌酐(Cr)的POF预测能力。
如图7(A)所示,在训练集(Training Cohort)中,与现有的评分系统APACHE II、BISAP和评价参数血尿素氮(BUN)、肌酐(Cr)相比,外泌体hsa-miRNA-4265和Cmi模型具有最优的POF预测能力;
如图7(B)、图7(C)和图7(D)所示,在第一验证集(Validation cohort1)、第二验证集(Validation cohort 2)和合并的V1+V2(Validation cohort 1+2)中,外泌体hsa-miRNA-4265比APACHE II、BISAP、BUN和Cr具有更优的POF预测能力,Cmi模型比BUN和Cr具有更优的POF预测能力,但在第一验证集和第二验证集中,Cmi模型的预测能力稍弱于APACHEII和BISAP。
综上所述,本发明通过血清外泌体hsa-miRNA检测和数学建模,为急性胰腺炎患者不良预后的预测提供了高灵敏度和高特异性的分子标志物组合,有利于临床医生更加精准地进行急性胰腺炎患者管理,并显著改善急性胰腺炎患者的预后。
申请人声明,本发明通过上述实施例来说明本发明的详细方法,但本发明并不局限于上述详细方法,即不意味着本发明必须依赖上述详细方法才能实施。所属技术领域的技术人员应该明了,对本发明的任何改进,对本发明产品各原料的等效替换及辅助成分的添加、具体方式的选择等,均落在本发明的保护范围和公开范围之内。