Science of The Total Environment: Yaqi Sheng, Lizhong Zhu
Science of The Total Environment: Yaqi Sheng, Lizhong Zhu
Science of The Total Environment: Yaqi Sheng, Lizhong Zhu
H I G H L I G H T S G R A P H I C A L A B S T R A C T
a r t i c l e i n f o a b s t r a c t
Article history: Biochar application to soil has been proposed for soil carbon sequestration and global warming mitigation. While
Received 5 October 2017 recent studies have demonstrated that soil pH was a main factor affecting soil microbial community and stability
Received in revised form 29 November 2017 of biochar, little information is available for the microbiome across different soil pH and the subsequently CO2
Accepted 29 November 2017
emission. To investigate soil microbial response and CO2 emission of biochar across different pH levels, compar-
Available online xxxx
ative incubation studies on CO2 emission, degradation of biochar, and microbial communities in a ferralsol
Editor: J Jay Gan (pH 5.19) and a phaeozems (pH 7.81) with 4 biochar addition rates (0.5%, 1.0%, 2.0%, 5.0%) were conducted. Bio-
char induced higher CO2 emission in acidic ferralsol, largely due to the higher biochar degradation, while the
Keywords: more drastic negative priming effect (PE) of SOC resulted in decreased total CO2 emission in alkaline phaeozems.
Soil carbon sequestration The higher bacteria diversity, especially the enrichment of copiotrophic bacteria such as Bacteroidetes,
Biochar Gemmatimonadetes, and decrease of oligotrophic bacteria such as Acidobacteria, were responsible for the in-
pH creased CO2 emission and initial positive PE of SOC in ferralsol, whereas biochar did not change the relative abun-
SOC dances of most bacteria at phylum level in phaeozems. The relative abundances of other bacterial taxa (i.e.
Microorganisms
Actinobacteria, Anaerolineae) known to degrade aromatic compounds were also elevated in both soils. Soil pH
was considered to be the dominant factor to affect CO2 emission by increasing the bioavailability of organic car-
bon and abundance of copiotrophic bacteria after biochar addition in ferralsol. However, the decreased bioavail-
ability of SOC via adsorption of biochar resulted in higher abundance of oligotrophic bacteria in phaeozems,
leading to the decrease in CO2 emission.
© 2017 Published by Elsevier B.V.
⁎ Corresponding author at: Department of Environmental Science, Zhejiang University, Hangzhou, Zhejiang 310058, China.
E-mail address: [email protected] (L. Zhu).
https://doi.org/10.1016/j.scitotenv.2017.11.337
0048-9697/© 2017 Published by Elsevier B.V.
1392 Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399
Table 1
Soil physical and chemical properties of the biochar-amended treatments and control after 90 days incubation.
pH NH+
4 -N NO−
3 -N TC WEOCc
(mg/kg) (mg/kg) (mg/g) (mg/L)
F 5.996 ± 0.015aa 0.639 ± 0.003b 19.65 ± 0.43b 11.39 ± 0.31a 28.82 ± 0.39b
FB 0.5 5.933 ± 0.018a 0.511 ± 0.007b 4.034 ± 0.38a 13.74 ± 0.14a 26.63 ± 0.26ab
FB 1.0 6.188 ± 0.016ab 0.286 ± 0.004ab 4.289 ± 0.23a 16.31 ± 0.26ab 25.45 ± 0.29ab
FB 2.0 6.207 ± 0.025ab 0.028 ± 0.001a 6.797 ± 0.17a 20.92 ± 0.75ab 24.11 ± 0.28ab
FB 5.0 6.758 ± 0.014b /b 7.817 ± 0.76a 34.76 ± 0.79b 21.67 ± 0.58a
P 8.237 ± 0.001c 0.055 ± 0.002a 6.247 ± 0.53a 37.68 ± 0.66bc 33.42 ± 0.07b
PB 0.5 8.274 ± 0.002c 0.044 ± 0.001a 5.253 ± 0.06a 39.02 ± 0.64bc 32.09 ± 0.04b
PB 1.0 8.325 ± 0.002c 0.040 ± 0.001a 4.762 ± 0.87a 40.67 ± 0.43bc 29.49 ± 0.46b
PB 2.0 8.334 ± 0.003c 0.034 ± 0.001a 4.725 ± 0.12a 45.19 ± 0.83c 27.41 ± 0.16ab
PB 5.0 8.407 ± 0.004c 0.021 ± 0.001a 4.362 ± 0.44a 56.46 ± 1.32c 26.80 ± 0.02ab
a
Values in the same column followed by different letters were statistically different at p b 0.05, n = 3.
b
Blow limit of detection.
c
Water extracted organic carbon.
was set with 6 replicates. Deionized water of 80% soil water holding ca- normalized, pooled in equimolar and paired-end sequenced (2
pacity was added to the vials, the vials were sealed and incubated in the × 300 bp) on an Illumina MiSeq platform (Guhe Information Technolo-
dark at 25 °C. Three empty vials were set to measure the concentration gy Co. Ltd., Hangzhou, China). Sequence reads processing was per-
of CO2 of surrounding atmosphere. The measurement of CO2 was de- formed using QIIME (version 1.9.0), including additional quality
tailed before (Sheng et al., 2016). The CO2 emission was measured trimming and demultiplexing. Operational taxonomic unit (OTU) was
with a gas-tight glass syringe (Agilent) through injecting 100 μL head- picked by Vsearch v1.11.1, followed by dereplication, cluster, detection
space CO2 from the incubation vials into the gas chromatograph (Fuli of chimeras. The barcoded 16S rRNA gene sequences were annotated
9790, Hangzhou, China). Limits of detection (LOD) of CO2 measurement into libraries by QIIME, resulting in operational taxonomic units
was 2 ppm, the relative standard deviation (RSD) were lower than 5%. (OTUs), bacterial diversity, and phenetic classification. Taxonomy was
The incubation experiments was carried over a period of 90 days, and assigned to an OTU based on its representative sequence by Greengene
CO2 evolution was measured periodically at 1, 2, 4, 7, 14, 21, 30, 40, 13.8. Calculations of alpha-diversity (including Shannon, Simpson,
50,60, 75, 90 days of incubation. Chao1) were performed by QIIME.
Phylogenetic investigation of PICRUSt was used to predict the Kyoto
2.3. Degradation of biochar and priming effects (PE) of soil organic carbon Encyclopedia of Genes and Genomes (KEGG). PICRUSt predicted
(SOC) metagenomes from 16S data by using evolutionary modeling and com-
paring with a reference genome database. Generally, OTUs of 16S rRNA
δ13C analysis by Precon-IRMS (Delta V Thermo Finnigan, Waltham, sequences were normalized by PICRUSt, and then the metagenomes
America) was used to calculate the amount of biochar mineralized predicted by PICRUSt algorithm were collapsed into clusters of
CO2 (CB) from biochar-amended soil, methods has been detailed before orthologous groups of proteins (COGs) and KEGG. Predictive COGs and
(Sheng et al., 2016). CB was calculated as: KEGGs were screened out and visualized by Statistical Analysis of
Metagenomic Profiles (STAMP) software package v 2.1.3. The KEGG
C T δT 13 CO2 −δS 13 CO2 profiles across all samples were pairwise compared between control
CB ¼ ð1Þ and biochar-amended treatments.
δB 13 CO2 −δS 13 CO2
2.5. Statistical analyses
CT is the total C mineralized from the biochar-amended soil, δ13T CO2
is the δ13T C value of the CO2-C evolved from the biochar-amended The regression and correlation analyses were performed on SPSS
soils, δ13 13
S CO2 is the δT C value of CO2-C evolved from the control soil, 20.0 (IBM, Armonk, USA). The differences in physiochemical properties
and δ13S CO 2 is the initial δ13
T C value of the biochar. The amount of the and relative abundances of dominant communities were tested by the
soil-C mineralized (CS,B) from the biochar-amended soil and the prim- analysis of variance (ANOVA). Kyoto Encyclopedia of Genes and Ge-
ing effects (PE) induced by the biochar on mineralization of the native nomes (KEGG) pathways were drawn by using STAMP (version 2.1.3).
organic carbon in the soils were calculated as: Pearson correlation analysis was used to assess the relationships
among physiochemical properties of soils, CO2 emission and relative
C S;B ¼ C T −C B ð2Þ abundances of dominant communities. The relationship between mi-
crobial community structure and the selected environmental variables
PE ¼ C S;B −C S ð3Þ was analyzed by the redundancy analysis (RDA) to determine the inde-
pendent contributions of these selected variables to the variation in
CS is the amount of CO2-C evolved from the control soil, CT is the total community composition.
C mineralized from the biochar-amended soil, CB is the amount of C
mineralized from biochar. The positive priming effect means accelerat- 3. Results
ing of SOC mineralization by biochar. In contrast, negative priming effect
means the suppression of SOC mineralization by biochar. 3.1. Soil chemical properties
2.4. High-throughput sequencing analysis Generally, biochar addition increased soil pH in both soils except
ferralsol amended with 0.5% biochar, with the greatest increase of
Metagenomic DNA was extracted from 500 mg freeze-dried soil 0.825 and 0.16 at the rate of 5.0% for ferralsol and phaeozems, respec-
Samples using the PowerMax (stool/soil) DNA isolation kit (MoBio Lab- tively (Table 1). The increase of pH in ferralsol was more drastic than
oratories). The V3-V4 hypervariable region of the 16S rRNA gene was that in phaeozems because of the higher buffer ability of phaeozems
amplified from genomic DNA by primers. Purified amplicons were whose initial content of total carbon was higher (Murray et al., 2015).
1394 Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399
−
The concentration of NH+ 4 and NO3 was highest in untreated ferralsol, More biochar C was mineralized in ferralsol (Fig. 3), the highest min-
NH+ 4 content decreased significantly (p b 0.05) as the biochar addition eralization rate of biochar appeared on the first day. Biochar degrada-
rate increased. NO− 3 content drastically decreased from 19.65 mg/kg in tion followed an order of FB 5.0 N FB 2.0 N FB 1.0 N FB 0.5 over the
untreated ferralsol to 4.03 mg/kg in FB 0.5, then arose as addition rate whole 90 days of incubation in ferralsol. However, degradation rate of
increased. In contrast, only a slightly decrease in the concentration of biochar showed an opposite trend on the first day in phaeozems, the
−
NH+ 4 and NO3 were observed after biochar application in phaeozems. degradation rate of biochar began to increase at higher addition rate
The total carbon (TC) content and water extractable organic carbon over time. Following sharp declines after the transient high degradation
(WEOC) was higher in phaeozems (Table 1). Total carbon (TC) content rate of biochar on the first day, second peak of biochar degradation was
was 37.68 mg/g in phaeozems, which was almost 2.31 folders higher observed on the 4-7th day of the incubation except the addition rate of
than ferralsol. Biochar significantly increased TC content from 11.39 to 5.0%. After the 30th day of the incubation, the degradation rate of bio-
34.76 mg/g and 37.68 to 56.46 mg/g for ferralsol and phaeozems at char decreased to a relative stable phase for both soils.
the rate of 5.0%. While WEOC decreased as the biochar addition rate in- Biochar addition initially accelerated the SOC mineralization in acid-
creased in both soils, biochar significantly (p b 0.05) decreased WEOC by ic ferralsol, this positive priming effect (PE) was facilitated at higher bio-
7.6%, 11.7%, 16.3% and 24.8% in ferralsol at the rate of 0.5%, 1.0%, 2.0% and char addition rate (Fig. 4a). The highest cumulative PE was 4.33, 7.41,
5.0%, respectively. For phaeozems, biochar decreased WEOC by 4.0%, 16.71 and 26.66 μg for 0.5%, 1.0% 2.0% and 5.0% biochar- amended
11.7%, 17.9% and 19.8% at the rate of 0.5%, 1.0%, 2.0% and 5.0%. ferralsol, respectively. The positive PE were transient and shifted to neg-
Biochar amendment altered the molecular-level composition and ative PE as time went by except ferralsol amended with 5.0% addition
distribution of native soil organic matter (Fig. 1). The proportion of rate of biochar. Negative cumulative PE were observed on the 4th, 7th
most aliphatic hydrocarbon and carbohydrates by solvent extraction and 20th day of the incubation for ferralsol amended with 0.5%, 1.0%
was higher in phaeozems (Tables S2, S3), resulted in higher bioavail- and 2.0% biochar. Compared with ferralsol, the treatments of
ability. After biochar amendment, the proportion of most aliphatic hy- phaeozems exhibited negative PE through the whole incubation. Soil
drocarbon and carbohydrates such as palmitic acid, lignoceric acid, amended with higher biochar addition rate (5.0%) led to a 2–3 fold
stearic acid, trehalose and glycerol monostearate significantly increased higher negative PE than the other treatments.
(p b 0.05) in ferralsol at the end of the incubation (Table S2). However,
as the addition rate increased, palmitic acid, stearic acid and trehalose 3.3. Shift of microbial community composition under biochar application
significantly decreased in phaeozems (Table S3). For example, the pro-
portion of palmitic acid decreased from 7.093% to 5.565% and 5.494% at Biochar amendment showed higher OTU richness and bacterial di-
the addition rate of 1.0% and 5.0%, respectively. More aromatic com- versity in acidic ferralsol, but showed slightly decrease in phaeozems.
pounds such as 2, 6-Bis (tert-butyl) phenol were shown in ferralsol The untreated phaeozems exhibited 2.38 folders higher OTUs than un-
(Table S4), aromatic compounds (i.e. benzenepropanoic acid, phthalic treated acidic ferralsol (Table 2). Changes in both bacterial community
acid and quinoline) significantly increased as the biochar addition rate structure and composition were found under different biochar addition
arose. rate by high-throughput sequencing. The three replicates of each treat-
ment were generally assembled together and distinct from each other
(Fig. S1). Furthermore, comparison of the sites effect prior to biochar
3.2. CO2 emission, degradation of biochar and priming effect (PE) of SOC treatment showed significant differences in community structure, em-
phasizing the importance of original soil properties on microorganism.
The cumulative amount of C mineralized from the biochar was sig- The predominant phyla in soils were Proteobacteria,
nificantly higher in ferralsol than that in phaeozems. Biochar signifi- Gemmatimonadetes, Firmicutes, Chloroflexi, Bacteroidetes, Actinobacteria,
cantly increased the total CO2 emission in the acidic ferralsol at all Acidobacteria, Cyanobacteria, Nitrospirae and TM7. These taxa accounted
addition rates (p b 0.05), the total C mineralized was significantly higher for N 89% of the bacterial sequences in all treatments (Table S5). As
in the 5.0% biochar treatment than other treatments (Fig. 2). At the end shown in Fig. 5, the relative abundance of Firmicutes decreased under
of the incubation, CO2 emission in ferralsol increased from 80.36 mg to biochar treatments in ferralsol, Acidobacteria decreased except the in
83.16, 84.92, 91.43 and 108.54 mg at the rate of 0.5%, 1.0%, 2.0% and the FB0.5, while Chloroflexi, Gemmatimonadetes, Bacteroidetes and
5.0%, respectively. In contrast, the total CO2 emission was suppressed Actinobacteria increased. Acidobacteria was the most sensitive phylum
in the alkaline phaeozems after biochar application, as the addition to biochar in ferralsol, with the relative abundances decreasing from
rate increase, the suppression of CO2 emission increased. At the end of 10.43% to 9.50%, 7.52% and 6.90% for FB 1.0, FB 2.0 and FB 5.0, respective-
the incubation, biochar decreased CO2 emission by 3.85%, 5.30%, 6.34% ly. Firmicutes and TM7 were also significantly decreased by biochar
and 14.9% at the rate of 0.5%, 1.0%, 2.0% and 5.0%. amendment, from 27.24% and 10.42% to 17.76% and 2.13% for 5.0%
14
(a) F
14 (b) P
12 F-1 12 P-1
F-5 P-5
10 10
8 8
Percent
Percent
6 6
4 4
2 2
0 0
0 100 200 300 400 500 600 700 800 900 0 100 200 300 400 500 600 700 800 900
Molecular weight Molecular weight
Fig. 1. Distribution of small molecular compounds (molecular weight b 1000) extracted by solvent in ferralsol (a) and phaeozems (b) amended with biochar across different addition rate
at the end of the incubation.
Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399 1395
60 60
40 40
20 20
0 0
0 10 20 30 40 50 60 70 80 90 0 10 20 30 40 50 60 70 80 90
Fig. 2. Total CO2 release (mg C kg−1 soil) in (a) ferralsol, (b) phaeozems across 0.5%, 1.0%, 2.0% and 5.0% biochar addition rate over 90 days of incubation. The ferralsol amend with 0%, 0.5%,
1.0%, 2.0% and 5.0% biochar were annotated as F, FB 0.5, FB 1.0, FB 2.0 and FB 5.0, the phaeozems amend with 0%, 0.5%, 1.0%, 2.0% and 5.0% biochar were annotated as P, PB 0.5, PB 1.0, PB 2.0
and PB 5.0.
addition rate, respectively (Table S5). On the other hand, no significant CO2 emission incorporated degradation of biochar and mineralization
differences were observed for biochar amended phaeozems except for of SOC, the effect of biochar on the total emission depended on the
Proteobacteria, which showed a slightly decrease in comparison to the sum of biochar mineralization and the PE of SOC induced by biochar.
control except PB 1.0. In addition, Actinobacteria and Chloroflexi also sig- Biochar degradation and the initial positive PE of SOC contributed to
nificantly increased for PB 5.0. the increased CO2 emission in acidic ferralsol. Over time, although the
The differences of diversity and richness of soil bacterial communi- positive PE shifted to negative except FB5.0, biochar degradation was
ties could primarily be explained by significant increase of pH in acidic larger and played a more important role. In contrast, the negative PE
ferralsol, whereas in phaeozems, TC is dominated factor for soil bacterial of SOC induced by biochar was higher than its own degradation in alka-
communities after biochar amendment. The RDA result indicated that line phaeozems.
biochar could strongly shift the soil microbial community by changing The cumulative amount of C mineralized from the biochar was sig-
soil chemical properties. The first two axes of redundancy analysis nificantly higher in soil with lower pH. The CO2 release from the carbon-
(RDA) together explain 47.6% and 34.9% of the total variation in the bac- ates of biochar was largely depended on soil pH. The transient high
terial communities (at the family level) of ferralsol and phaeozems, re- degradation rate of biochar followed by a sharp decline in acidic soil in-
−
spectively. In total, pH, TC, NH+4 , NO3 and WEOC could explain 53.8% dicated a faster abiotic release of carbonates. Besides, since acidic
and 46.1% of the soil microbial community composition variation for ferralsol exhibited lower proportion of aliphatic hydrocarbon and
ferralsol and phaeozems, respectively (Table S8). The importance of higher aromatic compounds, the bioavailability of SOC would be
variables on microorganism were followed as pH N TC N WEOC N NH+ 4 lower. Previous researches also showed the presence of significant
N NO− 3 in ferralsol, whereas in phaeozems TC and pH could explain char like aryl and O-aryl C of SOC in acidic soil (Farrell et al., 2013). As
10.4% and 7.40% of the variance in soil bacterial community composition a result, the labile fraction in biochar induced higher proportion of ali-
(Fig. 6). phatic hydrocarbon and carbohydrates in ferralsol, thus to increase
the bioavailability of sample and change the preferential substrate utili-
4. Discussion zation by soil microbiome, leading to higher microbial diversity and bio-
char degradation in ferralsol. Besides, the microbiome degrading char-
4.1. Higher CO2 emission, degradation of biochar and priming effect (PE) of like substance in soil with low pH would be more adaptive in the min-
SOC in acidic ferralsol eralization of recalcitrant fraction in biochar, also contributed to the
long-term biochar degradation. In contrast, the adsorption of WEOC
As the addition rate increased, biochar induced higher CO2 emission on biochar decreased the concentration of easily degradable substances
in acidic ferralsol, largely due to the higher biochar degradation, while such as stearic acid and trehalose in phaeozems, resulted in lower bio-
the more drastic negative priming effect (PE) of SOC resulted in de- availability of SOC, thus to decreased the microbial activity and degrada-
creased total CO2 emission in alkaline phaeozems. Generally, the total tion of biochar.
3.5 3.5
Respiration rate (mg CO2-C kg soil d )
-1
(a) (b)
3.0 FB0.5 3.0 PB0.5
-1
-1
FB1.0 PB1.0
2.5 FB2.0 2.5 PB2.0
FB5.0 PB5.0
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0 0.0
-10 0 10 20 30 40 50 60 70 80 90 -10 0 10 20 30 40 50 60 70 80 90
Incubation time (days) Incubation time (days)
Fig. 3. Respiration rate of biochar (mg CO2-C kg−1 soil) in (a) ferralsol, (b) phaeozems across 0.5%, 1.0%, 2.0% and 5.0% biochar addition rate over 90 days of incubation.
1396 Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399
120
(a) FB0.5 PB0.5
(b)
FB1.0 PB1.0
80
FB2.0 PB2.0
FB5.0 PB5.0
Primed soil-C ( g)
40
-40
-80
-120
-160
-10 0 10 20 30 40 50 60 70 80 90
Fig. 4. Primed soil organic carbon of ferralsol and phaeozems across 0.5%, 1.0%, 2.0% and 5.0% biochar addition rate over 90 days of incubation (a). Positive primed soil organic carbon (SOC)
indicated higher carbon mineralization of SOC in biochar-amended soil; negative primed SOC indicated decreased SOC mineralization. Silica-alumina mineral inserted to the pore space of
biochar through organo-mineral interactions in ferralsol at 5.0% biochar addition rate at the end of the incubation by scanning electron microscope (b).
Biochar addition initially accelerated the SOC mineralization in acid- 2012b) and stabilize SOC compared to relatively fresh biochar (Luo et
ic ferralsol, but induced drastic negative PE in alkaline phaeozems al., 2011). The occlusion of clay was indicated by in Fig. 4b. The silica-
through the whole incubation. Liming effect and co-metabolism was alumina mineral gradually inserted into the pore space of biochar at
supposed to be responsible for the initial positive PE. Since biochar sig- the 90th day of incubation, which would occur via adsorption (Singh
nificantly increased the pH of ferralsol, the growth of microbiome et al., 2012) and ligand exchange reactions (Joseph et al., 2010). The
would be greatly facilitated, thus to accelerate the mineralization of organo-mineral interaction was more important in the long-term soil
SOC. Whereas in phaeozems, the pH only met a slight change. The easily C sequestration by biochar.
mineralized labile fraction in biochar resulted in increased microbial In conclusion, the degradation of biochar was responsible for the en-
biomass and the simultaneous enzyme production (Kuzyakov et al., hanced CO2 emission in ferralsol, while the adsorption of native SOC at-
2009; Maestrini et al., 2015). The decreased cumulative positive priming tributed to the emission reduction of CO2 in phaeozems. pH was the
effect and correspondent decreased biochar degradation at the 7th day of main parameter in determining the degradation of biochar and PE of
incubation also indicated the co-metabolism. Although co-metabolism SOC through carbonates dissolution, affecting organo-mineral interac-
also existed in phaeozems, the decreased WEOC through adsorption tions, as well as changing the composition of SOC and microbiome.
and lower aliphatic hydrocarbon in biochar-amended phaeozems result- Since the key role of microbial community in regulating dynamics of
ed in the reduction in the accessibility of C by the microbes and enzymes CO2 emissions, studies on microbial community change under two dif-
(Lehmann et al., 2005; Lehmann et al., 2011; Zimmerman et al., 2011). ferent pH levels and its relation with CO2 emission should be conducted.
The toxicity effects of some compounds such as phthalic acid and PAHs
in biochar (Qian and Chen, 2014; Spokas, 2013) also led to the negative 4.2. pH induced microbial community composition change
PE by reducing soil microbial activity in phaeozems.
Due to the depletion of labile fraction in biochar and the stabilization The increase in bacterial diversity after biochar amendment could be
of SOC by organo-mineral interactions, the PE shifted from positive to prior attributed to the alkalinity of biochar (liming effect) in ferralsol. As
negative in ferralsol over time. This organo-mineral interaction was
also affected by soil pH. As previous studies indicated, complexation of TM7 Chloroflexi Nitrospirae Firmicutes
carboxyl/hydroxyl functional groups of organic matter and iron oxide Gemmatimonadetes Bacteroidetes Proteobacteria
surfaces via ligand exchange increases with decreasing pH (Gu et al., Cyanobacteria Actinobacteria Acidobacteria
1.0
1994; von Lutzow et al., 2006). This was proved by the earlier shift to
negative PE of SOC at lower addition rate. Furthermore, the aging of bio-
char developed negative charge through oxidation reactions would in- 0.8
crease sorption of SOC to biochar particles (Cao et al., 2009; Lin et al.,
Relative abundance
Table 2 0.6
Alpha diversity per sample of bacterial community at the end of the incubation.
1.0
1.0
(a) NH4+ (b) WEOC
FB0.5-1 NH4+ P-1
RDA2 (12.8%)
RDA2 (7.2%)
F-3
F-2 PB0.5-3
FB1.0-2 PB1.0-2
FB1.0-1 NO3-
FB1.0-3 PB1.0-1
PB2.0-1
FB5.0-1
PB5.0-1
FB2.0-3 FB5.0-2 PB2.0-2
PB2.0-3
FB2.0-1 PB5.0-2
pH
PB1.0-3 PB5.0-3
FB2.0-2 TC FB5.0-3
pH TC
-1.0
-1.0
the bacterial diversity was the highest in neutral soils and lower in acid- to degrade recalcitrant polymers (Kirby, 2005). Previous studies also in-
ic soils (Fierer and Jackson, 2006), soil with higher biochar addition rate dicated biochar amendment was characterized by an increase in the rel-
led to higher soil pH and bacterial diversity. What's more, since a more ative abundance of Actinobacteria (Kolton et al., 2011; O'Neill et al.,
gram-negative bacteria-dominated community in neural and alkaline 2009). The relative abundances of Anaerolineae, which were assumed
soils (Feng et al., 2010) was proposed, copiotrophic phylum like to be a potential degrader of relatively recalcitrant C compounds such
Chloroflexi, Gemmatimonadetes, Bacteroidetes significantly increased as phenol (Zheng et al., 2016), met significantly increases in biochar
after biochar amendment. In addition, biochar also led to a significant treatments. At the family level, Azoarcus group which could grow with
decline in abundance of Acidobacteria and TM7, which was reported to aromatic compounds as the only source and Sphingomonas also in-
dominate in low pH soils (Hanada et al., 2014; Verastegui et al., 2014). creased in phaeozems (Rabus, 2005; White et al., 1996). Furthermore,
The metabolic quotient and enzymatic activity would also be affected those acromatic degrading bacteria would directly affect the utilization
by the pH, thus to affect the utilization of the substrates and CO2 emis- of recalcitrant C in biochar, thus to alter the PE of SOC the total CO2 emis-
sion. Previous studies indicated the utilization of glucose and N-acetyl sion in the long-term.
glucosamine by Chloroflexi was low at low pH (Miura et al., 2007). Al- In conclusion, soil pH would directly alter the bacterial diversity as
though the results of ferralsol confirmed the importance of pH, the well as indirectly changed microbiome through affecting other physico-
slightly decrease of bacteria diversity in phaeozems after biochar addi- chemical variables such as bioavailability of SOC. The alternation of SOC
tion appeared to contradict with this, implying there may be other bioavailability incorporated by the compositional heterogeneity of bio-
mechanisms. char highlighted the importance of the relationship between soil pH
Soil pH would indirectly alter microbial communities by affecting and nature of SOC on carbon sequestration. Previous study also indicat-
the molecular-level composition of native soil organic matter. Com- ed high C/N ratio led to low microbial N availability and suppressed soil
pared to alkaline soil, acidic ferralsol exhibited lower proportion of ali- C mineralization (Liu et al., 2016). The observed difference in bacterial
phatic hydrocarbon and higher aromatic compounds, resulted in diversity by soil pH would in turn affect the stability of biochar and
different microbiome. Although biochar decreased water extracted or- soil as well as the CO2 emission. As the response patterns of bacterial di-
ganic carbon (WEOC) through adsorption (Fidel et al., 2017) in both versity to physicochemical factors was synergistic effects, figuring out
soils, the proportion of aliphatic hydrocarbon and carbohydrates in- the combined effect of variables would be crucial in large-scale biochar
creased in ferralsol after biochar amendment, leading to higher bioavail- application.
ability of the substrate. Previous researches also showed higher
substrate limitation of bacteria in the acidic arenosol (Sheng et al., 4.3. Combined effect of physicochemical properties on microbial communi-
2016). As a result, the phylum Proteobacteria which exhibited ties change and carbon sequestration
copiotrophic attributes and were found abundant in soils with high C
availability (Fierer et al., 2007; Verastegui et al., 2014) increased in Soil pH was dominated factor in evaluating the richness of bacteria,
ferralsol but significantly decreased in phaeozems. In contrast, CO2 emission, as well as the priming effect. Although significant rela-
Acidobacteria which was enriched in soils with low resource availability tionship had been found for other variables, it could be ascribed to the
and associated with lower SOC mineralization rate (Fierer et al., 2007) co-variation of pH since the significant correlation (p b 0.01) among
decreased in ferralsol after biochar amendment. This was in consistence each variable (Table S3). This was in consistent with previous findings
with previous study which indicated lower substrate limitation and which found SOC and C/N ratio was significant in shaping the soil bacte-
higher gram-negative bacteria proportion in biochar-amend acidic rial community composition due to their high correlations with soil pH
soils (Sheng et al., 2016). In contrast, decrease of several aliphatic hy- (Xu et al., 2014; Zheng et al., 2016). Previous studies also suggested that
drocarbons as well as the WEOC after biochar addition contributed to microbial biogeography is controlled primarily by soil pH (Fierer and
higher diversity of oligotrophic bacteria such as Actinobacteria and Jackson, 2006; Rousk et al., 2010). Higher pH by biochar amendment re-
−
Chloroflexi in phaeozems. sulted in lower NH+ 4 content and higher adsorption of NO3 in both soils,
Besides, recalcitrant aromatic compounds such as quinoline, benzoic as well as higher activity of nitrobacteria, resulted in a relative higher
acid and phthalic acid increased in both soils after biochar addition, concentration of NO− 3 in ferralsol when the addition rate increased. Fur-
resulting higher abundance of bacteria to degrade refractory and aro- thermore, the decreased correlation coefficient between pH and PE at
matic compound (Table S5). The higher presence level of Actinobacteria the end of the incubation indicated the relative short-term effect of
in both biochar-amended ferralsol and phaeozems was due to its ability soil pH acted as a prior factor stimulating CO2 of SOC. This was in
1398 Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399
accordance with the initial drastic CO2 release after biochar addition, es- Cheng, H.G., Hill, P.W., Bastami, M.S., Jones, D.L., 2017. Biochar stimulates the decomposi-
tion of simple organic matter and suppresses the decomposition of complex organic
pecially in ferralsol. The long-term carbon sequestration of biochar was matter in a sandy loam soil. Glob. Change Biol. Bioenergy. 9, 1110–1121.
combined effect that may be dominated by other factors such as TC and Ding, G.C., Pronk, G.J., Babin, D., Heuer, H., Heister, K., Kogel-Knabner, I., et al., 2013. Min-
clay content. eral composition and charcoal determine the bacterial community structure in artifi-
cial soils. FEMS Microbiol. Ecol. 86, 15–25.
The strong correlation between CO2 emission, priming effect (PE) Farrell, M., Kuhn, T.K., Macdonald, L.M., Maddern, T.M., Murphy, D.V., Hall, P.A., et al.,
and soil diversity proposed the microbial driven carbon sequestration 2013. Microbial utilisation of biochar-derived carbon. Sci. Total Environ. 465,
of biochar, indicating soils with high pH and TC tend to have negative 288–297.
Feng, Z.J., Zhu, L.Z., 2017. Impact of biochar on soil N2O emissions under different biochar-
PE. The higher pH and bioavailability of organic C in the biochar carbon/fertilizer-nitrogen ratios at a constant moisture condition on a silt loam soil.
amended ferralsol resulted in higher abundance of copiotrophic bacte- Sci. Total Environ. 584, 776–782.
ria such as Chloroflexi, Gemmatimonadetes and Bacteroidetes. Since Feng, X.J., Simpson, A.J., Schlesinger, W.H., Simpson, M.J., 2010. Altered microbial commu-
nity structure and organic matter composition under elevated CO2 and N fertilization
copiotrophic bacteria have higher growth rates, greater degree of vari-
in the duke forest. Glob. Chang. Biol. 16, 2104–2116.
ability in population size, metabolic quotient and lower substrate affin- Fidel, R.B., Laird, D.A., Parkin, T.B., 2017. Impact of biochar organic and inorganic carbon
ities than oligotrophic bacteria (Fierer et al., 2007), higher proportion of on soil CO2 and N2O emissions. J. Environ. Qual. 46, 505–513.
copiotrophic bacteria then enhanced the CO2 emission and co-metabo- Fierer, N., Jackson, R.B., 2006. The diversity and biogeography of soil bacterial communi-
ties. Proc. Natl. Acad. Sci. U. S. A. 103, 626–631.
lism of SOC. On the contrary, as the bioavailability of the SOC determines Fierer, N., Bradford, M.A., Jackson, R.B., 2007. Toward an ecological classification of soil
whether biochar is a more attractive substrate and alternative C source bacteria. Ecology 88, 1354–1364.
(Cheng et al., 2017; Whitman et al., 2014), the abundance of oligotro- Gu, B.H., Schmitt, J., Chen, Z.H., Liang, L.Y., Mccarthy, J.F., 1994. Adsorption and desorption
of natural organic-matter on iron-oxide - mechanisms and models. Environ. Sci.
phic bacteria such as Actinobacteria and Anaerolineae increased in bio- Technol. 28, 38–46.
char-amended phaeozems, combined with decreased copiotrophic Hanada, A., Kurogi, T., Giang, N.M., Yamada, T., Kamimoto, Y., Kiso, Y., et al., 2014. Bacteria
bacteria such as Proteobacteria, mainly due to the strong adsorption of the candidate phylum TM7 are prevalent in acidophilic nitrifying sequencing-
batch reactors. Microbes Environ. 29, 353–362.
and decreased bioavailability of SOC by biochar. The microbial commu- Hu, L., Cao, L.X., Zhang, R.D., 2014. Bacterial and fungal taxon changes in soil microbial
nity change could also be implied by Kyoto Encyclopedia of Genes and community composition induced by short-term biochar amendment in red oxidized
Genomes (KEGG) pathways mapping (Figs. S2, S3). In ferralsol the me- loam soil. World J. Microbiol. Biotechnol. 30, 1085–1092.
Jenkins, J.R., Viger, M., Arnold, E.C., Harris, Z.M., Ventura, M., Miglietta, F., et al., 2017. Bio-
tabolism of saccharides and amino acid such as glycine, serine and thre- char alters the soil microbiome and soil function: results of next-generation amplicon
onine significantly increased, as well as the transporters, while in sequencing across Europe. Glob. Change Biol. Bioenergy. 9, 591–612.
phaeozems the secretion system and bacterial mobility protein signifi- Joseph, S.D., Camps-Arbestain, M., Lin, Y., Munroe, P., Chia, C.H., Hook, J., et al., 2010. An
investigation into the reactions of biochar in soil. Aust. J. Soil Res. 48, 501–515.
cantly decreased after biochar addition.
Kan, T., Strezov, V., Evans, T.J., 2016. Lignocellulosic biomass pyrolysis: a review of product
properties and effects of pyrolysis parameters. Renew. Sust. Energ. Rev. 57,
5. Conclusions 1126–1140.
Keith, A., Singh, B., Singh, B.P., 2011. Interactive priming of biochar and labile organic mat-
ter mineralization in a smectite-rich soil. Environ. Sci. Technol. 45, 9611–9618.
Biochar can affect soil microbial community by altering the soil Khodadad, C.L.M., Zimmerman, A.R., Green, S.J., Uthandi, S., Foster, J.S., 2011. Taxa-specific
physicochemical properties, thus to change the CO2 emission and prim- changes in soil microbial community composition induced by pyrogenic carbon
ing effect (PE) of SOC. CO2 emission was accelerated in acidic ferralsol amendments. Soil Biol. Biochem. 43, 385–392.
Kirby, R., 2005. Actinomycetes and lignin degradation. Adv. Appl. Microbiol. 58C,
but suppressed in phaeozems after biochar amendment, the accelera- 125–168.
tion and suppression were increased at higher biochar addition rate. Al- Kolton, M., Harel, Y.M., Pasternak, Z., Graber, E.R., Elad, Y., Cytryn, E., 2011. Impact of bio-
though biochar addition initially exhibited positive PE of SOC, the long- char application to soil on the root-associated bacterial community structure of fully
developed greenhouse pepper plants. Appl. Environ. Microbiol. 77, 4924–4930.
term PE of two soils was negative, largely due to the organo-mineral in- Kolton, M., Graber, E.R., Tsehansky, L., Elad, Y., Cytryn, E., 2017. Biochar-stimulated plant
teraction. Among all the parameters, pH played a prior role in increasing performance is strongly linked to microbial diversity and metabolic potential in the
the bacteria diversity, especially in the abundance of copiotrophic bac- rhizosphere. New Phytol. 213, 1393–1404.
Kuzyakov, Y., Subbotina, I., Chen, H.Q., Bogomolova, I., Xu, X.L., 2009. Black carbon decom-
teria such as Chloroflexi, Gemmatimonadetes and Bacteroidetes after bio- position and incorporation into soil microbial biomass estimated by C-14 labeling.
char addition in ferralsol, resulting in more CO2 emission. Whereas Soil Biol. Biochem. 41, 210–219.
decreased content and bioavailability of SOC via adsorption of biochar Lehmann, J., Joseph, S., 2009. Biochar for Environmental Management: Science and Tech-
nology. Earthscan, London.
resulted in higher abundance of oligotrophic bacteria in phaeozems, Lehmann, J., Liang, B.Q., Solomon, D., Lerotic, M., Luizao, F., Kinyangi, J., et al., 2005. Near-
thus decreasing the CO2 emission. Further research is required to ex- edge X-ray absorption fine structure (NEXAFS) spectroscopy for mapping nano-scale
plore the relationship between soil physicochemical properties and mi- distribution of organic carbon forms in soil: application to black carbon particles.
Glob. Biogeochem. Cycles 19.
crobial function in the long-term carbon sequestration of biochar.
Lehmann, J., Rillig, M.C., Thies, J., Masiello, C.A., Hockaday, W.C., Crowley, D., 2011. Biochar
effects on soil biota - a review. Soil Biol. Biochem. 43, 1812–1836.
Acknowledgement Lin, Y., Munroe, P., Joseph, S., Henderson, R., Ziolkowski, A., 2012a. Water extractable or-
ganic carbon in untreated and chemical treated biochars. Chemosphere 87, 151–157.
Lin, Y., Munroe, P., Joseph, S., Kimber, S., Van Zwieten, L., 2012b. Nanoscale organo-min-
This work was supported by the National Natural Science Founda- eral reactions of biochars in ferrosol: an investigation using microscopy. Plant Soil
tion of China (21477111 and 21621005), and the National Key Technol- 357, 369–380.
Liu, S.W., Zhang, Y.J., Zong, Y.J., Hu, Z.Q., Wu, S., Zhou, J., et al., 2016. Response of
ogy Research and Development Program of the Ministry of Science and
soil carbon dioxide fluxes, soil organic carbon and microbial biomass carbon
Technology of China (2015BAC02B00). to biochar amendment: a meta-analysis. Glob. Change Biol. Bioenergy. 8,
392–406.
Appendix A. Supplementary data Luo, Y., Durenkamp, M., De Nobili, M., Lin, Q., Brookes, P.C., 2011. Short term soil priming
effects and the mineralisation of biochar following its incorporation to soils of differ-
ent pH. Soil Biol. Biochem. 43, 2304–2314.
Supplementary data to this article can be found online at https://doi. von Lutzow, M., Kogel-Knabner, I., Ekschmitt, K., Matzner, E., Guggenberger, G.,
org/10.1016/j.scitotenv.2017.11.337. Marschner, B., et al., 2006. Stabilization of organic matter in temperate soils: mecha-
nisms and their relevance under different soil conditions - a review. Eur. J. Soil Sci. 57,
426–445.
References Maestrini, B., Nannipieri, P., Abiven, S., 2015. A meta-analysis on pyrogenic organic matter
induced priming effect. Glob. Change Biol. Bioenergy. 7, 577–590.
Cao, X.D., Ma, L.N., Gao, B., Harris, W., 2009. Dairy-manure derived biochar effectively Mitchell, P.J., Simpson, A.J., Soong, R., Simpson, M.J., 2016. Biochar amendment altered the
sorbs lead and atrazine. Environ. Sci. Technol. 43, 3285–3291. molecular-level composition of native soil organic matter in a temperate forest soil.
Chen, J.H., Liu, X.Y., Zheng, J.W., Zhang, B., HF, Lu, Chi, Z.Z., et al., 2013. Biochar soil amend- Environ. Chem. 13, 854–866.
ment increased bacterial but decreased fungal gene abundance with shifts in com- Miura, Y., Watanabe, Y., Okabe, S., 2007. Significance of Chloroflexi in performance of sub-
munity structure in a slightly acid rice paddy from Southwest China. Appl. Soil Ecol. merged membrane bioreactors (MBR) treating municipal wastewater. Environ. Sci.
71, 33–44. Technol. 41, 7787–7794.
Y. Sheng, L. Zhu / Science of the Total Environment 622–623 (2018) 1391–1399 1399
Murray, J., Keith, A., Singh, B., 2015. The stability of low- and high-ash biochars in acidic transform ion cyclotron resonance mass spectrometry. Environ. Sci. Technol. 47,
soils of contrasting mineralogy. Soil Biol. Biochem. 89, 217–225. 13294–13302.
Nguyen, B.T., Lehmann, J., Kinyangi, J., Smernik, R., Riha, S.J., Engelhard, M.H., 2009. Long- Spokas, K.A., 2013. Impact of biochar field aging on laboratory greenhouse gas production
term black carbon dynamics in cultivated soil. Biogeochemistry 92, 163–176. potentials. Glob. Change Biol. Bioenergy. 5, 165–176.
O'Neill, B., Grossman, J., Tsai, M.T., Gomes, J.E., Lehmann, J., Peterson, J., et al., 2009. Bacte- Sun, D.Q., Meng, J., Liang, H., Yang, E., Huang, Y.W., Chen, W.F., et al., 2015. Effect of volatile
rial community composition in Brazilian anthrosols and adjacent soils characterized organic compounds absorbed to fresh biochar on survival of Bacillus mucilaginosus
using culturing and molecular identification. Microb. Ecol. 58, 23–35. and structure of soil microbial communities. J. Soils Sediments 15, 271–281.
Pietri, J.C.A., Brookes, P.C., 2009. Substrate inputs and pH as factors controlling microbial Verastegui, Y., Cheng, J., Engel, K., Kolczynski, D., Mortimer, S., Lavigne, J., et al., 2014.
biomass, activity and community structure in an arable soil. Soil Biol. Biochem. 41, Multisubstrate isotope labeling and metagenomic analysis of active soil bacterial
1396–1405. communities. MBio 5.
Qian, L.B., Chen, B.L., 2014. Interactions of aluminum with biochars and oxidized biochars: White, D.C., Sutton, S.D., Ringelberg, D.B., 1996. The genus Sphingomonas: physiology and
implications for the biochar aging process. J. Agric. Food Chem. 62, 373–380. ecology. Curr. Opin. Biotechnol. 7, 301–306.
Rabus, R., 2005. Functional genomics of an anaerobic aromatic-degrading denitrifying Whitman, T., Zhu, Z.H., Lehmann, J., 2014. Carbon mineralizability determines interactive
bacterium, strain EbN1. Appl. Microbiol. Biotechnol. 68, 580–587. effects on mineralization of pyrogenic organic matter and soil organic carbon. Envi-
Rousk, J., Brookes, P.C., Baath, E., 2009. Contrasting soil pH effects on fungal and bacterial ron. Sci. Technol. 48, 13727–13734.
growth suggest functional redundancy in carbon mineralization. Appl. Environ. Xu, H.J., Wang, X.H., Li, H., Yao, H.Y., Su, J.Q., Zhu, Y.G., 2014. Biochar impacts soil microbial
Microbiol. 75, 1589–1596. community composition and nitrogen cycling in an acidic soil planted with rape. En-
Rousk, J., Brookes, P.C., Baath, E., 2010. The microbial PLFA composition as affected by pH viron. Sci. Technol. 48, 9391–9399.
in an arable soil. Soil Biol. Biochem. 42, 516–520. Zheng, J.F., Chen, J.H., Pan, G.X., Liu, X.Y., Zhang, X.H., Li, L.Q., et al., 2016. Biochar decreased
Sheng, Y.Q., Zhan, Y., Zhu, L.Z., 2016. Reduced carbon sequestration potential of biochar in microbial metabolic quotient and shifted community composition four years after a
acidic soil. Sci. Total Environ. 572, 129–137. single incorporation in a slightly acid rice paddy from southwest China. Sci. Total En-
Singh, B.P., Cowie, A.L., 2014. Long-term influence of biochar on native organic carbon viron. 571, 206–217.
mineralisation in a low-carbon clayey soil. Sci. Rep. 4. Zimmerman, A.R., Gao, B., Ahn, M.Y., 2011. Positive and negative carbon mineralization
Singh, B.P., Cowie, A.L., Smernik, R.J., 2012. Biochar carbon stability in a clayey soil as a priming effects among a variety of biochar-amended soils. Soil Biol. Biochem. 43,
function of feedstock and pyrolysis temperature. Environ. Sci. Technol. 46, 1169–1179.
11770–11778.
Smith, C.R., Sleighter, R.L., Hatcher, P.G., Lee, J.W., 2013. Molecular characterization of
inhibiting biochar water-extractable substances using electrospray ionization Fourier