Structural Characteristics and Functional Implications of PM2.5 Bacterial Communities During Fall in Beijing and Shanghai, China

Air pollution characterized by fine particulate matter (PM_2.5) frequently has occurred in China, and has posed threats to human health. The physiochemical characteristics of airborne PM_2.5 have been extensively studied, but its bacterial structures and functions have not yet been well studied. Herein, we focused on the structural characteristics and functional implications of airborne bacteria under different pollution levels in Beijing and Shanghai. The α- and β-diversities showed no obvious difference in two cities (p > 0.05). The dominant phyla Proteobacteria, Firmicutes, and Actinobacteria with total abundance of over 92% were found in all PM_2.5 samples. The results of weighted unifrac non-metric multidimensional scaling (NMDS) suggested that air pollution was no obviously correlated with bacterial community but dispersed disorderly. Furthermore, canonical correlation analysis (CCA) and permutation test indicated that ${\text{NH}}_4^{+}$, ${\text{SO}}_4^{2-}$, and wind speed were the key factors that associated with airborne bacterial community structure. Chemical components of particulate matter played more important role in structuring bacterial community than meteorological conditions based on the result of partial CCA. In addition, the annotation of metabolic pathway suggested that the predominant genus Pseudomonas was obviously correlated with disease infections. Several dominant species might contribute to organic degradation, nitrogen cycles, and ice-nuclei activities in environments. Overall, this work enhanced our understanding of functions of airborne bacteria and highlighted their potential role in atmospheric chemical progresses.

Introduction

Severe urban air pollution caused by particulate matters is a hot environmental issue in recent years due to its great potential influences on atmospheric and terrestrial ecosystems (Delort et al., 2010; Pandey et al., 2014). Specially, it induces human health risks such as cardiovascular diseases, respiratory infections, and lung cancer (Esposito et al., 2012; Pope and Dockery, 2013; Maclntyre et al., 2014). Previously, extensive studies have been carried out in studying chemical composition of particulate matter, which are considered to be the key factors of these harmful influences (Sun et al., 2013; Chen et al., 2017). Actually, the airborne microorganisms also might pose threat to human, animals, and plants. Recent studies showed that microorganisms were one of main components in the airborne particulate matters, occupying about 10% of total aerosol particles (Wei et al., 2015). Among these microorganisms, bacteria accounted for >80% (Cao et al., 2014), with the abundance from 10⁵ to 10⁶ cells/m³ in near-surface atmosphere (Bowers et al., 2011; Wei et al., 2016). For one thing, the sizes of bacterial and fungal cells were reported to be ranged from 1.1 to 2.1 μm and from 2.1 to 3.3 μm aerodynamic diameters, respectively (Raisi et al., 2013), easily penetrating and depositing into tracheobronchial and alveolar regions (Brook et al., 2004). For another thing, several airborne bacteria could be the sources of pathogens and allergens for human or plant (Woo et al., 2013; Cao et al., 2014; Du et al., 2018). In addition, studies found that bacteria in clouds also could act as cloud condensation nuclei or ice nuclei and then influence the cloud and rain formation (Bauer et al., 2003; Amato et al., 2015; Hara et al., 2016). Even some of these bacteria were viable and could metabolize ammonia oxidization (Gao et al., 2016) or biotransformation of organic compounds (Delort et al., 2010; Vaïtilingom et al., 2010, 2013; Maki et al., 2018). Therefore, besides chemical composition, it is essential to understand the bacterial characteristics of airborne particles.

In regard to the relationships between pollution levels and airborne bacteria, most existing studies focused on the impact of pollution on bacterial cell concentration. Although both positive (Wei et al., 2016; Li et al., 2018; Lu et al., 2018) and negative correlations (Gao et al., 2015; Xu et al., 2017) between pollution levels and airborne bacterial concentrations were observed, it confirmed that air pollution influenced the airborne bacterial concentration. Besides, Zhong et al. (2016) found that microbial metabolic activities in particulate matter in hazy days were higher than that in sunny days. However, limited studies have been conducted for studying shifts of airborne bacterial community under different air pollution levels (Wei et al., 2016, 2017; Li et al., 2018). The results were inconsistent. For instance, Wei et al. (2016) indicated that air pollution posed no effect on bacteria structure due to the stable sources during sampling period, while Wei et al. (2017) and Li et al. (2018) found the opposite results. Wei et al. (2017) also found that potential pathogens such as Empedobacter brevis and Staphylococcus equorum were enriched in the polluted cloud. Thus, in order to understand effects of airborne bacteria on human health and ecosystem, more studies about effects of air pollution on airborne bacteria should be conducted.

In previous studies, meteorological parameters (relative humidity, wind speed, temperature) and gaseous components (NO₂, SO₂, O₃, CO) are often discussed as the main factors affecting the airborne bacterial diversity and structure (Woo et al., 2013; Lee et al., 2017; Xu et al., 2017). Nevertheless, airborne particulate matter included abundant organic matters, nitrate, and sulfate (Guo et al., 2014). Although Woo et al. (2013) showed the bioaerosol composition were determined by climate rather than airborne pollutants, Wei et al. (2017) found major ions were vital for shifting bacterial structures in clouds. It suggested that both climate and chemical components in PM_2.5 might affect the bacteria structure.

Beijing and Shanghai are two of the largest cities in China, and have more than 10 million people. Beijing is a hot city for particle pollution study because it has suffered from frequent polluted events (Cao et al., 2014; Wei et al., 2016), for example, about one-third of the fall of 2017 with the concentration of PM_2.5 > 75 μg/m³. Although the air pollution in Shanghai is relative mild, air pollution occasionally happened such as November 7–8, 2017 with an average PM_2.5 concentration of 98.5 μg/m³. Furthermore, few studies were conducted to study airborne bacteria in Shanghai. In this study, PM_2.5 samples were collected from these two megacities in polluted and unpolluted days in fall of 2017. The relation between air pollution and bacterial structure in atmospheric particulate matter was discussed; the dominant factors influencing airborne bacterial community were explored; and the key potential functions of these airborne bacteria were predicted.

Materials and Methods

Collection of PM_2.5 Samples

PM_2.5 samples were collected from the roof top of buildings at Chinese Research Academy of Environmental Sciences (Located in Beijing, 40.04° N, 116.41° E) and Shanghai Jiao Tong University (Located in Shanghai, 31.02° N, 121.43° E), respectively, where there were no major pollution sources nearby. Sampling was conducted by two high-volume samplers (KC-1000, Qingdao Laoshan Electric Inc., China), both equipped with PM_2.5 fractionating inlets. Sampling was operated in fall (2017.09–2017.11) with different pollution levels. The detailed sampling days were listed in Table 1. Each sample was collected at a flow rate of 1.05 m³/min for 23.5 h (9:00 am to 9:00 am the next day) onto a Tissuquartz membrane filter (PALL Life Sciences, 2500QAT-UP, 8 in × 10 in). Blank samples were collected similar to the sampling methods except without instrument operation. Prior to sampling, the membrane filters were sterilized at 450°C for 4 h in a Muffle furnace. One part of the sample each day was used to extract water soluble ions, and another part of sample was stored at −80°C until genomic DNA extraction.

TABLE 1

Table 1. Mean mass concentration of PM_2.5 gas pollutants and meteorological parameters during sampling period.

The environmental pollutant indexes (PM_2.5, SO₂, NO₂, CO, and O₃) and meteorological parameters (wind speed, temperature, and relative humility) per hour were obtained from the environmental monitoring stations at Beijing Olympic Sports Center and Shanghai Normal University, respectively, which were the nearest monitoring stations to sampling sites.

Extraction and Measurement of Water-Soluble Ions

One-eighth of membrane filter (4 in × 2.5 in) per day was cut into a 50-mL centrifuge tube and then submerged with 50 mL deionized water. After ultrasonication extraction for 15 min at room temperature, the extract was filtered with 0.45 μm Teflon filter. Dionex Aquion (ICS-90, Thermo Scientific) were used to analyze the concentrations of water soluble cations (${\text{NH}}_4^{+}$, Na⁺, K⁺, Ca²⁺, and Mg²⁺), anions (${\text{NO}}_3^{-}$, ${\text{SO}}_4^{2-}$, F^–, and Cl^–) and organic acids (HCOOH, CH₃COOH, C₂H₂O₄, C₄H₆O₄, C₅H₈O₄, CH₄O₃S) in the filtrate. The pH was detected with a precision pH meter (PHS-3C, Leici Instruments Co., Shanghai, China).

DNA Extraction and Sequencing

Considering the low DNA yield of airborne sample, one-eighth aerosol samples (about 185 m³ air) from several days were combined as one sample for the DNA extraction according to the method of Du et al. (2018) (see Table 1 for details). After the membrane filters were fully ground to powder in liquid nitrogen, 0.5 g sample was weighed to extract the total bacterial genomic DNA using the E.Z.N.A. Soil DNA Kit (Omega Bio-Tek, Norcross, GA, United States) according to the manufacturer’s instructions. To check sampling contamination, DNA of blank samples was extracted as samples. The extracted DNA was PCR amplified. No target fragments were observed for blank samples.

The bacterial 16S rRNA genes of aerosol samples were amplified using the PCR primers 341F (CCTACGGGNGGCW GCAG) and 806R (GGACTACHVGGGTATCTAAT) targeting the V3–V4 region (95°C for 2 min, followed by 27 cycles at 98°C for 10 s, 62°C for 30 s, 68°C for 30 s, and a final extension at 68°C for 10 min), where the eight-base barcodes were randomly added to the upstream of the universal primer to distinguish the different samples. PCR reactions were performed: 50 μL mixture containing 5 μL of 10× KOD buffer, 5 μL of 2.5 mM dNTPs, 1.5 μL of each primer (5 μM), 1 μL of KOD polymerase, and 100 ng of template DNA. After purification and quantification control, the V3–V4 tag PCR products were pooled with the other samples in equimolar and sequenced using 300 bp paired-end model on an Illumina HiSeq2500 platform in the GENE DENOVO Co., Ltd. (Guangzhou, China).

Sequence Data Analysis

After sequencing, raw data were filtered according to the following two rules: (1) removing reads containing >10% of unknown nucleotides (N) and (2) removing reads containing <80% of base with quality (Q-value > 20). The paired end clean reads were merged as raw tags by the FLASH software (v 1.2.11) (Magoč and Salzberg, 2011). Low quality fragments were filtered out using QIIME software (v 1.9.1). Chimera checking and removal was performed via UCHIME algorithm on mothur platform. The effective tags were clustered into OTUs of ≥97% similarity using UPARSE software (Edgar, 2013). Singletons were removed from the whole sequence data set. The α-diversities (Chao1, Simpson, Shannon, Ace, and goods_coverage) were calculated in QIIME (v 1.9.1). The dissimilarity test [non-metric multidimensional scaling (NMDS)] based on Bray–Curtis similarity distance matrices was performed by the Vegan package in R 3.4.1. To acquire bacterial community function, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt, v 1.0) was performed (Langille et al., 2013): firstly, predicting gene content for each OTU sequence based on Greengenes database; secondly, finding the reference genomes whose metabolic pathway has been annotated; and then summarizing their Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways according to the OTU abundance. The dominant OTU sequences were also blasted in the NCBI database to find the potentially corresponding strains and found out their ecological functions according to previous reports. And the canonical correlation analysis (CCA) was used to demonstrate the correlation between bacterial community and environmental factors. Monte Carlo permutation test was used to verify the significance of environmental factors related to bacterial community. Furthermore, the partial CCA (pCCA) was performed to compare the contributions of nutrients and meteorological conditions on the airborne bacterial structure (Liu et al., 2013). Raw 16S rRNA gene sequences reported here were deposited into the NCBI Sequence Read Archive (SRA¹) database (Accession Number: SRP167452).

Results

Definition and General Characteristics of Unpolluted and Polluted Days

In the present study, aerosols were sampled from two metropolis Beijing and Shanghai during fall. First, referring to the criteria of China Meteorological Administration (CMA), we divided samples into two levels: unpolluted (0 < PM_2.5 < 75 μg/m³) and polluted (PM_2.5 > 75 μg/m³) groups (Table 1). Next, combined with the water soluble ions, the classification of PM_2.5 samples was confirmed by cluster analysis in Figure 1. Therefore, we categorized samples SH-1, SH-2, SH-3, SH-4, BJ-1, and BJ-2 as unpolluted samples and samples SH-5, BJ-3, BJ-4, and BJ-5 as polluted samples.

FIGURE 1

Figure 1. Heatmap based on concentrations of water soluble components in all PM_2.5 samples. Red color indicated the higher concentration compared to other samples, while green color indicated the lower concentration compared to other samples. Cluster was also analyzed for all samples.

Similar chemical components were observed for all samples, and ${\text{NH}}_4^{+}$, ${\text{NO}}_3^{-}$, and ${\text{SO}}_4^{2-}$ (SNA) were dominant in water soluble extraction. From the point of city, there was no obvious difference for the concentrations of these water soluble ions in two cities due to the large concentration differences within group. However, air pollution levels influenced their concentrations (Figure 1). The sum of SNA increased as the PM_2.5 concentrations (R = 0.90, p < 0.01), which occupied 90.66 ± 0.03% of total measured water soluble matters in polluted days, significantly higher than that in unpolluted days (76.34 ± 0.07%). The concentrations of ${\text{NH}}_4^{+}$ (12.71 μg/m³), ${\text{NO}}_3^{-}$ (34.68 μg/m³), and ${\text{SO}}_4^{2-}$ (11.91 μg/m³) in polluted days were two to four times higher than that in unpolluted days (p < 0.01). Besides these three major ions, several organic acids were also found to be obviously enriched in polluted days such as C₂H₂O₄, C₄H₆O₄, C₅H₈O₄, and CH₄O₃S (Figure 1). Compared with unpolluted days, motor vehicle exhaust like CO, NO₂ in polluted days were obviously higher with 1.40 and 70.08 μg/m³, respectively (Table 1). Wind speed (1.71 m/s) in polluted days was lower than that (2.53 m/s) in unpolluted days.

Overview of Bacterial Communities in PM_2.5

More than 60,000 raw sequences were obtained for each sample, among which there were 80.71–85.62% effective sequences after quality control. According to Shannon rarefaction curve (Supplementary Figure S1), there were enough sequencing numbers for further analyses of airborne bacterial communities. These effective sequences resulted in 206–567 operational taxonomic units (OTUs) at 97% sequence identify cutoff (Table 2). T-test results showed that no obvious difference was observed for α-diversity between two cities (p > 0.05). Similarly, there was no significant difference for α-diversity between polluted days and unpolluted days (p > 0.05). Pearson correlation analysis also confirmed that PM_2.5 concentration was no obviously correlated with α-diversity (Supplementary Table S1).

TABLE 2

Table 2. Summary of sequence numbers and α-diversity of atmospheric bacterial community in PM_2.5 obtained by 16S rRNA genes sequencing on the illumina Hiseq platforms.

The NMDS result at OTU level further showed bacterial community structures of aerosol samples (Figure 2). The stress value was 0.04, which suggested the NMDS result was representative. As shown in Figure 2, all samples were not distributed by city or by pollution levels (Anosim, p > 0.05), but randomly dispersed and diverged from each other. It meant that bacterial compositions of aerosol samples were different, and other factors not pollution level influenced the bacterial community structure in PM_2.5.

FIGURE 2

Figure 2. Comparing the bacterial communities in unpolluted and polluted samples using NMDS analysis at the OTU level. Red color represented polluted samples and green color referred unpolluted samples.

Bacterial Structural Characteristics of PM_2.5

As many as 19 bacterial phyla and 87 orders were identified in the PM_2.5 samples. During the sampling time, the relative abundance of Proteobacteria in October and November in Beijing was higher than that in September, while Firmicutes was in contrast. No such trend was observed for Shanghai. The relative abundance of Pseudomonadales in Shanghai increased as the sampling time. However, based on the t-test results, no clear trend (p > 0.05) was observed for β-diversities of samples collected from Beijing and Shanghai, thus all samples were put together to analyze the bacterial compositions in polluted and unpolluted days. Among these phyla, Proteobacteria was the most abundant in all aerosol samples (Figure 3A), followed by Firmicutes and Actinobacteria, with total abundances of 92.61–99.87%. The phyla abundance of PM_2.5 samples was not strongly correlated with pollution levels. The average abundance of Proteobacteria (63.30 ± 8.23%) in polluted days was close to that in unpolluted days (64.87 ± 13.42%). In contrast, the bacterial abundances varied greatly among samples. Firmicutes (39.16%) in BJ-1 was about 2.5 times that in BJ-2, although both belonged to unpolluted days (Figure 3A). Similar phenomenon was observed at the order level. The dominant orders of all samples included Bacillales, Pseudomonadales, Enterobacteriales, and Burkholderales, while the proportion of each order was distinct from sample to sample (Figure 3B). Pearson correlation analysis demonstrated that Pseudomonadales were positively correlated with ${\text{NH}}_4^{+}$ (r = 0.73, p < 0.05). Overall, the bacterial structure of each sample had its own characteristics, independent of air pollution levels (PM_2.5 concentration), and environmental factors had obvious influences on airborne bacterial composition.

FIGURE 3

Figure 3. Bacterial community structure for PM_2.5 samples at the phylum (A) and order level (B). Taxonomic summary of the most abundant phylum (>1%) and order (>2%) was indicated in the figure.

Environmental Factors Shaping the Bacterial Community Structure

In order to further confirm the impact of environmental and meteorological conditions on airborne bacterial community, CCA and pCCA were operated at the OTU level. CCA biplot showed that 96.90% of the total variance of airborne bacterial community, calculated from the ratio of constrained inertia and total inertia, was totally explained by the environmental and meteorological factors (Figure 4). Based on the result of Monte Carlo permutation test (p = 0.04), this explanation was significant. Among these, the first two components (CCA1 and CCA2) together explained 53.62% of the total variation of OTU structure (Figure 4). Meanwhile, ${\text{NH}}_4^{+}$, wind speed, and ${\text{SO}}_4^{2-}$ were the predominant factors which shifted the bacterial community structure (r = 0.758, 0.552, and 0.498, respectively; p < 0.05). The results suggested that both ions and meteorological conditions together contributed to the structure of airborne microbial community. The pCCA result further showed that the explanation of chemical ions (60.35%) for the bacterial shift was higher than that of meteorological factors (38.57%).

FIGURE 4

Figure 4. Canonical correspondence analysis (CCA) of bacterial community structures (OTUs) from 10 PM_2.5 samples with respect to the eight environmental and meteorological variables. Arrows indicate the direction and magnitude of measurable variables associated with community structures.

Functional Estimation for PM_2.5 Bacterial Communities

Although the species and abundances of airborne microbes were gradually concerned, their functions were still not recognized, particularly their contribution for aerosol particle formation. PICRUSt software offered potential for gaining insight into the drivers of function dynamics during the formation process of particle pollution. In this study, weighted nearest sequenced taxon index (NSTI) of all samples were lower than 0.1, which suggested that the PICRUSt results were reliable. Overall, 232 reported KEGG pathways were obtained according to the annotated reference genomes for OTU sequences in PM_2.5 samples, and the predominant pathways were mainly related to metabolism, environmental information processing, genetic information processing, cellular processes, and human diseases. This study focused on the pathways of bacterial metabolism and human diseases. Firstly, the proportions of these pathways in polluted days were similar with that in unpolluted days. Secondly, the airborne bacterial community contained a versatile group (Figure 5). For example, Sphingomonas and Ralstonia had a significant correlation with xenobiotics biodegradation and metabolism. Pseudomonas were positively associated with infection disease (Figure 5). About the abundant OTUs with relative abundance >10%, previous studies demonstrated that their related species participated in multiple ecological functions such as N-cycling, forming atmospheric ice nuclei, organic degradation, and opportunistic pathogenic (Table 3).

FIGURE 5

Figure 5. Bacterial taxa are related to KEGG pathways. Pearson’s correlation coefficients were calculated for the relative abundances of genera and KEGG pathways. Red color referred to the positive correlation and green color indicated a negative correlation. Significant correlation was marked with ^∗p < 0.05, ^∗∗p < 0.01.

TABLE 3

Table 3. The dominant bacteria (relative abundance >10%) in PM_2.5 samples correlated with the potential ecological function.

Discussion

Effects of airborne particulate matter on human health have been fully illustrated in terms of chemical composition (Esposito et al., 2012; Pope and Dockery, 2013; Maclntyre et al., 2014). Besides the harm of chemical components, airborne bacteria are ubiquitous (10⁵–10⁶ cells/m³) and also may cause respiratory diseases such as allergy and asthma (Brook et al., 2004; Woo et al., 2013). Although the number of concentration of airborne bacteria has been confirmed to be closely correlated with the air pollution levels (Li et al., 2015; Wei et al., 2016; Du et al., 2018; Lu et al., 2018), an understanding of the structural and functional response of airborne bacteria to air pollution is still elusive. In the present study, PM_2.5 samples were collected to investigate the relationship of the structure and function implications of airborne bacterial community with air pollution.

Compared with unpolluted days, the concentrations of PM_2.5 and chemical species generally increased in polluted days (Zhang et al., 2016). On the basis of this aspect, in this study, six samples were classified into unpolluted samples and four samples were classified into polluted samples (Table 1 and Figure 1). However, despite several OTUs belonging to Acinetobacter and Massilia enriched in polluted days, the airborne bacterial communities showed no significant correlation with air pollution. On the one hand, air pollution posed no impact on all α-diversity indices (Shannon index, Simpson index, Chao1, and Ace), and these parameters showed low correlation with PM_2.5 concentration (r < 0. 2) (Supplementary Table S1). On the other hand, consistent with other studies about airborne bacteria in Beijing (Cao et al., 2014; Lee et al., 2017; Du et al., 2018; Li et al., 2018), the dominant bacteria in PM_2.5 samples belonged to Bacillales, Pseudomonadales, Enterobacteriales, Burkholderiales, and Sphingomonadales, which were similar in whether polluted or unpolluted days (Figure 3). Notably, airborne bacterial structure was not obviously related to air pollution (Figure 2). The proportions of dominant bacteria varied from each other, even from the same city and pollution level (Figure 3). The large differences within the group lowered the differences between groups. It suggested that there were some factors contributing to these differences.

Generally, meteorological conditions such as relative humidity, wind speed, the concentrations of SO₂ and NO₂, were considered to influence the airborne bacterial concentration and structure (Cao et al., 2014; Wei et al., 2017; Xu et al., 2017; Li et al., 2018). High humidity favored the survival of microbes in the air (Cao et al., 2014). SO₂ and NO₂ inhibited the growth and breeding of airborne bacteria (Xu et al., 2017). Wind speed conducted the structure of bacterial community in cloud (Evans et al., 2006; Wei et al., 2017). In this study, high explanation (53.62% by CCA1 and CCA2) was obtained from combining meteorological parameters and water soluble compounds in PM_2.5 (Figure 4). This explanation proportion was higher than the result (13.67%) of Li et al. (2018), while it was comparable with the result (65.9%) of Wei et al. (2017). In accordance with the results in this study, major water soluble ions and wind speed were the crucial environmental factor, different from Woo et al. (2013) who found climate not chemical compounds impacted airborne bacterial community. This might be due to the fact that soluble ions could act as the metabolic substrates. For example, Gao et al. (2016) confirmed the function of ammonia oxidation for the bacteria in PM_2.5. N-cycling-related OTUs were also found in this study (Table 3). Overall, air pollution showed no obvious correlation with bacterial community (Anosim, p > 0.05) as it reflected the combined influence of chemical compounds and meteorological conditions. Nevertheless, ${\text{SO}}_4^{2-}$ and ${\text{NO}}_3^{-}$ different from other compounds played an important role in conducting bacterial community of PM_2.5.

More to the point, part of airborne bacteria was viable (Hara and Zhang, 2012; Zhong et al., 2016; Li et al., 2018), but there was no systematic understanding about the functions of these bacteria. In the previous reports, pathogenic bacteria in air particulate matter were the focus of airborne bacterial function, and some pathogens has been identified in the air particulate matter such as Pseudomonas syringae and Pseudomonas cichorii (Pauwelyn et al., 2011; Vaïtilingom et al., 2013; Du et al., 2018). In the present study, the dominant genus Pseudomonas showed significant correlation with infectious diseases (Figure 5). Acinetobacter modestus (Nemec et al., 2016), an opportunistic pathogenic specie, was abundant in the PM_2.5 samples (Table 3). Other common genera of potential human pathogens including Streptococcus (Lin et al., 2003) and Prevotella (Ibrahim et al., 2017) were also detected in this study (Supplementary Figure S2). Besides the health risk, airborne bacteria may also play an important role in atmospheric biochemistry. The related species of several predominant OTUs, such as Exiguobacterium aurantiacum, Pseudomonas putida, Ralstonia pickettii, and Comamonas testosterone, have been reported to contribute to biodegrading hydrocarbons (Jeswani and Mukherji, 2013), phenolics, aromatic compounds (Hollender et al., 1997; Kahng et al., 2000; Gomes et al., 2005; Table 3). Pearson correlation analysis further confirmed that genera Ralstonia and Sphingomonas were significantly correlated with xenobiotics biodegradation and metabolism (p < 0.05) (Figure 5). Similar in the cloud, some airborne bacteria also act as an efficient ice nuclei because they could induce ice nucleation at a warmer temperature than mineral dusts (Amato et al., 2015). For example, the abundance related specie in this study, P. putida, was an ice nucleation active specie and showed the ice nucleation activity at −3°C (Hill et al., 2014). All these functions implied that airborne bacteria might be involved in changing meteorology and climate.

Conclusion

Airborne bacteria, including their structure response to different pollution levels, correlation with environmental factors, and potential role in atmospheric progresses, remain elusive. In the present study, airborne bacterial structure showed no significant correlation between polluted and unpolluted days, and between two cities in fall. The relationship of environmental factors and meteorological conditions with airborne bacterial community was explored, and found that ${\text{NH}}_4^{+}$, ${\text{SO}}_4^{2-}$, and wind speed significantly structured the bacterial community. Besides health risk from several potential pathogenic genera (Pseudomonas, Streptococcus, and Prevotella), some abundant species related to organic degradation and nitrification/denitrification may contribute to atmospheric biochemistry process. This knowledge helps for the comprehensive understanding of bacterial community biodiversity and their potential role in atmospheric progress. However, in the future study, more aerosol samples obtained through long term such as a whole year and high resolution sampling should be conducted to further disclose the correlation between airborne bacterial diversity and air quality. In addition, studies depending on the RNA technique and biochemical culture experiments need to be performed to confirm the functions in the atmosphere.

Data Availability Statement

Raw 16S rRNA gene sequences reported here were deposited into the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/sra) database (Accession Number: SRP167452).

Author Contributions

YP designed the study, performed the experiments, analyzed the data, and wrote the manuscript. XP performed the experiments and analyzed the data. HoX designed the experiments. HuX revised the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Nos. 31900106, 41425014, and 41867063) and the Doctoral Scientific Research Foundation of East China University of Technology (No. DHBK2017143).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02369/full#supplementary-material

Footnotes

1. ^ http://www.ncbi.nlm.nih.gov/sra

References

Amato, P., Joly, M., Schaupp, C., Attard, E., Möhler, O., Morris, C. E., et al. (2015). Survival and ice nucleation activity of bacteria as aerosols in a cloud simulation chamber. Atmos. Chem. Phys. 15, 6455–6465. doi: 10.5194/acp-15-6455-2015

CrossRef Full Text | Google Scholar

Bauer, H., Giebl, H., Hitzenberger, R., Kasper-Giebl, A., Reischl, G., Zibuschka, F., et al. (2003). Airborne bacteria as cloud condensation nuclei. J. Geophys. Res. 108:4658. doi: 10.1029/2003JD003545

CrossRef Full Text | Google Scholar

Bowers, R. M., McLetchie, S., Knight, R., and Fierer, N. (2011). Spatial variability in airborne bacterial communities across land-use types and their relationship to the bacterial communities of potential source environments. ISME J. 5, 601–612. doi: 10.1038/ismej.2010.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Brook, R. D., Franklin, B., Cascio, W., Hong, Y., Howard, G., Lipsett, M., et al. (2004). Air pollution and cardiovascular disease: a statement for healthcare professionals from the expert panel on population and prevention science of the American Heart Association. Circulation 109, 2655–2671. doi: 10.1161/01.CIR.0000128587.30041.C8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, T., Chen, L., Wu, S., Qian, L., and Ren, Q. (2010). Selection of denitrifying phosphorus-removing bacteria and its characteristic. Environ. Sci. 31, 2487–2492.

PubMed Abstract | Google Scholar

Cao, C., Jiang, W., Wang, B., Fang, J., Lang, J., Tian, G., et al. (2014). Inhalable microorganisms in Beijing’s PM2.5 and PM10 pollutants during a severe smog event. Environ. Sci. Technol. 48, 1499–1507. doi: 10.1021/es4048472

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Li, X., Zhu, T., Han, Y., and Lv, D. (2017). PM2.5-bound PAHs in three indoor and one outdoor air in Beijing: concentration, source and health risk assessment. Sci. Total Environ. 586, 255–264. doi: 10.1016/j.scitotenv.2017.01.214

PubMed Abstract | CrossRef Full Text | Google Scholar

Delort, A. M., Vaïtilingom, M., Amato, P., Sancelme, M., Parazols, M., Mailhot, G., et al. (2010). A short overview of the microbial population in clouds: potential roles in atmospheric chemistry and nucleation processes. Atmos. Res. 98, 249–260. doi: 10.1016/j.atmosres.2010.07.004

CrossRef Full Text | Google Scholar

Du, P., Du, R., Ren, W., Lu, Z., and Fu, P. (2018). Seasonal variation characteristic of inhalable microbial communities in PM2.5 in Beijing city. China. Sci. Total Environ. 61, 308–315. doi: 10.1016/j.scitotenv.2017.07.097

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/NMETH.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Esposito, V., Lucariello, A., Savarese, L., Cinelli, M. P., Ferraraccio, F., Bianco, A., et al. (2012). Morphology changes in human lung epithelial cells after exposure to diesel exhaust micron sub particels (PM1.0) and pollen allergens. Environ. Pollut. 171, 162–167. doi: 10.1016/j.envpol.2012.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, C. A., Coombes, P. J., and Dunstan, R. H. (2006). Wind, rain and bacteria: the effect of weather on the microbial composition of 16343588 roof-harvested rainwater. Water. Res. 40, 37–44. doi: 10.1016/j.watres.2005.10.034

CrossRef Full Text | Google Scholar

Gao, J., Fan, X., Pan, K., Li, H., and Sun, L. (2016). Diversity, abundance and activity of ammonia-oxidizing microorganisms in fine particulate matter. Sci. Rep. 6:38785. doi: 10.1038/srep38785

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, M., Qiu, T., Jia, R., Han, M., Song, Y., and Wang, X. (2015). Concentration and size distribution of viable bioaerosols during non-haze and haze days in Beijing. Environ. Sci. Pollut. Res. Int. 22, 4359–4368. doi: 10.1007/s11356-014-3675-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes, N. C., Kosheleva, I. A., Abraham, W. R., and Smalla, K. (2005). Effects of the inoculant strain Pseudomonas putida KT2442 (pNF142) and of naphthalene contamination on the soil bacterial community. FEMS Microbiol. Ecol. 54, 21–33. doi: 10.1016/j.femsec.2005.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, S., Hu, M., Zamora, M. L., Peng, J., Shang, D., Zheng, J., et al. (2014). Elucidating severe urban haze formation in China. Proc. Natl. Acad. Sci. U.S.A. 111, 17373–17378. doi: 10.1073/pnas.1419604111

PubMed Abstract | CrossRef Full Text | Google Scholar

Hara, K., Maki, T., Kobayashi, F., Kakikawa, M., Wada, M., and Matsuki, A. (2016). Variations of ice nuclei concentration induced by rain and snowfall within a local forested site in Japan. Atmos. Environ. 127, 1–5. doi: 10.1016/j.atmosenv.2015.12.009

CrossRef Full Text | Google Scholar

Hara, K., and Zhang, D. (2012). Bacterial abundance and viability in long-range transported dust. Atmos. Environ. 47, 20–25. doi: 10.1016/j.atmosenv.2011.11.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, T., Moffett, B. F., Demott, P. J., Georgakopoulos, D. G., Stump, W. L., and Franc, G. D. (2014). Measurement of ice nucleation-active bacteria on plants and in precipitation by quantitative PCR. Appl. Environ. Microbiol. 80, 1256–1267. doi: 10.1128/AEM.02967-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollender, J., Hopp, J., and Dott, W. (1997). Degradation of 4-chlorophenol via meta cleavage pathway by Comamonas testosteroni JH5. Appl. Environ. Microbiol. 63, 4567–4572. doi: 10.0000/PMID16535738

PubMed Abstract | CrossRef Full Text | Google Scholar

Ibrahim, M., Subramanian, A., and Anishetty, S. (2017). Comparative pan genome analysis of oral Prevotella species implicated in periodontitis. Funct. Integr. Genomics 17, 513–536. doi: 10.1007/s10142-017-0550-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeswani, H., and Mukherji, S. (2013). Batch studies with Exiguobacterium aurantiacum degrading structurally diverse organic compounds and its potential for treatment of biomass gasification wastewater. Int. Biodeter. Biodegr. 80, 1–9. doi: 10.1016/j.ibiod.2013.02.002

CrossRef Full Text | Google Scholar

Kahng, H., Byrne, A. M., Olsen, R. H., and Kukor, J. J. (2000). Characterization and role of tbuX in utilization of toluene by Ralstonia pickettii PKO1. J. Bacteriol. 182, 1232–1242. doi: 10.1128/JB.182.5.1232-1242.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Langille, M. G. I., Zaneveld, J., Caporaso, J. G., Mcdonald, D., Dan, K., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. Y., Park, E. H., Lee, S., Ko, G., Honda, Y., Hashizume, M., et al. (2017). Airborne bacterial communities in three east asian cities of China, South Korea, and Japan. Sci. Rep. 7:5545. doi: 10.1038/s41598-017-05862-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. R., Lee, R. E., Strong-Gunderson, J. M., and Minges, S. R. (1995). Isolation of ice-nucleating active bacteria from the freeze-tolerant frog Rana sylvatica. Cryobiology 32, 358–365. doi: 10.1116/1.2800327

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Pan, X., Zhang, D., Lee, D., Ai-Misned, F. A., and Mortuza, M. G. (2015). Anaerobic nitrate reduction with oxidation of Fe(II) by Citrobacter Freundii strain PXL1 – a potential candidate for simultaneous removal of As and nitrate from groundwater. Ecol. Eng. 77, 196–201. doi: 10.1016/j.atmosenv.2015.09.070

CrossRef Full Text | Google Scholar

Li, W., Yang, J., Zhang, D., Li, B., Wang, E., and Yuan, H. (2018). Concentration and community of airborne bacteria in response to cyclical haze events during the fall and midwinter in Beijing. China. Front. Microbiol. 9:1741. doi: 10.3389/fmicb.2018.01741

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, F. Y., Weisman, L. E., Troendle, J., and Adams, K. (2003). Prematurity is the major risk factor for late-onset group B streptococcus disease. J. Infect. Dis. 188, 267–271. doi: 10.1086/376457

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Yang, J., Yu, X., Chen, G., and Yu, Z. (2013). Patterns in the composition of microbial communities from a subtropical river: effects of environmental, spatial and temporal factors. PLoS One 8:e81232. doi: 10.1371/journal.pone.0081232

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, R., Li, Y., Li, W., Xie, Z., Fan, C., Liu, P., et al. (2018). Bacterial community structure in atmospheric particulate matters of different sizes during the haze days in Xi’an, China. Sci. Total Environ. 637-638, 244–252. doi: 10.1016/j.scitotenv.2018.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Maclntyre, E. A., Gehring, U., Mölter, A., Fuertes, E., Klümper, C., Krämer, U., et al. (2014). Air pollution and respiratory infections during early childhood: an analysis of 10 european birth cohorts within the ESCAPE project. Environ. Health Persp. 122, 107–113. doi: 10.1289/ehp.1306755

PubMed Abstract | CrossRef Full Text | Google Scholar

Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Maki, T., Furumoto, S., Asahi, Y., Lee, K. C., Watanabe, K., Aoki, K., et al. (2018). Long-range transported bioaerosols captured in snow cover on Mount Tateyama, Japan: impacts of Asian-dust events on airborne bacterial dynamics relating to ice-nucleation activities. Atmos. Chem. Phys. 18, 8155–8171. doi: 10.5194/acp-18-8155-2018

CrossRef Full Text | Google Scholar

Nemec, A., Radolfova-Krizova, L., Maixnerova, M., Vrestiakova, E., Jezek, P., and Sedo, O. (2016). Taxonomy of haemolytic and/or proteolytic strains of the genus Acinetobacter with the proposal of Acinetobacter courvalinii sp. nov. (genomic species 14 sensu Bouvet & Jeanjean), Acinetobacter dispersus sp. nov. (genomic species 17), Acinetobacter modestus sp. nov., Acinetobacter proteolyticus sp. nov. and Acinetobacter vivianii sp. nov. Int. J. Syst. Evol. Microbiol. 66, 1673–1685. doi: 10.1099/ijsem.0.000932

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, B., Agrawal, M., and Singh, S. (2014). Coal mining activities change plant community structure due to air pollution and soil degradation. Ecotoxicology 23, 1474–1483. doi: 10.1007/s10646-014-1289-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Pauwelyn, E., Vanhouteghem, K., Cottyn, B., Vos, P. D., Maes, M., Bleyaert, P., et al. (2011). Epidemiology of pseudomonas cichorii, the cause of lettuce midrib rot. J. Phytopathol. 159, 298–305. doi: 10.1111/j.1439-0434.2010.01764.x

CrossRef Full Text | Google Scholar

Pope, C. A., and Dockery, D. W. (2013). Air pollution and life expectancy in China and beyond. Proc. Natl. Acad. Sci. U.S.A. 110, 12861–12862. doi: 10.1073/pnas.1310925110

PubMed Abstract | CrossRef Full Text | Google Scholar

Raisi, L., Aleksandropoulou, V., Lazaridis, M., and Katsivela, E. (2013). Size distribution of viable, cultivable, airborne microbes and their relationship to particulate matter concentrations and meteorological conditions in a Mediterranean site. Aerobiologia 29, 233–248. doi: 10.1007/s10453-012-9276-9

CrossRef Full Text | Google Scholar

Ramette, A., Frapolli, M., Saux, M. F. L., Gruffaz, C., Meyer, J. M., Défago, G., et al. (2011). Pseudomonas protegens sp. nov., widespread plant-protecting bacteria producing the biocontrol compounds 2,4-diacetylphloroglucinol and pyoluteorin. Syst. Appl. Microbiol. 34, 180–188. doi: 10.1016/j.syapm.2010.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, Y., Yang, L., and Liang, X. (2014). The characteristics of a novel heterotrophic nitrifying and aerobic denitrifying bacterium, Acinetobacter junii YB. Bioresour. Technol. 171, 1–9. doi: 10.1016/j.biortech.2014.08.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Z., Mu, Y., Liu, Y., and Shao, L. (2013). A comparison study on airborne particles during haze days and non-haze days in Beijing. Sci. Total Environ. 456-457, 1–8. doi: 10.1016/j.scitotenv.2013.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaïtilingom, M., Amato, P., Sancelme, M., Laj, P., Leriche, M., and Delort, A. M. (2010). Contribution of microbial activity to carbon chemistry in clouds. Appl. Environ. Microbiol. 76, 23–29. doi: 10.1128/AEM.01127-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaïtilingom, M., Deguillaume, L., Vinatier, V., Sancelme, M., Amato, P., Chaumerliac, N., et al. (2013). Potential impact of microbial activity on the oxidant capacity and organic carbon budget in clouds. Proc. Natl. Acad. Sci. U.S.A. 110, 559–564. doi: 10.1073/pnas.1205743110

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, K., Zheng, Y., Li, J., Shen, F., Zou, Z., Fan, H., et al. (2015). Microbial aerosol characteristics in highly polluted and nearpristine environments featuring different climatic conditions. Sci. Bull. 60, 1439–1447. doi: 10.1007/s11434-015-0868-y

CrossRef Full Text | Google Scholar

Wei, K., Zou, Z., Zheng, Y., Li, J., Shen, F., Wu, C., et al. (2016). Ambient bioaerosol particle dynamics observed during haze and sunny days in Beijing. Sci. Total Environ. 550, 751–759. doi: 10.1016/j.scitotenv.2016.01.137

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, M., Xu, C., Chen, J., Zhu, C., Li, J., and Lv, G. (2017). Characteristics of bacterial community in cloud water at Mt Tai: similarity and disparity under polluted and non-polluted cloud episodes. Atmos. Chem. Phys. 17, 5253–5270. doi: 10.5194/acp-17-5253-2017

CrossRef Full Text | Google Scholar

Woo, A. D., Brar, M. S., Chan, Y., Lau, M. C. Y., Leung, F. C. C., Scott, J. A., et al. (2013). Temporal variation in airborne microbial populations and microbially-derived allergens in a tropical urban landscape. Atmos. Environ. 74, 291–300. doi: 10.1016/j.atmosenv.2013.03.047

CrossRef Full Text | Google Scholar

Xu, C., Wei, M., Chen, J., Wang, X., Zhu, C., Li, J., et al. (2017). Bacterial characterization in ambient submicron particles during severe haze episodes at Ji’nan, China. Sci. Total Environ. 580, 188–196. doi: 10.1016/j.scitotenv.2016.11.145

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Huang, W., Cai, T., Fang, D., Wang, Y., Song, J., et al. (2016). Concentrations and chemical compositions of fine particles (PM2.5) during haze and non-haze days in Beijing. Atmos. Res. 174, 62–69. doi: 10.1016/j.atmosres.2016.02.003

CrossRef Full Text | Google Scholar

Zhong, X., Qi, J., Li, H., Dong, L., and Gao, D. (2016). Seasonal distribution of microbial activity in bioaerosols in the outdoor environment of the Qingdao coastal region. Atmos. Environ. 140, 506–513. doi: 10.1016/j.atmosenv.2016.06.034

CrossRef Full Text | Google Scholar