图示：易基因组学技术研究内容

多组学关联分析方法

（1） 直接关联

一个基因的功能元件甲基化情况影响该基因的表达。

重叠分析

Pearson/Spearman 相关性分析

（2）模型关联

基于基因转录、蛋白质、代谢物等之间的上下游相互作用联系。

多元线性模型(multiple linear model)

（3）网络关联

基于分子功能和通路的富集性。

WGCNA module correlation

EMDN algorithm

SNF algorithm

多组学关联分析方法图例

（1）直接关联

① 重叠分析

特点：简单粗暴，也适用于样本量少的情况。

分析结果：韦恩图。

图例：WGBS + total RNA-seq，含DMR的差异表达基因、差异miRNA靶向的差异表达基因、差异siRNA靶向的差异表达基因三者之间的重叠关系分析

关联理论基础：DNA甲基化、miRNA和siRNA协同作用调控基因表达。

② 皮尔森/斯皮尔曼相关性分析

特点：准确计算相关性程度（R值），及其显著性（p值）。

分析结果：散点图（+拟合线）、相关性热图

图例：血液组织RRBS+RNA-seq，DMR的甲基化水平与差异表达的表达水平之间的皮尔森关联分析。

关联理论基础：DNA甲基化可直接调控基因表达。

图例：肠道宏基因组+宏病毒组，健康与炎症性肠病（UC）组分别计算噬菌体与细菌丰度之间的皮尔森相关性并发现差异。

关联理论基础：噬菌体可侵染细菌进而影响肠道菌群的结构。

图例：胚胎scRNA-seq+蛋白质组，两连续发育阶段之间的RNA和表蛋白质表达的相关性分析。

关联理论基础：RNA转录和蛋白质翻译具有上下游关系。

图例：肠道宏基因组+代谢组，不同管理状态下川金丝猴肠道中短链脂肪酸与肠道菌群的皮尔森相关性分析。

关联理论基础：肠道中有多类菌群可直接产生短链脂肪酸这类益生物质。

（2）模型关联

回归分析（regression analysis)是确定两组或两组以上变量间关系的统计方法。回归分析按照变量的数量分为一元回归和多元回归。两个变量使用一元回归，两个以上变量使用多元回归。

多元线性回归模型(multiple linear model)

特点：基于基因表达、蛋白质、代谢物等之间的直接和间接相互作用联系。

分析结果：关联网络图

图例：肠道菌群16S+血液代谢组+肝脏转录组，通过多元线性模型关联分析，筛选出若干优秀模型，组建低剂量抗生素饲喂促进仔猪快速生长的多组学调控网络。

关联理论基础：肠道菌群可通过产生代谢物进入血液，运输至肝脏影响肝脏细胞的基因表达。

（3）网络关联

细胞内所有大分子相互作用的集合，称为相互作用组（Interactome），是大多数基因型与表型关系的基础，可以用来指导解释组学技术检测到的变化如何干扰整个机体。

机体的分子响应和变化具有功能富集性、通路富集性。因此不同组学检测数据也具有相似的功能富集性和变化规律。

网络关联算法正是基于这些生物学理论基础。

基于WGCNA的共变关联网络分析（WGCNA module correlation)

基于表观模块的差异网络分析（EMDN algorithm）

相似性网络融合分析（SNF algorithm）

分析结果：关联网络图

①基于WGCNA的共变关联网络分析

利用组间差异基因鉴定共甲基化和共表达模块。

模块-模块相关性、模块表型相关性可以有效识别具有功能富集性的多组学变化模块。

图例：血液组织RRBS + RNA-seq，基因的差异甲基化模式与基因表达模式的共变关联网络分析。

前期直接关联得到的基因很少；

改变策略，采用基于WGCNA的共变关联网络分析，得到的共甲基化和共表达基因均富集于自噬相关通路。

关联理论基础：基因组DNA的甲基化与基因表达变化具有功能富集性。

图例：肠道宏基因组+代谢组，不同管理状态下川金丝猴粪便代谢物与肠道菌群基于WGCNA的共变关联网络分析。

鉴定了2对强正相关的物种和代谢物的共变化模块。

圈养条件的代表性模块中发现了潜在致病菌和相关代谢物。

关联理论基础：肠道菌群可响应环境变化改变肠道微环境中相关的代谢产物的浓度。

基于表观模块的差异网络分析（EMDN algorithm）利用组间差异基因鉴定共甲基化（共表达）模块。

差异共甲基化（差异共表达）网络筛选。

从多个差异共变网络中筛选共同网络。

图例：基于表观模块的差异网络分析（EMDN algorithm）

相似性网络融合分析 (SNF algorithm）

图例：SNF算法(图1d)使用了一种基于信息传递理论的非线性算法，该方法迭代更新每个网络，使其与其他网络更加相似。经过几次迭代之后，SNF收敛到单个网络中。

算法的优点是，弱相似性(低权重的边)消失，有助于降低噪声，而在一个或多个网络中存在的强相似性(高权重的边)被添加到其他网络中。

图例：融合网络三个cluster内部的连接性、紧密性和cluster之间相对较少的边界，说明该算法可以更清晰地显示多形性成胶质细胞瘤(GBM)患者的分型情况。

其他网络关联分析方法

从关联走向因果：组学分子实验验证

基因表达相关的组学：

基因敲除/抑制

基因过表达

甲基化组学：

甲基化酶基因的敲除与过表达

宏基因组（肠道菌群）：

无菌动物模型

粪菌移植

参考文献：

Yan H, Bombarely A, Xu B, et al.Autopolyploidization in switchgrass alters phenotype and flowering time viaepigenetic and tranion regulation[J]. Journal of experimental botany,2019

Zhang, D., Hu, Q., … Gao, F. (2019). Epigenetic and transcriptional signatures of ex situ conserved golden snub-nosed monkeys (Rhinopithecus roxellana). Biological Conservation, 237, 175–184.

Zuo, T., Lu, X. J., Zhang, Y. (2019). Gut mucosal virome alterations in ulcerative colitis. Gut, 68(7), 1169–1179.

Gao, Y., Liu, X.. (2017). Protein Expression Landscape of Mouse Embryos during Pre-implantation Development. Cell Reports, 21(13), 3957–3969.

Evaluating the influence of conservation activities on the gut microecosystem of Rhinopithecus roxellana based on establishment of a gut microbiome gene catalog. Unpublished.

Inter-correlated gut microbiota and SCFAs changes upon antibiotics exposure links with rapid body-mass gain in weaned piglet model. The Journal of nutritional biochemistry,2019, 74: 108246.

Multiple network algorithm for epigenetic modules via the integration of genome-wide DNA methylation and gene expression data. BMC Bioinformatics , 2017,18(1), 1–13.

Similarity network fusion for aggregating data types on a genomic scale. Nature Methods , 2014,11(3), 333–337.

Di Nanni, N., Bersanelli, M.(2020). Network Diffusion Promotes the Integrative Analysis of Multiple Omics. Frontiers in Genetics,11(February), 1–12.