CN107220525A - 基于rnn的基因调控网络构建与动态差异性分析方法 - Google Patents
基于rnn的基因调控网络构建与动态差异性分析方法 Download PDFInfo
- Publication number
- CN107220525A CN107220525A CN201710355357.5A CN201710355357A CN107220525A CN 107220525 A CN107220525 A CN 107220525A CN 201710355357 A CN201710355357 A CN 201710355357A CN 107220525 A CN107220525 A CN 107220525A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- gene
- network
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Bioethics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于RNN的基因调控网络构建与动态差异性分析方法，包括以下步骤：第一步、基于deepRNN的基因动态调控网络构建；第二步、基于亚型内动态调控网络的时序变化演化分析；第三步、基于亚型间动态调控网络的网络差异演化分析，不同亚型网络的演化分析包括动力学分析、差异性分析和扰动分析。本发明提供一种精确性较好的基于RNN的基因调控网络构建与动态差异性分析方法。
Description
技术领域
本发明涉及一种基于RNN的基因调控网络构建与动态差异性分析方法。
背景技术
基因调控网络的建模与演化分析能够很好的挖掘基因表达数据中的深层信 息，是当前生物信息学研究的重要领域和关键问题。二十世纪90年代以来，随着 基因芯片技术的发展和二代测序技术的兴起，基因调控网络建模的研究取得了巨 大进展。
基因调控网络建模主要根据基因表达数据推理网络中的调控关系，并表示为 拓扑结构，属于依靠数据挖掘进行的逆向工程研究。构建基因调控网络首先需要 确定网络模型，然后根据模型选择合适的建模算法。经典的网络模型包括布尔网 络、关联网络、微分方程、贝叶斯网络。
(a)布尔网络。布尔网络对基因状态做了相应简化，用布尔函数代替了微分 和导数描述基因间的相互关系。该模型的缺点在于不精确性，仅仅通过使用固定 的逻辑规则刻画和反映基因间相互作用，并不能准确描述真实的基因调控网络拓 扑结构，而且对基因数据进行离散化时不可避免的会造成很多重要的表达信息丢 失。Kauffman等人最先提出了布尔网络的分析框架模型，随后Akusu等人对布尔 网络在推理过程中的最少样本数进行证明。Liang等人设计了REVEAL算法，在原 有的离散化模型上尽可能少的减少信息损失。此外，Lyla等人提出了一种新的概 率布尔网络(PBN)，这是对传统布尔网络的拓展，同时量化基因间作用关系和灵 敏度从而解决模型选择过程中的不确定性，提高了模型的精确性。
(b)关联网络。关联网络的建模主要通过基因表达数据间的关联度实现。通 常使用互信息、皮尔森相关系数等测度计算基因间的相似度，若基因对间的相似 度高于某一阈值，则该基因对在网络中直接连通。Butte等人首先利用互信息计算 所有基因对之间的关联度，然后设置互信息阈值。后来发现，若基因对间具有相 同或相近的调控机制，则两个基因的关联度较高，尤其是同一转录因子的靶基因 或同一条生物通路上的基因。Margolin等人提出ARCANE方法，通过信息论构建 关联网络，该方法的优点是模型的建立简单易操作，但是构建的网络存在很多假 阳性的边。为降低所构建网络结构的假阳性率，得到接近真实拓扑的调控网络， 一般在计算基因对间的关联度时隔绝其它基因的影响。
(c)贝叶斯网络。贝叶斯网络(BN)通过局部概率的乘积来近似描述整体 网络结构的复杂概率分布，属于概率图模型，将节点之间的连边表示为节点间存 在的概率依赖关系。动态贝叶斯网络(DBN)是对静态贝叶斯网络模型的扩展， 通过引入时间因素形成动态变化网络，更加真实地表示随机系统的动态性。基因 调控网络本质上是一个复杂而连续的动态网络系统，所以在具体建模的时候，往 往对DBN进行简化从而降低计算复杂度。DBN克服了静态BN有向无环的不足， 更好地刻画了基因调控网络的动态特性，提高了模型的预测精度。Norbert为了能 够从基因扰动型实验数据中学习动态贝叶斯网络，利用离散化方法来对基因表达 数据进行预处理，结合基因调控的负反馈与时延因素提出新的数据整合模型，利 用并行算法加速构建基因调控网络。
随着2006年Hinton教授在《科学》上的一篇文章，深度学习拉开帷幕，并在 各个领域表现不俗。同时，学术界和行业都强调了深度学习的洞察力在生物信息 学中的应用，例如基于深度学习的蛋白质结构预测、基因调控码学习、基因表达 预测、癌症分类预测、复杂疾病分类、多平台癌症数据综合分析等。
Guillen等人设计基于多层感知器的深度学习算法捕获基因表达特征进行癌症分类，表明了神经网络可以高效率地对不同的样本进行分类，在最后的预测结果 中实现了较高精度。Bhat等人通过深度生成学习检测癌症，使用对抗性特征学习 过程挖掘数据特征，然后使用常规分类器进行分类。最终试验通过指定适当的超 参数，在两个不同数据集上执行得相当好。Danaee等人使用堆叠去噪自动编码器 (SDAE)从高维基因表达谱中提取深度功能特征，通过分析SDAE连接矩阵确定 了一组高度互动的基因用于癌症生物标志物检测。Chira等人使用基因表达值随时 间推移的模式开发基于形状的聚类模型，并且进一步结合基因表达水平与输出值 之间的相关关系，考虑共同表达模式与测量输出的关系，以指导结果的生物学解 释。Singh等人提供层叠特征选择与堆叠稀疏自动编码器(SSAE)从数据中学习 高级特征，每层执行特征选择是一种启发式的，可以在每个阶段获得相关特征，并且在调整过程中减少计算量，该算法在GEMLeR数据库的36个数据集上进行了 测试，其中35个数据集的效果超越了GEMLeR基准测试结果。Liang等人提出了一 种多峰深度信念网络(DBN)的新学习模型，从多平台观测数据对癌症患者进行 聚类，并为个性化癌症治疗提供了有效指导。同时应用对比度发散(CD)学习算 法，以无监督的方式推断多模态DBN模型参数。Xie等人基于多层感知器和堆叠 去噪自动编码器(MLP-SAE)的深度学习回归模型预测变异基因型的基因表达， 其中堆叠去噪自动编码器用于训练回归模型以提取有效特征，并利用多层感知器 进行反向传播，同时通过添加dropout防止过拟合。Chen等人设计了一种深度学习 方法(D-GEX)，充分捕捉基因表达间的非线性相关关系，利用大约1000个标记 基因推断剩余的靶基因表达，旨在降低基因表达谱测定成本。
发明内容
为了克服已有基因调控网络建模及差异性分析方法的精确性较差的不足，本 发明提供一种精确性较好的基于RNN的基因调控网络构建与动态差异性分析方 法。
本发明解决其技术问题所采用的技术方案是：
一种基于RNN的基因调控网络构建与动态差异性分析方法，包括以下步骤：
第一步、基于deepRNN的基因动态调控网络构建
基因表达数据表示为其中表示第j个 样本的第i个基因在tk时刻的表达量，在t0时刻输入的样本基因表达数据矩阵表示 为输入序列表示为X＝{Vi|i∈[1,n]}；输出的调控 权值矩阵包括了节点度值、连边权重及调控方向，输出矩阵 W＝{ωij|i∈[1,n],j∈[1,n],s＝sa,t＝tk}实质是一个上三角矩阵，包含样本sa的信息 基因在tk时刻的调控关系，在某一时刻的前后短时间段内，在基因A的表达作用 于基因B时不受基因B的反作用；若|ωij|＜threshold则表示该基因对之间无连边； 否则当ωij＞0表示基因i对基因j上调，当ωij＜0表示下调；若i＝j则表示基因的 自调控；
第二步、基于亚型内动态调控网络的时序变化演化分析
定义C1亚型在t0时刻的有向加权图拓扑结构表示为即将信息基 因抽象成节点集合V，将基因间的关系或作用描绘成边集合E，其中表示 基因i以权重ωij对基因j的作用，为+(上调)或者-(下调)；使用真阳率、假 阳率、阳性预测率、准确率对调控网络的性能进行定量评价；
第三步、基于亚型间动态调控网络的网络差异演化分析
不同亚型网络的演化分析包括动力学分析、差异性分析和扰动分析，
所述动力学分析使用差分方程对离散的网络动力学行为进行分析，对于不同 亚型的动态调控网络，分析同一时间段关联基因对的节点度值、连边权重、表达 变化量相对比率；通过提取不同亚型网络的关联特征，并以此为基础构建多网络 协同演化模型；
所述差异性分析对相同时间窗口内不同亚型间的两个基因调控网络作基于节 点局部结构特征的减法运算，检测网络结构间存在的差异边，根据差异网络鉴别 关键枢纽基因，然后利用GO信息和KEGG通路功能富集性分析检验所发现基因集 的显著性，得到癌症亚型相关控制基因作为进一步生物实验的检验标记；
所述扰动控制分析中，关键枢纽基因节点在细胞生化过程中具有以下特征： 同功能中心，即该节点附近的基因属于某类功能的基因集；同驱动中心，即受到 该节点表达调控的同距离区间内的基因具有类似的生化功能，对于关键枢纽节点 的调控输入一个随机扰动ΥPer，对不同网络在同距离区间内的同功能基因集取交 集，得到亚型网络间的动态调控差异节点。
进一步，所述第一步中，基于deepRNN的基因动态调控网络构建包括以下步 骤：
1.1预处理，首先，提取亚型网络之间的信息基因，然后，将同一亚型内部 的样本按照百分比随机分为训练集80％，验证集10％，测试集10％，进一步，将同 一样本的基因表达按照时间序列展开作为输入向量： 表示第cx类亚型中具有mcx个样本，提 取其中的ninf个信息基因，按照不同时刻k的展开；
1.2激活函数与损失函数，采用ReLU非饱和激活函数，值域为[0,+∞)，公式 如下：
其中为激活函数；
deepRNN由一个输入层、一个或多个循环体隐藏层和一个输出层组成，所有 隐藏的层都有相同数量的隐藏单元，将上一时刻的状态与当前时刻的输入拼接成 一个大的向量作为循环体中神经网络的输入，得到第l层的第j个单元的信号输出 为：
其中H是隐藏单元个数，表示与第j个隐藏单元uj相连接的需要被训 练的权重和偏差；当误差从输出层反向传播回来时，使用LSTM的记忆元保存信 息，在训练时，将每个输出单元的均方误差作为损失函数，即：
其中M'表示训练样本个数，N表示每个训练样本基因个数，ωm(i,j)表示在t时刻 样本m中的基因gi对基因gj的作用效果，即连边权重，表示预测值；
1.3dropout方法，在训练过程中，对于每个训练样本的隐藏单元及其边缘将 会以概率为p被暂时丢弃；因此前向传播和后向传播将在一个特别“薄”的稀疏网 络上进行；对于deepRNN，只在同一时刻的不同层循环体之间使用dropout，即仅 在同一时刻t中，从h1到hlast的不同层循环体之间使用dropout；将在区间[0％，25％] 之间比较不同程度的正则化效果，寻找最优dropout比率；
1.4加速梯度优化和权重初始化，拟采用动量法进行加速优化，即通过在迭 代过程中累积损失函数的梯度方向来代替梯度进行参数更新，对于神经网络参数 Θ的损失函数L(·)，动量计算公式如下：
其中，μ∈[0,1]是动量系数，η是学习率；
隐藏层单位的权重使用均匀分布进行采样，定义如下:
其中ni，no分别表示隐藏单元的扇入扇出个数；
1.5输出，在循环体中的神经网络供给当前时刻的输出后，将会使用另外一 个全连接神经网络实现将当前时刻的状态转化为最终的输出。
再进一步，所述第二步中，网络的拓扑属性是描述网络本身及其内部节点或 边结构特征的测度，包括：
聚类系数，体现部分节点间存在的密集连接性质，在有向网络中，标准化的 聚类系数被定义为：
其中kout表示节点v的出度，n表示所有v所指向的节点彼此存在的边数，
介数表明一个节点在其他节点彼此连接中所起的作用，标准化至[0,1]区间的 计算公式如下：
其中σij是节点i到节点j的最短路径条数，σivj表示σij中通过节点v的路径条数；
紧密度是描述一个节点到网络中其他所有节点平均距离的指标，定量衡量节 点接近网络“中心”的程度，节点v的紧密度Cv计算公式如下：
其中dvj表示节点v到节点j的最短距离(路径中所经过边的权重之和最小)。紧 密度越小，节点越接近中心。
基于网络结构的拓扑属性变化在时间序列上对时间窗口Δt进行微分展开，得 到动态调控网络的时空演化测度Γ'(·)的计算公式如下：
其中Θ表示函数参数，ωCC、ωB、ωC分别为对应指标的影响权重；
通过分析动态网络在不同时刻的节点指标(CCv、Bv、Cv)，挖掘在不同时 间窗口内的关键调控基因节点，解释其在生命活动过程中扮演的重要性。
所述第三步中，所述动力学分析过程中，动力差异计算公式如下：
其中Θ表示节点度值、连边权重、表达变化量相对比率三项指标，表示Cx类亚型样本sx在时间窗口ta内的对应指标值，ωΘ表示对应于Θ的权重，ζ表示网 络间的全局耦合强度。
所述第三步中，所述差异性分析过程中，检测网络结构间存在的差异边的计 算公式如下：
其中表示Cx类亚型的样本sx在ta时间窗口内的基因i对基因j的调控权 重。
所述第三步中，所述扰动控制分析中，得到亚型网络间的动态调控差异节点， 表达式为：
其中表示Cx类亚型样本sa在区间dis(ω)内的基因集。
本发明的技术构思为：分析同一癌症亚型和不同癌症亚型的基因调控差异， 针对表达数据中癌症基因间的高度非线性相关性，基于深层循环神经网络(deep recurrentneural network，deepRNN)对不同癌症亚型在用药后的连续时序变化下 的基因表达数据构建调控网络，分析亚型间的表达差异性。
在基因表达数据的癌症关联基因特征提取后，完成癌症亚型的聚类分析，针 对不同的亚型聚类结果分别构建对应的基因调控网络分析其差异性。本项目提出 基于深层循环神经网络(deepRNN)的基因调控网络建模方法，利用深层循环神 经网络的时序处理特性，预测基因动态调控网络的节点度值与连边权重。其次， 纵向分析不同时间窗口中相同亚型调控网络的节点与连边变化，挖掘相关基因在 癌症演化过程中的调控功能，以及对病症发展的后续阶段进行预测。最终，横向 分析不同亚型间的调控网络差异，并对时间序列下的协同演化过程中的差异变化 进行生物学意义上的解释，为个性化临床治疗方案提供科学合理的指导。
本发明的有益效果主要表现在：精确性较好。
附图说明
图1是基因表达动态时序网络及差异性演化分析示意图。
图2是基于deepRNN的基因调控时序网络构建框图。
图3是亚型内部时序展开动态调控网络构建示意图。
图4是不同亚型间的基因调控网络渐变演化示意图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1～图4，一种基于RNN的基因调控网络构建与动态差异性分析方法， 分析同一癌症亚型和不同癌症亚型的基因调控差异，针对表达数据中癌症基因间 的高度非线性相关性，基于深层循环神经网络(deep recurrent neural network， deepRNN)对不同癌症亚型在用药后的连续时序变化下的基因表达数据构建调控 网络，分析亚型间的表达差异性；
如图1所示，首先，在t0时刻的癌症样本被聚类为C1、C2、C3三种亚型， 其余三个黑点表示奇异样本。对于C1类簇，基于deepRNN的调控网络构建如蓝色 虚线框中的t0时刻网络，显示根据A-H的8个信息基因构建调控网络，并通过真阳 率、假阳率、阳性预测率、准确率对网络性能进行定量评价；然后，在后续的数 据流输入后，网络的节点度值、连边权重值及节点位置发生迁移，得到了诸如t1、 t2…tl的动态演化调控网络，从而设计基于多层次动力系统模型的分析方法揭示 基因间调控过程中的逻辑关系；最后，进行不同亚型在网络间的横向分析，对于 不同亚型间的两个基因调控网络作基于节点局部结构特征的减法运算，检测网络 结构间存在的差异边，得到差异网络进而鉴别关键枢纽基因，同时利用GO信息和KEGG通路功能富集性分析检验所发现基因集的显著性，最终识别出癌症亚型相 关控制基因作为进一步生物实验的检验标记。
为了验证本项目提出的算法在处理真实癌症基因表达数据的实时性、有效性 和可靠性，并且获得算法的优化参数，本项目将先对常用基因表达数据库(如 GEO、TCGA、SMD、GXD、GENT等)中的癌症表达标准数据进行有针对性的 分类与分析，验证算法性能。
所述基因调控网络构建与动态差异性分析方法包括以下步骤：
第一步、基于deepRNN的基因动态调控网络构建
基因的表达具有时空性，是基因与外界环境相互作用的结果，会根据当前的 表达状况决定未来的表达，所以适合利用deepRNN的历史记忆效应，学习训练隐 藏层参数，最终以矩阵形式输出调控网络权值。如图2所示为按时序展开的循环 神经网络的构建及调控权重的训练过程。
基因表达数据表示为其中表示第j个 样本的第i个基因在tk时刻的表达量。如上图所示，在t0时刻输入的样本基因表达 数据矩阵表示为输入序列表示为 X＝{Vi|i∈[1,n]}；输出的调控权值矩阵包括了节点度值、连边权重及调控方向， 输出矩阵W＝{ωij|i∈[1,n],j∈[1,n],s＝sa,t＝tk}实质是一个上三角矩阵，包含样本 sa的信息基因在tk时刻的调控关系，当然是基于以下前提：即在某一时刻的前后 短时间段内，在基因A的表达作用于基因B时不受基因B的反作用。若 |ωij|＜threshold则表示该基因对之间无连边；否则当ωij＞0表示基因i对基因j上 调，当ωij＜0表示下调；若i＝j则表示基因的自调控。下面介绍如何利用deepRNN对基因表达数据进行基因动态调控网络的构建：
1.1预处理。首先，提取亚型网络之间的信息基因，一方面是因为样本的过 长输入时间序列间隔会导致优化时的“梯度弥散”问题；另一方面也是因为在某一 调控过程中的无关基因相当于噪声，使用强有力的控制基因能够更好的挖掘调控 关系。然后，将同一亚型内部的样本按照百分比随机分为训练集80％，验证集10％， 测试集10％。进一步，将同一样本的基因表达按照时间序列展开作为输入向量： 表示第cx类亚型中具有mcx个样本，提 取其中的ninf个信息基因，按照不同时刻k的展开。
1.2激活函数与损失函数。激活函数作为非线性处理单元(如Sigmoid、tanh 函数)，实现的功能是将来自前一层的输入线性组合结果动态范围压缩到特定值 域。为了缓解深度神经网络的“梯度弥散”问题，加快训练收敛速度，拟采用ReLU 这类非饱和激活函数(值域为[0,+∞))，公式如下：
其中为激活函数。
deepRNN由一个输入层、一个或多个循环体隐藏层和一个输出层组成。所有 隐藏的层都有相同数量的隐藏单元，将上一时刻的状态与当前时刻的输入拼接成 一个大的向量作为循环体中神经网络的输入，得到第l层的第j个单元的信号输出 为：
其中H是隐藏单元个数，表示与第j个隐藏单元uj相连接的需要被训 练的权重和偏差。为了处理长期依赖问题，隐藏层单元使用三门结构(输入门、 遗忘门、输出门)的长短时记忆单元(Long Short Term Memory，LSTM)单元， 当误差从输出层反向传播回来时，可以使用LSTM的记忆元保存信息。在训练时， 将每个输出单元的均方误差作为损失函数，即：
其中M'表示训练样本个数，N表示每个训练样本基因个数，ωm(i,j)表示在t时刻 样本m中的基因gi对基因gj的作用效果，即连边权重，表示预测值。
1.3dropout方法。dropout是对神经网络进行模型平均和正则化的技术。在训 练过程中，对于每个训练样本的隐藏单元及其边缘将会以概率为p被暂时丢弃。 因此前向传播和后向传播将在一个特别“薄”的稀疏网络上进行。对于deepRNN， 一般只在同一时刻的不同层循环体之间使用dropout，即仅在同一时刻t中，从h1到 hlast的不同层循环体之间使用dropout，这样能够使得网络更加健壮。参考相关文 献，将在区间[0％，25％]之间比较不同程度的正则化效果，寻找最优dropout比率。
1.4加速梯度优化和权重初始化。拟采用动量法进行加速优化，即通过在迭 代过程中累积损失函数的梯度方向来代替梯度进行参数更新。对于神经网络参数 Θ的损失函数L(·)，动量计算公式如下：
其中，μ∈[0,1]是动量系数，η是学习率，在训练过程中随着错误率变化而不断减小，使用动量法在训练深度神经网络时能够提高收敛速度。深度网络的权重使用 归一法进行初始化，旨在稳定训练过程中的激活和反向传播梯度的差异。隐藏层 单位的权重使用均匀分布进行采样，定义如下:
其中ni，no分别表示隐藏单元的扇入扇出个数。
1.5输出。在循环体中的神经网络供给当前时刻的输出后，将会使用另外一 个全连接神经网络实现将当前时刻的状态转化为最终的输出。
第二步、基于亚型内动态调控网络的时序变化演化分析
通过基因表达数据构建动态基因调控网络来体现真实的动态调控过程能够更 精准地反映调控机理，理解基因之间的相互作用机制。如图3所示为某一亚型样本 内部的信息基因以时间序列展开后的动态调控网络构建示意图。
定义C1亚型在t0时刻的有向加权图拓扑结构表示为即将信息基 因抽象成节点集合V，将基因间的关系或作用描绘成边集合E，其中表示 基因i以权重ωij对基因j的作用，:可以为+(上调)或者-(下调)。可以使用 真阳率、假阳率、阳性预测率、准确率对调控网络的性能进行定量评价。
网络的拓扑属性是描述网络本身及其内部节点或边结构特征的测度。主要包 括以下几项：
聚类系数。聚类系数体现了部分节点间存在的密集连接性质，在有向网络中， 标准化的聚类系数被定义为：
其中kout表示节点v的出度，n表示所有v所指向的节点彼此存在的边数。
介数。介数表明了一个节点在其他节点彼此连接中所起的作用，标准化至 [0,1]区间的计算公式如下：
其中σij是节点i到节点j的最短路径条数，σivj表示σij中通过节点v的路径条数。介数越高，意味着节点在保持网络连接紧密性中越重要。
紧密度。紧密度是描述一个节点到网络中其他所有节点平均距离的指标，可 以定量衡量节点接近网络“中心”的程度。节点v的紧密度Cv计算公式如下：
其中dvj表示节点v到节点j的最短距离(路径中所经过边的权重之和最小)。紧 密度越小，节点越接近中心。
为了描述调控网络的动力学性质，基于网络结构的拓扑属性变化在时间序列 上对时间窗口Δt进行微分展开，得到动态调控网络的时空演化测度Γ'(·)的计算公 式如下：
其中Θ表示函数参数，ωCC、ωB、ωC分别为对应指标的影响权重。
通过分析动态网络在不同时刻的节点指标(CCv、Bv、Cv)，挖掘在不同时 间窗口内的关键调控基因节点，解释其在生命活动过程中扮演的重要性。
第三步、基于亚型间动态调控网络的网络差异演化分析
亚型间的网络分析是为了充分挖掘历史表达数据的时空特性、调控的变化规 律、以及调控网络中节点和连边的迁移演化，从而提高网络建模算法的准确度和 可靠性，并对基因表达表达变化和网络动态演化进行预测。如图4所示为不同亚型 间的基因调控网络渐变演化示意图，其中差异网络是由不同亚型间的两个基因调 控网络作基于节点局部结构特征的减法运算得到的。
不同亚型网络的演化分析包括动力学分析、差异性分析和扰动分析。
动力学分析。由于基因表达数据采样的时间间隔较长，使用差分方程对离散 的网络动力学行为进行分析。对于不同亚型的动态调控网络，分析同一时间段关 联基因对的节点度值、连边权重、表达变化量相对比率，动力差异计算公式如下：
其中Θ表示节点度值、连边权重、表达变化量相对比率三项指标，表示Cx类亚型样本sx在时间窗口ta内的对应指标值，ωΘ表示对应于Θ的权重，ζ表示网 络间的全局耦合强度。通过提取不同亚型网络的关联特征，并以此为基础构建多 网络协同演化模型，该模型将为后续网络差异的演化分析及精准亚型个体化治疗 提供理论基础和实践依据。
差异性分析。对相同时间窗口内不同亚型间的两个基因调控网络作基于节点 局部结构特征的减法运算，检测网络结构间存在的差异边，计算公式如下：
其中表示Cx类亚型的样本sx在ta时间窗口内的基因i对基因j的调控权 重。根据差异网络鉴别关键枢纽基因，然后利用GO信息和KEGG通路功能富集性 分析检验所发现基因集的显著性，得到癌症亚型相关控制基因作为进一步生物实 验的检验标记。
扰动控制分析。关键枢纽基因节点在细胞生化过程中具有以下特征：同功能 中心，即该节点附近的基因属于某类功能的基因集；同驱动中心，即受到该节点 表达调控的同距离区间内的基因具有类似的生化功能。对于关键枢纽节点的调控 输入一个随机扰动ΥPer，对不同网络在同距离区间内的同功能基因集取交集，得 到亚型网络间的动态调控差异节点。具体表达式为：
其中表示Cx类亚型样本sa在区间dis(ω)内的基因集。
Claims (6)
1.一种基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：包括以下步骤：
第一步、基于deepRNN的基因动态调控网络构建
基因表达数据表示为其中表示第j个样本的第i个基因在tk时刻的表达量，在t0时刻输入的样本基因表达数据矩阵表示为输入序列表示为X＝{Vi|i∈[1,n]}；输出的调控权值矩阵包括了节点度值、连边权重及调控方向，输出矩阵W＝{ωij|i∈[1,n],j∈[1,n],s＝sa,t＝tk}实质是一个上三角矩阵，包含样本sa的信息基因在tk时刻的调控关系，在某一时刻的前后短时间段内，在基因A的表达作用于基因B时不受基因B的反作用；若|ωij|＜threshold则表示该基因对之间无连边；否则当ωij＞0表示基因i对基因j上调，当ωij＜0表示下调；若i＝j则表示基因的自调控；
第二步、基于亚型内动态调控网络的时序变化演化分析
定义C1亚型在t0时刻的有向加权图拓扑结构表示为即将信息基因抽象成节点集合V，将基因间的关系或作用描绘成边集合E，其中表示基因i以权重ωij对基因j的作用，为+(上调)或者-(下调)；使用真阳率、假阳率、阳性预测率、准确率对调控网络的性能进行定量评价；
第三步、基于亚型间动态调控网络的网络差异演化分析
不同亚型网络的演化分析包括动力学分析、差异性分析和扰动分析，
所述动力学分析使用差分方程对离散的网络动力学行为进行分析，对于不同亚型的动态调控网络，分析同一时间段关联基因对的节点度值、连边权重、表达变化量相对比率；通过提取不同亚型网络的关联特征，并以此为基础构建多网络协同演化模型；
所述差异性分析对相同时间窗口内不同亚型间的两个基因调控网络作基于节点局部结构特征的减法运算，检测网络结构间存在的差异边，根据差异网络鉴别关键枢纽基因，然后利用GO信息和KEGG通路功能富集性分析检验所发现基因集的显著性，得到癌症亚型相关控制基因作为进一步生物实验的检验标记；
所述扰动控制分析中，关键枢纽基因节点在细胞生化过程中具有以下特征：同功能中心，即该节点附近的基因属于某类功能的基因集；同驱动中心，即受到该节点表达调控的同距离区间内的基因具有类似的生化功能，对于关键枢纽节点的调控输入一个随机扰动ΥPer，对不同网络在同距离区间内的同功能基因集取交集，得到亚型网络间的动态调控差异节点。
2.如权利要求1所述的基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：所述第一步中，基于deepRNN的基因动态调控网络构建包括以下步骤：
1.1预处理，首先，提取亚型网络之间的信息基因，然后，将同一亚型内部的样本按照百分比随机分为训练集80％，验证集10％，测试集10％，进一步，将同一样本的基因表达按照时间序列展开作为输入向量：
表示第cx类亚型中具有mcx个样本，提取其中的ninf个信息基因，按照不同时刻k的展开；
1.2激活函数与损失函数，采用ReLU非饱和激活函数，值域为[0,+∞)，公式如下：
其中为激活函数；
deepRNN由一个输入层、一个或多个循环体隐藏层和一个输出层组成，所有隐藏的层都有相同数量的隐藏单元，将上一时刻的状态与当前时刻的输入拼接成一个大的向量作为循环体中神经网络的输入，得到第l层的第j个单元的信号输出为：
其中H是隐藏单元个数，表示与第j个隐藏单元uj相连接的需要被训练的权重和偏差；当误差从输出层反向传播回来时，使用LSTM的记忆元保存信息，在训练时，将每个输出单元的均方误差作为损失函数，即：
<mrow>
<mi>L</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mi>i</mi>
</mrow>
<mi>N</mi>
</munderover>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>M</mi>
<mo>&prime;</mo>
</msup>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>M</mi>
<mo>&prime;</mo>
</msup>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>&omega;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中M'表示训练样本个数，N表示每个训练样本基因个数，ωm(i,j)表示在t时刻样本m中的基因gi对基因gj的作用效果，即连边权重，表示预测值；
1.3dropout方法，在训练过程中，对于每个训练样本的隐藏单元及其边缘将会以概率为p被暂时丢弃；因此前向传播和后向传播将在一个特别“薄”的稀疏网络上进行；对于deepRNN，只在同一时刻的不同层循环体之间使用dropout，即仅在同一时刻t中，从h1到hlast的不同层循环体之间使用dropout；将在区间[0％，25％]之间比较不同程度的正则化效果，寻找最优dropout比率；
1.4加速梯度优化和权重初始化，拟采用动量法进行加速优化，即通过在迭代过程中累积损失函数的梯度方向来代替梯度进行参数更新，对于神经网络参数Θ的损失函数L(·)，动量计算公式如下：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>&mu;V</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mi>&eta;</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&Theta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Theta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>&Theta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，μ∈[0,1]是动量系数，η是学习率；
隐藏层单位的权重使用均匀分布进行采样，定义如下:
<mrow>
<mi>W</mi>
<mo>~</mo>
<mi>U</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>6</mn>
</msqrt>
<msqrt>
<mrow>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>o</mi>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>,</mo>
<mfrac>
<msqrt>
<mn>6</mn>
</msqrt>
<msqrt>
<mrow>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>o</mi>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中ni，no分别表示隐藏单元的扇入扇出个数；
1.5输出，在循环体中的神经网络供给当前时刻的输出后，将会使用另外一个全连接神经网络实现将当前时刻的状态转化为最终的输出。
3.如权利要求2所述的基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：所述第二步中，网络的拓扑属性是描述网络本身及其内部节点或边结构特征的测度，包括：
聚类系数，体现部分节点间存在的密集连接性质，在有向网络中，标准化的聚类系数被定义为：
<mrow>
<msub>
<mi>CC</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>n</mi>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中kout表示节点v的出度，n表示所有v所指向的节点彼此存在的边数，
介数表明一个节点在其他节点彼此连接中所起的作用，标准化至[0,1]区间的计算公式如下：
<mrow>
<msub>
<mi>B</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>v</mi>
<mo>&Element;</mo>
<mi>V</mi>
</mrow>
</munder>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中σij是节点i到节点j的最短路径条数，σivj表示σij中通过节点v的路径条数；
紧密度是描述一个节点到网络中其他所有节点平均距离的指标，定量衡量节点接近网络“中心”的程度，节点v的紧密度Cv计算公式如下：
<mrow>
<msub>
<mi>C</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>v</mi>
<mo>&Element;</mo>
<mi>V</mi>
</mrow>
</munder>
<msub>
<mi>d</mi>
<mrow>
<mi>v</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中dvj表示节点v到节点j的最短距离(路径中所经过边的权重之和最小)。紧密度越小，节点越接近中心。
基于网络结构的拓扑属性变化在时间序列上对时间窗口Δt进行微分展开，得到动态调控网络的时空演化测度Γ'(·)的计算公式如下：
<mrow>
<msup>
<mi>&Gamma;</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>C</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&Delta;CC</mi>
<mi>v</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>B</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&Delta;B</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>C</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&Delta;C</mi>
<mi>v</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Θ表示函数参数，ωCC、ωB、ωC分别为对应指标的影响权重；
通过分析动态网络在不同时刻的节点指标(CCv、Bv、Cv)，挖掘在不同时间窗口内的关键调控基因节点，解释其在生命活动过程中扮演的重要性。
4.如权利要求3所述的基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：所述第三步中，所述动力学分析过程中，动力差异计算公式如下：
<mrow>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>d</mi>
<mi>y</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&zeta;</mi>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mn>3</mn>
</munderover>
<msub>
<mi>&omega;</mi>
<mi>&Theta;</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&Theta;</mi>
<mrow>
<mo><</mo>
<mi>i</mi>
<mi>j</mi>
<mo>,</mo>
<msub>
<mi>s</mi>
<mi>a</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>a</mi>
</msub>
<mo>></mo>
</mrow>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Theta;</mi>
<mrow>
<mo><</mo>
<mi>i</mi>
<mi>j</mi>
<mo>,</mo>
<msub>
<mi>s</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>b</mi>
</msub>
<mo>></mo>
</mrow>
<msub>
<mi>C</mi>
<mi>b</mi>
</msub>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Θ表示节点度值、连边权重、表达变化量相对比率三项指标，表示Cx类亚型样本sx在时间窗口ta内的对应指标值，ωΘ表示对应于Θ的权重，ζ表示网络间的全局耦合强度。
5.如权利要求4所述的基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：所述第三步中，所述差异性分析过程中，检测网络结构间存在的差异边的计算公式如下：
其中表示Cx类亚型的样本sx在ta时间窗口内的基因i对基因j的调控权重。
6.如权利要求4所述的基于RNN的基因调控网络构建与动态差异性分析方法，其特征在于：所述第三步中，所述扰动控制分析中，得到亚型网络间的动态调控差异节点，表达式为：
<mrow>
<msub>
<mi>gene</mi>
<mrow>
<mi>D</mi>
<mi>I</mi>
<mi>F</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>{</mo>
<msubsup>
<mi>gene</mi>
<mrow>
<mo><</mo>
<msub>
<mi>s</mi>
<mi>a</mi>
</msub>
<mo>,</mo>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
</mrow>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
</msubsup>
<mo>}</mo>
<mo>&cap;</mo>
<mo>{</mo>
<msubsup>
<mi>gene</mi>
<mrow>
<mo><</mo>
<msub>
<mi>s</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
<mi>d</mi>
<mi>i</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
</mrow>
<msub>
<mi>C</mi>
<mi>b</mi>
</msub>
</msubsup>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
3
其中表示Cx类亚型样本sa在区间dis(ω)内的基因集。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710355357.5A CN107220525B (zh) | 2017-05-19 | 2017-05-19 | 基于rnn的基因调控网络构建与动态差异性分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710355357.5A CN107220525B (zh) | 2017-05-19 | 2017-05-19 | 基于rnn的基因调控网络构建与动态差异性分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107220525A true CN107220525A (zh) | 2017-09-29 |
| CN107220525B CN107220525B (zh) | 2021-06-18 |
Family
ID=59944302
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710355357.5A Active CN107220525B (zh) | 2017-05-19 | 2017-05-19 | 基于rnn的基因调控网络构建与动态差异性分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107220525B (zh) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107992945A (zh) * | 2017-12-14 | 2018-05-04 | 浙江工业大学 | 基于深度学习和进化计算的特征基因选择方法 |
| CN108090326A (zh) * | 2018-02-09 | 2018-05-29 | 国家卫生计生委科学技术研究所 | 单细胞网络调控关系的构建方法 |
| CN108549832A (zh) * | 2018-01-21 | 2018-09-18 | 西安电子科技大学 | 基于全连接神经网络的低截获雷达信号分类方法 |
| CN108846261A (zh) * | 2018-04-26 | 2018-11-20 | 浙江工业大学 | 基于可视图算法的基因表达时序数据分类方法 |
| CN109215735A (zh) * | 2018-09-21 | 2019-01-15 | 西南民族大学 | 一种构建基因调控网络的方法 |
| CN109243523A (zh) * | 2018-08-24 | 2019-01-18 | 东北大学 | 一种基于乳腺癌疾病的调控网络构建及分析方法 |
| CN109360607A (zh) * | 2018-10-16 | 2019-02-19 | 武汉大学 | 一种动态基因调控网的网络演化分析方法及装置 |
| CN110310707A (zh) * | 2019-07-08 | 2019-10-08 | 西北工业大学 | 基于时序数据的表型本体构建方法 |
| CN110309528A (zh) * | 2019-04-15 | 2019-10-08 | 南京航空航天大学 | 一种基于机器学习的雷达方案设计方法 |
| CN110322930A (zh) * | 2019-06-06 | 2019-10-11 | 大连理工大学 | 基于水平关系的代谢组学网络标志物识别方法 |
| CN110941928A (zh) * | 2019-11-26 | 2020-03-31 | 哈尔滨理工大学 | 一种基于dropout-SAE和Bi-LSTM的滚动轴承剩余寿命预测方法 |
| CN111192631A (zh) * | 2020-01-02 | 2020-05-22 | 中国科学院计算技术研究所 | 用于构建用于预测蛋白质-rna相互作用结合位点模型的方法和系统 |
| CN111508555A (zh) * | 2020-04-15 | 2020-08-07 | 淮南师范学院 | 一组用于度量生物分子集的网络拓扑学特性的方法 |
| CN111816246A (zh) * | 2020-05-27 | 2020-10-23 | 上海大学 | 从差异网络中识别驱动基因的方法 |
| CN112037850A (zh) * | 2020-09-04 | 2020-12-04 | 中国科学院重庆绿色智能技术研究院 | 基于动量加速的缺失蛋白质间相互作用预测装置和方法 |
| CN112885404A (zh) * | 2021-03-29 | 2021-06-01 | 哈尔滨理工大学 | 一种多层布尔网络的模型辨识方法及系统 |
| CN113032776A (zh) * | 2021-02-08 | 2021-06-25 | 浙江工业大学 | 面向特征嵌入中毒攻击的检测方法、装置和系统 |
| EP4202941A1 (en) * | 2021-12-21 | 2023-06-28 | Dassault Systèmes | Inferrence of a gene expression profile |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093119A (zh) * | 2013-01-24 | 2013-05-08 | 南京大学 | 一种利用网络结构信息识别显著生物通路的方法 |
-
2017
- 2017-05-19 CN CN201710355357.5A patent/CN107220525B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093119A (zh) * | 2013-01-24 | 2013-05-08 | 南京大学 | 一种利用网络结构信息识别显著生物通路的方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘杰辉等: "基于ReLU激活函数的轧制力神经网络预报模型", 《锻压技术》 * |
| 张文兵: "基因调控网络的随机动力学分析", 《中国博士学位论文全文数据库基础科学辑》 * |
| 杨博: "基因调控网络构建及在癌症基因预测中的应用", 《中国博士学位论文全文数据库基础科学辑》 * |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107992945B (zh) * | 2017-12-14 | 2020-03-24 | 浙江工业大学 | 基于深度学习和进化计算的特征基因选择方法 |
| CN107992945A (zh) * | 2017-12-14 | 2018-05-04 | 浙江工业大学 | 基于深度学习和进化计算的特征基因选择方法 |
| CN108549832A (zh) * | 2018-01-21 | 2018-09-18 | 西安电子科技大学 | 基于全连接神经网络的低截获雷达信号分类方法 |
| CN108549832B (zh) * | 2018-01-21 | 2021-11-30 | 西安电子科技大学 | 基于全连接神经网络的低截获雷达信号分类方法 |
| CN108090326A (zh) * | 2018-02-09 | 2018-05-29 | 国家卫生计生委科学技术研究所 | 单细胞网络调控关系的构建方法 |
| CN108846261A (zh) * | 2018-04-26 | 2018-11-20 | 浙江工业大学 | 基于可视图算法的基因表达时序数据分类方法 |
| CN108846261B (zh) * | 2018-04-26 | 2020-05-29 | 浙江工业大学 | 基于可视图算法的基因表达时序数据分类方法 |
| CN109243523A (zh) * | 2018-08-24 | 2019-01-18 | 东北大学 | 一种基于乳腺癌疾病的调控网络构建及分析方法 |
| CN109243523B (zh) * | 2018-08-24 | 2021-06-11 | 东北大学 | 一种基于乳腺癌疾病的调控网络构建及分析方法 |
| CN109215735A (zh) * | 2018-09-21 | 2019-01-15 | 西南民族大学 | 一种构建基因调控网络的方法 |
| CN109360607A (zh) * | 2018-10-16 | 2019-02-19 | 武汉大学 | 一种动态基因调控网的网络演化分析方法及装置 |
| CN109360607B (zh) * | 2018-10-16 | 2021-05-04 | 武汉大学 | 一种动态基因调控网的网络演化分析方法及装置 |
| CN110309528A (zh) * | 2019-04-15 | 2019-10-08 | 南京航空航天大学 | 一种基于机器学习的雷达方案设计方法 |
| CN110322930B (zh) * | 2019-06-06 | 2021-12-03 | 大连理工大学 | 基于水平关系的代谢组学网络标志物识别方法 |
| CN110322930A (zh) * | 2019-06-06 | 2019-10-11 | 大连理工大学 | 基于水平关系的代谢组学网络标志物识别方法 |
| CN110310707A (zh) * | 2019-07-08 | 2019-10-08 | 西北工业大学 | 基于时序数据的表型本体构建方法 |
| CN110941928A (zh) * | 2019-11-26 | 2020-03-31 | 哈尔滨理工大学 | 一种基于dropout-SAE和Bi-LSTM的滚动轴承剩余寿命预测方法 |
| CN111192631A (zh) * | 2020-01-02 | 2020-05-22 | 中国科学院计算技术研究所 | 用于构建用于预测蛋白质-rna相互作用结合位点模型的方法和系统 |
| CN111192631B (zh) * | 2020-01-02 | 2023-07-21 | 中国科学院计算技术研究所 | 用于构建用于预测蛋白质-rna相互作用结合位点模型的方法和系统 |
| CN111508555A (zh) * | 2020-04-15 | 2020-08-07 | 淮南师范学院 | 一组用于度量生物分子集的网络拓扑学特性的方法 |
| CN111816246A (zh) * | 2020-05-27 | 2020-10-23 | 上海大学 | 从差异网络中识别驱动基因的方法 |
| CN111816246B (zh) * | 2020-05-27 | 2023-01-10 | 上海大学 | 从差异网络中识别驱动基因的方法 |
| CN112037850A (zh) * | 2020-09-04 | 2020-12-04 | 中国科学院重庆绿色智能技术研究院 | 基于动量加速的缺失蛋白质间相互作用预测装置和方法 |
| CN113032776A (zh) * | 2021-02-08 | 2021-06-25 | 浙江工业大学 | 面向特征嵌入中毒攻击的检测方法、装置和系统 |
| CN113032776B (zh) * | 2021-02-08 | 2022-08-05 | 浙江工业大学 | 面向特征嵌入中毒攻击的检测方法、装置和系统 |
| CN112885404A (zh) * | 2021-03-29 | 2021-06-01 | 哈尔滨理工大学 | 一种多层布尔网络的模型辨识方法及系统 |
| CN112885404B (zh) * | 2021-03-29 | 2023-11-21 | 哈尔滨理工大学 | 一种多层布尔网络的模型辨识方法及系统 |
| EP4202941A1 (en) * | 2021-12-21 | 2023-06-28 | Dassault Systèmes | Inferrence of a gene expression profile |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107220525B (zh) | 2021-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107220525A (zh) | 基于rnn的基因调控网络构建与动态差异性分析方法 | |
| CN107885853A (zh) | 一种基于深度学习的组合式文本分类方法 | |
| Dong et al. | Automatic age estimation based on deep learning algorithm | |
| CN112328767B (zh) | 基于bert模型和比较聚合框架的问答匹配方法 | |
| CN105975916B (zh) | 基于多输出卷积神经网络和有序回归的年龄估计方法 | |
| CN107153642A (zh) | 一种基于神经网络识别文本评论情感倾向的分析方法 | |
| CN103714261B (zh) | 二阶段混合模型的智能辅助医疗决策支持方法 | |
| CN108171209A (zh) | 一种基于卷积神经网络进行度量学习的人脸年龄估计方法 | |
| CN106919951A (zh) | 一种基于点击与视觉融合的弱监督双线性深度学习方法 | |
| CN109934261A (zh) | 一种知识驱动参数传播模型及其少样本学习方法 | |
| CN110321361A (zh) | 基于改进的lstm神经网络模型的试题推荐判定方法 | |
| CN111275113A (zh) | 基于代价敏感混合网络的偏斜类时间序列异常检测方法 | |
| Li et al. | A new method of identification of complex lithologies and reservoirs: task-driven data mining | |
| CN110674483B (zh) | 一种基于多模态信息的身份识别方法 | |
| CN105955959A (zh) | 一种情感分类方法及系统 | |
| CN111460157A (zh) | 用于多领域文本分类的循环卷积多任务学习方法 | |
| CN110021341A (zh) | 一种基于异构网络的gpcr药物和靶向通路的预测方法 | |
| CN109783629A (zh) | 一种融合全局事件关系信息的微博客事件谣言检测方法 | |
| CN107944049A (zh) | 一种基于深度学习的电影推荐方法 | |
| CN115269860A (zh) | 一种飞机维修领域的知识图谱构建方法 | |
| CN113868374B (zh) | 基于多头注意力机制的图卷积网络生物医学信息提取方法 | |
| CN110413781A (zh) | 一种基于LD和ANN-SoftMax Regressor的中医药短文本多分类方法 | |
| Liu et al. | Comparison and evaluation of activation functions in term of gradient instability in deep neural networks | |
| CN111652444A (zh) | 一种基于K-means和LSTM的日游客量预测方法 | |
| Puri et al. | Review on missing value imputation techniques in data mining |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |