Xue Zhou
Dario Fornara
Makoto Ikenaga4
Ruifu Zhang
Zhongjun Jia- 1State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Agri-Food and Biosciences Institute, Belfast, Ireland
- 4Research Field in Agriculture, Agriculture Fisheries and Veterinary Medicine Area, Kagoshima University, Kagoshima, Japan
- 5National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agriculture University, Nanjing, China
Forest soil ecosystems are associated with large pools and fluxes of carbon (C) and nitrogen (N), which could be strongly affected by variation in rainfall events under current climate change. Understanding how dry and wet cycle events might influence the metabolic state of indigenous soil microbes is crucial for predicting forest soil responses to environmental change. We used 454 pyrosequencing and quantitative PCR to address how present (DNA-based) and potentially active (RNA-based) soil bacterial communities might response to the changes in water availability across three different forest types located in two continents (Africa and Asia) under controlled drying and rewetting cycles. Sequencing of rRNA gene and transcript indicated that Proteobacteria, Actinobacteria, and Acidobacteria were the most responsive phyla to changes in water availability. We defined the ratio of rRNA transcript to rRNA gene abundance as a key indicator of potential microbial activity and we found that this ratio was increased following soil dry-down process whereas it decreased after soil rewetting. Following rewetting Crenarchaeota-like 16S rRNA gene transcript increased in some forest soils and this was linked to increases in soil nitrate levels suggesting greater nitrification rates under higher soil water availability. Changes in the relative abundance of (1) different microbial phyla and classes, and (2) 16S and amoA genes were found to be site- and taxa-specific and might have been driven by different life-strategies. Overall, we found that, after rewetting, the structure of the present and potentially active bacterial community structure as well as the abundance of bacterial (16S), archaeal (16S) and ammonia oxidizers (amoA), all returned to pre-dry-down levels. This suggests that microbial taxa have the ability to recover from desiccation, a critical response, which will contribute to maintaining microbial biodiversity in harsh ecosystems under environmental perturbations, such as significant changes in water availability.
Introduction
Climate change effects on seasonal precipitation regimes and on rates of soil evapotranspiration are predicted to strongly influence soil water availability across forest ecosystems worldwide through its effects on both seasonal precipitation regimes and soil evapotranspiration (Dale et al., 2001; Westerling et al., 2006; Malhi et al., 2008). Seasonal changes in rainfall events might lead to cooler-wetter winters and hotter-drier summers with important consequences for soil nutrient mobilization and soil microbiological processes. In particular large increases in water availability during intense rainfall events occurring at the end of prolonged dry seasons might greatly enhance soil CO2 efflux pulses. Such increase in soil respiration is primarily due to enhanced microbial mineralization of carbon (C) and nitrogen (N) organic substrates and is referred to as “Birch Effect” (Birch, 1958; Borken and Matzner, 2009; Inglima et al., 2009). This effect well explains large CO2 pulses, which result from rewetting events after dry summer periods and accounts for a significant proportion of annual C budgets in forest ecosystems (Xu and Qi, 2001). Despite soil microbial communities may play a crucial role in mediating the “Birch Effect” it is not clear how soil microbial abundance and community structure will change under drying and rewetting cycles in forest ecosystems.
Previous research studies have mainly focused on the chemical process of microbial mineralization, which is associated with high rates of C (Fierer and Schimel, 2002; Miller et al., 2005) and N release during the dry-wet transition in forest soil ecosystems (Franzluebbers, 1999; Fierer and Schimel, 2002; Mikha et al., 2005). Similar studies show that response patterns of soil microbial biomass to rainfall events could be site-dependent whereby microbial biomass can either increase (Basu et al., 1991; Jia and Insam, 1991; Luizao et al., 1998) or decrease during the rainy season (Singh et al., 1989; Barbhuiya et al., 2004). Microorganisms have an important role in the regulation of ecosystem processes such as nutrient cycling. The abundance and diversity of microbial communities has important implications for the stability and function of nature and semi-natural ecosystems (Hashsham et al., 2000; Bell et al., 2005). So far, field investigations addressing potential responses of soil microorganisms to dry-wet cycles have been limited by the use of techniques with relatively low taxonomic resolution. Recent findings show how the use of clone libraries can help better understanding the response of soil fungal and bacterial communities to climate change (Castro et al., 2010). By using 16S-based denaturing gradient gel electrophoresis (DGGE), Griffiths et al. (2003) have also demonstrated how bacterial communities in grassland soils can be highly resistant to water stress.
Desiccation will contribute to increasing nutrient-poor soil environments with negative effects on the metabolic state of microorganisms as it was shown in both natural communities and laboratory strains (Braeken et al., 2006). Survival of bacterial under desiccation conditions has been studied in different soil environments such as in soils under rhizoremediation (Vilchez and Manzanera, 2011) or in grassland soils (Barnard et al., 2013). In these studies, a large variation of survival times has been observed where for example, Gram-positive bacteria tend to better survive to desiccation than Gram-negative organisms. Changes in specific life strategies of microorganisms were also documented as key response to drought events. However, most of the experimental work so far has been culture-dependent often failing to address the effects of repeated cycles in water availability on microbial dormancy mechanisms in response to desiccation.
Despite changes in microbial community composition and structure under drying and rewetting cycles can greatly influence soil nutrient availability and productivity in forest ecosystems, the net response of soil microorganisms to multiple dry/wet cycles is still poorly understood. Here we suggest that high-taxonomic-resolution techniques can provide detailed phylogenetic-level data (Shi et al., 2015), which could help better understanding the genetic basis of dormancy regulation in soil ecosystems. The proportional change of certain phenotypes could be linked with changes in key variables of the complex soil matrix without the need to use specific biomarkers. High-pyrosequencing based on bacterial microbial DNA and RNA allows the phylogenetic characterization of bacterial groups (rRNA genes), which have the activity of synthetic proteins (rRNA). It has been then suggested that rRNA is a reliable indicator of metabolic state in microbial assemblages (Schippers et al., 2005; Jones and Lennon, 2010; Blazewicz et al., 2013). Thus the main aim of this study is to investigate the potential response of the soil microbial community to extreme desiccation events and subsequent rewetting across three different forest sites. First we used 454 pyrosequencing to analyze ribosomal RNA genes and transcripts from soil samples, which have been exposed to drying-rewetting cycles under different climatic conditions. Second we measured the abundance of selected phylogenetic marker genes (i.e., bacterial and archaeal 16S rRNA and amoA gene) and transcripts by performing quantitative PCR. By addressing how soil microbial communities might respond to drying-rewetting cycles, in forest ecosystems, our study contributes to improving our mechanistic understanding of how microorganisms can ultimately deal with significant fluctuations in soil water availability.
Materials and Methods
Site Description and Soil Sampling
In our study, we choose three typical forest soils with long distance and similar precipitation pattern, including: (i) African forest soil was collected from Egerton University (0°22′S, 35°55′E), Kenya. The annual mean temperature is 17.7°C, and the annual mean precipitation is 1049 mm. (ii) Chinese forest soil was sampled from the Tea Research Institute of the Chinese Academy of Agricultural Sciences (120°09′E, 30°14′N), Hangzhou City. The annual mean temperature is 17.0°C, and the annual mean precipitation is 1533 mm. (iii) Japanese forest soil was harvested from Minami–Kyushu city, Kagoshima (31°21′E, 130°26′N), Japan. The annual mean temperature is 18.8°C, and the annual mean precipitation is 2400 mm.
At each forest site, three large independent plots were randomly selected under same topography and vegetation type, and three replicate grids were then randomly established within each plot. Then within each grid three replicate soil samples were collected from the A horizon between 0 and 20 cm depth. One composite soil sample was calculated by homogenization of replicated samples after removal of plant residues, roots, stones, and obvious macrofauna. Immediately after transportation to the laboratory, soil samples were sieved through a 2 mm mesh before storage at 4°C until further use. Sub-samples were air-dried for physical and chemical analysis.
Laboratory Dry-down and Wet-up
Soils were exposed to drying and rewetting treatments during the incubation at 28°C. There were 12 replicates for each soil so that soils could be destructively sampled in triplicate after each dry-down and wet-up incubation. The treatment consisted of two dry/wet cycles (Figure 1A). One dry/wet cycle started with a 7-days dry-down period, and followed by a 7-days rewetting, which consisted of rewetting to ca. 40% soil maximum water holding capacity (WHC). Soils were incubated in sealed polythene bottles. The soils were dried to constant weight by removing the jar lid and blowing dry air across the soil for 7 days. Soils were rewetted by slowly adding sterile deionized water into the jar using a syringe until the soils reached to 40% WHC, and incubated for 7 days. The total incubation time for the 2-cycle treatment was 28 days.
FIGURE 1. (A) Dynamics of soil water content during dry-wet cycles in three experimental soils (D: no water input, F: weekly water inputs). (B) Dynamics of soil total soil organic carbon during dry-wet cycles. (C) Dynamics of total soil nitrogen during dry-wet cycles. (D) Changes in nitrate and ammonium concentration in soil microcosms in response to changes in water availability. Bars indicate standard deviations of triplicate soil samples. Abbreviations represent different treatments: fresh soil (F1), 1st dry-down (D1), 1st wet-up (F2), 2nd dry-down (D2), and 2nd wet-up (F3). The designation C1 represent the first cycle from fresh to dry condition and C2 represent the second cycle from fresh to dry condition.
Soil Properties
Soil properties were determined as previously described (Zhou et al., 2015). Soil water content was measured at 105°C for 8 h. Soil pH was determined in 2.5:1 (w:v) ratios of soil with distilled water using a DMP-2 mV/pH detector (Quark Ltd, Nanjing, China). Total soil C & N (%) were determined simultaneously by the Dumas method using a LECO CN 2000. NH4+-N and NO3--N were extracted with 2 M KCl at a soil/solution ratio of 1:5 by shaking at 200 rpm for 60 min and determined by a Continuous Flow Analyser (San++System, Skalar, Holland).
Soil Nucleic Acid Extraction
Soil nucleic acids were extracted from 0.5 g soil using the FastDNA spin kit for soil (Qbiogene, Inc., Irvine, CA, USA) according to the manufacturer’s instructions. DNA purification was conducted with 5.5 M guanidine thiocyanate solution by removing humic substance contamination. The quality and concetration of the extracted DNA was measured by gel electrophoresis (0.8% agarose) and Nanodrop ND-1000 UV-Vis spectrophotometer (NanoDrop Technologies Inc, Wilmington, DE, USA), and stored at -20°C.
RNA was extracted using the protocol of Griffiths et al. (2000) with the modification that glass beads was performed twice. The RNA extraction method is detailed as described in (Ding et al., 2015). rRNA was purified by the RNeasy mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions, and quantified using UV-vis Spectrophotometer (ND-1000 NanoDrop). otal RNA was converted to cDNA using PrimeScript 1st Strand cDNA Synthesis Kit (TaKaRa) using random hexamers, and stored at -20°C for amplicon sequencing.
Pyrosequencing Analysis of 16S rRNA Gene
Pyrosequencing of total 16S rRNA genes of the V4 regions was carried out on the Roche 454 GS FLX Tianium sequencer (Roche Diagnostics Corporation, Branford, CT, USA) using the putatively universal primer pairs of Tag-515F and Tag-907R (Supplementary Table S1), as described previously (Xia et al., 2011). The PCR reaction performed in a thermal cycler (Bio-Rad Laboratories, Hercules, CA, USA) under the following procedure: 94°C for 5 min, followed by 30 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 55 s, with a 10 min extension at 72°C (Supplementary Table S1). The reaction mixture consisted of 0.2 mM of each deoxynucleoside triphosphate, 1 × PCR buffer (Mg2+ Plus), 0.1 μM of each primer, 1.25 U of TaKaRa Taq HS polymerase (TaKaRa Biotech, Dalian, China), and 50 ng of soil DNA template in a final volume of 50 μl. PCR products were checked by 1.2% agarose gels, and purified using the Agarose Gel DNA Purification Kit (TaKaRa), and then the PCR amplicons were combined in equimolar ratios into a single tube in preparation for pyrosequencing analysis with soil samples from each treatment.
The pyrosequencing data of 16S rRNA genes were analyzed using the Quantitative Insights Into Microbial Ecology (QIIME) soft package (Caporaso et al., 2010), and sequences below quality score of 25 or 200 bp in length were excluded before further analysis. A total of 310, 308 high-quality sequence reads of 16S rRNA gene were obtained with an average length of ∼390 bp (Supplementary Table S2). All samples were rarefied to 2,300 sequences for statistical analysis including rarefaction curves of Shannon diversity (Supplementary Figure S1). The taxonomic identity of each phylotype was determined using ribosomal database project (RDP) Classifier with a confidence threshold of 0.80 (Cole et al., 2005). Finally, sequences were clustered into operational taxonomic units (OTUs) based on 97% sequence similarity with UCLUST (Edgar, 2010). Based on the sequences and/or OTUs obtained, microbial diversity, and richness analyses, principal component analysis, and canonical correspondence analysis were performed. The Crenarchaeota-like 16S rRNA genes were screened by RDP classifier at phyla level (Cole et al., 2009), while bacterial ammonia oxidizers of Nitrosococcus and Nitrosospira and nitrite oxidizers of Nitrobacter and Nitrospira were assigned at the genus level for further analysis.
Population Size
The abundance of genes and transcripts were assessed by quantitative polymerase chain reaction (qPCR) on a CFX96 Optical Real-Time Detection System (Bio-Rad Laboratories, Inc. Hercules, CA, USA). The population size of bacterial and archaeal communities was assessed by quantification of 16S rRNA gene copies in soil samples using the primer pairs 515F–907R (Stubner, 2002) and A364aF–A934bR (Kemnitz et al., 2005), respectively. AOA and AOB in soil samples were quantified with the primer pairs of Arch-amoAF/Arch-amoAR (Francis et al., 2005) and amoA-1F/amoA-2R (Rotthauwe et al., 1997), respectively.
The DNA and cDNA templates were 10-fold diluted to 1–10 ng in each reaction mixture. Quantitative PCR reactions were carried out in 20-μL reaction mixtures containing 10 μL 2x SYBR Premix Ex Taq (Takara Biotech, Dalian, China), 200 nM of each primer, and 20 ng DNA template. The thermal program for the real-time PCR assay was shown in Supplementary Table S1. Standard curves, spanning 107–101, were constructed by dilution series of plasmids harboring the aimed gene. The amplification efficiency ranges from 90 to 115% and R2 values were approximately 0.997–0.999 in each reaction. Special amplification of target genes was confirmed by melting curve analysis always resulting in a single peak.
Statistical Analysis
Statistical analyses of pyrosequencing were performed in R software (Version 3.0.2, vegan package). We performed analysis of variance (ANOVA) to test for potential significant effects of soil moisture on the abundance of 16S rRNA genes and transcripts across our microbial communities. We performed multiple Student–Newman–Keuls tests to identify significant differences between each treatment. We also used MRPP, ANOSIM, and ADONIS to test for group difference by using R software. All data analyses were performed using SPSS version 10.0 (IBM Co., Armonk, NY, USA). Differences at P < 0.05 were considered statistically significant. The pyrosequencing reads of the total 16S rRNA genes have been deposited in the European Nucleotide Archive under the accession number Hx2000053322.
Results
Soil C and N Pools under Drying-rewetting Cycles
We found that soil organic C and N content did not change during the dry-down and wet-up treatments across the three forest sites (Figures 1B,C). Soil NH4+ and soil NO3- pools significantly responded to changes in water availability (Figure 1D). In particular, NH4+ levels significantly increased with greater water availability in soils collected from Asian (Chinese and Japanese) forest sites (from 23.48 to 124.65 ppm and from 45.28 to 156.77 ppm, respectively). However, soil NO3- levels were not significantly affected by changes in water availability in Chinese soil and soil NO3- was slightly increased in the second wet-up in Japanese soil. Soil NH4+ and NO3- pools followed the opposite pattern in the African (AS) forest soils. Here soil NH4+ increased on average by 336% with repetitive dry-down treatments, and declined by 87.5% after repetitive wet-up treatments, whereas large amounts of NO3 accumulated during repetitive wet-up (+97.1%) but decreased after dry-down treatment (by average -17.4%).
Quantification of Bacterial and Archaeal 16S rRNA Genes and Transcripts
Population size of bacterial and archaeal 16S genes (DNA) and transcripts (RNA) differed among sites over dry-down and wet-up treatments (Figures 2A,B). In the African and Chinese forest soils repetitive dry-down treatments decreased the overall abundance of bacterial and archaeal 16S rRNA gene. In contrast, repetitive wet-up stimulated the abundance of bacterial and archaeal 16S rRNA gene back to the original level. However, neither the abundance of bacterial nor archaeal 16S rRNA genes was significantly affected by the change of water content in the Japanese soil.
FIGURE 2. Changes in the abundance of bacterial (A) and archaeal (B) 16S rRNA genes and transcripts in response to dry-down and wet-up treatments on the three forest soil types. Correlation of abundance of bacterial (C) and archaeal (D) 16S rRNA gene and transcript. The abundances were log-transformed to reduce skewness of the data. Red points represent dry-down treatment, and black points represent pre-dry-down and rewetting treatment. The designation AF1 represents the fresh African forest soil. The Africa D1 represents the 1st dry-down treatment of African forest soil. The designation Africa_F2 represents the 2nd wet-up treatment of African forest soil. The original data is shown in Supplementary Table S6.
In our study, we define rRNA transcipt: rRNA gene ratio as the indicator of potential activity of microorganisms (Figures 2C,D), as previous investigations (Kemp et al., 1993; Kerkhof and Ward, 1993; Poulsen et al., 1993; Muttray et al., 2001). In the African and Chinese soil, the bacterial rRNA transcript: bacterial rRNA gene ratio and archaeal rRNA transcript: archaeal rRNA gene ratio both stimulated by dry-down and declined by wet-up. However, neither the bacterial ratio nor the archaeal ratio was influenced by the change of water availability in Japanese soil.
Quantification of Ammonia-oxidizing Microbes
Real-time quantitative PCR was also employed to address changes in population size of AOA and AOB by quantifying archaeal and bacterial amoA genes (Figure 3A). The abundance of archaeal and bacterial amoA genes showed different response patterns to changes in water availability across the three forest sites. In African and Chinese soils, archaeal amoA genes significantly decreased with repetitive dry-down treatments (average -96.2 and -99%, respectively), and increased to the original level with the repetitive wet-up treatments. However, neither dry-down nor wet-up affected the abundance of AOA in Japanese soil (P > 0.05). Bacterial amoA gene was generally 2 to 3 orders magnitude lower than archaeal amoA genes across different soils. Bacterial amoA gene decreased with repetitive dry-down treatments (average -99%), and significantly increased to original levels under repetitive wet-up treatments (P < 0.05) in AS. In contrast, no detectable dry-down-related or wet-up-related changes in bacterial amoA gene were found in Chinese and Japanese soils, suggesting a marked resistance to changes in water availability.
FIGURE 3. (A) Changes in the abundance of bacterial and archaeal amoA genes in response to dry-down and wet-up treatments on the three forest soil types. Changes in the relative abundance of Crenarchaeota (AOA) (B), AOB (C), and NOB (D) in soil microcosms in response to changes in water availability. Bars indicate standard deviations of triplicate soil samples. All the designations are as the same as above.
Response of Microbial Community Structure
A total of 309,840 high-quality sequence reads of DNA and cDNA were obtained with length of ∼390 bp (Supplementary Table S1) with an average ± SE of 3443 ± 473 in each sample. Rarefaction curves of Shannon diversity index were used to estimate microbial diversity among dry-rewetting cycles in each plot (Supplementary Figure S1; Supplementary Table S3). All rarefaction curves reached a plateau indicating that sequence-derived diversity in this study was sufficient to characterize microbial species in each soil. Present (DNA-based) and potentially active (RNA-based) microbial diversities at each site responded differently to changes in water availability. In African and Japanese soil, present (DNA-based) microbial diversity remained stable under both dry-down and wet-up treatments. However, potentially (RNA-based) active microbial diversity decreased under both dry-down and rewetting treatments in these two soils. On the contrary, present (DNA-based) microbial diversity increased significantly as dry-down progressed in Chinese soil, then returned to the original diversity after rewetting. Potentially (RNA-based) active diversity displayed remarkable resistance to drying and rewetting cycles.
Figure 4 shows shifts in the microbial communities during the change of water availability in non-metric multidimensional scaling (NMDS) plots of the pairwise UniFrac distance ordinations. We found that F1, F2, F3 were clustered and D1, D2 were clustered in DNA of Chinese soil, and in RNA of African and Japanese soil (P < 0.05) (Supplementary Table S4).
FIGURE 4. (A) The compositional structure of 16S rRNA genes and transcripts among three forest soils described by non-metric multi-dimensional scaling (NMDS) plots. (B) The change pattern of bacterial 16S rRNA genes and transcripts in response to dry-down and wet-up treatments described by NMDS plots.
Response of Bacterial Groups and Nitrifying Communities to Dry-down and Wet-up
Bacterial 16S rRNA Gene
The relative abundance of bacteria groups was analyzed by aggregating all taxonomic classification at the phylum and class levels. Soil bacterial communities at the three forest sites were dominated in relative abundance by the Actinobacteria phylum (mostly Actinobacteria class), and the Proteobacteria phylum (primarily alpha, beta, gamma, and delta class). Although the bacterial community structure was not significantly influenced by dry-down and wet-up, the relative abundance of the different bacterial phyla and class followed different response patterns to changes in water availability (Table 1, Supplementary Table S5).
TABLE 1. Significant changes (%) in bacterial phylotypes at the taxonomic level of classes (P < 0.05) affected by the dry-down and wet-up.
Overall, there are 5, 4, and 2 phyla were significantly influenced by repeated dry-down and wet-up treatments in the African soil, Chinese soil and Japanese soil, respectively. Actinobacteria were the most responsive phylum to changes in water availability at African and Japanese sites. Actinobacteria increased with dry-down and decreased with wet-up. The relative abundance of Verrucomicrobia in African and Chinese soils followed a pattern similar to Actinobacteria. Additionally, Firmicutes (Africa), Acidobacteria (China), and Planctomycetes (China) phyla also increased with dry-down and decreased with wet-up treatments. The other main responsive phylum was Proteobacteria, which responded in different ways. Our data shows that dry-down decreased the relative abundance of Proteobacteria (mainly Betaproteobacteria and Gammaproteobacteria class) in Chinese and Japanese soils whereas wet-up increased their abundance. Proteobacteria (mainly Alphaproteobacteria class) in African soil showed the opposite response pattern to the dry-down and wet-up treatments. At the class level, Betaproteobacteria and Gammaproteobacteria were highly responsive to changes in water availability and decreased under dry-down and wet-up treatments at Chinese and Japanese sites.
Bacterial 16S rRNA Transcript
In our study, changes in the relative abundance of rRNA appeared to reflect changes in potential microbial activity. Each soil has 4 phyla which were significantly influenced by the change of water availability. The community-dominating Actinobacteria transcripts were similar in African and Chinese soils and increased with dry-down and decreased with wet-up treatments. Similarly, the Firmicutes phylum followed the Actinobacteria pattern in Japanese soil. In contrast, Acidobacteria showed the opposite change pattern in Chinese and Japanese soils, which declined with dry-down and accumulated with wet-up. Likewise, Chloroflexi transcripts followed patterns similar to that of Acidobacteria. At the class level, Alphaproteobacteria (China) and Gammaproteobacteria and Betaproteobacteria (Japan) were the most responsive class in Chinese and Japanese soils, respectively, and decreased with dry-down and increased with wet-up cycles. Interestingly, the Bacilli class significantly increased under dry-down and sharply decreased with wet-up treatment in Japanese soils.
Nitrifying Communities
Proportional changes in nitrifying population were assessed by screening 16S rRNA gene and transcript sequences of some dominant phylotypes in ammonia oxidizers and nitrite oxidizers (Figures 3B–D). Among crearchaeota, we observed a decline in the proportion of crearchaeota with dry-down treatments (from 8 to 5%, and 12 to 4%, respectively), and an increase under wet-up treatments (from 5 to 12%, and 4 to 10%, respectively) in AS, which is consistent with the result from quantitative PCR. However, in Chinese and Japanese soils, we were unable to detect a significant change in crearchaeota. Additionally, no detectable change in relative abundance of AOB and NOB were observed in the experimental soils. At the transcript level, the increase in the relative expression of crearchaeota-16S rRNA was generally consistent with trends at the DNA level, that the proportion of crearchaeota decreased with dry-down (from 3 to 1%, and 5 to 1%, respectively), and increased with wet-up (from 1 to 5%, and 1 to 3%, respectively) in the African soils. Similarly, in Chinese and Japanese soils, three nitrifying microbes (crearchaeota, AOB, and NOB), with extremely low transcript abundance, showed no response to changes in water availability.
Discussion
Overall our findings demonstrate that the microbial community of three different forest soils shows high desiccation-tolerance and high resilience to changes in soil water availability. To estimate the state of activity of our soil microorganisms we calculated the rRNA concentration per cell by computing rRNA:rRNA gene ratios (Jones and Lennon, 2010; DeAngelis et al., 2011). We found significantly higher bacterial and archaeal rRNA: rDNA ratios in response to desiccation in our African and Chinese. This suggests an increased accumulation of ribosome concentration when microorganisms face a drought stress. There are organisms that can survive desiccation by entering inactive metabolic states (Vilchez and Manzanera, 2011 #1854; Jones and Lennon, 2010 #1671; Oliver, 2010 #1851), which are associated with the reduced concentrations of nucleic acids (Choi et al., 1996; Dell’Anno et al., 1998; Lebaron et al., 2001; Suzina et al., 2004). Increased ribosome under water stress can provide microorganisms with a higher protein synthesis potential when environmental conditions improve again (Sukenik et al., 2012). In additionally, we found that the change in water availability did not influence the bacterial and archaeal rRNA: rDNA ratios at the forest sites in Japan. Such variability in microbial response to water availability among our study sites could be explained by the effects that different local factors may have on microbial life-strategies under environmental perturbations. For example, the organic carbon and nitrogen content was significantly higher in the Japanese soils and may provide large reservoirs of organic substrate, which can support energy requirements for microorganisms, and prevent resource starvation when facing the desiccation stress. Additionally, a significant proportion of organic carbon might be preferentially mineralized upon rewetting, and thus limit decreases in bacterial and archaeal 16S rRNA genes. Finally, it is possible that the different response patterns of microorganisms to water availability partly depend on changes in total microbial community structure and composition as shown in Figure 4.
The response patterns of bacterial phyla and classes are both site-specific and taxa-specific across our experimental forest soils. For example, Proteobacteria (the most abundant phylum of the present and potentially active community), displayed a water-related opportunistic pattern at the DNA level in soils collected from Chinese and Japanese forest sites. In particular, the relative abundance of Proteobacteria declined with desiccation while increased after rewetting. In contrast, Proteobacteria increased in relative abundance with desiccation, and decreased with rewetting at the DNA level in soils collected from African forest sites. Such variability could be explained by differences in response patterns at the class level whereby Betaproteobacteria were dominant at the Chinese and Japanese sites whereas Alphaproteobacteria were dominant in the African soils. Previous studies show how the abundance of Betaproteobacteria is positively correlated with C mineralization rates (Fierer et al., 2007), which increases with rewetting of dry soils (Birch, 1958; Borken and Matzner, 2009; Inglima et al., 2009). Another important phylum, Actinobacteria, was either stimulated by desiccation or reduced by rewetting at the DNA level in the African and Japanese soils. This agrees with the response of bacterial communities in a California grassland soil study (Barnard et al., 2013). Interestingly, Actinobacteria followed the same pattern at the RNA level in African and Chinese soils. Increases in Actinobacteria during dry-down periods may be attributed to the preparation for next nutrient acquisition. Previous studies have demonstrated that large members of the Gram-positive, high G+C content Actinobacteria phylum were drought resistant and able to grow under extrem dry conditions (Goodfellow and Williams, 1983; Chowdhury et al., 2009). Acidobacteria abundance at the RNA level across Chinese and Japanese soils showed an opposite pattern compared to Actinobacteria. This would agree with clear life-strategy differences between these two groups where Actinobacteria are considered copiotrophic (analog to r-strategist: fast growing and highly variable population size) and Acidobacteria are oligotrophic (analog to k-strategist: slow growing and stable population size). Additionally, beyond potential differences in life-resource strategies our findings clearly suggest that the response of key bacterial phyla and classes to changes in water availability was driven by niche differentiation. Microbial responses to environmental stress may be taxa-specific and dependent on changes in space dimension. Because relationships between rRNA concentration and growth rate can differ significantly among taxa (Barnard et al., 2013), relative rRNA abundance of each phylum may not provide convincing information about which taxa are more active following changes in water availability.
We found that nitrification processes were negligible under dry-down treatments, mainly because nitrifies remain “inactive” under dry soil conditions (Placella and Firestone, 2013). Following water additions, however, we observed high rates of gross N mineralization, which is consistent with previous findings (Birch, 1964; Hungate et al., 1997; Saetre and Stark, 2005). These N mineralization responses were different among all the three forest sites. African forest soils showed large NO3- pulses following rewetting as was found in previous studies (Holloway et al., 1998; Baron et al., 2009), and this increased nitrification was coupled with significant increases in crearchaeota 16S transcripts. On the contrary, little NO3 accumulation was observed in Chinese and Japanese forest soils. No detectable changes were observed either of the abundance of crearchaeota, AOB and NOB 16S transcripts. Such differences in response patterns of nitrifiers among forest soils, might be due to variation in soil gross nitrification rates and soil heterogeneity. The activity of AOA in the nitrification process of African soils confirms findings from the previous studies that AOA might prefer the ammonia released from mineralization (Gubry-Rangin et al., 2010; Zhang et al., 2010). The lower relative abundance of AOA transcripts in Chinese and Japanese soils might relate to negligible nitrification activity during the changes in water availability. The difference in the relative abundance of AOA transcripts might be attributed to variation in soil pH values, with African soils having soil pH (pH 6.0) significantly higher than the others (Chinese soil: pH 4.1; Japanese soil: pH 4.5). Finally, gene and transcript abundances of bacterial and archaeal 16S rRNA and the amoA gene abundance of AOA and AOB also returned to pre-dry-down abundance after water additions.
Despite the different response patterns of microbial taxa to the changes in water availability, our evidence is that upon rewetting, the present and potentially active bacterial community structure as well as the abundance of bacterial (16S), archaeal (16S), and ammonia oxidizers (amoA), all returned to pre-dry-down levels. This resilience ability may greatly contributes to the maintenance of microbial diversity in the forest ecosystem soils under future environmental change.
Author Contributions
XZ and ZJ contributed to the experimental design of the study. XZ collected conducted the analysis of results and wrote the manuscript. All authors contributed to writing and reviewing the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was financially supported by the National Key Basic Research Program of China (No. 2015CB150501), the Strategic Priority Research Program of the CAS (XDB15040000), and the National Science Foundation of China (41090281 and 31270147). We thank Ms. Rong Huang for excellent technical assistance. We greatly appreciate the constructive comments of reviewers.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.01101
References
Barbhuiya, A. R., Arunachalam, A., Pandey, H. N., Arunachalam, K., Khan, M. L., and Nath, P. C. (2004). Dynamics of soil microbial biomass C, N and P in disturbed and undisturbed stands of a tropical wet-evergreen forest. Eur. J. Soil Biol. 40, 113–121. doi: 10.1016/j.ejsobi.2005.02.003
Barnard, R. L., Osborne, C. A., and Firestone, M. K. (2013). Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 7, 2229–2241. doi: 10.1038/ismej.2013.104
Baron, J. S., Schmidt, T. M., and Hartman, M. D. (2009). Climate-induced changes in high elevation stream nitrate dynamics. Glob. Change Biol. 15, 1777–1789. doi: 10.1111/j.1365-2486.2009.01847.x
Basu, S., Joshi, S. K., Pati, D. P., and Behera, N. (1991). Soil respiration in relation to microbial biomass in a tropical deciduous forest floor from India. Rev. D Ecol. Et De Biol. Du Sol. 28, 377–386.
Bell, T., Newman, J. A., Silverman, B. W., Turner, S. L., and Lilley, A. K. (2005). The contribution of species richness and composition to bacterial services. Nature 436, 1157–1160. doi: 10.1038/nature03891
Birch, H. F. (1958). The effect of soil drying on humus decomposition and nitrogen availability. Plant Soil 10, 9–31. doi: 10.1007/BF01343734
Birch, H. F. (1964). Mineralisation of plant nitrogen following alternate wet and dry conditions. Plant Soil. 20, 43–49. doi: 10.1007/BF01378096
Blazewicz, S. J., Barnard, R. L., Daly, R. A., and Firestone, M. K. (2013). Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 7, 2061–2068. doi: 10.1038/ismej.2013.102
Borken, W., and Matzner, E. (2009). Reappraisal of drying and wetting effects on C and N mineralization and fluxes in soils. Glob. Change Biol. 15, 808–824. doi: 10.1111/j.1365-2486.2008.01681.x
Braeken, K., Moris, M., Daniels, R., Vanderleyden, J., and Michiels, J. (2006). New horizons for (p) ppGpp in bacterial and plant physiology. Trends Microbiol. 14, 45–54. doi: 10.1016/j.tim.2005.11.006
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
Castro, H. F., Classen, A. T., Austin, E. E., Norby, R. J., and Schadt, C. W. (2010). Soil microbial community responses to multiple experimental climate change drivers. Appl. Environ. Microbiol. 76, 999–1007. doi: 10.1128/AEM.02874-09
Choi, J. W., Sherr, E. B., and Sherr, B. F. (1996). Relation between presence-absence of a visible nucleoid and metabolic activity in bacterioplankton cells. Limnol. Oceanogr. 41, 1161–1168. doi: 10.4319/lo.1996.41.6.1161
Chowdhury, S. P., Schmid, M., Hartmann, A., and Tripathi, A. K. (2009). Diversity of 16S-rRNA and nifH genes derived from rhizosphere soil and roots of an endemic drought tolerant grass, Lasiurus sindicus. Eur. J. Soil Biol. 45, 114–122. doi: 10.1016/j.ejsobi.2008.06.005
Cole, J. R., Chai, B., Farris, R. J., Wang, Q., Kulam, S. A., McGarrell, D. M., et al. (2005). The ribosomal database project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res. 33, D294–D296. doi: 10.1093/nar/gki038
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–D145. doi: 10.1093/nar/gkn879
Dale, V. H., Joyce, L. A., McNulty, S., Neilson, R. P., Ayres, M. P., Flannigan, M. D., et al. (2001). Climate change and forest disturbances. BioScience 51, 723–734. doi: 10.1641/0006-3568(2001)051[0723:CCAFD]2.0.CO;2
DeAngelis, K. M., Wu, C. H., Beller, H. R., Brodie, E. L., Chakraborty, R., DeSantis, T. Z., et al. (2011). PCR amplification-independent methods for detection of microbial communities by the high-density microarray PhyloChip. Appl. Environ. Microbiol. 77, 6313–6322. doi: 10.1128/AEM.05262-11
Dell’Anno, A., Fabiano, M., Duineveld, G., Kok, A., and Danovaro, R. (1998). Nucleic acid (DNA., RNA) quantification and RNA/DNA ratio determination in marine sediments: comparison of spectrophotometric, fluorometric, and highperformance liquid chromatography methods and estimation of detrital DNA. Appl. Environ. Microbiol. 64, 3238–3245.
Ding, L.-J., Su, J.-Q., Xu, H.-J., Jia, Z.-J., and Zhu, Y.-G. (2015). Long-term nitrogen fertilization of paddy soil shifts iron-reducing microbial community revealed by RNA-13C-acetate probing coupled with pyrosequencing. ISME J. 9, 721–734. doi: 10.1038/ismej.2014.159
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Fierer, N., and Schimel, J. P. (2002). Effects of drying-rewetting frequency on soil carbon and nitrogen transformations. Soil Biol. Biochem. 34, 777–787. doi: 10.1016/s0038-0717(02)00007-x
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. U.S.A. 102, 14683–14688. doi: 10.1073/pnas.0506625102
Franzluebbers, A. J. (1999). Potential C and N mineralization and microbial biomass from intact and increasingly disturbed soils of varying texture. Soil Biol. Biochem. 31, 1083–1090. doi: 10.1016/s0038-0717(99)00022-x
Goodfellow, M., and Williams, S. (1983). Ecology of actinomycetes. Ann. Rev. Microbiol. 37, 189–216. doi: 10.1146/annurev.mi.37.100183.001201
Griffiths, R. I., Whiteley, A. S., O’Donnell, A. G., and Bailey, M. J. (2000). Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA-and rRNA-based microbial community composition. Appl. Environ. Microbiol. 66, 5488–5491. doi: 10.1128/AEM.66.12.5488-5491.2000
Griffiths, R. I., Whiteley, A. S., O’Donnell, A. G., and Bailey, M. J. (2003). Physiological and community responses of established grassland bacterial populations to water stress. Appl. Environ. Microbiol. 69, 6961–6968. doi: 10.1128/aem.69.12.6961-6968.2003
Gubry-Rangin, C., Nicol, G. W., and Prosser, J. I. (2010). Archaea rather than bacteria control nitrification in two agricultural acidic soils. FEMS Microbiol. Ecol. 74, 566–574. doi: 10.1111/j.1574-6941.2010.00971.x
Hashsham, S. A., Fernandez, A. S., Dollhopf, S. L., Dazzo, F. B., Hickey, R. F., Tiedje, J. M., et al. (2000). Parallel processing of substrate correlates with greater functional stability in methanogenic bioreactor communities perturbed by glucose. Appl. Environ. Microbiol. 66, 4050–4057. doi: 10.1128/AEM.66.9.4050-4057.2000
Holloway, J., Dahlgren, R., Hansen, B., and Casey, W. (1998). Contribution of bedrock nitrogen to high nitrate concentrations in stream water. Nature 395, 785–788. doi: 10.1038/27410
Hungate, B. A., Lund, C. P., Pearson, H. L., and Chapin, I. I. I. F. S. (1997). Elevated CO2 and nutrient addition after soil N cycling and N trace gas fluxes with early season wet-up in a California annual grassland. Biogeochemistry. 37, 89–109. doi: 10.1023/A:1005747123463
Inglima, I., Alberti, G., Bertolini, T., Vaccari, F. P., Gioli, B., Miglietta, F., et al. (2009). Precipitation pulses enhance respiration of Mediterranean ecosystems: the balance between organic and inorganic components of increased soil CO2 efflux. Glob. Change Biol. 15, 1289–1301. doi: 10.1111/j.1365-2486.2008.01793.x
Jia, C. Y., and Insam, H. (1991). Microbial biomass and relative contributions of bacteria and fungi in soil beneath tropical rain-forest, hainan island, china. J. Trop. Ecol. 7, 385–393. doi: 10.1017/S0266467400005678
Jones, S. E., and Lennon, J. T. (2010). Dormancy contributes to the maintenance of microbial diversity. Proc. Natl. Acad. Sci. U.S.A. 107, 5881–5886. doi: 10.1073/pnas.0912765107
Kemnitz, D., Kolb, S., and Conrad, R. (2005). Phenotypic characterization of rice cluster III archaea without prior isolation by applying quantitative polymerase chain reaction to an enrichment culture. Environ. Microbiol. 7, 553–565. doi: 10.1111/j.1462-2920.2005.00723.x
Kemp, P., Lee, S., and LaRoche, J. (1993). Estimating the growth rate of slowly growing marine bacteria from RNA content. Appl. Environ. Microbiol. 59, 2594–2601.
Kerkhof, L., and Ward, B. B. (1993). Comparison of nucleic acid hybridization and fluorometry for measurement of the relationship between RNA/DNA ratio and growth rate in a marine bacterium. Appl. Environ. Microbiol. 59, 1303–1309.
Lebaron, P., Servais, P., Agogué, H., Courties, C., and Joux, F. (2001). Does the high nucleic acid content of individual bacterial cells allow us to discriminate between active cells and inactive cells in aquatic systems? Appl. Environ. Microbiol. 67, 1775–1782. doi: 10.1128/AEM.67.4.1775-1782.2001
Luizao, F. J., Proctor, J., Thompson, J., Luizao, R. C. C., Marrs, R. H., Scott, D. A., et al. (1998). Rain forest on maraca island., roraima., brazil: soil and litter process response to artificial gaps. For. Ecol. Manage 102, 291–303. doi: 10.1016/S0378-1127(97)00170-9
Malhi, Y., Roberts, J. T., Betts, R. A., Killeen, T. J., Li, W., and Nobre, C. A. (2008). Climate change., deforestation., and the fate of the amazon. Science 319, 169–172. doi: 10.1126/science.1146961
Mikha, M. M., Rice, C. W., and Milliken, G. A. (2005). Carbon and nitrogen mineralization as affected by drying and wetting cycles. Soil Biol. Biochem. 37, 339–347. doi: 10.1016/j.soilbio.2004.08.003
Miller, A. E., Schimel, J. P., Meixner, T., Sickman, J. O., and Melack, J. M. (2005). Episodic rewetting enhances carbon and nitrogen release from chaparral soils. Soil Biol. Biochem. 37, 2195–2204. doi: 10.1016/j.soilbio.2005.03.021
Muttray, A. F., Yu, Z., and Mohn, W. W. (2001). Population dynamics and metabolic activity of Pseudomonas abietaniphila BKME-9 within pulp mill wastewater microbial communities assayed by competitive PCR and RT-PCR. FEMS Microbiol. Ecol. 38, 21–31. doi: 10.1111/j.1574-6941.2001.tb00878.x
Oliver, J. D. (2010). Recent findings on the viable but nonculturable state in pathogenic bacteria. FEMS Microbiol. Rev. 34, 415–425. doi: 10.1111/j.1574-6976.2009.00200.x
Placella, S. A., and Firestone, M. K. (2013). Transcriptional response of nitrifying communities to wetting of dry soil. Appl. Environ. Microbiol. 79, 3294–3302. doi: 10.1128/AEM.00404-13
Poulsen, L. K., Ballard, G., and Stahl, D. A. (1993). Use of rRNA fluorescence in situ hybridization for measuring the activity of single cells in young and established biofilms. Appl. Environ. Microbiol. 59, 1354–1360.
Rotthauwe, J. H., Witzel, K. P., and Liesack, W. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712.
Saetre, P., and Stark, J. M. (2005). Microbial dynamics and carbon and nitrogen cycling following re-wetting of soils beneath two semi-arid plant species. Oecologia 142, 247–260. doi: 10.1007/s00442-004-1718-9
Schippers, A., Neretin, L. N., Kallmeyer, J., Ferdelman, T. G., Cragg, B. A., Parkes, R. J., et al. (2005). Prokaryotic cells of the deep sub-seafloor biosphere identified as living bacteria. Nature 433, 861–864. doi: 10.1038/nature03302
Shi, Y., Xiang, X., Shen, C., Chu, H., Neufeld, J. D., Walker, V. K., et al. (2015). Vegetation-Associated impacts on arctic tundra bacterial and microeukaryotic communities. Appl. Environ. Microbiol. 81, 492–501. doi: 10.1128/AEM.03229-14
Singh, J. S., Raghubanshi, A. S., Singh, R. S., and Srivastava, S. C. (1989). Microbial biomass acts as a source of plant nutrients in dry tropical forest and savanna. Nature 338, 499–500. doi: 10.1038/338499a0
Stubner, S. (2002). Enumeration of 16S rDNA of Desulfotomaculum lineage 1 in rice field soil by real-time PCR with SybrGreen (TM) detection. J. Microbiol. Methods 50, 155–164. doi: 10.1016/S0167-7012(02)00024-6
Sukenik, A., Hadas, O., Kaplan, A., and Quesada, A. (2012). Invasion of Nostocales (cyanobacteria) to subtropical and temperate freshwater lakes–physiological, regional, and global driving forces. Front. Microbiol. 3:86. doi: 10.3389/fmicb.2012.00086
Suzina, N., Mulyukin, A., Kozlova, A., Shorokhova, A., Dmitriev, V., Barinova, E., et al. (2004). Ultrastructure of resting cells of some non-spore-forming bacteria. Microbiology 73, 435–447. doi: 10.1023/B:MICI.0000036990.94039.af
Vilchez, S., and Manzanera, M. (2011). Biotechnological uses of desiccation-tolerant microorganisms for the rhizoremediation of soils subjected to seasonal drought. Appl. Microbiol. Biotechnol. 91, 1297–1304. doi: 10.1007/s00253-011-3461-6
Westerling, A. L., Hidalgo, H. G., Cayan, D. R., and Swetnam, T. W. (2006). Warming and earlier spring increase western U.S. Forest Wildfire Activity. Science 313, 940–943. doi: 10.1126/science.1128834
Xia, W., Zhang, C., Zeng, X., Feng, Y., Weng, J., Lin, X., et al. (2011). Autotrophic growth of nitrifying community in an agricultural soil. ISME J. 5, 1226–1236. doi: 10.1038/ismej.2011.5
Xu, M., and Qi, Y. (2001). Soil-surface CO2 efflux and its spatial and temporal variations in a young ponderosa pine plantation in northern California. Glob. Change Biol. 7, 667–677. doi: 10.1046/j.1354-1013.2001.00435.x
Zhang, L.-M., Offre, P. R., He, J.-Z., Verhamme, D. T., Nicol, G. W., and Prosser, J. I. (2010). Autotrophic ammonia oxidation by soil thaumarchaea. Proc. Natl. Acad. Sci. U.S.A. 107, 17240–17245. doi: 10.1073/pnas.1004947107
Keywords: rRNA, rDNA, amoA, pyrosequencing, qPCR, forest, seed bank
Citation: Zhou X, Fornara D, Ikenaga M, Akagi I, Zhang R and Jia Z (2016) The Resilience of Microbial Community under Drying and Rewetting Cycles of Three Forest Soils. Front. Microbiol. 7:1101. doi: 10.3389/fmicb.2016.01101
Received: 02 June 2016; Accepted: 01 July 2016;
Published: 19 July 2016.
Edited by:
Hongchen Jiang, Miami University, USAReviewed by:
Qiuyuan Huang, Miami University, USAYe Deng, Research Center for Eco-Environmental Sciences, China
Copyright © 2016 Zhou, Fornara, Ikenaga, Akagi, Zhang and Jia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhongjun Jia, jia@issas.ac.cn