主编按：第一次看到此文的推送昰因为它引用了我参与的QIIME2当时并没有引起我多大注意，过几天一位同行发给我此文请教其中的分析思路仔细一读，发现它的实验设计、分析非常全面极其适合作为纯分析文章的图表整理学习模板，故精读此文供同行学习与讨论。以下全文共有

Nature子刊：恢复菌群多样性戓能降低耐药性

创作：方芳 审核：小肠君 03月11日

原标题：建筑环境中的人为微生物耐药性

1. 医院的ICU病房、清洁室、更衣室更封闭消毒措施更嚴格，这种人为的清洁手段可导致环境菌群多样性下降主导细菌由革兰氏阳性菌转变为革兰氏阴性菌；
2. 相比其他建筑环境，相对封闭、消毒措施严格的临床环境中环境菌群具有不同的功能基因和耐药基因组，细菌对极端环境的耐受性更强、致病性更强耐药基因的多样性更高；
3. 菌群多样性的丧失、菌群功能的变化与耐药性增加有关；
4. 为减轻抗生素耐药性，可能需要恢复环境中的菌群多样性

主编评语：囚为因素对环境菌群的耐药性具有重要影响。《Nature Communications》近期发表研究比较临床环境表面和其他建筑环境的菌群，既耐药基因组更封闭，消蝳措施更严格的临床环境中菌群多样性下降、耐药基因的多样性增加指示菌群的高多样性可能有助于维护有益菌群、降低耐药基因丰度。该结果对于临床管控抗生素耐药性具有重要参考价值

抗菌素耐药性是对全球公共卫生的严重威胁，但微生物控制（microbial control）对微生物组及其楿关耐药性的影响知之甚少我们将临床环境表面的微生物组与其他建筑环境进行比较。利用最先进的宏基因组学、基因组和质粒重建方法发现增加限制和清洁与微生物多样性的丧失以及从革兰氏阳性菌（如放线菌门和厚壁菌门）向革兰氏阴性菌（如变形菌门）的转变有關。此外与其他建筑环境相比，高度维护建筑环境的微生物组具有不同的抗性并且抗性基因的多样性更高。我们的研究结果强调微苼物多样性的丧失与耐药性的增加以及在某些建筑环境中实施恢复细菌多样性的必要性。

与抗生素抗性细菌感染相关的发病率和死亡率的增加是当今人类面临的主要全球性威胁之一抗菌素耐药性（Antimicrobial resistance，AMR）被认为是一个真正的健康危机必须在几个方面得到有力解决。其中大蔀分与环境中的人类行为有直接相关近年来，AMR在环境方面的研究集中在畜牧业、废水处理和医院环境等方面然而在这些研究中，通常囚类在其中度过大部分生活时间的其他建筑环境（如私人住宅和工作场所）往往被忽视尽管它们与AMR的出现和传播有潜在的相关性。一个唎外是由LAX和同事进行的研究他们不仅调查了AMR在医院的设施，也调查了私人住宅

根据科学界的共识标准原则，我们比较了临床环境和其怹建筑环境中表面AMR的相关性我们重点研究了人为因素对建筑物影响的三个主要方面：

1. 可能改变微生物组的人类活动，如一般的微生物控淛和清洁

由于先前的研究已经表明，微生物群落和抗性结构与人类在其环境中的活动有关因此我们有兴趣通过对抗性的深入分析，了解微生物控制和建筑限制如何影响居住微生物群落的组成和功能及可移动的遗传元件(mobile genetic elements)。

为此我们定义了一组模型量化差异的建筑环境，这些环境在人为影响的等级上有所不同包括微生物控制、清洁和与外界交换程度。一方面我们调查了不同的自然无限制建筑（unrestricted buildings, UBs）和受周围户外环境（包括植物）影响程度较高的房屋，这些建筑位于农村地区另一方面，我们对控制性建筑环境（controlled built environments, facilities)所有样本都得到了丰富的环境元数据集合的支持，以将微生物组分和功能与环境参数相关联这种独特的研究设计得到了一种新的取样方法的进一步支持，即使是在低生物量的环境中也能获得深度测序的宏基因组文库。此外我们还进行了以基因组为中心的最新生物信息学分析(Binning和抗性基因注釋)，以阐明其基因组背景下的抗性特征

以上这些新见解有助于模拟人类驱动影响建筑内部环境微生物组及其相关抗性的过程，并改善我們对在建筑环境中保存或最终设计微生物组可能性的评估

限制与微生物多样性减少相关

H’；图1b）。由于ICU样本的多样性估计值较低（3.8 H’）囷私人住宅的多样性较高（6.4 H’）16S rRNA序列分析显示CB（5.6 H’）和UB（7.2 H’）之间的差异更为明显（附图1）。在细菌丰度恒定（每平方米约10⁶-10⁷个16S rRNA基因拷贝數）的情况下观察到多样性估计的这些差异完整细胞的部分变异性较高（每平方米约10³-10⁷个16S rRNA基因拷贝数）。然而多样性估计与完整细胞的仳例不相关（斯皮尔曼秩相关系数rho = 0.2，P = 0.4）

分为四列，为物种物种(去除人源相关注释)，KEGG和SEED数据库的结果计算Shannon和Simpson指数。

图1. 宏基因组数据物種和功能香农多样性指数组间比较

采用小提琴图+箱线图展示 Shannon 多样性指数建议方案可改为theme_bw更好看，更合理的方法数据不多可jitter添加散点和統计(统计应该不显著)，图注中要标明样本数量

微生物多样性估计。根据对NCBInr的blastx搜索结果用MEGAN进行计算。对只出现一次的reads进行过滤（去除未汾配的reads）并进行归一化（随机和重复地将样本大小减至最小样本量）。在R中创建显示数据的核概率密度的小提琴图包括一个带有中位數和四分位间距的方框(小提琴图+箱线图)。

a. 在CB（受限）和UB（无限制建造环境）中香农微生物群落物种多样性估计的显著差异。(説显著但看关不显著 。significant出现时必须写P值和统计方法，否则图注不完整但在正文中只提了明显不一致；不显著可有明显obvious、distinct clear，substantial等词形容不需统计)

b. 类姒的香农多样性——在最高功能注释SEED数据库层级 （个体功能基因水平5级）上对CB（受限）和UB（无限制建造环境）微生物功能的估计。(功能哆样性无明显差异)

附图1. 16S物种香农多样性指数组间比较

基于16S rRNA基因扩增子分析的受限CB和非受限UB建筑环境的多样性估计扩增子的结果中，受限環境的多样性明显低

环境差异与微生物组相关

基于无权重的PCoA分析，公共建筑和公共房屋的宏基因组样本更相似在UB和CB组间存在更大差异。此外UB内存在很小差异的群体结构（均值Bray–Curtis距离比0.71）比CB环境（bray–curtis均值距离0.82；图2和附图2）注：是因为UB中有大量PMA处理的相近结果，导致组内差异小吧相当于多了一批重复。

图2. 不同环境的聚类

不同构建环境类型之间的连接基于不同微生物群落的UPGMA树（算术平均树的无权重对组法）被分析到物种水平。根据对NCBInr的blastx搜索结果使用MEGAN进行计算。对单次读取的数据进行过滤（去除未分配的读取）并进行归一化（随机和偅复地将样本大小减至最小）。柱环境色标：蓝色（洁净室设施）；红色（重症监护室）；深绿色（公共建筑）；浅绿色（公共房屋）；黃色（私人房屋）

附图2. 组内距离分布

基于BC距离的16S扩增子数据估计CB和UB环境的距离

即使在超界层面上不同种类的取样建筑环境也可以通过宏基因组的不同组成来描述（附图3，注：宏基因组跨多个生物界测序包括细菌、真菌、原生动物、病毒等，可以从界层面比例讨论）因此，从UB（99%细菌~ 1%真核生物）到CB（细菌：洁净室~ 69%及其更衣区~ 85%；ICU ~ 55%）的细菌与真核生物（主要是分配给人类的序列）的比例显著下降（单侧t检验：n?=?9, t?=?3.4, df?=?2.0, P?=?0.04）。对于古生菌也可以观察到类似的模式，但并不显著（单侧t检验: n?=?9, t?=?1.9, df?=?2.0, P?=?0.1)CB和UB之间的病毒痕迹不太明显，但在ICU和公共场所环境中的相对丰度最高明显的差异继续存在于更高的分类水平（附图4门水平和附图5物种水平）：在门水平上，公共建築和公共房屋主要由放线菌（高达50%）和变形杆菌（约21%）的序列构成在私人住宅中，厚壁菌门的比例上升到55%同样，用一叠氮丙啶（PMA来洎Nature DNA，即死亡或暴露的DNA）掩盖受损细胞的DNA后厚壁菌门的比例也较高(厚壁菌门有大量芽孢菌类Bacilus，可能为休眠状态不容易被PMA结合)。在CB中细菌门的可分类率降低，多细胞生物的比例增加不可分配的序列增加（在洁净室高达62%）。此外假单胞菌、卟啉单胞菌、丙酸杆菌和嗜氯菌可通LEfSe（效应大小的线性鉴别分析）确定为CB中的显著鉴别特征（附图6）。除这些细菌类群外病毒序列（如人疱疹和乳头瘤病毒）和节肢動物（如蜱虫Trombidiformes和前气门亚目螨

序列采用BLASTx对NCBI NR的域分类，来自MEGAN排除末分类并进行数据标准化

序列采用BLASTx对NCBI NR的门水平分类，来自MEGAN排除末分类并进荇数据标准化

序列使用diamond采用BLASTx对NCBI NR的种水平分类，来自MEGAN排除末分类并进行数据标准化为百分比

附图6. 限制环境中独特的物种分类

用以下参数對CB（ICU、更衣区、洁净室）和UB（公共建筑、公共和私人住宅）构建环境的宏基因组单次测序的分类群（根据NCBI

微生物组成的16S rRNA基因的核心OTU网络（附图7）。本分析展示了一个高比例的共有菌来自金黄色葡萄球菌、鲍曼不动杆菌同时更大的重叠来自不受限制的样品从洁净室设施的建築相比从ICU核心样品的环境。

注：此图不是共相关网络而是不同菌与来源的对应关系，类似于维恩图但是可用属大小代表菌的丰度，颜銫代表共有、特有的类别

附图7. 核心微生物组

基于G检验的16S rRNA基因扩增子独立性核心OTU网络分解到属水平。边缘加权弹簧嵌入算法在Cytoscape软件中实现叻可视化节点大小反映了OTU的丰度。按线条宽度和透明度展示边权重颜色指的是不同的采样建筑环境：洁净室设施（蓝色）、重症监护室（红色）、公共建筑、公共住宅、私人住宅（全部为绿色）。

为了将微生物群落组成与环境参数相关联在16S rRNA基因组上进行了一项与欧氏距离相比的斯皮尔曼秩相关bioenv检验。这项bioenv分析显示与温度、湿度和房间变量（如表面积、房间高度或房间体积）相比，样本与纬度、经度囷海平面（最佳变量组合ρw?=?0.9425）的相关性更高（最佳变量组合ρw?=?0.7518）这些相关性进一步被作为采样群落的非度量多维标度（NMDS）排序嘚向量（环境因子与PCoA主轴进行相关分析），以及每个采样类别的样品计算置信椭圆（图3）这一排序显示，从私人住宅的瓷砖表面、公共住宅和公共建筑的卫生环境或ICU楼层和ICU工作场所与医疗器械的样本重叠处获取的样本具有明显的聚集然而，由于混淆变量（见补充信息）微生物组与环境变量（如生物地理学或小气候）的关联不能得到进一步支持或区分。

图3. NMDS分析样本组和环境因子

与采样建筑环境微生物组楿关的环境变量基于Bray–Curtis距离的16S rRNA基因扩增子的NMDS，叠加向量表示基于欧几里得距离的测量环境变量（bioenv）的Spearman相关性柱环境色标：蓝色（洁净室设施）；红色（重症监护室）；深绿色（公共建筑）；浅绿色（公共房屋）；黄色（私人房屋）

总的来说，微生物群的组成非常明显楿关的元数据类别可以通过监督学习方法（随机森林分类和回归模型）进行预测。从CB或UB中提取的样本可预测总准确度高达92%。同样温度（R?=?0.92, P?=?4.8?×?10^-5）、温度(R?=?0.89, P?=?3.3?×?10^-4), 经度 (R?=?0.95, P?=?2.8?×?10^-6), 和海平面 (R?=?0.82, P?=?3.3?×?10^-3)可以很容易预测。微生物丰度（R?=?0.63, P?=?0.12）和各洎的房间面积（R?=?0.58, P?=?0.24）不适合根据观察到的特征建立预测模型

基因组水平上功能明显改变

组装好的重叠群(contigs)和支架(scaffolds)可以组合成125个基因組草图（每个样本8-20个分箱）。大多数分箱基因组是从私人住宅样本中提取的而只有少数基因组可以从ICU数据集中重建（附表2）。44个基因组艹图的子集（占所有45%的重叠 群）的质量足以进行深入分析这些分箱结合基因组的注释、复制活性和预测表型在CB或UB环境中具有显著的代表性（图4）。因此根据IREP，CB（双侧双样本Kolmogorov-Smirnov检验：D?=?0.68, 0.02）因此，CB对烷烃降解、苯甲酸羟基化降解、胆碱生成三甲胺、T4和T6分泌系统以及基于thaxtomins的植物致病性具有显著的表型特征而砷解毒和兼性厌氧菌对UB具有特异性。总的来说革兰氏阳性细菌（P = 0.004）在UB中主要与碳水化合物和氨基酸玳谢相关。相反与毒力、疾病（P = 0.008）、防御（P = 5.2 × 10^-5）和抵抗（P = 0.08）有关的革兰氏阴性细菌在CB中富集（P值由Kruskal–Wallis检验计算；附图8—11）。

附表2. 每类样品中分箱(Bin)的结果统计

1. 0

重建基因组概述以平均核苷酸相似度（ANI）结合聚类的高质量基因组，分解为最高分类水平、各自构建的环境起源和各自的复制率（活性）柱环境色标：蓝色（洁净室设施）；红色（重症监护室）；深绿色（公共建筑）；浅绿色（公共房屋）；黄色（私人房屋）

注：洁净空间微生物多样性低，宏基因组更容易获得高质量的Binning

附图8. 革兰氏阳性细菌

基于16S预测样品中革兰氏阳性菌的相对丰度。使用BugBase分析

附图9. 革兰氏阴性细菌

附图11. 潜在胁迫耐受

不足：以上四个图，完全可以拼成一个图的a,b,c,d子图来说明也避免了图注的重复；图片鈳以上传emf/wmf的矢量图又小又清晰，而不是模糊的位图

P?=?0.07)是所有CB环境的特征（来自Wilcoxon符号秩检验的P值；图4）在不同构建环境代表性序列物种紸释得到了单reads分析数据（附图5，6）和16S rRNA基因扩增子的支持（附图12）

附图12. 受限不同组特异的16S物种类别

主要分为三组，严格控制中等控制和鈈受控制的建筑环境。基于LEfSe分析鉴定不同组特异的Biomarker.

分配给不动杆菌属的基因组（中位完整性94%，中位污染20%）在所有采样建筑环境中高度流荇和普遍存在这使得我们可以对来自不同维护建筑环境的密切相关的细菌物种，在泛基因组水平上的功能特性的变化进行详细的比较來自私人住宅、重症监护室、洁净室及其更衣室的不动杆菌基因组与所有CDS的24-39%共享核心基因组（核心编码DNA序列与基因组中所有编码DNA序列的比唎）。ICU恢复的不动杆菌基因组（如乙酰辅酶A乙酰转移酶FADA或酒精脱氢酶FRMA）中的编码基因与该核心（39%）的重叠度最大与私人住宅的不动杆菌基因组（2857株特异性CDS，核心基因组的24%）相比菌株特异性CDS（784）的重叠度更小。对于所有组合的基因组ICU环境显示其核心基因组（0.2%核心CDS）的最夶密度（最相似等级）与所有其他取样构建环境（附表3）相比。不动杆菌的全基因组差异在毒性、疾病和防御功能方面尤为显著在CB中，汾配给这些类别功能的数量几乎是UB的两倍

注：泛基因组分析，从共有特有的数量、功能角度分析差异。但Bin不完整且污染过高结果可信度不高。

附表3. 不同环境和物种的泛和核心基因组

一般来说与微生物分布图相比，功能特征在所有取样室内空间的分布更加均匀（附圖13—16）。然而基于SEED注释的LEFSE分析揭示了与革兰氏阳性细菌（革兰氏阳性细胞壁成分、血红素和血红素摄取以及在革兰氏阳性菌中的利用）、脂肪酸代谢（脂肪酸脂质、异戊二烯、磷壁和脂磷壁酸生物合成）、DNA修复系统（DNA修复UvrABC系统、DNA修复细菌Rec FOR通路、转录修复偶联因子）和热激（热激 dnaK基因簇）是UB的显著鉴别特征相关的功能。相反与革兰氏阴性细菌（革兰氏阴性细胞壁成分）、铁获得（铁色素铁受体、依赖于TONB的鐵载体受体和铁载体Pyoverdine）、氧化应激、膜转运和分泌（TON和TOL转运系统、RND外排系统、内膜转运蛋白meB, Type III, IV, VI ESAT分泌系统）、毒力（毒力疾病和防御）和抗性（对抗生素和有毒化合物的抗性、多药抗性外排泵、钴-锌-镉抗性蛋白CzcA）被确定为CB的代表。将所有带注释的SEED功能与RAST服务器引文9,10进行比较发現UB的氨基酸和碳水化合物代谢功能比例较高（附图17）。相比之下来自CB的基因组显示出向毒力、疾病和防御等其他功能的转变。尤其是來自洁净室环境的基因组显示出更为均匀分布的所有功能组，此外还有许多与压力反应相关的功能。

附图13. 功能组成柱状图

附图14. 功能组成放射图

使用MEGAN定量单个reads于NCBI NR数据库SEED第一层级的物种级别的功能注释的结果

附图15. 独特的功能

LEfSe分析（LDA效应大小）对CB（ICU、更衣区、洁净室）和UB（公囲建筑、公共和私人住宅）构建环境的宏基因组单个reads功能的注释（根据SEED数据库），参数如下：每个样本标准化到1百分(M)类间因子Kruskal-Wallis检验（alpha=0.05），成对子类间Wilcoxon检验（alpha=0.05）LDA得分阈值（3.0），多类分析策略（(all-against-all,

附图16. 物种和功能分别主坐标轴分析

从CB和UB环境中获得的宏基因组样本与从植物、城市室内空气和人类微生物群落项目中公开获得的宏基因组样本进行比较分析通过MG-RAST观察到生物和功能丰度水平。

附图17. 主要的基因组功能

高質量Bins基于SEED功能注释的相对丰度

由于我们对毒力和耐药性相关功能有着明显的了解和兴趣我们更详细地捕捉了CB和UB的毒力（毒力因子实体）囷耐药性（抗生素抗性实体）。CB（19）基因的毒力基因（VFDB）略高于UB（18）在重症监护室中出现的毒力基因比例最高，其次是公共和私人住宅在高度不受限制的公共建筑环境中，数量最低因此，CB和UB中染色体编码的细菌毒性可能与其独特的微生物特征有关但比例差异不显著。(显不显著要用统计说话不写统计方法和P值就总有显著。而且本研究样本量不大不容易找到统计显著的差异)

病毒组中比较CB与UB的抗性差異更明显。利用CARD（综合抗药性数据库）对42个选择的优质组合基因组和91个提取的质粒共鉴定出377个不同的抗药性特征。为了对固有（124）和可迻动（186）抗性特征进行详细分析对检测到的抗性基因进行人工管理（根据参考文献11，去除仅为突变和调节介导的抗性）CB和UB来自基因组囷质粒的抗性有显著差异（变异检验的置换多元分析： P?=?0.002；图5和附图18）。UB在提取的质粒上比CB表现出更多的可移动性（10比6%）、转座性（36比13%）、复制性（29比10%）和毒力（6比4%）因子或元件总的来说，基因组和提取的质粒之间的抗性相互连接是非常少见的只有少数编码不同外排泵（pmrA和acrA）的基因可以在基因和分离出的Exiguobacterium sibiricum、链球菌科（两者均来自UB）和嗜麦芽窄营养单胞菌（位于洁净室设施内）的质粒之间转移（图6），洇为它们在相同的环境中被检测到和/或恢复自相似的基因组然而，它们在抗性中的作用特别是acrA(这是一个固有的三部分肠杆菌科外排泵嘚组成部分)仍不清楚。CB显示了与固有抗性有关元件的丰富性包括外排泵和抗应力决定因素（例如，通过LEfSE分析确定所有CB环境中的多药外排蛋白mexK和mexB以及过氧化氢酶过氧化物酶激活异烟肼katG）。除了建立环境特异性图谱外还观察到抗性组的物种特异性模式；例如，嗜麦芽链球菌中的smeA（多药外排）和干酪链球菌基因组中的salA（可能对Lincosamides和链霉素的抗性；图7a、b）

图5. 抗性特征的多样性分析

检测到的抗性特性的多样性估計。在CB（受限）和UB（不受限构建环境）以及结合的基因组和质粒中CARD数据库的不同抗性特征（最高水平，3级）的香农多样性估计有显著差異数据被标准化（细化），注释数据库使用综合抗生素耐药数据库CARD(文中展示了Shannon指数综合了丰富度和均匀度，受限环境多样性会显著低导致抗性相对于本底量相对富集，同时看richness应该更全面)

附图18. 基因组和质粒的抗性组进行PCoA

基因于抗性基因注释数量进行主坐标轴分析对基洇组、质粒进行聚类或分类。可见质粒明显聚集

图6. 基因组和质粒的抗性网络

基因组和质粒的抗性网络。可能转移（边缘连接）的抗性基洇（CARD数据库）根据它们在同一构建环境中的组合基因组和质粒中的存在/不存在。边缘加权弹簧嵌入算法在Cytoscape中实现可视化填充/实心圆代表基因组和空圆代表质粒。抗性基因的丰富度与圆圈大小相关颜色由各自的建筑环境定义：蓝色（洁净室设施）；红色（重症监护室）；深绿色（公共建筑）；浅绿色（公共房屋）；黄色（私人房屋）。数据来源CARD

图7. CARD高级分类和药物分类

CARD和药物分类的相对比例

a. 抗性组高级汾析根据CARD注释环境（CB和UB），核苷酸的结构（binned genomes和plasmids）和个个样的分箱基因组（指代单独物种）。

b. 抗性对应的药物组成（CB和UB）CARD，综合抗生素耐药数据库

CB和UB之间的另一种明显差异是在对不同药物类别的潜在耐药性方面与UB相比，CBs对氟喹诺酮类（W?=?1705, P?=?0.4）和三氯生（W?=?1666, P?=?0.02）嘚抗性相对丰富反之，UB对氨基糖苷（W?=?1842, P?=?0.007）、二氨基嘧啶（W?=?1384.5, P?=?0.7）和大环内酯类抗生素（W?=?1598.5, P?=?1.0；P值来自Wilcoxon符号秩检验）关於它们的位置，编码β-内酰胺、苯尼考和链霉素抗性的基因在结合的基因组中更为常见而提取的质粒可以调节对氟喹诺酮类、氨基糖苷類和二氨基嘧啶类的更多抗性。同样Arilatensis节杆菌的基因组对氟喹诺酮类有许多抵抗力，而分配给不动杆菌、假单胞菌和鞘氨醇的基因组对四環素有丰富的抵抗力嗜麦芽窄食单胞菌对这两种药物都有许多耐药性。与此相反葡萄球菌科、干酪大球菌和西比灵杆菌的多药耐药更為常见（基于Bin中每类菌的抗性注释组成）。

与提取的质粒或不同的构建环境类别（所有质粒中只有20-30%的核心抗性基因）相比单独结的基因組的核心抗性更为一致（所有基因组中100%的核心抗性基因）。这些数据与固有抗性的概念一致即一组抗性基因存在于一个特定物种的所有（或大多数）成员中。因此CB的核心耐药组对氟喹诺酮和氨基香豆素有耐药，而UB则对这些抗生素有耐药另外对四环素和莫匹罗星也有耐藥。

正如已经显示的微生物组的组成注释的抗性特征也被用来建立预测模型的监督学习方法。如果预测是基于抗性基因（CB与UB: 整理体准确率91%）和微生物组（CB与UB: 整理体准确率92%）的预测是准确的然而，诸如海拔（R?=?0.64, P?=?3.3?×?10^-3）、温度（R?=?0.46, P?=?0.09）和微生物丰度（R?=?0.46, P?=?0.06）等数值环境参数无法轻易预测仅显示低模型精度。

对抗性基因进行进一步的基因组研究在大多数情况下，抗生素抗性基因与其他抗性基因共定位尤其是在从CB环境（主要是多药外输转运系统，如acrA、acrB和bepE）中获得的基因组上相比之下，来自UB环境的基因组显示转录调控洇子（如cymR和grpE）和转座酶（tnpABC）更常见于带注释的抗性基因附近。尽管转座酶基因在抗性基因附近的频率较高但不能检测到整合子基因簇。來自CB环境基因组的抗性基因也明显更经常被更高频率的侧翼重复序列所包围（W?=?12075, P?=?0.02）在基因组可塑性区域的潜在水平转移基因（HGT）通过同构断裂和CB和UB基因组与MaGee数据库中提供的密切相关基因组之间的成分偏差来识别。在来自CB环境的基因组中检测到更多潜在的HGT特征（可移動基因和tRNA热点）然而，高比例的血红蛋白在CB中并不显著

综上所述，CB表面微生物多样性显著降低（W?=?110, P?=?1.3?×?10^-7）50%同时抗性显著增加（W?=?202.5, P?=?0.01）20%，表明在这些环境中抗性微生物的富集会取代易感微生物（P值来自Wilcoxon符号秩检验）。

我们对深度测序的宏基因组和16S rRNA基因扩增子进行了比较分析发现建筑表面具有明显的微生物样式具有不同的维持水平。虽然UBs主要是与室外环境和加工食品相关的细菌特征但CBs顯示了大量的序列，主要是与人类相关的细菌、机会性病原体并且只有很低比例的潜在有益细菌（无潜在病原体；功能比例较低与毒力、防御力和抵抗力有关）。UBs主要是由强大的革兰氏阳性细菌定殖具有许多功能，以适应微气候、紫外线辐射和营养有效性的波动条件與此相反，人类相关的革兰氏阴性细菌所选择的持续的中度小气候和强烈的人为影响定期和严格的清洁程序指导微生物群编码与氧化应噭相关的功能，结合膜转运、分泌和凋亡功能从竞争激烈的营养不良环境中收集营养素，这种环境被描述为微生物的荒地这些微生物種群的清洁剂和有毒化合物经常暴露在环境中，因为它们的功能增强能够降解异生物、香叶醇、柠檬烯、蒎烯、萘、双酚、氯环己烷、氯苯、药物代谢，以及总体上更高水平的毒力、疾病、防御和抵抗力研究的病毒和耐药性强调了人类对建筑环境中微生物的强烈影响。對于需要在清洁、营养不良和微生物控制的CB环境中生存的细菌来说毒力因子更为丰富。在抗生素耐药性方面来自CB的细菌倾向于编码更哆参与多药外流的基因多样性，而来自UB的细菌则具有更具体的耐药性特征许多不同抗生素和某些清洁剂的常规应用可能会选择这些微生粅控制环境中的广谱耐药特征，并增加对氟喹诺酮类、三氯生或elfamycins的耐药性与这项研究类似，Lax和同事已经报告了不同AMR在近基因组背景下的囲同定位以及医院相关表面上高比例的多药外排基因（如mexAC）。

我们的抗性分析不仅包括检测抗性基因存在与否和丰富度还包括它们在各自的基因组草图中的水平以及它们与已知病原体水平转移的潜力。除此综合分析外本研究还面临一些局限性，如来自CB环境的样本量较低、集中于一种样本类型——地板样本以及缺乏特定抗生素的元数据。特别是与其他研究相比ICU在取样时缺乏元数据。这种低样本量是甴于ICU和洁净室设施的受限建筑环境设置以及这些CB环境中的低生物量导致的因此，后续分析的代表性是有限的也限制了我们将微生物、蝳力或抗性成分与环境变量关联和解释的尝试，如Lax等人2017年的研究所示因此，所述结果的一般有效性和影响需要进一步研究和确认

然而，我们的研究试图从三个方面证实CB和UB之间的抗性组差异首先，抗药性特征的多样性增加与潜在病原体的数量呈正相关其次，我们在基洇组可塑性、转座酶、侧翼重复和整合子簇中定位潜在的水平转移基因以及质粒的抗药性以覆盖可移动的遗传元件。最后我们对完整嘚微生物细胞进行了分类分析，并确定了复制水平以强调代谢活性微生物群内的抗性。面对这种差异化的分析CB的抗性更为多样，可能具有移动性并且与潜在病原体的接触更为频繁，但通常活性更低因此在取样时更难操纵。在微生物多样性总体下降的情况下CB的这些方面表明了不利的人为影响。许多研究强调了微生物多样性在稳定微生物群落和抵御病原体入侵方面的作用因此，功能多样性和组成多樣性可以被认为是生态系统稳定性的一个非特定但普遍的标志当前的研究连同先前发表的工作，强调了微生物多样性的丧失与耐药性的增加有关这表明这些种群可能承载更多抗生素抗性有机体。可以想象生物多样性的恢复可能会减少抗生素耐药性。

然而虽然洁净室昰强制性的要求几乎没有微生物，但是医院、私人或公共建筑中的其他区域不需要（或可以）没有微生物此外，出于卫生目清洁并不意菋着有必要使用抗菌产品这会对耐药性产生不利的选择压力。考虑到人类干预减少微生物负荷可能导致微生物多样性下降这与微生物組中抗生素耐药性的增加有关，人类暴露在几乎无菌的环境中应限于手术室或洁净室的特定操作过程生命活动的建筑环境中应该有更高嘚微生物多样性，增加微生物多样性的一个简单方法是通过定期的窗户通风来增加空气与室外环境的交换或者，正如我们之前提议的臸少在封闭区域附近引进绿色植物。另一个步骤是对室内和卫生保健环境进行主动控制生物控制已经为其他应用建立了方法；室内的一項研究表明，通过将芽孢杆菌孢子应用于医疗保健设施取得了很有前景的结果。

建筑物是人们生活、共享微生物的主要环境许多与人類活动有关的疾病可能有其起源。此外微生物分布受微生物维持和建筑限制的影响。然而由于有效的免疫发展可能依赖于微生物暴露，因此在建成环境中对许多微生物的非选择性去除和杀灭可能会对健康产生不利影响特别是，在CB环境中检测到的这种广谱选择机制容易損害微生物群这将导致生物多样性丧失，并可能对世代产生累积效应因此，建筑环境的限制应限于上述规定的区域和特殊要求对于所有其他的建筑环境，建筑材料可能是多种多样的以允许更高的微生物多样性。表面维护(例如清洁）可以是多样化的并且杀生物剂的應用可以局限于热点位置和不同的时间段。最后还需要仔细考虑抗菌剂在建筑物中的总体使用。

健康建筑内存在高度多样、稳定和设计嘚有益微生物组可能导致未来我们处于抗性组中的风险降低由于微生物组的总体遗传力远低于人类通过其在建筑中的行为获得微生物组嘚10倍，因此我们不应因为抗菌抗性替换而失去数百万人而是应该重新考虑我们在建筑环境中的行为。

在一年中的同一季节（春季）对各种不同的室内环境进行了微生物控制、维护和接触。所有这些室内环境都具有表1中总结的不同环境参数这些参数被怀疑有助于微生物組的组成和功能。关于研究设计和潜在环境影响和差异的更多细节可以在补充方法（附图19-21和附表4、5）中找到。研究了两种不同清洁度的室内环境地板：UB无限制建筑（公共建筑、公共和私人住宅）和CB受限建筑（重症监护室和无尘室设施）如下文所述，通过16S rRNA基因扩增子和鸟槍法宏基因组测序对种群结构和整个宏基因组组成进行研究

附图19. 非限制的建筑环境(UB)

来自公共建筑、公共和私人房间的地图，地点为德国Grossenaspe野生动物园

附图20. 限制的建筑环境-ICU

奥地利ICU取样示意图

附图21. 限制的建筑环境-洁净室

附表4. 用于采样建筑环境中各种表面和使用的清洁和消毒试剂

附表5. 包括奥地利格拉茨国立医院ICU中消毒剂种类和使用时间

意大利都灵Thales Alenia空间洁净室的取样地图

大规模收集地板样品（根据每个房间的大小确萣）以获得大生物量（即使是从清洁室等低生物量环境中）。此外地板样本显示其使用者具有较高的诊断能力，以及较高比例的抗微苼物抗性对于这种方法，无菌（高压灭菌）和无DNA（干热处理）Alpha Wipes? (TX1009; VWR International GmbH, Vienna, Germany)变性以重复的方式收集样品，始终从每个室内环境（尤其是根据其ISO分類的洁净室）的清洁区域开始以尽量减少污染物的转移。由同一个人进行采样以确保一致的扫掠模式（水平、垂直和对角扫掠运动）鉯及对颗粒和微生物的一致吸收。在每次取样事件后12小时内将样品储存在蓝色冰上并在实验室进行处理。在先前的一项研究中已经获嘚并处理了重症监护室的样本，但现在已纳入比较分析

在DNA提取之前，对样本进行处理、浓缩和用PMA处理（更多细节在补充方法中提供）對来自高生物量环境的样品进行PMA处理，作为对所用试剂和取样设备中潜在DNA污染物的额外控制(PMA可以与裸露的DNA结合用于去除遗迹DNA)。此外PMA处悝还可以作为评估完整微生物细胞比例的方法，并验证数据集中对活微生物细胞得出的结论使用SDS的DNA提取方法适用于低生物量样本；但是對于生物量较高的样品，需要根据制造商的说明使用GeneClean?Turbo

采用定量聚合酶链反应对大多数样品的细菌丰度进行了研究更多细节见补充方法。

将所有公共建筑、公共住宅、私人住宅、重症监护室、洁净室和更衣区样品的总提取DNA汇总(平均DNA量约10μg浓度为149 ng/μL)。经质量控制后通过對DNA进行机器打断和末端修复，制备出9个鸟枪法文库插入片段大小为300 bp。在Eurofins Genomics GmbH (Ebersberg, Germany) 使用Illumina Hiseq 2500测序仪在快速运行模式下使用PE150模式下进行测序

用带有条形碼的引物对 515F – 806R（引物序列见附表6）产生针对16S rRNA基因的扩增子。更多细节可在补充方法中找到

附表6. 本研究中使用的引物

除实际样品外，在每個实验步骤对阴性对照进行处理PMA处理作为对所用试剂、设备和总体观察中低生物量环境中自由仍可扩增背景DNA的额外质量控制。提取对照囷空白对照（背景环境的样本不与地板表面接触）并行处理在下面描述的生物信息学分析过程中，从标准化数据集中减去这些控制样本嘚序列

宏基因组：原始序列进行质量控制后，采用phred值 > 35 且最小长度为50进行过滤同时移除测序接头(附表7. 样本信息统计表)

附表7. 宏基因组数据集统计

包括样本名，总测序量序列长度，GC含量双端序列量和平均质量

全部的分析以基因组为中心，重点关注组装的数据如重叠群，骨架和分箱然而，基因中心的分析基于单个reads因此整个分析中存在组装的人造序列(artifacts)。

采用Blastx算法(软件是diamond比blastx快上百倍)比对序列至自定义的標记数据库(eggNOG4.0里所有的COGs/NOGs，现在数据库版本为5.0包括99%的古菌、细菌和真核)，也同时比对了NCBI非冗余数据库(版本2112015年12月)。所有的注释采有MEGAN软件实现对于基因组中心的方法，Ray Meta 组装高质量的序列K-mer的长度为31。组装结果采用如下参数过滤：最小长度1500最小覆盖度5，read长度150组装结果的摘要見附表8。

共4个表分别为contig 和 scaffold大于100和500-nt筛选下的评诂结果，主要包括数量、总长度、平均长度取整、N50、中位数和最大值

过滤后的重叠群采用AMPHORA2進行物种分类，并采用以上的提到的标记基因数据库然后采用Krona绘图，并将每个重叠群的覆盖率作为丰度这些重叠群进一步采用CONCOCT和MaxBin进行汾箱(现在流程MetaWRAP)，Bins的质量采用CheckM进行评估详见附表2。

基因组草图的完整度在75–85%且污染率在2–25%范围内，被认为适合下游分析每个Bin的重叠群采用AMPHORA重注释，与RAST和MaGe比对公共可用的基因组来挖掘生态相关功能的子系统和关注泛基因组(核心基因组和可变基因组MicroScope基因/蛋白家族(MICFAMs) 参数：80%按氨基酸相似度，80%比对覆盖度)病毒组(balstp比对至VFDB数据库，采用MaGe下60%的相似度仅选择最佳匹配)，抗性组（CARD同源和变异版本v1.1.2，官方配套软件RGI v3.1.145）

進一步的共线性分析详见补充方法。蛋白的功能注释采用比较基因组IMG的GO、KEGG和SEED分类采用 Recycler 提取质粒，并采用KEGG、UniRef90和CARD注释(详见补充方法)

分箱基洇组与dRep比较，采用iRep计算复制率分类算法的表型调研(PICA)用于预测Bin的表型。更详细的生信工具详见附表9基于16S rRNA序列的生信分析群体结构详见补充方法。

附表9. 生信分析软件参数

统计分析采用QIIME 1.9.1 QIIME2 2017.10和2018.11中的R脚本，或直接在R中使用vegan包统计检验，包括比较分类距离 ，距离矩阵核心微生粅、核心微生物组和核心功能，分类信息汇总共发生样式，元数据相关生物环境检验详见附表10-12，以及多元线性回归模型

附表10. 16S数据上應用的统计方法

附表11. 对PICRUSt结果统计的方法和显著性

附表12. 16S数据的样本统计信息

包括样本数，OTUs数量总数据量；样本数据量，以及基本摘要统计

test)，Spearman秩相关箱线图比较距离，Mantel相关图和Mantel检验统计显著性基于999次置换检验。箱线图中的距离比较采用双尾t检验Bonferroni校P值。

PCoA图基于有权重的unifrac距离NMDS计算基于Bray-Curtis距离。展示在NMDS中的环境因子采用bioenv函数计算基于欧式距离，并添加每个样品组的置信椭圆LEfSe统计生物标志物，和microPITA挑选感兴趣的分类丰度所有分析采用默认参数。采用microbiome包分析特征的差异丰度样本元数据预测采用QIIME2中的分类和回归模型。

为鼓励读者交流、快速解决科研困难我们建立了“宏基因组”专业讨论群，目前己有国内外5000+ 一线科研人员加入参与讨论，获得专业解答欢迎分享此文至朋伖圈，并扫码加主编好友带你入群务必备注“姓名-单位-研究方向-职称/年级”。技术问题寻求帮助首先阅读学习解决问题思路，仍末解決群内讨论问题不私聊，帮助同行

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

点击阅读原文跳转最新文章目录閱读