基因调控网络 (Gene Regulatory Network) 01
本文为入门级的基因调控网络文章,主要介绍一些基本概念及常见的GRN模型。
概念:基因调控网络 (Gene Regulatory Network, GRN),简称调控网络,指细胞内或一个基因组内基因和基因之间的相互作用关系形成的网络,特指基因调控 (gene regulation) 导致基因之间的作用。
GRN是生物体内控制基因表达的机制,基因表达的主要过程是转录+翻译
GRN构建方法:
多数方法使用静态数据来分析基因网络,如基因表达矩阵,该矩阵是一个时刻的基因表达情况。实际上,我们需要考虑动态网络,这样才能逼近真实的GRN。
一些网络模型:
1、布尔网络
布尔网络是一种最简单的模型。布尔网络中,每个基因的状态只有“开”和“关”,“开”表示基因有表达,“关”表示基因未表达。基因间相互作用由布尔表达式表示:and, or, not,如 A and not B -> C。
该网络过于简化,存在局限。
2、线性模型
线性模型是一种连续的GRN模型。在线性模型中,一个基因的表达水平由若干其他基因表达水平的加权和表示,权是基因之间相互关系的定量化:正权表示基因激发,负权表示基因抑制,0权表示两个基因没有关系。
Xi(t+Δt)=∑wij Xj(t)+η
该网络是一种简单的数学模型,只能处理具有线性关系的基因表达数据,应用范围小。
相关模型:加权矩阵模型
3、马尔可夫模型
马尔可夫链是一种随机过程,适用于分析时间序列的基因表达数据。在马尔可夫模型中,马尔可夫链假设某一时刻的基因表达水平决定了下一时刻的基因表达水平,公式如下:
C(t)=J C(t-1)
构建GRN过程中,基于马尔可夫模型对gene expression profile的特征提取和聚类都表现出良好的适应性。
如果要提高模型的准确性,可提高马尔可夫模型的阶数。
4、微分方程模型
微分方程模型假设一个基因为一个变量,由n个基因组成的网络可以由如下n维微分方程表示:
dxi(t)/dt=fi(x1,x2,…,xn),xi(t)是第i个基因的表达水平,n表示网络中的基因数。
根据生物数据和建模的不同要求,函数fi根据实际确定。通常有如下两种形式:
fi(xj,Θij,α)=xjα/(xj+θijα)
fi(xj,α)=e-αxj/(1+e-αxj)
微分方程构建GRN的优点:强大灵活,有利于描述基因网络中的复杂关系。
5、贝叶斯网络模型
以贝叶斯定理和假设为理论基础,用有向无环图 (DAG) 的形式表示随机变量间的概率关系,网络中每个基因是一个节点,每个调控关系是一条边。
该模型可以处理随机事件,控制噪声,可以获得变量间的因果关系,在GRN模型中,贝叶斯网络比其他模型更有优势。
相关模型:神经网络模型,图解高斯模型
6、互信息关联模型
互信息关联模型用熵和互信息描述基因之间的关联。
一个基因表达模式A的熵,P(xi)是基因表达值出现在xi的频率,n为表达水平的区间数目。熵越大,基因表达水平越趋近随机分布。
两个基因表达模式之间的互信息MI(A,B)=H(A)+H(B)-H(A,B),若MI(A,B)=0,则两个基因不相关,若MI(A,B)越大,两个基因越非随机相关,之间的生物关系越密切。
7、随机方程模型
......
构建GRN模型时可参考的一些经验:GRN是稀疏的,也就是并不是所有基因之间都有相互作用关系 [3];基因调控网络的大多数变量是连续的,粗糙的离散化会丢失信息,精细的离散化参数太多,最好直接用连续变量 [4];
Ref:
[1]. https://baike.baidu.com/item/%E8%B0%83%E6%8E%A7%E7%BD%91%E7%BB%9C/5844691
[2]. https://wenku.baidu.com/view/34dff5ef19e8b8f67c1cb958.html
[3]. Maetschke S R, Ragan M A. Characterizing cancer subtypes as attractors of Hopfield networks[J]. Bioinformatics, 2014, 30(9): 1273-1279.
[4]. 雷耀山,史定华,王翼飞.基因调控网络的生物信息学研究[J].自然杂志,2004(01):7-12.
基因调控网络 (Gene Regulatory Network) 01的更多相关文章
- 基于基因调控网络(Hopfield network)构建沃丁顿表观遗传景观
基因调控网络的概念在之前已经简要介绍过:https://www.cnblogs.com/pear-linzhu/p/12313951.html 沃丁顿表观遗传景观(The Waddington's e ...
- 【转】理解Docker容器网络之Linux Network Namespace
原文:理解Docker容器网络之Linux Network Namespace 由于2016年年中调换工作的原因,对容器网络的研究中断过一段时间.随着当前项目对Kubernetes应用的深入,我感觉之 ...
- VMware虚拟机上网络连接(network type)的三种模式--bridged、host-only、NAT
VMware虚拟机上网络连接(network type)的三种模式--bridged.host-only.NAT VMWare提供了三种工作模式,它们是bridged(桥接模式).NAT(网络地址转换 ...
- [USACO08JAN]手机网络Cell Phone Network
[USACO08JAN]手机网络Cell Phone Network 题目描述 Farmer John has decided to give each of his cows a cell phon ...
- linux 网络虚拟化: network namespace 简介
linux 网络虚拟化: network namespace 简介 network namespace 是实现网络虚拟化的重要功能,它能创建多个隔离的网络空间,它们有独自的网络栈信息.不管是虚拟机还是 ...
- 洛谷 P2812 校园网络【[USACO]Network of Schools加强版】 解题报告
P2812 校园网络[[USACO]Network of Schools加强版] 题目背景 浙江省的几所OI强校的神犇发明了一种人工智能,可以AC任何题目,所以他们决定建立一个网络来共享这个软件.但是 ...
- Solaris11.1网络配置(Fixed Network)
Solaris11的网络配置与Solaris10有很大不同,Solaris11通过network configuration profiles(NCP)来管理网络配置. Solaris11网络配置分为 ...
- 洛谷P2899 [USACO08JAN]手机网络Cell Phone Network
P2899 [USACO08JAN]手机网络Cell Phone Network 题目描述 Farmer John has decided to give each of his cows a cel ...
- (58)zabbix网络拓扑图配置network map
zabbix网络地图介绍 “zabbix network map”可以简单的理解为动态网络拓扑图,可以针对业务来配置zabbix map, 通过map可以了解应用的整体状况:服务器是否异常.网络是否有 ...
随机推荐
- 爬虫(十七):Scrapy框架(四) 对接selenium爬取京东商品数据
1. Scrapy对接Selenium Scrapy抓取页面的方式和requests库类似,都是直接模拟HTTP请求,而Scrapy也不能抓取JavaScript动态谊染的页面.在前面的博客中抓取Ja ...
- 030-PHP日期查询函数
<?php , , ))//检查日期函数 { print("2,18,1970 :" . "这是一个正确的日期格式"); } else { print(& ...
- 【转】AutoMapper对象映射
什么是AutoMapper?AutoMapper是一个简单的小型库,用于解决一个看似复杂的问题 - 摆脱将一个对象映射到另一个对象的代码.这种类型的代码是相当沉闷和无聊的写,所以为什么不发明一个工具来 ...
- Sklearn 预处理数据
## 版权所有,转帖注明出处 章节 SciKit-Learn 加载数据集 SciKit-Learn 数据集基本信息 SciKit-Learn 使用matplotlib可视化数据 SciKit-Lear ...
- centos破解压缩文件密码
rarcrack是linux系统端的一款破解加密压缩包的工具,rarcrack使用的穷举法进行破解已经加密的rar.zip和7z压缩包,支持设置多线程,和文件类型进行破解. 1.rarcrack下载页 ...
- POJ 3295:Tautology
Tautology Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 10482 Accepted: 3982 Descri ...
- Python安装和虚拟环境创建以及外部库的安装
Python.虚拟环境.外部库的安装 一 安装Python 1 Windows 到官网下载对应的版本 下载地址 我选择的是Python3.6.8 下载完成后双击运行 !!!勾选Add Python 3 ...
- UVALive 4043 转化最佳完美匹配
首先黑点和白点是组成一个二分图这毫无疑问 关键是题目中要求的所有黑白配的线不能交叉...一开始我也没想到这个怎么转化为二分图里面的算法. 后来看书才知道,如果两两交叉,则可以把两根线当四边形的对角线, ...
- Aizu 2155 Magic Slayer 背包DP
这是上上次对抗赛的题目了 其实现在发现整个代码从头到尾,都是用了背包,怪我们背包没深入学好. 比赛的时候,聪哥提出的一种思路是,预处理一下,背包出 ALL攻击 和 single攻击的 血量对应的最小花 ...
- 浅谈__slots__
__slots__在python中是扮演属性声明(Attribute Declaration)的角色,当然属性本身不用声明,只需要在使用前赋值即可,但是实际上,属性的应用远非在使用前赋值就行,所以有a ...