Subscriber access provided by UNIV OF LOUISIANA
Agricultural and Environmental Chemistry
Evidences of N2O emissions in chloropicrin-fumigated soil wensheng fang, Dongdong Yan, xianli wang, bin huang, Zhaoxin Song, jie liu, Xiaoman Liu, Qiuxia Wang, Yuan Li, canbin ouyang, and Aocheng Cao J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b04351 • Publication Date (Web): 19 Oct 2018 Downloaded from http://pubs.acs.org on October 20, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 40
Journal of Agricultural and Food Chemistry
1 2 3 4
Evidences of N2O emissions in chloropicrin-fumigated soil
5
Wensheng Fang, Dongdong Yan, Xianli Wang, Bin Huang, Zhaoxin Song, Jie Liu,
6
Xiaoman Liu, Qiuxia Wang, Yuan Li, Canbin Ouyang, Aocheng Cao*
7
Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193,
8
China
9
*Corresponding
author email address:
[email protected], Tel: +86-62-815-904
10 11 12 13 14 15 16 17 18 19 20 21 1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
22 23
Abstract
24
The mechanism of N2O production following chloropicrin (CP) fumigation was
25
investigated in this study. Our results showed that CP fumigation increased N2O production
26
from 23 to 25 times compared with the control, and significantly decreased the abundance
27
of 16S rRNA and N-cycling functional genes. CP also decreased the soil bacterial diversity
28
and caused a shift in the community composition. The N2O emissions in fumigated soil
29
were significantly correlated with soil environmental factors (NH4+, dissolved amino acid,
30
microbial biomass nitrogen, and NO3−) but were not correlated with the abundance of
31
functional genes. Metatranscriptomes and dual-label 15N-18O isotope analysis revealed that
32
CP fumigation inhibited the expression of gene families involved in N2O production and
33
sink processes and shifted the main pathway of N2O production from nitrification to
34
denitrification. These results provided useful information for environmental safety
35
assessments of CP in China, for improving our understanding of the N-cycling pathways
36
in fumigated soils.
37
Keywords: Chloropicrin, N2O emissions, Functional genes, Nitrogen cycling,
38
Microorganisms
39
Introduction
40
Soil fumigation has been used to control soilborne pathogens, nematodes, and weeds
41
for several decades worldwide. Methyl bromide (MeBr) was the most used and best
42
fumigant for more than 50 years 1. However, the use of MeBr has been completely 2
ACS Paragon Plus Environment
Page 2 of 40
Page 3 of 40
Journal of Agricultural and Food Chemistry
43
prohibited in agriculture system since 2015 because of its detrimental effect on ozone 2. As
44
one of the most suitable alternatives to MeBr, chloropicrin (CP) does not deplete ozone
45
and exhibits high efficacy in killing board-spectrum soil-borne diseases by rapidly
46
immersing into the tissues of biota and killing cells 3. To date, CP has been widely
47
employed in crops, such as strawberry, cucumber, tomato, ginger, melon, and eggplant, to
48
control soil-borne pathogens and plant-parasitic nematodes 4-6.
49
Soil fumigation can effectively target organisms but also affects off-target
50
microorganisms, which are important to the soil microbial community. Studies reported
51
that CP causes a significant decrease in bacterial diversity and shift on the population size
52
and soil bacterial composition
53
transforming activities, including mineralization stimulation 9 and nitrification reduction
54
10,11.
55
of microbial community in increasing N2O emissions and the biochemical pathways used
56
by microorganisms to produce N2O emissions after CP fumigation have not been reported
57
yet.
7,8.
CP fumigation triggers significant effects on N-
Soil treated with CP increased N2O emissions by 12.6 times 12,13. However, the role
58
N2O is a potent greenhouse gas with a global warming potential that is 298-fold higher
59
than that of CO2 14. N2O is one of the most significant ozone depletes that affect the long-
60
term recovery of the ozone layer
61
primary sources of N2O. Nitrification starts with the oxidation of NH3 to NO2− and of NO2−
62
to NO3− 16. Under aerobic conditions, chemoautotrophs use NH3 as energy source 17,18 to
63
produce N2O through several pathways, such as nitrifier nitrification (NN), nitrifier
15.
Microbial nitrification and denitrification are the
3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 40
64
denitrification (ND), and nitrification-coupled denitrification (NCD). The biological
65
denitrification
66
NO3−→NO2−→NO→N2O→N2; these steps are catalyzed by NO3− reductase (nar/nap),
67
NO2− reductases (Nir), nitric oxide reductases (Nor), and N2O reductase (nos), respectively
68
19.
69
induced by N2O-reducing microorganisms indicates that they are the only microbial sink
70
for N2O
71
codenitrification
72
lead to N2O formation.
process
consists
of
four
reaction
steps,
namely,
N2O can result from a multiplicity of pathways, but the process of N2O reduction to N2
20.
Other microbial processes, such as heterotrophic nitrification 22,
21,
chemodenitrification, and dissimilatory NO3− reduction to NH3, can
73
In this study, we examined the production of N2O and the response of different groups
74
of N-cycling microbes living in Jiangxi lateritic red and Beijing Fluvo-aquic soils
75
fumigated with CP. At 59 days after fumigation, we quantified N2O emission rates, changes
76
in geochemical parameters NO3− and NH4+, microbial biomass nitrogen (MBN), and
77
dissolved amino acids (DAA). We also used real-time PCR to monitor the abundance of
78
key functional marker genes (nifH, AOA amoA, AOB amoA, nxrB, narG, napA, nirK, nirS,
79
cnorB, qnorB, and nosZ) involved in microbial N fixing, nitrification, and denitrification.
80
Moreover, 16S rRNA gene amplicon sequencing techniques were used to determine the
81
diversity and community structure of N-transfer microorganisms in the fumigated soils.
82
Finally, we analyzed dual-labeled 15N-18O isotopes and metatranscriptome to observe shifts
83
in the N2O pathway and determine the possible production mechanism of N2O in the
84
fumigated soils. We aimed to (i) discuss the biochemical pathways used by microorganisms 4
ACS Paragon Plus Environment
Page 5 of 40
Journal of Agricultural and Food Chemistry
85
to produce N2O emissions from soils fumigated with CP and (ii) determine the response
86
and recover rate of the N-cycling microorganisms in CP-fumigated soil at the post-
87
fumigation phase. Results are predicted to improve the understanding of N-cycling
88
mechanisms in fumigated soil ecosystems and provide users with valuable guidance
89
regarding the efficient and effective management of CP.
90
Methods and Materials
91
Soil sample collection
92
Soil samples were obtained from the top 20 cm of two agricultural fields in Beijing,
93
Shunyi and Jiangxi, Yudu. The two samples are typical alkaline Fluvo-aquic (Beijing) and
94
acidic lateritic red soils (Jiangxi). Detailed soil analytical data are presented in Table 1.
95
The soils were preincubated for 7 days at room temperature (25 ± 5 °C) in the dark and
96
then adjusted to 40% water-holding capacity (WHC) before fumigation.
97
Experiment I: Post-fumigation soil incubation experiment
98
Experimental setup
99
Microcosms were prepared with 220 g of sieved soil in 500 mL Duran® wide-neck
100
glass bottles (Schott AG, Mainz, Germany). CP (99.5% purity) was obtained from Dalian
101
Dyestuffs & Chemicals Co., Dalian, Liaoning, China. CP was added at a typical field
102
application rates of 65 mg kg−1 to both soil types as treatment groups. The control group
103
was treated with deionized water. Each treatment was prepared triplicates. The soil
104
microcosm bottles were closed with butyl rubber stoppers with outlet ports for syringe
105
sampling and incubated at the ambient temperature of 28 °C in daylight. 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
106
Soil sampling and geochemical analyses
107
N2O emission rates were measured using the method previously described 14,23. During
108
the incubation stage, the soil samples in the bottles were thoroughly stirred for 10-15 min
109
every day to create an aerobic incubation environment. At the designated time, a 10 mL
110
gas sample was withdrawn from each bottle using a gas-tight syringe (Agilent
111
Technologies, USA). The gas samples were transferred to a 21 mL headspace vial that was
112
flushed with helium and evacuated (10 mL) before use. An Agilent 7890A gas
113
chromatograph with an electron capture detector (63Ni-ECD) coupled with an Agilent
114
7694E headspace sampler were used to quantify the N2O concentrations. The gas
115
chromatograph conditions and gas emission calculation method were described previously
116
24.
117
any residual fumigant. After sampling, all bottles were pumped with fresh air and then
118
returned into the incubator. During incubation, water content was controlled
119
gravimetrically and adjusted each week to the initial WHC by adding deionized water by
120
a spray bottle. Gas samples were collected every week, and soil samples were collected at
121
days 10, 24, 38, and 59 for geochemical and molecular biological analyses.
After each gas sample, the microcosm bottles were opened in a ventilation hood to vent
122
Soil samples were extracted for mineral N (NH4+-N and NO3–-N) by using 2 M KCl.
123
Mineral N concentrations were determined with a continuous flow automated analyzer
124
(Futura Continuous Flow Analytical System, Alliance Instruments, France). MBN was
125
estimated by the chloroform fumigation method
126
ninhydrin reaction method 26.
25,
and DAA was measured using the
6
ACS Paragon Plus Environment
Page 6 of 40
Page 7 of 40
127
Journal of Agricultural and Food Chemistry
Real-time quantitative PCR and high-throughput sequencing
128
Soil total genomic DNA was extracted from 0.25 g of each soil sample by using a
129
MoBio Powersoil® DNA Isolation Kit (MoBio Laboratories, CA, USA) according to the
130
manufacturer’s protocol. Quantification of functional marker genes (16S rRNA gene
131
[bacteria], amoA [bacteria and archaea], nifH, napA, narG, nirK, nirS, cnorB, qnorB, and
132
nosZ) was carried out using the SsoFast EvaGreen Supermix (Bio-Rad Laboratories, CA,
133
USA) and gene-specific primers. Quantitative PCR was performed using CFX96 Real-time
134
PCR System (Bio-Rad, CA, USA). The details of gene-specific qPCR primers, reaction
135
mixtures, and thermal programs are listed in Tables S1, S2, and S3, respectively. To
136
analyze the abundance and community structure of N-transfer microbes, Majorbio Bio-
137
pharm Technology Co., Ltd. (Shanghai, China) conducted MiSeq sequencing of the 16S
138
rRNA genes in the V3–V4 regions by using the total DNA extracted from Beijing soil
139
microorganisms. The raw reads were deposited into the NCBI Sequence Read Archive
140
database (no. SRP124701). All the certified nitrogen cycle related functional microbes
141
were collected from the NCBI database (https://www.ncbi.nlm.nih.gov/) (showed in Table
142
S4 & Table S5).
143
Experiment II: Short-term fumigation experiment
144
Experimental setup
145
A total of 3 g (dry weight) incubated Beijing soil samples (the same soils used in
146
Experiment I) were placed into a 21 mL clean headspace vial. The following treatments
147
for each fumigant were established: (1) 18O-H2O + NO3− + NH4+; (2) H2O + 18O-NO3− + 7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 40
148
NH4+; (3) H2O + 15N-NO3− + NH4+; and (4) H2O + NO3− + 15N-NH4+. The isotopes 18O and
149
15N
150
treatments received the total mineral N contents in the soil of 50 mg NH4+-N and 50 mg
151
NO3−-N kg−1 soil. All treatments were incubated at 40% of WHC at 28 °C. Each treatment
152
was prepared in triplicates for each isotope and soil total RNA extraction. The gas and soil
153
samples were taken at the end of 70 h of incubation. Our preliminary experiment found
154
that the volume of N2O emitted from the soil peaked at 70 h after fumigation.
155
Isotope-labeled gas and soil sample analyses
(H2O, NO3−, NH4+) were enriched with 1.0%
18O
and 10.0%
15N,
respectively. All
156
The rate of N2O emissions and its isotopic signatures were recorded at the Institute of
157
Environment and Sustainable Development in Agriculture (CAAS) by using an IsoPrime
158
trace gas concentration system interfaced to an IsoPrime100 isotope ratio mass
159
spectrometer (IsoPrime Ltd., UK). Soil samples, which were extracted with 1 M KCl
160
followed by segmented flow analyses, were obtained after gas sampling (Skalar Analytical
161
B.V., Breda, The Netherlands). The concentrations of the 15N enrichments were analyzed
162
using an elemental analyzer interfaced to a continuous flow isotope ratio mass spectrometer
163
(Sercon 20-20, Sercon Ltd., UK). The 18O and 15N isotopic enrichments of the accumulated
164
N2O, the 18O signature in H2O and NO3−, and the 15N signatures in NO3− and NH4+ were
165
used to calculate the N2O produced by NN, ND, NCD, and heterotrophic denitrification
166
(HD) N2O pathways, respectively 27,28.
167
Metatranscriptomic analysis
168
Total RNAs were extracted from nine samples by using the EZNA® Soil RNA Midi 8
ACS Paragon Plus Environment
Page 9 of 40
Journal of Agricultural and Food Chemistry
169
Kit (Omega Bio-tek Inc., USA) according to manufacturer’s protocols. The RNA
170
concentration
171
spectrophotometer (Thermo Fisher Scientific, USA). RNA quality was examined using a
172
1% agarose gel electrophoresis system and assessed using RNA 6000 Nano Kit and
173
Reagents (total RNA) in an Agilent 2100 Bioanalyzer (Agilent Technologies, USA). The
174
total RNA contained in samples was based on an rRNA removal procedure by using the
175
Ribo-Zero rRNA Removal Kit (Illumina Inc., USA) according to the manufacturer’s
176
instructions. cDNA libraries were constructed using the TruSeq™ RNA Library Prep Kit
177
v2 (Illumina Inc., USA) according to the manufacturer’s instructions. The barcoded
178
libraries were paired end sequenced on the Illumina HiSeq 2500 System (Illumina Inc.,
179
USA) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China) by using HiSeq
180
3000/4000 PE Cluster Kits (Illumina Inc., USA) according to the manufacturer’s
181
instructions. A total of 325,931,416 raw reads were obtained. Quality control was carried
182
out by the SeqPrep (https://github.com/jstjohn/SeqPrep) to strip the 3′ and 5′ ends, as well
183
as by Sickle (https://github.com/najoshi/sickle) and SortMeRNA
184
reads (length < 50 bp, with a quality value < 20, or possessing N bases) and rRNA reads,
185
thereby resulting in 200,925,212 high-quality sequences. A total of 30,093 open reading
186
frames with an average length 463 bp from each sample were predicted using
187
TransGeneScan (http://sourceforge.net/projects/transgenescan/). All sequences with a 95%
188
sequence identity (90% coverage) were clustered as nonredundant gene catalog by the
189
program CD-HIT (http://www.bioinformatics.org/cd-hit/). Reads after quality control were
and
purity
were
quantified
using
9
ACS Paragon Plus Environment
a
NanoDrop™
2000/2000c
to remove low-quality
Journal of Agricultural and Food Chemistry
190
mapped to the representative genes with 95% identity and reads per kilobase of transcript
191
per million mapped reads and evaluated using RNA-Seq by expectation maximization
192
(http://deweylab.biostat.wisc.edu/rsem/). Functional and taxonomic assignments of the
193
predicted genes used the Basic Local Alignment Search Tool (BLAST) program to
194
compare nucleotide or protein sequences contained in the databases eggNOG, KEGG, and
195
NCBI nr. The BLAST E-value cut-off was set at 1.e−5. All metatranscriptomic data were
196
deposited in the NCBI database (no. SRP143409).
197
Statistical analysis
198
The major effects of fumigation on the biochemical parameters and abundance of
199
functional genes involved in N cycling was studied using univariate ANOVA with least
200
significant difference test by using the SPSS statistics software package version 18.0 (IBM,
201
USA). Univariate ANOVA was used to reveal differences between the control and
202
fumigated microcosms. All concentrations or gene copy number values from the control at
203
each time point of sampling were individually compared with the fumigated soil
204
microcosms. The α-diversity indices Chao1, ACE, Shannon, and Simpson were calculated
205
using the software Mothur to determine the diversity of the bacterial communities in
206
Beijing soil. Hierarchical clustering and a heat map were used to determine changes in the
207
relative abundance of the bacterial genera involved in N cycling. The Bray–Curtis
208
algorithm was used to show the relative abundance of these genera according to
209
hierarchical clustering. A heat map figure and statistical correlations were generated to
210
show the sequencing results by using pheatmap-package and vegan package in R, 10
ACS Paragon Plus Environment
Page 10 of 40
Page 11 of 40
Journal of Agricultural and Food Chemistry
211
respectively (version 2.15.3). Linear discriminant analysis coupled with effect size
212
measurement was used to identify the relative abundance of the bacterial genera in CP-
213
fumigated and unfumigated Beijing soils. Spearman’s rank correlation coefficient was used
214
to determine the relationship among N2O emission rate, physicochemical parameters, and
215
microbial functional genes.
216
Results
217
N2O production rates
218
The N2O production rates increased by 23.5- and 25.4-fold within 10 days of CP
219
fumigation in Beijing and Jiangxi soil compared with unfumigated soils, respectively
220
(Figure 1). However, N2O release rates gradually decreased at the post-fumigation phase
221
and returned to emission levels similar to that of the control at days 45 and 38 in Beijing
222
soil and Jiangxi soil.
223
Changes in physiochemical parameters
224
Compared with unfumigated soils, CP fumigation significantly increased the NH4+-N
225
and DAA contents in both soil types, whereas NO3−-N and MBN concentrations were
226
significantly decreased (Table 2). All these parameters recovered to levels similar to those
227
of the controls 59 days after fumigation in Beijing soil. However, the NH4+-N, NO3--N, and
228
MBN values were not recovered to the control levels until 59 days of fumigation in Jiangxi
229
soil. This result suggested that CP induced larger changes in these parameters that persisted
230
for longer in Jiangxi soil than in Beijing soil.
231
Changes in abundance of 16S rRNA and N-cycling functional marker genes 11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
232
CP fumigation induced a significant decrease in 16S rRNA gene abundance in Beijing
233
soil, and this inhibitory effect was evident in fumigated Jiangxi soil (Figure 2). The total
234
bacterial abundance of Beijing soil following CP fumigation recovered to the control level
235
38 days after treatment. However, in Jiangxi soil, CP fumigation inhibited the total
236
bacterial abundance for the entire incubation period. The copy number of the observed 11
237
N-cycling functional genes (nifH, AOA amoA, AOB amoA, nxrB, napA, narG, nirK, nirS,
238
cnorB, qnorB, and nosZ) were significantly decreased in CP-fumigated Beijing and Jiangxi
239
soils (Figure 2–4). In CP fumigated Jiangxi soil, only gene populations nxrB and nirS
240
recovered to the control level 59 days after treatment. In CP-fumigated Beijing soil,
241
functional genes, such as AOB amoA, nxrB, nirK, nirS, and cnorB recovered to levels
242
similar to that of the control at day 38.
243
Changes in soil bacterial diversity and community composition
244
The Shannon, ACE, and Chao1 diversity indices decreased significantly (p < 0.01)
245
relative to those in the control group in CP-fumigated Beijing soil samples during the entire
246
incubation phase (Table 3), but the Simpson diversity index significantly increased (p
0.69, p < 0.0002) and DAA (r > 0.83, p