《环境微生物群落多样性分析.docx》由会员分享,可在线阅读,更多相关《环境微生物群落多样性分析.docx(22页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、环境微生物群落多样性分析微生物群落多样性的根本概念环境中微生物的群落构造及多样性与微生物的功能及代谢机理是微生物生态学的研究热点。长期以来,由于受到技术限制,对微生物群落构造与多样性的认识还不全面,对微生物功能及代谢机理方面了解的也很少。但随着高通量测序、基因芯片等新技术的不断更新,微生物分子生态学的研究方法与研究途径也在不断变化。第二代高通量测序技术尤其是Roche 454高通量测序技术的成熟与普及,使我们能够对环境微生物进展深度测序,灵敏地探测出环境微生物群落构造随外界环境的改变而发生的极其微弱的变化,对于我们研究微生物及环境的关系、环境治理与微生物资源的利用以及人类医疗安康有着重要的理论
2、与现实意义。在国内,微生物多样性的研究涉及农业、土壤、林业、海洋、矿井、人体医学等诸多领域。以在医疗领域的应用为例,通过比拟正常与疾病状态下或疾病不同进程中人体微生物群落的构造与功能变化,可以对正常人群及某些疾病患者体内的微生物群体多样性进展比拟分析,研究获得人体微生物群落变化同疾病之间的关系;通过深度测序还可以快速地发现与检测常见病原及新发传染病病原微生物。研究方法进展环境微生物多样性的研究方法很多,从国内外目前采用的方法来看大致上包括以下四类:传统的微生物平板纯培养方法、微平板分析方法、磷脂脂肪酸法以及分子生物学方法等等。 近几年,随着分子生物学的开展,尤其是高通量测序技术的研发及应用,为
3、微生物分子生态学的研究策略注入了新的力量。 目前用于研究微生物多样性的分子生物学技术主要包括:DGGE/TGGE/TTGE、T-RFLP、SSCP、FISH、印记杂交、定量PCR、基因芯片等。DGGE等分子指纹图谱技术,在其实验结果中往往只含有数十条条带,只能反映出样品中少数优势菌的信息;另一方面,由于分辨率的误差,局部电泳条带中可能包含不只一种16SrDNA序列,因此要得悉电泳图谱中具体的菌种信息,还需对每一条带构建克隆文库,并筛选克隆进展测序,此实验操作相对繁琐;此外,采用这种方法无法对样品中的微生物做到绝对定量。生物芯片是通过固定在芯片上的探针来获得微生物多样性的信息,“只能验证,却无法
4、探索未知,此方法通过信号强弱判断微生物的丰度也不是非常的准确。 而近年来以454焦磷酸测序为代表的高通量测序技术凭借低本钱、高通量、流程自动化的优势为研究微生物群落构造提供了新的技术平台。Roche454高通量测序技术能同时对样品中的优势物种、稀有物种及一些未知的物种进展检测,获得样品中的微生物群落组成,并将其含量进展数字化。最近,美吉生物推出了新的测序平台MiSeq。MiSeq高通量测序平台集中了Roche454与IlluminaHiSeq2500的优点,不仅可实现对多样品的多个可变区同时测序,而且在测序速度与测序通量上都有进一步提升,目前此平台已在微生物多样性群落构造研究方面受到了广阔学者
5、的认可。 第二代高通量测序技术 产品优势 无需培养别离菌群: 直接从环境样本中扩增核糖体RNA高变区进展测序,解决了大局部菌株不可培养的难题。 客观复原菌群构造: 专业、成熟、稳定的样本制备流程,严格控制PCR循环数,客观复原样品本身的菌群构造及丰度比例。 痕量菌检测: 充分发挥高通量测序的大数据量优势,能检测出丰度低至万分之一的痕量菌。 效劳流程 样本要求 环境样品土壤:5-10g; 水体:2L水样或0.22m滤膜过滤; 粪便:3g; 黏膜:指甲大小; 植物内生菌:10-20g叶片;3-5g根系; 底泥:5-10g; 血液:10mL; 叶片:50-100g。 DNA浓度10ng/L总量500
6、ng的DNA,OD260/280介于1.8-2.0之间并确保DNA无降解。 PCR产物仅限Roche 454平台 PCR产物浓度5ng/L,总量100ng,OD260/280介于1.8-2.0之间并确保PCR产物无降解; PCR产物需经电泳切胶回收纯化; 送样管管口使用Parafilm封口膜密封; 样品保存期间切忌反复冻融,使用干冰运输。 生信分析1. 稀释性曲线Rarefaction Curve采用对测序序列进展随机抽样的方法,以抽到的序列数及它们所能代表OTU的数目构建曲线,即稀释性曲线。 当曲线趋于平坦时,说明测序数据量合理,更多的数据量对发现新OTU的边际奉献很小;反之那么说明继续测序
7、还可能产生较多新的OTU。 横轴:从某个样品中随机抽取的测序条数;Label 0.03 表示该分析是基于OTU 序列差异水平在0.03,即相似度为97% 的水平上进展运算的,客户可以选取其他不同的相似度水平。 纵轴:基于该测序条数能构建的OTU数量。 曲线解读:图1中每条曲线代表一个样品,用不同颜色标记; 随测序深度增加,被发现OTU的数量增加。当曲线趋于平缓时表示此时的测序数据量较为合理。 2. Shannon-Wiener 曲线反映样品中微生物多样性的指数,利用各样品的测序量在不同测序深度时的微生物多样性指数构建曲线,以此反映各样本在不同测序数量时的微生物多样性。 当曲线趋向平坦时,说明测
8、序数据量足够大,可以反映样品中绝大多数的微生物物种信息。 横轴:从某个样品中随机抽取的测序条数。 纵轴:Shannon-Wiener指数,用来估算群落多样性的上下。 Shannon 指数计算公式: 其中, Sobs=实际测量出的OTU数目; ni=含有i 条序列的OTU数目; N=所有的序列数。曲线解读: 图2每条曲线代表一个样品,用不同颜色标记,末端数字为实际测序条数; 起初曲线直线上升,是由于测序条数远缺乏覆盖样品导致; 数值升高直至平滑说明测序条数足以覆盖样品中的大局部微生物。 3.Rank-Abundance 曲线用于同时解释样品多样性的两个方面,即样品所含物种的丰富程度与均匀程度。
9、物种的丰富程度由曲线在横轴上的长度来反映,曲线越宽,表示物种的组成越丰富; 物种组成的均匀程度由曲线的形状来反映,曲线越平坦,表示物种组成的均匀程度越高。横轴:OTU 相对丰度含量等级降序排列。 纵轴:相对丰度比例。 曲线解读: 图3及图4中每条曲线对应一个样本参考右上角图标; 图3及图4中横坐标表示的是OTU物种丰度排列顺序,纵坐标对应的是OTU物种所占相对丰度比例图3为相对百分比例,图4为换算后Log值,曲线趋于水平那么表示样品中各物种所占比例相似;曲线整体斜率越大那么表示样品中各物种所占比例差异较大。 4. 样本群落组成分析:多样本柱状图/ 单样本饼状图根据分类学分析结果,可以得知一个或
10、多个样品在各分类水平上的物种组成比例情况,反映样品在不同分类学水平上的群落构造。 柱状图图5横轴:各样品的编号。纵轴:相对丰度比例。 图标解读:颜色对应此分类学水平下各物种名称,不同色块宽度表示不同物种相对丰度比例; 可以在不同分类学水平下作图分析。 饼状图图6 在某一分类学水平上,不同菌群所占的相对丰度比例。不同颜色代表不同的物种。 5. 样品OTU 分布Venn 图用于统计多个样品中共有或独有的OTU数目,可以比拟直观地表现各环境样品之间的OTU组成相似程度。 不同样品用不同颜色标记,各个数字代表了某个样品独有或几种样品共有的OTU数量,对应的OTU编号会以EXCEL表的形式在结题报告中呈
11、现。 分析要求单张分析图,样本分组至少两个,最多5个。 默认设置为97%相似度水平下以OTU为单位进展分析作图。 6.Heatmap图 用颜色变化来反映二维矩阵或表格中的数据信息,它可以直观地将数据值的大小以定义的颜色深浅表示出来。将高丰度与低丰度的物种分块聚集,通过颜色梯度及相似程度来反映多个样品在各分类水平上群落组成的相似性与差异性。 相对丰度比例: 热图图8中每小格代表其所在样品中某个OTU的相对丰度。以图8为例,红框高亮的小格所对应的信息为:样本R11-1Z中OTUOTU128的相对丰度比例大概为0.2%。 丰度比例计算公式BrayCurtis算法:其中,SA,i=表示A样品中第i个O
12、TU所含的序列数 SB,i=表示B样品中第i个OTU所含的序列数 样品间聚类关系树: 进化树表示在选用成图数据中,样本及样本间序列的进化关系差异关系。处于同一分支内的样品序列进化关系相近。 物种/OTU丰度相似性树: 丰度相似性树表示选用成图的数据中样品及样品中的OTU或序列在丰度上的相似程度。丰度最相近的会分配到同一分支上。 客户自定义分组:根据研究需求对菌群物种/OTU研究样本进展二级分组 二级物种/OTU分组:将下级分类学水平物种或OTU分配到对应的上级分类学水平,以不同颜色区分; 二级样品分组:根据研究需要,对样品进展人为的分组,以不同颜色区分。 7.主成分分析PCA(Principa
13、lComponentAnalysis) 在多元统计分析中,主成分分析是一种简化数据集的技术。主成分分析经常用于减少数据集的维数,同时保持数据集中对方差奉献最大的特征,从而有效地找出数据中最“主要的元素与构造,去除噪音与冗余,将原有的复杂数据降维,提醒隐藏在复杂数据背后的简单构造。 通过分析不同样品的OTU组成可以反映样品间的差异与距离,PCA运用方差分解,将多组数据的差异反映在二维坐标图上,坐标轴为能够最大程度反映方差的两个特征值。如样品组成越相似,反映在PCA图中的距离越近。 横轴与纵轴:以百分数的形式表达主成分主要影响程度。以图9为例,主成分1PC1与主成分2PC2是造成四组样品红色,蓝色
14、,黄色与绿色的两个最大差异特征,奉献率分别为41.1%与27.1%。 十字穿插线:在图9中作为0点基线存在,起到辅助分析的作用,本身没有意义。 图例解读: PCA分析图是基于每个样品中所含有的全部OTU完成的; 图9中每个点代表了一个样本;颜色那么代表不同的样品分组; 两点之间在横、纵坐标上的距离,代表了样品受主成分PC1或PC2影响下的相似性距离; 样本数量越多,该分析意义越大;反之样本数量过少,会产生个体差异,导致PCA分析成图后形成较大距离的分开,建议多组样品时,每组不少于5个,不分组时样品不少于10个; 图10中的圆圈为聚类分析结果,圆圈内的样品,其相似距离比拟接近。 8. RDA/
15、CCA 分析图 基于对应分析开展的一种排序方法,将对应分析及多元回归分析相结合,每一步计算均及环境因子进展回归,又称多元直接梯度分析。主要用来反映菌群及环境因子之间的关系。RDA是基于线性模型,CCA是基于单峰模型。分析可以检测环境因子、样品、菌群三者之间的关系或者两两之间的关系。 横轴与纵轴:RDA与CCA分析,模型不同,横纵坐标上的刻度为每个样品或者物种在及环境因子进展回归分析计算时产生的值,可以绘制于二维图形中。 图例解读: 冗余分析可以基于所有样品的OTU作图,也可以基于样品中优势物种作图; 箭头射线:图11中的箭头分别代表不同的环境因子即图中的碳酸氢根离子HCO3-,醋酸根离子AC-
16、等,图中的其它环境因子因研究不同代表的意义不同,因此不再赘述; 夹角:环境因子之间的夹角为锐角时表示两个环境因子之间呈正相关关系,钝角时呈负相关关系。环境因子的射线越长,说明该影响因子的影响程度越大; 图11中不同颜色的点表示不同组别的样品或者同一组别不同时期的样品,图中的拉丁文代表物种名称,可以将关注的优势物种也纳入图中; 环境因子数量要少于样本数量,同时在分析时,需要提供环境因子的数据,比方pH值,测定的温度值等。 9.单样品/多样品分类学系统组成树 根据NCBI提供的已有微生物物种的分类学信息数据库,将测序得到的物种丰度信息回归至数据库的分类学系统关系树中,从整个分类系统上全面了解样品中
17、所有微生物的进化关系与丰度差异。 单样品图图12:可以了解单样品中的序列在各个分类学水平上的分布情况。图例解读:图12中不同的层次反映不同的分类学水平; 分支处的圆面积说明了分布在该分类学水平,且无法继续往下级水平比对的序列数量,面积越大,说明此类序列越多; 每个分支上的名词后面的两组数字分别表示比对到该分支上的序列数与驻留在该节点上的序列数; 图13中为某单一水平物种分布情况,并非是序列分布。 多样品图图14:比对多个样品在不同分类学分支上序列数量差异。 图例解读: 比对不同样品在某分支上的序列数量差异,通过带颜色的饼状图呈现,饼状图的面积越大,说明在分支处的序列数量越多,不同的颜色代表不同
18、的样品。 某颜色的扇形面积越大,说明在该分支上,其对应样品的序列数比其他样品多。 多样品在做该分析时,建议样品数量控制在10个以内,或者将重复样本数据合并成一个样本后,总样品数在10个以内。 10.系统发生进化树 在分子进化研究中,基于系统发生的推断来提醒某一分类水平上序列间碱基的差异,进而构建进化树。 图例解读: 图15中表达的是序列进化差异情况,处在同一分支上的物种说明进化关系较近。 图15左下角的图例为距离标尺,分支距离越长,进化关系越远。 11.(un)WeightedUniFracPCoA/Tree分析 利用各样品序列间的进化信息来计算样品间距离,反映环境样品在进化树中是否有显著的微
19、生物群落差异。 PCoAprincipalco-ordinatesanalysis是一种研究数据相似性或差异性的可视化方法,通过一系列的特征值与特征向量进展排序后,选择主要排在前几位的特征值,PCoA可以找到距离矩阵中最主要的坐标,结果是数据矩阵的一个旋转,它没有改变样品点之间的相互位置关系,只是改变了坐标系统。通过PCoA可以观察个体或群体间的差异。 图例解读: 图16与图17中不同颜色代表不同分组; PCoA分析建议不分组时,样本数量不少于10个;多组样本时,每组样本数量不少于5个; 对于某一功能基因,进展进化树分析时,建议采用OTU数目控制在10,000以内,或者由客户指定分析优势OTU
20、个数。 12.NMDS分析 NMDSNonmetricMultidimensionalScaling常用于比对样本组之间的差异,可以基于进化关系或数量距离矩阵。 横轴与纵轴:表示基于进化或者数量距离矩阵的数值在二维表中成图。 图例解读: 图18中不同的颜色代表不同的分组; 建议不分组时,样本数量不少于10个;多组样本时,每组样本数量不少于5个; 图18中的点代表样本,点及点之间的距离表示差异程度。 13.含相似性树柱状图 根据样品中相似程度进展排布,并绘制对应样本树状图反映样本中群落构造。 图例解读: 图19中左侧是相似度树状图,样本之间的差异越小,样本便会处在相近的同一分支上; 右侧柱状图,
21、展示样本中微生物的群落构造。不同颜色代表不同物种。 14.Unifrac显著性差异分析 比拟样品间进化差异的显著性分析。 图例解读: 图20横坐标为两组样品; 纵轴坐标为unifrac进化距离序列差异。 15.单因素unweightedunifracPCA分析 在某单一因素上,进展unweightedunifracPCA分析。 图例解读: 图21横轴为不同变量本例为不同年龄阶段下的样本; 纵坐标为主成分,图21中显示同一年龄阶段内与不同年龄阶段间的由主成分导致的差异情况。 16. 个性化分析案例展示 案例描述一 Community composition of root-associated
22、fungi in a Quercus-dominated temperate forest:codominance of mycorrhizal and root-endophytic fungi. Ecol Evol. 2021 May; 3(5):1281-93. 样本来源 样本数量 高变区 测序平台 测序量reads 植物根系/泥土 159/38 ITS RocheGS 70,495 样本来自于以日本橡木为主的温带森林,使用454测序分析多个生态系统的真菌多样性。样品来源为植物根系样本,从12株植物中提取了159份根系。分析结果说明,外生真菌群落与根系内生真菌群落互相作用,维持了这种以日
23、本橡木为主的温带森林的生态环境。 图例解读: 图22构建了地下植物及真菌的相互关系网络; 图22中灰色圆点表示及植物共生的微生物物种,菱形代表真菌OTU,它们之间的关系用灰色连接线表示; 它们之间的密集程度越高表示它们之间相互作用被观察的次数越多; 互利共生真菌OTU用粉色菱形表示;寄生微生物OTU用橘色菱形表示;未知功能OTU用蓝色菱形表示; A:子囊真菌门;B:担子菌门;G:球囊菌门;U:表示门水平未知真菌。 案例描述二 Theignoreddiversity:complexbacterialcommunitiesinintensivecareunitsrevealedby16Spyros
24、equencing.SciRep.2021;3:1413. 室内微生物群落对日常生活中人类安康起着重要作用,尤其是医院的重症监护室。采用扩增焦磷酸测序研究ICU中微生物群落可以检测多种微生物序列,及现有的传统标准培养技术相比,有极大的优越性。 传统培养方式只能检测总细菌多样性的2.5。结合外部环境及物种系统发育谱分析发现,许多微生物及潜在的人类病原菌相关,当然也包含有益菌,一共7个门76个属。此外,丙酸杆菌属,假单孢菌属与伯克霍尔德氏菌被确定为感染的重要来源。在地板,医疗器械与工作间微生物组成有显著差异,但网络分析与一致性分子指纹印记分析发现该三个地点微生物组成也有一定的相似性。这些信息将帮助
25、加护病房进展新的公共卫生风险评估,帮助建立新的卫生协议,帮助深入了解医院获得性感染的情况。 样本来源 样本数量 高变区 测序平台 测序量reads 重症监护室ICU中地板,医疗设备及工作间外表样品 34 16SrRNA Roche454FLX+ 5000/样本 图23为丰度最高的40种菌(OTU)在不同采样区域的分布情况。不同颜色代表不同样本来源,绿色为ICU地板,红色为医疗设备,蓝色为工作间。结点面积表示该OTU在三个采样区域之间的相对丰度。 案例描述三 CMC与HIES均属于罕见的原发性免疫缺陷病。本文研究了CMC与HIES病人皮肤及口腔微生物及正常人群的微生物菌群差异。结果说明,相比正常
26、人体的微生物组成主要包括棒杆菌属与放线菌科等,病人的所携带的微生物中革兰氏阴性菌比例增加。而革兰氏阴性菌如不动杆菌属可以抑制机体对假丝酵母与金色葡萄球菌的免疫反响,结果很可能导致患者对这些菌感染的敏感性增加。本文的研究说明,此类免疫缺陷病人所携带的微生物菌群可以影响宿主的免疫防御系统,进而有可能通过基于微生物的辅助疗法来治疗免疫缺陷的患者。 样本来源 样本数量 高变区 测序平台 测序量reads 慢性皮肤黏膜念珠菌病CMC、高IgE综合症HIES患者皮肤及口腔微生物 60 16SrRNA V4区 MiSeq 6389/样本 图24为免疫缺陷疾病各亚型患者皮肤与口腔微生物群落构造的协方差。a图表
27、示35个皮肤样本的微生物多样性的差异;b图表示21个口腔样本微生物多样性差异。其中,绿色方块表示安康对照组,红色圆圈表示CMC,蓝色三角形表示HIES。安康人体皮肤外表的微生物主要包括葡萄球菌与棒状杆菌,而疾病组多为莫拉菌科。 培训文档境微生物多样性分析-根底 境微生物多样性分析-高级 常见问题Q1.能对哪些环境下的样品进展分析?A 目前已经成功运用Roche454技术发表的文章涵盖农业、土壤、林业、海洋、矿井石油等、人体医学等诸多领域,共计近1400篇,Paper目录可向本公司索取。其中有大量文献发表于国际顶级杂志上,包括Science、Nature、ISME、PNAS、AEM等。 Q2.D
28、GGE技术及Roche454高通量测序技术在环境微生物群落多样性研究中有什么区别? A DGGE等分子指纹图谱技术,在其实验结果中往往只含有数十条条带,只能反映出样品中的优势种群,也无法得到细菌种类及其绝对含量。而Roche454高通量测序技术能同时对样品中的优势种群及微量菌进展检测,获得样品中的微生物群落组成,并将其含量进展数字化。Q3.在医学领域有哪些应用? A 在人类的皮肤、口腔、呼吸道、胃肠道与尿道等处存在着大量及人体安康密切相关的正常菌群。它们能够合成并辅助机体吸收一些必需氨基酸与维生素;加工诸如植物多糖等人类饮食中一些难以消化的组分;占据人体的不同粘膜外表并产生天然的抗生素,抑制有
29、害菌的着落及生长。目前,人们已经采用传统方法对人人体微生物作了许多研究,但是这种方法需要对微生物进展纯培养,从而大大限制了我们对机体正常菌群,尤其是许多不可培养微生物的认识。第二代高通量测序技术能对人体菌群进展深度测序,从而准确检测到千分之一机体微生物群落的数量与组成构造的变化,突破了传统方法基于纯培养的限制,能使我们从整体上认识微生物群落及宿主之间的相互关系及其对人体安康的影响。例如,利用第二代高通量测序技术可以对处于疾病状态下的人体微生境进展研究,比拟分析正常与疾病状态下或疾病不同进程人体微生物群落的构造与功能变化;可以对正常人群及某些疾病患者体内的微生物群体多样性进展比拟分析,研究微生物
30、群落变化同疾病之间的关系;通过深度测序还可以快速地发现与检测常见病原及新发传染病病原微生物。 Q4.高通量测序是否需要做平行样? A 随着高通量测序的开展,提出研究过程中设置重复样的要求,也日益被大家所承受。在高通量研究中,设置重复样不仅仅是表达了科研的一种严谨态度,同时也表达了结果的真实性,防止了个体差异造成的影响。 目前的开展趋势:普通样品水体,土壤,处理样本内部设置重复等一般会要求设置重复样;对于大面积研究水体与土壤等样本时,可以采用多点采样混样研究。稀有样本稀有动物粪便,冰川极地冰样,过多采样会影响整体环境构成的样本等对于重复样设置的要求相比照拟宽松。 Q5.454一块PTP板最多可以
31、放多少样本?1个一样的样本分别研究细菌,真菌,古菌是否可以放在同一块PTP板上测序?A理论上,知道每个样本的测序量与一个454PTP板能够容纳的测序量,便能得知可以测序的样本量。实际上,考虑到测序质量的因素,不可如此计算上样量。一块PTP能够容纳多少样本往往应该考虑以下两个因素: 1. barcode数量:barcode是用于区分上样样品的标签序列,标签的个数限制着一块PTP能够容纳多少样本; 2. PTP板分区情况:PTP板分区越多可以增加上样的样品个数,但分区越多,会影响测序总量。 一块PTP板,在不分区的情况下,一般可以产出80万条序列,综合考虑barcode数量,分区情况以及单样本的测
32、序量,便可计算一块PTP板容纳的样本数量。 对于同一个样本,分别研究细菌,真菌,古菌一般是可以放到同一个板里进展测序,但要考虑到不同类型PCR产物的长度对测序结果的影响,如果PCR产物长度差异较大,那么不可以放入同一块PTP板进展测序。 参考文献AleksandarD.Kostic,DirkGevers,etal.GenomicanalysisidentifiesassociationofFusobacteriumwith colorectalcarcinoma.GenomeResearch,2021,22:292298. ChiachiHwang,FangqiongLing,etal.Mic
33、robialcommunitydynamicsofanurbandrinkingwater distributionsystemsubjectedtophasesofchloraminationandchlorinationtreatments.Microbiol, 2021,78(22):7856. CraigE.Nelson,StuartJ.Goldberg,etal.Coralandmacroalgalexudatesvaryinneutralsugar compositionanddifferentiallyenrichreefbacterioplanktonlineages.TheI
34、SMEJournal,2021,7:962979. ElieJami,AdiIsrael,etal.Exploringthebovinerumenbacterialcommunityfrombirthtoadulthood. TheISMEJournal,2021.2:111. ErikCvonRosenvinge,YangSong,etal.Immunestatus,antibioticmedicationandpHareassociated withchangesinthestomachfluidmicrobiota.TheISMEJournal,2021.33:113. KrissiM.
35、Hewitt,CharlesP.Gerba,etal.Officespacebacterialabundanceanddiversityinthreemetro politanAreas.PLoSONE,2021,7(5):e37849. SizhongYang,XiWen,etal.Pyrosequencinginvestigationintothebacterialcommunityinpermafrost soilsalongtheChina-Russiacrudeoilpipeline(CRCOP).PLoSONE,2021,7(12):e52730. VasviChaudhry,At
36、eequrRehman,etal.Changesinbacterialcommunitystructureofagriculturalland duetolong-termorganicandchemicalamendments.MicrobEcol,2021,64:450-460. WeiguangChen,FanlongLiu,etal.Humanintestinallumenandmucosa-associatedmicrobiotain patientswithcolorectalcancer.PLoSONE,2021,7(6):e39743. KatherineRAmato,Carl
37、JYeoman,etal.Habitatdegradationimpactsblackhowlermonkey(Alouatta pigra)gastrointestinalmicrobiomes.TheISMEJournal,2021,7:13441353. KabirGPeay,ChristopherBaraloto,etal.Strongcouplingofplantandfungalcommunitystructureacross westernAmazonianrainforests.TheISMEJournal,2021,7:18521861. DagneDuguma,PaulRu
38、gman-Jones,etal.Bacterialcommunitiesassociatedwithculexmosquitolarvaeand twoemergentaquaticplantsofbioremediationimportance.PLoSONE,2021,8(8):e72522. CoolonJD,JonesKL,etal.Long-termnitrogenamendmentaltersthediversityandassemblageofsoil bacterialcommunitiesintallgrassprairie.PLoSOne.2021,8(6):e67884.
39、 ZhaoJ,SchlossPD,etal.Decade-longbacterialcommunitydynamicsincysticfibrosisairways.ProcNatl AcadSciUSA,2021,109(15):5809-5814. BoLiu,LinaL.Faller,NielsKlitgord,etal.DeepSequencingoftheOralMicrobiomeRevealsSignaturesof PeriodontalDisease.PLoSONE,2021,7(6):e37919. BaldrianP,KolakM,StursovM,etal.Activeandtotalmicrobialcommunitiesinforestsoilarelargely differentandhighlystratifiedduringdecomposition.TheISMEJournal,2021,6(2):248-258. ScottTBates,JoseCClemente,GilbertoEFlores.Etal.Globalbiogeographyofhighlydiverseprotistan communitiesinsoil.TheISMEJournal,2021,7:652659. 第 22 页
限制150内