加权相关网络分析

✍ dations ◷ 2025-09-12 16:06:55 #加权相关网络分析

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

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

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

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

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

${displaystyle s_{ij}^{unsigned}=|cor(x_{i},x_{j})|}$ ${displaystyle s_{ij}^{unsigned}=|cor(x_{i},x_{j})|}$

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

${displaystyle s_{ij}^{signed}=0.5+0.5cor(x_{i},x_{j})}$ ${displaystyle s_{ij}^{signed}=0.5+0.5cor(x_{i},x_{j})}$

与无符号指标 ${displaystyle s_{ij}^{unsigned}}$ ${displaystyle s_{ij}^{unsigned}}$ 一样，有符号相似性 ${displaystyle s_{ij}^{signed}}$ ${displaystyle s_{ij}^{signed}}$ 取0和1之间的值。请注意，两个表达量相反的基因在无符号时等于1，有符号时等于0；两个表达量无相关性的基因，无符号时等于0，有符号时等于0.5。

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

${textstyle a_{ij}=(s_{ij})^{beta }}$ ${textstyle a_{ij}=(s_{ij})^{beta }}$ ,

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

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

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

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

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

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

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

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

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

相关

翠竹街道翠竹街道，是中华人民共和国广东省深圳市罗湖区下辖的一个乡镇级行政单位。翠竹街道下辖以下地区：百仕达社区、新村社区、翠岭社区、水贝社区、翠辉社区、民新社区、木头龙社区
兰亭兰亭古迹是中国浙江绍兴的重要文化古迹和旅游胜地，位于绍兴城区西南13公里的兰渚山麓。东晋永和九年（353年）三月初三，王羲之与谢安、王献之等40多位名士在此举办修禊集会，王羲之
区域菊石区域菊石（学名：），又名恩德莫菊石，是生存在白垩纪海洋中的一属菊石。它们是区域菊石超科的代表属，但其化石却相当地稀有。其化石被发现于德国及波兰等地。区域菊石和与其同科的其他
星震学星震学（英文：Asteroseismology，来自古希腊文 ἀστήρ，，恒星、σεισμός, ，振动、-λογία, ，研究。或称为stellar seismology）是借由分析恒星震动频谱研究恒星内部结构的
国际数学物理大会国际数学物理大会（International Congress on Mathematical Physics, ICMP）是数学物理学界最大的研究性学术会议。这一会议每三年举办一次，代表国际数学物理联合会（IAMP）。庞加莱
乌塔尔帕拉科特龙格乌塔尔帕拉科特龙格（Uttarpara Kotrung），是印度西孟加拉邦Hugli县的一个城镇。总人口150204（2001年）。该地2001年总人口150204人，其中男性78661人，女性71543人；0—6岁人口12292人，其
奇异建筑奇异建筑是后现代主义建筑中的一个流派。其特点是各不相同。有人说它的创作目的就是成名，把建筑当作游戏。
众议院 (波兰) 政治主题波兰共和国众议院（波兰语：Sejm Rzeczypospolitej Polskiej），通称音译瑟姆（Sejm，试听），是波兰议会的下议院。其自波兰立陶宛时代起存在，建立了波兰民主制的传统，并一路延续到现代。在20世纪前，词语“Sejm”指的是由三大议院（英语：Chambers of parliament）组成的整个波兰议会，包括下院（Izba Poselska）、上院（Senat）以及国王，是称呼三大阶层国会的常用名词。然而随着波兰第二共和国（1918年－1939年）的建立，“Sejm”
穗花澳大利亚杉穗花澳大利亚杉()为澳大利亚杉属()中唯一的物种，并且为太平洋地区新喀里多尼亚大特雷岛（英语：Grande Terre (New Caledonia)）的特有种。穗花澳大利亚杉为常绿乔木，高度可生长3至25米。树梢为密集分支，然其在拥有自由生长空间时，具有不会生长到如此高大及繁茂的习性（英语：Habit (biology)）。树皮为红棕色纤维状，会剥落成短带或鳞片状。树叶呈螺旋状排列于树枝上，叶柄长达8毫米。深绿色的树叶长17公分、宽0.8公分，叶型呈线形与披针形间，前端为尖状，叶缘略为卷起。在树叶的上侧
2016年10月台湾宣告 - 检阅