一、相关信息
原文题目：Lossofsoilmicrobialdiversityexacerbatesspreadofantibioticresistancerecruitingbeneficialmicrobes第一作者通讯作者：Qing-LinChenYong-GuanZhu期刊：SoilEcologyLetters
图1：文章信息二、文章摘要
生物多样性的丧失严重危害人类社会所依赖的生态系统过程。自然生态系统对外来物种入侵的抵抗能力不同，而这种抵抗能力可能取决于系统的多样性。然而，人们对微生物多样性为应对微生物入侵所提供的屏障知之甚少。耐药细菌的肆意横行是对21世纪公共卫生的严重威胁。因此，该研究探讨了减少土壤微生物多样性对抗生素耐药性传播的影响。
该研究使用稀释到灭绝，并结合高容量定量PCR的方法研究了微生物多样性和抗生素耐药性的入侵之间的关系。研究结果表明，微生物多样性与耐药基因丰度呈负相关。并且在考虑了孵育时间和微生物丰度等其他潜在驱动因素后，这种相关性得以维持。这说明高微生物多样性可以作为一个生物屏障，抵抗抗生素耐药性的传播。
研究意义：探讨了土壤微生物多样性在生态系统健康中的作用。
三、研究方法
1、实验设计
年8月于浙江省嘉兴市采集水稻农田表层土壤(0~20cm)。土壤过筛分至2mm，部分土壤经35kGy的γ射线辐照灭菌。将剩余的土壤作为接种剂进行梯度稀释接种实验。将0.5g灭菌土壤涂于含有色氨酸大豆肉汤(TSB)和马铃薯葡萄糖琼脂(PDA)的平板上，以检测辐照土壤的无菌性。6天后未观察到细菌或真菌的生长视为灭菌完全。
不同微生物多样性的土壤样品制造过程如下：将30g未灭菌土壤与50mL无菌脱盐水用紫外线灭菌混合器以最大速度搅拌5分钟配制成初始土壤悬浮液。将未灭菌土壤制成的土壤进行梯度稀释。选取4个浓度梯度，依次为未稀释（D0）、稀释倍（D1）、10-4（D2）、10-8（D3），分别取5ml悬液加入到含有30g灭菌土壤的烧瓶中。使未灭菌土壤接种量依次为10-1、10-2、10-3、10-4。阳性对照和阴性对照分别为非无菌土壤(S)和无菌土壤(SS)接种相同体积的无菌水。
每个处理组为4个重复。每个烧瓶加入1.2g猪粪作为耐药基因（ARGs）的来源。盖的无菌盖，在20°C、65%持水量（使用无菌脱盐水进行控制）条件下培养。每个处理组分别在第0、10、30、60和90天收集样品，用于分析细菌丰度、多样性，以及抗生素耐药性。
图2：实验设计2、分析内容
分别取0.5g土壤样品，提取DNA，并进行16SrDNA扩增子测序以分析细菌群落结构。采用高通量PCR的方法分析种ARGs、8种转座子、2种整合子的相对丰度。
四、图文结果
1、微生物群落的组成和多样性
从所有样本中共获得条高质量序列，每个样本的序列范围为~条。土壤微生物群落分布在17个门中。群落主要由放线菌门、厚壁菌门、变形菌门、拟杆菌门、普朗克菌门和绿菌门组成，这些是世界范围内最常见的土壤细菌。与放线菌门和拟杆菌门相比，厚壁菌门的相对丰度随时间显著下降。据报道，厚壁菌门丰度的下降预示着随着时间的推移，土壤环境中肥料衍生微生物群的减少。不同稀释处理组间优势菌差异不大。但经过大量稀释后，变形菌门显著富集(P0.01)，特别是在孵育第10天。
通过计算Good’scoverage发现测序深度在所有样本中具有高度可比性，从94.25%到98.51%，表明测序深度可以足够可靠地描述细菌微生物组。正如作者预估的那样，稀释处理中物种丰富度(观测到的OTUs)逐渐降低。表明稀释-灭绝方法可以有效地控制土壤等复杂自然生态系统的微生物多样性。物种丰富度也随着时间的推移而降低。
基于样本间加权UniFrac距离的主坐标分析(PCoA)进一步表明，稀释和培养时间显著改变了细菌组成的总体特征。沿着PC2轴，微生物群落通过初始稀释水平聚集在一起，解释了10.41%的总变化，而沿着PC1轴，微生物群落通过培养时间聚集在一起(解释了53.55%的总变化)。这一结果表明，时间对细菌组成的影响比稀释作用更大。
此外，作者还使用荧光量PCR检测了16SrRNA基因绝对丰度，以此检测了细菌丰度变化。结果表明，S和D0处理的细菌丰度在开始时(第0天)较高(因为原生菌群更丰富)，但随着时间的推移而下降。与S和D0相比，D1、D2、D3和SS(原生菌群多样性较低)的细菌丰度则更为稳定。
2、微生物抗药性组装和多样性
采用HT-qPCR检测了各样本ARGs的丰度和多样性。共检测到个基因，个ARGs，1类整合子整合酶基因和8个转座酶基因。每个样本中检测到的ARGs数目在47~之间，平均85个。这些ARGs可对大多数常见的人类使用的主要抗生素类产生耐药性，如氨基糖苷类、beta-内酰胺类、大环内酯林考胺-链霉素B(MLSB)、多药外排、磺胺类、四环素和万古霉素等。这些ARGs的丰度范围在5.53×~1.65×拷贝/g土壤之间。在开始时(第0天)，S组的ARGs丰度更高，这是可能与土著耐药微生物类群的高丰度有关。这一观察结果与细菌丰度一致。稀释处理对培养10天后ARGs丰度开始变化，随着稀释程度的增加，ARG丰度增加。
图3：OTUs及ARGs丰度3、微生物抗药性与多样性的关系
作者使用普通最小二乘(OLS)回归模型来确定物种多样性和ARGs丰度之间的关系。分析结果显示，物种丰富度(OTU丰富度)的减少与ARG丰度的增加呈线性相关(P=0.,R2=0.)。但相反地是，微生物Simpson多样性与ARGs的丰度呈高度负相关(P=0.,R2=0.)。相比之下，均匀度与ARG丰度无显著相关(P=0.,R2=0.)。
图4：OLS回归模型进一步分析发现，OTU丰富度与OTU多样性正相关(P0.,R2=0.)，而与OTU均匀度负相关(P=0.,R2=0.)。OLS回归模型显示，ARGs丰度随着时间的推移而下降(P0.,R2=0.)。这一结果表明，施肥后的ARG丰度并没有保持急剧上升的趋势。
图5：ARGs丰度与孵育时间的关系因此，作者根据已知的效应和关系生成了SEM模型，以探究当考虑到其他可能的驱动因素时，微生物多样性和耐药性之间的关系是否能够得到维持？SEM模型是一种可以分离多种影响途径的检验方法，对于探索生态系统中复杂的关系网络较为有用。作者将稀释倍数、孵育时间和MEGs（转座子和整合子）以及细菌丰度都纳入到该模型中。模拟结果显示此模型解释了ARGs模式中86%的方差。培养时间和稀释度可以直接影响ARGs的分布模式，也可以通过强烈影响细菌的多样性和丰度以及MGEs的丰度来间接影响ARGs。细菌多样性显著且直接影响ARGs丰度(r=-0.09,P0.05)，并通过强烈影响MGEs丰度间接影响ARGs(r=-0.07,P0.05)。然而，细菌丰度对ARG丰度并没有显著的直接影响(r=0.02,P0.05)。
图6：SEM模型