Evidences of N2O Emissions in Chloropicrin-Fumigated Soil - Journal

Oct 19, 2018 - ... of Agricultural Sciences , Beijing 100193 , People's Republic of China .... Continuous Flow Analytical System, Alliance Instruments...
1 downloads 0 Views 1MB Size
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