Molecular chemodiversity of dissolved organic matter in paddy soils

111. Soil (12 g) was placed in polypropylene centrifuge tubes and Milli-Q water (60 mL;. 112. i.e. 1:5 w/v) added. Samples were then shaken (170 rpm, ...
2 downloads 10 Views 9MB Size
Subscriber access provided by READING UNIV

Article

Molecular chemodiversity of dissolved organic matter in paddy soils Xiaoming Li, Guo-Xin Sun, Songcan Chen, Zhi Fang, Hai-Yan Yuan, Quan Shi, and Yong-Guan Zhu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b00377 • Publication Date (Web): 04 Jan 2018 Downloaded from http://pubs.acs.org on January 4, 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 free 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 accessible to all readers and 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.

Environmental Science & Technology 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

Environmental Science & Technology

1

Molecular chemodiversity of dissolved organic matter in

2

paddy soils

3 4

Xiao-Ming Li1,4; Guo-Xin Sun1; Song-Can Chen1,4; Zhi Fang2; Hai-Yan Yuan1,4; Quan

5

Shi2; Yong-Guan Zhu1,3*

6

7

1

8

Eco-Environmental Sciences, Chinese Academy of Sciences, Shuangqing Road, No.

9

18, Haidian District, Beijing 100085, People’s Republic of China

10

2

State Key Laboratory of Urban and Regional Ecology, Research Center for

State Key Laboratory of Heavy Oil Processing, China University of Petroleum, 18

11

Fuxue Road, Changping, Beijing 102249, China

12

3

13

Chinese Academy of Sciences, Jimei Road, No. 1799, Jimei District, Xiamen 361021,

14

People’s Republic of China

15

4

16

District, Beijing 100049, People’s Republic of China

Key Laboratory of Urban Environment and Health, Institute of Urban Environment,

University of Chinese Academy of Sciences, Yuquan Road, No. 19A, Shijingshan

17 1

ACS Paragon Plus Environment

Environmental Science & Technology

18

*Corresponding author. Phone: 86-10-62936940, E-mail: [email protected].

2

ACS Paragon Plus Environment

Page 2 of 40

Page 3 of 40

Environmental Science & Technology

19

Abstract:

20

Organic matter (OM), and dissolved organic matter (DOM), have a major

21

influence upon biogeochemical processes; most significantly, the carbon cycle. To

22

date, very few studies have examined the spatial heterogeneity of DOM in paddy soils.

23

Thus, very little is known about the DOM molecular profiles and the key

24

environmental factors that underpin DOM molecular chemodiversity in paddy soils.

25

Here, Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS)

26

was applied to unambiguously resolve 11,361 molecular formulae in 16 paddy soils;

27

thereby elucidating the molecular characteristics of paddy soil DOM. Soil pH, iron

28

complexing index (Fep/FeR) and C/N ratio were established to be key factors

29

controlling DOM profiles. Polycyclic aromatics (derived from combustion) and

30

polyphenols (derived from plants) increased with increasing pH, while polyphenols

31

molecules, pyrogenic
aromatics and carboxylic compounds decreased with increasing

32

iron complexing index. Patterns in molecular profiles indicated DOM in paddy soils

33

to become more recalcitrant at higher soil C/N ratio and higher pH. Furthermore,

34

plant-derived polyphenols and pyrogenic DOM were retained favorably by iron and

3

ACS Paragon Plus Environment

Environmental Science & Technology

35

the chemodiversity of DOM in paddy soil increased with increasing soil C/N ratios.

36

This study provides critical information about DOM characteristics at a molecular

37

level and will inform better global management of soil carbon in paddy soil

38

ecosystems.

4

ACS Paragon Plus Environment

Page 4 of 40

Page 5 of 40

Environmental Science & Technology

39 40

TOC / Abstract Art

41 42

5

ACS Paragon Plus Environment

Environmental Science & Technology

43

Page 6 of 40

Introduction

44

Dissolved organic matter (DOM) is ubiquitous in terrestrial ecosystems and the

45

cycling of this carbon is of pivotal importance in global biogeochemistry and

46

greenhouse gas emissions1, 2. Thus, characterizing DOM will allow a better

47

understanding of the global carbon cycling3; and, other biogeochemical cycles linked

48

to it4. Given the importance of this topic to environmental and ecological fields5, 6 the

49

molecular characterization of DOM has become a priority. While efforts have been

50

made to evaluate DOM quality7, its extreme complexity and heterogeneity presents

51

challenges to understanding its molecular characteristics and reactivity8, 9.

52

Recent progress in advanced mass spectrometry and data processing has enabled

53

detailed characterization of DOM molecular profiles and has allowed the influence of

54

broad-scale environmental gradients on DOM chemodiversity to be defined10,

55

Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) is a

56

robust approach through which to obtain detailed information on individual molecules

57

present in complex mixtures, such as natural DOM12-14. FT-ICR MS facilitates the

58

assignment of formulae to thousands of DOM moieties; thus providing a detailed

59

profiling of DOM chemodiversity15. Thus, the information generated by FT-ICR MS

60

provides new molecular understanding regarding sources and processing of DOM.

11

.

61

While this approach has recently been successfully used to characterize DOM

62

from several environments (e.g. ocean water11, 16-18, lake water12, soil pore water19,

63

sediment20, 21 and the atmosphere20), the molecular diversity of soil DOM has rarely

64

been characterized9, 22; this is particularly true of paddy soils. Paddy soils are of global 6

ACS Paragon Plus Environment

Page 7 of 40

Environmental Science & Technology

65

significance as they have been identified as one of the major sources of greenhouse

66

gases3 (emitting between 60 and 100 Mt methane per year4). This information, if

67

sustained, could be used to regulate and forecast the geochemical processes involving

68

DOM in paddy soils. About 30.3 Mha of paddy soil cover ~27% of the cropped area

69

in China (FAOSTAT, 2014)23. Thus, China presents an ideal model to establish the

70

influence of environmental gradients on chemodiversity profiles of DOM in paddy

71

soils.

72

We hypothesize that, the molecular composition and diversity of paddy soil

73

DOM are related to entrained soil chemical attributes. To elucidate the influence of

74

environmental gradients on DOM chemodiversity sixteen representative paddy soils

75

were collected from across China (from northeastern to southeastern regions - latitude

76

20.55 - 47.58 and longitude 105.29 - 132.74). Across this gradient fundamental soil

77

attributes varied: soil pH (5.0-8.1), soil C/N ratio (8.5-15.4) and dissolved organic

78

carbon (DOC) concentration (105-485 mg kg-1) (Table S1). Trough interrogation of

79

these samples using FT-ICR MS, we elucidate: 1) DOM chemodiversity profiles in

80

paddy soils and how it varies across a broad environmental gradient; 2) the main

81

factors driving DOM molecular chemodiversity in paddy soils, and; 3) specific

82

chemo-marker molecules that can be used to ascribe fundamental differences in DOM

83

composition.

84

Methods

85

Soil sampling

86

Samples were collected in October 2014 from 16 paddy fields distributed across 7

ACS Paragon Plus Environment

Environmental Science & Technology

87

China (Fig. S1) at 5 cm depth. Samples were transported to the laboratory, on ice, in

88

sterile plastic bags. Soil chemical properties were analyzed according to standard

89

methods24. Soil pH was measured using a digital pH meter (PHS-3C, Shanghai Lida

90

Instrument Company, China) at a soil:deionized-water ratio of 1:2.5. Total carbon and

91

total nitrogen (TC, TN) were determined via dry combustion in an element analyzer

92

(Vario EL III - Elementar, Germany)25. The concentration of DOC was measured

93

using a TOC analyzer (Liquic TOC - Elementar, Germany). Nitrate and ammonium

94

were extracted from soil using KCl (2 mol·L-1) and measured using a continuous flow

95

analyzer (SAN++, Skalar, Holand)26.

96

Quantitative mineralogical analysis of the main mineral elements in soil samples

97

(such as, Fe, Al, Mn, Ca, Mg and Si) were carried out by selective dissolution

98

procedures and quantified by Inductively Coupled Plasma-Optical Emission

99

Spectrometer (ICP-OES; Optimal 2000 DV, Perkin-Elmer, USA)27. To determine the

100

total mineral elements (Fetot, Altot, Mntot, Catot, Mgtot and Sitot) in the soils, samples were

101

wet-digested using methods described previously28. Quantitation of reactive Fe

102

minerals (FeR) was performed using a citrate-bicarbonate-dithionite (CBD) solution

103

(0.1 mol·L-1)29. Short-range ordered (SRO) minerals were determined following

104

extraction (4 h) with ammonium oxalate solution (0.2 mol·L-1 at pH 3)22. The

105

organically complexed minerals (Fep, Alp, Mnp, Cap, Mgp and Sip) were prepared by

106

extraction with sodium-pyrophosphate (0.1 mol·L-1) for 16 h and analyzed by

107

ICP-OES24.

108

DOM sample preparation for FT-ICR MS using solid-phase extraction 8

ACS Paragon Plus Environment

Page 8 of 40

Page 9 of 40

Environmental Science & Technology

109

Water was used to extract DOM that was, thereafter, isolated with cartridges

110

obtained from Agilent Technologies (Bond Elut PPL, 500 mg, 6 mL)30. Milli-Q 18

111

mΩ water (DDI; Millipore, Bedford, MA, USA) was used throughout all experiments.

112

Soil (12 g) was placed in polypropylene centrifuge tubes and Milli-Q water (60 mL;

113

i.e. 1:5 w/v) added. Samples were then shaken (170 rpm, at 25 ± 1°C) for 8 h31.

114

Thereafter, the tubes were centrifuged (at 2,800 × g) for 10 min32. The supernatant

115

was filtered through a 0.45 µm mixed cellulose ester membrane.

116

HPLC methanol (20 mL) was passed through PPL cartridges to clean them; these

117

were then conditioned with acidified (pH 2) Milli-Q water (20 mL)33. To increase the

118

subsequent extraction efficiency of organic molecules loaded onto the PPL cartridges,

119

soil derived DOM water extracts were acidified to pH 2 (with HCl, p.a. grade, Merck,

120

Germany)34. Sufficient volume of sample (to deliver 100 µg C; informed by the

121

predetermined DOC content of the sample) was then loaded on to the PPL Cartridges

122

under gravity11. The cartridges were then rinsed with 20 mL of acidified Milli-Q water

123

(pH = 2)5 and ultrapure N2 gas blown onto the cartridges to completely dried them.

124

DOM was collected from the cartridges by eluting them using HPLC grade methanol

125

(5 mL; Merck, Germany). The elutes were stored in acid washed (and combusted)

126

glass vials and kept in the dark in the refrigerator (-20°C) prior to FT-ICR MS

127

measurement. The scheme for DOM sample preparation is provided in Figure S3.

128

FT-ICR MS analysis

129

Bruker Apex Ultra FT-ICR MS (equipped with a 9.4 T superconducting magnet5

130

interfaced with negative-ion mode electrospray ionization35) was used to perform 9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 40

131

ultra-high resolution mass spectrometry analysis of DOM samples. Diluted Suwannee

132

River Natural Organic Matter solution (50 mg L-1) (obtained from IHSS, USA) was

133

used for calibration. To compare the relative intensities of MS spectra in various

134

samples deuterated stearic acid (C18D35H1O2; Sigma-Aldrich, USA) was added as an

135

internal standard to the samples13. Ammonium hydroxide was added prior to

136

electrospray to increase the ionization efficiency34. Samples (180 µL h-1) were

137

injected into the electrospray source (voltages of the capillary and spray shield were,

138

respectively, 4 kV and 3.5 kV). Spectra were acquired over 128 scans (with an ion

139

accumulation time of 0.6 s; ion flight time into the ICR cell of 0.0012 s, and a m/z

140

range of 150 - 1,200). PPL extraction blanks and solvent blanks were prepared and

141

analyzed to check for contamination and carry-over5; peaks found in these blank were

142

removed the profiles obtained for the DOM samples33.

143

Molecular formulas of DOM were calculated using a custom software; the data

144

processing is described elsewhere13, 35. Briefly, detected mass peaks with S/N greater

145

than 6 were exported to datasheets for processing. Elemental combinations were

146

limited to molecular formulas containing

147

mass measurement errors less than 1 13. Each assigned class species was extended in

148

double bond equivalence (DBE) values (with a mass interval of 2 Da) and carbon

149

numbers (CH2 unit, with a mass interval of 14 Da) by an automatic retrieval within a

150

tolerance of ± 0.001 Kendrick mass defect (KMD)33. Through the application of these

151

criteria formulae were ascribed with a high level of confidence; and, only formulae

152

that were physically possible were permitted35. Element P were excluded because no

12

C0-100, 1H0-200,

10

ACS Paragon Plus Environment

14

N0-4,

16

O0-30,

32

S0-2 and

Page 11 of 40

Environmental Science & Technology

153

systematic trends in P-containing compounds were observed across the samples.

154

Where mass peaks were detected in more than two of the three replicates for

155

each sample they were retained for further data analysis17. Relative abundance of

156

moieties was determined through normalization of the sum of all FT-ICR MS signal

157

intensities of each mass of the respective DOM formulae. Numerous studies have

158

evidenced quantitative reproducibility of FT-ICR MS approaches36. Compound

159

groups were delineated by elemental ratios, aromaticity index (AI)11,

160

bond equivalence (DBE)40 and H/C cutoffs37: polycyclic aromatics originating from

161

combustion (AI > 0.66), polyphenols derived from plants (0.66 ≥ AI > 0.50), highly

162

unsaturated and phenolic compounds (AI ≤ 0.50 and H/C < 1.5), aliphatic compounds

163

(2.0 ≥ H/C ≥ 1.5)12, lipids (O/C = 0−0.3, H/ C = 1.5−2.0), protein/amino sugars (O/C

164

= 0.3−0.67, H/C = 1.5−2.2, N/C ≥ 0.05), carbohydrates (O/C = 0.67−1.2; H/C =

165

1.5−2), unsaturated hydrocarbons (O/C = 0−0.1, H/C = 0.7−1), lignins (O/C =

166

0.1−0.67, H/C = 0.7−1.5, AI < 0.67), condensed aromatics (O/C = 0−0.67, H/C =

167

0.2−0.7, AI ≥ 0.67) and tannins (O/C = 0.67−1.2, H/C = 0.5−1.5, AI < 0.67)41, 42.

168

Statistical analyses of DOM chemodiversity

37-39

, double

169

Moiety accumulation curve, ranked abundance and redundancy analysis were

170

performed in R software using the normalized intensities of all compounds12. Moiety

171

accumulation curve was calculated as the number of unique compounds by adding the

172

samples in a random order using the R package vegan. Confidence interval (95%) of

173

the molecule accumulation curve was estimated by 1,000 permutations. By this

174

approach a confidence interval for the number of unique compounds detected with 11

ACS Paragon Plus Environment

Environmental Science & Technology

175

each additional sample was obtained. Regarding the ranked abundance curves color

176

was used to depict the proportion (%) of samples in which a compound was acquired.

177

The rank abundance curve was obtained using the R package BiodiversityR43.

178

RDA was carried out using relative peak abundances (matrix: Bray-Curtis

179

dissimilarities) in R using package vegan44. Geography and soil chemistry variables

180

were fitted to each ordination using ‘envfit’ function implemented in package vegan,

181

with P-values reckoned over 999 permutations.

182

Correlation of DOM chemodiversity and environmental attributes

183

Soil chemistry associated molecular patterns were evaluated using Spearman’s

184

rank correlation coefficients. Where the formulae could be unambiguously assigned,

185

spearman correlation was performed between environmental parameters and the

186

sum-normalized intensity of peaks, following the method described previously12. The

187

environmental parameters covered soil chemistry (DOC, pH, total carbon, total

188

nitrogen, mineral elements (Fe, Mn, Al, Mg, Ca, Si)) and geography (elevation,

189

latitude and longitude). Van Krevelen45 diagrams were then plotted for each

190

environmental parameter, major patterns of chemodiversity were visually identified,

191

and variables were grouped based on patterns of similar correlation.

192

Reproducibility and duplicate analysis

193

The same amount of DOM was extracted from each sample in triplicates. All

194

samples were solid phase extracted following the same protocol in a random order

195

and analyzed on the same FT-ICR MS with the same settings, which has previously

196

shown high reproducibility of replicates12. Soil from MY and WC as representatives 12

ACS Paragon Plus Environment

Page 12 of 40

Page 13 of 40

Environmental Science & Technology

197

were extracted in triplicate for quality control (Table S6). The variation in the

198

replicates was negligible. Compared with applying presence/absence-based methods,

199

using normalized peak intensities broadened the range of statistical tools available for

200

interpreting molecular patterns of DOM46.

201

Index of chemodiversity

202

The chemodiversity index was calculated using the Chao 1 metric47 in R,

203

wherein the individual molecule was assigned as the “species” and the relative

204

intensity of the peaks used to ascribe species abundance12. Spearman’s correlation was

205

calculated between environmental parameters and Chao 1 index establish possible

206

associations between DOM chemodiversity, soil chemistry attributes and geography

207

variables. Automatic linear modeling was performed (95% confidence level) using

208

IBM SPSS Statistics for Mac (21.0)48. Finally, Spearman’s correlation was carried out

209

to identify molecules that were associated with Chao 1. These associations were

210

plotted in a van Krevelen diagrams to support further interpretation of the

211

relationships between DOM molecules and chemodiversity.

212

Linear discriminant analysis effect size (LEfSe)

213

DOM molecules established to be most likely to explain differences between

214

samples were determined by linear discriminant analysis effect size (LEfSe)17, 49. The

215

threshold on the LDA score for discriminative features was set at 2.0. Kruskal–Wallis

216

tests were performed at a significance level of P < 0.05

217

LEfSe was performed in the Galaxy framework49 with online analysis tools

218

(http://huttenhower.sph.harvard.edu/lefse/). 13

ACS Paragon Plus Environment

17

. An implementation of

Environmental Science & Technology

219

Results

220

Dominant DOM molecular compounds in paddy soils.

221

Individual molecules differed among samples with respect to their presence and

222

signal intensity indicating that paddy soil DOM molecules were observed to be very

223

complex and very diverse. Across the 16 paddy soils a total of 11,361 molecular

224

formulae were assigned. The identified formulae ranged in mass from 152.07 to

225

799.39 Da and had weighed average masses per sample between 312.1 and 355.3 Da

226

(average 320.8 Da). Statistical analyses of the classified groups for the formulae

227

showed that CHO (44.3%) and CHON (42.9%) molecules were dominant, followed

228

by CHOS (12.4%) (Table S4). N-containing compounds were present in paddy soil

229

DOM in large amounts; and, most moieties contained two or more N atoms. Van

230

Krevelen analysis revealed that paddy soil DOM has the greatest proportion of lignin

231

components and aromatic component; these can be ascribed as soil-derived

232

“humics”38. The relative abundance of lignin and aromatic components increased with

233

soil pH (Fig. 1).

234

Distribution of organic molecules in paddy soils across China.

235

Moiety accumulation curve and rank abundance were derived to characterize the

236

molecular distribution across the dataset. A saturation of compounds was observed in

237

molecular accumulation across the samples, with 95% of the molecular richness being

238

reached in 13 out of 16 paddy soils (Fig. 2a). The chemodiversity of paddy soil DOM

239

was captured within a relatively few number of sample sites. Accordingly, fewer

240

previously undetected compounds would be identified with each additional sample. 14

ACS Paragon Plus Environment

Page 14 of 40

Page 15 of 40

Environmental Science & Technology

241

Thus, we are confident that the variability of paddy soil DOM molecular composition

242

was well defined in this research.

243

Rank-abundance analysis is a well-known and informative method of analyzing

244

patterns of diversity43. Rank abundance for the DOM molecular profiles, for all

245

samples, showed that compounds with normalized intensities in the highest 30% were

246

present in at least 80% of the samples (Fig. 2b). This observation suggests that there

247

exists a fundamental foundation of DOM compounds that, while they varying in their

248

abundance, are ubiquitous across the large environmental gradient captured in this

249

research.

250

Linking DOM molecular profiles to environmental factors

251

Redundancy analysis (RDA) was used to evaluate relationships between the

252

intensities of all chemical moieties and environmental attributes (Fig. 3). Our results

253

indicated that DOM molecular profile composition depended, to a large extent, upon

254

soil parameters such as pH, soil minerals, carbon and nitrogen contents (Fig. 3). pH

255

emerged as the most important driver influencing the DOM compound profiles (Fig.

256

4a). Our analysis highlighted that soil pH resulted in divergent chemical composition.

257

Similarly, Fep/FeR and (Alp+Fep) also appeared as important drivers. This was

258

congruent with the opposing correlation observed with soil C/N ratio and NO3¯-N

259

concentration (Fig. 3).

260

Van Krevelen diagrams45, that cross-plot hydrogen to carbon (H/C) ratio with

261

oxygen to carbon (O/C) ratio of DOM moieties assist in the visualization of DOM

262

profile chemodiversity. The van Krevelen diagrams also facilitate the assignment of 15

ACS Paragon Plus Environment

Environmental Science & Technology

263

molecular formulae into different broad molecular groups, for example: aliphatics,

264

polyphenols and polycyclic aromatics. Our analysis, of van Krevelen plots (Fig. 4),

265

suggested that all DOM moieties were significantly correlated (P < 0.05) with either

266

pH, iron complexing index (Fep/FeR), or soil C/N ratio. Spearman rank,

267

sum-normalized, compound intensities were observed to correlate with environmental

268

variables, including pH, Fep/FeR, and C/N ratio. Across pH, Fep/FeR and C/N ratio,

269

clear trends were observed in the DOM moiety profiles (Fig. 4a-c). In general, the

270

relative abundance of polyphenols (derived from plants) and polycyclic aromatics

271

(derived from combustion) increased in prevalence in paddy soils with higher pH (Fig.

272

4a). Compounds with higher H/C ratio and O/C ratio were observed to be represented

273

in greater proportions in paddy soil with lower pH; a similar distribution pattern was

274

observed for soil C/N ratio, and; an inverse pattern observed for Fep/FeR (Fig. 4b-c).

275

DOM moieties that were positively correlated with pH had low H/C ratios and low

276

O/C ratios, and were clearly differentiated from moieties that were negatively

277

correlated with pH (Fig. 4a).

278

DOM molecular associations with soil pH

279

Spearman correlations between the intensity of individual moieties and

280

environmental factors (Fig. 4) highlighted that variability in the DOM chemodiversity

281

profiles was driven mainly by soil pH. Roughly half of all detected molecular

282

formulas (~4,904) varied significantly across paddy soils. Chemical moieties

283

positively correlated to pH included polycyclic aromatics (derived from combustion),

284

polyphenols (derived from plants) and carboxyl-containing aliphatic molecule. Seven 16

ACS Paragon Plus Environment

Page 16 of 40

Page 17 of 40

Environmental Science & Technology

285

different magnitude-weighted parameters (H/C)w, (O/C)w, (MW)w, (DBE)w,

286

(DBE/C)w, (DBE/O)w and (DBE-O)w18 were used to obtain more detailed

287

characterization of molecular chemodiversity profiles. The highest positive

288

correlation of pH with (DBE)w, (DBE/C)w, (DBE/O)w and (DBE-O)w, and the highest

289

negative correlation with (H/C)w and (O/C)W, indicated that at higher pH, the DOM

290

molecules were more aromatic (Table S5). The positive correlation between pH and

291

(DBE-O)w suggests that more C=C unsaturation existed with higher pH values.

292

Consequently, in more acidic paddy soils DOM displayed more carbon saturated and

293

N-containing compounds (e.g. fatty acids and proteinaceous materials), which were

294

more labile and inherently degradable.

295

Marker molecules of DOM in paddy soils

296

LEfSe was performed to determine markers in DOM composition under different

297

soil pH. LEfSe analysis pin-pointed 564 marker molecules that defined differences in

298

DOM composition at contrasting soil pH (Fig. 5). DOM in paddy soils of pH > 5.63

299

had 340 marker components compared to 224 for DOM in paddy soils of pH < 5.63

300

(Fig. 5). In general, at elemental composition level CHO compounds were dominant

301

in high pH group, whereas the CHON compounds are major components in low pH

302

group. At the compound level, DOM were marked by recalcitrant component (shown

303

in black) in high pH soil, such as polycyclic aromatics (derived from combustion) and

304

polyphenols (derived from plants), while labile compounds (shown in red), such as

305

proteins, amino sugars, carbohydrates, and phenolic compounds were marker

306

components in low pH soil. This result suggested that soil pH has a strong effect on 17

ACS Paragon Plus Environment

Environmental Science & Technology

307

the DOM molecular distribution. Most significantly, our findings suggest that DOM

308

in paddy soils will become more labile with further soil acidification.

309

DOM molecular associations with iron complexing index (Fep/FeR)

310

About 32.4% of the molecules that have been unambiguously identified showed

311

a significant correlation with iron complex index (Fep/FeR). For all there was

312

consistency in molecular-level chemodiversity patterns within DOM moieties, highly

313

unsaturated compounds and phenolic compounds (with formulae typical of aromatic

314

and condensed aromatic compounds), lignin and tannin, show a significant shift in

315

chemodiversity profiles associated with pH and Fep/FeR (Fig. 4a-b). Polycyclic

316

aromatics (derived from combustion) and polyphenols (derived from plants)

317

compounds were negatively correlated to Fep/FeR, with an inverse relationship with

318

pH. In addition, a large number of N-containing compounds, most containing two or

319

more N-atoms, were significantly correlated to iron complexing index.

320

DOM molecular associations with soil C/N ratio

321

Soil C/N ratio is a recognized indicator of soil quality56. Our results highlight

322

that soil C/N ratio has a close relationship with the DOM compounds and is also a

323

good predictor of the molecular-level DOM distribution (Fig. 4c). Molecules that

324

were positively correlated with C/N ratio formed a cluster and mainly plot in the

325

aromatic and condensed aromatics regions. Low soil C/N ratio was associated with

326

molecules with higher H/C ratios. Higher H/C ratios, indicative of greater saturation,

327

are characteristic for compounds that are labile (e.g. proteinaceous materials and

328

carbohydrates)37. A majority (69.4%) of nitrogenous compounds in our dataset were 18

ACS Paragon Plus Environment

Page 18 of 40

Page 19 of 40

Environmental Science & Technology

329

positively correlated with soil C/N ratio. This suggests that the inputs of organic

330

matter with high C/N ratios (e.g. charred residues, organic fertilizer and returned

331

straw) under intensive management favor accumulation of recalcitrant compounds,

332

and aromatic structures may prevail in the DOM of old paddies.

333

Chemodiversity of DOM in paddy soil

334

Chao 1 index was used to qualify DOM molecular chemodiversity in paddy soils

335

(Table S3). The chemodiversity ranged from 3,258 to 3,853, with an average diversity

336

index of 3,607, whereas the mean chemodiversity of sediment DOM ranged from

337

1,848

338

While DOM composition reflects variations in pH, Fep/FeR, and soil C/N ratio,

339

chemodiversity was primarily predicted by soil C/N ratio. Automatic linear modeling

340

indicated soil C/N ratio to be the most important predictor of DOM chemodiversity

341

(Fig. S2). Paddy soils established to have relatively low DOM chemodiversity were

342

noted to contained less aromatic but more aliphatic compounds (Fig. 4d). Compounds

343

established to have with a H/C ratio of less than 1 (primarily polyphenols and

344

polycyclic aromatics) were synonymous with higher chemodiversity. Through the

345

application of numerical analysis, more traditionally applied to biological species

346

biodiversity assessment (e.g. Chao 1 index)11, 12, we highlight the opportunity to apply

347

these methods to evaluate the chemodiversity profiles of DOM in soils across; and to

348

use this analysis to qualify the influence of broad environmental gradients on DOM

349

transformation.

350

Discussion

14

to 1,919 50, indicating a larger spatial heterogeneity of DOM in paddy soils.

19

ACS Paragon Plus Environment

Environmental Science & Technology

351

Here we present an original characterization of DOM chemodiversity in paddy

352

soil. The number of molecular formulae (11,361) was observed to be much higher in

353

paddy soil than in other environments, for example, agricultural soil DOM (1,237

354

formulae)9, sediment (1,848

355

and boreal rivers (4,109 formulae)51. With combination of multivariate analysis, we

356

demonstrated that the molecular distribution of paddy soil DOM was closely related

357

to soil pH, C/N ratio and iron complexing index (Table S7).

358

The effects of pH on paddy soil DOM distribution

14

to 1,919

50

formulae), boreal lakes (7,122 formulae)12,

359

It has been well documented that rice cultivation and increasing inputs of

360

chemical nitrogen fertilizers generally lead to more acidic and oligotrophic soil

361

conditions52-54, with changes in pH affecting organically complexed iron55, 56, along

362

with changes in soil C/N ratio56, ultimately affecting the distribution of DOM

363

chemical diversity. It has been reported that soil pH has decreased significantly in

364

major Chinese croplands over the last 25 years54, and specifically paddy soil has been

365

reported to have acidified (with an average pH decline of 0.59) over this period57.

366

Thus, with the pH decrease in paddy soils, the abundance of polyphenols (derived

367

from plant) and polycyclic aromatics (derived from combustion) decreased, while the

368

abundance of labile components such as aliphatic, highly unsaturated and phenolic

369

compounds have increased (Fig. 4a). Our results highlight the vulnerability of paddy

370

soil DOM to degradation, and, in particular, the antagonism acidification has with

371

respect to DOM degradation.

372

Iron-DOM complexing effects on paddy soil DOM distribution 20

ACS Paragon Plus Environment

Page 20 of 40

Page 21 of 40

Environmental Science & Technology

373

The importance of solid iron in the cycling of organic carbon has long been

374

known29, 58. However, direct field evidence linking the biogeochemical cycles of

375

DOM with iron, at the molecular level, has to date been lacking. Using FT-ICR MS,

376

we observed chemodiversity of DOM to be reduced as iron complexing index

377

(Fep/FeR) increased; our results further suggesting that organically complexed iron

378

selectively traps polyphenols (derived from plants) and polycyclic aromatics (derived

379

from combustion) (Fig. 4b). This observation could be explained by the fractionation

380

of DOM induced by iron oxyhydroxides. Iron oxyhydroxides complex with carboxyl

381

groups of DOM compounds21, 59 (Table S5, Spearman correlation: -0.347, P < 0.05),

382

resulting in the precipitation of phenolic and aromatic compounds with iron40 and

383

consequential release of hydrocarbons and amino acids into soil solutions. Iron is

384

ubiquitous in paddy soils61 and our results highlight its importance for SOM storage

385

as iron associated DOM. Our results also reveal that 32.4% of the DOM molecules in

386

paddy soils were significantly correlated with iron; this finding may have far-reaching

387

implications for the understanding of global cycling and the long-term storage of

388

SOM in paddy soils.

389

C/N ratio effects on paddy soil DOM distribution

390

Intensive inputs of chemical N fertilizer and simultaneous decrease in

391

applications of organic fertilizer have caused the decrease in soil C/N ratio in paddy

392

soils, as well as soil acidification. Paddy soils with high C/N ratio contain the most

393

diverse mixture of DOM molecules. DOM composition is more recalcitrant and

394

diverse at high soil C/N ratio. The abundance of polyphenols (derived from plants) 21

ACS Paragon Plus Environment

Environmental Science & Technology

395

and polycyclic aromatics (derived from combustion) decreased with decreasing soil

396

C/N ratio. The decreased abundance of polyphenols and polycyclic aromatics would

397

reduce the recalcitrant soil organic carbon and nitrogen pool; potentially leading to

398

soil degradation40 and a reduction in food security52. Furthermore, decreased

399

recalcitrant pyrogenic compounds could reduce the stability of soil organic carbon

400

pool; thus, potentially increase greenhouse gas emissions.

401

Paddy soils are of prominent importance to the global carbon cycling as they have

402

a pivotal role to play in carbon fixation, its respiration and methane emissions. Our

403

results have revealed that qualitative changes in DOM were mainly driven by soil pH,

404

iron bonding and C/N ratio in paddy soils. The applications of organic and chemical N

405

fertilizers significantly affect these factors, especially excess inputs of N fertilizer to

406

soil. These results present a new framework for understanding the distribution of

407

DOM in paddy soils and will facilitate predictions on how paddy soil DOM

408

chemodiversity and recalcitrance responds to soil management and environmental

409

gradients.

410

Notes

411

Declaration: We, the authors, have no competing financial interest linked to this

412

research and its findings.

413 414

Acknowledgements

415

This research was financially supported by the State Key Program of Natural Science

416

Foundation of China (No. 41430858), the National Natural Science Foundation of

417

China (No. 41371459) and the National Basic Research (973) Program of China (No. 22

ACS Paragon Plus Environment

Page 22 of 40

Page 23 of 40

Environmental Science & Technology

418

2014CB441102).

419 420

Supporting Information

421

Details on supportive methods and discussion; sampling sites, predictive importance

422

of chemical variables, schematic diagram for the preparation of SPE-DOM,

423

magnitude-weighted parameters, molecular characteristics, spearman correlation

424

analysis.

425 426

References

427

(1) van der Wal, A.; de Boer, W. Dinner in the dark: Illuminating drivers of soil

428

organic matter decomposition. Soil Biol. Biochem. 2017, 105, 45-48;

429

DOI:10.1016/j.soilbio.2016.11.006.

430

(2) Kalbitz, K.; Kaiser, K. Ecological aspects of dissolved organic matter in soils.

431

Geoderma 2003, 113 (3–4), 177-178; DOI:10.1016/S0016-7061(02)00359-2.

432

(3) Khan, S.; Chao, C.; Waqas, M.; Arp, H. P.; Zhu, Y. G. Sewage sludge biochar

433

influence upon rice (Oryza sativa L) yield, metal bioaccumulation and greenhouse gas

434

emissions from acidic paddy soil. Environ. Sci. Technol. 2013, 47 (15), 8624-8632;

435

DOI:10.1021/es400554x.

436

(4) Yan, X.; Akiyama, H.; Yagi, K.; Akimoto, H. Global estimations of the inventory

437

and mitigation potential of methane emissions from rice cultivation conducted using

438

the 2006 intergovernmental panel on climate change guidelines. Global Biogeochem.

439

Cy. 2009, 23 (2), 125-151; DOI:10.1029/2008gb003299. 23

ACS Paragon Plus Environment

Environmental Science & Technology

440

(5) Li, Y.; Xu, C.; Chung, K. H.; Shi, Q. Molecular characterization of dissolved

441

organic matter and its subfractions in refinery process water by Fourier transform ion

442

cyclotron resonance mass spectrometry. Energ. Fuel. 2015, 29 (5), 2923-2930;

443

DOI:10.1021/acs.energyfuels.5b00333.

444

(6) Nebbioso, A.; Piccolo, A. Molecular characterization of dissolved organic matter

445

(DOM): a critical review. Anal. Bioanal. Chem. 2013, 405 (1), 109-124;

446

DOI:10.1007/s00216-012-6363-2.

447

(7) Li, Y.; Fang, Z.; He, C.; Zhang, Y.; Xu, C.; Chung, K. H.; Shi, Q. Molecular

448

characterization and transformation of dissolved organic matter in refinery wastewater

449

from water treatment processes: characterization by Fourier transform ion cyclotron

450

resonance mass spectrometry. Energ. Fuel. 2015, 29 (11), 6956-6963;

451

DOI:10.1021/acs.energyfuels.5b01446.

452

(8) Hassett, J. P. Dissolved natural organic matter as a microreactor. Science 2006,

453

311 (5768), 1723-1724; DOI:10.1126/science.1123389.

454

(9) Ohno, T.; Ohno, P. E. Influence of heteroatom pre-selection on the molecular

455

formula assignment of soil organic matter components determined by ultrahigh

456

resolution mass spectrometry. Anal. Bioanal. Chem. 2013, 405 (10), 3299-3306;

457

DOI:10.1007/s00216-013-6734-3.

458

(10) Stubbins, A.; Spencer, R. G. M.; Chen, H.; Hatcher, P. G.; Mopper, K.; Hernes, P.

459

J.; Mwamba, V. L.; Mangangu, A. M.; Wabakanghanzi, J. N.; Six, J. Illuminated

460

darkness: Molecular signatures of Congo River dissolved organic matter and its

461

photochemical alteration as revealed by ultrahigh precision mass spectrometry. 24

ACS Paragon Plus Environment

Page 24 of 40

Page 25 of 40

Environmental Science & Technology

462

Limnol. Oceanogr. 2010, 55 (4), 1467-1477; DOI:10.4319/lo.2010.55.4.1467.

463

(11) Kellerman, A. M.; Kothawala, D. N.; Dittmar, T.; Tranvik, L. J. Persistence of

464

dissolved organic matter in lakes related to its molecular characteristics. Nat. Geosci.

465

2015, 8 (6), 454-457; DOI:10.1038/ngeo2440.

466

(12) Kellerman, A. M.; Dittmar, T.; Kothawala, D. N.; Tranvik, L. J. Chemodiversity

467

of dissolved organic matter in lakes driven by climate and hydrology. Nat. Commun.

468

2014, 5, 3804; DOI:10.1038/ncomms4804.

469

(13) Fang, Z.; He, C.; Li, Y.; Chung, K. H.; Xu, C.; Shi, Q. Fractionation and

470

characterization of dissolved organic matter (DOM) in refinery wastewater by revised

471

phase retention and ion-exchange adsorption solid phase extraction followed by ESI

472

FT-ICR MS. Talanta 2017, 162, 466-473; DOI:10.1016/j.talanta.2016.10.064.

473

(14) Zhang, L.; Wang, S.; Xu, Y.; Shi, Q.; Zhao, H.; Jiang, B.; Yang, J. Molecular

474

characterization of lake sediment WEON by Fourier transform ion cyclotron

475

resonance mass spectrometry and its environmental implications. Water Res. 2016,

476

106, 196-203; DOI:10.1016/j.watres.2016.09.059.

477

(15) Robert W. Kramer; Elizabeth B. Kujawinski; Hatcher, P. G. Identification of

478

black carbon derived structures in a volcanic ash soil humic acid by Fourier transform

479

ion cyclotron resonance mass spectrometry. Environ. Sci. Technol. 2004, 38 (12),

480

3387-3395; DOI:10.1021/es030124m.

481

(16) Osterholz, H.; Singer, G.; Wemheuer, B.; Daniel, R.; Simon, M.; Niggemann, J.;

482

Dittmar, T. Deciphering associations between dissolved organic molecules and

483

bacterial communities in a pelagic marine system. ISME J. 2016, 10 (7), 1717–1730; 25

ACS Paragon Plus Environment

Environmental Science & Technology

484

DOI:10.1038/ismej.2015.231.

485

(17) Lucas, J.; Koester, I.; Wichels, A.; Niggemann, J.; Dittmar, T.; Callies, U.;

486

Wiltshire, K. H.; Gerdts, G. Short-term dynamics of north sea

487

bacterioplankton-dissolved organic matter coherence on molecular level. Front.

488

Microbiol. 2016, 7, 321; DOI:10.3389/fmicb.2016.00321.

489

(18) D'Andrilli, J.; Dittmar, T.; Koch, B. P.; Purcell, J. M.; Marshall, A. G.; Cooper, W.

490

T. Comprehensive characterization of marine dissolved organic matter by Fourier

491

transform ion cyclotron resonance mass spectrometry with electrospray and

492

atmospheric pressure photoionization. Rapid Commun. Mass Spectrom. 2010, 24 (5),

493

643-650; DOI: 10.1002/rcm.4421.

494

(19) Roth, V. N.; Dittmar, T.; Gaupp, R.; Gleixner, G. The molecular composition of

495

dissolved organic matter in forest soils as a function of pH and temperature. PLoS

496

One 2015, 10 (3), e0119188; DOI:10.1371/journal.pone.0119188.

497

(20) Hertkorn, N.; Harir, M.; Koch, B. P.; Michalke, B.; Schmitt-Kopplin, P.

498

High-field NMR spectroscopy and FTICR mass spectrometry: powerful discovery

499

tools for the molecular level characterization of marine dissolved organic matter.

500

Biogeosciences 2013, 10 (3), 1583-1624; DOI:10.5194/bg-10-1583-2013.

501

(21) Riedel, T.; Zak, D.; Biester, H.; Dittmar, T. Iron traps terrestrially derived

502

dissolved organic matter at redox interfaces. Proc. Natl. Acad. Sci. USA 2013, 110

503

(25), 10101-10105; DOI:10.1073/pnas.1221487110.

504

(22) Ohno, T.; Parr, T. B.; Gruselle, M. C.; Fernandez, I. J.; Sleighter, R. L.; Hatcher,

505

P. G. Molecular composition and biodegradability of soil organic matter: a case study 26

ACS Paragon Plus Environment

Page 26 of 40

Page 27 of 40

Environmental Science & Technology

506

comparing two new England forest types. Environ. Sci. Technol. 2014, 48 (13),

507

7229-7236; DOI:10.1021/es405570c.

508

(23) FAOSTAT. FAO Statistical databases. Food and Agriculture Organization (FAO)

509

of the United Nations. http://www.fao.org/ (accessed 15 September, 2017).

510

(24) Gregorich, M. R. C. a. E. G., Soil Sampling and Methods of Analysis, Second

511

Edition. CRC Press: 2007.

512

(25) Wang, N.; Chen, Z.; Li, H. B.; Su, J. Q.; Zhao, F.; Zhu, Y. G. Bacterial

513

community composition at anodes of microbial fuel cells for paddy soils: the effects

514

of soil properties. J. Soils Sed. 2015, 15 (4), 926-936;

515

DOI:10.1007/s11368-014-1056-4.

516

(26) Yang, X. R.; Li, H.; Nie, S. A.; Su, J. Q.; Weng, B. S.; Zhu, G. B.; Yao, H. Y.;

517

Gilbert, J. A.; Zhu, Y. G.; Nojiri, H. Potential contribution of anammox to nitrogen

518

loss from paddy soils in southern China. Appl. Environ. Microbiol. 2015, 81 (3),

519

938-947; DOI:10.1128/aem.02664-14.

520

(27) Wang, H. Y.; Wen, S. L.; Chen, P.; Zhang, L.; Cen, K.; Sun, G. X. Mitigation of

521

cadmium and arsenic in rice grain by applying different silicon fertilizers in

522

contaminated fields. Environ. Sci. Pollut. Res. Int. 2016, 23 (4), 3781-3788;

523

DOI:10.1007/s11356-015-5638-5.

524

(28) Sun, G. X.; Liu, X.; Williams, P. N.; Zhu, Y. G. Distribution and translocation of

525

selenium from soil to grain and its speciation in paddy rice (Oryza sativa L.). Environ.

526

Sci. Technol. 2010, 44 (17), 6706-6711; DOI:10.1021/es101843x.

527

(29) Lalonde, K.; Mucci, A.; Ouellet, A.; Gelinas, Y. Preservation of organic matter in 27

ACS Paragon Plus Environment

Environmental Science & Technology

528

sediments promoted by iron. Nature 2012, 483 (7388), 198-200;

529

DOI:10.1038/nature10855.

530

(30) Dittmar, T.; Koch, B.; Hertkorn, N.; Kattner, G. A simple and efficient method for

531

the solid-phase extraction of dissolved organic matter (SPE-DOM) from seawater.

532

Limnol. Oceanogr. Methods 2008, 6 (6), 230-235; DOI:10.4319/lom.2008.6.230.

533

(31) Guigue, J.; Mathieu, O.; Lévêque, J.; Mounier, S.; Laffont, R.; Maron, P. A.;

534

Navarro, N.; Chateau, C.; Amiotte-Suchet, P.; Lucas, Y. A comparison of extraction

535

procedures for water-extractable organic matter in soils. Eur. J. Soil Sci. 2014, 65 (4),

536

520-530; DOI:10.1111/ejss.12156.

537

(32) Xiao, J.; He, X.; Hao, J.; Zhou, Y.; Zheng, L.; Ran, W.; Shen, Q.; Yu, G. New

538

strategies for submicron characterization the carbon binding of reactive minerals in

539

long-term contrasting fertilized soils: implications for soil carbon storage.

540

Biogeosciences 2016, 13 (12), 3607-3618; DOI:10.5194/bg-13-3607-2016.

541

(33) Zhang, H.; Zhang, Y.; Shi, Q.; Ren, S.; Yu, J.; Ji, F.; Luo, W.; Yang, M.

542

Characterization of low molecular weight dissolved natural organic matter along the

543

treatment trait of a waterworks using Fourier transform ion cyclotron resonance mass

544

spectrometry. Water Res. 2012, 46 (16), 5197-5204;

545

DOI:10.1016/j.watres.2012.07.004.

546

(34) Ohno, T.; He, Z.; Sleighter, R. L.; Honeycutt, C. W.; Hatcher, P. G. Ultrahigh

547

resolution mass spectrometry and indicator species analysis to identify marker

548

components of soil-and plant biomass-derived organic matter fractions. Environ. Sci.

549

Technol. 2010, 44 (22), 8594-8600; DOI:10.1021/es101089t. 28

ACS Paragon Plus Environment

Page 28 of 40

Page 29 of 40

Environmental Science & Technology

550

(35) Wang, X.; Wang, J.; Zhang, Y.; Shi, Q.; Zhang, H.; Zhang, Y.; Yang, M.

551

Characterization of unknown iodinated disinfection byproducts during

552

chlorination/chloramination using ultrahigh resolution mass spectrometry. Sci. Total

553

Environ. 2016, 554-555, 83-88; DOI:10.1016/j.scitotenv.2016.02.157.

554

(36) Seidel, M.; Yager, P. L.; Ward, N. D.; Carpenter, E. J.; Gomes, H. R.; Krusche, A.

555

V.; Richey, J. E.; Dittmar, T.; Medeiros, P. M. Molecular-level changes of dissolved

556

organic matter along the Amazon River-to-ocean continuum. Mar. Chem. 2015, 177,

557

218-231; DOI:10.1016/j.marchem.2015.06.019.

558

(37) Koch, B. P.; Dittmar, T. From mass to structure: an aromaticity index for

559

high-resolution mass data of natural organic matter. Rapid Commun. Mass Spectrom.

560

2006, 20 (5), 926-932; DOI:10.1002/rcm.2386.

561

(38) Santl-Temkiv, T.; Finster, K.; Dittmar, T.; Hansen, B. M.; Thyrhaug, R.; Nielsen,

562

N. W.; Karlson, U. G. Hailstones: a window into the microbial and chemical inventory

563

of a storm cloud. PLoS One 2013, 8 (1), e53550; DOI:10.1371/journal.pone.0053550.

564

(39) Dittmar, T.; Koch, B. P. Thermogenic organic matter dissolved in the abyssal

565

ocean. Mar. Chem. 2006, 102 (3-4), 208-217; DOI:10.1016/j.marchem.2006.04.003.

566

(40) Lv, J.; Zhang, S.; Wang, S.; Luo, L.; Cao, D.; Christie, P. Molecular-scale

567

investigation with ESI-FT-ICR-MS on fractionation of dissolved organic matter

568

induced by adsorption on iron oxyhydroxides. Environ. Sci. Technol. 2016, 50 (5),

569

2328-2336; DOI:10.1021/acs.est.5b04996.

570

(41) Feng, L.; Xu, J.; Kang, S.; Li, X.; Li, Y.; Jiang, B.; Shi, Q. Chemical composition

571

of microbe-derived dissolved organic matter in cryoconite in Tibetan Plateau glaciers: 29

ACS Paragon Plus Environment

Environmental Science & Technology

572

insights from Fourier transform ion cyclotron resonance mass spectrometry analysis.

573

Environ. Sci. Technol. 2016, 50 (24), 13215-13223; DOI:10.1021/acs.est.6b03971.

574

(42) Antony, R.; Grannas, A. M.; Willoughby, A. S.; Sleighter, R. L.; Thamban, M.;

575

Hatcher, P. G. Origin and sources of dissolved organic matter in snow on the East

576

Antarctic ice sheet. Environ. Sci. Technol. 2014, 48 (11), 6151-6159;

577

DOI:10.1021/es405246a.

578

(43) Kindt, R.; Coe, R., Tree diversity analysis: a manual and software for common

579

statistical methods for ecological and biodiversity studies. World Agroforestry Centre:

580

2005; Vol. 235, p 368-374.

581

(44) Legendre, P.; Gallagher, E. D. Ecologically meaningful transformations for

582

ordination of species data. Oecologia 2001, 129 (2), 271-280;

583

DOI:10.1007/s004420100716.

584

(45) Kim, S.; Kramer, R. W.; Hatcher, P. G. Graphical method for analysis of

585

ultrahigh-resolution broadband mass spectra of natural organic matter, the van

586

Krevelen diagram. Anal. Chem. 2003, 75 (20), 5336-5344; DOI:10.1021/ac034415p.

587

(46) Roth, V.-N.; Dittmar, T.; Gaupp, R.; Gleixner, G. Latitude and pH driven trends in

588

the molecular composition of DOM across a north south transect along the Yenisei

589

River. Geochim. Cosmochim. Acta 2013, 123, 93-105;

590

DOI:10.1016/j.gca.2013.09.002.

591

(47) Chao, A., and Tsung-Jen Shen. Nonparametric estimation of Shannon’s index of

592

diversity when there are unseen species in sample. Environ. Ecol. Stat. 2003, 10 (4),

593

429-443; DOI:10.1023/A:1026096204727. 30

ACS Paragon Plus Environment

Page 30 of 40

Page 31 of 40

Environmental Science & Technology

594

(48) Sun, R.; Zhang, X.-X.; Guo, X.; Wang, D.; Chu, H. Bacterial diversity in soils

595

subjected to long-term chemical fertilization can be more stably maintained with the

596

addition of livestock manure than wheat straw. Soil Biol. Biochem. 2015, 88, 9-18;

597

DOI:10.1016/j.soilbio.2015.05.007.

598

(49) Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W. S.;

599

Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome. Biol.

600

2011, 12 (6), R60; DOI:10.1186/gb-2011-12-6-r60.

601

(50) He, W.; Chen, M.; Park, J.-E.; Hur, J. Molecular diversity of riverine

602

alkaline-extractable sediment organic matter and its linkages with spectral indicators

603

and molecular size distributions. Water Res. 2016, 100, 222-231;

604

DOI:10.1016/j.watres.2016.05.023.

605

(51) Stubbins, A.; Lapierre, J. F.; Berggren, M.; Prairie, Y. T.; Dittmar, T.; del Giorgio,

606

P. A. What's in an EEM? Molecular signatures associated with dissolved organic

607

fluorescence in boreal Canada. Environ. Sci. Technol. 2014, 48 (18), 10598-10606;

608

DOI:10.1021/es502086e.

609

(52) Wang, C.; Li, W.; Yang, Z.; Chen, Y.; Shao, W.; Ji, J. An invisible soil

610

acidification: Critical role of soil carbonate and its impact on heavy metal

611

bioavailability. Sci. Rep. 2015, 5, 12735; DOI:10.1038/srep12735.

612

(53) Tian, D.; Niu, S. A global analysis of soil acidification caused by nitrogen

613

addition. Environ. Res. Lett. 2015, 10 (2), 024019;

614

DOI:10.1088/1748-9326/10/2/024019.

615

(54) Guo, J. H.; Liu, X. J.; Zhang, Y.; Shen, J. L.; Han, W. X.; Zhang, W. F.; Christie, 31

ACS Paragon Plus Environment

Environmental Science & Technology

616

P.; Goulding, K. W. T.; Vitousek, P. M.; Zhang, F. S. Significant acidification in major

617

Chinese croplands. Science 2010, 327 (5968), 1008-1010;

618

DOI:10.1126/science.1182570.

619

(55) Kogel-Knabner, I.; Amelung, W.; Cao, Z. H.; Fiedler, S.; Frenzel, P.; Jahn, R.;

620

Kalbitz, K.; Kolbl, A.; Schloter, M. Biogeochemistry of paddy soils. Geoderma 2010,

621

157 (1-2), 1-14; DOI:10.1016/j.geoderma.2010.03.009.

622

(56) Kalbitz, K.; Solinger, S.; Park, J.-H.; Michalzik, B.; Matzner, E. Controls on the

623

dynamics of dissolved organic matter in soils: a review. Soil. Sci. 2000, 165 (4),

624

277-304; DOI:10.1097/00010694-200004000-00001.

625

(57) Zhou, X. Y.; Zhou, S. W.; Xu, M. G.; Gilles, C. Evolution characteristics and

626

influence factors of acidification in paddy soil of southern China. Scientia Agricultura

627

Sinica 2015, 48 (23), 4811-4817; DOI:10.3864/j.issn.0578-1752.2015.23.022.

628

(58) Kaiser, K.; Guggenberger, G. The role of DOM sorption to mineral surfaces in the

629

preservation of organic matter in soils. Org. Geochem. 2000, 31 (7-8), 711-725;

630

DOI:10.1016/s0146-6380(00)00046-2.

631

(59) Henneberry, Y. K.; Kraus, T. E. C.; Nico, P. S.; Horwath, W. R. Structural stability

632

of coprecipitated natural organic matter and ferric iron under reducing conditions.

633

Org. Geochem. 2012, 48, 81-89; DOI:10.1016/j.orggeochem.2012.04.005.

634

(60) Ding, L. J.; Su, J. Q.; Xu, H. J.; Jia, Z. J.; Zhu, Y. G. Long-term nitrogen

635

fertilization of paddy soil shifts iron-reducing microbial community revealed by

636

RNA-13C-acetate probing coupled with pyrosequencing. ISME J. 2014, 9 (3),

637

721-734; DOI:10.1038/ismej.2014.159. 32

ACS Paragon Plus Environment

Page 32 of 40

Page 33 of 40

Environmental Science & Technology

638

(61) Ding, L. J.; An, X. L.; Li, S.; Zhang, G. L.; Zhu, Y. G. Nitrogen loss through

639

anaerobic ammonium oxidation coupled to iron reduction from paddy soils in a

640

chronosequence. Environ. Sci. Technol. 2014, 48 (18), 10641-10647;

641

DOI:10.1021/es503113s.

642

33

ACS Paragon Plus Environment

Environmental Science & Technology

643

Figure legends

644

Figure 1. Comparison of relative abundance (%) of van Krevelen diagram-derived

645

classification classes from the FT-ICR MS analysis of the DOM components for 16

646

paddy soils.

647

Figure 2. Molecular distributions of FT-ICR MS assigned molecules in sixteen DOM

648

samples across China. (a) Number of unique molecules with each added sample. The

649

black dotted line indicates 95% of compounds. (b) Rank abundance of the compounds

650

across all paddy soil DOM samples. Molecular compounds are color coded by the

651

percentage of samples in which they occurred.

652

Figure 3. Multivariate analysis of molecular data and drivers using RDA. Ordinations

653

are based on Bray-Curtis, which utilizes relative compound abundances information.

654

Environmental variables were fit to the ordination. Grey-shaded circles are DOM

655

compounds, whereas black circles indicate the site. Variables with a significance level

656

of < 0.05 (blue), < 0.01 (light green) and < 0.001 (purple) are shown. FeR, the

657

citrate-bicarbonate-dithionite extracted reactive Fe. FeP, the sodium pyrophosphate

658

extracted organic complexed Fe. FeP/FeR, this ratio represents the iron complexing

659

index. (Alp+ Fep), the sum of organic complexed Al and Fe.

660

Figure 4. Molecular-level DOM patterns across 16 paddy fields. Significant

661

Spearman rank correlation coefficients (P < 0.05) of individual molecules with (a) pH,

662

(b) Fep/FeR, (c) C/N and (d) the chemodiversity index (Chao 1). The color scale

663

indicates Spearman correlations between the intensity of individual molecules and pH,

664

Fep/FeR, C/N and chemodiversity index (red, positive; blue, negative). Circles indicate 34

ACS Paragon Plus Environment

Page 34 of 40

Page 35 of 40

Environmental Science & Technology

665

compounds without N and diamonds indicate N-containing compounds. Compound

666

groups include combustion-derived polycyclic aromatics (Aromaticity index (AI) >

667

0.66), vascular plant-derived polyphenols (0.66 ≥ AI > 0.50), highly unsaturated and

668

phenolic compounds (AI ≤ 0.50 and H/C < 1.5, and aliphatic compounds (2.0 ≥ H/C ≥

669

1.5). Compound category labels for delineation in panels (d) also apply to delineated

670

regions in (a) (c) and (d). Lines separating compound categories on van Krevelen

671

diagrams are for visualization only and exact categorization may slightly differ. The

672

number of significant correlations is (a) n=4,904; (b) n=3,676; (c) n=2,498; (d)

673

n=2,714

674

Figure 5. Linear discriminant effect size analysis (LEfSe) results on DOM molecules.

675

van Krevelen of linear discriminant analysis (LDA) scores computed for DOM

676

molecules, differently abundant across paddy soils. Data are displayed as molecular

677

hydrogen-to-carbon versus oxygen-to-carbon ratio. Circles indicate individual

678

molecule compounds. Red indicate the molecules that are enriched in paddy soils (pH

679

< 5.63) and black indicate the molecules that are enriched in paddy soils (pH > 5.63).

680

Molecule compounds with no significant differences are not shown.

681 682

35

ACS Paragon Plus Environment

Relative abundance(%)

100

Environmental Science & Technology

Page 36 of 40

80 60 40 20 0 TY AH YT

LZ JSJ JX WC JS

Tannins Ligins Carbohydrates Lipids

HL SY

BJ DH MY AX

JL

Condensed Aromatics Unsaturated Hydrocarbons Protein ACS Paragon Plus Environment

BC

Page 37 of 40

Environmental Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

ACS Paragon Plus Environment

Page 38 of 40

Page 39 of 40

Environmental Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

F ig u r e 5

Page 40 of 40

2 .5

p H < 5 .6 3 n = 2 2 4 p H > 5 .6 3 n = 3 4 0

H /C

2 .0 A lip h a tic

1 .5

H ig h ly u n s a tu r a te d a n d p h e n o lic c o m p o u n d s

1 .0

V a s c u la r p la n t- d e r iv e d p C o m b u o ly p h e n s tio n - d o ls e r iv e d p o ly c y c lic a r o m a tic s

0 .5

0 .0 0 .0

0 .2

0 .4

0 .6

O /C

ACS Paragon Plus Environment

0 .8

1 .0

1 .2