发明内容
在本发明中,发明人从TCGA公共数据库内结肠腺癌病人的RNA-seq的基因表达数据及其临床预后信息出发,运用两步的有监督的机器学习方法,无偏地从整个转录组中获得一种12-基因表达分类器,准确预测结肠腺癌病人的预后。
一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,包括:数据训练阶段和验证阶段,所述训练阶段包括第一阶段和第二阶段,所述第一阶段使用有监督的机器学习方法建立能预测结肠腺癌病人预后的基因表达分类器雏形,所述第二阶段进一步使用机器学习的方法获得预测结肠腺癌病人预后的基因表达分类器。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述第一阶段通过使用结肠腺癌病人的基因表达信息和临床信息,采用有监督的机器学习方法挑选与真实预后情况高度相关的基因。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中对所有基因的表达和预后情况的Pearson系数的绝对值由大到小排序,获得与真实三年内复发情况高度相关的基因。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中使用留一交叉验证方法获得最佳的基因类别和数目,从而构建所述基因表达分类器雏形。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中使用基因表达分类器雏形计算用于验证的肿瘤样本的风险系数,预测其复发风险,通过比较真实复发风险和预测复发风险的一致性,验证所述基因表达分类器雏形的效能。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述第二阶段包括特征(基因)排序,所述特征(基因)排序基于单因素的Cox比例风险回归模型所得p值进行。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述第二阶段采用KM生存分析计算高风险组和低风险组之间的时序检验p值,时序检验p值最小的Cox模型所包含的基因种类和数目就是所构建的基因表达分类器。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中验证阶段使用基因表达芯片数据集验证所获得的基因表达分类器。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器及其构建方法。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器及其构建方法,实现对结肠腺癌病人的无复发生存期和总体生存期进行精准预测。
另一方面,本发明的目的是提供一种预测结肠腺癌病人预后的基因表达分类器及其构建方法,将结肠腺癌病人分为差预后和预后良好的两个亚群,为其辅助治疗提供判断依据。
另一方面,本发明的目的是提供一种在多个结肠腺癌数据集中预测结肠腺癌病人预后的基因表达分类器及其构建方法。
另一方面,本发明的目的是提供一种在结肠腺癌数据集中预测结肠腺癌病人预后的基因表达分类器及其构建方法,所述多个结肠腺癌数据集选自GSE39582、GSE17538。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其包括:数据训练阶段和验证阶段。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述训练阶段包括第一阶段和第二阶段。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述第一阶段包括:使用TCGA结肠腺癌病人的基因表达信息和临床信息,使用有监督的机器学习方法建立能预测结肠腺癌病人预后的基因表达分类器雏形。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中有监督的机器学习的方法如下:
分组:对满足条件的结肠腺癌病人的基因表达数据在分组前进行标准化处理,然后分别从临床上差预后和预后良好的病人中随机选取第一数量病人和第二数量作为训练组,剩余病人作为验证组;
筛选与真实预后情况高度相关的基因:从训练组数据着手,差预后的病人真实预后状态赋值1,预后良好的病人真实预后状态赋值0。在训练组队列中,计算每个基因的表达值与真实预后状态的Pearson相关系数,选择绝对Pearson相关系数大于或等于0.3的基因,为了排除获得的回归系数分布是随机概率产生的可能性,运用置换检验方法产生了10,000个Monte-Carlo模拟试验,从而对训练组的病人的基因表达数据和预后分组信息进行随机化处理;
有监督的分类方法:将挑选出来的基因按照相关系数的绝对值从大到小进行排序,得到一个排行榜,从排在最前面的两个基因开始,每次从排行榜再依次添加两个基因建立一个分类器,如此循环,直到排行榜中所有的基因都被作为报告子用尽,建立基因表达分类器雏形。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中使用留一交叉验证法检查所述基因表达分类器雏形的效能,步骤如下:
第一步,将一名病人留出,以备后续检验用;
第二步,通过取基因表达Z-分数的平均值计算在差预后组该分类器中所涉及的所有基因的表达模式;类似地计算在预后良好组的基因表达模式;接着,定义一个风险系数(risk-coef),所述风险系数即指该肿瘤内预后良好组的基因表达模式的相关系数减去差预后组基因表达模式的相关系数;
第三步,计算剩余的训练肿瘤样本以及第一步中留存的肿瘤样本的风险系数,将这些样本按照风险系数从小到大排序,第一数量肿瘤病人被划分为高基因组风险组,第二数量病人被划分为低基因组风险,检查每个病人真实的临床预后情况和预测的基因组风险的一致性;
循环第一到第三步,直到所有训练病人样本都被留出过一次,每次当被留出的样本的基因组风险和实际的临床预后情况相背离时,错误计数器加1次;
对于每一个分类器,设置风险系数的门槛值,将病人分为第一数量高风险和第二数量低风险的病人;同时,计算独立于训练样本的验证样本的风险系数;通过前面设置的风险系数的门槛值,可判定验证样本的基因组风险的高低;同样地,也统计验证组病人的基因组风险和实际临床预后出现不一致的次数;最终发现分别包含排行榜的前16、36、40、42、44、46、48、50、56、58、60、62、64、66、68、70、72、74、76、78、80、82、84、86或118基因时的分类器预测错误次数为最低,运用所有病人的风险系数和真实预后状态进行接受者操作特性曲线分析,比较这25个分类器的准确度和敏感度。结果显示这25个分类器的准确度和敏感度没有显著差异。为了尽量不遗漏与预后高度相关的基因,排行榜的前118个基因都被选择进入第二阶段的训练步骤,即选择118-基因表达分类器作为雏形。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中运用一种改良版的留一交叉验证法来解决信息遗露(information leakage)的问题:
第一步,留出一个样本待验证;
第二步,运用剩余样本计算所有基因的表达和预后情况的Pearson系数,过滤得到|coef|≥0.3的基因;
第三步,运用第二步中过滤得到的基因构建分类器,并据此预测被留出的那个样本的基因组风险。
第四步,重复第一到第三步,直到所有肿瘤病人都被留出过一次,从而获得基因表达分类器。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述第二阶段包括:
基于训练的第一阶段获得的118-基因表达分类器,进一步运用机器学习的方法获得更加简洁的风险评分***来预测结肠腺癌病人的预后,同样是运用TCGA的病人数据,但是这次包括了在三年内删失的病例,在建模过程中,基因被称为特征;
随机分组:将这些样本随机划分为:组1和组2。
特征排序:以组1为训练数据,通过单因素的Cox比例风险回归模型计算单个特征的回归系数和p值,按照p值从小到大将118个特征重新排序,排在越前面的特征,Cox回归p值越小,与预后的相关性越大;
特征数目的优化:从排序后的第一个特征开始,从前往后每次加一个特征,运用多因素的Cox回归分析获得各个特征的Cox回归系数;运用组2病人进行交叉验证,评估此Cox模型的好坏:将组2中每个病人的相关特征的表达值与多因素Cox回归系数相乘并累加得到一个分数值,分数值高低表示病人死亡或者复发风险的高低;接着采用KM生存分析计算组2的高风险和低风险亚组之间的时序检验p值;如此循环,直到所有的特征都被纳入Cox回归模型;时序检验p值最小的Cox模型所包含的特征种类和数目就是最优的,获得基因表达分类器。
另一方面,本发明提供一种预测结肠腺癌病人预后的基因表达分类器的构建方法,其中所述验证阶段包括:
使用GEO数据库内的结肠腺癌病人的基因表达芯片数据集验证以上所获得基因表达分类器;在一个独立的数据集中,那些分数值高于群体分数值中位数的病人被划为高风险组,而其余为低风险组;KM分析用来比较高风险组和低风险组的生存曲线;时序检验p值<0.05表示有统计学差异。
在另一实施方式中,本发明提供一种12-基因表达分类器。
在另一实施方式中,本发明提供一种12-基因表达分类器,其中12基因表达分类器的基因选自由TREML2、PADI4、NCKIPSD、PTPRN、PGLYRP1、C5orf53、TREML3、NOG、VIP、FAM171B、NKAIN4和RIMKLB组成的组。
另一方面,本发明提供一种基因表达分类器,其能作为一种有效的诊断手段将结肠腺癌病人中有较高风险发展成差预后的亚群分离出来。
另一方面,本发明提供一种基因表达分类器,其能作为一种有效的诊断手段将结肠腺癌、肾癌、肺鳞癌、眼底和皮肤黑色素瘤、脑癌以及胰腺癌病人中有较高风险发展成差预后的亚群分离出来。
有益效果
本发明运用有监督的机器学习的方法建立了一种基因表达分类器来精准预测结肠腺癌预后。该基因表达分类器具有以下几种明显的优点:
1.能预测结肠腺癌病人的预后。可为结肠腺癌病人术后辅助治疗提供精准判断依据。
2.不仅可预测结肠腺癌病人的预后,在包括肾癌、肺鳞癌、眼底和皮肤黑色素瘤、脑癌以及胰腺癌中都有较为显著的预后预测功效。
具体实施方式
本发明将会参照下面的实施例进行阐述,但本发明将不限于下面的实施例。
1研究材料和方法
1.1TCGA(The Cancer Genome Altas)和GEO(Gene
Expression Omnibus)数据集的下载
从TCGA RNA-seq数据库(基因https://cancergenome.nih.gov/)下载多种癌症的RNA-seq转录组数据机器临床信息。结肠癌基因表达芯片及其临床信息从GEO数据库下载(https://www.ncbi.nlm.nih.gov/geo/)。
1.2基因表达分类器的构建
基因表达分类器的构建方法分为训练和验证两个步骤。
1.2.1训练步骤
1.2.1.1第一阶段
分组:TCGA的结肠腺癌(COAD)的基因表达数据和临床信息被用来进行第一阶段的机器学习,获得基因表达分类器雏形。满足条件的结肠腺癌病人共91位,其中42位病人在3年内(<36个月)发生了复发事件,为差预后;49位病人在等于或长于3年(>=36个月)的期间内未发生复发事件,为预后良好。这些病人的基因表达值在分组前需进行标准化(取Z-分数)。分别从差预后和预后良好的病人中随机取得32和39位病人,作为训练组。余下的病人作为验证组。
与预后高度相关的基因的筛选:RNA-seq转录组总共有20530个基因。从训练组数据着手,差预后的病人真实预后状态赋值1,预后良好的病人真实预后状态赋值0。在训练组队列中,计算每个基因的表达值与真实预后状态的Pearson相关系数。总计1510个基因的绝对Pearson相关系数大于或等于0.3。为了排除获得的回归系数分布是随机概率产生的可能性,运用置换检验方法产生了10,000个Monte-Carlo模拟试验,从而对训练组的71个病人的基因表达数据和预后分组信息进行随机化处理。在10,000个Monte-Carlo试验中,|相关系数|≥0.3的基因数目的试验的频数分布见图1B。获得|相关系数|≥0.3的基因数目≥1510的概率为0.0019(p<0.05),据此可拒绝零假设。
有监督的分类方法:很显然,1510个基因对于分类器来说数目太过庞大,基因数目需要优化。将这1510个基因按照相关系数的绝对值从大到小进行排序,得到一个排行榜(Top list)。从排在最前面的两个基因开始,每次从排行榜再依次添加两个基因建立一个分类器,如此循环,直到排行榜中所有的1510个基因都被作为报告子(reporter)用尽。因此总共建立了755个分类器。
一种称为留一交叉验证(leave-one-out cross-validation,LOOCV)的方法被用来检验这些分类器的效能,步骤如下:
第一步,将一名病人留出,以备后续检验用。
第二步,计算在差预后组该分类器中所涉及的所有基因的表达模式(poor-prognosis expression template)(取表达Z分数的平均值);类似地,计算在预后良好组的基因表达模式(good-prognosis expression template)。接着,定义一个风险系数(risk-coef)。对于一个肿瘤来说,风险系数即指该肿瘤内预后良好组的基因表达模式的Pearson相关系数减去差预后组基因表达模式的Pearson相关系数:
风险系数=预后良好组的基因表达模式的pearson相关系数-差预后组基因表达模式的pearson相关系数
第三步,计算剩余的70个训练肿瘤样本以及第一步中留存的肿瘤样本的风险系数。将这71个样本按照风险系数从小到大排序,前面32个肿瘤病人被划分为高基因组风险(high genomic risk)组,而剩余的39个病人被划分为低基因组风险(low genomic risk)组。检查每个病人真实的临床预后情况和预测的基因组风险的一致性。
循环第一到第三步,直到所有的71个训练病人样本都被留出过一次。每次当被留出的样本的基因组风险和实际的临床预后情况相背离时,错误计数器加1次。
在留一交叉验证过程中,最好的分类器应是,错误计数器收集的预测错误次数应该最少。当基因表达分类器包含排行榜的前16、36、40、42、44、46、48、50、56、58、60、62、64、66、68、70、72、74、76、78、80、82、84、86或118基因时,留一交叉验证过程中的预测错误次数最低。还需要进一步判断这25个分类器中哪一个最好。对于每一个分类器,设置风险系数的门槛值,将71个病人分为32个高风险和39个低风险的病人。同时,计算独立于训练样本的20个验证样本的风险系数。通过前面设置的风险系数的门槛值,我们可判定这20个验证样本的基因组风险的高低。同样地,我们也统计验证组病人的基因组风险和实际临床预后出现不一致的次数。对于每个分类器,运用这91个TCGA病人的风险系数和真实预后状态进行接受者操作特性曲线分析(Receiver Operating Characteristic Curve,ROC),比较这25个分类器的准确度和敏感度。结果显示这25个分类器的准确度和敏感度没有显著差异(图2和表1)。为了尽量不遗漏与预后高度相关的基因,排行榜的前118个基因都被选择进入第二阶段的训练步骤。
表1 ROC分析的统计表
备注:AUC:曲线下的面积SE:标准误差95%CI:95%置信区间
无信息遗漏(information leak)的交叉验证:由于以上1510个基因是基于所有的71个训练样本获得的(包括被留出的用来验证的那个样本),因此存在信息遗漏(information leakage)所产生的过度拟合(Over-fitting)的可能性。为此,发明人运用一种改良版的留一交叉验证来评估过度拟合的程度:
第一步,留出一个样本待验证。
第二步,运用剩余的70个样本计算所有基因的表达和预后情况的Pearson系数。过滤得到|相关系数|≥0.3的基因。
第三步,运用第二步中过滤得到的基因构建分类器,并据此预测被留出的那个样本的基因组风险。
第四步,重复第一到第三步,直到所有的71个肿瘤病人都被留出过一次。
因此又获得71个分类器。结果显示原初的1510个基因中绝大多数的基因都存在于这71个分类器中(图3)。据此,可判定前面的训练过程中所引入的过度拟合相当有限。
1.2.1.2第二阶段
基于训练的第一阶段获得的118个基因,进一步运用机器学***均数)。
随机分组
将这些样本随机划分为(随机化后两组病人各临床信息无显著差异):组1和组2(表2)。
表2.240名TCGA COAD病人临床病理特征在组1和组2数据集中的比较
t检验
Fisher精确检验
§卡方检验
特征排序
以组1为训练数据,通过单因素的Cox比例风险回归模型(Cox’s proportionalhazards regression model,CPH),计算单个特征的回归系数和p值。按照p值从小到大将118个特征重新排序。排在越前面的特征,Cox回归p值越小,与预后的相关性越大。
特征数目的优化:这是一种迭代优化的过程。从排序后的第一个特征开始,从前往后每次加一个特征,运用多因素的Cox回归分析获得各个特征的Cox回归系数。运用组2病人进行交叉验证,评估此Cox模型的好坏:将组2中每个病人的相关特征的表达值与多因素Cox回归系数相乘并累加得到一个分数值(score),分数值高低表示病人死亡或者复发风险的高低。接着采用KM(Kaplan-Meier)生存分析计算组2的高风险和低风险亚组之间的时序检验p值(log Rank p-value)。如此循环,直到所有的特征都被纳入Cox回归模型。时序检验p值最小的Cox模型所包含的特征种类和数目就是最优的。
基因表达分类器的确立:
通过以上步骤确定基因表达分类器的基因数目和种类(表3),及其表达加权系数(多因素Cox回归系数),建立基因表达评分体系:
Ei:基因i的表达水平;βi:基因i的多变量Cox回归系数
1.2.2验证阶段
GEO数据库内的结肠腺癌病人的基因表达芯片数据集被用来验证以上所获得12-基因表达分类器。因为每个基因可能有多个探针,选择单因素CPH p值最小的探针来代表该基因。类似的,也需要对每个基因的探针信号进行均一化处理(将每个基因的探针信号值除以TFRC、GUSB、和RPLP0的探针信号值的几何平均数):在一个独立的数据集中,那些分数值高于群体分数值中位数的病人被划为高风险组,而其余为低风险组。KM分析用来比较高风险组和低风险组的生存曲线。时序检验p值<0.05表示有统计学差异。
另外,TCGA数据库内其他癌种的数据集也被用来验证该12-基因模型。
2.研究结果
2.1分类器雏形的建立
运用TCGA结肠腺癌的数据,采用一种无偏的筛选方法获得了1510个与真实预后状态高度相关的基因(Pearson相关系数的绝对值大于等于0.3)。在后续的留一交叉验证过程中,构建了755个分类器。当分类器包含排行榜的前16、36、40、42、44、46、48、50、56、58、60、62、64、66、68、70、72、74、76、78、80、82、84、86或118基因时,预测的复发风险与真实的预后状态不相符合的累计次数最少。在这71个训练样本中,包含有前118个基因的分类器的预测准确性如图4。另外,这些分类器的准确度和敏感度在20个训练病人中也进行了验证(图2和图4,表1)。ROC分析发现,这25个分类器在预测这91个TCGA结肠腺癌病人(71个训练样本,20个验证样本)的准确度和敏感度都较为理想,且相互之间没有明显差异(图2和表1)。
2.2 12-基因表达分类器的建立
为了建立更加简洁有效的分类器,通过进一步的机器学习方法对从第一阶段的训练过程中获得的118个基因进行进一步筛选,获得了12-基因表达分类器模型。每个基因的表达加权系数如表3。每个病人都可计算出相应的12-基因表达风险值。KM生存分析的结果发现,在TCGA的训练组和验证组中,低12-基因表达风险值亚组的预后明显好于高风险值亚组(图5:DFS训练数据集:KM时序检验p=0.0001;验证数据集:KM时序检验p=0.0005)。
表3最终12-基因表达分类器模型中的基因及其表达加权系数。
2.3 12-基因表达分类器诊断预测效能验证
GSE39582和GSE17538被用来验证12-基因表达分类器的诊断预测效能。
GSE39582包含有566个结肠腺癌样本以及19个肺粘膜组织样本。Marisa等人通过这些样本开发了基因表达分类器,并定义了6种结肠腺癌的分子亚型,这些分子亚型具有不同的分子特征、临床以及预后特征(Marisa等,2013)。对于GSE39582中pMMR的病人,本发明构建的12-基因表达分类器可高效准确地将预后良好和预后较差的组别区分开[图6A和B,OS:KM时序p=0.005;无复发生存期(RFS):KM时序p=0.022]。
GSE17538中,Smith和Freeman等人通过运用结肠腺癌患者的基因芯片表达和临床预后信息,构建了一种特异表征癌症转移的基因表达模式,去预测结肠腺癌病人的复发和死亡风险(Smith等,2010;Freeman等,2012)。本发明构建的12-基因表达分类器可将GSE17538内232个结肠腺癌病人中预后良好和差预后的病人有效区分开[图6C-E,无疾病生存期(DFS):KM时序p=0.0004;疾病特异性存活(DSS):KM时序p=0.0034;总体生存期(OS):KM时序p=0.0336]。另外,晚期病人在高12-基因分值组得到显著富集(表4)。
表4高评分组和低评分组晚期病人的分布使用Fisher精确检验进行统计学分析
3.4 12-基因分类器预测其他癌症预后的能力
从TCGA数据库下载24种其他癌症的RNA-seq和临床信息,运用类似的方法来验证12-基因表达分类器在这些癌种中的预后预测能力。KM生存分析结果显示12-基因表达分类器可将结肠腺癌以外的多种癌症中的差预后和预后良好的病人区分开,这些癌种包括KIPAN(图7A,OS:KM时序p=6.815e-6)、KIRC(图7B,DFS:KM时序p=0.0480)、KIRP(图7C,DFS:KM时序p=0.0027;图7D,OS:时序p=0.0129)、LUSC(图7E,DFS;时序p=0.0071)、SKCM(图7F,DFS:时序p=0.01117)、LGG(图8A,OS:时序p=0.0031)、UVM(图8B,OS:时序p=0.0054)、GBM(图8C,OS:时序p=0.0074)、CESC(图8D,OS:时序p=0.0090)、PAAD(图8E,OS:时序p=0.0127)和STAD(图8F,OS:时序p=0.0456)。
4.讨论
为了准确预测结肠癌病人的预后,人们已经做了很多尝试去构建各种基因表达分类器。一篇荟萃分析对这些已经发表的基因表达分类器进行了总结和评估(Sanz-Pamplona等,2012),结果显示虽然这些分类器在各自的数据集中与预后高度相关,然而在其他独立的数据集中预测预后的准确度不高。因此需要更加强大更加准确的基因表达分类器在多种其他独立的数据集中也能准确预测结肠腺癌预后。通过独特的两步有监机器学习方法,我们获得一种12-基因表达分类器,该分类器的预测预后的准确度在两个独立的大型数据集中都得到了验证(GSE39582,N=459;GSE17538,N=232)。
结直肠癌术后辅助治疗的选择在临床和学术圈是近二十年来研究和讨论的热点(Dotan等,2011;Meropol,2011;Vachani等,2013)。苏格兰校际协作指南网(the ScottishIntercollegiate Guidelines Network,SIGN)、美国临床肿瘤学会(American Society ofClinical Oncology,ASCO)、美国国家综合癌症网络(National Comprehensive CancerNetwork,NCCN)遵循的临床指南也不尽相同(Gao等,2016)。就国内临床来讲,行手术的pMMR的结肠腺癌病人一般会进行5-FU单药的术后辅助化疗(Alex等,2016),由于缺少有效的临床检测方法进一步将肿瘤患者精细分群,因此很多结肠腺癌病患很可能被过度治疗。本发明开发的12-基因分类器可以将GSE39582数据集内pMMR病人预后良好和预后较差的病人有效区分开,说明该12-基因表达分类器可以开发成为一种体外诊断产品,为pMMR结肠腺癌病人术后辅助治疗的选择提供精准判断依据。
该分类器所包含的12个基因中,有7个编码的蛋白是免疫***调控相关,说明肿瘤内免疫相关因子的不正常表达往往预示着该肿瘤预后较差,这也可以解释我们构建的基因表达分类器也能预测多种其他癌症的预后。
该发明构建的12-基因表达分类器提供了一种新型的评分***,其中基因的表达值是相对表达量。相对表达量是采用TFRC,GUSB和RPLP0表达值的几何平均数来做均一化处理获得的。为了最大程度地保留不同肿瘤间的异质性,β-肌动蛋白(β-Actin,ACTB)和甘油醛-3-磷酸脱氢酶(glyceraldehyde-3-phosphate dehydrogenase,GAPDH)没有被用来作为内参基因进行基因的均一化处理,这是因为细胞骨架重塑和能量代谢的异常调控在肿瘤发展和转移过程中起着非常重要的功能,这两种关键途径的信号水平在不同的癌症患者之间往往存在很大的差异。基因表达的均一化处理非常重要,这使得该基因表达评分***灵活适用于不同的基因表达检测***,如qPCR、RNA-seq以及QuantiGene Plex等。另外本发明的12-基因表达分类器不仅可以预测结肠腺癌患者的预后,在包括肾癌、肺鳞癌、眼底和皮肤黑色素瘤、脑癌以及胰腺癌中都有较为显著的预后预测功效。
综上所述,本发明构建了一种12基因表达分类器***可准确预测结肠腺癌病人的预后,可为结肠腺癌患者术后辅助治疗的选择提供精准判断依据。
应当强调,本发明的上述实施例仅仅是可能的示例实施方式,其仅仅是为了清楚地理解本公开的原理而提出的。在不脱离本公开的精神和原理的情况下,可以对本公开的上述实施例进行许多变化和修改。所有这些修改和变化旨在被包括在本发明的范围内并由所附权利要求保护。
参考文献
1.Alex AK,Siqueira S,Coudry R,Santos J,Alves M,Hoff PM,RiechelmannRP.Response to Chemotherapy and Prognosis in Metastatic Colorectal CancerWith DNA Deficient Mismatch Repair.Clin Colorectal Cancer.2016;pii:S1533-0028(16)30255-9.
2.BrychtováV, R,Hrstka R,VídeňskáP,BencsikováB,HanákováB,DubskáL,Nenutil R,BudinskáE.Molecular Pathology of ColorectalCancer,Microsatellite Instability-the Detection,the Relationship to thePathophysiology and Prognosis.Klin Onkol.2017Spring;30(Supplementum1):153-155.Czech.
3.Cunningham D,Atkin W,Lenz HJ,Lynch HT,Minsky B,Nordlinger B,Starling N.Colorectal cancer.Lancet.2010;375:1030-47.
4.De Sousa E,Melo F,Wang X,Jansen M,et al..Poor-prognosis coloncancer is defined by a molecularly distinct subtype and develops fromserrated precursor lesions.NAT MED 2013;19:614-8.
5.Dotan E,Cohen SJ.Challenges in the management of stage II coloncancer.Semin Oncol.2011;38(4):511-20.
6.Ebata T,Hirata H,Kawauchi K.Functions of the Tumor Suppressorsp53and Rb in Actin Cytoskeleton Remodeling.Biomed Res Int.2016;2016:9231057.doi:10.1155/2016/9231057.
7.Freeman TJ,Smith JJ,Chen X,Washington MK,Roland JT,Means AL,Eschrich SA,Yeatman TJ,Deane NG,Beauchamp RD.Smad4-mediated signalinginhibits intestinal neoplasia by inhibiting expression ofβ-catenin.Gastroenterology.2012 Mar;142(3):562-71.
8.Gao S,Tibiche C,Zou J,Zaman N,Trifiro M,O'Connor-McCourt M,WangE.Identification and construction of combinatory cancer hallmark-based genesignature sets to predict recurrence and chemotherapy benefit in stage IIcolorectal cancer.JAMA Oncol.2016;2(1):37-45.
9.Marisa L,de Reyniès A,Duval A,Selves J,Gaub MP,Vescovo L,Etienne-Grimaldi MC,Schiappa R,Guenot D,Ayadi M,Kirzin S,Chazal M,Fléjou JF,BenchimolD,Berger A,Lagarde A,Pencreach E,Piard F,Elias D,Parc Y,Olschwang S,Milano G,Laurent-Puig P,Boige V.Gene expression classification of colon cancer intomolecular subtypes:characterization,validation,and prognostic value.PLoSMed.2013;10(5):e1001453.
10.Meropol N,Lyman GH,Chien R.Use of a multigene prognostic assay forselection of adjuvant chemotherapy in patients with stage II colon cancer:Impact on quality-adjusted life expectancy and costs.Poster presented at theAmerican Society Clinical Oncology Gastrointestinal Cancers Symposium;SanFrancisco,CA;January 2011.
11.Meropol NJ.Ongoing challenge of stage II colon cancer.J ClinOncol.2011;29(25):3346-3348.
12.Moloney JN,Cotter TG.ROS signalling in the biology of cancer.SeminCell Dev Biol.2017 Jun 3.pii:S1084-9521(16)30383-4.doi:10.1016/j.semcdb.2017.05.023.
13.Quasar Collaborative Group,Gray R,Barnwell J,McConkey C,Hills RK,Williams NS,Kerr DJ.Adjuvant chemotherapy versus observation in patients withcolorectal cancer:a randomised study.Lancet.2007,370(9604):2020-9.
14.Sanz-Pamplona R,Berenguer A,Cordero D,Riccadonna S,SoléX,Crous-BouM,GuinóE,Sanjuan X,Biondo S,Soriano A,Jurman G,Capella G,Furlanello C,MorenoV.Clinical value of prognosis gene expression signatures in colorectalcancer:a systematic review.PloS One.2012;7(11):e48877.
15.Smith JJ,Deane NG,Wu F,Merchant NB,Zhang B,Jiang A,Lu P,JohnsonJC,Schmidt C,Bailey CE,Eschrich S,Kis C,Levy S,Washington MK,Heslin MJ,CoffeyRJ,Yeatman TJ,Shyr Y,Beauchamp RD.Experimentally derived metastasis geneexpression profile predicts recurrence and death in patients with coloncancer.Gastroenterology.2010 Mar;138(3):958-68.
16.Stewart BW,Wild CP,editors.World Cancer Report 2014.2014.Lyon,France:International Agency for Research on Cancer.
17.Vachani C,Giantonio B.Stage II colon cancer:To treat or not totreat?2013(http://www.oncolink.org/types/article.cfm?c=124&id=9621).
18.van‘t Veer LJ,Dai H,van de Vijver MJ,Friend SH,Hart AAM,Mao M,etal.,Gene expression profiling predicts clinical outcome of breast cancer,Nature 415(2002)530-536.
19.Venook A,Niedzwiecki D,Lopatin M,et al.Validation of a 12-genecolon cancer recurrence score(RS)in patients(pts)with stage II colon cancer(CC)from CALGB 9581.Presented at:American Society of Clinical Oncology;June2011;Chicago,IL.