Original Research ARTICLE
Diversity and Composition of Airborne Fungal Community Associated with Particulate Matters in Beijing during Haze and Non-haze Days
- China Pharmaceutical Culture Collection, Institute of Medicinal Biotechnology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
To assess the diversity and composition of airborne fungi associated with particulate matters (PMs) in Beijing, China, a total of 81 PM samples were collected, which were derived from PM2.5, PM10 fractions, and total suspended particles during haze and non-haze days. The airborne fungal community in these samples was analyzed using the Illumina Miseq platform with fungi-specific primers targeting the internal transcribed spacer 1 region of the large subunit rRNA gene. A total of 797,040 reads belonging to 1633 operational taxonomic units were observed. Of these, 1102 belonged to Ascomycota, 502 to Basidiomycota, 24 to Zygomycota, and 5 to Chytridiomycota. The dominant orders were Pleosporales (29.39%), Capnodiales (27.96%), Eurotiales (10.64%), and Hypocreales (9.01%). The dominant genera were Cladosporium, Alternaria, Fusarium, Penicillium, Sporisorium, and Aspergilus. Analysis of similarities revealed that both particulate matter sizes (R = 0.175, p = 0.001) and air quality levels (R = 0.076, p = 0.006) significantly affected the airborne fungal community composition. The relative abundance of many fungal genera was found to significantly differ among various PM types and air quality levels. Alternaria and Epicoccum were more abundant in total suspended particles samples, Aspergillus in heavy-haze days and PM2.5 samples, and Malassezia in PM2.5 samples and heavy-haze days. Canonical correspondence analysis and permutation tests showed that temperature (p < 0.01), NO2 (p < 0.01), PM10 (p < 0.01), SO2(p < 0.01), CO (p < 0.01), and relative humidity (p < 0.05) were significant factors that determine airborne fungal community composition. The results suggest that diverse airborne fungal communities are associated with particulate matters and may provide reliable data for studying the responses of human body to the increasing level of air pollution in Beijing.
In recent years, particulate matters (PMs) have become a major air pollutant in densely populated cities or urban areas, and their impact to the public health may be profound (Zhang et al., 2010; Cheng et al., 2013; Evans et al., 2013; Beelen et al., 2014; DeFranco et al., 2015). PM pollutants are usually highly toxic to humans and can cause severe sickness or even death. Different PMs may affect various parts of the human body. For example, total suspended particles (TSP, <100 μm) cause skin allergy, PM10 (<10 μm) deposit mainly in the head airways, and PM2.5 (<2.5 μm) are more likely to penetrate and deposit deeper in the tracheobronchial and alveolar regions (Brook et al., 2004). Historical data suggest that exposure to high level of atmospheric PM pollutants, particularly PM2.5, increases morbidity and mortality, and thereby shortening the life span (Hunt et al., 2003; Pope et al., 2009).
While the physical and chemical properties of PM pollutants have been extensively studied (Putaud et al., 2004; Kim et al., 2005; Kuhn et al., 2005; Bressi et al., 2013; Li et al., 2013; Huang et al., 2014; Yang et al., 2016), relatively little is known about microbes associated with PMs, especially airborne fungi. Fungal spores can account for large proportions of air PMs, and some fungi are major pathogens or allergens for humans, animals, and plants. Moreover, air is the primary medium for fungal dispersal (Adhikari et al., 2004). Hitherto, a few studies carried out provide only limited insight into airborne fungi associated with PMs. Airborne fungi associated with PMs have been reported from different regions of the world, including Beijing (China; Cao et al., 2014; Gao et al., 2015), Xi'an (China; Li et al., 2015), Hong Kong (China; Woo et al., 2013), Denver and Greeley (USA; Bowers et al., 2013), New Haven (USA; Yamamoto et al., 2012), Berkeley (USA; Adams et al., 2013), Rehovot (Israel; Dannemiller et al., 2014), Graz (Austria; Haas et al., 2013), Manaus (Brazil; Womack et al., 2015), and Mainz (Germany; Fröhlich-Nowoisky et al., 2009).
Most previous studies focused on the fungal diversity associated with PMs using traditional isolation methods (Awad et al., 2013; Haas et al., 2013; Almaguer et al., 2014; Gao et al., 2015; Li et al., 2015). However, the cultured method cannot reflect the actual diversity of airborne fungi because of its selectivity. In addition, a few studies surveyed PMs-associated airborne fungi using conventional DNA-based molecular methods (e.g., cloning approaches, q-PCR, DGGE; Fröhlich-Nowoisky et al., 2009; Li et al., 2010; Rittenour et al., 2014). In recent years, diversity and community composition of PMs-associated airborne fungi have been revealed by high-throughput sequencing (Yamamoto et al., 2012, 2014; Woo et al., 2013; Cao et al., 2014; Dannemiller et al., 2014; Emerson et al., 2015; Womack et al., 2015), which can significantly enhance the characterization of fungal diversity compared to traditional methods.
Beijing, the capital of China with a population of over 20 million, in recent years, has been suffering from frequent smog events caused mainly by PMs (Zhang et al., 2010, 2012; Cheng et al., 2013; Cao et al., 2014; Gao et al., 2015). While previous findings have improved our understanding of air fungal diversity in urban area of Beijing (Cao et al., 2014; Gao et al., 2015), the fungal community associated with PMs during haze and non-haze days are not well-understood. The aim of this study was to use high-throughput sequencing to investigate the PMs-associated air fungal community in Beijing to address the following questions: (1) what are air fungal diversity and community composition associated with PMs, (2) does the fungal community composition differ among various PMs (PM2.5, PM10, and TSP), and (3) does the fungal community composition vary among different air quality levels?
Materials and Methods
Samples for PM2.5, PM10, and TSP were collected from the roof top of the Conference Building at Institute of Medicinal Biotechnology, Chinese Academy of Medical Sciences (39°52′43″N, 116°23′21″E, ~8 m above the ground, ~400 m from Temple of Heaven Park), an area without major pollution sources nearby. Sampling was conducted by three portable ambient air samplers (Air Metrics, USA), one of which, impactors were removed from the filter holder for TSP samples, the other one was assembled with PM10 impactor for PM10 samples, and another one was assembled with both PM10 and PM2.5 impactors for PM2.5 samples. Ambient air was drawn at an average flow rate of 5 L/min for 24 h (12:00 p.m. to 12:00 p.m. at next day) per sampling day. PM samples were collected on 47-mm quartz aerosol collection filters (Pall, USA). Prior to sampling, all the filters were sterilized by autoclaving at 121°C for 20 min. The filter holders were cleaned with 75% ethanol and all the tools used for changing new filters were autoclaved every day to avoid contamination. After sampling, the filters were kept in a 2-ml centrifuge tube and stored at −20°C.
A total of 81 PM samples were collected for 27 days (June–November 2014) during different air quality levels (Table 1). The air quality index (AQI), which is an index for reporting daily air quality (Wang et al., 2006), was used to indicate the pollutant level. We define a day with AQI lower than 100 as a non-haze day, that with AQI in the range of 100–200 as a light-haze day, and that with AQI higher than 200 as a heavy-haze day. The environmental variables such as PM2.5, PM10, CO, SO2, and NO2 were recorded from the monitoring data of Temple of Heaven Park site (~800 m from sampling site) of Beijing Municipal Environmental Monitoring Center (http://zx.bjmemc.com.cn/). Temperature (Temp) and relative humidity (RH) were recorded according to the reports of Chinese National Meteorological Center (http://www.nmc.cn/). The air environmental variables are listed in Table S1.
Table 1. Summary data for sequencing data from the 81 particulate matter samples in the present study.
DNA Extraction and PCR Amplification
To obtain high total DNA yield, a half filter was cut into small pieces by scissor. Genomic DNA was extracted using PowerSoil DNA isolation kit (MoBio, USA) according to manufacturer's protocols. The internal transcribed spacer 1 (ITS1) region of the fungal ribosomal RNA gene was amplified by PCR (95°C for 3 min, followed by 35 cycles at 95°C for 30 s, 55°C for 40 s, and 72°C for 45 s, and a final extension at 72°C for 10 min) using primers ITS5 (5′-barcode-GGAAGTAAAAGTCGTAACAAGG-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′; Bellemain et al., 2010), where barcode is an 8-base sequence unique to each sample. PCR amplification was conducted using high fidelity TransStart Fastpfu DNA Polymerase (Transgen, China). PCR reactions were performed in triplicate in 20 μL mixtures containing 4 μL of 5 × FastPfu Buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.4 μL of FastPfu Polymerase, and 10 ng of template DNA.
Illumina MiSeq Sequencing
Amplicons were purified using AxyPrep DNA Gel Extraction Kit (Axygen, USA) according to the manufacturer's instructions and quantified using QuantiFluor™-ST (Promega, USA). Purified amplicons were pooled in equimolar and paired-end sequenced (2 × 250 bp) on an Illumina MiSeq platform according to the standard protocols. The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database under accession number SRP068265.
Processing of Sequencing Data
Raw files were demultiplexed, quality-filtered using QIIME (version 1.18; Caporaso et al., 2010) with the following criteria: (i) The 300 bp reads were truncated at any site receiving an average quality score <20 over a 50-bp sliding window, discarding the truncated reads that were shorter than 50 bp; (ii) exact barcode matching, two nucleotide mismatch in primer matching, reads containing ambiguous characters were removed; (iii) only sequences that overlap longer than 10 bp were assembled according to their overlap sequence. Reads which could not be assembled were discarded.
All samples were normalized at the same sequence depth (9840 reads). Operational taxonomic units (OTUs) were clustered with 97% similarity cutoff using UPARSE (version 7.1, http://drive5.com/uparse/); chimeric sequences were identified and removed using UCHIME. The OTUs were used as a basis for calculating alpha diversity and beta diversity metrics. The taxonomy of each ITS1 sequence was analyzed by RDP Classifier (Wang et al., 2007) against the Unite database (Release 7.0, http://unite.ut.ee/index.php; Kõljalg et al., 2013) using confidence threshold of 70%. Statistical analyses of the OTU diversity of each air sample via Chao1, Good's coverage estimator, and Shannon's index (H′) were performed using QIIME 1.8.0 software (Caporaso et al., 2010). Statistical comparisons of Chao1 and Shannon indices among different samples were made by one-way analysis of variance (ANOVA) using Statistical package for the social sciences (SPSS, version 23.0). An analysis of similarities (ANOSIM) was performed using QIIME 1.8.0 software (Caporaso et al., 2010) to determine whether various types of air samples had significantly different fungal communities. Metastats analysis (White et al., 2009) was used to identify the significantly different fungal genera in different types of air samples using Mothur software (Version 1.35.1). The program Metastats can use count data from annotated sequences to compare two populations in order to detect differentially abundant features and then produce a tab-delimited table displaying the mean relative abundance of a feature, variance, and standard error together with a p-value and q-value to describe significance of the detected variations (http://metastats.cbcb.umd.edu/).
Airborne Fungal Richness and Diversity
A total of 797,040 reads with 9840 reads per sample passed quality filtering, subsampled, and were clustered into 1633 unique OTUs (97% sequence similarity). The alpha diversity in each sample was estimated by Chao1, Good's coverage estimator, and Shannon's indices (Table 1). The rarefaction curves and Shannon index curves indicated that high-throughput sequencing captured the dominant phylotypes (Figure S1). On average, there were 94 (from 34 to 216), 155 (from 59 to 285), and 153 (from 80 to 244) OTUs in samples of PM2.5, PM10, and TSP, respectively.
ANOVA showed that TSP and PM10 samples had higher Chao1-values as compared with PM2.5 samples (p < 0.05), which indicated higher richness in TSP and PM10 samples. A relatively lower Shannon index in PM10 and TSP samples was also observed, but there was no significant difference (p > 0.05). Samples in heavy-haze days showed a higher value of Shannon index, which indicated a higher diversity than non-haze and light-haze days (p < 0.05), whereas there was no significant difference between non-haze and light-haze days (Figure 1).
Figure 1. Statistical comparisons of Chao1 and Shannon indices among three PM types (PM2.5, PM10, and TSP) and three air quality levels (non-haze, light-haze, and heavy-haze). (A) Values of Chao1 in three particulate matter samples; (B) Values of Shannon index in three particulate matter samples; (C) Values of Chao1 in three air quality levels; (D) Values of Shannon index in three air quality levels.
Taxonomic Composition of Airborne Fungal Community
A total of four phyla, 71 orders, and 368 genera were identified in the present study. The airborne fungal community was dominated by the Ascomycota, followed by the Basidiomycota, Zygomycota, and Chytridiomycota (Figure 2). The dominant orders were Pleosporales (29.39%), Capnodiales (27.96%), Eurotiales (10.64%), and Hypocreales (9.01%; Figure S2), which belonged to Ascomycota. The abundant orders in Basidiomycota included the Ustilaginales (4.99%), Polyporales (3.57%), and Agaricales (2.78%; Figure S2). Two genera Cladosporium and Alternaria were the most abundant in fungal community, followed by Fusarium, Penicillium, Sporisorium, and Aspergilus (Figure 3).
Cladosporium and Alternaria were two dominant genera in all types of PM samples. Unlike in PM2.5 and PM10 samples, Alternaria was more abundant than Cladosporium in TSP samples. Fusarium was the third abundant genus in PM2.5 and PM10 samples and the fourth abundant genus in TSP samples. Epicoccum was one of the dominant genera in TSP samples, whereas there was only a few Epicoccum in PM2.5 and PM10 samples. Sporisorium showed more abundant in PM2.5 samples. The top three genera in PM2.5 and PM10 samples were same (i.e., Cladosporium, Alternaria, and Fusarium; Table S2).
In addition, Cladosporium and Alternaria were two dominant genera in non-haze and light-haze days, whereas, Alternaria was most abundant in heavy-haze days. Sporisorium showed more abundant in none-haze and light-haze days but not (<1%) in heavy-haze days. Trametes was the fourth abundant in non-haze days but not in both light-haze and heavy-haze days. Fifteen dominant genera accounted for 82% of total reads and more diverse genera were observed in heavy-haze days (Table S3).
Airborne Fungal Communities Associated with Different PMs and Air Quality Levels
ANOSIM (R = 0.175, p = 0.001) indicated that there was a significant effect of particle size on airborne fungal community composition and a Venn diagram showed that 462 OTUs were shared among PM2.5, PM10, and TSP (Figure 4A). Moreover, there was a significant difference in fungal communities among non-haze, light-haze, and heavy-haze samples (R = 0.076, p = 0.006). In addition, the Venn diagram showed that there were 433 OTUs shared among non-haze, light-haze, and heavy-haze days (Figure 4B).
Figure 4. Venn diagrams illustrating the number of unique and shared OTUs among (A) three PM types (PM2.5, PM10, and TSP) and (B) three air quality levels (non-haze, light-haze, and heavy-haze).
A heatmap diagram illustrated the distributions of abundant fungal genera in PMs collected from various air quality days (Figure 5). For examples, Aspergillus increased as the haze got heavy and decreased as the particle got bigger; Penicillium and Malassezia were more abundant in PM2.5 samples and increased during heavy-haze days, especially in PM2.5 and PM10 samples; Alternaria was more abundant in TSP samples when compared with PM10 and PM2.5 samples; Cladosporium was less during heavy-haze days; Fusarium was more abundant in PM2.5 and TSP samples during non-haze days, whereas it is most abundant in PM10 samples during heavy-haze days; Epicoccum showed significant abundance in TSP samples.
Figure 5. A heatmap diagram showing the distribution of fungal genera in three PM types (PM2.5, PM10, and TSP) and three air quality levels (non-haze, light-haze, and heavy-haze).
Metastats analysis demonstrated that some fungal genera showed significant differences among samples of PM2.5, PM10, and TSP. For example, many fungal genera, including Alternaria, Penicillium, Aspergillus, Epicoccum, Malassezia, Nigrospora, Talaromyces, Funalia, Emericella, Neurospora, and Oedocephalum, showed significant differences between samples of PM2.5 and TSP (Table 2). In addition, some fungal genera showed significant differences among non-haze, light-haze, and heavy-haze days. For example, Cladosporium, Aspergillus, Trametes, Davidiella, Schizophyllum, Nigrospora, Talaromyces, Flammulina, Emericella, Ustilago, Irpex, Phaeoacremonium, Periconia, Chaetomium, Fimetariella, Torula, Mycosphaerella, Bjerkandera, Tilletiopsis, Trichoderma, and Hyphodiscus showed significant differences between non-haze and heavy-haze days (Table 3).
Table 2. Metastats analysis showing the fungal genera which are significantly different among three PM types.
Table 3. Metastats analysis showing the fungal genera which are significantly different among three air quality levels.
The Effects of Environmental Factors on the Airborne Fungal Community
A Pearson correlation analysis was used to examine the relationships between the environmental factors and fungal diversity (Table S4). Values of Shannon index were significantly correlated with RH (r = 0.350, p < 0.01), PM2.5 (r = 0.339, p < 0.01), CO (r = 0.331, p < 0.01), NO2 (r = 0.264, p < 0.05), PM10 (r = 0.239, p < 0.05), and SO2 (r = 0.224, p < 0.05).
Canonical correspondence analysis (CCA, Figure 6) and permutation tests (Table S5) were performed to examine the relationships between the environmental factors and fungal community composition. Six environmental factors, including temperature (r2 = 0.2792, p < 0.01), NO2 (r2 = 0.2173, p < 0.01), PM10 (r2 = 0.1692, p < 0.01), SO2 (r2 = 0.1382, p < 0.01), CO (r2 = 0.1522, p < 0.01), and RH (r2 = 0.0977, p < 0.05), showed significant correlations with airborne fungal community composition.
Figure 6. Canonical correspondence analysis showing the relationships between environmental factors and airborne fungal community composition.
The aim of the present study was to analyze the airborne fungal community associated with PMs in Beijing using high-throughput sequencing. This study demonstrated the diversity and composition of airborne fungal community from different PM types and air quality levels, which was not clarified before in Beijing.
In the present study, 1633 OTUs (34–285 OTUs per sample) were obtained, which indicated the high richness of airborne fungi associated with PMs. These results are somewhat consistent with a previous study that the fungal OTU numbers associated with the PMs in Rehovot (Israel) were 121–178 OTUs (Dannemiller et al., 2014). Furthermore, the number of OTUs associated with PM10 and TSP was higher than those that were associated with PM2.5 samples. Raisi et al. (2013) found that the highest concentrations of the airborne fungi were determined at aerodynamic diameters between 2.1 and 3.3 μm. In addition, during heavy-haze days, the values of Shannon index were significantly higher than those in non-haze and light-haze days, which are consistent with the concentration of airborne fungi in a culture-dependent study (Li et al., 2015). A high number of particulate matters were observed in haze days as compared with non-haze and light-haze days, which would provide a number of favorable substrates for fungal colonization. Therefore, the fungal diversity in haze days was increased.
The phyla Ascomycota and Basidiomycota dominated the PMs-associated fungal community in Beijing. This result is similar to the previous studies of fungi in the atmosphere (Bowers et al., 2013; Dannemiller et al., 2014; Oh et al., 2014; Womack et al., 2015). Moreover, this result makes sense in light of feature of the Ascomycota, which have single-celled or filamentous vegetative growth forms that are small to become aerosolized, whereas most of the Basidiomycota are too large to be easily aerosolized (Womack et al., 2015). However, some other studies showed that the abundance of Basidiomycota was more than Ascomycota (Barnes et al., 2001; Yamamoto et al., 2014).
In the phylum Ascomycota, Capnodiales and Pleosporales were two abundant orders. Capnodiales have been previously detected in air samples (Oh et al., 2014; Rittenour et al., 2014; Emerson et al., 2015). Pleosporales also have been reported in Kansas City as a broadly cross-reactive and important allergenic fungal order to local residents (Barnes et al., 2001; Sáenz-de-Santamaría et al., 2006; Rittenour et al., 2014). Alternaria and Cladosporium were the most abundant fungi in our study which is consistent with a previous study (Yamamoto et al., 2012).
Different PMs may affect various regions of the human body: TSP cause skin allergy, PM10 deposit mainly in the head airways, and PM2.5 are more likely to penetrate and deposit deeper in the tracheobronchial and alveolar regions (Brook et al., 2004). Therefore, it is possible to estimate certain fungal infections based on the PM type and haze severity. In the present study, most of the pathogenic and allergenic fungi increased with the air quality levels. For example, the relative abundance of Aspergillus significantly increased during heavy-haze days, making human body more susceptible to this pathogenic and allergenic genus in haze days. In addition, some fungal genera which were dominant genera (>1%) and showed significant differences between heavy-haze and non-haze samples may serve as indicators for assessment to the severity of the haze days with other environmental variables. For example, Davidiella and Trametes showed remarkable differences between heavy-haze and non-haze samples, especially in PM10 and TSP samples. The fungal genus Schizophyllum also dramatically changed between heavy-haze and non-haze samples in PM2.5 samples. Previous studies also demonstrated that these three genera were prevalent in the atmosphere (Yamamoto et al., 2012, 2014; Oh et al., 2014; Rittenour et al., 2014).
Fungal DNA detected in airborne filter samples may be from spores or fragments of fungal hyphae. The genera Alternaria and Epicoccum were much abundant in TSP samples as compared with PM2.5 and PM10, which may owe to their large spore diameters (>10 μm; McCartney et al., 1993; Yamamoto et al., 2012). In addition, the airborne fungal composition may be related to the release of fungal spores in the air and the release of spores would depend on environmental factors, such as relative humidity, and wind speed (Pasanen et al., 1991; Kildesø et al., 2003). For example, spores of Aspergillus and Penicillium were released at air velocity of 0.5 ms−1, whereas the release of Cladosporium spores required at least a velocity of 1.0 ms−1 (Pasanen et al., 1991). Therefore, Cladosporium showed less abundant and Aspergillus and Penicillium were more abundant in heavy-haze days, which is always associated with less wind as compared with non-haze and light-haze days.
Multiple environmental variables may determine the fungal community composition in the atmosphere. PMs were significantly correlated with fungal diversity and community composition in our study, which indicated that PMs might influence the colonization of different fungi. These results are consistent with previous studies that both PM2.5 and PM10 showed positive correlation with culturable fungi concentration (Haas et al., 2013; Liu et al., 2014). In the present study, RH and temperature were also significantly correlated with fungal diversity and community composition. These environmental factors might impact the fungal spores concentration and fungal growth (Troutt and Levetin, 2001; Rodríguez-Rajo et al., 2005; Erkara et al., 2008; Almaguer et al., 2014), and they might show spatial and temporal patterns. For example, the distinct differences of Alternaria and Cladosporium were observed among the summer and autumn samples, indicating the temporal changes of fungal groups. In addition, some low levels of genera in most of the samples, in our study, showed abundant in a few samples and these irregular changes may owe to the accidental event, such as rainfalls or strong winds (Barnes et al., 2001). Similarly, the airborne fungal communities in Hong Kong were best explained by spatial and temporal factors (Woo et al., 2013). Woo et al. (2013) found that the airborne fungal communities were obviously different among locations (i.e., mid-western urban US locations and Hong Kong locations) and months. In addition, there are many other factors influencing the fungal community in the atmosphere such as wind speed and time of sunlight. Therefore, airborne fungal communities should be surveyed for a long time and their relationships with spatial–temporal factors require further study.
In summary, we have provided an integrated characterization of the airborne fungal community associated with different sizes of particles in Beijing using high-throughput sequencing. Nevertheless, the ecological roles of airborne fungi and their physiological functions remain poorly understood. In the further study, a combination of different technologies, such as traditional culture-based method and metatranscriptomics, may help to answer the pending questions.
DY, JS, LLZ, HW, and XMF sampled. DY performed the laboratory work and wrote the manuscript. TZ performed the part of laboratory work and revised the manuscript. LYY designed the research and revised the manuscript. DY, YQZ, and HYL analyzed the data for the work. All authors contributed significantly in the preparation of the manuscript. All authors approve of the submission of this manuscript.
Conflict of Interest Statement
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.
This research was supported by the National Infrastructure of Microbial Resources (No. NIMR-2015-3), the National Natural Science Foundation of China (NSFC; Nos. 81373452, 31300115, and 81321004), Program (No. 2014AA021504) and National S&T Major Special Project on Major New Drug Innovation (No. 2012ZX09301002-001). LYY is supported by Xiehe Scholar.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.00487
Adams, R. I., Miletto, M., Taylor, J. W., and Bruns, T. D. (2013). Dispersal in microbes: fungi in indoor air are dominated by outdoor air and show dispersal limitation at short distances. ISME J. 7, 1262–1273. doi: 10.1038/ismej.2013.28
Adhikari, A., Sen, M. M., Gupta-Bhattacharya, S., and Chanda, S. (2004). Airborne viable, non-viable, and allergenic fungi in a rural agricultural area of India: a 2-year study at five outdoor sampling stations. Sci. Total Environ. 326, 123–141. doi: 10.1016/j.scitotenv.2003.12.007
Almaguer, M., Aira, M. J., Rodríguez-Rajo, F. J., and Rojas, T. I. (2014). Temporal dynamics of airborne fungi in Havana (Cuba) during dry and rainy seasons: influence of meteorological parameters. Int. J. Biometeorol. 58, 1459–1470. doi: 10.1007/s00484-013-0748-6
Awad, A. H. A., Gibbs, S. G., Tarwater, P. M., and Green, C. F. (2013). Coarse and fine culturable fungal air concentrations in urban and rural homes in Egypt. Int. J. Environ. Res. Public Health. 10, 936–949. doi: 10.3390/ijerph10030936
Barnes, C., Tuck, J., Simon, S., Pacheco, F., Hu, F., and Portnoy, J. (2001). Allergenic materials in the house dust of allergy clinic patients. Ann. Allergy Asthma Immunol. 86, 517–523. doi: 10.1016/S1081-1206(10)62899-2
Beelen, R., Raaschou-Nielsen, O., Stafoggia, M., Andersen, Z. J., Weinmayr, G., Hoffmann, B., et al. (2014). Effects of long-term exposure to air pollution on natural-cause mortality: an analysis of 22 European cohorts within the multicentre ESCAPE project. Lancet 383, 785–795. doi: 10.1016/S0140-6736(13)62158-3
Bellemain, E., Carlsen, T., Brochmann, C., Coissac, E., Taberlet, P., and Kauserud, H. (2010). ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiol. 10:189. doi: 10.1186/1471-2180-10-189
Bowers, R. M., Clements, N., Emerson, J. B., Wiedinmyer, C., Hannigan, M. P., and Fierer, N. (2013). Seasonal variability in bacterial and fungal diversity of the near-surface atmosphere. Environ. Sci. Technol. 47, 12097–12106. doi: 10.1021/es402970s
Bressi, M., Sciare, J., Ghersi, V., Bonnaire, N., Nicolas, J. B., Petit, J. E., et al. (2013). A one-year comprehensive chemical characterisation of fine aerosol (PM2.5) at urban, suburban and rural background sites in the region of Paris (France). Atmos. Chem. Phys. 13, 7825–7844. doi: 10.5194/acp-13-7825-2013
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
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
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cheng, Z., Jiang, J., Fajardo, O., Wang, S., and Hao, J. (2013). Characteristics and health impacts of particulate matter pollution in China (2001–2011). Atmos. Environ. 65, 186–194. doi: 10.1016/j.atmosenv.2012.10.022
Dannemiller, K. C., Lang-Yona, N., Yamamoto, N., Rudich, Y., and Peccia, J. (2014). Combining real-time PCR and next-generation DNA sequencing to provide quantitative comparisons of fungal aerosol populations. Atmos. Environ. 84, 113–121. doi: 10.1016/j.atmosenv.2013.11.036
DeFranco, E., Hall, E., Hossain, M., Chen, A., Haynes, E. N., Jones, D., et al. (2015). Air pollution and stillbirth risk: exposure to airborne particulate matter during pregnancy is associated with fetal death. PLoS ONE 10:e0120594. doi: 10.1371/journal.pone.0120594
Emerson, J. B., Keady, P. B., Brewer, T. E., Clements, N., Morgan, E. E., Awerbuch, J., et al. (2015). Impacts of flood damage on airborne bacteria and fungi in homes after the 2013 Colorado Front Range flood. Environ. Sci. Technol. 49, 2675–2684. doi: 10.1021/es503845j
Erkara, I. P., Asan, A., Yilmaz, V., Pehlivan, S., and Okten, S. S. (2008). Airborne Alternaria and Cladosporium species and relationship with meteorological conditions in Eskisehir City, Turkey. Environ. Monit. Assess. 144, 31–41. doi: 10.1007/s10661-007-9939-0
Evans, J., van Donkelaar, A., Martin, R. V., Burnett, R., Rainham, D. G., Birkett, N. J., et al. (2013). Estimates of global mortality attributable to particulate air pollution using satellite imagery. Environ. Res. 120, 33–42. doi: 10.1016/j.envres.2012.08.005
Fröhlich-Nowoisky, J., Pickersgill, D. A., Després, V. R., and Pöschl, U. (2009). High diversity of fungi in air particulate matter. Proc. Natl. Acad. Sci. U.S.A. 106, 12814–12819. doi: 10.1073/pnas.0811003106
Gao, M., Jia, R., Qiu, T., Han, M., Song, Y., and Wang, X. (2015). Seasonal size distribution of airborne culturable bacteria and fungi and preliminary estimation of their deposition in human lungs during non-haze and haze days. Atmos. Environ. 118, 203–210. doi: 10.1016/j.atmosenv.2015.08.004
Haas, D., Galler, H., Luxner, J., Zarfel, G., Buzina, W., Friedl, H., et al. (2013). The concentrations of culturable microorganisms in relation to particulate matter in urban air. Atmos. Environ. 65, 215–222. doi: 10.1016/j.atmosenv.2012.10.031
Huang, X. H., Bian, Q., Ng, W. M., Louie, P. K., and Yu, J. Z. (2014). Characterization of PM2.5 major components and source investigation in suburban Hong Kong: a one year monitoring study. Aerosol. Air Qual. Res. 14, 237–250. doi: 10.4209/aaqr.2013.01.0020
Hunt, A., Abraham, J. L., Judson, B., and Berry, C. L. (2003). Toxicologic and epidemiologic clues from the characterization of the 1952 London smog fine particulate matter in archival autopsy lung tissues. Environ. Health Persp. 111, 1209. doi: 10.1289/ehp.6114
Kildesø, J., Würtz, H., Nielsen, K. F., Kruse, P., Wilkins, K., Thrane, U., et al. (2003). Determination of fungal spore release from wet building materials. Indoor Air 13, 148–155. doi: 10.1034/j.1600-0668.2003.00172.x
Kim, S. W., Yoon, S. C., Jefferson, A., Ogren, J. A., Dutton, E. G., Won, J. G., et al. (2005). Aerosol optical, chemical and physical properties at Gosan, Korea during Asian dust and pollution episodes in 2001. Atmos. Environ. 39, 39–50. doi: 10.1016/j.atmosenv.2004.09.056
Kõljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F., Bahram, M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481
Kuhn, T., Biswas, S., Fine, P. M., Geller, M., and Sioutas, C. (2005). Physical and chemical characteristics and volatility of PM in the proximity of a light-duty vehicle freeway. Aerosol. Sci. Technol. 39, 347–357. doi: 10.1080/027868290930024
Li, K., Dong, S., Wu, Y., and Yao, M. (2010). Comparison of the biological content of air samples collected at ground level and at higher elevation. Aerobiologia (Bologna) 26, 233–244. doi: 10.1007/s10453-010-9159-x
Li, Y., Fu, H., Wang, W., Liu, J., Meng, Q., and Wang, W. (2015). Characteristics of bacterial and fungal aerosols during the autumn haze days in Xi'an, China. Atmos. Environ. 122, 439–447. doi: 10.1016/j.atmosenv.2015.09.070
Li, Z., Gu, X., Wang, L., Li, D., Xie, Y., Li, K., et al. (2013). Aerosol physical and chemical properties retrieved from ground-based remote sensing measurements during heavy haze days in Beijing winter. Atmos. Chem. Phys. 13, 10171–10183. doi: 10.5194/acp-13-10171-2013
Liu, Z., Li, A., Hu, Z., and Sun, H. (2014). Study on the potential relationships between indoor culturable fungi, particle load and children respiratory health in Xi'an, China. Build. Environ. 80, 105–114. doi: 10.1016/j.buildenv.2014.05.029
Pasanen, A. L., Pasanen, P., Jantunen, M. J., and Kalliokoski, P. (1991). Significance of air humidity and air velocity for fungal spore release into the air. Atmos. Environ. A Gen. Top. 25, 459–462. doi: 10.1016/0960-1686(91)90316-Y
Putaud, J. P., Raes, F., Van Dingenen, R., Brüggemann, E., Facchini, M. C., Decesari, S., et al. (2004). A European aerosol phenomenology—2: chemical characteristics of particulate matter at kerbside, urban, rural and background sites in Europe. Atmos. Environ. 38, 2579–2595. doi: 10.1016/j.atmosenv.2004.01.041
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 (Bologna) 29, 233–248. doi: 10.1007/s10453-012-9276-9
Rittenour, W. R., Ciaccio, C. E., Barnes, C. S., Kashon, M. L., Lemons, A. R., Beezhold, D. H., et al. (2014). Internal transcribed spacer rRNA gene sequencing analysis of fungal diversity in Kansas City indoor environments. Environ. Sci. Process. Impact. 16, 33–43. doi: 10.1039/C3EM00441D
Rodríguez-Rajo, F. J., Iglesias, I., and Jato, V. (2005). Variation assessment of airborne Alternaria and Cladosporium spores at different bioclimatical conditions. Mycol. Res. 109, 497–507. doi: 10.1017/S0953756204001777
Sáenz-de-Santamaría, M., Postigo, I., Gutierrez-Rodríguez, A., Cardona, G., Guisantes, J. A., Asturias, J., et al. (2006). The major allergen of Alternaria alternata (Alt a 1) is expressed in other members of the Pleosporaceae family. Mycoses 49, 91–95. doi: 10.1111/j.1439-0507.2006.01195.x
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, Y., Zhuang, G., Sun, Y., and An, Z. (2006). The variation of characteristics and formation mechanisms of aerosols in dust, haze, and clear days in Beijing. Atmos. Environ. 40, 6579–6591. doi: 10.1016/j.atmosenv.2006.05.066
White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 5:e1000352. doi: 10.1371/journal.pcbi.1000352
Womack, A. M., Artaxo, P. E., Ishida, F. Y., Mueller, R. C., Saleska, S. R., Wiedemann, K. T., et al. (2015). Characterization of active and total fungal communities in the atmosphere over the Amazon rainforest. Biogeosci. Discuss. 12, 7177–7207. doi: 10.5194/bgd-12-7177-2015
Woo, A. C., Brar, M. S., Chan, Y., Lau, M. C., Leung, F. 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
Yamamoto, N., Bibby, K., Qian, J., Hospodsky, D., Rismani-Yazdi, H., Nazaroff, W. W., et al. (2012). Particle-size distributions and seasonal diversity of allergenic and pathogenic fungi in outdoor air. ISME J. 6, 1801–1811. doi: 10.1038/ismej.2012.30
Yamamoto, N., Nazaroff, W. W., and Peccia, J. (2014). Assessing the aerodynamic diameters of taxon-specific fungal bioaerosols by quantitative PCR and next-generation DNA sequencing. J. Aerosol Sci. 78, 1–10. doi: 10.1016/j.jaerosci.2014.08.007
Yang, F., Kawamura, K., Chen, J., Ho, K., Lee, S., Gao, Y., et al. (2016). Anthropogenic and biogenic organic compounds in summertime fine aerosols (PM2.5) in Beijing, China. Atmos. Environ. 124, 166–175. doi: 10.1016/j.atmosenv.2015.08.095
Zhang, J., Mauzerall, D. L., Zhu, T., Liang, S., Ezzati, M., and Remais, J. V. (2010). Environmental health in China: progress towards clean air and safe water. Lancet 375, 1110–1119. doi: 10.1016/S0140-6736(10)60062-1
Keywords: airborne fungi, fungal community composition, PM2.5, PM10, TSP, haze, high-throughput sequencing
Citation: Yan D, Zhang T, Su J, Zhao L-L, Wang H, Fang X-M, Zhang Y-Q, Liu H-Y and Yu L-Y (2016) Diversity and Composition of Airborne Fungal Community Associated with Particulate Matters in Beijing during Haze and Non-haze Days. Front. Microbiol. 7:487. doi: 10.3389/fmicb.2016.00487
Received: 14 January 2016; Accepted: 23 March 2016;
Published: 14 April 2016.
Edited by:Steve Brian Pointing, Auckland University of Technology, New Zealand
Reviewed by:Donnabella Castillo Lacap-Bugler, Auckland University of Technology, New Zealand
Teruya Maki, Kanazawa University, Japan
Copyright © 2016 Yan, Zhang, Su, Zhao, Wang, Fang, Zhang, Liu and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Li-Yan Yu, firstname.lastname@example.org
†These authors have contributed equally to this work.