Shu-Yi-Dan Zhou
Hu Li
Madeline Giles
Roy Neilson
Xiao-ru Yang
Jian-qiang Su- 1Key Laboratory of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Center for Excellence in Regional Atmospheric Environment, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen, China
- 4Ecological Sciences, The James Hutton Institute, Dundee, United Kingdom
The phyllosphere is populated by numerous microorganisms. Microbes from the wider environment, i.e., air and soil, are considered key contributors to phyllosphere microbial communities, but their contribution is unclear. This study seeks to address this knowledge gap by controlling the movement of microbes along the air-phyllosphere-soil continuum. Customized equipment with dual chambers was constructed that permitted airflow to enter the first chamber while the second chamber recruited filtered microbe-free air from the initial chamber. Allium schoenoprasum (chive) and Sonchus oleraceus (sow thistle) were cultivated in both chambers, and the microbial communities from air, phyllosphere, and soil samples were characterized. Shares of microbial OTUs in the equipment suggested a potential interconnection between the air, phyllosphere, and soil system. Fast expectation-maximization microbial source tracking (FEAST) suggested that soil was the major source of airborne microbial communities. In contrast, the contribution of airborne and soil microbes to phyllosphere microbial communities of either A. schoenoprasum or S. oleraceus was limited. Notably, the soilborne microbes were the only environmental sources to phyllosphere in the second chamber and could affect the composition of phyllosphere microbiota indirectly by air flow. The current study demonstrated the possible sources of phyllosphere microbes by controlling external airborne microbes in a designed microcosm system and provided a potential strategy for recruitment for phyllosphere recruitment.
Introduction
The phyllosphere represents one of the most important reservoirs of microorganisms on the planet (Lindow and Brandl, 2003; Monier and Lindow, 2004). On a global scale, it has been estimated that phyllosphere bacteria could comprise 1026 cells, with microbial density up to 106–107 cells per cm2 (Morris et al., 2002; Lindow and Brandl, 2003; Monier and Lindow, 2004). Phyllosphere bacteria are a component of the plant microbiome and have an important role in facilitating plant growth, protecting crops from external pathogens (Rasche et al., 2006), and mediating carbon and nitrogen cycles (Furnkranz et al., 2008; Redford and Fierer, 2009). Potential immigrants to the phyllosphere such as Escherichia coli pose a risk to human health through the transfer of human-pathogenic bacteria from ready-to-eat food (Teplitski et al., 2011; Chen et al., 2019; Zhou et al., 2020). Selection pressure is essential for structuring communities in the phyllosphere (Yang et al., 2001; Kim et al., 2012) and is typically strong as oligotrophic conditions provide low levels of nutrients in most part of the leaves (Lindow and Brandl, 2003). Additionally, the dynamic external environment such as fluctuations in UV light irradiation, humidity, and temperature is an additional pressure on phyllosphere microorganisms (Hirano and Upper, 2000; Kadivar and Stapleton, 2003; Redford and Fierer, 2009). Moreover, plant genotypes may contribute to the selection of phyllosphere microorganisms (Bodenhausen et al., 2014; Zarraonaindia et al., 2015; Agler et al., 2016).
The diversity of microorganisms within the phyllosphere is greater than previously thought but is lower than those in either rhizosphere or bulk soil (Delmotte et al., 2009; Knief et al., 2012). The phyllosphere microbial community does not exist in isolation and can recruit members from rainfall, irrigation, atmospheric deposition, and transfer from soils and other plant organs (Lopez-Velasco et al., 2013; Vacher et al., 2016). Soil has been widely considered a microbial reservoir for the plant microbiome (Vorholt, 2012; Gopal and Gupta, 2016). Emerging evidence suggests that soil microbial communities are important sources for the plant phyllosphere. For example, shared taxa have been found between the soil and the phyllosphere of grape, switchgrass, and perennial mustard, suggesting that soil is a potential contributor of phyllosphere microbes (Martins et al., 2013; Zarraonaindia et al., 2015; Grady et al., 2019). Furthermore, airborne microorganisms are another important source of phyllosphere microorganisms (Yan et al., 2020). Phyllosphere microbial communities also have the potential to affect the composition of the microbiome in surrounding air (Wei et al., 2017). In addition, the microbial composition of the phyllosphere could also be determined from seed hereditary (Kumar et al., 2005; Vacher et al., 2016). Furthermore, human intervention, such as irrigation and fertilization, could also be a driver of phyllosphere microbial communities (Zhu et al., 2017; Gekenidis et al., 2018). Meanwhile, human activities have the potential to introduce antibiotic-resistance genes and antibiotic-resistant pathogens and bacteria (ARB) into the plant system, posing a potential risk to human health (Chen et al., 2019).
While it is recognized that there are multiple sources to the phyllosphere microbial community including pathogenic bacteria or ARB, the contribution of each source is unknown (Vorholt, 2012). While previous studies have focused on the source-sink relationships between soil and phyllosphere or between air and phyllosphere (Lindemann and Upper, 1985; Martins et al., 2013), a more comprehensive evaluation is needed to understand the movement of bacteria within three sources (air, phyllosphere, and soil).
To determine the possible sources of phyllosphere bacteria, we designed a controlled system for plant cultivation that would prevent the recruitment of bacteria from external uncontrolled sources such as irrigation and aerosol deposition. By using a controlled environment, the external airborne microbes could be recruited as a potential source for the phyllosphere and at the same time facilitate evaluation of the contribution of possible other sources (air or soilborne microbes) to the phyllosphere. The second chamber used filtered air, enabling the evaluation of the contribution of soil solely. Thus, in this study we could characterize (1) the phyllosphere microbial community of two plant species and the impact of external airborne microbes to those phyllosphere communities and (2) determine the possible sources (air or soil-borne microbes) of the phyllosphere microbial communities.
Materials and Methods
Microcosm Design
A connected two-chamber device made of Acrylic plate was constructed for plant cultivation (Figure 1). One chamber permitted external air to flow through holes in one wall of the chamber, while the other chamber received filtered air from the previous chamber. A total of six filters (0.22 μm) with a diameter of 80 mm, located at the intersection of the two chambers and the back wall of chamber 2, were used for collection and filtering of airborne microbes. As chamber 1 had access to external air flow, it was considered for the purpose of this study as “microbial air” (M). As the air of chamber 2 was filtered, it was considered to be “non-microbial air” (NM).
Figure 1. Schematic of the equipment used in the current study for plant cultivation and identifying the source of microbial communities.
A vacuum pump (flow rate was 200 L min–1) was used to induce the air flow within both chambers. To facilitate manipulation of plants while maintaining a sealed environment, sterile gloves were built into the chamber sides (Figure 1). Sterile water was used for irrigation to avoid the introduction of external microbes and entered through two drilled holes along each long chamber side. Built-in pots could be arranged in the culturing system. Except for the air holes in chamber 1, the system was hermetically sealed for the duration of the experiment. The specific design parameters are listed in supplementary materials (Supplementary Figure 1).
Experimental Design and Sampling
Soil was collected from a vegetable field in Xiamen city, China (24°38′26.9″N 118°02′28.8″E) and sieved to 5 mm prior to use. Allium schoenoprasum and Sonchus oleraceus were chosen to evaluate the phyllosphere microbial community as they can be ingested raw as either a garnish or ready-to-eat salad and represent a major pathway for human health impact by phyllosphere microbiomes. Plants were cultivated in both chambers from seed and irrigated with sterile water. Plant pots and the surface of seeds were sterilized before cultivation. Plants were harvested after 30 days, and DNA was extracted from collected leaves. Airborne microbes were captured on filters (n = 6) for DNA extraction. Sample codes were the following: soils for cultivation of A. schoenoprasum and S. oleraceus were SCM and SSOM (in chamber 1), SCNM and SSONM (in chamber 2), respectively; the phyllosphere samples PCM and PSOM (in chamber 1) and PCNM and PSONM (in chamber 2); AM and ANM represented air samples with or without outdoor airborne microbes. Destructive sampling was conducted to collect the soil samples. Except for air samples (n = 3), there were four replicates for each sample type. All samples were stored at −20°C for DNA extraction.
DNA Extraction From Soil, Phyllosphere, and Air Samples
The collection of microbes from leaves and subsequent DNA extraction was conducted as previously reported (Zhou et al., 2019). Collected soil samples were contained in a sterilized bag, and 500 mg soil was used for DNA extraction. Filters that captured airborne microbes from chambers 1 and 2 were cut into small pieces by sterile scissors prior to DNA extraction. DNA extractions of air, phyllosphere, and soil samples used a FastDNA Spin Kit for Soil (MP Bio, United States) following the manufacturer’s instructions. The quality of the extracted DNA was evaluated by a Qubit 3.0 Fluorometer (Invitrogen, Ghent, Belgium) (Zheng et al., 2020). Extracted DNA was stored at −20°C until used.
Amplicon Sequencing of 16S RNA Genes and Analysis
The 515F/907R primer set was used to amplify the V4–V5 region of the 16S rRNA gene from air, phyllosphere, and soil samples (Turner et al., 1999). PCR protocols and conditions were as previously reported (Chen et al., 2017). To distinguish samples, unique barcodes were used for each sample. Prior to sequencing them on the Illumina 2,500 platform (Novogene, Beijing, China), the concentration of PCR products was quantified using a Qubit 3.0 Fluorometer. Quantitative Insight Into Microbiology Ecology (QIIME, version 1.9.1) was used for sequence analysis (Caporaso et al., 2010). Operational taxonomic units (OTU) were determined by UCLUST clustering (Edgar, 2010) with the similarity set at 97%. OTUs with only a single sequence were discarded from the final OTU table (Caporaso et al., 2010). Taxonomic classification of OTUs was performed using the RDP classifier with the Greengenes database version 13.5 (McDonald et al., 2012). Chao1, Observed species, PD whole tree, and Shannon diversity were calculated using QIIME. Sequences generated in this study were submitted to the National Center for Biotechnology Information (NCBI), with accession number PRJNA643678.
Statistical Analysis
Excel 2016 was used for mathematical calculations of the raw data (for example, means, standard errors, and sum). Analysis of variance was conducted in SPSS 21, and p < 0.05 was considered statistically significant. Principal Coordinate Analysis (PCoA) based on Bray-Curtis distances and permutational multivariate analysis of variance (PERMANOVA) were performed using vegan package (Dixon, 2003; Oksanen et al., 2019) and visualized using ggplot2 package (version 3.1) (Wickham et al., 2020) in R. The microbial source was tracked using fast expectation-maximization microbial source tracking (FEAST), following protocols provided by the authors of the R package (Shenhav et al., 2019). Bar charts were created in OriginPro 2018.
Results
Composition of Microbial Communities
A total of 3,947,497 high-quality sequences, which ranged from 48,053 (air) to 245,792 (soil) were recorded across all samples. Proteobacteria (25.5%), Chloroflexi (10.6%), Acidobacteria (9.3%), Firmicutes (7.7%), and Nitrospirae (7.5%) made up more than 65% of airborne microbial OTUs (Figure 2A). In contrast, Proteobacteria (74.8%) dominated phyllosphere samples, followed by Firmicutes (12.4%), Bacteroidetes (4.7%), and Actinobacteria (4.3%). Compared with both air and phyllosphere microbial communities, soil microbial communities were relatively homogeneous with the similar relative abundance of Proteobacteria, Acidobacteria, Firmicutes, and Crenarchaeota across all soil samples. The relative abundance of Proteobacteria in the phyllosphere was significantly greater than that in either air or soil samples (p < 0.01, ANOVA). Accordingly, the relative abundance of Acidobacteria, Chloroflexi, Planctomycetes, Nitrospirae, and Crenarchaeota in the phyllosphere were significantly lower than those in air and soil samples (p < 0.05, ANOVA). ANOVA results at phylum level are presented in the supplementary materials (Supplementary Tables 1, 2).
Figure 2. The percentage of microbial communities at phylum (A) and family (B) level, calculated by the average of replicates. “M” represents the chamber 1 samples (with extra outdoor microbes) and “NM” represents chamber 2 samples (without extra microbes). AM, PCM, PSOM, SCM, and SSOM represents the air, phyllosphere, and soil of Allium schoenoprasum and Sonchus oleraceus samples, respectively, in chamber 1 (with external airborne microbes). The ANM, PCNM, PSONM, SCNM, and SSONM represents the air, phyllosphere, and soil of A. schoenoprasum and S. oleraceus samples, respectively, in chamber 2 (without external airborne microbes).
The dominant families with relative abundance >1% in were Nitrososphaeraceae (7.3%), Bacillaceae (3.6%), Syntrophobacteraceae (3.1%), Pirellulaceae (2.6%), and Hyphomicrobiaceae (1.8%) (Figure 2B) in air samples. Rhizobiaceae (38.2%), Pseudomonadaceae (15.8%), and Bacillaceae (9.0%) represented 63.3% of the total microbial community in the phyllosphere. The abundance of Rhizobiaceae and Pseudomonadaceae was significantly greater in phyllosphere (p < 0.05, ANOVA) than either air or soil samples (Supplementary Tables 3, 4).
Comparative Analysis of Air, Phyllosphere and Soil Microbial Communities
The diversity of microbial communities in the phyllosphere was significantly lower than those of air and soil samples (p < 0.05, ANOVA, Supplementary Figure 2). Average linkage clustering of OTU level data (Figure 3A) revealed a clear separation between air, phyllosphere, and soil samples (p = 0.001, ANOVA). The different plant species, Sonchus oleraceus and A. schoenoprasum, separated into different clusters (p = 0.001, PERMANOVA). The phyllosphere samples of A. schoenoprasum from different chambers (with or without external air microbiota) were separated into two clusters; however, the phyllosphere samples of S. oleraceus were clustered together (Figure 3A).
Figure 3. Dissimilarity analysis of microbial OTUs, phyla, and families for all samples. “M” represents chamber 1 samples (with extra outdoor microbes) and “NM” represents chamber 2 samples (without extra microbes). (A) Cluster diagram at OTU level. (B,C) The Principal Coordinate Analysis (PCoA) is based on the Bray-Curtis distance for microbial phylum and family, respectively. Different colors represent different habitats, while different shapes indicate each chamber.
At the phylum level, the composition of microbes from air, phyllosphere, and soil samples clearly separated along PCo1 (p < 0.001, PERMANOVA, Figure 3B), which explained 84.7% of the total variance. At the family level, a similarly significant separation existed along PCo1, explaining 67.0% of the total variance (Figure 3C). The phyllosphere microbial communities were significantly different between the two plant species in both chambers 1 and 2 (p < 0.001, PERMANOVA, Supplementary Figures 3A,B). Within the same plant species, the A. schoenoprasum phyllosphere microbial communities were separated between chamber 1 and 2 along the PCo1 which explained 42% of the variance (p < 0.05, PERMANOVA, Supplementary Figure 3C). In contrast, no significant differences were observed within S. oleraceus samples (p > 0.05, PERMANOVA, Supplementary Figure 3D).
The Shared Microbiota and Source Tracking
The shared microbiota among air, phyllosphere, and soil samples in each chamber were analyzed at the OTU level. In chamber 1, air (AM), A. schoenoprasum phyllosphere (PCM) and A. schoenoprasum soil (SCM) samples shared 1099 OTUs. AM and SCM shared 3130 OTUs, which was greater than those between either AM and PCM (1265) or PCM and SCM (1199). A total of 115 and 90 unique microbial OTUs were found in the phyllosphere of A. schoenoprasum (PCM) and S. oleraceus (PSOM), respectively (Figures 4A,C). A similar pattern of shared OTUs was observed in chamber 2 (Figures 4B,D). 391 (chamber 1) and 294 (chamber 2) unique OTUs were found in air samples, which represents 10.7% and 8.1% of the total OTUs, respectively (Figure 4E). There were 2964 OTUs were shared between AM and ANM, which accounted for 88.3% and 91.0% of the total OTUs in each chamber. The PCN and PCNM, PSOM and PSONM shared a total 859 and 278 OTUs in chambers 1 and 2, respectively (Figures 4F,G).
Figure 4. Venn diagram at OTU level for the two chambers. “M” represents chamber 1 samples (with extra outdoor microbes) and “NM” represents chamber 2 samples (without extra microbes). AM and ANM, PCN and PCNM, PSOM and PSONM represent the air, phyllosphere, and soil samples of Allium schoenoprasum and Sonchus oleraceus in chambers 1 and 2 (with or without extra airborne microbes), respectively. (A,C) represent the overview of shared OTUs between air, phyllosphere, and soil samples in chamber 1; (B,D) represent shared OTUs between air, phyllosphere, and soil samples in chamber 2. (E–G) represent shared microbial OTUs between AM and ANM, PCN and PCNM, PSOM and PSONM, respectively.
A source tracking method, FEAST, was used for tracking the origin of the phyllosphere microbiota based on the OTU data. Contributions from each source to sink were calculated and represented as a percentage. The airborne microbial community had approximately 19.5% of the community sourced from soil and 3% (S. oleraceus) and 4% (A. schoenoprasum) from the respective phyllosphere (Figure 5A). A total of 5.0% of the microbiota in A. schoenoprasum phyllosphere originated from the air (2.3%) and A. schoenoprasum soil (2.7%). Soil associated with S. oleraceus and the air each accounted for <1% of the S. oleraceus phyllosphere communities. In general, the relative abundance of each microbial source to each of the three (air, phyllosphere, and soil) microbial communities in chamber 2 was similar to that for chamber 1 (Figure 5B).
Figure 5. Fast expectation-maximization microbial source tracking (FEAST) analysis for chamber 1 (A) and chamber 2 (B) based on OTUs level. Direction of the arrows represents the source-sink relationships, and percentages represent the contribution that each source provides.
Discussion
The Composition of Microbes
In this study, we created a purpose-built system that controlled air movement for comparison of microbial communities between air, phyllosphere, and soil samples. Air from the external environment entered chamber 1 and then chamber 2 after filtration. Results demonstrated that the composition of the microbial communities was significantly different in the three compartments (p < 0.05, PERMANOVA), which is consistent with previous studies (Vokou et al., 2012; Bulgarelli et al., 2013). At the phylum level, the archaea Crenarchaeota was sharply lower (<1%) in the phyllosphere but higher in both air and soil samples, which concurred with the previous study (Vorholt, 2012). The shared of Crenarchaeota within air, phyllosphere, and soil suggested that soil may be a possible source of archaea affecting adjacent air (Wehking et al., 2018). Previous studies have found that the bacterial families of Rhizobiaceae and Pseudomonadaceae are commonly observed in plants, which would help create differences in composition of the phyllosphere, air, and soil samples (Yanni et al., 1997; Hunter et al., 2010). A significantly higher abundance (p < 0.001, ANOVA) of Rhizobiaceae in the phyllosphere than the other habitats, which indicated that the phyllosphere may be the ideal habitat for surviving of Rhizobiaceae. It has been reported that the Rhizobiaceae can migrate from root to the above-ground plant parts through endophytic system in rice or tobacco, which suggested that phyllosphere may acquire certain microbes from soil and finally to the leaves (Chi et al., 2005; Ji et al., 2010).
The composition of soil microbial communities between the two chambers was similar, and there was no difference in the soil microbes associated with the two plant species used in this study. Although the same soil was used for growing both plant species, the phyllosphere microbial community composition was clearly distinct for each of the two plant species in both chambers. For instance, the phyla of Acidobacteria, Crenarchaeota, and Chloroflexi were significantly higher in PCM than PSOM in chamber 1. Hyphomicrobiaceae, Pirellulaceae, and Syntrophobacteraceae were significantly higher in PCNM than PSONM in chamber 2, which suggested plant species may therefore be one of the drivers to shape the microbial community in the phyllosphere (Whipps et al., 2008; Redford et al., 2010).
The differentiated composition of dominant taxa between the phyllosphere and other habitats suggested that rare species from the soil and air may colonize and grow on plant leaves then become dominant taxa in phyllosphere communities. The depth of the sequencing may be insufficient to effectively characterize rare species in microbial communities in air, phyllosphere, and soil sample.
Source of the Phyllosphere Microbial Community
As environmental factors were controlled in a sealed system in this study, FEAST analysis has the potential to identify the sources of the microbial communities in the air and phyllosphere. Unlike with previous studies, in this system we can precisely identify the source of the microbial communities in the air, phyllosphere, and soil samples (Zarraonaindia et al., 2015; Chen et al., 2017). In addition, the microcosm system allowed us to evaluate the contribution of soil solely by excluding the airborne microbes in chamber 2.
To understand the interactions between the air, phyllosphere, and soil samples, shared OTUs were investigated. The number of shared microbial OTUs between air and soil as well as between phyllosphere and soil samples was considerable, which suggested a potential interconnection of microbes from these sources. The number of unique OTUs in AM (391) and ANM (294) were fewer than those shared between AM and ANM (2,964). This difference in OTU number between air samples may have been affected by the microbial composition of soil and phyllosphere sources in each chamber. As the air was continuously in motion, the communities in the soil and phyllosphere may have been distributed into the air within each chamber, as perhaps highlighted with the noted contribution of phyllosphere and soil microbes to the air microbial communities (Lindemann and Upper, 1985; Lymperopoulou et al., 2016). Soil and phyllosphere microbes were the main sources of airborne microbial communities in both chambers. It is likely that the air circulation in the chamber may have resulted in the suspension of soil and phyllosphere microbes in the air, allowing them to become part of the airborne microbial community (Bowers et al., 2013; Smets et al., 2016). To judge from this system, the contribution of the phyllosphere to the airborne microbial community was greater than that from the air to the phyllosphere, which concurred with Lymperopoulou et al. (2016).
Soil can affect the microbial communities in the phyllosphere either directly or indirectly. It has been reported that microbes originating from the soil can colonize plant roots and then reach the phyllosphere (Beattie and Lindow, 1999; Bodenhausen et al., 2013) through endophytic transport (Chi et al., 2005), which would directly affect the phyllosphere microbiome. Previous studies had been mainly focused on this direct effect. However, there are few studies that consider the indirect effect of the soil microbial communities on the phyllosphere, such as through airborne distribution. In the system, the composition of airborne microbes in chamber 2 was predominately acquired from soil microbes through air currents with no external microbes input from environment, which implied a potential pathway for microbes transportation from soil to air and finally to phyllosphere (Lindemann and Upper, 1985; Martins et al., 2013; Grady et al., 2019). Thereby, soil microbes may indirectly contribute to phyllosphere microbes. In this study, the contributions of the soil and air to the phyllosphere microbial communities were low, whereas a considerable proportion of unknown sources contributed to the phyllosphere microbial communities according to FEAST analysis. Although there was a shared component of the microbial community in the soil, air, and phyllosphere, none of the shared microbes were dominant microbes, which suggests that plant microbes may be acquired inherently from seed (Wulff et al., 2003; Vorholt, 2012). The results from a recent study (preprint) supports our point (Abdelfattah et al., 2020). Thus, inheritance of microbes first from seeds and parent material and supplemented from the surrounding environment, such as the soil and air, may be the potential pathway to explain the origination of phyllosphere microbes (Mukhopadhyay et al., 1996; Truyens et al., 2015; Cope-Selby et al., 2017). Future studies including the seed microbiome with a time-series sampling strategy are needed for more effective characterization of community change in each habitat.
In summary, the identified shared taxa highlighted the potential exchange of microbes within the air-phyllosphere-soil continuum. Part of the airborne microbes originated from the soil and phyllosphere in both chambers. However, source tracking analysis indicated the soil and air may not be the major sources of the A. schoenoprasum and S. oleraceus phyllosphere microbial communities, although we found considerable taxa overlap within three habitats.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
S-Y-DZ contributed to the conception, analysis, depict, and drafted the manuscript. HL contributed to the conception, design of the microcosm, experiment, and review of the manuscript. MG, RN, X-rY, and J-qS contributed to the review of the manuscript. All authors gave the final approval and agreed to be accountable for all aspects of the work.
Funding
This work was funded by the National Natural Science Foundation of China (41977201 and 41701280), and The James Hutton Institute receives financial support from Scottish Government Rural and Environment Science and Analytical Services (RESAS).
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.
The reviewer DJ declared a shared affiliation with several of the authors S-Y-DZ, HL, and X-rY to the handling editor at the time of review.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.615481/full#supplementary-material
Supplementary Figure 1 | The specific design parameters of the equipment.
Supplementary Figure 2 | Chao1, Observed species, PD whole tree, and Shannon analysis of Alpha diversity.
Supplementary Figure 3 | Principal Coordinate Analysis (PCoA) based on Bray-Curtis distances of microbial OTU.
References
Abdelfattah, A., Wisniewski, M., Schena, L., and Tack, A. J. (2020). Experimental Evidence of Microbial Inheritance in Plants and Transmission Routes from seed to Phyllosphere and Root. Durham, NC: Research Square, doi: 10.21203/rs.3.rs-27656/v1
Agler, M. T., Ruhe, J., Kroll, S., Morhenn, C., Kim, S. T., Weigel, D., et al. (2016). Microbial hub taxa link host and abiotic factors to plant microbiome Variation. PLoS Biol. 14:e1002352. doi: 10.1371/journal.pbio.1002352
Beattie, G. A., and Lindow, S. E. (1999). Bacterial colonization of leaves: a spectrum of strategies. Phytopathology 89, 353–359. doi: 10.1094/phyto.1999.89.5.353
Bodenhausen, N., Bortfeld-Miller, M., Ackermann, M., and Vorholt, J. A. (2014). A synthetic community approach reveals plant genotypes affecting the phyllosphere microbiota. PLoS Genet. 10:e1004283. doi: 10.1371/journal.pgen.1004283
Bodenhausen, N., Horton, M. W., and Bergelson, J. (2013). Bacterial communities associated with the leaves and the roots of Arabidopsis thaliana. PLoS One 8:e56329. doi: 10.1371/journal.pone.0056329
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
Bulgarelli, D., Schlaeppi, K., Spaepen, S., van Themaat, E. V. L., and Schulze-Lefert, P. (2013). Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 64, 807–838. doi: 10.1146/annurev-arplant-050312-120106
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
Chen, Q. L., An, X. L., Zhu, Y. G., Su, J. Q., Gillings, M. R., Ye, Z. L., et al. (2017). Application of struvite alters the antibiotic resistome in soil, rhizosphere, and phyllosphere. Environ. Sci. Technol. 51, 8149–8157. doi: 10.1021/acs.est.7b01420
Chen, Q.-L., Cui, H.-L., Su, J.-Q., Penuelas, J., and Zhu, Y.-G. (2019). Antibiotic resistomes in plant microbiomes. Trends Plant Sci. 24, 530–541. doi: 10.1016/j.tplants.2019.02.010
Chi, F., Shen, S. H., Cheng, H. P., Jing, Y. X., Yanni, Y. G., and Dazzo, F. B. (2005). Ascending migration of endophytic rhizobia, from roots to leaves, inside rice plants and assessment of benefits to rice growth physiology. Appl. Environ. Microbiol. 71, 7271–7278. doi: 10.1128/aem.71.11.7271-7278.2005
Cope-Selby, N., Cookson, A., Squance, M., Donnison, I., Flavell, R., and Farrar, K. (2017). Endophytic bacteria in Miscanthus seed: implications for germination, vertical inheritance of endophytes, plant evolution and breeding. Glob. Chang. Biol. Bioenergy 9, 57–77. doi: 10.1111/gcbb.12364
Delmotte, N., Knief, C., Chaffron, S., Innerebner, G., Roschitzki, B., Schlapbach, R., et al. (2009). Community proteogenomics reveals insights into the physiology of phyllosphere bacteria. Proc. Natl. Acad. Sci. U.S.A. 106, 16428–16433. doi: 10.1073/pnas.0905240106
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Furnkranz, M., Wanek, W., Richter, A., Abell, G., Rasche, F., and Sessitsch, A. (2008). Nitrogen fixation by phyllosphere bacteria associated with higher plants and their colonizing epiphytes of a tropical lowland rainforest of Costa Rica. ISME J. 2, 561–570. doi: 10.1038/ismej.2008.14
Gekenidis, M. T., Schoner, U., von Ah, U., Schmelcher, M., Walsh, F., and Drissner, D. (2018). Tracing back multidrug-resistant bacteria in fresh herb production: from chive to source through the irrigation water chain. FEMS Microbiol. Ecol. 94:fiy149. doi: 10.1093/femsec/fiy149
Gopal, M., and Gupta, A. (2016). Microbiome selection could spur next-generation plant breeding strategies. Front. Microbiol. 7:1971. doi: 10.3389/fmicb.2016.01971
Grady, K. L., Sorensen, J. W., Stopnisek, N., Guittar, J., and Shade, A. (2019). Assembly and seasonality of core phyllosphere microbiota on perennial biofuel crops. Nat. Commun. 10:4135. doi: 10.1038/s41467-019-11974-4
Hirano, S. S., and Upper, C. D. (2000). Bacteria in the leaf ecosystem with emphasis on Pseudomonas syringae - a pathogen, ice nucleus, and epiphyte. Microbiol. Mol. Biol. Rev. 64, 624–653. doi: 10.1128/mmbr.64.3.624-653.2000
Hunter, P. J., Hand, P., Pink, D., Whipps, J. M., and Bending, G. D. (2010). Both leaf properties and microbe-microbe interactions influence within-species variation in bacterial population diversity and structure in the Lettuce (Lactuca Species) Phyllosphere. Appl. Environ. Microbiol. 76, 8117–8125. doi: 10.1128/aem.01321-10
Ji, K.-X., Chi, F., Yang, M.-F., Shen, S.-H., Jing, Y.-X., Dazzo, F. B., et al. (2010). Movement of Rhizobia Inside Tobacco and lifestyle alternation from endophytes to free-living Rhizobia on leaves. J. Microbiol. Biotechnol. 20, 238–244. doi: 10.4014/jmb.0906.06042
Kadivar, H., and Stapleton, A. E. (2003). Ultraviolet radiation alters maize phyllosphere bacterial diversity. Microb. Ecol. 45, 353–361. doi: 10.1007/s00248-002-1065-5
Kim, M., Singh, D., Lai-Hoe, A., Go, R., Rahim, R. A., Ainuddin, A. N., et al. (2012). Distinctive phyllosphere bacterial communities in tropical trees. Microb. Ecol. 63, 674–681. doi: 10.1007/s00248-011-9953-1
Knief, C., Delmotte, N., Chaffron, S., Stark, M., Innerebner, G., Wassmann, R., et al. (2012). Metaproteogenomic analysis of microbial communities in the phyllosphere and rhizosphere of rice. ISME J. 6, 1378–1390. doi: 10.1038/ismej.2011.192
Kumar, K., Gupta, S. C., Baidoo, S. K., Chander, Y., and Rosen, C. J. (2005). Antibiotic uptake by plants from soil fertilized with animal manure. J. Environ. Qual. 34, 2082–2085. doi: 10.2134/jeq2005.0026
Lindemann, J., and Upper, C. D. (1985). Aerial dispersal of epiphytic bacteria over bean-plants. Appl. Environ. Microbiol. 50, 1229–1232. doi: 10.1128/aem.50.5.1229-1232.1985
Lindow, S. E., and Brandl, M. T. (2003). Microbiology of the phyllosphere. Appl. Environ. Microbiol. 69, 1875–1883. doi: 10.1128/aem.69.4.1875-1883.2003
Lopez-Velasco, G., Carder, P. A., Welbaum, G. E., and Ponder, M. A. (2013). Diversity of the spinach (Spinacia oleracea) spermosphere and phyllosphere bacterial communities. FEMS Microbiol. Lett. 346, 146–154. doi: 10.1111/1574-6968.12216
Lymperopoulou, D. S., Adams, R. I., and Lindow, S. E. (2016). Contribution of vegetation to the microbial composition of nearby outdoor Air. Appl. Environ. Microbiol. 82, 3822–3833. doi: 10.1128/aem.00610-16
Martins, G., Lauga, B., Miot-Sertier, C., Mercier, A., Lonvaud, A., Soulas, M. L., et al. (2013). Characterization of epiphytic bacterial communities from grapes, leaves, bark and soil of grapevine plants grown, and their relations. PLoS One 8:e73013. doi: 10.1371/journal.pone.0073013
McDonald, D., Price, M. N., Goodrich, J., Nawrocki, E. P., DeSantis, T. Z., Probst, A., et al. (2012). An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618. doi: 10.1038/ismej.2011.139
Monier, J. M., and Lindow, S. E. (2004). Frequency, size, and localization of bacterial aggregates on bean leaf surfaces. Appl. Environ. Microbiol. 70, 346–355. doi: 10.1128/aem.70.1.346-355.2004
Morris, C. E., Barnes, M. B., and McLean, R. J. C. (2002). Biofilms on Leaf Surfaces: Implications for the Biology, Ecology and Management of Populations of Epiphytic Bacteria. St Paul: Amer Phytopathological Soc.
Mukhopadhyay, K., Garrison, N. K., Hinton, D. M., Bacon, C. W., Khush, G. S., Peck, H. D., et al. (1996). Identification and characterization of bacterial endophytes of rice. Mycopathologia 134, 151–159. doi: 10.1007/bf00436723
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., Mcglinn, D., et al. (2019). vegan: Community Ecology Package. Available online at: https://cran.r-project.org/web/packages/vegan/index.html (accessed September 1, 2019).
Rasche, F., Marco-Noales, E., Velvis, H., van Overbeek, L. S., Lopez, M. M., van Elsas, J. D., et al. (2006). Structural characteristics and plant-beneficial effects of bacteria colonizing the shoots of field grown conventional and genetically modified T4-lysozyme producing potatoes. Plant Soil 289, 123–140. doi: 10.1007/s11104-006-9103-6
Redford, A. J., Bowers, R. M., Knight, R., Linhart, Y., and Fierer, N. (2010). The ecology of the phyllosphere: geographic and phylogenetic variability in the distribution of bacteria on tree leaves. Environ. Microbiol. 12, 2885–2893. doi: 10.1111/j.1462-2920.2010.02258.x
Redford, A. J., and Fierer, N. (2009). Bacterial succession on the leaf surface: a novel system for studying successional dynamics. Microb. Ecol. 58, 189–198. doi: 10.1007/s00248-009-9495-y
Shenhav, L., Thompson, M., Joseph, T. A., Briscoe, L., Furman, O., Bogumil, D., et al. (2019). FEAST: fast expectation-maximization for microbial source tracking. Nat. Methods 16, 627–632. doi: 10.1038/s41592-019-0431-x
Smets, W., Wuyts, K., Oerlemans, E., Wuyts, S., Denys, S., Samson, R., et al. (2016). Impact of urban land use on the bacterial phyllosphere of ivy (Hedera sp.). Atmos. Environ. 147, 376–383. doi: 10.1016/j.atmosenv.2016.10.017
Teplitski, M., Warriner, K., Bartz, J., and Schneider, K. R. (2011). Untangling metabolic and communication networks: interactions of enterics with phytobacteria and their implications in produce safety. Trends Microbiol. 19, 121–127. doi: 10.1016/j.tim.2010.11.007
Truyens, S., Weyens, N., Cuypers, A., and Vangronsveld, J. (2015). Bacterial seed endophytes: genera, vertical transmission and interaction with plants. Environ. Microbiol. Rep. 7, 40–50. doi: 10.1111/1758-2229.12181
Turner, S., Pryer, K. M., Miao, V. P. W., and Palmer, J. D. (1999). Investigating deep phylogenetic relationships among Cyanobacteria and plastids by small submit rRNA sequence analysis. J. Eukaryot. Microbiol. 46, 327–338. doi: 10.1111/j.1550-7408.1999.tb04612.x
Vacher, C., Hampe, A., Porte, A. J., Sauer, U., Compant, S., and Morris, C. E. (2016). “The phyllosphere: microbial jungle at the plant-climate interface,” in Annual Review of Ecology, Evolution, and Systematics, Vol. 47, ed. D. J. Futuyma (Palo Alto: Annual Reviews), 1–24. doi: 10.1146/annurev-ecolsys-121415-032238
Vokou, D., Vareli, K., Zarali, E., Karamanoli, K., Constantinidou, H. I. A., Monokrousos, N., et al. (2012). Exploring biodiversity in the bacterial community of the mediterranean phyllosphere and its relationship with airborne bacteria. Microb. Ecol. 64, 714–724. doi: 10.1007/s00248-012-0053-7
Vorholt, J. A. (2012). Microbial life in the phyllosphere. Nat. Rev. Microbiol. 10, 828–840. doi: 10.1038/nrmicro2910
Wehking, J., Pickersgill, D. A., Bowers, R. M., Teschner, D., Poschl, U., Frohlich-Nowoisky, J., et al. (2018). Community composition and seasonal changes of archaea in coarse and fine air particulate matter. Biogeosciences 15, 4205–4214. doi: 10.5194/bg-15-4205-2018
Wei, X. Y., Lyu, S. H., Yu, Y., Wang, Z. H., Liu, H., Pan, D. M., et al. (2017). Phylloremediation of air pollutants: exploiting the potential of plant leaves and leaf-associated microbes. Front. Plant Sci. 8:1318. doi: 10.3389/fpls.2017.01318
Whipps, J. M., Hand, P., Pink, D., and Bending, G. D. (2008). Phyllosphere microbiology with special reference to diversity and plant genotype. J. Appl. Microbiol. 105, 1744–1755. doi: 10.1111/j.1365-2672.2008.03906.x
Wickham, H., Chang, W., Henry, L., Pedersen, L., Takahashi, K., Wilke, C., et al. (2020). ggplot2: Create Elegant Data Visualisations using the Grammar of Graphics. Available online at: https://cran.r-project.org/web/packages/ggplot2/ (accessed June 19, 2020).
Wulff, E. G., van Vuurde, J. W. L., and Hockenhull, J. (2003). The ability of the biological control agent Bacillus subtilis, strain BB, to colonise vegetable brassicas endophytically following seed inoculation. Plant Soil 255, 463–474. doi: 10.1023/a:1026026917355
Yan, Z. Z., Chen, Q. L., Zhang, Y. J., He, J. Z., and Hu, H. W. (2020). Industrial development as a key factor explaining variances in soil and grass phyllosphere microbiomes in urban green spaces. Environ. Pollut. 261:114201. doi: 10.1016/j.envpol.2020.114201
Yang, C. H., Crowley, D. E., Borneman, J., and Keen, N. T. (2001). Microbial phyllosphere populations are more complex than previously realized. Proc. Natl. Acad. Sci. U.S.A. 98, 3889–3894. doi: 10.1073/pnas.051633898
Yanni, Y. G., Rizk, R. Y., Corich, V., Squartini, A., Ninke, K., PhilipHollingsworth, S., et al. (1997). Natural endophytic association between Rhizobium leguminosarum bv. trifolii and rice roots and assessment of its potential to promote rice growth. Plant Soil 194, 99–114. doi: 10.1023/a:1004269902246
Zarraonaindia, I., Owens, S. M., Weisenhorn, P., West, K., Hampton-Marcell, J., Lax, S., et al. (2015). The soil microbiome influences gapevine-associated microbiota. mBio 6:e2527-14. doi: 10.1128/mBio.02527-14
Zheng, F., Zhu, D., Chen, Q. L., Bi, Q. F., Yang, X. R., O’Connor, P., et al. (2020). The driving factors of nematode gut microbiota under long-term fertilization. FEMS Microbiol. Ecol. 96:fiaa037. doi: 10.1093/femsec/fiaa037
Zhou, S.-Y.-D., Wei, M.-Y., Giles, M., Neilson, R., Zheng, F., Zhang, Q., et al. (2020). Prevalence of antibiotic resistome in ready-to-eat salad. Front. Public Health 8:92. doi: 10.3389/fpubh.2020.00092
Zhou, S.-Y.-D., Zhu, D., Giles, M., Yang, X.-R., Daniell, T., Neilson, R., et al. (2019). Phyllosphere of staple crops under pig manure fertilization, a reservoir of antibiotic resistance genes. Environ. Pollut. 252, 227–235. doi: 10.1016/j.envpol.2019.05.098
Keywords: phyllosphere, leaf microbiota, source tracking, microcosm, airborne microbial community
Citation: Zhou S-Y-D, Li H, Giles M, Neilson R, Yang X-r and Su J-q (2021) Microbial Flow Within an Air-Phyllosphere-Soil Continuum. Front. Microbiol. 11:615481. doi: 10.3389/fmicb.2020.615481
Received: 09 October 2020; Accepted: 02 December 2020;
Published: 12 January 2021.
Edited by:
Yurong Liu, Huazhong Agricultural University, ChinaReviewed by:
Rachel Adams, University of California, Berkeley, United StatesCuijing Zhang, Shenzhen University, China
Decai Jin, Research Center for Eco-Environmental Sciences (CAS), China
Copyright © 2021 Zhou, Li, Giles, Neilson, Yang and Su. 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: Hu Li, hli@iue.ac.cn
†These authors have contributed equally to this work