加权相关网络分析

✍ dations ◷ 2025-11-04 07:31:50 #加权相关网络分析

加权相关网络分析也称为加权基因共表达网络分析,是一种广泛使用的数据挖掘方法,它用两两变量间相关系数研究生物网络。它适用于高维度数据集,在基因组学领域应用的最为广泛。

此分析假定基因网络服从无尺度分布,定义基因共表达相关矩阵,计算生成邻接矩阵,然后计算各个节点的相异系数,建立分层聚类树。此树不同分支代表不同的基因模块,模块内基因有较高的共表达水平,从中可以鉴定出枢纽基因。可通过基因网络与表型之间联系,最终找到表型的靶点基因、基因网络。此分析还可以用于研究共表达模块之间的关系、比较不同网络的网络拓扑学(差异网络分析)。加权基因共表达网络分析可以进行数据提炼(与斜交因子分析相关)、聚类分析(模糊聚类)、特征选择方法(比如,进行基因筛选)、综合互补(基因组)数据的框架(根据数量变量之间的加权相关性)和探索性数据分析。

虽然加权基因共表达网络分析结合了很多传统数据探索技术,易于理解的网络语言和分析框架使它优于任何标准的分析方法。由于使用网络方法,同时适合整合互补基因组的数据,它可认为是一种系统生物学或系统遗传数据的分析方法。加权基因共表达网络分析通过在相似的模块间选择模块间枢纽基因,也使基于网络的元分析技术成为可能。

加权基因共表达网络分析方法是加州大学洛杉矶分校大卫·格芬医学院人类遗传学教授、加州大学洛杉矶分校菲尔丁公共卫生学院生物统计学教授Steve Horvath、他在加州大学洛杉矶分校的同事、实验室成员(主要有Peter Langfelder、Bin Zhang、Jun Dong)联合创立的。研究的很大部分是应用领域的研究者合作产生的。加权相关网络是与癌症研究者Paul Mischel、Stanley F. Nelson和神经科学家Daniel H. Geschwind、Michael C. Oldham的讨论后发明的(根据里的致谢章节)。存在很多关于依赖网络、自由尺度网络和共表达网络的研究文献。

第一,指定基因共表达相似性指标来定义网络。我们用 s i j {displaystyle s_{ij}} 表示基因i和j之间的基因共表达相似性。很多共表达研究都使用相关性的绝对值,作为无符号共表达相似指标,

s i j u n s i g n e d = | c o r ( x i , x j ) | {displaystyle s_{ij}^{unsigned}=|cor(x_{i},x_{j})|}

公式里,多个样本中基因表达量构成了基因表达谱 x i {displaystyle x_{i}} x j {displaystyle x_{j}} 。然而,使用相关性的绝对值就不区分基因的抑制和激活,造成了生物学信息的损失。而在有符号网络中,基因之间的相似性表达了其表达谱的相关方向。为定义基因表达谱 x i {displaystyle x_{i}} x j {displaystyle x_{j}} 之间一个有符号共表达指标,我们可以使用对相关性的简单转换:

s i j s i g n e d = 0.5 + 0.5 c o r ( x i , x j ) {displaystyle s_{ij}^{signed}=0.5+0.5cor(x_{i},x_{j})}

与无符号指标 s i j u n s i g n e d {displaystyle s_{ij}^{unsigned}} 一样,有符号相似性 s i j s i g n e d {displaystyle s_{ij}^{signed}} 取0和1之间的值。请注意,两个表达量相反的基因在无符号时等于1,有符号时等于0;两个表达量无相关性的基因,无符号时等于0,有符号时等于0.5。

下一步,使用邻接矩阵(网络) A = {displaystyle A=} 定量确定两个基因相互连接的程度。将共表达相似性矩阵 S = {displaystyle S=} 设定阈值定义为 A {displaystyle A} 。对相似性指标 S {displaystyle S} 设立固定的阈值会将数据分为两组,生成无权重基因共表达网络。如果 s i j > τ {displaystyle s_{ij}>tau } ,无权重网络邻接矩阵设定为1,否则设定为0。因为设定固定的阈值会将基因连接分成两组,对阈值的选择较为敏感,造成共表达信息的损失。如果不设固定阈值,共表达信息连续的特性就得以保留。加权基因共表达网络分析使用以下的幂函数评估连接强度:

a i j = ( s i j ) β {textstyle a_{ij}=(s_{ij})^{beta }} ,

β {displaystyle beta } 是不固定阈值的参数。无符号和有符号网络分别使用 β = 6 {displaystyle beta =6} β = 12 {displaystyle beta =12} 作为默认值。 β {displaystyle beta } 也可以使用无尺度拓扑标准进行选择,即选择 β {displaystyle beta } 可达成无尺度网络的最小值 。

由于 l o g ( a i j ) = β l o g ( s i j ) {displaystyle log(a_{ij})=beta log(s_{ij})} ,加权网络邻接矩阵与共表达相似矩阵线性的对数值相关。 β {displaystyle beta } 取高值,会将高度相似性转换为高度邻接矩阵、低相似度推向0。由于这种应用于成对相关矩阵的不固定阈值的步骤会产生加权邻接矩阵,随后的分析称为加权基因共表达网络分析。

模块中心分析的重要一步是:用网络拓扑指标将基因集聚为网络模块。简单来说,一对基因如果相互连接,距离会较近。按照惯例,两个基因的距离在0到1之间。通常来说,加权基因共表达网络分析使用拓扑重叠指标作为距离。

拓扑重叠指标将两个基因之间的邻接矩阵和两个基因与其他基因之间的连接强度混合,是网络相互连接性稳健的指标。将此指标输入平均连锁层级聚类。动态分支切割法生成的聚类结果中的分支,就成了模块。之后,在特定模块内的基因与模块特征向量一起,可以被视为标准模块表达数据最好的总结。标准化表达谱的第一个主成分定义为模块特征向量。为了寻找与感兴趣的临床性状相关的模块,模块特征向量之后与其进行相关性比较,得出特征向量显著指标。人们可以在模块特征向量之间建立共表达网络,即以节点为模块的网络。

为了判断特定模块的模块间枢纽基因,可以使用两种连接指标。第一种称为 k M E i = c o r ( x i , M E ) {displaystyle kME_{i}=cor(x_{i},ME)} ,是将基因与相关模块的特征向量进行相关定义的。第二种称为kIN,定义为邻接矩阵相关模块基因的和。实际运用中,两种指标是等同的。为了测试一个模块是否存在于另一个数据集,可以使用多种网络统计量,如 Z s u m m a r y {displaystyle Zsummary}

可以将加权相关网络理解成特殊形式的加权网络、依赖网络或者相关网络。加权相关网络分析由于以下几点,很具吸引力:

加权基因共表达网络分析被广泛应用于在分析表达数据(即转录数据)。如,寻找模块间枢纽基因。

在系统遗传学应用中,此分析常用于数据提炼步骤,模块由“模块特征向量”表示。可用模块特征向量将模块与临床性状相关联。特征向量网络是模块特征向量之间的共表达网络(即节点为模块的网络)。 神经科学应用加权基因共表达网络分析,分析微阵列数据、单细胞RNA测序数据、DNA甲基化数据、miRNA数据、肽计数、人类微生物组计划数据(16S rRNA基因序列)的基因组数据。其他的应用包括:脑影像数据,如功能性磁共振成像数据。

加权基因共表达网络分析的R软件包提供了加权网络分析所有功能(模块建立、节点基因选择、模块保存统计、差异网络分析、网络统计)的函数。可以通过CRAN下载加权基因共表达网络分析包。

相关

  • 股疝气股疝是一种腹外疝,是指腹腔内的器官或组织连同腹膜壁层形成的疝囊通过股环经股管向卵圆窝突出。其发病率约占腹外疝的3~5%,多见于中老年妇女。根据传统定义,股疝不属于腹股沟疝,
  • 药物经济学药物经济学是一门应用经济学原理和方法来研究和评估药物治疗的成本与效果及其关系的边缘学科。药物经济学的研究任务主要是通过成本分析对比不同的药物治疗方案或药物治疗方
  • 科曼多尔群岛坐标:55°N 167°E / 55°N 167°E / 55; 167科曼多尔群岛(俄语:Командорские острова)是俄罗斯在太平洋沿岸的一个火山岛群,位于堪察加半岛以东200公里的太
  • 甬剧甬剧,中国戏曲剧种之一,流行于宁波地区。舞台语言采用宁波方言。甬剧曲调包含基本调、快二簧、慢二簧等。甬剧剧目有《田螺姑娘》《两兄弟》《王鲲》《亮眼哥》《典妻》等。目
  • 格拉塞伦巴赫格拉塞伦巴赫(德语:Grasellenbach)是德国黑森州的一个市镇。总面积22.87平方公里,总人口3772人,其中男性1939人,女性1833人(2011年12月31日),人口密度165人/平方公里。
  • 洪弃生洪弃生(1866年12月17日-1928年),原名攀桂,字月樵,台湾清末、日治时期重要文学家。彰化鹿港人,台湾割让后改名洪�,字弃生,以字行。清朝光绪十年(1884年)以第一名入学生员,光绪廿一年(1895年
  • 兰菁轩兰菁轩 (1999年12月16日-),浙江人,中国足球运动员,司职中后卫,目前效力于天津泰达。最后更新:2020年8月7日
  • 在线考试在线考试是指通过局域网、城域网进行考试的一种考试形式,也可以说成是通过网络媒体进行的考试。  
  • 马尔玛杜克·比克达尔穆罕默德·马尔玛杜克·比克达尔 (英语:Muhammad Marmaduke Pickthall,1875年4月7日-1936年5月19日) ,出生时名为马尔玛杜克·威廉·比克达尔(英语:Marmaduke William Pickthall),是
  • 网络赌博在线赌博又称网络赌博、网络赌博、网上赌博、网上博彩,是指人们在互联网上进行的赌博。种类包括网上扑克、在线赌场和在线体育赌博。第一个能公开访问的在线赌博场所,是1994年10月列支敦士登发行的国际彩票。许多国家限制或禁止网上赌博。然而在美国的一些州、加拿大的一些省份、欧洲联盟的大多数国家和加勒比地区的一些国家,网上赌博是合法的。