基因网络 python

转载

岁月静好呀 2024-09-12 15:49:03

文章标签 基因网络 python 互信息线性模型数据 文章分类 Python 后端开发

本文为入门级的基因调控网络文章，主要介绍一些基本概念及常见的GRN模型。

概念：基因调控网络 (Gene Regulatory Network, GRN)，简称调控网络，指细胞内或一个基因组内基因和基因之间的相互作用关系形成的网络，特指基因调控 (gene regulation) 导致基因之间的作用。
GRN是生物体内控制基因表达的机制，基因表达的主要过程是转录+翻译

基因网络 python_互信息

GRN构建方法：
多数方法使用静态数据来分析基因网络，如基因表达矩阵，该矩阵是一个时刻的基因表达情况。实际上，我们需要考虑动态网络，这样才能逼近真实的GRN。
一些网络模型：
1、布尔网络
布尔网络是一种最简单的模型。布尔网络中，每个基因的状态只有“开”和“关”，“开”表示基因有表达，“关”表示基因未表达。基因间相互作用由布尔表达式表示：and, or, not，如 A and not B -> C。
该网络过于简化，存在局限。
2、线性模型
线性模型是一种连续的GRN模型。在线性模型中，一个基因的表达水平由若干其他基因表达水平的加权和表示，权是基因之间相互关系的定量化：正权表示基因激发，负权表示基因抑制，0权表示两个基因没有关系。
X_i(t+Δt)=∑w_ij X_j(t)+η
该网络是一种简单的数学模型，只能处理具有线性关系的基因表达数据，应用范围小。

相关模型：加权矩阵模型

3、马尔可夫模型

马尔可夫链是一种随机过程，适用于分析时间序列的基因表达数据。在马尔可夫模型中，马尔可夫链假设某一时刻的基因表达水平决定了下一时刻的基因表达水平，公式如下：

C(t)=J C(t-1)

构建GRN过程中，基于马尔可夫模型对gene expression profile的特征提取和聚类都表现出良好的适应性。

如果要提高模型的准确性，可提高马尔可夫模型的阶数。

4、微分方程模型

微分方程模型假设一个基因为一个变量，由n个基因组成的网络可以由如下n维微分方程表示：

dx_i(t)/dt=f_i(x₁,x₂,…,x_n)，x_i(t)是第i个基因的表达水平，n表示网络中的基因数。

根据生物数据和建模的不同要求，函数fi根据实际确定。通常有如下两种形式：

f_i(x_j,Θ_ij,α)=x_jα/(x_j+θ_ijα)

fi(xj,α)=e^-αx_j/(1+e^-αx_j)