中国海洋大学学报自然科学版  2021, Vol. 56 Issue (6): 88-98  DOI: 10.16441/j.cnki.hdxb.20200135

引用本文  

汪瑶, 祁建华, 张婷, 等. 青岛多种天气下生物气溶胶中细菌群落特征[J]. 中国海洋大学学报(自然科学版), 2021, 56(6): 88-98.
WANG Yao, QI Jian-Hua, ZHANG Ting, et al. Bacterial Community Characteristics in Different Weather Conditions of Atmospheric Bioaerosols in Qingdao[J]. Periodical of Ocean University of China, 2021, 56(6): 88-98.

基金项目

国家自然科学基金项目(41775148);中央高校基本科研业务费专项(201762006)资助
Supported by the National Natural Science Foundation of China (41775148);the Fundamental Research Funds for the Central Universities (201762006)

通讯作者

张大海, E-mail:dahaizhang@ouc.edu.cn

作者简介

汪瑶(1994-),女,硕士生。E-mail:1280314744@qq.com

文章历史

收稿日期:2020-04-28
修订日期:2020-06-05
青岛多种天气下生物气溶胶中细菌群落特征
汪瑶1,2 , 祁建华1,2 , 张婷3 , 张大海3     
1. 中国海洋大学海洋环境与生态教育部重点实验室, 山东 青岛 266100;
2. 青岛海洋科学与技术试点国家实验 室海洋生态与环境科学功能实验室, 山东 青岛 266237;
3. 中国海洋大学化学化工学院, 山东 青岛 266100
摘要:为探究雾、霾、沙尘等多种天气下生物气溶胶中细菌群落结构特征,于2017年12月—2018年4月采集了青岛晴天和多种特殊天气下的大气生物气溶胶样品,运用16S rRNA高通量测序方法分析了样品中细菌群落结构。结果表明,在雾、霾、沙尘天均检测到对人体健康存在潜在危害的条件致病菌,其中变形菌门(Proteobacteria)均是优势细菌门,不受天气条件变化影响;霾天时其相对丰度最大。除此以外,沙尘和霾发生时,细菌群落结构显著变化,丰富度和多样性高于晴天,但随霾污染程度的加剧,群落丰富度和多样性下降;雾天多样性下降,丰富度变化不显著。在霾天,不同气团来源中细菌群落结构也存在显著差异。属水平上,霾天不动杆菌(Acinetobacter)和鞘氨醇单胞菌(Sphingomonas)占主导地位;雾天微球菌(Micrococcus)、副球菌(Paracoccus)和假单胞菌(Pseudomonas)占主导地位;沙尘天的优势菌为芽孢杆菌(Bacillus)。RDA分析显示O3是影响菌群的主要环境因子。
关键词生物气溶胶    细菌    16S rRNA    群落结构    健康风险    

大气中生物气溶胶是由含有微生物或生物大分子等生命活性物质的微粒组成,包括细菌、真菌、病毒、花粉、动植物碎裂分解体等[1]。细菌是生物气溶胶的重要组成部分,在近地表大气层中无处不在。它们可在大气中存活较长时间并且借助介质进行传播,一般浓度在104~108 cell/m3之间变化[2-5]。有研究表明,在PM2.5和PM10的可识别生物成分中,细菌的相对百分比占80%以上[6-7]。这些在空气中传播的细菌有些属于人类病原体,可引发过敏和哮喘反应,急性、慢性疾病以及动植物疾病等[8-9]。有些细菌还可以作为云凝结核,如假单胞菌(Pseudomonas)、根瘤菌(Rhizobium)和芽孢杆菌(Bacillus)[10]。细菌群落反映生物气溶胶中细菌的种类、数量和组成,群落结构在不同天气条件下的改变可以揭示雾、霾和沙尘事件对大气中生物活动环境的影响[6, 11]。同时,细菌群落的变化,特别是致病菌的增加还对人体健康构成威胁,一定程度会影响空气质量[12]。空气受到污染时,细菌群落结构随污染类型和污染水平而变化。与清洁大气相比,污染天大气中微生物浓度和致病菌都会增加[13],群落结构和优势菌也会发生变化。近年来,霾成为国内常见的污染天气,对北京一次霾事件的研究发现,重度霾天中存在1 000多种可识别微生物,致敏性和引起呼吸系统疾病的细菌丰度显著增加[6]。在空气质量为优及轻度污染时,济南大气中优势细菌属为蓝藻(Cyanobacteria);在中度污染时葡萄球菌(Staphylococcus)为优势菌;重度污染时异常球菌(Deinococcus)和芽孢杆菌(Bacillus)为优势细菌[16]。但也有研究发现,北京霾天大气中细菌群落组成在门水平和属水平上无显著差异[17]。由于生物气溶胶的研究仍非常有限,现阶段尚不能对这种空间和时间上的差异有深入的剖析。

研究也发现,沙尘天气下细菌群落会有显著变化,如2016年非洲一次强沙尘事件中,厚壁菌门和放线菌门为主要细菌类群,厚壁菌门在粗颗粒物中的相对丰度高于在细颗粒,放线菌门则相反[14],而亚洲沙尘期间,孢子类微生物成为主要的细菌类群[15]。同样,这种差异缺乏深入的科学分析。雾天相对湿度增加,悬浮着的大量小液滴可以作为微生物的载体,液滴所含的营养物质更易被微生物吸收[18]。Dong等[18]发现青岛雾天粗颗粒上的总微生物浓度增加;微生物群落结构较晴天变得更加复杂,不动杆菌(Acinetobacter)和假单胞菌(Pseudomonas)为雾天的优势菌[19]

生物气溶胶受排放源、地理环境和气象条件等因素的共同影响,不同地区存在较大空间差异[11]。由于生物气溶胶来源和时空分布的差异,不同地区特殊天气过程对空气中细菌群落的影响较难进行系统的对比研究。本文以沿海城市青岛作为研究区域,于2017年12月—2018年4月采集晴天、雾、霾及沙尘天的生物气溶胶样品,采用16S rRNA高通量测序方法,对比分析不同天气条件下生物气溶胶中细菌群落结构特征及其影响因素,初步评估污染天气发生时细菌对人体的潜在健康风险。

1 材料与方法 1.1 样品采集

采样点位于青岛市中国海洋大学崂山校区教学楼顶(36°09′N,120°29′E,见图 1)。采样高度为1.5 m,采样点距离地面垂直高度约为9 m,距离海边约为7 km,海拔约为90 m。采样点位于教学区,校园绿地面积占总面积大约50%,附近无明显人为污染排放源。于2017年12月—2018年4月,运用FA-1撞击式生物采样器(辽宁康洁仪器研究所)将大气气溶胶样品采集于0.22 μm的聚碳酸酯膜,采样流量为28.3 L/min。采样器共有6个粒径范围,为了确保不同天气下的样品能够分析测定出细菌群落,将6级合为1个样品进行分析。从出现特殊天气(雾、霾和沙尘)现象时开始采样,特殊天气结束时停止采样,采样平均时长为6 h。晴天样品从早上8点开始采集6 h。所有实验用具(采样膜、镊子、铝箔等)在121 ℃灭菌30 min,烘干后置于超净台待用。在进行DNA提取前,将采集的样品放置在约-80 ℃冰箱中储存。

图 1 大气生物气溶胶采样点示意图 Fig. 1 Location of the sampling site for atmospheric bioaerosols

采样期间PM2.5、PM10和CO、NO2、SO2、O3的浓度数据来源于青岛市环境保护局(http://www.qepb.gov.cn),温度和相对湿度(RH)数据来源于青岛市气象局(http://qdqx.qingdao.gov.cn/),数据为采样时段1 h内的平均值。采样期间气象数据如表 1所示。共采集9个生物气溶胶样品,分别为晴天样品1个(DQ1)、霾天样品6个(DQ2、DQ8、DQ4、DQ6、DQ5、DQ9)、沙尘天样品1个(DQ7)、雾天样品1个(DQ3)。其中,参照中华人民共和国气象行业标准《霾的观测和预报等级》(QX/T113-2010),将6个霾天样品划分为轻微霾天(DQ2)、轻度霾天(DQ8、DQ4、DQ6、DQ5)和中度霾天(DQ9)(能见度5~10 km,为轻微霾天;3~5 km,为轻度霾天;2~3 km,为中度霾天;小于2 km,为重度霾天)。

表 1 2017年12月—2018年4月青岛大气中颗粒物浓度和气象参数统计表 Table 1 The related parameters in atmosphere bioaerosols in Qingdao from Dec. 2017 to Apr., 2018
1.2 样品处理 1.2.1 DNA提取

吸取1 000 μL CTAB裂解液(十六烷基三甲基溴化铵)转移至2 mL离心管中,加入溶菌酶,取适量的样品加入裂解液中,在65 ℃水浴条件下裂解样品。取上清液加入离心管中,加入体积比为25∶24∶1的酚(pH=8.0)、氯仿和异戊醇溶液,颠倒混匀后12 000 r/min离心10 min。继续取上清液,加入体积比为24∶1的氯仿和异戊醇溶液,充分混匀后12 000 r/min离心10 min。移取上清液至1.5 mL离心管中,加入异丙醇,摇晃均匀,-20 ℃的环境中静置沉淀。然后12 000 r/min下离心10 min,倒出液体。用1 mL 75%乙醇洗涤2次加入双蒸水(ddH2O)溶解DNA样品,于55~60 ℃下孵育10 min助溶,然后加入1 μL核酸内切酶(RNase A)消化RNA,37 ℃静置15 min,样品置于-80 ℃的冰箱中保存。

1.2.2 基于16S rDNA的测序

以样品的DNA作为PCR反应模板,对细菌16S rRNA基因的V4-V5区进行扩增,正反向引物分别为515F:5’-GTGCCAGCMGCCGCGGTAA-3’和907R:5’-CCGTCAATTCC TTTGAGTTT-3’。PCR反应体系(30 μL):2×PhusionMasterMix 15 μL,Primer(2 μL)3 μL(6 μM),gDNA(1 ng/μL)10 μL(5~10 ng),H2O2 μL。PCR反应程序:98 ℃预变性1 min,98 ℃10 s,50 ℃30 s,72 ℃30 s扩增30个循环,72 ℃延伸5 min。PCR产物使用2%琼脂糖凝胶进行电泳检测。扩增后得到450 bp的扩增片段,检验合格后添加测序接头,IonS5TMXL测序平台进行高通量测序。

1.3 测序数据及处理

使用Cutadapt(V1.9.1)过滤处理原始序列(Raw Reads),去除嵌合体,得到有效序列。利用Uparse软件(Uparse v7.0.1001)对所有样品的全部有效序列进行聚类,默认以97%的相似性将序列聚类为OTUs(Operational Taxonomic Units)。根据OTUs结果进行物种注释和丰度、Alpha多样性计算、Venn图和花瓣图等分析,同时对OTUs进行多序列比对并构建系统发育树,探究不同样本或组间群落结构的差异。

2 结果与讨论 2.1 不同天气生物气溶胶样品中OTUs和多样性分析

对9个大气样品经高通量测序后将得到的原始数据进行拼接、质控和去嵌合体,共得到80 157条序列。根据Alpha多样性指数稀释曲线(见图 2),对每个样品随机抽取62 186条序列进行均一化处理,以97%相似度对序列进行聚类,得到1 387个OTUs。OTUs数目最多的为沙尘天样品(DQ7),有920个OTUs;其次为霾天样品,平均有(432±94)个(n=6);雾天(DQ3)有323个。有研究[20]指出,沙尘是生物气溶胶中微生物的重要来源之一,沙尘会携带大量外源微生物使下风向地区微生物浓度增加,群落结构发生变化,且粗颗粒物上的微生物更丰富[21],粗颗粒可更好的为微生物提供营养和基质,抵抗紫外线对微生物的伤害等。因此OTUs在沙尘天最高;霾天相对湿度较高、风速和压强较低,气象条件不利于污染物扩散,造成颗粒物累积,以颗粒物为载体的微生物也有所增加;而雾天颗粒物浓度低,较高的相对湿度使革兰氏阳性菌[22]外的微生物生存受到限制,雾天较低的颗粒物浓度和较高的相对湿度可能是导致雾天OTUs较低的原因。

图 2 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌Alpha多样性稀释曲线 Fig. 2 Alpha diversity index dilution curve in bioaerosol samples in different weather conditions in Qingdao from Dec. 2017 to Apr., 2018

不同程度霾天OTUs表现为:轻度霾天((479±73)个)>轻微霾天(361个)>中度霾天(315个)。OTUs在轻度霾天时最高,表明OTU并非与霾污染程度呈线性关系。在轻微和轻度污染时,硫酸盐和硝酸盐的促进作用占主导[23],可为微生物的生长提供营养。随霾污染加重,细颗粒物浓度升高,研究表明细颗粒物中的细菌比粗颗粒物中的更易受到有毒有害物质和环境的影响[24]。细颗粒物中有害物质(重金属和多环芳烃)随PM2.5浓度的升高而显著增加[23, 25],对微生物的毒性作用增强,使中度霾天OTUs减少。

所有样品共有OTUs为85个(见图 3),占OTUs总数的6.1%,表明不同天气条件下生物气溶胶样品中细菌OTUs差异较大。沙尘天(DQ7)独有的OTUs最多,为300个。沙尘传输时会携带大量沙尘源区及途径地区的微生物物种,一些高耐受性的细菌在长距离传输中可长期存活,能被沙尘有效的分散到下风向地区[18];雾天(DQ3)独有OTUs最少,雾天大气环境(颗粒物浓度低及相对湿度高)对微生物的选择性较高,新出现的物种最少;6个霾天独有OTUs居中,介于8~58个之间。4个轻度霾天独有的OTUs在28~58个之间;轻微霾天(DQ2)独有OTUs略高于中度霾天(DQ9),分别为15和8个。随霾污染的加剧,独有OTUs减少,与其他对青岛霾天OTUs的研究结果一致[7],霾发生后有适应霾天环境的新物种出现,轻度霾天新出现的物种最多,中度霾天时减少。

图 3 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌群落花瓣图 Fig. 3 Veen OTUs of bacterial community in bioaerosol samples in different weather conditions in Qingdao from Dec., 2017 to Apr., 2018

Alpha多样性指数如表 2所示。所有样品的测序覆盖度均在99.9%以上,测序深度基本覆盖全部物种。不同天气条件的Chao1指数表现为沙尘天最高,为971;4个轻度霾天次之,介于460~615之间,平均为536±76;轻微霾天高于中度霾天,分别为423和357;雾天和晴天的Chao1指数相差不大,分别为397和390。Shannon指数表现为沙尘天>轻度霾天>晴天>轻微霾天>中度霾天>雾天。沙尘天(DQ7)的Chao1指数和Shannon指数最大,表明沙尘天丰富度和多样性最高,且细菌群落最为均匀。沙尘天微生物种类增加,每个种类数量较多。一般来说,长距离传输过程中空气微生物会受到紫外辐射、营养缺乏和缺水等不利环境的考验导致其存活率降低,但微生物可通过各种防御机制抵抗不利环境,如形成休眠孢子[15],附着在颗粒物上进行传输,粗颗粒可保护细胞抵抗恶劣环境的伤害。研究表明,沙尘期间细菌丰度的增加与粗颗粒物数浓度的增加有密切关系[26]。此外,沙尘还是矿物颗粒的重要来源可为微生物提供碳源等能量源[27],故沙尘天物种丰富度和多样性都增加。雾天时颗粒物浓度低,PM2.5和PM10的浓度分别为54和99 μg/m3,微生物附着的基质减少;雾天相对湿度高达95%,对微生物的影响较为复杂。一方面,相对湿度较高时,革兰氏阳性菌容易存活[22];另一方面,较高的相对湿度可以促进颗粒物的吸湿性生长,从而促使气溶胶颗粒的干沉降[28]。雾天较低的颗粒物浓度及较高的相对湿度可能是造成细菌群落丰富度和多样性低的原因。

表 2 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌群落Alpha多样性指数 Table 2 Alpha diversity index of bacterial community in bioserosol samples in different weather conditions in Qingdao from Dec., 2017 to Apr., 2018

总体而言,霾发生后物种丰富度和多样性增加,与武丽婧[19]和邵聪聪等[29]运用变形梯度凝胶电泳(DGGE)研究青岛霾天大气细菌群落结构多样性高于晴天的结果一致。从轻微霾天到轻度霾天再到中度霾天的过程中,大气细菌群落丰富度和多样性先增加后降低,在轻度霾天时最高,细菌丰富度和多样性并不是随着霾污染的加剧而增加。低污染时,颗粒物累积导致微生物增加[18],一些化合物(硫酸盐、硝酸盐和铵)可为微生物的生长代谢提供营养[28];而随霾污染加重,二次气溶胶的化学反应增强,有助于细模态颗粒物的增加[16, 30],细颗粒物吸附的重金属、挥发性有机化合物和多环芳烃等物质会对微生物产生毒害作用。有研究表明,冬季重度霾天细颗粒物中多环芳烃和部分重金属的浓度分别为非霾天的3倍和3~8倍[31],SO42-和NO3-的浓度分别约为非霾天的17和6倍[32]。这些高浓度的有毒有害物种会导致细胞成分损伤[33],对微生物生长不利。生物气溶胶颗粒表面吸附的化学物质所引起的促生长效应和毒性效应之间的竞争关系可能是不同程度霾天细菌丰富度和多样性发生变化的原因。随霾污染的加重,化学物质的毒性效应强于促生长效应,导致中度霾天大量物种无法存活,物种多样性和丰富度均最低。

2.2 不同天气生物气溶胶中细菌群落结构组成分析

不同天气生物气溶胶样品测序共得1 387个OTUs,注释到33门,54纲,113目,204科,389属,166种。选取前十门做物种相对丰度图如图 4所示。变形菌门(Proteobacteria)、厚壁菌门(Firmicutes)、放线菌门(Actinobacteria)和拟杆菌门(Bacteroidetes)在所有样品中的相对丰度均超过1%,为相对丰度最高的4个细菌类群,相对丰度之和在56.2%~80.9%之间,这与青岛地区[19]和北疆地区[34](乌鲁木齐和石河子)细菌群落的研究一致,说明门水平上细菌群落结构具有很高的相似性,可能是由于这些细菌类群比较容易气溶胶化,能在相对不利的环境中存活[17]。但西安地区[35]秋冬霾天细菌群落优势菌为变形菌门和厚壁菌门,二者相对丰度之和在80%以上,说明不同地区微生物种类的相对占比差异较大。

图 4 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌前十门相对丰度 Fig. 4 The relative abundance of top ten in the level of phylum classification in bioserosol samples in different weather conditions in Qingdao from Dec., 2017 to Apr., 2018

任何天气下变形菌门都是相对丰度最高的细菌类群,在22.4%~55.7%之间。总体来讲,变形菌门在霾天时的相对丰度增加,其中在轻度霾天(DQ6)和中度霾天(DQ9)均达50%以上,其他天气条件下相对丰度相差不大,在22.9%~25.6%之间。变形菌门为革兰氏阴性菌,是细菌中最大的门类,分布广泛,对温度、辐射等具有较强的抵抗力[36]。厚壁菌门在9个样品中的相对丰度次之,为1.9%~33.7%。晴天时厚壁菌门处于较高水平(29.0%),在特殊天气时相对丰度有所下降,沙尘天的相对丰度比霾天和雾天的高,这可能是由于厚壁菌门以细菌聚合体形式存在,使其与大颗粒物浓度更相关[36]。放线菌门在不同类型天气的相对丰度在9.2%~18.8%之间,沙尘天和雾天的相对丰度最高,均为18.8%。有研究表明,放线菌与温度呈正相关[37]。本研究中,霾天样品采集时的平均温度为(7.8±6.0) ℃,而沙尘和雾天的平均温度分别为15.3和12.9 ℃,低温可能是霾天放线菌相对丰度较低的原因。沙尘天显著增加的有酸杆菌门(Acidobacteria),酸杆菌门在沙尘天的相对丰度为6.9%,而在其他天气下不足1%。酸杆菌门(Acidobacteria)广泛存在于各种生境中,特别是在土壤环境中,约占土壤细菌的20%~50%[38],本文的研究结果表明沙尘期间土壤来源的细菌相对丰度增加。

注释到属分类水平上,前十属物种相对丰度图如图 5所示。所有样品中前十属的相对丰度在16%~47%之间,不同天气前十属的相对丰度相差较大。仅不动杆菌(Acinetobacter)在所有样品中的相对丰度超过1%,在霾天相对丰度增加且在轻度霾天DQ4(14.1%)、DQ6(19.3%)和中度霾天DQ9(11.9%)处于较高水平,其他样品中为1.2%~2.5%。不动杆菌广泛存在于水体和土壤环境中,具有很强的适应能力,当大气中含有较多颗粒,尤其是相对湿度较高时可在颗粒物上迅速繁殖[39]。因此,不动杆菌能在霾天较为恶劣的环境中良好生存。鞘氨醇单胞菌(Sphingomonas)也是霾天的优势细菌,在轻度霾天的相对丰度较高,为1%~13.8%。马曼曼等[7]同样发现青岛地区鞘氨醇单胞菌在霾天相对丰度明显升高。鞘氨醇单胞菌具有特殊的代谢调控机制,可利用多种复杂的有机物进行生长代谢以适应环境的改变[7, 40-41],其某些菌株能够降解高分子有机污染物(尤其是多环芳烃)[41]。沙尘天的优势菌为芽孢杆菌(Bacillus),相对丰度为12.2%,约为晴天的两倍,在轻微霾天、中度霾天及雾天的相对丰度在1%左右。芽孢杆菌是大气中广泛存在的细菌,能够产生对不利条件具有特殊抵抗力的芽孢,抵抗脱水和极端环境,在很多地区都作为优势菌存在[42-45]。雾天的优势菌为微球菌(Micrococcus)和副球菌(Paracoccus),相对丰度分别为9.4%和7.6%。雾天相对丰度增加的还有假单胞菌,假单胞菌具有增强自身气溶胶化的能力,并且能在恶劣的环境中保持其完整性[46]。结果显示假单胞菌(Pseudomonas)在雾天的相对丰度从其他天气的不足1%增加到5.8%,这与武丽婧[19]使用DGGE方法的研究结果一致。总体而言,相比晴天,霾天和雾天检测出更多的致病菌,其中不动杆菌、鞘氨醇单胞菌都是条件致病菌,可通过呼吸道进入人体引发感染。

图 5 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌前十属相对丰度 Fig. 5 The relative abundance of top ten in the level of genus classification in bioserosol samples in different weather conditions in Qingdao from Dec., 2017 to Apr., 2018
2.3 不同天气生物气溶胶中细菌群落结构相似性分析

对不同天气的样品进行非加权组平均法(UPGMA)聚类分析(见图 6)。所有样品总体分为两支,第一支由轻度霾天(DQ4和DQ6)和中度霾天(DQ9)构成。反映在细菌群落结构上,这三个样品表现为一种或两种优势菌占主导的现象,细菌群落组成单一。不动杆菌在这三个样品中的相对丰度均超过12%,而在其他样品中不足3%。第二支可分为沙尘天样品(DQ7)和其他天气样品。沙尘天的样品聚类位置区别于其他特殊天样品,表明沙尘的到来可以显著改变大气细菌群落组成。而晴天DQ1和轻微霾天DQ2聚为一支,说明轻微霾天时微生物群落没有发生显著变化,微生物群落在霾天的变化与霾污染强度有关。6个霾天样品分别聚在两支,轻度霾天(DQ5和DQ8)与轻微霾天(DQ2)聚为一支。这3个霾天样品细菌群落结构均表现为物种组成丰富而相对丰度较低。另外,许多研究表明气体污染物对细菌群落结构的影响较大,其中O3与微生物呈负相关[47-48]。本研究中RDA(冗余分析)结果显示,O3是对微生物群落影响最大的因素,其次是NO2。比较霾天气体污染物浓度发现,DQ8、DQ5和DQ2的O3浓度分别为23.4、19.9和33.0 μg/m3,处于较低浓度水平。NO2的浓度也较为接近,在63.0~89.4 μg/m3之间,相较其他霾天处于较高水平。而轻度霾天DQ4和DQ6的O3浓度分别是89.6和120 μg/m3,处于较高水平。

图 6 2017年12月—2018年4月青岛不同天气生物气溶胶中细菌群落相似性聚类 Fig. 6 Cluster of bacterial community in bioserosol samples in different weather conditions in Qingdao in April from Dec., 2017 to Apr., 2018

尽管DQ8、DQ5、DQ4和DQ6同为轻度霾天样品,但在Alpha多样性、物种相对丰度和聚类结果之间存在差异。DQ8和DQ5的多样性高于DQ4和DQ6;DQ4和DQ6中优势菌占比较大,使得物种组成单一;聚类结果显示DQ4和DQ6在同一支,而DQ8和DQ5在另一支。从后向轨迹分析(距离地面1 km,72 h)。从图 7可以看出,这4个轻度霾天的气团来源可分为两类:DQ4和DQ6来自南方陆地源;DQ8和DQ5来自北方陆地源,这可能是4个轻度霾天细菌群落结构存在差异的原因。可见,霾天细菌群落结构也受到不同气团来源的影响。综上,不同程度霾天细菌群落结构的相似性是多种因素(气体污染物、气象因素和气团来源等)共同作用的结果,影响霾天细菌群落结构变化的因素还有待进一步的研究。

图 7 2017年12月—2018年4月青岛大气生物气溶胶样品72 h后向轨迹聚类 Fig. 7 Cluster results of 72 h air mass back trajectories for bioaerosol samples in Qingdao in April from Dec., 2017 to Apr., 2018
2.4 主要菌群与环境因子的关系

为了解主要菌群与大气环境因子的变化关系,对前十门细菌与大气颗粒物及污染物因子进行冗余分析(RDA),结果表明(见图 8),前两个排序轴的特征值分别为0.400和0.270。环境因子AQI(-0.106)、PM2.5(-0.120)、PM10(-0.105)、NO2(-0.166)、O3(0.290)、CO(-0.130)主要贡献于第一排序轴。两轴中细菌与环境因子的累计百分比为67.0%。所有环境因子中,O3是主要菌群的第一影响因子,其余依次为NO2、PM2.5、CO、AQI和PM10。主要细菌群中,变形菌门与CO、PM2.5和NO2呈正相关,厚壁菌门与AQI和PM10呈正相关,酸杆菌门和拟杆菌门与O3呈正相关。

图 8 环境因子与微生物群落(门水平)的冗余分析 Fig. 8 Redundancy analysis to explore relationship between microorganism community (phylum) and environmental factors

本研究中变形菌门与PM2.5呈正相关,微生物附着在颗粒物上,颗粒物增加会使微生物浓度增加[18]。Qi等[24]研究表明变形菌门的相对丰度在霾天显著改变,尤其是在细颗粒物上。Yu等[49]发现霾天变形菌门在细颗粒物上的相对丰度为88.9%,而在粗颗粒物上的相对丰度为35.7%。与细颗粒的相关性可能是变形菌门在霾天相对丰度较大,且在中度霾天时依然能维持较高相对丰度的原因;另一优势菌厚壁菌门与PM10的相关性更大。Yu等[49]也发现霾天粗粒径上厚壁菌门相对丰度为16.5%,细粒径上为4.4%。这可能是由于厚壁菌门以细菌聚合体形式存在,使其与大颗粒物更相关[36]。放线菌门与O3有很好的正相关性。高浓度的O3会使细菌细胞氧化破裂,细胞膜结构紊乱,渗透性增强使膜活性受阻,最终导致细胞死亡[50],但适宜浓度的O3可能会与颗粒物上的其他物质发生化学反应,从而减少了对微生物的伤害[51],尤其放线菌可以在极端温度和pH、营养贫瘠的环境中存活[52]。但目前气体污染物对微生物的影响和机制尚不清楚,还需更深入的研究。

2.5 健康风险分析

研究表明致病菌比其他细菌更能承受空气中的毒素,当空气污染严重时病原菌仍可以生存,而其他细菌则可能死亡,因此可能导致致病菌比例的增加[53]。相比晴天,特殊天(沙尘天、霾天和雾天)能检测到更多致病菌或条件致病菌,可能会对人体健康造成潜在影响。霾天出现的优势菌不动杆菌是医院感染的主要细菌,广泛分布于水体和土壤中,也存在于人体皮肤、口腔和呼吸道中,是唯一能在人体皮肤表面生存的革兰氏阴性杆菌[54],主要引起人体呼吸道感染,还能引起中耳炎、脑膜炎、烧伤感染、败血症和泌尿系统感染[55]。霾天的另一优势菌鞘氨醇单胞菌广泛分布于自然界中,能引起植物根部和动物伤口感染,会通过感染的伤口进入人体而引起各种炎症、败血症等疾病[39]。雾天假单胞菌和微球菌的相对丰度显著增加。假单胞菌是条件致病菌,是引起医院感染的常见病原菌。假单胞菌包含很多类别,其中的铜绿假单胞菌(Pseudomonas aeruginosa)是引起医院感染的三大病原菌之一,可在潮湿环境中大量繁殖,能引起皮肤感染、呼吸道感染、泌尿道感染和囊性纤维变性等[56]。铜绿假单胞菌感染的发病率正逐年升高,医院感染日益严重[57];微球菌广泛存在于土壤、水体及人体皮肤和呼吸道中,一般不致病,但在机体抵抗力低下时则可致病,但近年来对发热病例的长期观察发现,微球菌的发病率和致病性有逐年增高的趋势[58]

3 结语

对青岛2017年12月—2018年4月期间不同天气条件大气细菌群落结构的16S rRNA基因高通量分析结果表明,特殊天气(沙尘天、霾天和雾天)发生时细菌群落组成和多样性变化显著。沙尘天物种丰富度和多样性最高,雾天最低,霾天丰富度和多样性与霾污染程度有关,污染加重,群落多样性降低。对细菌群落结构影响最大的环境因子是O3,同时还受到颗粒物浓度和气团来源的影响。沙尘天、霾天和雾天优势菌属差异较大,出现不同的条件致病菌,会对人体健康造成潜在风险。由于样品有限,不同天气对细菌群落结构的影响机制并不清楚,有待进一步研究。

参考文献
[1]
Ariya P A, Amyot M. New directions: The role of bioaerosols in atmospheric chemistry and physics[J]. Atmospheric Environment, 2004, 38(8): 1231-1232. DOI:10.1016/j.atmosenv.2003.12.006 (0)
[2]
Lange J L, Thorne P S, Lynch N L. Application of flow cytometry and fluorescent in situ hybridization for assessment of exposures to airborne bacteria[J]. Applied and Environmental Microbiology, 1997, 63(4): 1557-1563. DOI:10.1128/AEM.63.4.1557-1563.1997 (0)
[3]
Bauer H, Heidi G. Airborne bacteria as cloud condensation nuclei[J]. Journal of Geophysical Research, 2003, 108(21): 4658-4661. (0)
[4]
Albrecht A, Witzenberger R, Bernzen U, et al. Detection of airborne microbes in a composting facility by cultivation based and cultivation-independent methods[J]. Annals of Agricultural and Environmental Medicine: AAEM, 2007, 14(1): 81-85. (0)
[5]
Bowers R M, Lauber C L, Wiedinmyer C, et al. Characterization of airborne microbial communities at a high-elevation site and their potential to act as atmospheric ice nuclei[J]. Applied & Environmental Microbiology, 2009, 75(15): 5121-5130. (0)
[6]
Cao C, Jiang W, Wang B, et al. Inhalable microorganisms in Beijing's PM2.5 and PM10 pollutants during a severe smog event[J]. Environmental Science & Technology, 2014, 48(3): 1499-1507. (0)
[7]
马曼曼, 甄毓, 米铁柱, 等. 青岛冬季霾天不同粒径生物气溶胶中细菌群落特征研究[J]. 中国环境科学, 2017, 37(8): 2855-2865.
Ma M M, Zhen Y, Mi T Z, et al. Bacterial community characteristics in different particle sizes of bioaerosols in winter haze days in Qingdao[J]. China Environment Science, 2017, 37(8): 2855-2865. DOI:10.3969/j.issn.1000-6923.2017.08.007 (0)
[8]
Schwartz D A, Thorne P S, Yagla S J, et al. The role of endotoxin in grain dust-induced lung disease[J]. American Journal of Respiratory & Critical Care Medicine, 1995, 152(2): 603-608. (0)
[9]
Husman T. Health effects of indoor-air microorganisms[J]. Scandinavian Journal of Work, Environment & Health, 1996, 22(1): 5-13. (0)
[10]
Ariya A P, Amyot M. New directions: The role of bioaerosols in atmospheric chemistry and physics[J]. Atmospheric Environment, 2004, 38(8): 1231-1232. DOI:10.1016/j.atmosenv.2003.12.006 (0)
[11]
Li Y P, Fu H L, Wang W, et al. Characteristics of bacterial and fungal aerosols during the autumn haze days in Xi'an, China[J]. Atmospheric Environment, 2015, 122: 439-447. DOI:10.1016/j.atmosenv.2015.09.070 (0)
[12]
Gao M, Jia R Z, Qiu T L, et al. Seasonal size distribution of airborne culturable bacteria and fungi and preliminary estimation of their deposition in human lungs during non-haze and haze days[J]. Atmospheric Environment, 2015, 18(4): 203-210. (0)
[13]
Liu H, Zhang X, Zhang H, et al. Effect of air pollution on the total bacteria and pathogenic bacteria in different sizes of particulate matter[J]. Environmental Pollution, 2018, 233: 483-493. DOI:10.1016/j.envpol.2017.10.070 (0)
[14]
Tang K, Huang Z, Huang J, et al. Characterization of atmospheric bioaerosols along the transport pathway of Asian dust during the Dust-Bioaerosol 2016 Campaign[J]. Atmospheric Chemistry & Physics, 2018, 18(10): 7131-7148. (0)
[15]
Hara K, Zhang D. Bacterial abundance and viability in long-range transported dust[J]. Atmospheric Environment, 2012, 47: 20-25. DOI:10.1016/j.atmosenv.2011.11.050 (0)
[16]
Wei M, Liu H, Chen J, et al. Effects of aerosol pollution on PM2.5-associated bacteria in typical inland and coastal cities of northern China during the winter heating season[J]. Environmental Pollution, 2020, 262: 114-118. (0)
[17]
王步英, 郎继东, 张丽娜, 等. 基于16S rRNA基因测序法分析北京霾污染过程中PM2.5和PM10细菌群落特征[J]. 环境科学, 2015, 36(8): 2727-2734.
Wang B Y, Lang J D, Zhang L N, et al. Characterizing Beijing's communities in PM2.5 and PM10 samples during haze pollution episodes using 16S rRNA gene analysis method[J]. Environmental Science, 2015, 36(8): 2727-2734. (0)
[18]
Dong L, Qi J H, Shao C C, et al. Concentration and size distribution of total airborne microbes in hazy and foggy weather[J]. Science of the Total Environment, 2015, 541: 1011-1018. (0)
[19]
武丽婧. 青岛及黄海生物气溶胶中微生物群落多样性研究[D]. 青岛: 中国海洋大学, 2014.
Wu L J. Microbial Community Diversity of Bioaerosols in Qingdao and the Yellow Sea[D]. Qingdao: Ocean University of China, 2014. (0)
[20]
Jeon E M, Kim H J, Jung K, et al. Impact of Asian dust events on airborne bacterial community assessed by molecular analyses[J]. Atmospheric Environment, 2011, 45(25): 4313-4321. DOI:10.1016/j.atmosenv.2010.11.054 (0)
[21]
Bowers R M, Clements N, Emerson J B, et al. Seasonal variability in bacterial and fungal diversity of the near- surface atmosphere[J]. Environmental Science & Technology, 2013, 47(21): 12097-12106. (0)
[22]
Theunissen H J H, Toom L D, Burggraaf A, et al. Influence of temperature and relative humidity on the survival of Chlamydia pneumonia in aerosols[J]. Applied & Environmental Microbiology, 1993, 59(8): 2589-2593. (0)
[23]
Chen J, Qiu S S, Shang J, et al. Impact of relative humidity and water soluble constituents of PM2.5 on visibility impairment in Beijing, China[J]. Aerosol and Air Quality Research, 2014, 14(1): 260-268. DOI:10.4209/aaqr.2012.12.0360 (0)
[24]
Qi J H, Li M Z, Zhen Y, et al. Characterization of bioaerosol bacterial communities during hazy and foggy weather in Qingdao, China[J]. Journal of Ocean University of China, 2018, 17(3): 71-81. (0)
[25]
Kong S, Li X, Li L, et al. Variation of polycyclic aromatic hydrocarbons in atmospheric PM2.5 during winter haze period around 2014 Chinese Spring Festival at Nanjing: Insights of source changes, air mass direction and firework particle injection[J]. Science of the Total Environment, 2015, 520: 59-72. DOI:10.1016/j.scitotenv.2015.03.001 (0)
[26]
Lee S, Choi B, Yi S M, et al. Characterization of microbial community during Asian dust events in Korea[J]. Science of the Total Environment, 2009, 407(20): 5308-5314. DOI:10.1016/j.scitotenv.2009.06.052 (0)
[27]
Kim K Y, Kim C N. Airborne microbiological characteristics in public buildings of Korea[J]. Building & Environment, 2007, 42(5): 2188-2196. (0)
[28]
Liu Z, Hu B, Wang L, et al. Seasonal and diurnal variation in particulate matter (PM10 and PM2.5) at an urban site of Beijing: analyses from a 9-year study[J]. Environmental Science and Pollution Research, 2015, 22(1): 627-642. DOI:10.1007/s11356-014-3347-0 (0)
[29]
邵聪聪, 甄毓, 米铁柱, 等. 青岛冬季霾天与非霾天大气细菌群落结构特征研究[J]. 城市环境与城市生态, 2016, 29(3): 7-13.
Shao C C, Zhen Y, Mi T Z, et al. Characterization of airborne bacteria communities in Qingdao during haze and non-haze episodes in winter[J]. Urban Environment and Urban Ecology, 2016, 29(3): 7-13. (0)
[30]
Zhang T, Li X Y, Chen H X, et al. Microbial aerosol chemistry characteristics in highly polluted air[J]. Science China Chemistry, 2019, 62(8): 1051-1063. DOI:10.1007/s11426-019-9488-3 (0)
[31]
Xie Z, Li Y, Lu R, et al. Characteristics of total airborne microbes at various air quality levels[J]. Journal of Aerosol Science, 2018, 116. (0)
[32]
Guo Z, Wang Z, Zhao Z, et al. Biological and chemical compositions of atmospheric particulate matter during hazardous haze days in Beijing[J]. Environmental Science and Pollution Research, 2018, 25(34): 34540-34549. DOI:10.1007/s11356-018-3355-6 (0)
[33]
Xu W, Chen H, Li D, et al. A case study of aerosol characteristics during a haze episode over Beijing[J]. Procedia Environmental Sciences, 2013, 18: 404-411. DOI:10.1016/j.proenv.2013.04.054 (0)
[34]
Wei H Y, Zhang Y Y, Song S J, et al. Alveolar macrophage reaction to PM2.5 of hazy day in vitro: Evaluation methods and mitochondrial screening to determine mechanisms of biological effect[J]. Ecotoxicology and Environmental Safety, 2019, 174: 566-573. DOI:10.1016/j.ecoenv.2019.02.093 (0)
[35]
苟欢歌, 谢纯斌, 童延斌, 等. 石河子市春季大气颗粒物TSP和PM10中微生物群落特征[J]. 环境工程学报, 2017, 11(4): 2342-2349.
Gou H G, Xie C B, Tong Y B, et al. Assessment of microbial communities in TSP and PM10 of Shihezi during spring[J]. Chinese Journal of Enivironmental Engineering, 2017, 11(4): 2342-2349. (0)
[36]
宋颖. 西安市秋冬季微生物气溶胶中细菌群落结构特征研究[D]. 西安: 长安大学, 2017.
Song Y. Bacterial Community Structure of Atmospheric Aerosol in autumn and winter in Xi'an[D]. Xi'an: Chang'an University, 2017. (0)
[37]
肖桃. 变形杆菌比较基因组学分析及耐药基因岛的研究[D]. 北京: 中国疾病预防中心, 2019.
Xiao T. Comparative Genomics Analysis of Proteobacteria and Study on Drug Resistance Gene Island[D]. Beijing: Chinese Center for Disease Control and Prevention, 2019. (0)
[38]
Franzetti A, Gandolfi I, Gaspari E, et al. Seasonal variability of bacteria in fine and coarse urban air particulate matter[J]. Applied Microbiology & Biotechnology, 2011, 90(2): 745-753. DOI:10.1007/s00253-010-3048-7 (0)
[39]
王光华, 刘俊杰, 于镇华, 等. 土壤酸杆菌门细菌生态学研究进展[J]. 生物技术通报, 2016, 32(2): 14-20.
Wang G H, Liu J J, Yu Z H, et al. Research progress of Acidobacteria ecology in soils[J]. Biotechnnology Bulletin, 2016, 32(2): 14-20. (0)
[40]
Doughari H J, Ndakidemi P A, Human I S, et al. The ecology, biology and pathogenesis of Acinetobacter spp.: An overview[J]. Microbes & Environments, 2011, 26(2): 101-112. (0)
[41]
胡杰, 何晓红, 李大平, 等. 鞘氨醇单胞菌研究进展[J]. 应用与环境生物学报, 2007, 13(3): 431-437.
Hu J, He X H, Li D P, et al. Progress in research of Sphingomonas[J]. Chinese Journal of Applied and Enivironmental Biology, 2007, 13(3): 431-437. DOI:10.3321/j.issn:1006-687X.2007.03.030 (0)
[42]
Fegatella F, Cavicchioli R. Physiological responses to starvation in the marine oligotrophic ultramicro bacterium Sphingomonas sp. Strain RB2256[J]. Applied and Environmental Microbiology, 2020, 66(5): 2037-2044. (0)
[43]
Kim KY, Kim C N. Airborne microbiological characteristics in public buildings of Korea[J]. Building and Environment, 2007, 42(5): 2188-2196. DOI:10.1016/j.buildenv.2006.04.013 (0)
[44]
Pastuszka J S, Paw U K T, Lis D O, et al. Bacterial and fungal aerosol in indoor environment in Upper Silesia, Poland[J]. Atmospheric Environment, 2000, 34(22): 3833-3842. DOI:10.1016/S1352-2310(99)00527-0 (0)
[45]
Smith D J, Griffin D W, Jaffe D A. The high life: Transport of microbes in the atmosphere[J]. Eos Transactions American Geophysical Union, 2011, 92(30): 249-250. DOI:10.1029/2011EO300001 (0)
[46]
方治国, 欧阳志云, 胡利锋, 等. 城市生态系统空气微生物群落研究进展[J]. 生态学报, 2004, 24(2): 315-322.
Fang Z Z, Ouyang Z Y, Hu L F, et al. Progresses of airborne microbial communities in urban ecosysytem[J]. Acta Ecologica Sinica, 2004, 24(2): 315-322. DOI:10.3321/j.issn:1000-0933.2004.02.023 (0)
[47]
韩晨, 谢绵测, 祁建华, 等. 青岛市不同空气质量下可培养生物气溶胶分布特征及影响因素[J]. 环境科学研究, 2016, 29(9): 1264-1271.
Han C, Xie M C, Qi J H, et al. Size distribution characteristics of culturable bioaerosols in relation to air quality levels in Qingdao[J]. Research of Environmental Sciences, 2016, 29(9): 1264-1271. (0)
[48]
Du P, Du R, Ren W, et al. Variations Variations of bacteria and fungi in PM2. 5 in Beijing, China[J]. Atmospheric Environment, 2018, 172: 55-64. DOI:10.1016/j.atmosenv.2017.10.048 (0)
[49]
Yu R L, Wang S K, Wu X L, et al. Community structure variation associated with airborne particulate matter at central south of China during hazy and nonhazy days[J]. Atmospheric Pollution Research, 2020. (0)
[50]
Bowers R M, Clements N, Emerson J B, et al. Seasonal variability in bacterial and fungal diversity of the near-surface atmosphere[J]. Environmental Science&Technology, 2013, 47(21): 12097-12106. (0)
[51]
孙东平, 李红召. 臭氧对细菌杀灭作用的初步研究[J]. 生物学杂志, 1999, 16(5): 19-20.
Sun D P, Li H Z. Study on the sterilization effects of bacteria by the trement of ozone[J]. Joural of Biology, 1999, 16(5): 19-20. (0)
[52]
Araujo R, Gupta V V S R, Reith F, et al. Biogeography and emerging significance of Actinobacteria in Australia and Northern Antarctica soils[J]. Soil Biology and Biochemistry, 2020, 146. (0)
[53]
Liu H, Zhang X, Zhang H, et al. Effect of air pollution on the total bacteria and pathogenic bacteria in different sizes of particulate matter[J]. Environmental Pollution, 2018, 233: 483-493. DOI:10.1016/j.envpol.2017.10.070 (0)
[54]
陈静, 张才成, 邹海虹, 等. 烧伤患者不动杆菌感染及耐药性分析[J]. 实用医学杂志, 2012, 28(1): 138-140.
Chen J, Zhang C C, Zhou H H, et al. Analysis of Acinetobacter infection and drug resistance in urban patients[J]. Journal of Practial Medicine, 2012, 28(1): 138-140. DOI:10.3969/j.issn.1006-5725.2012.01.059 (0)
[55]
王潭枫, 符振华. 医院感染不动杆菌和肠球菌的分离及耐药性调查分析[J]. 中国微生态学杂志, 2011, 23(11): 138-140.
Wang T F, Fu Z H. Analysis on distribution and antibiotic resistance of nosocomial infection of Acinetobacter and Enterococcus[J]. Chinese Journal of Micorecology, 2011, 23(11): 138-140. (0)
[56]
Fahlgren C, Bratbak G, Sandaa R A, et al. Diversity of airborne bacteria in samples collected using different devices for aerosol collection[J]. Aerobiologia, 2011, 27(2): 107-120. DOI:10.1007/s10453-010-9181-z (0)
[57]
魏金花. 铜绿假单胞菌PA2800在耐药性和致病性中的作用研究[D]. 西安: 西北大学, 2010.
Wei J H. The Roles of PA2800 in Antibiotic Resistande and Virulence in Pseudomonas aeruginosa[D]. Xi'an: Northwest University, 2010. (0)
[58]
李钧, 杜生祥. 小儿微球菌性败血症36例临床分析[J]. 新乡医学院学报, 1989(4): 264-265.
Li J, Du S X. Clinical analysis of 36 cases of Microcococcal septicemia in children[J]. Journal of Xinxiang Medical College, 1989(4): 264-265. (0)
Bacterial Community Characteristics in Different Weather Conditions of Atmospheric Bioaerosols in Qingdao
WANG Yao1,2 , QI Jian-Hua1,2 , ZHANG Ting3 , ZHANG Da-Hai3     
1. Key Laboratory of Marine Environmental Science and Ecology, Ministry of Education, Ocean University of China, Qingdao 266100, China;
2. Laboratory for Marine Ecology and Environmental Science, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266237, China;
3. College of Chemical Engineering, Ocean University of China, Qingdao 266100, China
Abstract: In order to evaluate the impact of foggy, hazy and dusty weather on the bacterial communities in atmospheric bioaerosols, bioaerosol samples were collected in the coastal region of Qingdao on different weather from Dec. 2017 to Apr. 2018. Bioaerosol bacterial community compositions were measured using 16S rRNA high-through sequencing. The results showed that the compositions of bacterial community were significant varied on the hazy and dusty weather, as indicated by the higher richness and diversity of community than those on the sunny weather, and have a negative correlation with the aggravation of haze pollution, but the diversity decreased in foggy days while the richness almost kept the same. Air mass source had an effect on bacterial communities in hazy days. Proteobacteria was the dominant bacteria on any weather conditions, with the highest relative abundance on hazy days. At the genus level, Actinetobacter and Sphingomonas were the dominant genera during the period of haze, and the relative abundance of Micrococcus, Paracoccus and Pseudomonas were increased on foggy weather, while Bacillus was the dominant genus on dusty days. RDA analysis showed that O3 was the major environmental factor that could affect bacterial community. Opportunistic pathogens could be detected in foggy, hazy and dusty weather, which had potential harm to human health.
Key words: bioaerosols    bacteria    16S rRNA    community structure    health risks