为进一步提高《微生物组实验手册》稿件质量，本项目新增大众评审环节。文章在通过同行评审后，采用公众号推送方式分享全文，任何人均可在线提交修改意见。公众号格式显示略有问题，建议电脑端点击文末阅读原文下载PDF审稿。在线文档(https://kdocs.cn/l/cL8RRqHIL)大众评审页面登记姓名、单位和行号索引的修改建议。修改意见的征集截止时间为推文发布后的72小时，文章将会结合有建设性的修改意见进一步修改后获得DOI在线发表，同时根据贡献程度列为审稿人或致谢。感谢广大同行提出宝贵意见。

土壤微生物响应环境变化的系统发育保守性和环境阈值

Phylogenetic conservation and environmental threshold of soil microbes in response to environmental changes

陈贝贝，焦硕^*，韦革宏

西北农林科技大学生命科学学院，旱区作物逆境生物学国家重要实验室，陕西省农业与环境微生物重点实验室，陕西杨凌

*通讯作者邮箱：shuojiao@nwafu.edu.cn

摘要：土壤微生物在促进土壤养分循环和调节植物矿质营养等方面起着重要的生态作用。但由于其对环境变化的敏感性，外界环境变化通常会导致其多样性和群落结构的变化。对于土壤微生物响应外界环境变化的系统发育保守性和环境阈值知之甚少。本文以中国东部农田土壤微生物群落数据为示例，基于微生物群落对环境梯度的响应，计算Blomberg’s K值和环境阈值，以评估系统发育信号强度和环境适应性宽度，给出了计算微生物系统发育保守性和环境阈值的标准化流程，为土壤微生物响应外界环境变化的研究提供了新思路。

关键词：土壤微生物，系统发育保守性，环境阈值

仪器设备

1.普通Windows系统个人电脑，内存8G，需求内存空间约500M

软件

1.R (v3.6.3)

注：本教程是基于已经在个人电脑上安装好的相关软件和依赖包进行的。

实验步骤

1.数据准备

本分析用到两个数据，即土壤微生物群落数据，环境因子数据(与采样点一一对应)。

2.数据准备

2.1导入已抽平的微生物群落数据及环境因子数据（图1和图2）

Micro<-< span="">read.csv("otu_data.csv", header = T,row.names=1)

Env<-read.csv("env.csv",header=t,row.names=1)< span="">

图1. 微生物群落数据集

图2. 环境因子数据集

3.系统发育信号的计算

3.1微生物群落特征数据集的获取

为了获得潜在性状信息，通过微生物类群与环境变量之间Spearman相关来定义每个OTU的生态偏好。例如，OTU与pH的正负相关被定义为acid-preferred或alkaline-preferred，相应的在性状数据集中：正相关为1，负相关为-1，没有显著相关的为0。

相应的计算过程如下：

library(Hmisc) #所需依赖包：lattice、survival、ggplot2和Formula

Micro<-< span="">Micro/11730 #11730为每个样本的最小抽平序列数

Micro<-t(micro)< span="">

cor<-rcorr(as.matrix(cbind(micro,env)),type="spear")< span="">

r<-cor$r< span="">

q<-cor$p< span="">

q[q<0.05]<-1< span="">

q[q>0.05&q<1]<-0< span="">

g<-r*q< span="">

heat<-g[1:ncol(micro),-(1:ncol(micro))]< span="">

traits<-replace(heat, heat<0,-1)< span="">

traits<-replace(traits, traits="">0,1)

write.csv(traits,"traits.csv") #写出相应的性状数据集

图3. 微生物群落性状数据集

3.2Blomberg’s K的计算

library(picante)  #所需依赖包：vegan、permute和nlme

phy<-read.tree("Micro.tre") #读取微生物类群对应的系统发育树

phy<-multi2di(phy) #无根树转为有根树

traits <- traits[phy$tip.label, ]

traits<-as.matrix(traits)< span="">

Micro_Blombergk<-multiPhylosignal(traits,phy) #计算Blomberg’s K

write.csv(Micro_Blombergk," Micro_Blombergk.csv")

写出计算结果（图4）

图4. 布朗K值数据集

K值描述了一个分类单元与系统发育的关系，正如布朗运动的系统发育信号强度度量所期望的那样。K值接近于0表示进化是随机的或收敛的，而K值大于1则表示强烈的系统发育信号和性状保守性。

3.3绘制系统发育信号柱状图

data<-read.csv("micro_blombergk.csv",header=t,row.names=1)< span="">

data<-t(data)< span="">

data<-data[1,]< span="">

data<-as.matrix(data)< span="">

barplot(data,axes=T,width=2,axis.lty=1,ylab="Blomberg’s K ",

legend=rownames(data),beside=TRUE,col=c("#366CA8"))

图5. 系统发育信号柱状图

4.环境阈值的计算

4.1对输入的微生物群落数据集Micro进行筛选，标准为：每个OTU至少在3个样本中存在，获得的数据集为Micro1

Micro1<- micro="">= 1)>= 3,]

Micro1<- Micro1/11730 #11730为每个样本的最小抽平序列数

Micro1<-t(micro1)

library(TITAN2)

pH_titan<-titan(Env[,2], Micro1,numPerm=250,

boot=TRUE,nBoot=500,imax=FALSE,ivTot=FALSE,

pur.cut=0.95,rel.cut=0.95,memory=TRUE)

plot_sumz(pH_titan,filter=TRUE,xlab=expression("pH"),pch1=20,pch2=20, col1="#fdbb2d",col2="#1E9600")

图6. 微生物群落负响应种（Z^—）和正响应种（Z⁺）指示总分沿pH梯度突变点的响应曲线

plot_taxa(pH_titan, z.med=F,leg.x=.8, leg.y=5, xlab=expression("pH"), cex.taxa=.9,cex.axis=1.35,cex=1,fil1="white",fil2="white",pch1=20, pch2=20,

col1="#fdbb2d",col2="#1E9600",pur.cut=0.95, pval.cut=0.05, rel.cut=0.95)

图7. 微生物群落沿pH梯度负响应种和正响应种分布图

4.2重复第4.1条步骤，分别探究微生物群落沿其他环境梯度的环境阈值。

致谢

本方案主要来源于课题组先前发表的相关文章(Jiao et al.,2020; Zhang et al.,2020)。相关研究得到了国家自然科学基金(41807030)、国家重点研发计划项目(2016YFD0200306) 和中国博士后科学基金(2018M630041) 等项目的资助。

参考文献

1.Baker, M. E. and R. S. King (2010). A new method for detecting and interpreting biodiversity and ecological community thresholds. Methods Ecol Evol 1: 25-37

2.Isobe, K., S. D. Allison, B. Khalili, A. C. Martiny and J. B. H. Martiny (2019). Phylogenetic conservation of bacterial responses to soil nitrogen addition across continents. Nat Commun 10(1): 2499.

3.Jiao, S. and Y. Lu (2020). Abundant fungi adapt to broader environmental gradients than rare fungi in agricultural fields. Glob Chang Biol 26(8).

4.Oliverio, A. M., M. A. Bradford and N. Fierer (2016). Identifying the microbial taxa that consistently respond tosoil warming across time and space. Glob Chang Biol. 23(5)

5.Zhang, Z., J. Zhang and S. Jiao (2020). Fungi show broader environmental thresholds in wet than dry agricultural soils with distinct biogeographic patterns. Sci Total Environ 750: 141761.

猜你喜欢

10000+：菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature Cell专刊 肠道指挥大脑

系列教程：微生物组入门 Biostar 微生物组  宏基因组

专业技能：学术图表 高分文章 生信宝典 不可或缺的人

一文读懂：宏基因组 寄生虫益处 进化树

必备技能：提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

扩增子分析：图表解读 分析流程 统计绘图

16S功能预测   PICRUSt  FAPROTAX  Bugbase Tax4Fun

在线工具：16S预测培养基 生信绘图

科研经验：云笔记  云协作 公众号

编程模板: Shell  R Perl

生物科普:  肠道细菌 人体上的生命 生命大跃进  细胞暗战 人体奥秘

写在后面

为鼓励读者交流、快速解决科研困难，我们建立了“宏基因组”专业讨论群，目前己有国内外5000+ 一线科研人员加入。参与讨论，获得专业解答，欢迎分享此文至朋友圈，并扫码加主编好友带你入群，务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份，另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助，首先阅读《如何优雅的提问》学习解决问题思路，仍未解决群内讨论，问题不私聊，帮助同行。

学习16S扩增子、宏基因组科研思路和分析实战，关注“宏基因组”

点击阅读原文下载PDF审稿，或浏览器直接访问下载链接：http://210.75.224.110/github/MicrobiomeProtocol/04Review/210429/2104057ShuoJiao1000264/protocol2104057R1.pdf