Original Research ARTICLE
The Arbuscular Mycorrhizal Fungus Funneliformis mosseae Alters Bacterial Communities in Subtropical Forest Soils during Litter Decomposition
- 1Key laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 2World Agroforestry Centre, East and Central Asia, Kunming, China
- 3Centre of Excellence in Fungal Research, Mae Fah Luang University, Chiang Rai, Thailand
- 4School of Science, Mae Fah Luang University, Chiang Rai, Thailand
- 5Department of Soil Ecology, UFZ-Helmholtz Centre for Environmental Research, Halle, Germany
Bacterial communities and arbuscular mycorrhizal fungi (AMF) co-occur in the soil, however, the interaction between these two groups during litter decomposition remains largely unexplored. In order to investigate the effect of AMF on soil bacterial communities, we designed dual compartment microcosms, where AMF (Funneliformis mosseae) was allowed access (AM) to, or excluded (NM) from, a compartment containing forest soil and litterbags. Soil samples from this compartment were analyzed at 0, 90, 120, 150, and 180 days. For each sample, Illumina sequencing was used to assess any changes in the soil bacterial communities. We found that most of the obtained operational taxonomic units (OTUs) from both treatments belonged to the phylum of Proteobacteria, Acidobacteria, and Actinobacteria. The community composition of bacteria at phylum and class levels was slightly influenced by both time and AMF. In addition, time and AMF significantly affected bacterial genera (e.g., Candidatus Solibacter, Dyella, Phenylobacterium) involved in litter decomposition. Opposite to the bacterial community composition, we found that overall soil bacterial OTU richness and diversity are relatively stable and were not significantly influenced by either time or AMF inoculation. OTU richness at phylum and class levels also showed consistent results with overall bacterial OTU richness. Our study provides new insight into the influence of AMF on soil bacterial communities at the genus level.
In terrestrial ecosystems, litter decomposition plays a crucial role in carbon and nutrient cycling, which consists of complex physical, chemical and biological processes (Berg, 2000). These processes are driven primarily by microorganisms in the soil using various enzymes to degrade different components of the litter (e.g., lignin, cellulose, and hemicellulose) (Lynd et al., 2002). Among these microorganisms, although soil saprophytic fungi are usually considered the major decomposers (Hanackova et al., 2015), soil bacteria also play a major role (Greaves, 1971; Purahong et al., 2015). For example, DeAngelis et al. (2013) found that many fast growing bacteria, including Alpha-, Beta-, and Gamma-proteobacteria, Clostridia and Bacteroidetes were the initial decomposers in tropical forest soils. In addition, Actinobacteria are able to degrade cellulose (Schellenberger et al., 2010; Stursova et al., 2012), and many other bacterial groups have also been found to be engaged in lignin decomposition (Masai et al., 2007; Bugg et al., 2011). A recent study provides evidence that complex litter decay processes are the result of a dynamic cross-kingdom, functional succession between fungi and bacteria (Purahong et al., 2016).
Previous studies on litter decomposition have focused on the role of soil fungal communities, especially saprophytic fungi (Kuramae et al., 2013; Persoh et al., 2013; Treseder et al., 2014; Hanackova et al., 2015). However, despite the fact that bacterial communities and non-saprophytic fungi co-exist in the soil, the interactions between these two groups and their effects on litter decomposition has been largely ignored. Arbuscular mycorrhizal fungi (AMF), from the phylum of Glomeromycota, is the most common mycorrhizal association between soil fungi and terrestrial plants, with about 70% of plants forming mycorrhizal symbioses with AMF (Smith et al., 1997). The essential function of this symbiosis is acquiring nutrients (e.g., N, P) through its external mycelium in exchange for carbohydrate from host plants’ photosynthate (Smith and Smith, 2011). Furthermore, AMF are known to have no saprophytic capacity, but they can affect litter decomposition, both positively (Hodge, 2001; Cheng et al., 2012; Herman et al., 2012; Koller et al., 2013), and negatively (Leifheit et al., 2015). Although the mechanism underlying this effect is unclear, AMF might impact litter decomposition through effects on other saprophytic microorganisms in the soil (Nuccio et al., 2013) or on soil structure (Rillig and Mummey, 2006). Soil bacteria can obtain C from the exudates released by mycorrhizal hyphae or by using the hyphae themselves as substrate (Qin et al., 2014). The plant–AMF–soil bacteria interaction likely occurs in the mycorrhizosphere: as defined by Oswald and Ferchau (1968), this is the surface area of the mycelium and consists of two components, the rhizosphere and hyphosphere, which provide niches for the growth of bacteria (Toljander et al., 2006). Additionally, AM fungi are able to enhance the activity of bacteria associated with soil N availability (for plants, fungi, and bacteria), either using raw organic material as a substrate or by influencing the plant N cycle (Saia et al., 2014). In a microcosm experiment, Hodge et al. (2001) reported that the presence of AMF increased the decomposition of complex organic matter and that AMF were able to influence litter decomposition by altering the activity of hyphosphere bacteria. Furthermore, Herman et al. (2012) also reported that AM fungal inoculation altered the C flow as a result of altering the soil microbial community.
However, using the same experimental design as Herman et al. (2012), Nuccio et al. (2013) found that the presence of AMF resulted in a change of about 10% of the relative abundance of the bacterial community. In addition, Welc et al. (2010) reported that AMF mycelium could suppress the soil bacterial community. Furthermore, many specific bacterial taxa are highly associated with AMF as they colonize and live on the surface of the hyphae (Cruz and Ishii, 2011; Iffis et al., 2014). For example, certain Proteobacteria and Firmicutes taxa have been documented as co-existing with AMF (Bonfante and Anca, 2009; Scheublin et al., 2010; Lecomte et al., 2011).
In this study, a microcosm with two units separated by an air-gap was used to study the interaction between the hyphae of Funneliformis mosseae and soil bacterial communities during litter decomposition. We applied high-throughput sequencing techniques to characterize the soil bacterial community during litter decomposition. Specifically, we aimed to investigate how F. mosseae temporally altered soil bacterial diversity (i) and community structure (ii) during litter decomposition. Based on past studies investigating the interaction between AMF and soil bacteria (Toljander et al., 2006; Herman et al., 2012; Saia et al., 2014), we hypothesize that F. mosseae will significantly alter the bacterial community by stimulating key bacterial groups associated with litter decomposition.
Materials and Methods
Materials Used in the Experiment
(1) AMF Inoculum
Arbuscular mycorrhizal fungi inoculum (F. mosseae) was provided by The Institute of Plant Nutrition and Resources, Beijing Academy of Agriculture and Forestry Sciences (Beijing, China). The inoculum contained approximately 1000 spores per 20 g application and the sterile F. mosseae spores were contained (pre-mixed) in a rock flour material (20 g).
The soil (pH = 4.1) used in the microcosms was collected from a subtropical forest located in southwestern China (N 21° 31′ 42.13″, E 100° 29′ 41.87″). The top 5 cm of soil was collected after first removing the litter layer. The soil was sieved using a 2 mm mesh in order to remove any stones or root material. Soil properties have been published in Gui et al. (2017).
Litterbags were made of 200 μm nylon mesh. The litter comprised dried leaves of Calophyllum polyanthum Wall. ex Choisy, an indigenous tree to Yunnan Province, and one of the dominant species from the forest used to collect the soil. The leaves were collected from nursery grown C. polyanthum saplings, and were oven dried at 65°C to a constant weight. Then the dried leaves (Total N = 1.44%) were cut into small pieces (ca. 5 mm∗5 mm), 2 g of which were put into the litterbags. Note that the purpose of the litterbags was to allow for the monitoring of the impact of AMF on litter decomposition, the results of which have been published separately (Gui et al., 2017).
(4) The Host Plant
Trifolium repens L. cv. Milkanova was selected as the host plant. The seeds of T. repens were sterilized before use.
Experimental Setup and Validation
Our experiments were conducted in an acrylic microcosm unit, the design of which has been described in Gui et al. (2017). Briefly, the microcosm unit consisted of two compartments. The first compartment (Host), which was designed to pot the host plant, was filled with sterilized vermiculite and fine gravel (ca. 0.3 cm diameter), which was evenly mixed in a 1:1 ratio. In order to test how F. mosseae interacts with soils from a tropical forest, the soil was placed in the second compartment (Litter). A litterbag (5 cm∗5 cm) was buried in the soil at a depth of 5 cm so that the soil microbial community could fully interact with the surface of the litter. Even though this burying set-up could not exactly simulate what happens in the field, it did allow for the litter to be exposed to the microbial communities in the soil. The two compartments of the microcosm were separated by a plate, which was drilled with evenly spaced holes (4 mm in diameter) and covered by 20 μm nylon mesh on both sides, the mesh allows the hyphae of F. mosseae to pass through, but not plant roots (Supplementary Figure S1).
The microcosms were divided into two treatments, based on the host plant being inoculated with F. mosseae (AM) or not (NM). For NM control, the same amount of sterilized rock flour (20 g) was added into the host compartment. Four replicates were set for each treatment. Additionally, four time-phase samplings were taken at monthly intervals. Thus, in total, 32 microcosms were set up for the experiment. All the microcosms were randomly placed in a greenhouse with daily temperature ranging from 20 to 25°C. Plants received natural light only and no rainwater.
Plant Growth and Percentage AMF Colonization
The 0.2 g of T. repens seeds were planted in the Host compartment of the microcosm (Supplementary Figure S1). The Host compartment received 10 ml of distilled water twice a week and 10 ml of modified Long Ashton nutrient solution once a week (Hewitt and Bureaux, 1966). The Litter compartment received 10 ml of distilled water once a week in order to maintain the moisture levels. After 2 weeks the N and P concentrations in the Long Ashton solution was diluted to 1/10 of the original concentration (34 mg⋅L-1 NaNO3 + 21.4 mg⋅L-1 NH4Cl, 29.2 mg⋅L-1 NaH2PO4⋅2H2O + 4.7 mg⋅L-1 Na2HPO4⋅12H2O) according to Leigh et al. (2009) and adjusted to pH 7.0 with NaOH.
The percentage of T. repens root length colonized by AMF was determined using fresh root samples. Root pieces (2 cm in length) were washed in distilled water and then rinsed with 10% KOH, which was stained with pen ink according to the methods of Vierheilig et al. (1998). The percentage colonization was calculated by a modified line intersection method (McGonigle et al., 1990).
Soil Sampling and Nutrient Analyses
The original soil was collected from the forest as four subsamples which were bulked into one composite sample and then preserved at -20°C for the later analysis and this original soil was marked as “O.” The first time of sampling from the microcosm was conducted 90 days (T90) after planting, subsequent sampling times were at 120 days (T120), 150 days (T150), and 180 days (T180) after planting. Five grams of the soil around the litterbag was collected and preserved at -20°C for further DNA analysis.
Total soil organic carbon was determined by Dumas combustion (White et al., 1997), and total N using a semi-micro Kjeldahl apparatus (Yuen and Pollard, 1953). Total phosphorus (P) and potassium (K) were measured spectrophotometrically after digesting a mixture of concentrated H2SO4 and H2O2 (Murphy and Riley, 1986). Hydrolysable N was analyzed through a reaction with iron (II) sulfate and sodium hydroxide in a diffusion procedure (Mulvaney and Khan, 2001). Available P and K were determined using ammonium fluoride and ammonium acetate (Hedley et al., 1982). All the chemical analyses were conducted in Yunnan Agriculture Academy, Yunnan Province, China.
DNA Extraction, Illumina Sequencing Analysis, and Data Processing
Total soil genomic DNA was extracted from 2 g of fresh soil using OMEGA Soil DNA kit, following the manufacturer’s instructions and stored at -80°C until PCR amplification. The DNA extracts were used to partially amplify the 16S rDNA genes using barcoded primers 519F and 907R (Biddle et al., 2008) on the Illumina Miseq platform. Before amplification, DNA concentration and purity was monitored on 1% agarose gels. DNA was diluted to 1 ng/μL using sterile water according to the concentration. The PCR was carried out in 50-μl reaction mixtures with the following components: 4 μl (2.5 mM) of deoxynucleoside triphosphates, 2 μl (10 μM) of forward and reverse primers, 2U of Taq DNA polymerase with 0.4 μl (TaKaRa), and 1 μl of template containing approximately 50 ng of genomic community DNA as a template. Thirty-five cycles (95°C for 45 s, 56°C for 45 s, and 72°C for 60 s) were performed with a final extension at 72°C for 7 min. Triplicate reaction mixtures per sample were pooled, purified using the QIA quick PCR Purification kit (QIAGEN, Germany) and quantified using a NanoDrop ND-1000 (Thermo Scientific, United States).
Sequencing libraries were generated using TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States) following the manufacturer’s instructions, and index codes were added. Library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. The library was sequenced on an Illumina HiSeq 2500 platform and 250 bp paired-end reads were generated.
Paired-end sequencing reads were assigned to samples based on their unique barcode and truncated by cutting off the barcode and primer sequence and then merged using FLASH (V1.2.71; Magoc and Salzberg, 2011). Quality filtering on the raw tags (splicing sequences) were performed under specific filtering conditions to obtain the high-quality clean tags (Bokulich et al., 2013) using the QIIME (V1.7.02; Caporaso et al., 2010) quality control process with the default parameters that sequences were quality trimmed (>25 quality score and 200 bp in length), and assigned to soil samples based on unique 5-bp barcodes. Furthermore, UCHIME algorithm (UCHIME Algorithm3; Edgar et al., 2011) was used to compare the tags with reference database (Gold database4) to detect and remove chimera sequences and then to obtain the effective tags according to Haas et al. (2011). Sequences with ≥97% similarity were assigned to the same operational taxonomic units (OTUs). Rare OTUs (less than five sequences), which could potentially have originated from sequencing errors were removed from the dataset. The most abundant sequence from each OTU was then selected as a representative sequence for that OTU. For each representative sequence, the Green Gene Database5 was applied to annotate taxonomic information based on RDP 3 classifier algorithm (Version 2.26). At the end, OTUs abundance information was normalized using a sequence standard corresponding to the sample with the least sequences. Subsequent statistical analysis was performed basing on this output normalized data.
Nucleotide Accession Number
The bacterial 16S rDNA genes Illumina sequencing data are deposited in the Sequence Read Archive (SRA) of National Center for Biotechnology Information (NCBI) under the BioProject number PRJNA353568.
Richness estimator (Chao1) and Shannon diversity index were calculated to test whether AMF inoculation and time phase influenced the bacterial OTU richness and diversity. We calculated these two indices using QIIME (Version 1.7.0) and displayed with R software (Version 2.15.3). One-way ANOVA with a general linear model (GLM) procedure (SPSS 18.0) was applied to determine the significance of different treatments (AM fungal inoculation or non-inoculation and sampling times) and their influence on the composition of the soil bacterial community. The results were expressed as mean values with standard error, and compared using Duncan’s multiple range tests. Statistical significance was determined at P < 0.05.
Two-way permutational multivariate analysis of variance (PERMANOVA) was carried out using the software PAST to investigate the effect of sampling time and AMF inoculation treatment on bacterial community structure (Hammer et al., 2001). Principal Component Analysis (PCA) was performed to get principal coordinates and visualizations of complex multidimensional data. A distance matrix of weighted or unweighted unifrac among samples was transformed to a new set of orthogonal axes, by which the maximum variation factor is demonstrated by the first principal coordinate, the second maximum variation factor by the second principal coordinate, and so on. PCA analysis was displayed by WGCNA package, stat packages and ggplot2 package in R software (Version 2.15.3). Goodness-of-fit statistics (R2) of measured factors fitted to the PCA ordination of the soil bacterial community were calculated using the “envfit” function in the “vegan package,” with P-values based on 999 permutations (Oksanen, 2007). We determined the gradient length of the soil bacterial community using Detrended Correspondence Analysis (DCA) and the results show that this bacterial community has short gradient (less than 2.7 SD) (Ramette, 2007). Thus, linear methods such as PCA are appropriate for analysis of soil bacterial community in this study.
AM Fungal Root Colonization and Hyphae Development
Rates of root colonization by F. mosseae in the AMF inoculation treatment increased overtime from 24.9% at first month to 72% at fourth month and were significantly higher than control for all sampling times (Supplementary Figure S2). The control treatment had low levels of root colonization (ranged from 1.5 to 1.75%) that were not significantly different across sampling dates. Furthermore, the roots of T. repens remained un-nodulated for the duration of the experiment.
As a separate study we determined the soil fungal community profile in the Litter compartment using Illumina sequencing (unpublished data). These results show that the relative abundance of Glomeromycota in AMF inoculation treatment was significantly higher than control from T90 to T150 (Supplementary Figure S3). We present this data here as evidence of F. mosseae activity in the Litter compartment.
Soil Bacterial Community Composition
On average, 41,385 bacterial sequences per sample survived from quality trimming and chimera removal, which were subsequently normalized to 27310 sequences per sample and assembled into 4956 OTUs with >97% similarity. We found that 729 and 217 OTUs were detected specifically in the control and AMF inoculation treatments, respectively, and 4002 OTUs were shared between treatments (for detailed OTU distribution, see Supplementary Table S1). There were eight bacterial OTUs detected in the original soil collected from the forest that were not detected in both the control and AMF inoculation treatments. Most of the obtained OTUs belonged to the phyla Proteobacteria, Acidobacteria, and Actinobacteria with a relative abundance of more than 70% per sample on average. Chloroflexi was the fourth most abundant phyla (5.2% relative abundance on average) followed by Planctomycetes, Gemmatimonadete, Verrucomicrobia, and Bacteroidetes (Figure 1B). Detailed information showing how the relative abundance of the 10 most abundant phyla changed over time and treatment has been given in Supplementary Table S2.
FIGURE 1. The relative abundance of the dominant bacteria class (A) and phyla (B) in soil samples derived from the different treatments and sampling times used in our study. Abbreviations: AM represents the treatment inoculated with arbuscular mycorrhizal fungi, NM represents the uninoculated treatment, while O stands for the original soil collected from the forest. Different numbers (1, 2, 3, and 4) represent the different sampling times (T90, T120, T150, and T180, respectively). Relative abundance (>1%) is based on the proportional frequencies of those DNA sequences that could be classified at the phylum or class levels. Phylogenic groups that account for less then 1% of all classified sequences are divided into the artificial group “others.”
The relative proportion of the phyla and classes differed slightly according to sampling time and F. mosseae inoculation. However, the effects of F. mosseae inoculation on the bacterial communities were detected at a fine taxonomic resolution (genera level). The presence of AMF slightly decreased the relative abundance of Acidobacteria and increased that of Actinobacteria at each harvest, compared to the non-AM treatment (Figure 1B). At the class level, the 10 most abundant bacterial classes, representing over 75% of total abundance, was also slightly influenced by time and by F. mosseae inoculation. In each treatment, the class composition of the bacterial community showed successive changes over time (Figure 1A). Detailed information showing how the relative abundance of the 10 most abundant classes changed over time and treatment has been given in Supplementary Table S3. F. mosseae inoculation consistently increased the relative abundance of Actinobacteria and Thermoleophilia and decreased the relative abundance of Gammaproteobacteria and Acidobacteria. Additionally, the relative abundance of the most dominant class, Alphaproteobacteria, increased in the early stages (T90 and T120) of litter decomposition but decreased in the later stages (T150 and T180) (Figure 1A). While the relative abundance of the other major classes of the bacterial community fluctuated at different times and with different treatments, no consistent effects due to F. mosseae inoculation during litter decomposition were observed.
Analysis of the bacterial community at fine taxonomic resolution across all samples resulted in the detection of 283 genera. When conducting a significance test of the top 30 genera by relative abundance, significant changes either by time or AM fungal inoculation were found for 18 of these 30 (P < 0.05) (Figure 2). We found that F. mosseae inoculation significantly changed the bacterial community dynamics overtime. Specifically, in the NM treatment, the following genera demonstrated significant changes over time: Bradyrhizobium, Burkholderia, Salinispora, Candidatus Xiphinematobacter, Dyella, Kaistobacter, Phenylobacterium, Alicyclobacillus, Halomonas, and Stenotrophomonas. Furthermore, the observed changes occurred in the later stage of litter decomposition for all these genera except for Burkholderia, Salinispora, and Halomonas. Whereas changes in the relative abundance of Salinispora and Halomonas were only noted in the early stage of decomposition. On the other hand, in the AMF treatment, the relative abundance of Burkholderia, Candidatus Xiphinematobacter, Phenylobacterium, Alicyclobacillus, Lactococcus, Telmatospirillum, Candidatus Nitrososphaera, and Rhodoplanes were found to change significantly over the sampling time. All the noted changes occurred in the later stages of litter decomposition, while only three genera also exhibited significant changes in the early (Burkholderia) and both early and later stages (Lactococcus and Rhodoplanes). Bradyrhizobium, Salinispora, Dyella, Kaistobacter, Halomonas, and Stenotrophomonas, which showed significant changes overtime in NM treatment, exhibited no changes overtime in AM fungal inoculation treatment.
FIGURE 2. Heatmap and accompanying cluster analysis (x-axis) of the relative abundance of dominant bacterial genera in all the soil samples with different treatments over different sampling times. Abbreviations: AM represents the treatment inoculated with arbuscular mycorrhizal fungi, NM represents the uninoculated treatment, while O stands for the original soil collected from the forest. Different numbers (1, 2, 3, and 4) represent the different sampling times (T90, T120, T150, and T180, respectively). Mean relative abundance is shown for each samples. For each genus, a significant difference either between different sampling times or treatments is noticed by an asterisk (P < 0.05). Treatments with a letter in common are not different at P < 0.05 according to the Duncan’s multiple range test. Differences were tested among treatment combinations (harvesting times (0 is not included) and arbuscular mycorrhizal inoculation or not) for each genus (per row). The relative abundance for each genus in each soil sample is colored in the shades of yellow (low relative abundance) to red (high relative abundance).
When comparing the difference between the two treatments, F. mosseae inoculation significantly changed the relative abundance of the following genera: Burkholderia, Cryocola, Salinispora, Dyella, Kaistobacter, Phenylobacterium, Lactococcus, Telmatospirillum, Candidatus Nitrososphaera, Stenotrophomonas, and Halomonas (Figure 2).
The Effect of AMF on the Soil Bacterial Diversity and Richness
Although we obtained high number of sequence reads per sample (27310), the rarefaction curves showed that the number of OTUs still increased with the number of sequences, without reaching a plateau (Figure 3). Thus, we used both Shannon diversity index and Chao1 richness estimator as proxy for bacterial diversity and richness. The presence of F. mosseae did not significantly change the alpha diversity and richness of the soil bacterial community, during litter decomposition. The Shannon diversity remained stable over time and across the different treatments, with an average value of 6.18 (Supplementary Figure S4A). The Chao1 richness estimator showed a similar trend with that of Shannon diversity (Supplementary Figure S4B).
FIGURE 3. Rarefaction curve of soil bacterial OTUs clustered at 97% sequence identity across different soil samples. Abbreviations: AM represents the treatment inoculated with arbuscular mycorrhizal fungi, NM represents the uninoculated treatment, while O stands for the original soil collected from the forest. Different numbers (1, 2, 3, and 4) represent the different sampling times (T90, T120, T150, and T180, respectively).
We assessed the soil bacterial OTU richness at the phyla and class levels. The results indicated that F. mosseae did not significantly affect OTU richness at both these levels. At phylum level, the OTU-rich phyla for the NM control were phylum Proteobacteria (612 OTUs) followed by, Chloroflexi (278 OTUs), Acidobacteria (259 OTUs), Actinobacteria (222 OTUs), candidate division AD3 (38 OTUs), candidate division WPS-2 (43 OTUs). A similar pattern was also found for the OTU-rich phyla found in the AMF inoculation treatment: Proteobacteria (620 OTUs) followed by Chloroflexi (289 OTUs), Acidobacteria (259 OTUs), Actinobacteria (222 OTUs), candidate division WPS-2 (43 OTUs), and candidate division AD3 (38 OTUs). At the class level, we found that the bacterial communities from both the control and AMF inoculation treatments had the highest number of OTUs assigned to Alphaproteobacteria (control = 266 OTUs; AMF = 263 OTUs) followed by Ktedonobacteria (control = 201 OTUs; AMF = 205 OTUs), Planctomycetia (control = 194 OTUs; AMF = 196 OTUs), Gammaproteobacteria (control = 123 OTUs; AMF = 124 OTUs) and, Actinobacteria (control = 98 OTUs; AMF = 95 OTUs). Detailed changes in the richness of the 10 most abundant OTUs, for each harvest, is shown in Supplementary Table S4 (phylum level) and Supplementary Table S5 (class level).
The Effect of AMF on the Soil Bacterial Community Structure
The PCA and PERMANOVA, based on the relative abundance of all detected OTUs, showed that sampling time and AMF treatment significantly affected bacterial community composition (Pseudo-Fsampling time = 1.70, P < 0.004; Pseudo-FAMF inoculation = 1.99, P < 0.022). The interaction between sampling time and AMF treatment was not significant (Pseudo-Fsampling time x AMF inoculation= 1.09, P < 0.299). The PCA plot also showed that the first two canonical axes explained 28.7 and 18.1% of the total variability (Figure 4). The strongest effect of AMF inoculation on bacterial community structure was observed in T150 (Figure 4). Correlation analysis confirmed that AM fungal inoculation and sampling time were significantly correlated with the changes observed in the soil bacterial communities during litter decomposition. Additionally, none of the nutrients were correlated with bacterial community structure in this system (Table 1).
FIGURE 4. Principal component analyses (PCA) of soil bacterial community composition structure under different treatments and over different sampling times based on all the OTUs. Abbreviations: AM represents the treatment inoculated with arbuscular mycorrhizal fungi, NM represents the uninoculated treatment, while O stands for the original soil collected from the forest. Different numbers (1, 2, 3, and 4) represent the different sampling times (T90, T120, T150, and T180, respectively). The data is shown at mean value for each axis (±SE, n = 4).
Previous studies have used a variety of methods to investigate the effects of AMF on soil bacterial communities. These include PLFA analysis (Welc et al., 2010; Leifheit et al., 2015; Gui et al., 2017) and PCR-denaturing gradient gel electrophoresis (DGGE) (Marschner and Baumann, 2003; Roesti et al., 2006; Solis-Dominguez et al., 2011), as well as high-throughput methods such as 454 Pyrosequencing and Illumina sequencing (Nuccio et al., 2013; Qin et al., 2014). Our use of Illumina sequencing to test the effect of F. mosseae on soil bacterial communities during litter decomposition showed that although bacterial diversity and richness remained unchanged, F. mosseae altered the bacterial community composition and structure at different taxonomical levels during the litter decomposition process, including some key genera related to C and N cycling. We also found that the influence of F. mosseae on soil bacterial communities changed over time.
Although the effects of AMF on soil bacterial communities and litter decomposition have been documented in several studies (Hodge et al., 2001; Welc et al., 2010; Cheng et al., 2012; Leifheit et al., 2015), the results are not consistent between studies and the mechanisms by which AMF influences these processes remain unclear. Specifically, the relationship between AMF and bacterial richness and diversity are still largely unknown. In our study, we show that bacterial richness and diversity are not significantly affected by AMF inoculation and remain relatively stable during the litter decomposition processes. Despite a lack of difference in bacterial richness and diversity, 217 bacterial OTUs were detected in the AM treatment only, and 729 were detected in the NM treatment only. This is likely as a result of competition between the bacteria associated with AMF and other general soil bacteria, as reported in several studies (Welc et al., 2010; Nuccio et al., 2013; Mechri et al., 2014).
The bacterial community structure investigated in our study was not strongly influenced by F. mosseae at coarse taxonomic resolution, such as at the phylum or class levels. However, at the genus level we found that AMF had significant effects on bacterial community structure. This result is in agreement with past studies, which reported changes in the composition of soil bacterial communities as a result of AMF (e.g., Welc et al., 2010; Nuccio et al., 2013; Mechri et al., 2014). Furthermore, these observed changes were found at each sampling time, but became more pronounced over time.
The bacterial communities from the forest soils used in our study were mainly composed of Proteobacteria, Acidobacteria, Actinobacteria, and Choroflexi, regardless of the sampling time or treatment. This phylum level profile was found to be similar to the results of other studies using soil samples from different environments (Lauber et al., 2009; Li et al., 2014). Furthermore, for each sample, Proteobacteria was the most common and abundant phylum in the soil, a result that is also in line with past studies (Bastian et al., 2009; Li et al., 2014). However, our results indicated that F. mosseae did influence the relative abundance of these groups. For example, F. mosseae inoculation increased the relative abundance of Actinobacteria, and decreased the relative abundance of Acidobacteria. Actinobacteria is commonly positively correlated with AMF (Carpenter-Boggs et al., 1995; Nuccio et al., 2013), and known to inhibit the growth of other microorganisms (Duraipandiyan et al., 2010). This suppression of other organisms by Actinobacteria is the likely cause of the observed decline in the relative abundance of Acidobacteria. Furthermore, members of Proteobacteria (e.g., Alpha-, Beta-, and Gammaproteobacteria) have been identified as fast-growing bacteria that respond positively to and utilize AM fungal exudates. This could explain the observed increase in the relative abundance of Alphaproteobacteria during the initial harvests. However, a decrease in the abundance of Betaproteobacteria and Gammaproteobacteria also indicated that there was rhizosphere competition from other fast-growing bacteria such as Actinobacteria. This result agrees with that of Nuccio et al. (2013), which showed a similar decrease in the Betaproteobacterial family due to AM fungal inoculation.
At the genus level we observed that the relative abundance of certain bacterial genera were affected by F. mosseae inoculation. These genera shared the characteristic ability to degrade organic compounds in the soil (Xie and Yokota, 2005; Pearce et al., 2012) and included Candidatus Solibacter, Dyella, and Phenylobacterium. This result indicated that F. mosseae could influence litter decomposition through its effects on different bacteria in the soil; however, its effects were not uniform across genera. For example, Candidatus Solibacter can produce enzymes that break down organic carbon (Pearce et al., 2012). F. mosseae inoculation resulted in a higher relative abundance of C. Solibacter. Similarly, although Dyella and Phenylobacterium are both able to degrade organic compounds, F. mosseae inoculation positively affected the relative abundance of Dyella only, while the relative abundance of Phenylobacterium was negatively affected during the later stages of decomposition. Furthermore, our results indicate that bacteria from the genus Burkholderia, a rhizosphere-colonizing bacteria with saprophytic abilities (Roberts et al., 1997), were suppressed at T90. We also found that Geobacter, which is able to oxidize organic compounds (Childers et al., 2002) was abundant throughout the decomposition process, although its relative abundance was not affected by F. mosseae inoculation.
Our findings agree with those of Ravnskov et al. (2002), who reported that Rhizophagus irregularis suppressed the growth and development of Burkholderia cepacia populations. However, Mansfeld-Giese et al. (2002) showed that the external mycelia of F. intraradices increased the population of Pseudomonas chlororaphis in the soil, whereas in our experiment Pseudomonas abundance did not change significantly with F. mosseae inoculation. Our work supports past findings showing that Pseudomonas is frequently detected in litter decomposition processes (Noll et al., 2010).
Our study confirms previous reports on the influence of AMF on soil bacterial communities and provides detailed insight into the changes that occur in bacterial communities due to the presence of AMF. In addition, we were able to show how the influence of AMF on bacterial communities changes over time. AMF significantly altered the bacterial community dynamics and structures at fine taxonomic resolution, with no noticeable changes occurring at the phylum and class levels. Our work demonstrates that analysis of bacterial communities using high-resolution culture independent methods can significantly improve our understanding of bacterial richness, diversity and community dynamics in complex soil environments.
HG carried out the experiment, performed the data analysis, and took lead on writing the manuscript. PM developed the original ideas, oversaw the research work, and contributed toward writing the manuscript. WP contributed toward data analysis and manuscript revision. JX and KH contributed toward writing the manuscript.
This research was supported by the 973 key project of the National Natural Science Foundation of China (Y413631231) and the Chinese Ministry of Science and Technology, under the 12th 5-year National Key Technology Support Program (NKTSP) 2013BAB07B06 integration and comprehensive demonstration of key technologies on Green Phosphate-mountain Construction. We would like to thank the support from CGIAR Research Program 6: Forest, Trees and Agroforestry, for partially funding this work. We also would like to thank the help from Biological technology open platform, Kunming Institute of Botany, Chinese Academy of Sciences.
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.
- ^ http://ccb.jhu.edu/software/FLASH/
- ^ http://qiime.org/index.html
- ^ http://www.drive5.com/usearch/manual/uchime_algo.html
- ^ http://drive5.com/uchime/uchime_download.html
- ^ http://greengenes.lbl.gov/cgi-bin/nph-index.cgi
- ^ http://sourceforge.net/projects/rdp-classifier/
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01120/full#supplementary-material
Bastian, F., Bouziri, L., Nicolardot, B., and Ranjard, L. (2009). Impact of wheat straw decomposition on successional patterns of soil microbial community structure. Soil Biol. Biochem. 41, 262–275. doi: 10.1016/j.soilbio.2008.10.024
Biddle, J. F., Fitz-Gibbon, S., Schuster, S. C., Brenchley, J. E., and House, C. H. (2008). Metagenomic signatures of the Peru Margin subseafloor biosphere show a genetically distinct environment. Proc. Natl. Acad. Sci. U.S.A. 105, 10583–10588. doi: 10.1073/pnas.0709942105
Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59. doi: 10.1038/NMETH.2276
Bugg, T. D. H., Ahmad, M., Hardiman, E. M., and Singh, R. (2011). The emerging role for bacteria in lignin degradation and bio-product formation. Curr. Opin. Biotechnol. 22, 394–400. doi: 10.1016/j.copbio.2010.10.009
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
Carpenter-Boggs, L., Loynachan, T. E., and Stahl, P. D. (1995). Spore germination of Gigaspora margarita stimulated by volatiles of soil-isolated actinomycetes. Soil Biol. Biochem. 27, 1445–1451. doi: 10.1016/0038-0717(95)00075-P
Cheng, L., Booker, F. L., Tu, C., Burkey, K. O., Zhou, L. S., Shew, H. D., et al. (2012). Arbuscular mycorrhizal fungi increase organic carbon decomposition under elevated CO2. Science 337, 1084–1087. doi: 10.1126/Science.1224304
Cruz, A. F., and Ishii, T. (2011). Arbuscular mycorrhizal fungal spores host bacteria that affect nutrient biodynamics and biocontrol of soil-borne plant pathogens. Biol. Open 1, 52–57. doi: 10.1242/bio.2011014
DeAngelis, K. M., Chivian, D., Fortney, J. L., Arkin, A. P., Simmons, B., Hazen, T. C., et al. (2013). Changes in microbial dynamics during long-term decomposition in tropical forests. Soil Biol. Biochem. 66, 60–68. doi: 10.1016/j.soilbio.2013.06.010
Duraipandiyan, V., Sasi, A. H., Islam, V. I. H., Valanarasu, M., and Ignacimuthu, S. (2010). Antimicrobial properties of actinomycetes from the soil of Himalaya. J. Mycol. Med. 20, 15–20. doi: 10.1016/j.mycmed.2009.11.002
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Gui, H., Hyde, K., Xu, J., and Mortimer, P. (2017). Arbuscular mycorrhiza enhance the rate of litter decomposition while inhibiting soil microbial community development. Sci. Rep. 7:42184. doi: 10.1038/srep42184
Haas, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110
Hanackova, Z., Koukol, O., Stursova, M., Kolarik, M., and Baldrian, P. (2015). Fungal succession in the needle litter of a montane Picea abies forest investigated through strain isolation and molecular fingerprinting. Fungal Ecol. 13, 157–166. doi: 10.1016/j.funeco.2014.09.007
Hedley, M. J., Stewart, J. W. B., and Chauhan, B. S. (1982). Changes in inorganic and organic soil-phosphorus fractions induced by cultivation practices and by laboratory incubations. Soil Sci. Soc. Am. J. 46, 970–976. doi: 10.2136/sssaj1982.03615995004600050017x
Herman, D. J., Firestone, M. K., Nuccio, E., and Hodge, A. (2012). Interactions between an arbuscular mycorrhizal fungus and a soil microbial community mediating litter decomposition. FEMS Microbiol. Ecol. 80, 236–247. doi: 10.1111/J.1574-6941.2011.01292.X
Hodge, A. (2001). Arbuscular mycorrhizal fungi influence decomposition of, but not plant nutrient capture from, glycine patches in soil. New Phytol. 151, 725–734. doi: 10.1046/j.0028-646x.2001.00200.x
Hodge, A., Campbell, C. D., and Fitter, A. H. (2001). An arbuscular mycorrhizal fungus accelerates decomposition and acquires nitrogen directly from organic material. Nature 413, 297–299. doi: 10.1038/35095041
Iffis, B., St-Arnaud, M., and Hijri, M. (2014). Bacteria associated with arbuscular mycorrhizal fungi within roots of plants growing in a soil highly contaminated with aliphatic and aromatic petroleum hydrocarbons. FEMS Microbiol. Lett. 358, 44–54. doi: 10.1111/1574-6968.12533
Koller, R., Rodriguez, A., Robin, C., Scheu, S., and Bonkowski, M. (2013). Protozoa enhance foraging efficiency of arbuscular mycorrhizal fungi for mineral nitrogen from organic matter in soil to the benefit of host plants. New Phytol. 199, 203–211. doi: 10.1111/nph.12249
Kuramae, E. E., Hillekens, R. H. E., de Hollander, M., van der Heijden, M. G. A., van den Berg, M., van Straalen, N. M., et al. (2013). Structural and functional variation in soil fungal communities associated with litter bags containing maize leaf. FEMS Microbiol. Ecol. 84, 519–531. doi: 10.1111/1574-6941.12080
Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. doi: 10.1128/AEM.00335-09
Lecomte, J., St-Arnaud, M., and Hijri, M. (2011). Isolation and identification of soil bacteria growing at the expense of arbuscular mycorrhizal fungi. FEMS Microbiol. Lett. 317, 43–51. doi: 10.1111/j.1574-6968.2011.02209.x
Leifheit, E. F., Verbruggen, E., and Rillig, M. C. (2015). Arbuscular mycorrhizal fungi reduce decomposition of woody plant litter while increasing soil aggregation. Soil Biol. Biochem. 81, 323–328. doi: 10.1016/J.Soilbio.2014.12.003
Leigh, J., Hodge, A., and Fitter, A. H. (2009). Arbuscular mycorrhizal fungi can transfer substantial amounts of nitrogen to their host plant from organic material. New Phytol. 181, 199–207. doi: 10.1111/J.1469-8137.2008.02630.X
Li, Y. Y., Wen, H. Y., Chen, L. Q., and Yin, T. T. (2014). Succession of bacterial community structure and diversity in soil along a chronosequence of reclamation and re-vegetation on coal mine spoils in China. PLoS ONE 9:e115024. doi: 10.1371/journal.pone.0115024
Lynd, L. R., Weimer, P. J., van Zyl, W. H., and Pretorius, I. S. (2002). Microbial cellulose utilization: fundamentals and biotechnology. Microbiol. Mol. Biol. Rev. 66, 739–739. doi: 10.1128/MMBR.66.4.739.2002
Mansfeld-Giese, K., Larsen, J., and Bodker, L. (2002). Bacterial populations associated with mycelium of the arbuscular mycorrhizal fungus Glomus intraradices. FEMS Microbiol. Ecol. 41, 133–140. doi: 10.1016/S0168-6496(02)00265-9
Masai, E., Katayama, Y., and Fukuda, M. (2007). Genetic and biochemical investigations on bacterial catabolic pathways for lignin-derived aromatic compounds. Biosci. Biotechnol. Biochem. 71, 1–15. doi: 10.1271/bbb.60437
McGonigle, T., Miller, M., Evans, D., Fairchild, G., and Swan, J. (1990). A new method which gives an objective measure of colonization of roots by vesicular-arbuscular mycorrhizal fungi. New Phytol. 115, 495–501. doi: 10.1111/j.1469-8137.1990.tb00476.x
Mechri, B., Manga, A. G. B., Tekaya, M., Attia, F., Cheheb, H., Ben Meriem, F., et al. (2014). Changes in microbial communities and carbohydrate profiles induced by the mycorrhizal fungus (Glomus intraradices) in rhizosphere of olive trees (Olea europaea L.). Appl. Soil Ecol. 75, 124–133. doi: 10.1016/j.apsoil.2013.11.001
Noll, M., Naumann, A., Ferrero, F., and Malow, M. (2010). Exothermic processes in industrial-scale piles of chipped pine-wood are linked to shifts in gamma-, alphaproteobacterial and fungal ascomycete communities. Int. Biodeterior. Biodegradation 64, 629–637. doi: 10.1016/j.ibiod.2010.07.002
Nuccio, E. E., Hodge, A., Pett-Ridge, J., Herman, D. J., Weber, P. K., and Firestone, M. K. (2013). An arbuscular mycorrhizal fungus significantly modifies the soil bacterial community and nitrogen cycling during litter decomposition. Environ. Microbiol. 15, 1870–1881. doi: 10.1111/1462-2920.12081
Pearce, D. A., Newsham, K. K., Thorne, M. A., Calvo-Bado, L., Krsek, M., Laskaris, P., et al. (2012). Metagenomic analysis of a southern maritime Antarctic soil. Front. Microbiol. 3:403. doi: 10.3389/fmicb.2012.00403
Persoh, D., Segert, J., Zigan, A., and Rambold, G. (2013). Fungal community composition shifts along a leaf degradation gradient in a European beech forest. Plant Soil 362, 175–186. doi: 10.1007/s11104-012-1271-y
Purahong, W., Kapturska, D., Pecyna, M. J., Jariyavidyanont, K., Kaunzner, J., Juncheed, K., et al. (2015). Effects of forest management practices in temperate beech forests on bacterial and fungal communities involved in leaf litter degradation. Microb. Ecol. 69, 905–913. doi: 10.1007/s00248-015-0585-8
Purahong, W., Wubet, T., Lentendu, G., Schloter, M., Pecyna, M. J., Kapturska, D., et al. (2016). Life in leaf litter: novel insights into community dynamics of bacteria and fungi during litter decomposition. Mol. Ecol. 25, 4059–4074. doi: 10.1111/mec.13739
Qin, H., Brookes, P. C., Xu, J. M., and Feng, Y. Z. (2014). Bacterial degradation of Aroclor 1242 in the mycorrhizosphere soils of zucchini (Cucurbita pepo L.) inoculated with arbuscular mycorrhizal fungi. Environ. Sci. Pollut. Res. 21, 12790–12799. doi: 10.1007/s11356-014-3231-y
Ravnskov, S., Larsen, J., and Jakobsen, I. (2002). Phosphorus uptake of an arbuscular mycorrhizal fungus is not effected by the biocontrol bacterium Burkholderia cepacia. Soil Biol. Biochem. 34, 1875–1881. doi: 10.1016/S0038-0717(02)00201-8
Roberts, D., Dery, P., Mao, W., and Hebbar, P. (1997). Use of a colonization-deficient strain of Escherichia coli in strain combinations for enhanced biocontrol of cucumber seedling diseases. J. Phytopathol. 145, 461–463. doi: 10.1111/j.1439-0434.1997.tb00350.x
Roesti, D., Gaur, R., Johri, B. N., Imfeld, G., Sharma, S., Kawaljeet, K., et al. (2006). Plant growth stage, fertiliser management and bio-inoculation of arbuscular mycorrhizal fungi and plant growth promoting rhizobacteria affect the rhizobacterial community structure in rain-fed wheat fields. Soil Biol. Biochem. 38, 1111–1120. doi: 10.1016/j.soilbio.2005.09.010
Saia, S., Benitez, E., Garcia-Garrido, J. M., Settanni, L., Amato, G., and Giambalvo, D. (2014). The effect of arbuscular mycorrhizal fungi on total plant nitrogen uptake and nitrogen recovery from soil organic material. J. Agric. Sci. 152, 370–378. doi: 10.1017/S002185961300004X
Schellenberger, S., Kolb, S., and Drake, H. L. (2010). Metabolic responses of novel cellulolytic and saccharolytic agricultural soil Bacteria to oxygen. Environ. Microbiol. 12, 845–861. doi: 10.1111/j.1462-2920.2009.02128.x
Scheublin, T. R., Sanders, I. R., Keel, C., and van der Meer, J. R. (2010). Characterisation of microbial communities colonising the hyphal surfaces of arbuscular mycorrhizal fungi. ISME J. 4, 752–763. doi: 10.1038/ismej.2010.5
Smith, S. E., and Smith, F. A. (2011). Roles of arbuscular mycorrhizas in plant nutrition and growth: new paradigms from cellular to ecosystem scales. Annu. Rev. Plant Biol. 62, 227–250. doi: 10.1146/annurev-arplant-042110-103846
Solis-Dominguez, F. A., Valentin-Vargas, A., Chorover, J., and Maier, R. M. (2011). Effect of arbuscular mycorrhizal fungi on plant biomass and the rhizosphere microbial community structure of mesquite grown in acidic lead/zinc mine tailings. Sci. Total Environ. 409, 1009–1016. doi: 10.1016/j.scitotenv.2010.11.020
Stursova, M., Zifcakova, L., Leigh, M. B., Burgess, R., and Baldrian, P. (2012). Cellulose utilization in forest litter and soil: identification of bacterial and fungal decomposers. FEMS Microbiol. Ecol. 80, 735–746. doi: 10.1111/j.1574-6941.2012.01343.x
Toljander, J. F., Artursson, V., Paul, L. R., Jansson, J. K., and Finlay, R. D. (2006). Attachment of different soil bacteria to arbuscular mycorrhizal fungal extraradical hyphae is determined by hyphal vitality and fungal species. FEMS Microbiol. Lett. 254, 34–40. doi: 10.1111/j.1574-6968.2005.00003.x
Treseder, K. K., Bent, E., Borneman, J., and Mcguire, K. L. (2014). Shifts in fungal communities during decomposition of boreal forest litter. Fungal Ecol. 10, 58–69. doi: 10.1016/j.funeco.2013.02.002
Welc, M., Ravnskov, S., Kieliszewska-Rokicka, B., and Larsen, J. (2010). Suppression of other soil microorganisms by mycelium of arbuscular mycorrhizal fungi in root-free soil. Soil Biol. Biochem. 42, 1534–1540. doi: 10.1016/j.soilbio.2010.05.024
White, D., MacNaughton, S., Pankhurst, C., Doube, B., and Gupta, V. (1997). Chemical and molecular approaches for rapid assessment of the biological status of soils. Biol. Indic. Soil Health 15, 371–396.
Keywords: arbuscular mycorrhizal fungi, soil bacteria, Illumina sequencing, litter decomposition, soil microbial community
Citation: Gui H, Purahong W, Hyde KD, Xu J and Mortimer PE (2017) The Arbuscular Mycorrhizal Fungus Funneliformis mosseae Alters Bacterial Communities in Subtropical Forest Soils during Litter Decomposition. Front. Microbiol. 8:1120. doi: 10.3389/fmicb.2017.01120
Received: 22 December 2016; Accepted: 01 June 2017;
Published: 20 June 2017.
Edited by:Pierre-Emmanuel Courty, Institut National de la Recherche Agronomique (INRA), France
Reviewed by:Sergio Saia, Council for Agricultural Research and Agricultural Economy Analysis, Italy
Erika Buscardo, University of Coimbra, Portugal
Copyright © 2017 Gui, Purahong, Hyde, Xu and Mortimer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.