CN109243523B - 一种基于乳腺癌疾病的调控网络构建及分析方法 - Google Patents

一种基于乳腺癌疾病的调控网络构建及分析方法 Download PDF

Info

Publication number
CN109243523B
CN109243523B CN201810972336.2A CN201810972336A CN109243523B CN 109243523 B CN109243523 B CN 109243523B CN 201810972336 A CN201810972336 A CN 201810972336A CN 109243523 B CN109243523 B CN 109243523B
Authority
CN
China
Prior art keywords
genes
centrality
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.)
Active
Application number
CN201810972336.2A
Other languages
English (en)
Other versions
CN109243523A (zh
Inventor
王之琼
曲璐渲
郭上慧
霍岳阳
高笑宇
钱唯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN201810972336.2A priority Critical patent/CN109243523B/zh
Publication of CN109243523A publication Critical patent/CN109243523A/zh
Application granted granted Critical
Publication of CN109243523B publication Critical patent/CN109243523B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出一种基于乳腺癌疾病的基因调控网络构建及分析方法，流程包括：乳腺癌相关基因初步筛选；基因调控网络的构建；根据基因调控网络结果，进行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性；取网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据作为为筛选出的相关基因；在人类基因中有效地筛选出了乳腺癌相关基因，从而建立与乳腺癌疾病相关的基因调控网络，并通过节点中心性分析得到重要基因，可以促进从基因学的角度对乳腺癌疾病的研究，为找到干预乳腺癌疾病发生的有效途径奠定基础。

Description

一种基于乳腺癌疾病的调控网络构建及分析方法
技术领域
本发明属于医学信息学领域，涉及一种基于乳腺癌疾病的基因调控网络构建及分析方法。
背景技术
遗传表达确定了细胞和组织的表型和发育状态。人类常见的疾病，如恶性肿瘤和神经退行性疾病都可以追溯到异常基因表达。通常，基因不是单独表达的，一个基因表达可以影响其他基因的表达，同时也受其他基因表达的影响。基因之间的相互作用和相互关联的调节关系形成了基因调控网络。基因调控网络的建模和分析为乳腺癌疾病提供了遗传视角的治疗和诊断。因此，构建与疾病相关的基因调控网络对探索人体机制和疾病治疗具有重要意义。
现在的网络建模方法有效地解决了基因调控网络构建的问题，但也存在一些缺点。例如，贝叶斯网络模型可以准确地建立监管网络，但仅适用于构建小规模网络。互信息可以支持大规模网络建模，但不能描述基因之间的调节关系。多模型组合可以改善上述两个问题，但仍不能有效满足基因组规模调节网络建模的需求。
发明内容
针对现有技术的不足，本发明的目的是提出基于乳腺癌疾病的基因调控网络构建及分析方法，利用基因对之间的互信息值进行逐轮筛选，在全基因组数据中筛选出与乳腺癌疾病相关的200个基因进行基因调控网络构建并进行网络分析，为乳腺癌疾病的研究人员提供有价值的参考。
本发明提出一种基于乳腺癌疾病的基因调控网络构建及分析方法，包括以下步骤：
步骤1：乳腺癌相关基因初步筛选：将乳腺癌基因分为基准基因、相关基因与筛选出的相关基因，基准基因指现有技术中已经认定与乳腺癌疾病有关系的基因，相关基因指除现有技术中以外的所有基因，筛选出的相关基因指用本申请的方法筛选出与乳腺癌疾病关系紧密的基因，包括步骤1.1～步骤1.4：
步骤1.1：将基因表达数据分为两个矩阵，X矩阵为乳腺癌基准基因的表达数据，Y矩阵是待筛选的其余乳腺癌相关基因的表达数据，在X矩阵与Y矩阵中，每一行为一个基因，每一列为一个样本；
步骤1.2：矩阵Y中的向量Yj依次与X中的向量Xi计算每个基因之间的互信息值，Yj和Xi代表第j个和第i个基因的数据；
步骤1.3：人工选择阈值Isn，若互信息值大于等于Isn，则向量Yj对应的基因添加到X矩阵中，并在矩阵Y中删除，阈值Isn为动态选择的，根据每一轮的互信息结果，对每轮的互信息的值大小进行排序，Isn即为当前轮中第sn个阈值，筛选出来互信息值最大的前sn个基因；
步骤1.4：重复步骤1.2至1.3，直到矩阵X中的基因数量达到预设定的值M个；
步骤2：基因调控网络的构建：利用步骤1中获得基因的表达数据矩阵X，通过BNFinder2工具进行调控网络构建，设置调控节点集合的数量，评分函数为BDE评分函数，得到一个概率网络，网络节点代表基因，有向边代表基因之间的调控关系；具体包括步骤2.1～步骤2.3：
步骤2.1：数据处理：根据X矩阵，确定基因个数及样本数量，并为每个基因设定潜在调控基因；
步骤2.2：设定基因调控网络结构的结构学习及参数学习，采用基于评分搜索的结构学习方法，遍历所有可能的结构，然后用BDE评分函数衡量各个结构，进而找出最好的结构；具体包括步骤2.2.1～步骤2.2.4：
步骤2.2.1：对于某个基因，首先获取该基因的潜在调控基因及其表达数据；
步骤2.2.2：为基因的每个潜在调控节点赋予权重；
步骤2.2.3：获取调控基因为空集的基因的BDE分数；
步骤2.2.4：遍历基因的潜在调控基因的集合，每次遍历结果即为该基因的一种调控结构，计算该结构的BDE分数，存入到一个堆栈中，遍历结束后，堆栈中留下的为评分最高的对应的调控基因的集合，堆栈堆的大小即为设定的调控节点集合的数量；
步骤2.3：根据网络的学习结果整合调控关系，得到基因调控网络结果。；
步骤3：根据基因调控网络结果，进行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性，将BNFinder2得到的基因调控网络，通过Python的扩展包networkx计算各个基因的度中心性，接近中心性，中介中心性和特征向量中心性；
步骤4：取网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据作为为筛选出的相关基因。
所述网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据，N≤10，且矩阵X中的基因数量M>N。
所述步骤1.3中人工选择阈值Isn，其中sn等于需要筛选出来的基因数量除以需要的轮数。
所述基准基因包括基因名称如下：AKT1、AKT2、APC、ARID1A、ARID1B、ARID2、ASXL1、BAP1、BRCA1、BRCA2、CASP8、CDH1、CDKN1B、CDKN2A、CCND1、MDM2、99Mb、ERBB2、ZNF217、FGFR1/ZNF703、MYC、GATA3、KRAS、MAP2K4、MAP3K1、MAP3K13、MIL2、MIL3、NCOR1、NF1、PIK3CA、PTEN、RB1、SETD2、SF3B1、SMAD4、SMARCD1、STK11、TBX3和TP53。
所述互信息的计算方法，如下公式：
Figure BDA0001776525790000031
其中，I(Xi,Yj)为向量Xi和Yj的互信息，|C(Xi)|是向量Xi的协方差矩阵的行列式，|C(Yj)|是向量Yj的协方差矩阵的行列式，|C(Xi,Yj)|是向量Xi和Yj的协方差矩阵的行列式。
所述节点度中心性(Degree Centrality)：在一个网络图G＝(V,E)中，节点v的度dv指的是与v相连的E中边的数量，公式为：
Figure BDA0001776525790000032
其中，CB(v)为节点度中心性，m为该节点的度，n为基因的数量；
所述接近中心性(Closeness Centrality)：考量每个结点到其它结点的最短路的平均长度。距离其它结点越近，那么它的中心度越高，公式为：
Figure BDA0001776525790000033
其中，Ccl(v)为接近中心性，n为基因的数量，d(v,u)为节点v和u之间的最短路径。
所述中介中心性(Betweenness Centrality)：以经过某个节点的最短路径数目来刻画节点重要性的指标。一个结点充当“中介”的次数越高，它的中介中心度就越大，公式为：
Figure BDA0001776525790000034
其中，CB(v)为中介中心性，σ(s,t|v)表示经过节点v的s到t的最短路径条数，σ(s,t)表示s到t的最短路径条数。
特征向量中心性(Eigenvector Centrality)：取决于其相邻节点的数量，即该节点的度，也取决于其相邻节点的重要性。特征向量中心性基于其相邻节点的中心性来计算节点的中心性。节点i的特征向量中心性是由以下等式定义的向量y的第i个元素：
Ay＝λy
其中，A是网络的邻接矩阵，特征值为λ，y为解上述方程得出的特征向量。
有益技术效果：
本发明是一种基于乳腺癌疾病的基因调控网络构建及分析方法，在人类基因中有效地筛选出了乳腺癌相关基因，从而建立与乳腺癌疾病相关的基因调控网络，并通过节点中心性分析得到重要基因，可以促进从基因学的角度对乳腺癌疾病的研究，为找到干预乳腺癌疾病发生的有效途径奠定基础。
附图说明
图1为本发明实施例的基于乳腺癌疾病的基因调控网络的构建及分析方法流程图；
图2为本发明实施例的乳腺癌相关基因逐轮筛选策略单轮筛选过程流程图；
图3为本发明实施例的乳腺癌相关基因的基因调控网络构建过程流程图；
图4为本发明实施例的乳腺癌相关基因的基因调控网络节点分析过程流程图；
图5为本发明实施例的调控基因的网络；
图6为本发明实施例的互信息计算结果截图
图7为本发明实施例的乳腺癌相关基因的基因节点中心性分析计算结果截图。
具体实施方式
下面结合附图和具体实施实例对发明做进一步说明：一种基于乳腺癌疾病的基因调控网络的构建及分析方法，其流程如图1所示，包括乳腺癌相关基因初步筛选过程、基因调控网络建模过程和基因调控网络分析过程，具体包括：
步骤1：乳腺癌相关基因初步筛选：将乳腺癌基因分为基准基因、相关基因与筛选出的相关基因，基准基因指现有技术中已经认定与乳腺癌疾病有关系的基因，一共40个基因，相关基因指除现有技术以外的所有基因，本实施例取24236个基因，筛选出的相关基因指用本申请的方法筛选出与乳腺癌疾病关系紧密的基因，包括步骤1.1～步骤1.4：
步骤1.1：将基因表达数据分为两个矩阵，X矩阵为乳腺癌基准基因的表达数据，Y矩阵是待筛选的其余乳腺癌相关基因的表达数据，在X矩阵与Y矩阵中，每一行为一个基因，每一列为一个样本；
步骤1.2：矩阵Y中的向量Yj依次与X中的向量Xi计算每个基因之间的互信息值，Yj和Xi代表第j个和第i个基因的数据；
基因筛选过程中采用的互信息的方法来衡量基因之间的关系。互信息通常用于测量两个变量X和Y之间的可靠性，因此两个基因之间的相关性可以通过对基因表达数据的计算得出。两个离散变量X和Y之间的互信息可以通过熵的形式表现如下：
I(X,Y)＝H(X)+H(Y)-H(X,Y)
其中，H(X)和H(Y)分别为变量X和变量Y的熵，H(X,Y)为变量X和Y的联合熵，离散变量X的熵H(X)可以表示如下：
Figure BDA0001776525790000051
其中，p(x)是向量X的边缘概率分布函数。向量X和Y的联合熵可以表示为
Figure BDA0001776525790000052
其中，p(x,y)为X和Y的联合概率分布函数。熵通过高斯核概率密度估计可以表示如下：
Figure BDA0001776525790000053
其中，C是变量X的协方差矩阵，n为C中的变量(基因)数量，由此，向量X和Y的互信息可以表示为协方差的形式：
Figure BDA0001776525790000054
其中，|C(X)|是变量X的协方差矩阵的行列式，|C(Y)|是变量Y的协方差矩阵的行列式，|C(X,Y)|是变量X和Y的协方差矩阵的行列式。将互信息计算公式转化成协方差的形式可以很大程度上提高互信息计算的效率。
应用到本实施例中具体公式为：
Figure BDA0001776525790000055
其中，I(Xi,Yj)为向量Xi和Yj的互信息，|C(Xi)|是向量Xi的协方差矩阵的行列式，|C(Yj)|是向量Yj的协方差矩阵的行列式，|C(Xi,Yj)|是向量Xi和Yj的协方差矩阵的行列式。
如图6所示，为本发明实施例的互信息计算结果截图，探针编号是原始的基因表达数据中基因名称对应的探针编号，一个基因可能对应多个探针，最后一列数据代表基准基因来自第几轮的筛选出来的基因，0则代表是来自最开始的40个基因中的。
步骤1.3：人工选择阈值Isn，若互信息值大于等于Isn，则向量Yj对应的基因添加到X矩阵中，并在矩阵Y中删除，阈值Isn为动态选择的，根据每一轮的互信息结果，对每轮的互信息的值进行排序，Isn即为当前轮中第sn个阈值，筛选出来互信息值最大的前sn个基因，其中sn等于需要初步筛选出来的基因数量除以需要的轮数，需要初步筛选出来的基因数为200，本实施例轮数设置为10；
步骤1.4：重复步骤1.2至1.3，直到矩阵X中的基因数量达到预设定的值M个，M＝240；
图2为本发明具体实施方式中乳腺癌相关基因逐轮筛选策略单轮筛选过程流程图，每轮过程一样，直到筛选到需要初步筛选出来的基因数量，本实施例设置为200个；
所述基准基因包括基因名称如下：AKT1、AKT2、APC、ARID1A、ARID1B、ARID2、ASXL1、BAP1、BRCA1、BRCA2、CASP8、CDH1、CDKN1B、CDKN2A、CCND1、MDM2、99Mb、ERBB2、ZNF217、FGFR1/ZNF703、MYC、GATA3、KRAS、MAP2K4、MAP3K1、MAP3K13、MIL2、MIL3、NCOR1、NF1、PIK3CA、PTEN、RB1、SETD2、SF3B1、SMAD4、SMARCD1、STK11、TBX3和TP53。
步骤2：基因调控网络的构建：利用步骤1中获得基因的表达数据矩阵X，包括240个基因数据，通过BNFinder2工具进行调控网络构建，设置调控节点集合的数量，评分函数为BDE评分函数，得到一个概率网络，网络节点代表基因，有向边代表基因之间的调控关系；具体包括步骤2.1～步骤2.3：
步骤2.1：数据处理：根据X矩阵，确定基因个数及样本数量，并为每个基因设定潜在调控基因，如图5所示，为本发明实施例的调控基因的网络，箭头表示调控方向，虚线表示负调控，实线是正调控；
步骤2.2：设定基因调控网络结构的结构学习及参数学习，采用基于评分搜索的结构学习方法，遍历所有可能的结构，然后用BDE评分函数衡量各个结构，进而找出最好的结构；具体包括步骤2.2.1～步骤2.2.4：
步骤2.2.1：对于某个基因，首先获取该基因的潜在调控基因及其表达数据；
步骤2.2.2：为基因的每个潜在调控节点赋予权重；
步骤2.2.3：获取调控基因为空集的基因的BDE分数；
步骤2.2.4：遍历基因的潜在调控基因的集合，每次遍历结果即为该基因的一种调控结构，计算该结构的BDE分数，存入到一个堆栈中，遍历结束后，堆栈中留下的为评分最高的对应的调控基因的集合，堆栈堆的大小即为设定的调控节点集合的数量；
步骤2.3：根据网络的学习结果整合调控关系，得到基因调控网络结果；
图3为本发明具体实施方式中乳腺癌相关基因的基因调控网络构建过程流程图。利用上一步骤获得的240个基因的表达数据通过BNFinder2工具进行调控网络构建，调控节点集合的数量设置为5，评分函数为BDE评分函数，可以得到一个概率网络，网络节点代表基因，有向边表示基因之间的调控关系。在此过程中，依次基于每个基因学习它的调控网络，待全部基因的子网络都构建完后合并成一个全局网络，即最终的基因调控网络。
本实施方式中，基因调控网络的构建采用的BNFinder2工具，该工具的具体内容如下：
BNFinder2是一个构建贝叶斯网络准确并有效的工具，静态贝叶斯和动态贝叶斯网络(DBN)都适用。一个贝叶斯网络定义包括一个有向无环图(DAG)和一个条件概率表集合。DAG中每一个节点表示一个随机变量{X1，X2，……，Xn}，可以是可直接观测变量或隐藏变量，而有向边表示随机变量间的条件依赖；条件概率表中的每一个元素对应DAG中唯一的节点，存储此节点对于其所有直接前驱节点的联合条件概率节点X的联合概率可以表示成：
Figure BDA0001776525790000071
其中，pa(i)表示节点i之“因”，或称pa(i)是i的parents(父母)。
BNFinder2采用基于评分搜索的结构学习方法，遍历所有可能的结构并采用合适的评分函数衡量结构好坏，进而找出最好的结构。采用的评分函数为BDe评分函数(BayesianDirichlet)，利用先验知识和数据寻求概率最大的网络结构，函数如下：
Figure BDA0001776525790000072
其中，Γ为Gamma函数，mi'jk是数据D中节点Vi状态为k且父节点状态组合为j的样本数，m′ij＝∑km′ijk，αijk为服从Dirichlet分布的函数的参数，
Figure BDA0001776525790000073
步骤3：根据基因调控网络结果，进行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性，将BNFinder2得到的基因调控网络，通过Python的扩展包networkx计算各个基因的度中心性，接近中心性，中介中心性和特征向量中心性；
所述节点度中心性(Degree Centrality)：在一个网络图G＝(V,E)中，节点v的度dv指的是与v相连的E中边的数量，公式为：
Figure BDA0001776525790000074
其中，CB(v)为节点度中心性，m为该节点的度，n为基因的数量；
所述接近中心性(Closeness Centrality)：考量每个结点到其它结点的最短路的平均长度。距离其它结点越近，那么它的中心度越高，公式为：
Figure BDA0001776525790000075
其中，Ccl(v)为接近中心性，n为基因的数量，d(v,u)为节点v和u之间的最短路径。
所述中介中心性(Betweenness Centrality)：以经过某个节点的最短路径数目来刻画节点重要性的指标。一个结点充当“中介”的次数越高，它的中介中心度就越大，公式为：
Figure BDA0001776525790000081
其中，CB(v)为中介中心性，σ(s,t|v)表示经过节点v的s到t的最短路径条数，σ(s,t)表示s到t的最短路径条数。
特征向量中心性(Eigenvector Centrality)：取决于其相邻节点的数量(即该节点的度)，也取决于其相邻节点的重要性。特征向量中心性基于其相邻节点的中心性来计算节点的中心性。节点i的特征向量中心性是由以下等式定义的向量x的第i个元素：
Ay＝λy
其中，A是网络的邻接矩阵，特征值为λ，y为解上述方程得出的特征向量。
如图7所示截图为根据因调控网络结果，行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性计算结果。
步骤4：取网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据作为筛选出的相关基因,，N＝10。
图4为本发明具体实施方式中乳腺癌相关基因的基因调控网络节点分析过程流程图。对生成的基因调控网络进行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性。将BNFinder得到的基因调控网络通过Python的扩展包networkx计算各个基因的度中心性，接近中心性，中介中心性和特征向量中心性。依据节点中心性的分析结果，挑选出四种度量指标中值最高的10个节点，作为筛选出的相关基因的结果，如表1所示。
表1筛选出的相关基因的结果
Figure BDA0001776525790000082

Claims (6)

1.一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，包括以下步骤：
步骤1：乳腺癌相关基因初步筛选：将乳腺癌基因分为基准基因、相关基因与筛选出的相关基因，基准基因指与乳腺癌疾病有关系的基因，相关基因指除基准基因以外所有基因，筛选出的相关基因指用如下方法筛选出与乳腺癌疾病关系紧密的基因，包括步骤1.1～步骤1.4：
步骤1.1：将基因表达数据分为两个矩阵，X矩阵为乳腺癌基准基因的表达数据，Y矩阵是待筛选的其余乳腺癌相关基因的表达数据，在X矩阵与Y矩阵中，每一行为一个基因，每一列为一个样本；
步骤1.2：矩阵Y中的向量Yj依次与X中的向量Xi计算每个基因之间的互信息值，Yj和Xi代表第j个和第i个基因的数据；
步骤1.3：人工选择阈值Isn，若互信息值大于等于Isn，则向量Yj对应的基因添加到X矩阵中，并在矩阵Y中删除，阈值Isn为动态选择的数值，根据每一轮的互信息结果，对每轮的互信息的值大小进行排序，Isn即为当前轮中第sn个阈值，筛选出来互信息值最大的前sn个基因；
步骤1.4：重复步骤1.2至1.3，直到矩阵X中的基因数量达到预设定的值M个；
步骤2：基因调控网络的构建：利用步骤1中获得基因的表达数据矩阵X，通过BNFinder2工具进行调控网络构建，设置调控节点集合的数量，评分函数为BDE评分函数，得到一个概率网络，网络节点代表基因，有向边代表基因之间的调控关系；具体包括步骤2.1～步骤2.3：
步骤2.1：数据处理：根据X矩阵，确定基因个数及样本数量，并为每个基因设定潜在调控基因；
步骤2.2：设定基因调控网络结构的结构学习及参数学习，采用基于评分搜索的结构学习方法，遍历所有可能的结构，然后用BDE评分函数衡量各个结构，进而找出最好的结构；具体包括步骤2.2.1～步骤2.2.4：
步骤2.2.1：对于某个基因，首先获取该基因的潜在调控基因及其表达数据；
步骤2.2.2：为基因的每个潜在调控节点赋予权重；
步骤2.2.3：获取调控基因为空集的基因的BDE分数；
步骤2.2.4：遍历基因的潜在调控基因的集合，每次遍历结果即为该基因的一种调控结构，计算该结构的BDE分数，存入到一个堆栈中，遍历结束后，堆栈中留下的为评分最高的对应的调控基因的集合，堆栈堆的大小即为设定的调控节点集合的数量；
步骤2.3：根据网络的学习结果整合调控关系，得到基因调控网络结果；
步骤3：根据基因调控网络结果，进行节点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性，将BNFinder2得到的基因调控网络，通过Python的扩展包networkx计算各个基因的度中心性，接近中心性，中介中心性和特征向量中心性；
步骤4：取网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据作为为筛选出的相关基因。
2.根据权利要求1所述一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，所述网络节点的度中心性，接近中心性，中介中心性和特征向量中心性的前N个数据，N≤10，且矩阵X中的基因数量M>N。
3.根据权利要求1所述一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，所述步骤1.3中人工选择阈值Isn，其中sn等于需要筛选出来的基因数量除以需要的轮数。
4.根据权利要求1所述一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，所述基准基因包括基因名称如下：AKT1、AKT2、APC、ARID1A、ARID1B、ARID2、ASXL1、BAP1、BRCA1、BRCA2、CASP8、CDH1、CDKN1B、CDKN2A、CCND1、MDM2、99Mb、ERBB2、ZNF217、FGFR1/ZNF703、MYC、GATA3、KRAS、MAP2K4、MAP3K1、MAP3K13、MIL2、MIL3、NCOR1、NF1、PIK3CA、PTEN、RB1、SETD2、SF3B1、SMAD4、SMARCD1、STK11、TBX3和TP53。
5.根据权利要求1所述一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，所述互信息的计算方法，如下公式：
Figure FDA0003009548640000021
其中，I(Xi,Yj)为向量Xi和Yj的互信息，|C(Xi)|是向量Xi的协方差矩阵的行列式，|C(Yj)|是向量Yj的协方差矩阵的行列式，|C(Xi,Yj)|是向量Xi和Yj的协方差矩阵的行列式。
6.根据权利要求1所述一种基于乳腺癌疾病的基因调控网络构建及分析方法，其特征在于，所述点中心性分析，包括节点的度中心性，接近中心性，中介中心性和特征向量中心性，用如下公式进行计算：
节点度中心性：在一个网络图G＝(V,E)中，节点v的度dv指的是与v相连的E中边的数量，公式为：
Figure FDA0003009548640000022
其中，Cd(v)为节点度中心性，dv为该节点的度，n为基因的数量；
所述接近中心性：
Figure FDA0003009548640000031
其中，Ccl(v)为接近中心性，n为基因的数量，d(v,u)为节点v和u之间的最短路径；
所述中介中心性：
Figure FDA0003009548640000032
其中，CB(v)为中介中心性，σ(s,t|v)表示经过节点v的s到t的最短路径条数，σ(s,t)表示s到t的最短路径条数；
特征向量中心性：节点i的特征向量中心性是由以下等式定义的向量y的第i个元素：
Ay＝λy
其中，A是网络的邻接矩阵，特征值为λ，y为得到的特征向量。
CN201810972336.2A 2018-08-24 2018-08-24 一种基于乳腺癌疾病的调控网络构建及分析方法 Active CN109243523B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810972336.2A CN109243523B (zh) 2018-08-24 2018-08-24 一种基于乳腺癌疾病的调控网络构建及分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810972336.2A CN109243523B (zh) 2018-08-24 2018-08-24 一种基于乳腺癌疾病的调控网络构建及分析方法

Publications (2)

Publication Number Publication Date
CN109243523A CN109243523A (zh) 2019-01-18
CN109243523B true CN109243523B (zh) 2021-06-11

Family

ID=65067892

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810972336.2A Active CN109243523B (zh) 2018-08-24 2018-08-24 一种基于乳腺癌疾病的调控网络构建及分析方法

Country Status (1)

Country Link
CN (1) CN109243523B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109994151B (zh) * 2019-01-23 2020-10-16 杭州师范大学 基于复杂网络与机器学习方法的肿瘤驱动基因预测系统
CN110555530B (zh) * 2019-09-02 2022-11-08 东北大学 一种基于分布式的大规模基因调控网络构建方法
CN110675912B (zh) * 2019-09-17 2022-11-08 东北大学 一种基于结构预测的基因调控网络构建方法
CN110734977A (zh) * 2019-11-12 2020-01-31 山西医科大学 以sf3b1为靶点在制备预防或治疗乳腺癌的药物中的应用
CN112086187B (zh) * 2020-09-16 2022-04-19 电子科技大学 一种基于复杂网络的疾病进展路径挖掘方法
CN112802543B (zh) * 2021-01-14 2024-04-30 东北大学 一种基于概率图的基因调控网络分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101719195A (zh) * 2009-12-03 2010-06-02 上海大学 一种逐步回归基因调控网络的推断方法
CN101763528A (zh) * 2009-12-25 2010-06-30 深圳大学 基于贝叶斯网络的基因调控网络构建方法
CN105740651A (zh) * 2016-03-07 2016-07-06 吉林大学 一种特定癌症差异表达基因调控网络的构建方法
CN107016260A (zh) * 2017-03-30 2017-08-04 广东工业大学 一种基于跨平台基因表达数据的基因调控网络重建方法
CN107220525A (zh) * 2017-05-19 2017-09-29 浙江工业大学 基于rnn的基因调控网络构建与动态差异性分析方法
CN107358062A (zh) * 2017-06-02 2017-11-17 西安电子科技大学 一种双层基因调控网络的构建方法
CN107368702A (zh) * 2017-06-02 2017-11-21 西安电子科技大学 一种构建miRNA调控网络的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101719195A (zh) * 2009-12-03 2010-06-02 上海大学 一种逐步回归基因调控网络的推断方法
CN101763528A (zh) * 2009-12-25 2010-06-30 深圳大学 基于贝叶斯网络的基因调控网络构建方法
CN105740651A (zh) * 2016-03-07 2016-07-06 吉林大学 一种特定癌症差异表达基因调控网络的构建方法
CN107016260A (zh) * 2017-03-30 2017-08-04 广东工业大学 一种基于跨平台基因表达数据的基因调控网络重建方法
CN107220525A (zh) * 2017-05-19 2017-09-29 浙江工业大学 基于rnn的基因调控网络构建与动态差异性分析方法
CN107358062A (zh) * 2017-06-02 2017-11-17 西安电子科技大学 一种双层基因调控网络的构建方法
CN107368702A (zh) * 2017-06-02 2017-11-21 西安电子科技大学 一种构建miRNA调控网络的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Identification of breast cancer mechanism based on weighted gene coexpression network analysis;X Guo 等;《Nature》;20170811;全文 *
Network analysis of EMT and MET micro-RNA regulation in breast cancer;Diana Drago-García 等;《Nature》;20171019;全文 *
基于信息论的基因调控网络结构推测算法研究;刘炜;《中国博士学位论文全文数据库基础科技辑》;20180615(第6期);全文 *
基于转录调控网络研究的乳腺癌与系统性红斑狼疮免疫系统发病机理探寻;孔薇 等;《青岛科技大学学报（自然科学版）》;20160229;第37卷(第1期);全文 *

Also Published As

Publication number Publication date
CN109243523A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN109243523B (zh) 一种基于乳腺癌疾病的调控网络构建及分析方法
CN108846526A (zh) 一种二氧化碳排放量预测方法
CN108595916B (zh) 基于生成对抗网络的基因表达全谱推断方法
CN109994200A (zh) 一种基于相似度融合的多组学癌症数据整合分析方法
Mu et al. Multi-objective ant colony optimization algorithm based on decomposition for community detection in complex networks
CN109754122A (zh) 一种基于随机森林特征提取的bp神经网络的数值预测方法
CN107609583A (zh) 分类模型的参数优化方法、装置、计算机设备及存储介质
CN107169871B (zh) 一种基于关系组合优化和种子扩张的多关系社区发现方法
Davarynejad et al. A novel general framework for evolutionary optimization: Adaptive fuzzy fitness granulation
CN113032902A (zh) 一种基于机器学习优化的高速列车气动头部外形设计方法
CN110990718A (zh) 一种公司形象提升系统的社会网络模型构建模块
CN107527110A (zh) 基于布谷鸟搜索优化mgm（1，n）的年用电量预测方法
Cui et al. Improved genetic algorithm to optimize the Wi-Fi indoor positioning based on artificial neural network
CN103637800B (zh) 基于八段阻抗模型的人体成分分析的方法
CN111599406A (zh) 结合网络聚类方法的全局多网络比对方法
CN112633591A (zh) 一种基于深度强化学习的空间搜索方法及设备
CN113177675A (zh) 一种基于天牛群算法优化神经网络的空调冷负荷预测方法
CN116187835A (zh) 一种基于数据驱动的台区理论线损区间估算方法及系统
CN111429970B (zh) 基于极端梯度提升方法进行特征选择来获取多基因风险评分的方法及系统
CN108388769B (zh) 基于边驱动的标签传播算法的蛋白质功能模块识别方法
CN108460491B (zh) 一种基于启发式算法的时变随机环境下景区线路设计方法
CN113919729A (zh) 一种区域三生空间影响与协作水平评估方法及系统
CN113762370A (zh) 一种结合高斯随机场的深度网络集合生成方法
CN113762602A (zh) 一种风电场短期风速预测方法
CN116614860A (zh) 一种基于卡尔曼滤波的无线传感器网络分簇路由优化方法

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