CN113345592B - 一种急性髓细胞样白血病预后风险模型的构建及诊断设备 - Google Patents

Abstract

本发明涉及一种急性髓细胞样白血病预后风险模型的构建及诊断设备。模型的构建方法包括：获取缺氧相关基因；获取急性髓细胞样白血病基因表达数据作为训练数据集；利用回归模型或机器学习，在训练数据集中，对所述缺氧相关基因进行特征选择，得到预后风险相关基因；利用所述预后风险相关基因构建急性髓细胞样白血病预后风险模型。本发明构建的模型预测急性髓细胞样白血病患者的预后准确率高，为临床治疗急性髓细胞样白血病指明了方向。

Description

一种急性髓细胞样白血病预后风险模型的构建及诊断设备

技术领域

本发明涉及生物信息技术领域，更具体地，涉及一种急性髓细胞样白血病预后风险模型的构建、诊断设备和计算机可读存储介质。

背景技术

急性髓系白血病(AML)是一种侵袭性的髓系恶性肿瘤，是成人中最常见的急性白血病类型，约占病例的70％，其特点是细胞快速增殖和高死亡率。由于对急性髓系白血病缺乏完整的认识，其临床结果仍不令人满意。经典治疗策略从化疗到靶向治疗效果都差强人意，对于大多数超过60岁的病人术后五年生存率只有10-20％，复发的患者很少能存活超过5年。为了提高AML患者的生存率，需要开发新的AML治疗方法和预测其预后的方法。

肿瘤缺氧是恶性肿瘤实体生长的标志，深刻影响恶性肿瘤的进展，并有助于发展治疗耐药性。然而，到目前为止，还没有关于AML的缺氧相关研究。

发明内容

本研究首次开发了AML缺氧预后风险模型，该模型可用于预测AML组的预后，模型所涉及的相关分子指标会帮助临床医生做出重要的治疗决策。

本发明提供一种急性髓细胞样白血病预后风险模型的构建方法，包括：

获取缺氧相关基因；

获取急性髓细胞样白血病基因表达数据作为训练数据集；

利用回归模型或机器学习，在训练数据集中，对所述缺氧相关基因进行特征选择，得到预后风险相关基因；

利用所述预后风险相关基因构建急性髓细胞样白血病预后风险模型。

进一步，所述获取缺氧相关基因是获取GSEA数据库中缺氧相关基因。

进一步，获取GSEA数据库(hallmark-hypoxia)200个缺氧相关基因。

当组织得不到充足的氧，或不能充分利用氧时，组织的代谢、机能、甚至形态结构都可能发生异常变化，这一病理过程称为缺氧(hypoxia)。细胞对缺氧的反应以一定数量的基因表达改变为特征，这些表达改变的基因被称之为缺氧相关基因。

进一步，获取缺氧相关基因还包括如下处理：通过蛋白质相互网络、基因共表达网络、转录因子调控网络、pathway网络中的一种或几种基于网络分析方法筛选出排名靠前的缺氧相关基因。

蛋白质-蛋白质相互作用(protein-protein interaction，PPI)是指两个或两个以上的蛋白质分子通过非共价键形成蛋白质复合体(protein complex)的过程。

可选的，通过PPI网络分析方法筛选出排名靠前的缺氧相关基因；更可选的，通过MCC、DMNC、MNC、Degree、EPC、EcCentricity、Closeness、Radiality中的一种或几种分析方法筛选出排名靠前的缺氧相关基因。

所述排名靠前可以是排名前50、40、30、20、10、5、3等任选整数，优选排名前5、前10或排名前20。

进一步，获取TCGA或GEO等公开数据库中的急性髓细胞样白血病基因表达数据作为训练数据集；可选的，将所述TCGA或GEO等公开数据库中的急性髓细胞样白血病基因表达数据分为训练集和验证集。

进一步，所述利用回归模型是通过单因素COX回归模型进行特征选择，得到预后风险相关基因，利用所述预后风险相关基因，采用LASSO COX回归分析或多因素COX回归分析，构建急性髓细胞样白血病预后风险模型。

Cox回归模型，又称“比例风险回归模型(proportional hazards model，简称Cox模型)”，是由英国统计学家提出的一种半参数回归模型。该模型以生存结局和生存时间为应变量，可同时分析众多因素对生存期的影响，能分析带有截尾生存时间的资料，且不要求估计资料的生存分布类型。

基因表达是指将来自基因的遗传信息合成功能性基因产物的过程，基因表达数据包括mRNA表达水平数据和/或蛋白表达水平数据。

进一步，将所述构建急性髓细胞样白血病预后风险模型作为特征向量，选取临床病理特征，利用所述特征向量和临床病理特征构建急性髓细胞样白血病预后综合风险模型。

可选的，临床病理特征由临床医生挑选出与急性髓细胞样白血病可能相关的候选临床病理因素；更可选的，临床病理特征包括年龄、性别等；更可选的，临床病理特征为年龄。

进一步，所述选取临床病理特征包括：采用回归模型或机器学习筛选出与急性髓细胞样白血病生存相关的临床病理特征。可选的，将具体某一项临床病理特征进行单因素Cox独立预后分析，筛选出与急性髓细胞样白血病生存相关的临床病理特征。

可选的，所述与急性髓细胞样白血病生存相关的临床病理特征是指进行单因素Cox独立预后分析时，pvalue值小于0.05。更可选的，pvalue值小于0.01。

进一步，输入年龄等临床病理因素与所述构建的急性髓细胞样白血病预后风险模型，采用多因素回归模型，构建急性髓细胞样白血病预后综合风险模型。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器；

所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的基因表达数据或基因表达数据及临床病理特征；

将所述待测样本的基因表达数据输入所述的急性髓细胞样白血病预后风险模型中，或将所述待测样本的基因表达数据及临床病理特征输入所述的急性髓细胞样白血病预后综合风险模型中；

获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：获取待测样本的EGFR、PYGM和GPC3基因表达数据；

根据EGFR、PYGM和GPC3基因表达情况，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：获取待测样本的EGFR、PYGM和GPC3基因表达数据和临床病理特征；

根据EGFR、PYGM和GPC3基因表达情况和临床病理特征，获得待测样本急性髓细胞样白血病预后风险预测结果。

可选的，所述临床病理特征为年龄。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的基因表达数据或基因表达数据及临床病理特征；

处理单元，用于将所述待测样本的基因表达数据输入所述的急性髓细胞样白血病预后风险模型中，或将所述待测样本的基因表达数据及临床病理特征输入所述的急性髓细胞样白血病预后综合风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据；

处理单元，用于根据EGFR、PYGM和GPC3基因表达情况，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据及临床病理特征；

处理单元，用于根据EGFR、PYGM和GPC3基因表达数据及临床病理特征，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的急性髓细胞样白血病预后风险预测。

本申请的优点：本申请利用数据库中的AML数据，基于缺氧相关基因构建了一个高效急性髓细胞样白血病预后风险模型，继而将该急性髓细胞样白血病预后风险模型与临床病理特征结合，获得更加准确的急性髓细胞样白血病预后综合风险模型，用临床病理特征和多个基因特征来预测急性髓细胞样白血病患者的预后，准确率高，本申请为临床治疗急性髓细胞样白血病指明了方向。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获取其他的附图。

图1是本发明实施例提供的一种基于急性髓细胞样白血病预后风险模型的构建示意流程图；

图2是训练集缺氧相关基因单因素Cox回归分析结果；

图3是训练集缺氧相关基因单因素Cox分析的森林图；

图4是训练集缺氧相关基因多因素Cox分析结果；

图5是训练集缺氧相关基因多因素Cox分析的森林图；

图6是TCGA-LAML训练集中高低风险分组的风险曲线，A图显示患者从左到右风险得分依次增大；B图横坐标为根据风险得分高低进行排序的患者样本，纵坐标为生存时间，C图为高低风险组的基因表达热图；

图7是TCGA-LAML训练集中高低风险分组的生存曲线；

图8是TCGA-LAML训练集中1/3/5年的ROC曲线；

图9是TCGA-LAML验证集中高低风险组生存曲线；

图10是TCGA-LAML验证集中1/3/5年的ROC曲线；

图11是本发明实施例提供的一种基于急性髓细胞样白血病预后综合风险模型的构建示意流程图；

图12是TCGA-SKCM不同临床信息下的生存和表达分析，A为年龄组，B为性别组；

图13是风险基因的表达与临床特征的相关性；

图14是单因素Cox独立预后分析；

图15是多因素Cox独立预后分析；

图16是列线图预测患者1、3和5年生存率；

图17是列线图1、3和5年矫正曲线；

图18是本发明实施例提供的一种急性髓细胞样白血病预后风险诊断设备；

具体实施方式

为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。

在本发明的说明书和权利要求书及上述附图中的描述的一些流程中，包含了按照特定顺序出现的多个操作，但是应该清楚了解，这些操作可以不按照其在本文中出现的顺序来执行或并行执行，操作的序号如101、102等，仅仅是用于区分开各个不同的操作，序号本身不代表任何的执行顺序。另外，这些流程可以包括更多或更少的操作，并且这些操作可以按顺序执行或并行执行。需要说明的是，本文中的“第一”、“第二”等描述，是用于区分不同的消息、设备、模块等，不代表先后顺序，也不限定“第一”和“第二”是不同的类型。

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获取的所有其他实施例，都属于本发明保护的范围。

图1是本发明实施例提供的一种基于急性髓细胞样白血病预后风险模型的构建示意流程图，具体地，所述方法包括如下步骤：

101：获取缺氧相关基因；

102：获取急性髓细胞样白血病基因表达数据作为训练数据集；

103：利用回归模型或机器学习，在训练数据集中，对所述缺氧相关基因进行特征选择，得到预后风险相关基因；

104：利用所述预后风险相关基因构建急性髓细胞样白血病预后风险模型。

在一个实施例中，101和102步骤可进行调换，先获取急性髓细胞样白血病基因表达数据作为训练数据集，后获取缺氧相关基因，或者同时获取。

在一个实施例中，获取GSEA数据库中缺氧相关基因，从“基因集合富集分析”(hallmark-hypoxia)中下载了与缺氧相关的基因集，其中包含200个对缺氧水平反应上调的基因。

在一个实施例中，获取缺氧相关基因还包括如下处理：通过蛋白质相互网络、基因共表达网络、转录因子调控网络、pathway网络中的一种或几种基于网络分析方法筛选出排名靠前的缺氧相关基因。

在一个实施例中，通过PPI网络分析方法筛选出排名靠前的缺氧相关基因；在一个实施例中，通过MCC、DMNC、MNC、Degree、EPC、EcCentricity、Closeness、Radiality中的一种或几种分析方法筛选出排名靠前的缺氧相关基因。

在一个实施例中，从TCGA数据库中获取140例AML患者的RNA-seq转录组数据和相应的临床病理信息。数据按7:3的比例分为训练集和验证集。

在一个实施例中，对获取140例AML患者的RNA-seq转录组数据进行归一化处理。

在一个实施例中，利用回归模型中的单因素Cox和/或多因素Cox回归分析，在训练数据集中，对所述GSEA数据库中缺氧相关基因进行特征选择，特征选择后的缺氧相关基因为危险因素。

在一个实施例中，利用单因素Cox回归分析发现有统计学意义的缺氧相关基因，得到候选预后风险相关基因，然后利用多因素Cox回归分析候选预后风险相关基因，得到预后风险相关基因及其回归系数，构建急性髓细胞样白血病预后风险模型：

其中，N＝缺氧基因的数量，Expi为缺氧基因的各自表达量，Coei为相应的多因素Cox回归系数。

可选的，N＝3-5，Expi为缺氧基因的各自表达量，Coei为相应的多因素Cox回归系数。

在一个实施例中，利用单因素Cox回归分析，在TCGA训练数据集中进行了回归分析，获得6个缺氧相关基因(ERRFI1、EGFR、SDC2、PPP1R3C、PYGM和GPC3)与患者的5年总生存期(OS)显著相关(见图2和图3)。然后对6个缺氧相关基因进行多因素Cox回归分析，选择P<0.05的3个缺氧相关基因(EGFR、PYGM和GPC3)构建急性髓细胞样白血病预后风险模型(图4和图5)，公式如下：

风险评分(riskscore)＝(-6.87*EGFR基因表达数据)+(0.38*PYGM基因表达数据)+(0.24*GPC3基因表达数据)

在一个实施例中，基因表达数据是向量，具有大小和正负。

在一个实施例中，根据构建的急性髓细胞样白血病预后风险模型，对急性髓细胞样白血病患者进行打分，得到每个患者的风险值，以风险评分中位值为界，将急性髓细胞样白血病患者分为高风险组和低风险组。根据风险模型绘制风险曲线图(图6)，风险曲线图由A、B、C三部分图组成，A图和B图中横坐标一致，为根据风险得分高低进行排序的患者样本，从左到右风险得分依次增大；纵坐标分别为风险得分和生存时间。虚线为风险得分中位值及其对应的病人数量。C图为高低风险组的基因表达热图。根据高低风险分组绘制OS(Overall survival)生存曲线，结果见图7。从图中可以看到高低风险组间病人的生存情况有显著差异(p<0.05)，高风险组的病人生存率更低。为进一步评估急性髓细胞样白血病预后风险模型的有效性，采用ROC计算模型的AUC面积以评估急性髓细胞样白血病预后风险模型的有效性。ROC曲线是一种显示分类模型在所有分类阈值下的效果的图表。该曲线绘制了以下两个参数：真阳性(True positive rate)和假阳性(False positive rate)。曲线下方部分的面积被称为AUC(Area Under Curve)，用来表示预测准确性，AUC值越高，也就是曲线下方面积越大，表示预测准确率越高。根据多因素风险回归得到的风险模型，以1、3、5年为生存时间节点，绘制ROC曲线，结果见图8。可以看到ROC曲线的AUC均大于0.6，急性髓细胞样白血病预后风险模型的效能较好。

在一个实施例中，我们使用验证集对急性髓细胞样白血病预后风险模型进行验证，我们同样对验证数据绘制风险曲线、高低风险组生存曲线(图9)和ROC曲线(图10)，可以看到验证集数据中同样高风险组在生存率低，且ROC曲线中1、3、5年节点的AUC均大于0.6，与训练集结果一致。

图11是本发明实施例提供的一种基于急性髓细胞样白血病预后综合风险模型的构建示意流程图，具体地，所述方法包括如下步骤：

201：获取缺氧相关基因；

202：获取急性髓细胞样白血病基因表达数据作为训练数据集；

203：利用回归模型或机器学习，在训练数据集中，对所述缺氧相关基因进行特征选择，得到预后风险相关基因；

204：利用所述预后风险相关基因构建急性髓细胞样白血病预后风险模型，将所述构建急性髓细胞样白血病预后风险模型作为特征向量；

205：选取临床病理特征，利用所述特征向量和临床病理特征构建急性髓细胞样白血病预后综合风险模型。

在一个实施例中，所述选取临床病理特征包括：采用回归模型或机器学习筛选出与急性髓细胞样白血病生存相关的临床病理特征。

在一个实施例中，机器学习的方法包括递归特征消除、CART、随机森林、线性回归、朴素贝叶斯、自定义的训练模型中的一种或几种。

在一个实施例中，将age(年龄)、gender(性别)临床病理因素与急性髓细胞样白血病预后风险模型进行生存(图12)和表达相关性分析(图13)，比较不同临床信息下的生存和表达是否有显著差异，我们发现风险评分在年龄、性别的不同分组之间无显著差异，但是年龄不同分组之间生存存在显著差异，而性别不同分组之间的生存无显著差异。

在一个实施例中，将age(年龄)、gender(性别)这2个临床病理因素进行单因素Cox独立预后分析(图14)，结果显示age和riskScore的pvalue均小0.05。因此选取age和riskScore纳作为特征向量进行多因素Cox独立预后分析(图15)。结果显示，age和riskScore的pvalue均小于0.05。说明构建的急性髓细胞样白血病预后综合风险模型有很好的预后价值。

在一个实施例中，选取年龄(年龄大于55或年龄小于等于55)作为临床病理特征，利用所述构建急性髓细胞样白血病预后风险模型(riskscore＝(-6.87*EGFR基因表达量)+(0.38*PYGM基因表达量)+(0.24*GPC3基因表达量))作为特征向量，构建急性髓细胞样白血病预后综合风险模型。

在一个实施例中，进行急性髓细胞样白血病患者生存列线图的构建及验证，将多因素Cox独立预后分析中pvalue小于0.05因素用于构建列线图，列线图对年龄和riskScore因素分别进行打分，每一种因素对应一个评分，各因素总评分相加对应总评分，然后根据总评分预测1、3和5年生存率(图16)，分数越高，生存率越低。基于上述预后综合风险模型绘制列线图绘制校准曲线(图17)，斜率越接近1，说明预测越准确。结果得知，本申请构建的预后综合风险模型对急性髓细胞样白血病患者1年、3年和5年生存率的预测准确度较高，说明构建的预后综合风险模型为有效模型。

在一个实施例中，我们以riskScore中位数将急性髓细胞样白血病样本分为高、低风险两组进行GSEA富集分析。GSEA对HALLMARK定义的基因集h.all.v7.2.symbols.gmt富集分析，我们发现低风险组富集到了HALLMARK HYPOXIA、HALLMARK ANGIOGENESIS和HALLMARKPI3K AKT MTOR SIGNALING等缺氧促进分化相关的通路；GSEA对KEGG定义的基因集c2.cp.kegg.v7.2.symbols.gmt富集分析，发现高风险组富集到了KEGG RIG ILIKERECEPTOR SIGNALING PATHWAY、KEGG PRIMARY IMMUNODEFICIENCY、KEGG ACUTE MYELOIDLEUKEMIA和KEGG BASAL CELL CARCINOMA等疾病免疫缺陷相关的通路。

图18是本发明实施例提供的一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器；

所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的基因表达数据；

将所述待测样本的基因表达数据输入所述的急性髓细胞样白血病预后风险模型中；

获得待测样本急性髓细胞样白血病预后风险预测结果。

在一个实施例中，本发明提供的一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器；

所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的基因表达数据及临床病理特征；

将所述待测样本的基因表达数据及临床病理特征输入所述的急性髓细胞样白血病预后综合风险模型中；

获得待测样本急性髓细胞样白血病预后风险预测结果。

在一个实施例中，本发明提供一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：获取待测样本的EGFR、PYGM和GPC3基因表达数据；

根据EGFR、PYGM和GPC3基因表达情况，获得待测样本急性髓细胞样白血病预后风险预测结果。

在一个实施例中，本发明提供一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：获取待测样本的EGFR、PYGM和GPC3基因表达数据和临床病理特征；

根据EGFR、PYGM和GPC3基因表达情况和临床病理特征，获得待测样本急性髓细胞样白血病预后风险预测结果。

可选的，所述临床病理特征为年龄。

本发明实施例提供的一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的基因表达数据；

处理单元，用于将所述待测样本的基因表达数据输入所述的急性髓细胞样白血病预后风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

在一个实施例中，本发明提供的一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的基因表达数据及临床病理特征；

处理单元，用于将所述待测样本的基因表达数据及临床病理特征输入所述的急性髓细胞样白血病预后综合风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据；

处理单元，用于根据EGFR、PYGM和GPC3基因表达情况，获得待测样本急性髓细胞样白血病预后风险预测结果。

在一个实施例中，本发明提供一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据及临床病理特征；

处理单元，用于根据EGFR、PYGM和GPC3基因表达情况及临床病理特征，获得待测样本急性髓细胞样白血病预后风险预测结果。

本发明的一个目的在于提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的急性髓细胞样白血病预后风险预测。

术语“样本”可以是从受试者分离的任何生物样本。例如，样本可以包括但不限于体液、全血、血小板、血清、血浆、粪便、红细胞、白细胞或白血球、内皮细胞、组织活组织检查、滑液、淋巴液、腹水、间质或细胞外液、细胞间空间的液体，包括龈沟液、骨髓、脑脊液、唾液、粘液、痰、精液、汗液、尿液、鼻刷液、巴氏涂片液或任何其他体液。体液可以包括唾液、血液或血清。例如，多核苷酸可以是从体液例如血液或血清分离的无细胞DNA。样本也可以是肿瘤样本，肿瘤样本可以通过各种方法从受试者获得，所述方法包括但不限于静脉穿刺、排泄、射精、按摩、活组织检查、针抽吸、灌洗、刮擦、手术切口或介入或其他方法。样本可以是无细胞样本(例如，不包含任何细胞)。

本验证实施例的验证结果表明，为适应症分配固有权重相对于默认设置来说可以适度改善本方法的性能。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器(ROM，Read Only Memory)、随机存取存储器(RAM，RandomAccess Memory)、磁盘或光盘等。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上对本发明所提供的一种计算机设备进行了详细介绍，对于本领域的一般技术人员，依据本发明实施例的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (14)

1.一种急性髓细胞样白血病预后风险模型的构建方法，包括：

获取缺氧相关基因；

获取急性髓细胞样白血病基因表达数据作为训练数据集；

利用回归模型或机器学习，在训练数据集中，对所述缺氧相关基因进行特征选择，得到预后风险相关基因，所述利用回归模型是通过单因素COX回归模型对所述缺氧相关基因进行特征选择，得到预后风险相关基因，所述预后风险相关基因为EGFR、PYGM和GPC3基因；

利用所述预后风险相关基因，采用LASSOCOX回归分析或多因素COX回归分析，构建急性髓细胞样白血病预后风险模型；

将所述构建急性髓细胞样白血病预后风险模型作为特征向量，选取临床病理特征，利用所述特征向量和临床病理特征构建急性髓细胞样白血病预后综合风险模型。

2.根据权利要求1所述的构建方法，其特征在于，获取缺氧相关基因还包括如下处理：通过蛋白质相互网络、基因共表达网络、转录因子调控网络、pathway网络中的一种或几种基于网络分析方法筛选出排名靠前的缺氧相关基因。

3.根据权利要求1所述的构建方法，其特征在于，所述获取缺氧相关基因是获取GSEA数据库中缺氧相关基因。

4.根据权利要求1所述的构建方法，其特征在于，通过PPI网络分析方法筛选出排名靠前的缺氧相关基因。

5.根据权利要求1所述的构建方法，其特征在于，所述选取临床病理特征包括：采用回归模型或机器学习筛选出与急性髓细胞样白血病生存相关的临床病理特征。

6.一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器；

所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的基因表达数据及临床病理特征；

将所述待测样本的基因表达数据及临床病理特征输入权利要求1-5任意一项所述的急性髓细胞样白血病预后综合风险模型中；

获得待测样本急性髓细胞样白血病预后风险预测结果。

7.一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器；

所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的基因表达数据；

将所述待测样本的基因表达数据输入权利要求1-5任意一项所述的急性髓细胞样白血病预后风险模型中；

获得待测样本急性髓细胞样白血病预后风险预测结果。

8.一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的EGFR、PYGM和GPC3基因表达数据；

将所述EGFR、PYGM和GPC3基因表达数据输入权利要求1-5任意一项所述的急性髓细胞样白血病预后风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

9.一种急性髓细胞样白血病预后风险诊断设备，所述设备包括：存储器和处理器，所述存储器用于存储程序指令；

所述处理器用于调用程序指令，当程序指令被执行时，用于执行以下操作：

获取待测样本的EGFR、PYGM和GPC3基因表达数据和临床病理特征；

将所述EGFR、PYGM和GPC3基因表达情况和临床病理特征输入权利要求1-5任意一项所述的急性髓细胞样白血病预后综合风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

10.一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的基因表达数据及临床病理特征；

处理单元，用于将所述待测样本的基因表达数据及临床病理特征输入权利要求1-5任意一项所述的急性髓细胞样白血病预后综合风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

11.一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的基因表达数据；

处理单元，用于将所述待测样本的基因表达数据输入权利要求1-5任意一项所述的急性髓细胞样白血病预后风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

12.一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据及临床病理特征；

处理单元，用于将所述EGFR、PYGM和GPC3基因表达数据及临床病理特征输入权利要求1-5任意一项所述的急性髓细胞样白血病预后综合风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

13.一种急性髓细胞样白血病预后风险诊断系统，包括：

获取单元，用于获取待测样本的EGFR、PYGM和GPC3基因表达数据；

处理单元，用于将所述EGFR、PYGM和GPC3基因表达数据输入权利要求1-5任意一项所述的急性髓细胞样白血病预后风险模型中，获得待测样本急性髓细胞样白血病预后风险预测结果。

14.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-5任意一项所述的急性髓细胞样白血病预后风险模型的构建方法。

Images