Name: 基于加权相关网络的不同栖息地根微生物群的差异 (视频) | JoVE Journal
Uploaded: 2021-09-25T13:30:00.000+00:00
Duration: 9 min 49 s
Description: 学习使用WGCNA算法构建水稻根系不同生态位微生物共现网络，识别核心属和非保守模块。该方法可揭示根际土壤、根表和根内微生物群落的差异及其生态适应机制。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本论文的撰写得到了中国国家自然科学基金-贵州省人民政府喀斯特科学研究中心项目（U1812401）、贵州师范大学博士研究项目（GZNUD[2017]1）、贵州省科技支撑计划项目（QKHZC[2021]YB459）以及贵阳市科技计划项目（[2019]2-8）的资助。
作者感谢 Edwards J.A 等人提供了水稻微生物组数据至公共数据库，并感谢 TopEdit（www.topeditsci.com）在本手稿撰写过程中提供的语言学协助。

1. 数据下载
<ol>
	<li>从NCBI数据库下载登录号PRJNA386367的数据。从该登录号的数据中，选取2014年在美国加利福尼亚州阿布克尔（Arbuckle）的淹水水稻田中种植14周的水稻植株的根际、根表和内生微生物组数据。 
	&#8203;注：根际、根表和内生微生物组数据在登录号PRJNA386367中以OTU表形式提供。</li>
</ol>
2. 最佳功率值的确定
注意：WGCNA 软件包包含以下所有功能参数。WGCNA 是一个用于加权相关网络分析的 R 软件包。主要命令行参见补充材料 S1。
<ol>
	<li>在 R 语言环境中，打开 Rstudio 软件并安装 WGCNA 软件包。</li>
	<li>加载数据，并使用 goodSamplesGenes 函数检查数据的正确性。执行以下命令行： 
	&#34;gsg = goodSamplesGenes(datExpr0, verbose = 3) 
	gsg$allOK &#34; 
	点击 Run。</li>
	<li>检查是否存在离群值，并保存符合要求的样本。当检查结果为 TRUE 时，继续下一步操作。保存结果。</li>
	<li>使用 PickSoftThreshold 函数计算两组数据在不同幂值下的无标度指数 R2。执行以下命令行： 
	&#34;sft = pickSoftThreshold(datExpr0, powerVector = powers, verbose = 5)&#34; 
	点击 Run。</li>
	<li>可视化结果（图 1）。执行以下命令行： 
	&#34;plot(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2], 
	&#160; &#160;xlab=&#34;软阈值（power）&#34;,ylab=&#34;无标度拓扑模型拟合度，有符号 R^2&#34;,type=&#34;n&#34;, 
	&#160; &#160;main = paste(&#34;ES_无标度性&#34;)); 
	text(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2], 
	&#160; &#160;labels=powers,cex=cex1,col=&#34;red&#34;); 
	abline(h=0.9,col=&#34;red&#34;) 
	plot(sft$fitIndices[,1], sft$fitIndices[,5], 
	&#160; &#160;xlab=&#34;软阈值（power）&#34;,ylab=&#34;平均连接度&#34;, type=&#34;n&#34;, 
	&#160; &#160;main = paste(&#34;ES_平均连接度&#34;)) 
	text(sft$fitIndices[,1], sft$fitIndices[,5], labels=powers, cex=cex1,col=&#34;red&#34;)&#34; 
	点击 Run。 
	注：加权相关网络算法的前提是所构建的共表达网络结构符合无标度拓扑标准，从而提高其鲁棒性。无标度指数越接近 1，表明网络结构越接近无标度网络。</li>
	<li>选择无标度指数 R2 大于 0.9 时对应的幂值，并进入下一步分析。 
	&#8203;注：当无标度指数接近 1 时，网络结构更接近无标度网络。在分析两个或多个网络时，需选择使每个网络均接近无标度网络的幂值，以满足共表达网络之间的可比性。</li>
</ol>
3. 共表达网络的构建与模块识别
注意：根据上述计算得到的功率值构建共现网络。关键命令行参见补充材料 S2。
<ol>
	<li>使用 WGCNA 软件包中的 adjacency 函数为构建符号共现网络添加有符号参数。执行命令行： 
	&#34;adjacency = adjacency(datExpr0, power = softPower)&#34; 
	点击 Run。</li>
	<li>应用 TOM-similarity 函数构建拓扑重叠网络，并计算不相似性网络。执行命令行： 
	&#34;TOM = TOMsimilarity(adjacency); 
	dissTOM = 1-TOM&#34; 
	点击 Run。 
	注意：添加有符号参数用于设定拓扑重叠网络的类型。</li>
	<li>使用 hclust 函数选择平均连接层次聚类方法进行层次聚类。执行命令行： 
	&#34;geneTree = hclust(as.dist(dissTOM), method = &#34;average&#34;);&#34; 
	点击 Run。</li>
	<li>使用 cutreeDynamic 函数执行动态分支剪切，并将 minClusterSize 参数设置为 30，以获得模块识别结果。执行命令行： 
	&#34;dynamicMods = cutreeDynamic(dendro = geneTree, distM = dissTOM, deepSplit = 2, pamRespectsDendro = FALSE, minClusterSize = minModuleSize);&#34; 
	点击 Run。 
	注意：最小模块大小不得低于 30。</li>
	<li>通过 moduleEigengenes 函数计算每个 OTU 模块的模块特征值。执行命令行： 
	&#34;MEList = moduleEigengenes(datExpr0, colors = dynamicColors) 
	MEs = MEList$eigengenes&#34; 
	点击 Run。 
	注意：模块特征值代表该模块中 OTU 的整体表达水平，它并非某个特定 OTU，而是通过奇异网络值分解获得的每个聚类的第一主成分。</li>
	<li>基于模块特征值的相关系数执行聚类函数。使用 mergeCloseModules 函数合并相关性低于 0.25 的模块。执行命令行： 
	&#34;merge = mergeCloseModules(datExpr0, dynamicColors, cutHeight = MEDissThres, verbose = 3)&#34; 
	点击 Run。</li>
	<li>最后，使用 plotDendroAndColors 函数进行可视化，以获得每个共表达网络的模块分配示意图（图 2）。使用 table 函数从模块分配表中提取每个 OTU 对应的模块归属信息。执行命令行： 
	&#34;plotDendroAndColors(geneTree, mergedColors, &#34;Merged dynamic&#34;,dendroLabels = FALSE, 
	&#160; &#160;hang = 0.03,addGuide = TRUE, guideHang = 0.05, 
	&#160; &#160;main = &#34;ES_Gene dendrogram and module colors&#34;)&#34; 
	点击 Run。 
	&#8203;注意：在共表达网络的模块分配图中，不同颜色代表不同的模块，灰色代表无法归入任何模块的 OTU。灰色模块中 OTU 数量越多，表明表达矩阵前期预处理质量较差。</li>
</ol>
4. 模块比较
注意：该方法可用于比较两个生态微生物群落的网络模块。本文中，比较内生环境与根面、内生环境与根际、根际与根面之间的微生物网络模块差异。
<ol>
	<li>保存性检验
	<ol>
		<li>加载前几步中保存的两组数据集的参数和结果。</li>
		<li>将一组微生物数据的网络模块划分结果设为参考组，另一组设为测试组。</li>
		<li>使用 modulePreservation 函数计算保守性统计参数 Z_summary 和 medianRank 的值。执行命令行： 
		&#34;system.time({mp=modulePreservation(multiExpr, 
		multiColor,referenceNetworks=1, 
		nPermutation=100, randomSeed=1,quickCor=0,verbose=3)})&#34; 
		点击 Run。 
		注意：该结果可用于量化模块之间的保守性。Z_summary&#62;10 表示两个模块高度保守，而 Z_summary&#60;2 表示模块未被保守。medianRank 通过排序反映模块的相对保守性，medianRank 值越高表示模块越未被保守。（关键命令行参见 Supplement S3。）</li>
		<li>使用 plot 函数可视化结果（图 3），获取 Z_summary 和 medianRank 参数值（表 1）。 
		注意：同时满足 Z_summary 值小于 2 且 medianRank 值处于最高水平的网络模块，是两个生态微生物群落中最未被保守的模块。</li>
		<li>根据上述两个统计参数的结果，确定两个网络中最为未被保守的模块。</li>
	</ol>
	</li>
	<li>模块成员关系的相关性分析
	<ol>
		<li>将两个网络的模块划分结果分别设为参考组和测试组。 
		注意：设置需与保存性检验一致。</li>
		<li>使用 corPvalueStudent 函数提取若干候选模块中每个 OTU 的 kME（模块成员度）值。 
		执行命令行： 
		&#34;Pvalue = as.data.frame(corPvalueStudent(as.matrix 
		(ModuleMembership), Samples))&#34; 
		点击 Run。 
		注意：kME 表示模块成员度，ME 表示模块特征向量（module eigen），代表模块中 OTU 表达的整体水平。kME 是每个 OTU 与 ME 之间的相关系数，通过 OTU 的 kME 值可量化其在网络中的重要性。（关键命令行参见 Supplement S4。）</li>
		<li>随后，使用 verboseScatterplot 函数计算两个网络中对应 OTU 的 kME 值之间的相关系数，并绘制相关性分析图（图 4）。 
		执行命令行： 
		&#34;verboseScatterplot(abs(TModuleMembership 
		[TmoduleGenes, Tcolumn]), 
		&#160; &#160;abs(NModuleMembership[NmoduleGenes, Ncolumn]), 
		&#160; &#160;xlab = paste(&#34;kME in&#34;, &#34;ES&#34;), 
		&#160; &#160;ylab = paste(&#34;kME in&#34;, &#34;RP&#34;), 
		&#160; &#160;main = paste(&#34;lightyellow&#34;), 
		&#160; &#160;cex.main = 1.7, cex.lab = 1.6, cex.axis = 1.6, col =&#160; &#160; &#160; &#160; &#160; modulecolor)&#34; 
		点击 Run。</li>
		<li>选择两个网络中 OTU 的 kME 值相关系数最小的模块，认为该模块在两个网络间的差异最大。</li>
	</ol>
	</li>
</ol>
5. 微生物差异网络模块的分析
<ol>
	<li>通过统计分析差异最大模块的OTU序列集，获得优势细菌门的数据。 
	注意：差异最大模块的OTU序列集按门水平分类进行汇总，优势细菌门所占比例超过10%。</li>
	<li>然后，使用exportNetworkToCytoscape函数获取最大差异模块中OTU相互作用关系信息的文件。 
	执行命令行： 
	&#34;cyt = exportNetworkToCytoscape(modTOM, 
	edgeFile = paste(&#34;NEW-ES_CytoscapeInput-edges-&#34;, modules , &#34;.txt&#34;, sep=""), 
	nodeFile = paste(&#34;NEW-ES_CytoscapeInput-nodes-&#34;, modules, &#34;.txt&#34;, sep=""), 
	weighted = TRUE,threshold = 0.5, nodeNames = modProbes, 
	altNodeNames = modGenes, nodeAttr = moduleColors[inModule])&#34; 
	点击 运行。</li>
	<li>将文件导入Cytoscape，设置阈值为0.5，并根据需要调整其他参数。</li>
	<li>构建差异微生物的共现网络（图5）。</li>
	<li>获得在网络中具有最重要调控作用的核心属的信息。 
	注意：根据OTU的kME值，可定义核心属。</li>
	<li>最后，评估核心属的功能，并分析其对整个差异网络的影响。</li>
</ol>

<ol><li>Philippot, L., Raaijmakers, J. M., Lemanceau, P., vander Putten, W. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Going+back+to+the+roots%3A+the+microbial+ecology+of+the+rhizosphere%5BTitle%5D)+AND+%22Nature+Reviews+Microbiology%22%5BJournal%5D)">Going back to the roots: the microbial ecology of the rhizosphere.</a> Nature Reviews Microbiology. 11, 789-799 (2013).</li>
<li>Fierer, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Embracing+the+unknown%3A+disentangling+the+complexities+of+the+soil+microbiome%5BTitle%5D)+AND+%22Nature+Review+Microbiology%22%5BJournal%5D)">Embracing the unknown: disentangling the complexities of the soil microbiome.</a> Nature Review Microbiology. 15 (10), 579-590 (2017).</li>
<li>Jin, J., Wang, G. H., Liu, X. B., Liu, J. D., Chen, X. L., Herbert, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Temporal+and+spatial+dynamics+of+bacterial+community+in+the+rhizosphere+of+soybean+genotypes+grown+in+a+black+soil%5BTitle%5D)+AND+%22Pedosphere%22%5BJournal%5D)">Temporal and spatial dynamics of bacterial community in the rhizosphere of soybean genotypes grown in a black soil.</a> Pedosphere. 19 (6), 808-816 (2009).</li>
<li>Ma, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genetic+correlation+network+prediction+of+forest+soil+microbial+functional+organization%5BTitle%5D)+AND+%22ISME+J%22%5BJournal%5D)">Genetic correlation network prediction of forest soil microbial functional organization.</a> ISME J. 12 (10), 2492-2505 (2018).</li>
<li>de Vries, F. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Soil+bacterial+networks+are+less+stable+under+drought+than+fungal+networks%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">Soil bacterial networks are less stable under drought than fungal networks.</a> Nature Communications. 9 (1), 3033(2018).</li>
<li>Colin, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Correlating+transcriptional+networks+to+breast+cancer+survival%3A+a+large-scale+coexpression+analysis%5BTitle%5D)+AND+%22Carcinogenesis%22%5BJournal%5D)">Correlating transcriptional networks to breast cancer survival: a large-scale coexpression analysis.</a> Carcinogenesis. (10), 2300-2308 (2013).</li>
<li>Ma, B., Zhao, K., Lv, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genetic+correlation+network+prediction+of+forest+soil+microbial+functional+organization%5BTitle%5D)+AND+%22ISME+J%22%5BJournal%5D)">Genetic correlation network prediction of forest soil microbial functional organization.</a> ISME J. 12, 2492-2505 (2018).</li>
<li>Zhang, B., Horvath, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+general+framework+for+weighted+gene+co-expression+network+analysis%5BTitle%5D)+AND+%22Statistical+applications+in+genetics+and+molecular+biology%22%5BJournal%5D)">A general framework for weighted gene co-expression network analysis.</a> Statistical applications in genetics and molecular biology. 4 (1), (2005).</li>
<li>Edwards, J. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Compositional+shifts+in+root-associated+bacterial+and+archaeal+microbiota+track+the+plant+life+cycle+in+field-grown+rice%5BTitle%5D)+AND+%22PLoS+Biology%22%5BJournal%5D)">Compositional shifts in root-associated bacterial and archaeal microbiota track the plant life cycle in field-grown rice.</a> PLoS Biology. 16 (2), 2003862(2018).</li>
<li>Bashan, Y., De-Bashan, L. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((How+the+Plant+Growth-Promoting+Bacterium+Azospirillum+Promotes+Plant+Growth-A+Critical+Assessment%5BTitle%5D)+AND+%22Advances+in+Agronomy%22%5BJournal%5D)">How the Plant Growth-Promoting Bacterium Azospirillum Promotes Plant Growth-A Critical Assessment.</a> Advances in Agronomy. 108, 77-136 (2010).</li>
<li>Lovley, D. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Geobacter%3A+The+Microbe+Electric%27s+Physiology%2C+Ecology%2C+and+Practical+Applications%5BTitle%5D)+AND+%22Advances+in+Microbial+Physiology%22%5BJournal%5D)">Geobacter: The Microbe Electric's Physiology, Ecology, and Practical Applications.</a> Advances in Microbial Physiology. 59, 1(2011).</li>
<li>Langfelder, P., Horvath, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((WGCNA%3A+an+R+package+for+weighted+correlation+network+analysis%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">WGCNA: an R package for weighted correlation network analysis.</a> BMC Bioinformatics. 9, 559(2008).</li>
<li>Zhang, B., Horvath, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+General+Framework+for+Weighted+Gene+Co-expression+Network+Analysis%5BTitle%5D)+AND+%22Statistical+Applications+in+Genetics+and+Molecular+Biology%22%5BJournal%5D)">A General Framework for Weighted Gene Co-expression Network Analysis.</a> Statistical Applications in Genetics and Molecular Biology. 4 (1), 17(2005).</li>
<li>Horvath, S., Dong, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Geometric+interpretation+of+Gene+Co-expression+Network+Analysis%5BTitle%5D)+AND+%22PLoS+Computational+Biology%22%5BJournal%5D)">Geometric interpretation of Gene Co-expression Network Analysis.</a> PLoS Computational Biology. 4 (8), 1000117(2008).</li>
<li>Langfelder, P., Horvath, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Eigengene+networks+for+studying+the+relationships+between+co-expression+modules%5BTitle%5D)+AND+%22BMC+Systems+Biology%22%5BJournal%5D)">Eigengene networks for studying the relationships between co-expression modules.</a> BMC Systems Biology. 1, 54(2007).</li>
<li>Horvath, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Analysis+of+Oncogenic+Signaling+Networks+in+Glioblastoma+Identifies+ASPM+as+a+Novel+Molecular+Target%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Analysis of Oncogenic Signaling Networks in Glioblastoma Identifies ASPM as a Novel Molecular Target.</a> Proceedings of the National Academy of Sciences of the United States of America. 103 (46), 17402-17407 (2006).</li>
<li>Carlson, M. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((and+Sequence+Conservation%3A+Predictions+from+Modular+Yeast+Co-expression+Networks%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">and Sequence Conservation: Predictions from Modular Yeast Co-expression Networks.</a> BMC Genomics. 7 (1), 40(2006).</li>
<li>Fuller, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Weighted+Gene+Co-expression+Network+Analysis+Strategies+Applied+to+Mouse+Weight%5BTitle%5D)+AND+%22Mammalian+Genome%22%5BJournal%5D)">Weighted Gene Co-expression Network Analysis Strategies Applied to Mouse Weight.</a> Mammalian Genome. 6 (18), 463-472 (2007).</li>
<li>Yin, L., Wang, Y., Lin, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Explorative+analysis+of+the+gene+expression+profile+during+liver+regeneration+of+mouse%3A+a+microarray-based+study%5BJ%5D%5BTitle%5D)+AND+%22Artificial+Cells+Nanomedicine+%26+Biotechnology%22%5BJournal%5D)">Explorative analysis of the gene expression profile during liver regeneration of mouse: a microarray-based study[J].</a> Artificial Cells Nanomedicine & Biotechnology. 47 (1), 1113-1121 (2019).</li>
<li>Oldham, M., Horvath, S., Geschwind, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Conservation+and+Evolution+of+Gene+Co-expression+Networks+in+Human+and+Chimpanzee+Brains%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Conservation and Evolution of Gene Co-expression Networks in Human and Chimpanzee Brains.</a> Proceedings of the National Academy of Sciences of the United States of America. 103 (47), 17973-17978 (2006).</li>
<li>Oldham, M. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+organization+of+the+transcriptome+in+human+brain%5BTitle%5D)+AND+%22Nature+Neuroscience%22%5BJournal%5D)">Functional organization of the transcriptome in human brain.</a> Nature Neuroscience. 11 (11), 1271-1282 (2008).</li>
<li>Keller, M. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+gene+expression+network+model+of+type+2+diabetes+links+cell+cycle+regulation+in+islets+with+diabetes+susceptibility%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">A gene expression network model of type 2 diabetes links cell cycle regulation in islets with diabetes susceptibility.</a> Genome Research. 18 (5), 706-716 (2008).</li>
<li>Weston, D., Gunter, L., Rogers, A., Wullschleger, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Connecting+genes%2C+coexpression+modules%2C+and+molecular+signatures+to+environmental+stress+phenotypes+in+plants%5BTitle%5D)+AND+%22BMC+Systems+Biology%22%5BJournal%5D)">Connecting genes, coexpression modules, and molecular signatures to environmental stress phenotypes in plants.</a> BMC Systems Biology. 2 (1), 16(2008).</li>
<li>Jorda´n, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Keystone+species+and+food+webs%5BTitle%5D)+AND+%22Biological+Sciences%22%5BJournal%5D)">Keystone species and food webs.</a> Biological Sciences. 364, 1733-1741 (2009).</li>
<li>Backhed, F., Ley, R. E., Sonnenburg, J. L., Peterson, D. A., Gordon, J. I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Host-Bacterial+Mutualism+in+the+Human+Intestine%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Host-Bacterial Mutualism in the Human Intestine.</a> Science. 307, 1915-1920 (2009).</li>
</ol>

作者无任何利益冲突需要披露。

相关性网络在生物信息学应用中正被越来越多地使用。WGCNA 是一种系统生物学方法，用于描述生物系统中各个元素之间关系的分析12。早期关于 WGCNA 的研究使用了 R 软件包13,14,15。该软件包包含用于网络构建、模块检测、拓扑性质计算、数据模拟、可视化以及与外部软件交互的功能。WGCNA 已被广泛应用于分析脑癌16、酵母细胞周期17、小鼠遗传学18,19、灵长类动物脑组织20,21、糖尿病22以及植物<sup class='x...

微生物组研究对于理解与调控生态系统过程具有重要意义1,4。微生物种群通过相互作用的生态网络彼此关联，这些网络的特征可影响微生物对环境变化的响应3,4。此外，这些网络的特性影响着微生物群落的稳定性，并与土壤功能密切相关5。加权基因共表达网络分析（WGCNA）目前已广泛应用于基因与微生物群落关系的研究6。以往的研究主要关注不同基因或种群网络与外部环境之间的关联7。然而，在不同环境条件下由微生物种群形成的共现网络差异却鲜有研究。本文所呈现的研究旨在提供深入见解和具体细节，介绍如何快速实施WGCNA算法，以构建在不同环境条件下采集的微生物组样本的共现网络。基于分析结果，我们评估了种群的组成及差异，并进一步探讨了不同微生物种群之间的关系。本研究应用了如下加权相关性网络算法的基本流程8：首先，通过计算操作分类单元（OTU）表达谱之间的皮尔逊相关系数构建相似性矩阵；随后，依据无标度拓扑标准选择邻接函数参数（幂函数或S型函数），将相似性矩阵转化为邻接矩阵，每个共现网络对应一个邻接矩阵；接着，采用基于TOM的不相似性指标结合平均连接层次聚类方法，将表达模式一致的OTU聚类为模块；进一步计算保守性统计量与相关参数分析模块之间的关系，最终识别出模块中的核心OTU。这些方法特别适用于分析不同环境条件下多种微生物种群间网络结构的差异。在本文中，我们详细描述了共表达网络构建的方法、模块间差异的分析过程，并简要概述了用于识别不同模块网络中核心物种的操作步骤。

<table><tbody><tr><td>R</td><td>The University of Auckland</td><td>version 4.0.2</td><td>R 是一个用于统计计算和图形的自由软件环境，可在多种 UNIX 平台、Windows 和 MacOS 系统上编译和运行。</td></tr><tr><td>RStdio</td><td>JJ Allaire</td><td>version 1.4.1103</td><td>RStudio 集成开发环境是一组集成工具，旨在帮助用户更高效地使用 R 和 Python。</td></tr><tr><td>Cytoscape</td><td></td><td>version 3.7.1</td><td>Cytoscape 是一个开源软件平台，用于可视化复杂网络，并将其与任意类型的属性数据进行整合。</td></tr><tr><td>NCBI database</td><td></td><td></td><td>美国国家生物技术信息中心通过提供生物医学和基因组信息，推动科学与健康事业的发展。</td></tr></tbody></table>

基于加权相关网络的不同生境下根系微生物群落的分化

Divergence of Root Microbiota in Different Habitats based on Weighted Correlation Networks

根部微生物组在植物生长和环境适应中发挥着重要作用。网络分析是研究群落的重要工具，能够有效探究不同环境中不同微生物物种之间的相互作用关系或共现模式。本文旨在详细介绍如何利用加权相关网络算法，分析因不同生态环境而在微生物群落中可能出现的各类共现网络。本实验的所有分析均在WGCNA软件包中进行。WGCNA是一个用于加权相关网络分析的R语言软件包。用于演示这些方法的实验数据来自NCBI（国家生物技术信息中心）数据库中的水稻（Oryza sativa）根系三个生态位的微生物群落数据。我们采用加权相关网络算法，分别构建了这三个生态位中微生物群落的共丰度网络，并鉴定了内生圈、根表层和根际土壤之间的差异共丰度网络。此外，通过“WGCNA”软件包获得了网络中的核心属，这些核心属在网络功能中发挥重要的调控作用。这些方法使研究人员能够分析微生物网络对环境扰动的响应，并验证不同的微生物生态响应理论。结果表明，在水稻的内生圈、根表层和根际土壤中均鉴定出显著差异的微生物网络。

本文中的代表性结果下载自NCBI数据库中2014年加利福尼亚Abaker水稻根系微生物组数据（PRJNA386367）9。该数据包括在淹水水稻田中生长14周的水稻植株的根际、根表和内生圈微生物组样本。我们采用WGCNA算法选择满足接近无尺度网络的三个网络的幂值（图1），并构建了三个共表达网络（图2）。在内生圈、根表和根际土壤微生物共表达网络中，分别鉴定了23、22和21个模块。这些结果表明，这三个生态位中的微生物相互作用网络数量基本相等。
我们进一步比较了内生圈、根表层和根际土壤中微生物网络模块的差异。获得了三个生态位组模块的以下保真性检验结果。在根际土壤与根表层之间存在三个极不保守的模块（图3a，表1）；在根际土壤与内生圈之间存在九个极不保守的模块（图3b<str...

采用网络分析方法评估不同生态微生物群落（如土壤、水体和根际）之间的关联。本文介绍了一种利用WGCNA算法分析因不同生态环境而产生的微生物群落中可能存在的各种共现网络的实验方案。

根际微生物组在植物生长和环境适应中发挥着重要作用。网络分析是研究群落的重要工具，……

该方法可有效探究不同环境中多种微生物物种之间的相互作用关系或共现网络。本文详细介绍了如何利用WGCNA算法构建微生物群落的共现网络。此外，基于分析结果，该方法可评估不同微生物群落之间微生物关系及组成的差异性。以下是该方法的基本流程。从NCBI数据库下载微生物群的组成和丰度数据，或使用您样本的测序数据。首先，打开RStudio软件并安装WGCNA包。然后，加载数据并使用“好样本基因”功能检查数据的正确性。检查是否存在离群值，并保存符合我们要求的样本。当检查结果为真时，继续下一步。保存结果。使用“剔除阈值”功能计算数据在不同幂值下的无标度指数 R 平方。可视化结果。当无标度指数接近1时，网络结构更接近于无标度网络。当无标度指数R平方大于0.9时，选择该幂律值。最后，使用相同的方法分析其余的微生物组数据。首先，使用邻接函数构建符号共现网络。此外，使用TOM相似性函数构建拓扑重叠网络。其次，使用Hclust函数进行层次聚类，并绘制所得的聚类树。第三，使用 cutree 动态函数进行动态分支切割，并通过最小聚类大小参数设置最小模块尺寸。最小模块尺寸通常设为大于 30。第四，计算每个模块的微生物特征。根据相关系数进行层次聚类，并将高度小于0.25的模块合并，以获得各模块的分布情况。第五步，使用plot dendro和colors函数可视化结果，得到共现网络模块的分配展示图。然后，重命名模块颜色，并构建与这些颜色对应的数字标签。将它们保存下来，以供后续部分使用。最后，对其他数据集重复上述过程。在此部分，比较并分析两组数据，并进行保存性测试。首先，加载之前步骤中保存的两个数据集的参数和结果。然后，将其中一个数据集的模块结果设为对照组，另一个设为测试组，并执行模块比较。接下来，计算保守性值、统计参数、Z 汇总值和中位秩，以量化模块之间的保守性。最后，对结果进行可视化。确定同时满足 Z 汇总值小于 2 且中位秩值位于前列的网络模块。该模块是两种微生物组数据中保守性最低的模块。对模块成员关系进行相关性分析。将两个网络的模块分配结果分别设为参考组和测试组。设置需与保存性检验相同。首先，根据保存性检验的结果，在多个候选模块中计算每个OTU的KME值。以黄色模块为例。计算两个黄色模块中KME值的相关系数，然后绘制模块内KME值的相关性分析图。最后，根据图中的相关系数判断该模块在两个数据集中的保守性，选择相关系数最小的模块。首先，使用“导出网络至 Cytoscape”功能，将网络导出为 Cytoscape 可读的边和节点列表文件。然后，将该文件导入 Cytoscape。将阈值设置为 0.5，并根据需要调整其他参数。最后，获得不同微生物的共现网络。本文采用WGCNA算法分析水稻根系三个生态位之间的差异。选择满足三个网络均接近无标度网络的幂值。在内生层、根表层和根际微生物共现网络中，分别鉴定出23、22和21个模块。利用保守性检验和相关性分析，找出水稻根系各两个生态位之间极不保守的模块。使用不同颜色代表不同的微生物门，构建这三个模块的共现网络。变形菌门、放线菌门、拟杆菌门、厚壁菌门和疣微菌门主导了这三个不同的微生物网络。此外，17个核心属主要调控这些网络。观看本视频后，您应能够充分理解如何执行一系列步骤，利用WGCNA算法分析因不同生态环境而在微生物群落中可能出现的不同共现网络。

divergence-root-microbiota-different-habitats-based-on-weighted

Research

科研

JoVE Journal

Biology

生物学

基于加权相关网络的不同栖息地根微生物群的差异

本文介绍了一种获取土壤、根际和根内微生物组扩增子序列数据的实验方案。该信息可用于研究与植物相关的微生物群落的组成和多样性，适用于多种植物物种。...

探索根系微生物组：从土壤、根际和根内圈提取细菌群落数据

本文详细描述了从田间挖掘植物根系以及将样品处理为内生圈、根际和土壤的方法，包括DNA提取和数据分析方法。本文旨在使其他实验室能够采用这些技术来研究土壤、内生圈和根际微生物组。...

多年生草本植物实验中土壤、根际及根部微生物群落的分离与分析

流式细胞术分析已被证明在研究纯培养物和监测微生物群落动态方面具有重要价值。我们分别提供了从取样到数据分析的三个完整工作流程，适用于透明培养基中的纯培养物和复杂群落，以及具有挑战性的基质环境。...

表征微生物组动态——从纯培养到自然群落的基于流式细胞术的工作流程

本文介绍了一种在体外培养结肠肠道微生物群的实验方案,...

应用先进的体外培养技术研究人类肠道微生物群

群落中的大多数植物可能都通过丛枝菌根（AM）真菌相互连接，但关于这些真菌对植物相互作用的调控作用，以往研究主要通过比较植物在有无菌根条件下的生长情况来开展。本文介绍一种在菌根植物之间调控共有菌根网络的方法，以研究其对植物相互作用的影响。...

利用旋转土芯法研究通过共同菌根网络连接的植物间相互作用

利用两步PCR和下一代16S rRNA基因测序进行微生物群分析

本实验方案描述了两种物种——Zea mays 和 Festuca rubra 根部丛枝菌根定殖模式与策略的评估方法。采用 MycoPatt...

菌根图谱作为探索红羊茅（Festuca rubra）和玉米（Zea mays）根部定殖模式与真菌策略的工具

堆肥微宇宙实验将自然界中存在的微生物多样性引入实验室，以促进对秀丽隐杆线虫（Caenorhabditis...

堆肥微宇宙作为微生物组研究中微生物多样性高且接近自然环境的实验体系 秀丽隐杆线虫

我们提出了一种优化的扫描电子显微镜实验方案，用于观察包含植物和真菌生物量以及微生物群和生物膜的高度异质且脆弱的样品。该方案可用于描述微生物群空间结构的维度特征。...

切叶蚁菌圃的微生物群落：通过扫描电子显微镜可视化微生物景观

本方案提供了在微电极阵列（MEA）上建立、维护、记录和电刺激神经元培养物所需的全部信息。In vitro...

如何在微电极阵列（MEAs）上培养、记录和刺激神经元网络

学习使用WGCNA算法构建水稻根系不同生态位微生物共现网络，识别核心属和非保守模块。该方法可揭示根际土壤、根表和根内微生物群落的差异及其生态适应机制。

基于加权相关网络的不同生境下根系微生物群落的分化

探索更多视频