Original Research ARTICLE
Seasonal Changes of Airborne Bacterial Communities Over Tokyo and Influence of Local Meteorology
- 1Department of Atmospheric Science, Colorado State University, Fort Collins, CO, United States
- 2National Institute of Polar Research, Tachikawa, Japan
- 3Department of Polar Science, School of Multidisciplinary Sciences, SOKENDAI (The Graduate University for Advanced Studies), Tachikawa, Japan
- 4National Research Institute for Earth Science and Disaster Resilience, Storm, Flood and Landslide Research Division, Tsukuba, Japan
In order to study airborne bacterial community dynamics over Tokyo, including fine-scale correlations between airborne microorganisms and meteorological conditions, and the influence of local versus long-range transport of microbes, air samples were collected on filters for periods ranging from 48 to 72 h. The diversity of the microbial community was assessed by next generation sequencing. Predicted source regions of airborne particles, from back trajectory analyses, changed abruptly from the Pacific Ocean to the Eurasian Continent in the beginning of October. However, the microbial community composition and the alpha and beta diversities were not affected by this shift in meteorological regime, suggesting that long-range transport from oceanic or continental sources was not the principal determinant controlling the local airborne microbiome. By contrast, we found a significant correlation between the local meteorology, especially relative humidity and wind speed, and both alpha diversity and beta diversity. Among four potential local source categories (soil, bay seawater, river, and pond), bay seawater and soil were identified as constant and predominant sources. Statistical analyses point toward humidity as the most influential meteorological factor, most likely because it is correlated with soil moisture and hence negatively correlated with the dispersal of particles from the land surface. In this study, we have demonstrated the benefits of fine-scale temporal analyses for understanding the sources and relationships with the meteorology of Tokyo’s “aerobiome.”
Microorganisms are abundant in the atmosphere, with numbers ranging from 104–105 cells m−3 in ambient air on mountain peaks (Bowers et al., 2012) to 106–107 m−3 in desert dust storms (Maki et al., 2017b). The atmosphere is widely recognized as a highly diverse microbiome (Franzetti et al., 2011; Bowers et al., 2011a,b, 2012, 2013; Bertolini et al., 2013; DeLeon-Rodriguez et al., 2013; Woo et al., 2013; Yooseph et al., 2013; Gandolfi et al., 2015; Seifried et al., 2015; Genitsaris et al., 2017; Šantl-Temkiv et al., 2018; Archer et al., 2019). Airborne bacteria and fungi also have the potential to cause human disease, as reviewed by Salem and Gardner (1994). For example, airborne fungi have been identified as the cause of respiratory problems such as asthma after thunderstorms (Grinn-Gofroń and Strzelczak, 2013), and also the lymph node syndrome, Kawasaki disease (Rodó et al., 2014). Endotoxins from airborne bacteria are also associated with health issues (Mueller-Annelling et al., 2004). Airborne microorganisms have impacts not only on human health, but also on climate and microbial biogeography. For example, the plant pathogen Pseudomonas syringae and related phylloplane bacteria have strong ice nucleation ability at 33°C warmer (−5°C) than the homogeneous freezing temperature of cloud droplets composed of pure water (Maki et al., 1974; Lindow et al., 1978; Morris et al., 2008; Hill et al., 2014). Thus, airborne ice nucleating bacteria may promote the formation of ice in clouds, which may modify cloud radiative forcing and promote precipitation. Until recently, most airborne microorganism studies depended on culture-based methods (e.g., Lindemann and Upper, 1985; Shaffer and Lighthart, 1997). However, culture-independent methods using next generation sequencing (NGS) are now widely applied (Franzetti et al., 2011; Bowers et al., 2011a,b, 2012, 2013; Bertolini et al., 2013; DeLeon-Rodriguez et al., 2013; Woo et al., 2013; Yooseph et al., 2013; Gandolfi et al., 2015; Seifried et al., 2015; Genitsaris et al., 2017; Šantl-Temkiv et al., 2018). This technology transforms our ability to describe spatial and temporal variability of many airborne microorganisms, and enables us to compare their temporal variation under different meteorological conditions. However, the relationship between aerial microbial community composition and meteorology is still not well characterized (Gunthe et al., 2016). For example, Bowers et al. (2011a) showed different land-use types, rather than local weather including relative humidity and wind speed, control the airborne bacterial composition in northern Colorado, US, during the early summer season. On the other hand, Gandolfi et al. (2015) showed wind speed and relative humidity affected airborne microbial community structure in two urban sites, in northern Italy. Šantl-Temkiv et al. (2018) showed that the alpha diversity of bacterial communities is positively correlated with air temperature and negatively correlated with relative humidity in western Greenland during mid-summer. Studies are often focused on urban areas, because these contain sources of human pathogens, such as water treatment facilities and densely populated areas. While spatial and temporal variations of bioaerosols in urban environments have been frequently studied in many cities (Brodie et al., 2007; Frohlich-Nowoisky et al., 2009; Bowers et al., 2013; Woo et al., 2013; Gandolfi et al., 2015; Genitsaris et al., 2017), few studies (Bowers et al., 2013) have used consecutive short-period samples to explore correlations between variation of the microbiome and meteorology. Additionally, samples from these studies were taken near ground level (generally several meters above) or on roof tops of lower buildings, with the result that specific local sources very near to sampling sites likely predominate over microorganisms originating regionally (e.g., Archer et al., 2019).
In this study, we investigated: (1) airborne bacterial communities over Tokyo, (2) fine-scale correlations between airborne microorganisms and meteorological conditions, and (3) the influence of local versus long-range transport of microbes, via analysis of a series of consecutive 48–72 h samples collected at a height of 458 m from the tallest communication/observation tower in Japan from summer 2016 to winter 2017.
Materials and Methods
Sampling of Air and Reference Samples
Airborne microorganisms were collected at a 458-m level measurement site (Misumi et al., 2018) on the western side of the Tokyo Skytree (Supplementary Table S1) from August 2016 to February 2017. Eight sets of 2-m-long conductive silicon tubes were fixed to the outside wall of the Tokyo Skytree with their inlets placed under a wind and rain shield, and sterilized inline NILU filter holders (Norwegian Institute for Air Research) attached to the inner ends. These were fitted with pre-combusted (500°C for 2 h) 47-mm-diameter quartz filters with 99.9% retention rate for 0.3 μm particles (Tissuquartz™ Filters, 2500 QAT-UP, Pall). The vacuum sides of the NILU filter holders were connected to low volume samplers (LV-40BW, Sibata Scientific Technology) programmed to sample from one unit at a time. Air samples were collected at a flow rate of 15–30 L min−1 for 48–72 h for each filter, and then the sampler switched to a new filter unit; thus, sampling was continuous throughout the study period. Sampling details are described in Supplementary Table S2. In order to avoid degradation of filter samples during and after sampling, all filter holders were placed in a 4°C refrigerator (JR-N40G, Haier) until retrieval (usually at 2–3 week intervals). Reference environmental (potential source) samples for source estimation were taken around Tokyo Skytree on Jan 3, 2017, and from around the building of the National Institute of Polar Research (NIPR) on June 6, 2016 (Supplementary Table S1, Supplementary Figure S1). We selected these sites because there are very limited open spaces in the Tokyo metropolis. Soil samples were placed in a 5-ml sterile plastic tube containing 2 ml of RNAlater (Thermo Fisher Scientific, MA, USA) using a pre-cleaned stainless-steel spoon. Seawater, river water, and pond water were sampled directly into a 50-ml sterile plastic bottle. All soil and water samples were kept at −20°C prior to DNA extraction.
DNA Extraction, Polymerase Chain Reaction, and DNA Sequencing
In order to avoid contamination, all processes prior to the polymerase chain reaction (PCR) amplification were done in a laminar flow clean bench (PCV-1305BNG3-AG, Hitachi). The clean bench was sanitized with a UV lamp overnight, and pipettes were sterilized in a DNA cross linker (CL-1000, UVP) box inside the clean bench. Genomic DNA in bioaerosols captured on quartz filters was extracted using the FastDNA™ SPIN Kit for Soil (MP Biomedicals, Santa Ana, CA). The quartz filter was initially pulverized during the bead beating step, but in order to maximize the yield of DNA (DNA adsorbs to quartz fibers), all fragments of the filter were carried over until the final elution step. Partial 16S rRNA gene sequences, spanning the V3 and V4 regions, were amplified using the primers Bakt_341F (5´ CCTACGGGNGGCWGCAG 3′) and Bakt_805R (5´ GACTACHVGGGTATCTAATCC 3′), with Illumina overhang adaptor sequences attached to their 5′ ends, by KAPA HiFi HotStart ReadyMix (KAPA Biosystems, MA, USA). PCR conditions comprised 35 cycles of denaturation at 95°C for 30 s, annealing at 55°C for 30 s, and elongation at 72°C for 30 s and an additional final elongation at 72°C for 5 min using a GeneAmp PCR System 9700 (Applied Biosystems, CA, USA). Subsequent clean-up and indexing of PCR amplicons were performed following Illumina methods for 16S metagenomic sequencing library preparation1. All samples were pooled into one flow cell of a MiSeq sequencer (Illumina, San Diego, CA) and sequenced at NIPR using a MiSeq V3 reagent kit.
Exact Sequence Variant Analysis
In order to study sequence differences in greater detail than the conventional 97% operational taxonomic unit (OTU) approach, we used a newly developed DADA2 (1.4) method which, by incorporating an error model, is able to infer sequence reads with single nucleotide resolution (Callahan et al., 2016). Since airborne bacterial DNA is relatively low in concentration and contamination during processing of samples is a known concern (e.g., Salter et al., 2014; de Goffau et al., 2018), air samples and three types of control sequences were analyzed together: (1) an extraction-PCR control using blank filters (07 W4, 14 W4, 16 W3); (2) a PCR control (Negative1, Negative2); and (3) laboratory amplicon contamination (sequences from previous studies analyzed at NIPR: DRA004099–DRA004103). Next, Amplicon Sequence Variants (ASVs) in all contamination controls were removed from further analysis. Taxonomy was assigned by the Naive Bayes Classifier method in the Ribosomal Database Project (RDP) Classifier (Wang et al., 2007) implemented in the DADA2 program against the customized Silva 128 database (data are available at https://doi.org/10.6084/m9.figshare.7772864) (Quast et al., 2013). All potential chimera sequences were removed by DADA2 and all chloroplast and mitochondria were removed manually. Alpha diversities, rarefaction curve of alpha diversities (Chao 1, the Shannon index and the reciprocal of Simpson’s index), unweighted UniFrac distances analysis (Lozupone and Knight, 2005) and unweighted Pair Group Method with Arithmetic Means (UPGMA) were analyzed by QIIME (Caporaso et al., 2012). Chao 1 estimates the total ASV richness; the Shannon index is a general diversity measure that is positively correlated with both diversity and evenness, being sensitive to differences in abundance of rare ASVs; and Simpson’s reciprocal is a measure of evenness, which has a lower bound of 1 for a community composed of only one ASV and, for example, a value of 10 for a community containing 10 equally abundant ASVs (Hill et al., 2003). A phylogenetic tree for unweighted UniFrac distances analysis was constructed using FastTree (Price et al., 2009) implemented in Geneious software, version 10.0.92. Based on rarefaction curve results and number of filtered sequences, sequences were rarefied to 970 for alpha diversity and 870 for beta diversity estimation. All codes and data are available from https://github.com/JUetake/TokyoSkytree-Bacteria. A taxonomy heatmap was created by the R Package “superheat” (Barter and Yu, 2018).
“SourceTracker” is a Bayesian approach program to estimate the proportion of exogenous sequences in a given community that come from possible source environments (Knights et al., 2011). We used the latest version, SourceTracker23, for source estimation of sampled bioaerosols sampled. Twenty-three reference samples collected around Tokyo Skytree (eight sites: Supplementary Figure S1) were analyzed as possible sources of airborne bacteria. In order to focus on local bacteria around the Tokyo area, six Amplicon Sequence Variants (ASVs), which have significant differences (p < 0.001, Mann-Whitney) between Pacific and continental periods (explained in section “Results and Discussion”), were removed from analysis.
Statistical Data Analysis
Correlations between alpha diversities, source contributions, and meteorological data were analyzed by Spearman’s correlation using the R package “ggcorrplot.” Mann-Whitney U tests between Pacific and continental periods were performed with XLSTAT software4. Analysis of similarities (ANOSIM) and permutational multivariate analysis of variance (PERMANOVA) for unweighted UniFrac dissimilarity matrices were performed using PRIMER 7 (Plymouth, UK).
All meteorological data taken by the Automated Meteorological Data Acquisition System (AMeDAS), which is managed by the Japan Meteorological Agency (JMA), are available from http://www.data.jma.go.jp/obd/stats/etrn/index.php. Air temperature, wind speed, sunshine duration (total duration of solar radiation above 0.12 kW m−2), and humidity are taken at the “Edogawa rinkai” station and precipitation measurements were taken at the “Tokyo” station (Supplementary Table S1). Wave height data in Tokyo Bay, which is managed by the Bureau of Port and Harbor, Tokyo Metropolitan Government, are available from http://www.kouwan.metro.tokyo.jp/yakuwari/choui/kako1-index.html. All daily data were transformed to 2–3 day averages to correspond to sampling durations (Supplementary Table S2, Figure 1). The trajectories of air masses arriving at the sampling location were calculated by HYSPLIT (Cohen et al., 2015) using data from the GDAS 0.5 model (lat = 35.7101, lon = 139.8107, height = 458, duration = 72) and were started from 12:00 (UTC) of each day (Figure 2). All samples were categorized as either continental, meaning of Eurasian continent origin, or Pacific, based on their latitude/longitude location 72 h prior to arriving at Skytree.
Figure 1. Seasonal changes of meteorological factors (temperature: pink, sunshine duration: orange, precipitation: blue, humidity: light blue, and wind speed: light green) from nearest automatic weather stations, wave height (purple), and alpha diversities (gray) during sampling periods.
Figure 2. 72 h HYSPLIT back trajectories during sampling period. Height above 1,600 m is shown in blue.
Results and Discussion
Long-Range Effects of Air Mass Shifts Upon Bacterial Composition
Bacteria, which have round shapes and typical diameters/lengths of ~1–2 μm, are readily emitted from surfaces under certain conditions, such as by wind, mechanical disturbances, and bubble bursting, and may have aerial residence times from seconds to months (Burrows et al., 2009; Arnold et al., 2016), during which time they may be transported over great distances before their removal by both wet and dry deposition (Prospero et al., 2005; Barberán et al., 2014; Maki et al., 2017a). On non-dusty days, bacterial concentrations near ground level (10 m a.g.l.) are 10–100 times higher than in air at higher levels (Maki et al., 2017a), indicating that near-ground levels are strongly influenced by local sources. The average daytime atmospheric boundary layer height in Tokyo generally exceeds 500 m, even in wintertime, but is generally just below 500 m at nighttime over all seasons (Dien et al., 2017). Therefore, the Tokyo Skytree sampling site would likely sample boundary layer air primarily during daytime, and air from the residual layer and potentially the lower free troposphere during nighttime. Since samples were obtained continuously over several days, the deposited particulate matter in any one sample will be derived from this mix of air masses.
In this study, an average of 5,862 sequences per sample (maximum: 19,913, minimum: 1,068) were used for analyses (Supplementary Table S2). Adequate coverage is required for reliable estimation by the indices (Hill et al., 2003), and indeed the depth of sequencing was enough for all three to stabilize (by 1,000 sequences; Supplementary Figure S2). Temporal variation of phylum Proteobacteria (mean 51.4% of relative abundances) was by far the most common followed by Firmicutes (13.6%), Cyanobacteria (7.9%), Actinobacteria (7.7%), Bacteroidetes (5.2%), and Acidobacteria (3.0%) (Figure 3). ANOSIM analyses on unweighted UniFrac distances (Supplementary Figure S3) clearly clustered into Pacific and continental sources (ANOSIM Global R = 0.337, p < 0.001: Supplementary Table S3), as predicted using HYSPLIT back trajectories.
Of all the phyla, only Parcubacteria (2.6%) exhibited seasonal changes, accounting for a larger fraction during August and September. A finer resolution of taxonomy (Supplementary Figure S2) shows that Paraburkholderia (mean 6.3%) had the highest mean relative abundance throughout the entire sampling period, followed by Sphingomonas (3.5%), Chroococcidiopsis (3.3%), and Bacillus (3.2%). Relative abundances of some genera were significantly (p < 0.05, Mann-Whitney) different in early October (Supplementary Table S4), when the HYSPLIT model indicated that the source region was shifted from Pacific to continental between October 2nd and 4th. At the genus level, the relative abundance of Paraburkholderia is significantly higher in the continental period, and that of an unidentified genus in the phylum Parcubacteria was higher in Pacific period (Supplementary Figure S4). Therefore, these genera are presumed to have originated from non-local sources, indicating that synoptic air mass movement influenced bacterial distributions.
To obtain more detail of these potentially long-range-transported bacteria, we investigated their ASV level taxonomy. Six ASVs (mostly members of the Parcubacteria) out of a total of 3,257 ASVs had significant differences in their abundances (p < 0.001, Mann-Whitney) between Pacific and continental periods (Supplementary Table S4, Supplementary Figure S5). The estimated sources of five of these ASVs (all the Parcubacteria and an Opitutus sp. from the Verrucomicrobia) that were abundant in the Pacific period, and hence were expected to be pelagic, were actually likely from freshwater environments according to their closest matches using BLAST searches (95.1–99.8% identity: Supplementary Table S5) and metametaDB analysis (25–78% MHI: Supplementary Figure S6). The sixth ASV (Paraburkholderia sp.) that was abundant in the period affected by continental air, and is mostly from soils (Supplementary Figure S6).
Their non-oceanic origin is also supported by the absence of ocean bacteria, such as bacteria belonging to the SAR group or typical marine groups (e.g., Oceanospirillales). Oceanospirillales are known to be readily aerosolized and carried long distances, since it was recorded in rainfall in the central Pyrenees mountains (Spain) by Cáliz et al. (2018), and at the Mt. Bachelor Observatory in North America (2.8 km above sea level) by Smith et al. (2013). This group seems to readily aerosolize and be transported long distances.
Thus, despite the HYSPLIT results, BLAST and metametaDB analysis for Parcubacteria and Opitutus indicate minimal inputs from the Pacific Ocean. This can be explained by a combination of low emissions from the ocean surface and higher inputs from regional freshwater. In previous studies, similarly lower contributions of ocean bacteria were reported. For example, the number of cultivable bacteria found in near-surface air during non-dust events over the ocean is generally very low (Prospero et al., 2005; Griffin et al., 2006), and the occurrence of marine bacteria in the near-surface air of a coastal city in Greece was found to be rare (Genitsaris et al., 2017).
One of the six ASVs, belonging to the genus Paraburkholderia (ASV137) and abundant in the continental period, was indicated to be a soil bacterium based on closest relative matches in BLAST searches (100% of identity: Supplementary Table S5) and metametaDB analysis (44% MHI: Supplementary Figure S6). Exactly the same sequences have been found in soil samples in many countries including China (e.g., KU323602, KX351056, KY427125). Since these were transported by westerly winds, which prevailed in the later sampling period when air masses mostly came from the Eurasian Continent, this ASV possibly originates from Chinese soil. Although some bacteria are known to be transported to Japan from arid regions in China during Asian dust events in March and April in Japan (Iwasaka et al., 2009), the genus Paraburkholderia (and its higher taxonomy order Burkholderiales) was not found in ground observations during Asian dust days (Park et al., 2016), in the upper Asian dust layer (Maki et al., 2017a) or in dust coming from the Gobi Desert, one of their potential source areas (Maki et al., 2017b). Therefore, ASV137 is not likely to be of Asian dust origin. Furthermore, Bacillus subtilis, which is common in Asian dust studies (Maki et al., 2010, 2017a), was found only once over the entire sampling period. Some Bacillus ASVs are relatively similar to uncultured bacteria from the Taklimakan Desert (AB696509 and AB696498); however, these were not perfectly matched by BLAST (ASV399: 99.1%, ASV6088: 98.6%, ASV1582: 98.1%). While six ASVs had long-range origins, the majority of ASVs were of local origin, and therefore were influenced by local meteorological conditions.
Local Meteorological Effects Upon Sources of Bacteria
Alpha diversities were relatively higher in August, October, and at the end of January and in the beginning of February. In contrast, there was also an extended period of consistently low alpha diversity during September (Supplementary Table S1, Figure 1). Correlation coefficients with the six local meteorological and ocean factors (Figure 4) show that, among meteorological and ocean factors, humidity and precipitation had the most significant correlations (both negative) with all diversity measures, while wind speed and sunshine duration had significant positive correlations with the alpha diversities. The same effects of local meteorology were observed in the PERMANOVA test of beta diversity on the effect of meteorological data (PERMANOVA, p < 0.005, Supplementary Table S6). While many previous studies did not analyze the cross-correlated meteorological variables (Burrows et al., 2009), multivariate statistics enable these relationships to be revealed (Bowers et al., 2013). We found that both simple correlations and multivariate analyses showed the same relationships, with humidity and wind speed being effective meteorological factors correlated with bacterial composition and diversity.
Figure 4. Spearman’s correlation between alpha diversities and meteorological and ocean factors. A red/orange circle shows positive correlation and a blue circle shows negative correlation. A cross “X” on the circle indicates no significance (p < 0.05).
The effect of relative humidity changes on the composition of the airborne bacterial community has been reported previously. For example, canonical correspondence analyses of bacterial community structure in urban bioaerosols in Italy showed that both relative humidity and wind speed affected airborne bacterial community structure (Gandolfi et al., 2015), and a study in western Oregon, U.S.A, showed that the concentrations of airborne bacteria were positively correlated with temperature but negatively correlated with relative humidity (Tong and Lighthart, 2000). In the present study, relative humidity increased during precipitation events and remained high, before decreasing with an improvement in the weather (Figure 1). Relative humidity was especially high in late August – September, when it was strongly affected by a passing typhoon and an autumn rain event; alpha diversities were, correspondingly, low in the same period. High humidity tends to keep the surface of the ground wet, and the bonding force of surface tension will keep particles attached to the surface. This binding effect reduces with drying (Jones and Harrison, 2004; Šantl-Temkiv et al., 2018).
Some potential source areas were also affected by humidity. Sourcetracker2 analysis shows the estimated contributions to Tokyo Skytree air samples from each potential source (Figure 5 for each category, Supplementary Figure S7). Total percentage contributions from the tested sources had a maximum of 52%, a minimum of 0.1%, and a mean of 6.8%. Correlations between source contributions and environmental factors (Figure 6) show that, among the six environmental factors, humidity had a significant negative correlation for all soil sites (Kiyari, Ueno, Meiji and Tachikawa). This correlation was highest in the middle of winter, the driest season during the sampling period (Figures 1, 6).
Figure 5. Seasonal change of estimated contribution from potential source types (bay, soil, pond, and river) by source tracking analysis.
Figure 6. Spearman’s correlation between estimated source contribution and meteorological factors and wave height in Tokyo Bay. A red/orange circle shows positive correlation and a blue circle shows negative correlation. Cross “X” on the circle indicates no significance (p < 0.05).
Precipitation, which obviously strongly affects the humidity and surface bonding forces, was, however, not found to have a consistent or strong relationship with airborne bacterial diversity (correlation coefficient: −0.3 to −0.4) and composition (Psuedo-F: 1.48, p = 0.023). While alpha diversities tended to decrease during precipitation events (Supplementary Figure S8: blue broken lines), exceptions occurred during heavy precipitation events resulting from a typhoon that set historical precipitation records and during some stationary fronts (Supplementary Figure S8: red broken lines). The common expectation, and what occurred during the majority of lighter precipitation events, is that airborne bacterial communities would be reduced by precipitation due to scavenging (wet deposition) (Tong and Lighthart, 2000; Yue et al., 2016) and wetting of the soil surface. However, alpha diversities slightly increased during the historical heavy rains on Aug. 21–22, Sep. 12–13, and Sep. 22–23 (Supplementary Figure S8: red broken line). For example, rain rates of 107.5 mm h−1 were reported in Ome (30 km west from the site) on Aug. 22 during typhoon passage, and 63.5 mm h−1 in Yokoshibahikari (80 km east) under a stationary front condition. Diversity increase during such heavy rain events may be explained by bioaerosol generation caused by the impact of large rain drops on plant, soil, and built surfaces (Tong and Lighthart, 2000; Wang et al., 2016; Joung et al., 2017; Kim et al., 2019), and in this case potential sources would be near the sampling site (e.g., even the outside walkway and wall of the tower).
Wind speed was positively correlated with airborne microbial concentrations in previous studies (Bovallius et al., 1978; Lighthart, 1997; Mouli et al., 2005). As well, dynamic water surfaces enhance aerosol generation into the atmosphere (Aller et al., 2005; Dueker et al., 2011; DeLeon-Rodriguez et al., 2013; Maki et al., 2017a; Cáliz et al., 2018). Wind speed is highest in late August, and was significantly correlated with wave height in Tokyo Bay (Figure 1). High wind speeds over the bay’s surface resulted in rougher waters, which would be expected to correlate with more aerosol generation via the bubble bursting process (Burrows et al., 2009) during wave-breaking, especially at the shore. Among the three potential sources, the bay sample had the most significant positive correlation of source contribution with wind speed and wave height. Wind speed was also weakly correlated with two of the soil sites (Meiji and Tachikawa), suggesting that dispersal of bacteria from the soil surface was enhanced by stronger winds.
Many previous studies showed temperature, which follows seasonal cycles in temperate regions, is also a factor controlling airborne bacterial emissions and community composition (Tong and Lighthart, 2000; Jones and Harrison, 2004; Bowers et al., 2013; Genitsaris et al., 2017). In contrast, temperature was not found to be a significant factor in this study.
The effects of local meteorological conditions upon airborne microorganisms are complicated, because recent studies have suggested that airborne bacterial communities are more affected by the nature of local sources than by changes in local meteorological conditions (Bowers et al., 2011a, 2013). This study, and other works, found raised relative humidity and rainfall were the principal drivers of primary bioaerosols and ice nucleating particle emissions (Huffman et al., 2013; Prenni et al., 2013; Tobo et al., 2013; Wright et al., 2014; Mason et al., 2015). Bacterial diversity and community structure are clearly influenced by both local ecotypes and meteorological factors.
In this study, air mass shifts detected by HYSPLIT analysis indicated limited long-range contributions. Estimated contributions from local sources, from seawater in Tokyo Bay and soil, were consistently found and predominated over other known sources over the sampling period (Figure 5). Since the samples used for source estimation were limited in type, number of sites, and season, our reference samples only explained, on average, <6.9% of the total relative contribution from multiple sources. This implies that we should include more diverse potential source types in future studies, such as plant surfaces (Šantl-Temkiv et al., 2012) and sites with exposed animal (cow, dog, and human) feces (Bowers et al., 2013). Adding more, and more varied, reference samples will likely enable more accurate and comprehensive source estimations.
Since a high proportion of airborne bacteria has been shown to remain viable after continental scale long-range transport (Hara and Zhang, 2012), deposition of airborne bacterial communities plays a crucial role in the local microbiome (Barberán et al., 2014; Honeyman et al., 2018). All Tokyo sites would be influenced and perhaps modified by the deposition of such airborne bacteria. However, deposited microbes will only persist if their new niche is favorable. For example, the genus Paraburkholderia, which dominated in the Continental period, was found only in one out of four soil sites (at NIPR).
In order to obtain enough DNA for analysis, we obtained samples that were integrated over 2–3 days. However, shorter sampling intervals may be helpful for understanding the effect of radical meteorological changes upon microbial communities, such as heavy rain enhancement of microbial emissions versus precipitation scavenging in other rain events. Recently, much higher flow rate and lower cost sampling devices have been proposed (Šantl-Temkiv et al., 2017) for use in the field (Šantl-Temkiv et al., 2018). Such capabilities may offer the ability to understand not only hourly changes of specific microbes, but also offer sufficient sample sizes to permit metagenomic analyses of potential microbiome function.
Consecutive 48–72 h samples from the tallest tower in Japan, combined with NGS analyses, were used to study airborne microbial communities from summer to winter. In this study, air mass shifts detected by HYSPLIT analysis indicated limited long-range effects on microbial populations. Our results showed only a limited number of ASVs that could potentially be associated with long-distance transport of bacteria. Further, relative abundances of most of the ASVs were not affected by abrupt air mass changes, indicating they were likely to have been from local sources. Source tracking analysis indicated that local inputs, from soil and the seawater in Tokyo Bay, were consistently found and were the predominant identified sources over the study period. Correlation and PERMANOVA analyses indicated that humidity and wind speed were key factors affecting bacterial alpha and beta diversity, and hence, these appear to be the controlling factors on emissions of bacteria from bay seawater and soil around the Tokyo Skytree. The combination of NGS, reliable metadata, and powerful statistical tools allowed us to study the airborne microbiome near the Tokyo Skytree with unprecedented resolution, and the accumulation of such knowledge from many environments across the world could provide a more comprehensive understanding of factors determining local, regional, and global airborne microbiomes.
Publicly available datasets were analyzed in this study. This data can be found here: https://doi.org/10.6084/m9.figshare.7772864.
JU designed and managed the research. JU, YT, YU, and RM set up the laboratory on the Tokyo Skytree. JU was responsible for sampling, with support from YT and YU. JU performed the bioinformatics and statistical analyses of the sequencing data. JU, YT, TH, PD, and SK wrote the manuscript. All authors read, edited, and approved the final manuscript.
This study was partially supported by crowdfunding operated by academist Co., Ltd. JU was supported for further analyses of data under the National Science Foundation grant AGS1660486.
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.
The authors thank Mr. Kenichi Watanabe, Ms. Ayumi Akiyoshi, and Ms. Mizuho Mori for assistance in DNA extraction, PCR, and sequencing. They also thank the two reviewers for their helpful suggestions and gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and the READY website (https://www.ready.noaa.gov/index.php) used in this publication. This manuscript has been released as a Pre-Print at bioRxiv: https://www.biorxiv.org/content/10.1101/542001v1.article-metrics. The authors specially thank all supporters of JU’s project (https://academist-cf.com/projects/36).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01572/full#supplementary-material
Aller, J. Y., Kuznetsova, M. R., Jahns, C. J., and Kemp, P. F. (2005). The sea surface microlayer as a source of viral and bacterial enrichment in marine aerosols. J. Aerosol Sci. 36, 801–812. doi: 10.1016/j.jaerosci.2004.10.012
Archer, S., Lee, K., Caruso, T., Lee, C., Maki, T., Cowan, D., et al. (2019). Microbial dispersal limitation to isolated soil microbial habitats of Antarctica. Nat. Microbiol. 4, 925–932. doi: 10.1038/s41564-019-0370-4
Arnold, S. R., Law, K. S., Brock, C. A., Thomas, J. L., Starkweather, S. M., von Salzen, K., et al. (2016). Arctic air pollution: challenges and opportunities for the next decade. Elem. Sci. Anth. 4:000104. doi: 10.12952/journal.elementa.000104
Barberán, A., Henley, J., Fierer, N., and Casamayor, E. O. (2014). Structure, inter-annual recurrence, and global-scale connectivity of airborne microbial communities. Sci. Total Environ. 487, 187–195. doi: 10.1016/j.scitotenv.2014.04.030
Barter, R. L., and Yu, B. (2018). Superheat: an R package for creating beautiful and extendable heatmaps for visualizing complex data. J. Comput. Graph. Stat. 27, 910–922. doi: 10.1080/10618600.2018.1473780
Bertolini, V., Gandolfi, I., Ambrosini, R., Bestetti, G., Innocente, E., Rampazzo, G., et al. (2013). Temporal variability and effect of environmental variables on airborne bacterial communities in an urban area of northern Italy. Appl. Microbiol. Biotechnol. 97, 6561–6570. doi: 10.1007/s00253-012-4450-0
Bovallius, A., Bucht, B., Roffey, R., and Anas, P. (1978). Three year investigation of the natural airborne bacterial flora at four localities in Sweden. Appl. Environ. Microbiol. 35, 847–852. doi: 10.5194/acp-9-9263-2009
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
Bowers, R. M., McCubbin, I. B., Hallar, A. G., and Fierer, N. (2012). Seasonal variability in airborne bacterial communities at a high-elevation site. Atmos. Environ. 50, 41–49. doi: 10.1016/j.atmosenv.2012.01.005
Bowers, R. M., McLetchie, S., Knight, R., and Fierer, N. (2011a). 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
Bowers, R. M., Sullivan, A. P., Costello, E. K., Collett, J. L., Knight, R., and Fierer, N. (2011b). Sources of bacteria in outdoor air across cities in the midwestern United States. Appl. Environ. Microbiol. 77, 6350–6356. doi: 10.1128/AEM.05498-11
Brodie, E. L., DeSantis, T. Z., Parker, J. P. M., Zubietta, I. X., Piceno, Y. M., and Andersen, G. L. (2007). Urban aerosols harbor diverse and dynamic bacterial populations. Proc. Natl. Acad. Sci. USA 104, 299–304. doi: 10.1073/pnas.0608255104
Burrows, S. M., Elbert, W., Lawrence, M. G., and Pöschl, U. (2009). Bacteria in the global atmosphere - part 1: review and synthesis of literature data for different ecosystems. Atmos. Chem. Phys. 9, 9263–9280. doi: 10.5194/acp-9-9263-2009
Cáliz, J., Triadó-margarit, X., Camarero, L., and Casamayor, E. O. (2018). A long-term survey unveils strong seasonal patterns in the airborne microbiome coupled to general and regional atmospheric circulations. Proc. Natl. Acad. Sci. USA 115, 1–6. doi: 10.1073/pnas.1812826115
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
Cohen, M. D., Ngan, F., Rolph, G. D., Draxler, R. R., Stunder, B. J. B., and Stein, A. F. (2015). NOAA’s HYSPLIT atmospheric transport and dispersion modeling system. Bull. Am. Meteorol. Soc. 96, 2059–2077. doi: 10.1175/bams-d-14-00110.1
de Goffau, M. C., Lager, S., Salter, S. J., Wagner, J., Kronbichler, A., Charnock-Jones, D. S., et al. (2018). Recognizing the reagent microbiome. Nat. Microbiol. 3, 851–853. doi: 10.1038/s41564-018-0202-y
DeLeon-Rodriguez, N., Lathem, T. L., Rodriguez-R, L. M., Barazesh, J. M., Anderson, B. E., Beyersdorf, A. J., et al. (2013). Microbiome of the upper troposphere: species composition and prevalence, effects of tropical storms, and atmospheric implications. Proc. Natl. Acad. Sci. USA 110, 2575–2580. doi: 10.1073/pnas.1212089110
Dien, N. T., Hirai, Y., and Sakai, S. I. (2017). Correlation between atmospheric boundary layer height and polybrominated diphenyl ether concentrations in air. Environ. Sci. Technol. 51, 356–364. doi: 10.1021/acs.est.6b03004
Dueker, M. E., Weathers, K. C., O’Mullan, G. D., Juhl, A. R., and Uriarte, M. (2011). Environmental controls on coastal coarse aerosols: implications for microbial content and deposition in the near-shore environment. Environ. Sci. Technol. 45, 3386–3392. doi: 10.1021/es1035128
Franzetti, A., Gandolfi, I., Gaspari, E., Ambrosini, R., and Bestetti, G. (2011). Seasonal variability of bacteria in fine and coarse urban air particulate matter. Appl. Microbiol. Biotechnol. 90, 745–753. doi: 10.1007/s00253-010-3048-7
Frohlich-Nowoisky, J., Pickersgill, D. a., Despres, V. R., and Poschl, U. (2009). High diversity of fungi in air particulate matter. Proc. Natl. Acad. Sci. USA 106, 12814–12819. doi: 10.1073/pnas.0811003106
Gandolfi, I., Bertolini, V., Bestetti, G., Ambrosini, R., Innocente, E., Rampazzo, G., et al. (2015). Spatio-temporal variability of airborne bacterial communities and their correlation with particulate matter chemical composition across two urban areas. Appl. Microbiol. Biotechnol. 99, 4867–4877. doi: 10.1007/s00253-014-6348-5
Genitsaris, S., Stefanidou, N., Katsiapi, M., Kormas, K. A., Sommer, U., and Moustaka-Gouni, M. (2017). Variability of airborne bacteria in an urban Mediterranean area (Thessaloniki, Greece). Atmos. Environ. 157, 101–110. doi: 10.1016/j.atmosenv.2017.03.018
Griffin, D. W., Westphal, D. L., and Gray, M. A. (2006). Airborne microorganisms in the African desert dust corridor over the mid-Atlantic ridge, ocean drilling program, leg 209. Aerobiologia 22, 211–226. doi: 10.1007/s10453-006-9033-z
Gunthe, S. S., Hoor, P., Burrows, S. M., Kampf, C. J., Lang-Yona, N., Després, V. R., et al. (2016). Bioaerosols in the earth system: climate, health, and ecosystem interactions. Atmos. Res. 182, 346–376. doi: 10.1016/j.atmosres.2016.07.018
Hill, T. C. J., 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
Hill, T. C. J., Walsh, K. A., Harris, J. A., and Moffett, B. F. (2003). Using ecological diversity measures with bacterial communities. FEMS Microbiol. Ecol. 43, 1–11. doi: 10.1111/j.1574-6941.2003.tb01040.x
Honeyman, A. S., Day, M. L., and Spear, J. R. (2018). Regional fresh snowfall microbiology and chemistry are driven by geography in storm-tracked events, Colorado, USA. PeerJ. 6:e5961. doi: 10.7717/peerj.5961
Huffman, J. A., Prenni, A. J., Demott, P. J., Pöhlker, C., Mason, R. H., Robinson, N. H., et al. (2013). High concentrations of biological aerosol particles and ice nuclei during and after rain. Atmos. Chem. Phys. 13, 6151–6164. doi: 10.5194/acp-13-6151-2013
Iwasaka, Y., Shi, G. Y., Yamada, M., Kobayashi, F., Kakikawa, M., Maki, T., et al. (2009). Mixture of Kosa (Asian dust) and bioaerosols detected in the atmosphere over the Kosa particles source regions with balloon-borne measurements: possibility of long-range transport. Air Qual. Atmos. Health 2, 29–38. doi: 10.1007/s11869-009-0031-5
Jones, A. M., and Harrison, R. M. (2004). The effects of meteorological factors on atmospheric bioaerosol concentrations - a review. Sci. Total Environ. 326, 151–180. doi: 10.1016/j.scitotenv.2003.11.021
Kim, S., Park, H., Gruszewski, H. A., Schmale, D. G., and Jung, S. (2019). Vortex-induced dispersal of a plant pathogen by raindrop impact. Proc. Natl. Acad. Sci. USA 116, 4917–4922. doi: 10.1073/pnas.1820318116
Knights, D., Kuczynski, J., Charlson, E. S., Zaneveld, J., Mozer, M. C., Collman, R. G., et al. (2011). Bayesian community-wide culture-independent microbial source tracking. Nat. Methods 8, 761–765. doi: 10.1038/nmeth.1650
Maki, T., Hara, K., Iwata, A., Lee, K. C., Kawai, K., Kai, K., et al. (2017a). Variations in airborne bacterial communities at high altitudes over the Noto Peninsula (Japan) in response to Asian dust events. Atmos. Chem. Phys. 17, 11877–11897. doi: 10.5194/acp-17-11877-2017
Maki, T., Kurosaki, Y., Onishi, K., Lee, K. C., Pointing, S. B., Jugder, D., et al. (2017b). Variations in the structure of airborne bacterial communities in Tsogt-Ovoo of Gobi desert area during dust events. Air Qual. Atmos. Health 10, 249–260. doi: 10.1007/s11869-016-0430-3
Maki, T., Susuki, S., Kobayashi, F., Kakikawa, M., Tobo, Y., Yamada, M., et al. (2010). Phylogenetic analysis of atmospheric halotolerant bacterial communities at high altitude in an Asian dust (KOSA) arrival region, Suzu City. Sci. Total Environ. 408, 4556–4562. doi: 10.1016/j.scitotenv.2010.04.002
Mason, R. H., Si, M., Li, J., Chou, C., Dickie, R., Toom-Sauntry, D., et al. (2015). Ice nucleating particles at a coastal marine boundary layer site: correlations with aerosol type and meteorological conditions. Atmos. Chem. Phys. 15, 12547–12566. doi: 10.5194/acp-15-12547-2015
Misumi, R., Miura, K., Tobo, Y., Iwanami, K., Uetake, J., Uji, Y., et al. (2018). Characteristics of droplet size distributions in low-level stratiform clouds observed from Tokyo Skytree. J. Meteorol. Soc. Japan. Ser. II 96, 405–413. doi: 10.2151/jmsj.2018-040
Morris, C. E., Sands, D. C., Vinatzer, B. A., Glaux, C., Guilbaud, C., Buffière, A., et al. (2008). The life history of the plant pathogen pseudomonas syringae is linked to the water cycle. ISME J. 2, 321–334. doi: 10.1038/ismej.2007.113
Mouli, P., Mohan, S., and Reddy, S. (2005). Assessment of microbial(bacteria) concentrations of ambient air at semi-arid urban region: influence of meteorological factors. Appl. Ecol. Environ. Res. 3, 139–149. doi: 10.15666/aeer/0302_139149
Mueller-Annelling, L., Avol, E., Peters, J. M., and Thorne, P. S. (2004). Ambient endotoxin concentrations in PM10 from Southern California. Environ. Health Perspect. 112, 583–588. doi: 10.1289/ehp.6552
Park, J., Ichijo, T., Nasu, M., and Yamaguchi, N. (2016). Investigation of bacterial effects of Asian dust events through comparison with seasonal variability in outdoor airborne bacterial community. Sci. Rep. 6, 1–3. doi: 10.1038/srep35706
Prenni, A. J., Tobo, Y., Garcia, E., DeMott, P. J., Huffman, J. A., McCluskey, C. S., et al. (2013). The impact of rain on ice nuclei populations at a forested site in Colorado. Geophys. Res. Lett. 40, 227–231. doi: 10.1029/2012GL053953
Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). Fasttree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077
Prospero, J. M., Blades, E., Mathison, G., and Naidu, R. (2005). Interhemispheric transport of viable fungi and bacteria from Africa to the Caribbean with soil dust. Aerobiologia 21, 1–19. doi: 10.1007/s10453-004-5872-7
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596. doi: 10.1093/nar/gks1219
Rodó, X., Curcoll, R., Robinson, M., Ballester, J., Burns, J. C., Cayan, D. R., et al. (2014). Tropospheric winds from northeastern China carry the etiologic agent of Kawasaki disease from its source to Japan. Proc. Natl. Acad. Sci. USA 111, 7952–7957. doi: 10.1073/pnas.1400380111
Salter, S. J., Cox, M. J., Turek, E. M., Calus, S. T., Cookson, W. O., Moffatt, M. F., et al. (2014). Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 12, 1–12. doi: 10.1186/s12915-014-0087-z
Šantl-Temkiv, T., Amato, P., Gosewinkel, U., Thyrhaug, R., Charton, A., Chicot, B., et al. (2017). High-flow-rate impinger for the study of concentration, viability, metabolic activity, and ice-nucleation activity of airborne bacteria. Environ. Sci. Technol. 51, 11224–11234. doi: 10.1021/acs.est.7b01480
Šantl-Temkiv, T., Finster, K., Hansen, B. M., Nielsen, N. W., and Karlson, U. G. (2012). The microbial diversity of a storm cloud as assessed by hailstones. FEMS Microbiol. Ecol. 81, 684–695. doi: 10.1111/j.1574-6941.2012.01402.x
Šantl-Temkiv, T., Gosewinkel, U., Starnawski, P., Lever, M., and Finster, K. (2018). Aeolian dispersal of bacteria in Southwest Greenland: their sources, abundance, diversity and physiological states. FEMS Microbiol. Ecol. 94, 1–10. doi: 10.1093/femsec/fiy031
Shaffer, B. T., and Lighthart, B. (1997). Survey of culturable airborne bacteria at four diverse locations in Oregon: urban, rural, forest, and coastal. Microb. Ecol. 34, 167–177. doi: 10.1007/s002489900046
Smith, D. J., Timonen, H. J., Jaffe, D. A., Griffin, D. W., Birmele, M. N., Perry, K. D., et al. (2013). Intercontinental dispersal of bacteria and archaea by transpacific winds. Appl. Environ. Microbiol. 79, 1134–1139. doi: 10.1128/AEM.03029-12
Tobo, Y., Prenni, A. J., DeMott, P. J., Huffman, J. A., McCluskey, C. S., Tian, G., et al. (2013). Biological aerosol particles as a key determinant of ice nuclei populations in a forest ecosystem. J. Geophys. Res. Atmos. 118, 10100–10110. doi: 10.1002/jgrd.50801
Tong, Y., and Lighthart, B. (2000). The annual bacterial particle concentration and size distribution in the ambient atmosphere in a rural area of the Willamette Valley, Oregon. Aerosol Sci. Technol. 32, 393–403. doi: 10.1080/027868200303533
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve 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
Woo, A. C., 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
Wright, T. P., Hader, J. D., McMeeking, G. R., and Petters, M. D. (2014). High relative humidity as a trigger for widespread release of ice nuclei. Aerosol Sci. Technol. 48, i–v. doi: 10.1080/02786826.2014.968244
Yooseph, S., Andrews-Pfannkoch, C., Tenney, A., McQuaid, J., Williamson, S., Thiagarajan, M., et al. (2013). A metagenomic framework for the study of airborne microbial communities. PLoS One 8:e81862. doi: 10.1371/journal.pone.0081862
Keywords: bioaerosol, air, DNA, airborne microbiome, urban microbiome, next generation sequencing
Citation: Uetake J, Tobo Y, Uji Y, Hill TCJ, DeMott PJ, Kreidenweis SM and Misumi R (2019) Seasonal Changes of Airborne Bacterial Communities Over Tokyo and Influence of Local Meteorology. Front. Microbiol. 10:1572. doi: 10.3389/fmicb.2019.01572
Edited by:Mark Alexander Lever, ETH Zürich, Switzerland
Reviewed by:Kai Waldemar Finster, Aarhus University, Denmark
Eyal Rahav, Israel Oceanographic and Limnological Research, Israel
Copyright © 2019 Uetake, Tobo, Uji, Hill, DeMott, Kreidenweis and Misumi. 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) and the copyright owner(s) 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: Jun Uetake, firstname.lastname@example.org