Effect of bioaugmentation and biostimulation on hydrocarbon degradation and microbial community composition in petroleum-contaminated loessal soil

Effect of bioaugmentation and biostimulation on hydrocarbon degradation and microbial community composition in petroleum-contaminated loessal soil

Journal Pre-proof Effect of bioaugmentation and biostimulation on hydrocarbon degradation and microbial community composition in petroleum-contaminate...

969KB Sizes 3 Downloads 112 Views

Journal Pre-proof Effect of bioaugmentation and biostimulation on hydrocarbon degradation and microbial community composition in petroleum-contaminated loessal soil

Manli Wu, Jialuo Wu, Xiaohui Zhang, Xiqiong Ye PII:

S0045-6535(19)31680-7

DOI:

https://doi.org/10.1016/j.chemosphere.2019.124456

Article Number:

124456

Reference:

CHEM 124456

To appear in:

Chemosphere

Received Date:

07 January 2019

Accepted Date:

25 July 2019

Please cite this article as: Manli Wu, Jialuo Wu, Xiaohui Zhang, Xiqiong Ye, Effect of bioaugmentation and biostimulation on hydrocarbon degradation and microbial community composition in petroleum-contaminated loessal soil, Chemosphere (2019), https://doi.org/10.1016/j. chemosphere.2019.124456

This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2019 Published by Elsevier.

Journal Pre-proof 1

Code number: CHEM59490

2 3

Effect of bioaugmentation and biostimulation on hydrocarbon

4

degradation

5

petroleum-contaminated loessal soil

6

Manli Wua,b, Jialuo Wua,b, Xiaohui Zhangc, Xiqiong Yea,b

7

(a Shaanxi Key Laboratory of Environmental Engineering, Xi’an University of Architecture and

8

Technology, Xi’an 710055, People’s Republic of China;

9

b Key Laboratory of Northwest Water Resource, Environment and Ecology, MOE, Xi’an

and

microbial

community

composition

in

10

University of Architecture and Technology, Xi’an 710055, People’s Republic of China;

11

c Exploration and Development Research Institute of Changqing Oilfield Company, PetroChina,

12

Xi’an 710018, People’s Republic of China)

13

Abstract: This study assessed the benefits of biostimulation with nitrogen and phosphorous

14

(BS) versus bioaugmentation with native petroleum degrading flora (BA) in terms of petroleum

15

hydrocarbon removal and the microbial community structure shift in petroleum-polluted loessal

16

soil. After 12 weeks of remediation, the TPH degradation efficiencies were 28.3% and 13.9% in

17

BS and BA treated soils, respectively. Biostimulation was more effective than bioaugmentation

18

for petroleum hydrocarbon degradation. Soil microbial community composition changed while

19

microbial diversity decreased greatly by bioaugmentation treatment. The inoculum could

20

survive, grow up quickly and become the predominant microorganisms after one week of

21

inoculation. In the biostimulation treatment, microbial community composition is more

22

evenness and richness than in the bioaugmented remediation. The strong positive correlations of

23

the nitrogen and phosphorus with the petroleum hydrocarbon suggest the importance of

24

nutrients for petroleum biodegradation processes in the contaminated loessal soil. The results 1

Journal Pre-proof 25

indicate that the stabilization and variety of the microbial community structure are essential for

26

the petroleum biodegradation performance. Further engineering is suggested to improve the

27

evenness and richness of the soil microbial community since an abundance of nitrogen and

28

phosphorus nutrients ensures the degraders’ activity in the petroleum polluted environment.

29

Keywords: Petroleum-contaminated soil; Bioaugmentation; Biostimulation; TPH; Microbial

30

community; α-diversity

31 32

Introduction

33

Petroleum pollution has become a global environmental concern. Soil contamination with

34

petroleum hydrocarbons is caused by accidental leakage from reservoirs and petroleum

35

refineries and spillage from transport pipelines. The presence of petroleum hydrocarbons affects

36

the physical, chemical and ecological properties of soil (Ramadass et al., 2015; Petrov et al.,

37

2016). Among a variety of the remediation strategies, bioremediation, based on the use of

38

microorganisms to degrade contaminants, has been widely used for treating oil-contaminated

39

sites and is an efficient, cost-effective and environmentally sound remediation approach.

40

Bioaugmentation and biostimulation are the two main bioremediation strategies for the

41

decontamination of petroleum polluted soil. Biostimulation enhances the metabolic activity of

42

the indigenous microbial community through nutrients amendment. Nitrogen and phosphorus

43

are the essential growth-limiting nutrients for microorganism growth. Many studies suggest that

44

the supplement of nitrogen and phosphorus is important to enhance the degradation of petroleum

45

in polluted soil (Leys et al., 2005; John et al., 2011; Shahi et al., 2016). Studies have indicated

46

that when the C/N/P ratio is regulated to 100/5/1, 100/10/1, and 100/15/1, the bacterial activity

47

in biostimulation practice is tolerant to various hydrocarbons and can utilize hydrocarbons as

48

carbon sources for their growth (Shahi et al., 2016). Since different properties of contaminated

49

soil has a diverse microbial community, the optimum C/N/P ratio may be different for 2

Journal Pre-proof 50

remediation. Generally, a ratio of C/N/P of 100/10/1 is widely accepted as a standard formula

51

for biostimulation strategy (Leys et al., 2005; Shahi et al., 2016).

52

Bioaugmentation, the inoculation of exogenous microorganisms into soil, is the best option

53

when contaminated soils have low indigenous populations of hydrocarbon degraders (Xu and

54

Lu, 2010; Castiglione et al., 2016). In addition, bioaugmentation gives an appropriate result

55

when the pollutants have a toxic effect on the native microorganisms (Tyagi et al., 2011; Wu et

56

al., 2013). One of the advantages of bioaugmentation is that the degradation process can start

57

immediately when specific microbial degraders are introduced.

58

Although bioaugmentation technology has been reported to be a sufficient method to

59

improve the degradation of petroleum in contaminated soil, the effects of bioaugmentation are

60

case specific and controversial. Some studies indicate that bioaugmentation only enhances

61

degradation efficiencies temporarily and that biostimulation represents a more ideal remediation

62

strategy to achieve effective petroleum decontamination (Kauppi et al., 2011; Liu et al., 2011;

63

Megharaj et al., 2011; Sayara et al., 2011; Abed et al., 2014; Yang et al., 2015; Wu et al., 2016;

64

Yulia et al., 2017). However, the reason for the inefficiency of bioaugmentation treatment in

65

some cases is still not apparent. There is a need to better understand the causes that lead to

66

bioaugmentation failure.

67

Microbial communities, as the essential component in soil ecological systems, are vital for

68

the remediation of petroleum hydrocarbon and other pollutants (Yan et al., 2016; Jia et al.,

69

2017). Since petroleum hydrocarbon remediation depends on soil microorganisms that have

70

suitable metabolic capabilities, understanding how bioaugmentation and biostimulation

71

influence the soil microbial community structure are fundamental to ensuring effective

72

remediation of oil-contaminated soil.

73

Molecular technologies provide a good understanding of the microbial community

74

composition in remediation systems (Zhang, et al., 2018). Various studies have attempted to 3

Journal Pre-proof 75

analyze the microbial structure that is responsible for the degradation of petroleum

76

hydrocarbons. Studies documenting the remediation of petroleum contaminated soil have

77

applied denaturing gel gradient electrophoresis (DGGE) (Jia et al., 2005; Pacwa-Płociniczak et

78

al., 2016; Wu et al., 2016), phospholipid fatty acid (PLFA) (Bundy et al., 2004; Margesin et al.,

79

2007), terminal-restriction fragment length polymorphism percentage(T-RFLP) (Grant et al.,

80

2007; Liu et al., 2011; Liu et al., 2013) and high throughput sequencing (Peng et al., 2015; Bao

81

et al., 2017; Wu et al., 2017a) to investigate bacterial community shifts. These

82

culture-independent molecular biology techniques can provide a good understanding of the

83

microbial community shift during remediation.

84

This paper presents a study of the effects of bioaugmentation and biostimulation strategies

85

applied to loessal soil polluted with 20,000 mg kg-1 total petroleum hydrocarbons (TPH). The

86

objectives were 1) to assess the petroleum hydrocarbon degradation efficiencies of the

87

bioaugmentation and biostimulation strategies, 2) to investigate the changes in the microbial

88

community structures as a function of time and type of treatment using high throughput

89

sequencing analysis, and 3) to evaluate the potential benefits of biostimulation with nutrients

90

versus bioaugmentation with hydrocarbon degrading microorganisms. We present evidence

91

regarding the superiority of the biostimulation strategy for the remediation of petroleum

92

contaminated loessal soil.

93

2. Materials and methods

94

2.1. Soil for the bioremediation study

95

A petroleum polluted soil sample was collected from an oil well in north Shaanxi

96

province, China. The site has a history of oil pollution over many years. According to the

97

sampling method “The Technical Specification for Soil Environmental monitoring (HJ/T 166

98

-2004, China)”, eight soil samples were sampled using a sterilized shovel from the top 10-cm

99

soil layer around the oil well. Then, these samples were well-mixed and homogenized into a 4

Journal Pre-proof 100

single soil sample. The fresh samples were loaded in polyethylene bag and immediately

101

transported to the laboratory by using a cooler. Table 1 summarizes selected physicochemical

102

and biological characteristics of this contaminated soil (IS).

103

2.2. Enrichment and identification of petroleum degrading flora from petroleum polluted soils

104

Five grams of petroleum-contaminated soil was suspended in 50 mL of MBS culture

105

medium. The medium contained 1% petroleum as the sole carbon and energy source. Incubation

106

was cultured at 25 ℃ by shaking on an oscillating incubator (150 rpm) for 7 days, 6% of the

107

cultures were then transferred into 50 mL of fresh MBS for further enriched. After five cycles of

108

repeated enrichment, the petroleum hydrocarbon degrading microbial flora PMC was obtained

109

by centrifugation (5000 rpm) for 15 min and suspended in PBS buffer. Identification of the

110

hydrocarbon degrading flora was performed using high throughput sequencing analysis by

111

Sangon Biotech Co., Ltd. China (ftp://ftp.sangon.com:21148).

112

2.3. Bioremediation experiment design

113

Two different remediation microcosms of BA and BS were prepared. (1) BA:

114

bioaugmentation was performed by adding the PMC hydrocarbon degrading flora (described in

115

2.2) into the petroleum-contaminated soil to obtain a density of 108 cfu g-1 soil as reported by

116

Abalos et al. (2004). (2) BS: biostimulation with NH4NO3 and KH2PO4, which were amended to

117

soil by a C/N/P ratio of 100/10/1.

118

Each microcosm was prepared in triplicate with 0.8 kg of soil. The microcosms were then

119

incubated at 24 °C for 12 weeks. Distilled water was amended to the soil to keep 15% soil

120

moisture. The soil was stirred weekly for sufficient aeration.

121

2.4. Physicochemical characteristics analysis

122

The soil pH, available phosphorus (AP), total nitrogen (TN), ammonia nitrogen (NH4+-N),

123

nitrate nitrogen (NO3--N), and oxidation-reduction potential (ORP) were determined according

124

to the technical specifications for soil analysis (2014). Total petroleum hydrocarbons (TPH) 5

Journal Pre-proof 125

were determined by the gravimetric method as described in Mishra et al. (2001).

126

2.5. Illumina sequencing

127

The initial petroleum contaminated soil (IS) and soil samples from the BA and BS

128

microcosms were collected at weeks 1 and 12 for MiSeq sequencing. DNA was extracted

129

according to a Power Soil DNA extraction kit (MoBio Laboratories, USA). The quality of DNA

130

was determined on a 0.5 % agarose gel. For bacterial community structure analysis, primers

131

341F

132

(GACTGGAGTTCCTTGGCACCCGAGAATTCCAGACTACHVGGGTATCTAATCC)

133

used to amplify the the V3-V4 region as reported by Qu et al. (2016). The PCR products then

134

loaded on Illumina-Miseq device according to the manufacturer’s protocols. After sequencing,

135

data were collected and processed by Sangon Biotech Co., Ltd (ftp://ftp.sangon.com:21148).

136

Operational taxonomic unit (OTU) clustering was performed using Usearch (version 8.1.1831)

137

at 97% sequence similarity threshold.

138

2.6. Statistical analysis

(CCCTACACGACGCTCTTCCGATCTGCCTACGGGNGGCWGCAG)

and

805R was

139

All experiments were conducted with three replications. The experimental results of the

140

TPH content, soil pH, AP, TN, NH4+-N, NO3--N, and ORP are presented as the mean ±standard

141

deviation (SD). To study the relationship between petroleum hydrocarbon degradation and soil

142

physicochemical characteristics, the Pearson correlation coefficient was calculated using

143

statistical software SPSS 25 (Statistical Package for the Social Sciences, China). A p value less

144

than 0.05 was considered to be statistically significant.

145

3. Results

146

3.1. Analysis of petroleum-degrading flora

147

The composition of petroleum-degrading flora PMC isolated from petroleum contaminated

148

soil was analyzed according to the high throughput sequencing method. The petroleum

149

degrading consortia were classified into two phyla and sixteen different genera. Among them, 6

Journal Pre-proof 150

the Proteobacteria phylum was predominant and accounted for 96.3%, and the Firmicutes was

151

the second dominant phylum, with a relative abundance of 3.70%. At the genus level,

152

Pseudomonas, Achromobacter, Bacillus, and Azomonas were dominant and accounted for

153

87.2%, 6.12%, 3.51 and 1.80%, respectively (Fig. 1).

154

3.2. Effects of bioremediation on TPH degradation

155

The concentration of total petroleum hydrocarbon (TPH) in the contaminated soil sample

156

was 20,000 mg kg-1 soil (Table 1). After 12 weeks of remediation, TPH was 17230 and 14340

157

mg kg-1 soil in the BA and BS microcosms (Fig. 2), which represented TPH removal

158

efficiencies of 13.9% and 28.3%, respectively. Thus, the BS treatment was more effective than

159

the BA treatment for petroleum degradation in contaminated loessal soil.

160

3.3. Physicochemical properties of petroleum-polluted soil

161

The ORP value of the initial oil-contaminated soil was -48 mV (Fig.3-a). The ORP values

162

changed to -52~-55 mV and -40~-47 mV in the BA and BS treatments, respectively. Soil ORP

163

was enhanced by biostimulation and reduced by bioaugmentation treatment.

164

The initial oil-contaminated soil had pH values of 7.88 (Table 1). During the remediation,

165

the pH values were approximately 7.82 and 7.76 in the BA and BS treatments, respectively

166

(Fig.3-b). The soil pH value was lower in the BS treatment than in the BA treatment during the

167

remediation.

168

Total nitrogen (TN), ammonia nitrogen (NH4+-N), and nitrate nitrogen (NO3--N) were

169

1159.7, 0.38 and 12.39 mg kg-1 in the initial contaminated soil, respectively (Table 1). In the BA

170

treatment, the concentrations of TN, NH4+-N, and NO3--N were 1159.6~1170.4, 4.21~1.38, and

171

12.1~18.5 mg kg-1, respectively (Fig.3-c-d, Fig.3-f). The various forms of nitrogen remained

172

relatively stable during the remediation. In the BS treatment, the contents of TN, NH4+-N, and

173

NO3--N respectively increased to 1650.3, 210.4 and 292.9 mg kg-1 after the amendment of

174

NH4NO3 (zero weeks) for the biostimulation treatment. After 12 weeks of biostimulation 7

Journal Pre-proof 175

treatment, the contents of TN, NH4+-N, and NO3--N decreased to 1523.1, 128.7 and 225.6 mg

176

kg-1, respectively (Fig 3-c-d, Fig.3-f).

177

The content of available phosphorus (AP) in the initial oil-contaminated soil was 15.9 mg

178

kg-1 (Table 1). An increase in soil AP at zero weeks was observed when the bioaugmentation

179

and biostimulation treatments were applied due to the addition of KH2PO4 in the BS treatment

180

or inoculation of PMC flora (which suspended in PBS buffer) in the BA treatment. After one

181

week of incubation, the AP concentration dropped from 175.7 to 93.5 mg kg-1 in the BS sample

182

and thereafter showed a substantial decrease. The AP content decreased from 35.0 to 22.2 mg

183

kg-1 and then increased in the BA treatment (Fig.3-e).

184

3.4. Correlation between the petroleum hydrocarbon and soil physicochemical properties

185

The correlations between the soil physicochemical properties and TPH are given in Table 2.

186

In the BS treatment, the TPH content had a strong positive relationship with TN and AP and a

187

correlation with pH and NH4+-N. In the BA treatment, the TPH content had a positive

188

correlation with NH4+-N and AP. Correlation between the petroleum hydrocarbon degradation

189

and nutrients in the BA treatment are different from BS treatment, which may be attributed to

190

different nitrogen utilization ways in the bioaugmentation and biostimulation strategies.

191

3.5. Evaluation of soil microbial diversity and community composition

192

In this study, we collected bioaugmentation (BA) and biostimulation (BS) soil samples at

193

weeks 1 and 12 for Illumina MiSeq sequencing analysis. They are hereafter referred to as BA1

194

and BS1 for the first week and BA12 and BS12 for the twelfth week of remediation,

195

respectively. Additionally, the initial petroleum contaminated soil (IS) was analyzed by Illumina

196

MiSeq sequencing analysis.

197

3.5.1. Sequencing quality control

198

The initial number of sequences was 45501, 43274, 32449, 42716, and 42269 in the IS,

199

BA1, BS1, BA12 and BS12 samples, respectively. After dislodging low-quality reads using 8

Journal Pre-proof 200

Usearch software, the valid sequence number was 41369, 41031, 30489, 40858, and 40,418,

201

respectively. The sequencing optimization and coverage rates were more than 90% and 98% for

202

all samples, indicating that the optimization procedures were reasonable and the sequencing

203

depths were adequate.

204

The valid sequences were then clustered into OTUs at a threshold of 97%. The results

205

indicated the OTUs number obtained in the IS (1669) and BS1 (1631) samples were the highest,

206

next to the BS12 sample (1217), followed by the BA1 (1027) and the BA12 (765) samples

207

(Table 3). The details of the bacterial taxonomy (phylum, class, order, family, and genus) in

208

different soil samples are shown in Table 3.

209

3.5.2. Soil microbial diversity

210

Table 3 summarize the biodiversity index of different soil samples. The species diversity

211

indices, including Chao 1, Ace, and Shannon indices, produced the highest values for the IS,

212

followed by BS1 and BS12, and the lowest value for the BA12 samples among the treatments

213

(Table 3), which indicated that the bacterial diversity was higher in the IS and BS treatments

214

than in the BA treatment. The BA treatment greatly decreased the biodiversity of the soil

215

microbial community. Compared to the BA treatment, the microbial community composition is

216

stable and diverse in the BS treatment. In addition, the α-diversity decreased along with the

217

remediation time compared with those of the same treatment at different remediation periods

218

(BA1 vs. BA12; BS1 vs. BS12).

219

3.5.3. Microbial community composition

220

The 10 most abundant phyla and top 20 predominant genera in each sample and their

221

relative abundances are shown in Fig. 4. In the initial contaminated soil (IS), Actinobacteria and

222

Proteobacteria were the dominant phyla, with the relative abundances of 47.3% and 37.4%,

223

respectively. The Firmicutes phylum was the subordinate dominant, with a relative abundance

224

of 9.16% (Fig.4-a). At the genus level, Promicromonospora and Exiguobacterium were the

225

dominant genera, with the relative abundances of 19.0% and 8.49%, respectively (Fig. 4-b). 9

Journal Pre-proof 226

Compared to the IS, the microbial community compositions significantly changed in the

227

BA treatment. At the phylum level, the Proteobacteria phylum was the most dominant in the

228

BA1 and BA12 samples, with the relative abundance increased to 96.3% and 87.4%,

229

respectively. Firmicutes (3.71% in the BA1 vs. 2.64% in the BA12) and Actinobacteria (2.46%

230

in the BA1 vs. 5.46% in the BA12) decreased and became the subordinate dominant phyla

231

(Fig.4-a). At the genus level, the Pseudomonas genus was predominant and accounted for 75.2%

232

in the BA1 sample and 76.4% in the BA12 sample. The Achromobacter (8.81% in the BA1 vs.

233

7.56% in the BA12) genus increased and became the secondary dominant genus. There was no

234

discernible change in the BA12 sample compared to the BA1sample (Fig. 4-b).

235

The microbial community composition changed in the BS12 sample compared to that in

236

the BS1 sample. The relative abundance of the Actinobacteria (17.8% in the BS1 vs. 50.4% in

237

the BS12) and Bacteroidetes phyla (0.64% in the BS1 vs. 2.16% in the BS12) increased

238

substantially, while that of the Proteobacteria (35.5% in the BS1 vs. 10.1% in the BS12) and

239

Firmicutes phyla (58.4% in the BS1 vs. 35.5% in the BS12) were significantly reduced with the

240

remediation time. At the genus level, the relative abundances of Nocardioides (1.66% in the

241

BS1 vs. 28.9% in the BS12) and Promicromonospora (6.34% in the BS1 vs. 14.8% in the BS12)

242

increased, while those Bacillus (31.6% in the BS1 vs. 22.7% in the BS12) and Planococcaceae

243

(6.64% in the BS1 vs. 2.07% in the BS12) were substantially reduced during the remediation.

244

4. Discussion

245

Strains isolated from oil contaminated soil that can effectively degrade hydrocarbons often

246

include the Proteobacteria phylum and Pseudomonas genus (Bento et al., 2005; Wang et al.,

247

2005; Han et al., 2009). In this study, the hydrocarbon degrading flora PMC obtained from

248

petroleum contaminated loessal soil primarily consisted of the Proteobacteria phylum and

249

Pseudomonas genus. These results indicate that the classification homology of petroleum

250

degraders isolated from the different soil types and textures. 10

Journal Pre-proof 251

The initial physiochemical parameters of the contaminated soil are valuable to the success

252

of remediation. It has been widely accepted that the optimum C/N/P ratio in soil for

253

biostimulation strategy is 100/10/1, which is equivalent to the nutrients required for active

254

bacteria (US EPA, 2002; Wu et al., 2016). Additionally, some studies suggested that the

255

appropriate degrading population in soil for bioaugmentation strategy is approximately 108

256

cfug-1of soil (Abaloset al., 2004; Masy et al., 2016).

257

The physiochemical characteristics of petroleum-contaminated loessal soil suggested that

258

the C/N/P ratio in the polluted soil was 660/10/0.14. Additionally, the initial contaminated soil

259

had lower populations of TPH degraders (Table 1). Theoretically, the addition of a suitable

260

amount of nitrogen and phosphorus for biostimulation or inoculation with the appropriate

261

hydrocarbon degraders for bioaugmentation are both useful strategies for the bioremediation of

262

the petroleum-contaminated soil.

263

Studies have demonstrated that the introduction of microbial flora enriched from

264

contaminated soil and inoculated back to the same soil was more effective in degrading

265

hydrocarbons than flora obtained from other soils (Wu et al., 2013). In this study, we isolated

266

the indigenous petroleum hydrocarbon degrading flora PMC from contaminated soil that was

267

actively under bioremediation and added it back to the same polluted soil. The results of high

268

throughput sequencing revealed the main compositions of flora PMC belonged to the

269

Proteobacteria phylum (96.3%) and Pseudomonas genus (87.2%) (Fig. 1). After one week of

270

inoculation, the relative abundance of the Proteobacteria phylum increased from 37.4% to

271

87.4% (Fig. 4-a), and that of the Pseudomonas genus increased from 2.99% to 76.4% in the BA

272

treatment (Fig. 4-b). After 12 weeks of incubation, the Proteobacteria phylum and

273

Pseudomonas genus were still the most dominant microorganisms in the BA12 sample. These

274

results suggest that the inoculated oil-degrading flora PMC survived and grew quickly in

275

contaminated soil. However, the TPH degradation efficiency was lower in the bioaugmentation

276

treatment (13.9%) than that in the biostimulation treatment (28.3%). It seems that the TPH 11

Journal Pre-proof 277

degrader population is not a prerequisite for petroleum hydrocarbon degradation.

278

The contents of TN, NH4+-N NO3--N, and AP are high in the BS treatment because of the

279

addition of NH4NO3 and KH2PO4 for the biostimulation treatment. The NH4+-N, NO3--N, and

280

AP contents decreased rapidly after one week of remediation, which can be explained by

281

microbial utilization and petroleum hydrocarbon degradation in soil under biostimulation (John

282

et al., 2011; Urakawa et al., 2012; Marchand et al., 2017). In addition, the soil pH value was

283

lower and ORP was higher in the biostimulation treatment than those in the bioaugmentation

284

treatment. The strong positive correlations of the nutrient contents with the petroleum

285

concentration in the BA and BS treatments indicated the importance of nutrients for the

286

bioremediation processes.

287

The similarity of the microbial community compositions in BA1 and BA12 suggested that

288

soil microbial community structure achieved balance rapidly in the bioaugmentation treatment

289

after one week of remediation. According to Fig. 4, soil microbial community structure in the

290

first week of bioaugmentation (BA1) greatly changed compared with the initial soil (IS), but

291

then remained basically unchanged in the twelfth week (BA12). Also, TPH degradation is the

292

highest in the first week during the whole remediation (shown Fig. 2). The highest TPH

293

degradation efficiencies in the first week indicated that the microbial community structure shift

294

may be associated with the petroleum hydrocarbon degradation.

295

Proteobacteria, Firmicutes and Actinobacteria were the dominant phyla across all soil

296

samples. These phyla are widely distributed in petroleum-contaminated loessal soil (Liu et al.,

297

2013). The relative abundance of these strains changed in different bioremediation strategies.

298

In the bioaugmentation treatment, the relative abundance of the Proteobacteria phylum and

299

Pseudomonas genus significantly increased, which can be attributed to the inoculation of the

300

petroleum-degrading flora PMC in contaminated soil. The inoculants become the most dominant

301

microorganisms in the contaminated soil, which restrained indigenous microorganism growth

302

during the remediation. 12

Journal Pre-proof 303

In the BS treatment, the relative abundance of the Proteobacteria phylum decreased, while

304

that of the Firmicutes phylum increased and became dominant compared to the IS. At the genus

305

level, the relative abundances of Exiguobacterium and Promicromonosporawere reduced, while

306

those of Nocardioides and Bacillus enhanced and became dominant.

307

According to the literature, the Firmicutes phylum and Bacillus genera are primary

308

petroleum degraders in oil-contaminated soil (Yeh, 1985; Ortega-González et al., 2015; Yang et

309

al., 2015). The increase of the Firmicutes phylum and Bacillus genus indicate the correlativity

310

between these strains and TPH degradation in the biostimulation treatment.

311

Biodiversity is dependent on both the richness and evenness of the microbial composition.

312

The diversity indices of Chao 1 and Ace represent the richness, while Shannon and Simpson

313

indices represent the evenness of a microbial community (Sengupta and Dick, 2015).

314

The bioaugmentation treatment led to the lowest microbial biodiversity among all soil

315

samples (Table 3). Soil microbial species were relatively single and inhomogeneous in the

316

bioaugmentation treatment, which constrained their degradation capability. On the contrary, the

317

microbial community composition was more suitable and diverse in the biostimulation

318

treatment. Thus, the advantage of biostimulation was ensuring that the soil microbial community

319

was well-matched to the specific soil environment, which led to improved petroleum

320

hydrocarbon degradation.

321

Previous studies have revealed that the success of oil-contaminated bioremediation depends

322

on the presence of sufficient microorganisms with the appropriate metabolic capabilities (Krutz

323

et al., 2005; Wu et al., 2017b). Our study found that TPH degrader enrichment was not a

324

prerequisite for enhancing hydrocarbon degradation and that the balance and diversity of the

325

microbial community structure were essential for petroleum hydrocarbon remediation. We also

326

demonstrated that the evenness and richness of species were vital for contaminant remediation

327

in the polluted soil. 13

Journal Pre-proof 328

5. Conclusion

329

Biostimulation with nitrogen and phosphorous represents a more effective remediation

330

strategy for TPH degradation than bioaugmentation with petroleum degrading flora in

331

contaminated loessal soil. In the bioaugmentation treatment, the inoculum was able to survive,

332

flourish quickly and become the most dominant microorganisms. Microbial species

333

was inhomogeneous by the bioaugmentation treatment. The soil microbial community

334

changed considerably and reached stability after one week of remediation. In the biostimulation

335

treatment, the microbial community composition was more diverse than that in the

336

bioaugmentation treatment. The balance and diversity of the microbial community structure

337

were beneficial to petroleum hydrocarbon decontamination.

338

Acknowledgments

339

This work was supported by the National Natural Science Foundation of China (No.

340

21577109), the Program for Innovative Research Team in Shaanxi (PIRT) (Grant No.

341

2013KCT-13), the Natural Science Foundation of Shaanxi Province (2015JM5163), and the Key

342

Laboratory Project of the Shaanxi Provincial Education Department (13JS048).

343

References

344

Abalos, A., Vinas, M., Sabate, J., Manresa, M.A., Solanas, A.M., 2004. Enhanced

345

biodegradation of Casablanca crude oil by a microbial consortium in presence of a

346

rhamnolipid produced by Pseudomonas aeruginosa AT10. Biodegradation 15, 249-260.

347

Abed, R.M.M., Al-Sabahi, J., Al-Maqrashi, F., Al-Habsi, A., Al-Hinai, M., 2014. Characteri

348

-zation of hydrocarbon-degrading bacteria isolated from oil-contaminated sediments in the

349

Sultanate of Oman and evaluation of bioaugmentation and biostimulation approaches in

350

microcosm experiments. Int. Biodeterior. Biodegrad. 89, 58-66.

351

Bento, F. M., Camargo, F. A. D. O., Okeke, B. C., Jr, W. T. F., 2005. Diversity of biosurfactant

352

producing microorganisms isolated from soils contaminated with diesel oil. Microbiol. Res. 14

Journal Pre-proof 353

160(3), 249-255.

354

Bao, Y. J., Xu, Z., Yang, L., Zhi, Y., Sun, J., Hui, S., 2017. High-throughput metagenomic

355

analysis of petroleum-contaminated soil microbiome reveals the versatility in xenobiotic

356

aromatics metabolism. J. Environ. Sci. 56 (6), 25-35.

357

Bundy, J. G., Paton, G. I., Campbell, C. D., 2004. Combined microbial community level and

358

single species biosensor responses to monitor recovery of oil polluted soil. Soil Biol.

359

Biochem. 36(7), 1149-1159.

360

Castiglione, M. R., Giorgetti, L., Becarelli, S., Siracusa, G., Lorenzi, R., Gregorio, S. D., 2016.

361

Polycyclic aromatic hydrocarbon-contaminated soils: bioaugmentationof autochthonous

362

bacteria and toxicological assessment of the bioremediation process by means of vicia faba,

363

l. Environ. Sci. Pollut. Res. 23(8), 1-12.

364

Han, P., Zheng, L., Cui, Z. S., Guo, X. C., & Tian, L., 2009. Isolation, identification and

365

diversity analysis of petroleum-degrading bacteria in shengli oil field wetland soil. Chinese

366

J. Appl. Ecol. 20(5), 1202-1208.

367 368 369

Jia, J., Li, G., Zhang, X., Lu, X., Dai, D., Zhong, Y., 2005. Microbial communities in petroleum polluted soils by pcr-dgge. Journal of Tsinghua University 45(9), 1217-1220. Jia, J., Zong, S., Hu, L., Shi, S., Zhai, X., Wang, B., 2017. The dynamic change ofmicrobial

370

communities in crude oil-contaminated soils from oil fields in china. J.Soil Contam. 26(2),

371

171-183.

372 373

John, R. C., Itah, A. Y., Essien, J. P., Ikpe, D. I., 2011. Fate of nitrogen-fixing bacteriain crude oil contaminated wetland ultisol. B. Environ. Contam.Tox. 87(3), 343-353.

374

Kauppi, S., Sinkkonen, A., Romantschuk, M., 2011. Enhancing bioremediation diesel fuel

375

contaminated soil in a boreal climate: comparison of biostimulation and bioaugmentation.

376

Int. Biodeterior. Biodegrad. 65, 359-368.

377

Krutz, L.J., Beyrouty, C.A., Gentry, T.J., 2005. Selective enrichment of a pyrene degrader

378

population and enhanced pyrene degradation in Bermuda grass rhizosphere. Biol. Fert. Soils. 15

Journal Pre-proof 379

41, 359-364.

380

Leys, N. M., Bastiaens, L., Verstraete, W., & Springael, D., 2005. Influence of the

381

carbon/nitrogen/phosphorus ratio on polycyclic aromatic hydrocarbon degradation by

382

mycobacterium, and sphingomonasin soil. Appl. Microbiol. Biotech. 66(6), 726-736.

383

Liu, P. W. G., Chang, T. C., Whang, L. M., Kao, C. H., Pan, P. T.,Cheng, S. S., 2011.

384

Bioremediation of petroleum hydrocarbon contaminated soil: effects of strategies and

385

microbial community shift. Int. Biodeterior. Biodegrad. 65(8), 1119-1127.

386 387

Liu, P. W. G., Chang, T. C., Chen, C. H., Wang, M. Z., Hsu, H. W., 2013. Effects of soil organic

388

matter and bacterial community shift on bioremediation of diesel-contaminated soil.

389

Int. Biodeterior. Biodegrad. 85(11),661-670.

390

Marchand, C., St-Arnaud, M., Hogland, W., Bell, T.H., Hijri, M., 2017. Petroleum

391

biodegradation capacity of bacteria and fungi isolated from petroleum-contaminated soil.

392

Int. Biodeterior. Biodegrad. 116, 48–57.

393

Margesin, R., Hämmerle, M., Tscherko, D., 2007. Microbial activity and community

394

composition during bioremediation of diesel-oil-contaminated soil:effects of hydrocarbon

395

concentration, fertilizers, and incubation time. Microb. Ecol. 53(2), 259-269.

396

Masy, T., Demanèche, S., Tromme, O., Thonart, P., Jacques, P., Hiligsmann, S., 2016.

397

Hydrocarbon biostimulation and bioaugmentation in organic carbon and clay-rich

398

soils. Soil Biol. Biochem. 99, 66-74.

399

Megharaj, M., Ramakrishnan, B., Venkateswarlu, K., Sethunathan, N., Naidu, R., 2011.

400

Bioremediation approaches for organic pollutants: a critical perspective. Environ.

401

Int. 37(8), 1362-1375.

402 403 404

Mishra, S., Jyot, J., Kuhad, R. C., Lal, B., 2001. In situ bioremediation potential ofan oily sludge-degrading bacterial consortium. Curr. Microbiol. 43(5), 328-335. National Agricultural TechnologyExtension Service Center. 2014. Technicalspecifications for 16

Journal Pre-proof 405 406

soil analysis[M]. 2th Edition. Beijing: China Agricultural Press. Ortega-González, D. K., Martínez-González, G., Flores, C. M., Zaragoza, D., Cancino-Diaz, J.

407

C., Cruz-Maya, J. A., 2015. Amycolatopsis, sp. poz14 isolated from oil-contaminated soil

408

degrades polycyclic aromatic hydrocarbons. Int. Biodeterior. Biodegrad. 99, 165-173.

409

Pacwa-Płociniczak, M., Płaza, G. A., Piotrowska-Seget, Z., 2016. Monitoring the changes in a

410

bacterial community in petroleum-polluted soil bioaugmented with hydrocarbon-degrading

411

strains. Appl. Soil Ecol. 105, 76-85.

412

Peng, M., Zi, X., Wang, Q., 2015. Bacterial community diversity of oil-contaminated soils

413

assessed by high throughput sequencing of 16s rRNA genes. Int.J.Env. Res. Pub. He. 12(10),

414

12002.

415

Petrov, A. M., Versioning, A. A., Karimullin, L. K., Akaikin, D. V., Tarasov, O. Y., 2016.

416

Dynamics of ecological and biological characteristics of soddy-podzolic soils under

417

long-term oil pollution. Eurasian Soil Sci. 49(7), 784-791.

418

Qu, Y.Y., Zhang, X.W., Shen, W.L., Ma, Q., You, S.N., Pei, X.F., Li, S.Z., Ma, F., Zhou, J.T.,

419

2016. Illumina MiSeq sequencing reveals long-term impacts of single-walled carbon

420

nanotubes on microbial communities of wastewater treatment systems. Bioresour. Technol.

421

211, 209-215.

422 423 424

Ramadass, K., Megharaj, M., Venkateswarlu, K., Naidu, R., 2015. Ecologicalimplications of motor oil pollution: earthworm survival and soil health. Soil Biol. Biochem. 85, 72-81. Sayara, T., Borras, E., Caminal, G., Sarr a, M., S anchez, A., 2011. Bioremediation of PAHs

425

contaminated soil through composting: Influence of bioaugmentation and biostimulation on

426

contaminant biodegradation. Int. Biodeterior. Biodegrad. 65, 859-865.

427

Sengupta, A., Dick, W. A., 2015. Bacterial community diversity in soil under two tillage

428

practices as determined by pyrosequencing. Microb. Ecol.70(3), 853-859.

429 430

Shahi, A., Aydin, S., Ince, B., Ince, O., 2016. Evaluation of microbial population and functional genes during the bioremediation of petroleum-contaminated soil as an effective monitoring 17

Journal Pre-proof 431

approach. Ecotox. Environ. Safe. 125,153-160.

432

Tyagi, M., Da, F. M., de Carvalho, C. C., 2011. Bioaugmentation and biostimulation strategies

433

to improve the effectiveness of bioremediation processes. Biodegradation 22(2), 231-41.

434

Urakawa, H., Garcia, J.C., Barreto, P.D., Molina, G.A., Barreto, J.C., 2012. A sensitive crude

435

oil bioassay indicates that oil spills potentially induce a change of major nitrifying

436

prokaryotes from the Archaea to the Bacteria. Environ. Pollut. 164, 42-45.

437

U.S. Environmental Protection Agency, 2002. Application, performance, and costs of

438

biotreatment technologies for contaminated soils, EPA/600/R-03/037, Battelle Contract No.

439

68-C-00-185.

440

Wang, Z. Q., Wu, Q., Ye, S. Y., Li, F. Q., Xie, H. L., Jin, X. L., 2005. Study on the ground

441

water

442

petroleum contaminant biodegradation by high efficient microorganism. Environ. Sci. 26(6),

443

61-64.

444

Wu, M. L., Chen, L. M., Tian, Y. Q., Ding, Y., Dick, W.A., 2013. Degradation of polycyclic

445

aromatic hydrocarbons by microbial consortia enriched from three soils using two different

446

culture media. Environ. Pollut. 178, 152-158.

447

Wu, M. L., Dick, W. A., Li, W., Wang, X. C., Yang, Q., Wang, T. T., 2016. Bioaugmentation

448

and biostimulation of hydrocarbon degradation and the microbial community in a

449

petroleum-contaminated soil. Int. Biodeterior. Biodegrad. 107(3), 158-164.

450

Wu, M. L., Ye, X. Q., Chen, K. L., Li, W.,Yuan, J.,Jiang, X., 2017a. Bacterial community shift

451

and hydrocarbon transformation during bioremediation of short-term petroleum

452

contaminated soil. Environ. Pollut. 223, 657-664.

453

Wu, M. L, W., Dick, W. A., Ye, X.Q., Chen, K. L., Kost, D., Chen L. M., 2017b.

454

Bioremediation of hydrocarbon degradation in a petroleum-contaminated soil and

455

microbial population and activity determination. Chemosphere169, 124-130.

456

Xu, Y. H., Lu, M., 2010. Bioremediation of crude oil-contaminated soil: comparison of different 18

Journal Pre-proof 457

biostimulation and bioaugmentation treatments. J. Hazard. Mater. 183(1), 395-401.

458

Yan, L., Sinkko, H., Penttinen, P., Lindström, K., 2016. Characterization of successional

459

changes in bacterial community composition during bioremediationof used motor

460

oil-contaminated soil in a boreal climate. Sci.Total Environ. 542(Pt A), 817-825.

461

Yang, Q., Wu, M. L., Nie, M. Q., Wang, T. T., Zhang, M. H., 2015. Effects and biological

462

response on bioremediation of petroleum contaminated soil. Environ. Sci.36(5), 1856-1863.

463

Yeh, G. C., 1985. Treatment and separation of petroleum and related materials.US, US4515685.

464

Yulia, M., Polyak, Ludmila, G., Bakina, Marina, V., Chugunova, M. V., 2017. Effectof remediation

465

strategies on biological activity of oil-contaminated soil-a field study. Int. Biodeterior.

466

Biodegrad. 126, 57-68.

467

Zhang, H.H., Feng, J., Chen, S.N., Zhao, Z. F., Li, B. Q., Wang, Y., Jia, J. Y., Li, S. L., Wang,

468

Y., Yan, M. M., Lu, K. Y., Hao, H. Y., 2018. Geographical Patterns of nirS Gene

469

Abundance and nirS-Type Denitrifying Bacterial Community Associated with Activated

470

Sludge

471

10.1007/s00248-018-1236-7.

from

Different

Wastewater

Treatment

472 473 474

19

Plants.

Microb.

Ecol.

DOI:

Journal Pre-proof Fig.1. hydrocarbon degrading flora isolated from petroleum contaminated soil. Fig.2. Degradation of TPH by biostimulation and bioaugmentation in oil-contaminated soil. Errors bars indicate ±SD of triplicate samples. Different letters in the same week represent a significant difference at P<0.05. Fig.3. Changes of physicochemical properties of soil. (3-a ORP; phosphorus;

3-b pH;

3-c NH4+-N;

3-d NO3-N;

3-e Available

3-f Total nitrogen).

Fig.4. Microbial community composition in the BA and BS treatments from Illumina sequencing (4-a The most abundant bacterial phyla; 4-b The most abundant bacterial genus).

Journal Pre-proof

Fig.1

Journal Pre-proof

BA

19000

a

a

BS

a a

-1

TPH concentrations(mg kg )

20000 a

18000

a

a

17000

b b

16000

a

a b b

b

b

15000

b

14000

b

b

13000 12000 0

1

2

3

4

5

6

7

Time(weeks) Fig.2

8

9

10 11 12 13

Journal Pre-proof

BA -30

8.6

(a)

-35 -40

8.2

-45

8.0

-50 -55

7.4

-65

7.2

-70 0

2

3

4

5

6

7

8

+

7.0

9 10 11 12

-1

200

1

2

3

4

5

6

7

8

9 10 11 12

2

3

4

5

6

7

8

9 10 11 12

2

3

4

5

6

7

8

9 10 11 12

(d)

300 240

160 120

-

180

80

120

40

60

0 0 180

1

2

3

4

5

6

7

8

0

9 10 11 12

0

1800 -1

TN (mg kg )

(e)

160 -1

0

NO3 -N (mg kg )

(c)

-1

NH4 -N (mg kg )

1

360

240

AP (mg kg )

7.8 7.6

-60

140 120 100 80

1

(f)

1700 1600 1500 1400

60 40

1300

20

1100

0

(b)

8.4

pH

ORP (mv)

BS

1200

0

1

2

3

4

5

6

7

8

9 10 11 12

1000

Time(weeks)

0

1

Time(weeks)

Fig.3

Journal Pre-proof

100 Others Verrucomicrobia

Relative abundance(%)

80

Planctomycetes Gemmatimonadetes

60

Chloroflexi Candidatus Saccharibacteria

40

Bacteroidetes Acidobacteria

20

Proteobacteria Actinobacteria

0

Firmicutes

IS

BA1

BA12

BS1

BS12

100

(b)

Others Olivibacter Achromobacter Planococcaceae_incertae_sedis Fictibacillus Paenisporosarcina Paenibacillus Arthrobacter Blastococcus Altererythrobacter Pseudoxanthomonas Acinetobacter Leifsonia Mycobacterium Nocardioides Unclassified Citrobacter Exiguobacterium Bacillus Pseudomonas Promicromonospora

80

Relative abundance(%)

(a)

60 40 20 0

IS

BA1

BA12

Treatments

Fig.4

BS1

BS12

Journal Pre-proof 

Biostimulation is better than bioaugmentation for petroleum hydrocarbon degradation.



Bioaugmentation decrease soil microbial diversity.



Soil microbial are more diverse in the biostimulation than in the bioaugmentation.



The advantage of biostimulation is ensuring soil microbial evenness and diversity.

Journal Pre-proof Table 1 Physicochemical and microbiological properties of the petroleum-contaminated soil (mean ± standard deviation (SD), n = 3). Main characteristics

Values

TPHa(g kg-1)

19.8±0.38

Total carbon(%)

7.65 ±0.05

Organic matter(%)

0.67± 0.03

Moisture content(%)

5.40±0.11

pH

7.88± 0.11

Total nitrogen(mg kg-1)

1159.7± 95.2

Available phosphorus(mg kg-1)

171.6 ± 18.4

Ammonia nitrogen(mg kg-1)

0.38 ± 0.08

Nitrate nitrogen(mg kg-1)

12.39± 0.87

Total bacterial numbers(cells g-1)

(6.8 ± 0.3)×106

TPH degraders(MPNb g-1)

(3.30 ± 0.6)×104

aTPH:

total petroleum hydrocarbon most probably number

bMPN:

Journal Pre-proof Table 2 Correlation matrix (coefficients and significance levels) for parameters of BS and BA treatments. pH ORP N-NH4+ N-NO3TN AP Biostimulation treatment TPH 0.784* -0.303 pH 1 -0.682* ORP 1 + N-NH4 N-NO3TN AP Bioaugmentation treatment TPH -0.185 -0.29 pH 1 0.034 ORP 1 + N-NH4 N-NO3TN AP

0.786* 0.567* 0.197 1

0.036 0.311 -0.529* -0.287 1

0.806** 0.746* -0.111 0.831** 0.294 1

0.908** 0.922** -0.364 0.825** 0.157 0.915** 1

0.574* -0.745* -0.213 1

-0.329 0.709* -0.5* -0.171 1

-0.467 -0.561* -0.351 0.964** 0.098 1

0.508* -0.499 -0.787* 0.456 0.061 0.477 1

*,p> 0.05;**,p> 0.01.

TPH: total petroleum hydrocarbon; ORP: oxidation-reduction potential; nitrogen; AP: Available phosphorus.

TN: total

Journal Pre-proof

Table 3 Bacterial taxonomy and diversity indices of the BA and BS treatments at the 1st and 12th weeks of remediation. Treatments

IS

BA1

BS1

BA12

BS12

1669

1027

1631

765

1217

Phylum

21

17

22

20

19

Class

42

37

49

43

46

Order

67

51

70

60

64

Family

136

100

139

115

132

Genus

353

241

380

248

322

4.38

2.00

4.10

1.62

3.25

Ace index

3327.1

2674.4

3180.4

2699.6

3009.2

Chao1 index

2534.1

1815.4

2481.2

1584.2

2088.9

0.05

0.39

0.07

0.5

0.12

sequencing analysis OTUs num

Diversity indices Shannon index

Simpson index