一种克服信号失真的肿瘤细胞异常通路识别方法技术

本发明专利技术提供一种克服信号失真的肿瘤细胞异常通路识别方法，包括：利用非负矩阵分解算法对基因表达谱图像矩阵进行通路解构，得到通路表达谱矩阵；根据接受者操作特性曲线和所述通路表达谱矩阵，绘制出每条通路的活性异常曲线；根据所述每条通路的活性异常曲线，计算出每条通路的活性异常值；将每条通路的活性异常值与预设阈值进行比较，当其活性异常值高于预设阈值时，则该通路为肿瘤细胞活性异常通路。本发明专利技术克服了传统方法对通路活性建模的失真缺陷，消除了通路关联噪声对通路活性评估的影响，具有更好的肿瘤细胞异常通路检测精度，稳定性与可靠性。本发明专利技术可用于药物研发、个性化用药中的药物靶标的检测与识别。

【技术实现步骤摘要】

本专利技术涉及肿瘤细胞标志物检测
，具体是一种克服信号失真的肿瘤细胞异常通路识别方法。
技术介绍
肿瘤已经成为人类的“头号杀手”，严重威胁着人类的健康。由于生命系统的复杂性以及肿瘤细胞的高度异质性，如何有效诊断和治疗癌症是一项世界性的难题。随着现代生物技术的发展，在基因水平上识别和发现肿瘤标志物已经成为当前医学界的共识。细胞中的基因往往以一种通路结构(由多个基因做构成)形式来执行特定的生物学功能，最终影响和决定细胞的命运(表型)，比如正常细胞的癌变事件等。这些通路模块本质上构成一种复杂的多基因网络。越来越多的医学研究证实，是异常的通路活性而非单个基因的异常表达构成了细胞发生癌变的根本原因。因此，有效识别肿瘤细胞中的异常通路在肿瘤的早期诊断与治疗以及靶向药物的研制等方面都具有重要的应用价值。国内外的研究者针对肿瘤细胞异常通路的识别已经进行了一定的研究和探讨。其中，经典的方法有基因集合富集法(GSEA)、全局检验法(Globaltest)、PathOlogist、信号影响分析法(SPIA)和差异表达分析法(DEAP)。这些方法的思路可概括为三类：1)基于通路中差异基因的富集性，比如GSEA方法；2)基于通路中基因对表型的统计回归能力，比如Globaltest方法；3)基于通路的网络拓扑特性，比如PathOlogist、SPIA和DEAP等方法。然而，这些方法都未能考虑通路间的相互作用对通路活性信号检测的影响。不同的通路行使不同的生物学功能，但这些通路可能共享同样的基因，在功能上相互关联，从而在更高层次上构成复杂的生物分子系统，产生复杂多样的细胞表型。从系统生物学角度来看，只有有效去除通路相互作用对通路活性检测的噪声影响才能可靠识别肿瘤细胞异常通路。上述现有方法简单假设通路间是独立的，结果导致通路活性信号检测失真。因此，这些现有的方法仅仅在理想的孤立通路活动情况下是有效的，而并不适合在真实细胞情况下(即通路间具有复杂关联关系时)的异常通路活性检测与识别。本专利技术采用一种非负矩阵分解算法实现肿瘤细胞异常通路识别中的通路关联噪声的提取与去除，从而提高了肿瘤细胞异常通路的识别精度和稳定性。
技术实现思路
本专利技术的目的在于提供一种克服信号失真的肿瘤细胞异常通路识别方法，避免现有肿瘤细胞异常通路识别方法的不足，利用非负矩阵分解来实现通路间关联噪声的提取与去除，获得真实的通路活性，进而更加精确、可靠地识别肿瘤细胞中的异常通路。本专利技术的技术方案为：一种克服信号失真的肿瘤细胞异常通路识别方法，包括以下步骤：(1)利用非负矩阵分解算法对基因表达谱图像矩阵进行通路解构，得到通路表达谱矩阵；(2)根据接受者操作特性曲线和所述通路表达谱矩阵，绘制出每条通路的活性异常曲线；(3)根据所述每条通路的活性异常曲线，计算出每条通路的活性异常值；(4)将每条通路的活性异常值与预设阈值进行比较，当其活性异常值高于预设阈值时，则该通路为肿瘤细胞活性异常通路。所述的克服信号失真的肿瘤细胞异常通路识别方法，所述步骤(1)采用以下公式实现：minW,H(||E-WH||2+λ|W|1),s.t.W≥0,H≥0]]>其中，E表示基因表达谱图像矩阵，每一行表示一个基因，每一列表示一个细胞样本；W表示通路-基因关系矩阵，其中行表示基因，列表示通路，每一个元素表示对应基因与通路间的关联系数；H表示通路表达谱矩阵，每一行表示一个通路，每一列表示一个细胞样本；λ表示稀疏约束系数。所述的克服信号失真的肿瘤细胞异常通路识别方法，所述公式的求解包括以下步骤：a、根据通路基因成员表对通路-基因关系矩阵W进行0-1初始化，对通路表达谱矩阵H进行随机初始化；b、利用乘性更新算法更新通路-基因关系矩阵W和通路表达谱矩阵H：Wij←Wij(EHT)ij(WHHT+λW)ij]]>Hjk←Hjk(WTE)jk(WTWH)jk]]>其中，Wij表示通路-基因关系矩阵W中位于第i行、第j列的元素，其含义是基因i与通路j之间的关联系数，i＝1，2，…，L，j＝1，2，…，M；Hjk表示通路表达谱矩阵H中位于第j行、第k列的元素，其含义是通路j在样本细胞k中的活性水平，j＝1，2，…，M，k＝1，2，…，N；c、根据准则|E-WtHt|＜10-4，其中Wt、Ht分别为第t步迭代所获得的W和H，判断迭代是否收敛，若是，则执行步骤d，若否，则返回步骤b；d、得到通路表达谱矩阵H。所述的克服信号失真的肿瘤细胞异常通路识别方法，所述步骤(2)具体包括以下步骤：a、根据所述通路表达谱矩阵，找出每条通路在各个样本细胞中的活性水平最小值和最大值：Bu=mink=1N(Hjk),Bw=maxk=1N(Hjk)]]>其中，Hjk表示通路j在样本细胞k中的活性水平，Bu表示通路j在各个样本细胞中的活性水平最小值，Bw表示通路j在各个样本细胞中的活性水平最大值；b、在Bu与Bw之间均匀选择m个阈值常数h，采用以下公式对每个h分别计算其对应的fpr和tpr值，Bu＜h＜Bw，m≥100：fpr=FPFP+TP]]>tpr=TPFP+TP]]>其中，Lk表示样本细胞k的标签，Lk＝0表示样本细胞k为正常细胞，Lk＝1表示样本细胞k为肿瘤细胞，当Hjk＞h&Lk＝0为真时，I(Hjk＞h&Lk＝0)＝1，当Hjk＞h&Lk＝0为假时，I(Hjk＞h&Lk＝0)＝0，当Hjk＞h&Lk＝1为真时，I(Hjk＞h&Lk＝1)＝1，当Hjk＞h&Lk＝1为假时，I(Hjk＞h&Lk＝1)＝0；c、以计算出的fpr为横坐标、tpr为纵坐标绘制每条通路的活性异常曲线(fpr，tpr)。所述的克服信号失真的肿瘤细胞异常通路识别方法，采用以下公式计算每条通路的活性异常值Y：Y=∫BuBwtpr(h)fpr′(h)dh]]>其中，fpr′(h)表示fpr(h)的一阶导数。所述的克服信号失真的肿瘤细胞异常通路识别方法，通过以下步骤近似求解每条通路的活性异常值Y：a、对于每条通路，根据m个fpr的升序方式绘制新的通路活性异常曲线(fq，tq)，其中fq表示按照升序排列的第q(本文档来自技高网...

【技术保护点】
一种克服信号失真的肿瘤细胞异常通路识别方法，其特征在于，包括以下步骤：(1)利用非负矩阵分解算法对基因表达谱图像矩阵进行通路解构，得到通路表达谱矩阵；(2)根据接受者操作特性曲线和所述通路表达谱矩阵，绘制出每条通路的活性异常曲线；(3)根据所述每条通路的活性异常曲线，计算出每条通路的活性异常值；(4)将每条通路的活性异常值与预设阈值进行比较，当其活性异常值高于预设阈值时，则该通路为肿瘤细胞活性异常通路。

【技术特征摘要】
1.一种克服信号失真的肿瘤细胞异常通路识别方法，其特征在
于，包括以下步骤：
(1)利用非负矩阵分解算法对基因表达谱图像矩阵进行通路解
构，得到通路表达谱矩阵；
(2)根据接受者操作特性曲线和所述通路表达谱矩阵，绘制出
每条通路的活性异常曲线；
(3)根据所述每条通路的活性异常曲线，计算出每条通路的活
性异常值；
(4)将每条通路的活性异常值与预设阈值进行比较，当其活性
异常值高于预设阈值时，则该通路为肿瘤细胞活性异常通路。
2.根据权利要求1所述的克服信号失真的肿瘤细胞异常通路识
别方法，其特征在于，所述步骤(1)采用以下公式实现：
minW,H(||E-WH||2+λ|W|1),s.t.W≥0,H≥0]]>其中，E表示基因表达谱图像矩阵，每一行表示一个基因，每一列表
示一个细胞样本；W表示通路-基因关系矩阵，其中行表示基因，列
表示通路，每一个元素表示对应基因与通路间的关联系数；H表示通
路表达谱矩阵，每一行表示一个通路，每一列表示一个细胞样本；λ表
示稀疏约束系数。
3.根据权利要求2所述的克服信号失真的肿瘤细胞异常通路识
别方法，其特征在于，所述公式的求解包括以下步骤：
a、根据通路基因成员表对通路-基因关系矩阵W进行0-1初始化，
对通路表达谱矩阵H进行随机初始化；
b、利用乘性更新算法更新通路-基因关系矩阵W和通路表达谱矩
阵H：
Wij←Wij(EHT)ij(WHHT+λW)ij]]>Hjk←Hjk(WTE)jk(WTWH)jk]]>其中，Wij表示通路-基因关系矩阵W中位于第i行、第j列的元素，其
含义是基因i与通路j之间的关联系数，i＝1，2，…，L，j＝1，2，…，M；Hjk表示通路表达谱矩阵H中位于第j行、第k列的元素，其含义是通路j在
样本细胞k中的活性水平，j＝1，2，…，M，k＝1，2，…，N；
c、根据准则|E-WtHt|＜10-4，其中Wt、Ht分别为第t步迭代
所获得的W和H，判断迭代是否收敛，若是，则执行步骤d，若否，
则返回步骤b；
d、得到通路表达谱矩阵H。
4.根据权利要求1所述的克服信号失真的肿瘤细胞...

【专利技术属性】
技术研发人员：王红强，
申请(专利权)人：中国科学院合肥物质科学研究院，
类型：发明
国别省市：安徽;34