Chemotaxis Increases the Retention of Bacteria in Porous Media with

Publication Date (Web): June 1, 2018 ... partitioned into the aqueous phase to create localized chemoattractant gradients throughout the column. A pul...
0 downloads 0 Views 930KB Size
Subscriber access provided by University of Winnipeg Library

Environmental Processes

Chemotaxis increases the retention of bacteria in porous media with residual NAPL entrapment Joanna S.T. Adadevoh, C. Andrew Ramsburg, and Roseanne Marie Ford Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01172 • Publication Date (Web): 01 Jun 2018 Downloaded from http://pubs.acs.org on June 1, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

Environmental Science & Technology

1

Chemotaxis increases the retention of bacteria in

2

porous media with residual NAPL entrapment

3

Joanna S. T. Adadevoh,† C. Andrew Ramsburg,‡ and Roseanne M. Ford*,†

4

†Department of Chemical Engineering, University of Virginia, Charlottesville, Virginia 22904

5

‡Department of Civil and Environmental Engineering, Tufts University, Medford, Massachusetts

6

02155

7

Keywords: chemotaxis, bioremediation, porous media, NAPL, transport phenomena,

8

Pseudomonas putida

9 10

ABSTRACT

11

Chemotaxis has the potential to decrease the persistence of non-aqueous phase liquid (NAPL)

12

contaminants in aquifers by allowing pollutant-degrading bacteria to move toward sources of

13

contamination and thus influence dissolution. This experimental study investigated the migratory

14

response of chemotactic bacteria to a distribution of residual NAPL ganglia entrapped within a

15

laboratory-scale sand column under continuous-flow at a superficial velocity of 0.05 cm/min.

16

Naphthalene dissolved in a model NAPL 2,2,4,4,6,8,8-heptamethylnonane partitioned into the

17

aqueous phase to create localized chemoattractant gradients throughout the column. A pulse

18

mixture of equal concentrations of Pseudomonas putida G7, a strain chemotactic to naphthalene,

19

and Pseudomonas putida G7 Y1, a non-chemotactic mutant, was introduced to the column and

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 23

20

effluent bacterial concentrations were measured with time. Breakthrough curves (BTCs) for the

21

two strains were noticeably different upon visual inspection. Differences in BTCs (compared to

22

non-chemotactic controls) were quantified in terms of percent recovery and were statistically

23

significant (p< 0.01). Chemotaxis reduced percent recovery in the effluent by 45% thereby

24

increasing the population of bacteria that were retained within the column in the vicinity of

25

residual NAPL contaminants. An increase in flow rate to a superficial velocity of 0.25 cm/min

26

did not diminish cell retention associated with the chemotactic effect.

27 28

INTRODUCTION

29

Non-aqueous phase liquids (NAPLs) typically comprise persistent organic pollutants that pose

30

long-term threats to groundwater quality.1–5 In fact, complete dissolution of entrapped NAPL

31

contaminants may take many decades.2,6,7 Chemotaxis-aided bioremediation has been recently

32

highlighted as a promising technology for more rapid aquifer restoration.8–10 Chemotaxis is a

33

phenomenon in which motile cells, such as bacteria, can detect the concentration gradient of a

34

chemical and move preferentially towards a region of higher chemical/chemoattractant

35

concentration if it is advantageous to do so.11–13 Chemotaxis has the potential to improve

36

bioremediation strategies by augmenting the mass transfer of pollutant-degrading bacteria to the

37

source of NAPL contamination.

38

Bacterial chemotaxis in porous media is well studied. In microfluidic devices, and laboratory-

39

scale sand columns and microcosms, researchers documented enhanced migration of chemotactic

40

bacteria towards attractant sources in directions transverse to flow in both physically

41

homogeneous and heterogeneous porous media.14–16 Controlled field studies showed that the

42

coupled effects of cell chemotaxis and proliferation greatly enhanced bacterial transport towards

ACS Paragon Plus Environment

2

Page 3 of 23

Environmental Science & Technology

43

a chemical source.17 In all these studies, bacteria exhibited chemotaxis towards macroscopic

44

contaminant plumes. NAPLs entrapped in porous media typically form discrete ganglia, each

45

serving as a source for contaminant dissolution.2–4 In essence, what results is a distribution of

46

multiple, localized, concentration gradients within a NAPL source zone. A few studies have

47

investigated chemotactic bacterial response to residual NAPL entrapment at the micro-scale.5,18

48

However, this phenomenon is yet to be explored in the meso- and macro-scales. This work with

49

NAPLs builds on a recent study in which localized concentration gradients in a sand column

50

were created by randomly distributing solid naphthalene crystals within the packed bed.19 While

51

similar in design, solid contaminants present different physical interfaces and chemical

52

characteristics than are present in porous media containing entrapped NAPL.

53

Hence, the purpose of this study was to investigate the migratory response of chemotactic

54

bacteria to a distribution of residual NAPL ganglia entrapped within a lab-scale sand-packed

55

column. This work is important for assessing the role of bacterial chemotaxis on the

56

bioremediation and biorestoration of NAPL source zones.

57

MATERIALS AND METHODS

58

Bacteria and Culture Conditions. Two motile bacterial strains were used in this study:

59

Pseudomonas putida G7 (PpG7), which is chemotactic to naphthalene, and Pseudomonas putida

60

G7 Y1 (PpG7 Y1), a nahY::Km non-chemotactic mutant strain.20 A detailed description of the

61

bacterial culture conditions is provided in a previous article by Adadevoh and co-workers.19

62

Before proceeding with transport experiments, cell motility was verified under oil immersion at

63

100× with a Zeiss microscope (F100/1.25 oil). According to a previous study no difference in

64

bacteria swimming speed was found between the two strains. 19

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 23

65

Column Assembly. Sand column experiments were performed in duplicates for each type of

66

column set-up (i.e., NAPL experiments, and No NAPL experiments); for each replicate, the sand

67

column was re-assembled as follows. A glass chromatography column (diameter 4.8 cm, length

68

15.5 cm) was dry packed with 30-40 mesh quartz sand (VWR item BDH9274) in 1-2 cm

69

increments under mixing and vibration. The average gravimetrically estimated porosity for

70

replicate sand columns was 0.38. To aid subsequent buffer saturation, air was displaced from the

71

sand-packed column by flushing the column with several pore volumes of CO2 as CO2 is more

72

soluble in the aqueous phase than air. The sand column was then saturated with 10% (v/v)

73

random motility buffer (RMB; 11.2 g/L K2HPO4 (Fisher, 7758-11-4), 4.8 g/L KH2PO4

74

(Amresco, 7778-77-0), 0.029 g/L EDTA (Sigma-Aldrich, 60-00-4)) at a superficial flow rate of 5

75

mL/min, against gravity, using an ÄKTAexplorer (Amersham Pharmacia Biotech, 18-1300-00).

76

After fully saturating the sand column with 10% RMB, a conservative, non-reactive, solute tracer

77

test (10 mL pulse of 0.05 M NaNO3 (Sigma-Aldrich, 221341) at a superficial rate of 0.905

78

mL/min) was conducted to quantify the dispersivity within the column. To capture tracer

79

breakthrough, 5 mL effluent samples were collected continuously via a fraction collector

80

(Pharmacia Biotech, Frac-900) for 1.6 pore volumes (i.e., 190 minutes). Nitrate concentrations in

81

the effluent samples were quantified via absorbance at 300 nm using a Beckman Coulter

82

spectrophotometer (DU® 640).

83

2,2,4,4,6,8,8-Heptamethylnonane (HMN; Acros Organics, 4390-04-9) was selected as a model

84

NAPL based on its prior use when examining the influence of chemotaxis on desorption and

85

degradation of naphthalene.8,22 HMN is not known to elicit a chemotactic response in PpG7 and

86

was not toxic to the bacteria at the concentrations used in this study. The HMN NAPL contained

87

33 g/L naphthalene (Fisher, 91-20-3) for chemotaxis experiments and no naphthalene for HMN

ACS Paragon Plus Environment

4

Page 5 of 23

Environmental Science & Technology

88

control experiments. The aqueous-phase concentration of naphthalene in equilibrium with the

89

HMN NAPL containing 33 g/L naphthalene was 15.4 mg/L (see the Supporting Information).

90

For experiments containing NAPL, a residual NAPL saturation was produced following the

91

tracer test using well established protocols.21 In brief, NAPL was introduced into the column at a

92

constant flow rate of 1 mL/min using the ÄKTAexplorer. Upon reaching the maximum HMN

93

NAPL saturation of 75%, aqueous flow (10% RMB buffer) was re-established and continued

94

until reaching a residual, entrapped NAPL saturation of 19%. The NAPL saturation was

95

determined gravimetrically. After establishing the residual NAPL saturation, a second tracer test

96

was conducted, in a manner similar to that of the first tracer test, to determine the resulting

97

aqueous-phase dispersivity.3 All tracer tests were conducted at a superficial flow rate of 0.905

98

mL/min.

99

Bacterial Transport Experiment and Sample Analysis. Bacterial transport experiments

100

were performed in duplicate and for each run, the column packing and NAPL entrapment

101

process was repeated, as previously described. A 10 mL mixture (~10% of the pore volume) of

102

equal concentrations of PpG7 and PpG7 Y1 (~2 × 108 cells/mL in 10% RMB) was introduced

103

into the sand column containing entrapped globules of naphthalene dissolved in HMN at a

104

superficial velocity of 0.905 mL/min against gravity. Prior to bacterial transport experiments,

105

PpG7 and PpG7 Y1 were stained with red (FM4-64; Molecular Probes, T3166) and green

106

(calcein AM; Molecular Probes, C1430) fluorescent dyes, respectively, to aid cell differentiation.

107

FM4-64 and calcein AM have excitation/emission wavelengths of ~515/640 nm and ~495/516

108

nm, respectively. A previous study showed that these stains had no effect on the cell size, swim

109

speed, or zeta potential of PpG7 and PpG7 Y1.19 The same study also showed that the transport

110

of either of these strains is not influenced by the presence of the other. Bacterial effluent samples

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 23

111

(i.e., 5 mL) were collected continuously from the column via a fraction collector (Pharmacia

112

Biotech, Frac-900) for 1.7 pore volumes (i.e., 200 min). Once 1.7 pore volumes was reached, the

113

flow rate was increased by a factor of 5, while 5 mL effluent samples were still being collected,

114

to observe the effect of increased flow rate on bacterial transport within the column. Bacterial

115

concentrations in effluent samples were enumerated via a flow cytometer (BD Accuri™ C6).

116

The flow cytometry procedure utilized in this study is the same as that previously provided by

117

Adadevoh and co-workers.19 Concentrations as low as 105 cells/mL were detectable using this

118

method. After each transport experiment, cell motility, and hence viability, was verified via

119

visual inspection under oil immersion at 100× with a Zeiss microscope (F100/1.25 oil).

120

Naphthalene concentration in the column effluent was also measured through the course of the

121

experiment once residual NAPL entrapment was achieved. Samples collected during the

122

bacterial experiments were first filtered using 0.22 µm PTFE syringe filters (Celltreat Scientific

123

Products, 229757) to remove bacterial cells. Naphthalene was quantified via absorbance at 220

124

nm using a Shimadzu Prominence UFLC equipped with a UV detector. Isocratic (85%

125

acetonitrile and 15% DI H2O) separation was accomplished on a C-18 column at a flow rate of

126

0.6 mL/min.22 To serve as control experiments, bacterial transport was also observed in columns

127

containing only sand (i.e., no NAPL) and in columns containing entrapped HMN with no

128

dissolved naphthalene. Interstitial velocities for columns containing entrapped NAPL and for

129

columns containing no NAPL were 2.33 ± 0.04 m/d and 1.88 ± 0.02 m/d, respectively (note that

130

± represents 1 standard deviation about a mean for replicate experiments, unless noted

131

otherwise). For columns containing no NAPL, bacterial experiments were conducted for 2 pore

132

volumes (i.e., 240 min) before an increase in flow rate. Cell proliferation on naphthalene over the

133

course of these experiments was negligible.19,23

ACS Paragon Plus Environment

6

Page 7 of 23

Environmental Science & Technology

134

Quantitative Assessment of Breakthrough Curves. Previous studies have used a one-

135

dimensional advection-dispersion equation to quantify differences in experimentally acquired

136

bacterial breakthrough curves (BTCs).15,19,24,25 One such study employed fitted parameters of the

137

1D equation as apparent values used to empirically assess differences in transport between

138

chemotactic and non-chemotactic bacteria.19 We have adopted a similar approach in this study.

139

The following equation was used to quantify differences in shape and magnitude of BTCs that

140

correspond to certain features (e.g. spread in the peak, position of the peak, and area under the

141

peak) for chemotactic and non-chemotactic bacterial populations      =  −  −   (1)   

142

where b is the dimensionless species concentration in the aqueous phase [-], t is time [T], z is the

143

longitudinal position [L], R is used to quantify the reversible retention of bacterial mass [-], Dbz

144

is the longitudinal hydrodynamic dispersion coefficient [L2T-1], vf is the interstitial pore water

145

velocity [LT-1], and km represents irreversible, first-order retention of bacteria within the medium

146

[T-1]. The nonlinear least-squares parameter optimization method in CXTFIT was employed to

147

fit the 1D equation to the experimental BTC data.26 Specifically, R, Dbz, and km parameter values

148

were fit. Longitudinal dispersivity, αz [L], was calculated via  =  +   where Deff

149

[L2T-1] is the effective motility coefficient for bacteria in porous media.9 For P. putida strains in

150

30-40 mesh quartz sand, Deff was reported to be 1.3 × 10-5 cm2/s.11,15,17,27

151

Note that the parameters in Equation 1 do not explicitly represent chemotaxis – the 1D

152

advection-dispersion equation was employed only to show that there are quantitative differences

153

in BTCs for chemotactic bacteria by comparing apparent or effective parameters obtained by

154

fitting Equation 1. Thus, the reader is cautioned against attributing these parameter values to the

155

description of any particular physical process. Also note that Equation 1 does not include kinetic

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 23

156

sorption/desorption terms, so it cannot capture the tailing phenomenon that is evident in the

157

experimental BTCs. Our goal in this study was not to develop a mathematical model for bacterial

158

transport, rather we used Equation 1 to quantify differences between various characteristics of

159

BTCs that were evident upon visual inspection.

160

Cell normalized mean travel time (τ) and percent recovery from the sand column were also

161

calculated via the first and zeroth moment of the bacterial BTCs, respectively.28 Additionally, the

162

1D advection-dispersion equation was fit to tracer BTCs in a manner similar to that described for

163

bacterial BTCs with the exception that only Dbz was fit while R and km were prescribed at 1.0 and

164

0, respectively.

165

RESULTS AND DISCUSSION

166

Sand Column Characterization. Three types of sand column experiments were conducted in

167

this study: (1) water-saturated sand columns with no HMN NAPL and no naphthalene (i.e., No

168

NAPL); (2) sand columns with entrapped HMN NAPL, but no naphthalene (i.e., HMN NAPL);

169

and (3) sand columns with entrapped HMN NAPL containing 33 g/L naphthalene (i.e., HMN-

170

NAP NAPL). Table 1 lists characteristic parameters for two replicates of each type of sand

171

column experiment. Prior to NAPL entrapment, the average pore water velocity and standard

172

deviation within the sand column was determined to be 1.88 ± 0.02 m/d for six experimental runs

173

(i.e., duplicate experiments for HMN-NAP NAPL, HMN NAPL, and No NAPL). After

174

entrapment, residual NAPL saturation within the column was gravimetrically estimated to be 19

175

± 2% and the resulting interstitial velocity was 2.33 ± 0.04 m/d. For HMN-NAP NAPL

176

experiments, effluent naphthalene concentrations were observed to remain at a steady-state value

177

of 10.2 ± 0.1 mg/L. Longitudinal dispersivity assessed by the nitrate tracer tests was 0.39 ± 0.12

ACS Paragon Plus Environment

8

Page 9 of 23

Environmental Science & Technology

178

mm. Values for the dispersivity parameter assessed from bacterial BTCs are reported in Table 2

179

with an average of 1.8 ± 0.75 mm.

180

Influence of NAPL on Bacterial Transport. Chemotactic (PpG7) and non-chemotactic

181

(PpG7 Y1) bacterial transport was observed in duplicate experiments in each of the previously

182

described types of sand columns: (1) water-saturated sand columns with no HMN NAPL and no

183

naphthalene; (2) sand columns with entrapped HMN NAPL, but no naphthalene; and (3) sand

184

columns with entrapped HMN NAPL containing 33 g/L naphthalene. BTCs in Figure 1 c-d were

185

from columns that contained HMN, but no naphthalene. Because there was no chemoattractant,

186

the BTCs were expected to look the same for both chemotactic and non-chemotactic species and

187

this was the case. The open circles in panels c and d were consistent with what we expected to

188

observe. Bacteria did not exhibit a chemotactic response in the absence of naphthalene. The open

189

diamond symbols also showed no difference between chemotactic and nonchemotactic species.

190

Although we see differences in the qualitative aspects of the BTCs (comparing diamonds to

191

circles), the differences between the percent recoveries – one of the main quantitative descriptors

192

– is similar to that observed between the two BTCs in Figure 1b. There is the possibility that the

193

column packing yielded differences in the qualitative trend of the BTCs, such as the peak values.

194

This, however, should not detract from the main message that we obtain from the BTCs –

195

bacteria did not exhibit a chemotactic response in the absence of naphthalene.

196

To focus on the effect of residual NAPL entrapment on bacterial transport, we direct our

197

attention to bacterial migration in columns containing only HMN and in those with no NAPL

198

(Figure 1c-f). The presence of NAPL within the column led to an earlier breakthrough of cells

199

and resulted in an average of 17% decrease in normalized mean travel time, τ, (Table 2),

200

compared to results obtained in the absence of NAPL for both PpG7 and PpG7 Y1. The earlier

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 23

201

cell breakthrough is a result of the greater interstitial velocity arising in those experiments

202

occurring in the presence of NAPL (i.e., 19% lower water saturation with a similar superficial

203

velocity). For experiments with no NAPL, bacterial normalized mean travel time was greater

204

than 1 due to reversible retention of cells in the sandy medium; values of τ are comparable to

205

those previously reported for these bacterial strains in 30-40 mesh quartz sand (e.g., τ = 1.26).19

206

Bacterial cell recovery was found to be approximately 50% lower in the presence of entrapped

207

NAPL (38-46% in the absence of NAPL versus 18-25% in the presence of the HMN NAPL).

208

Correspondingly, the parameter associated with irreversible retention (km) was found to be

209

approximately double, while reversible retention (R) remained relatively similar across all

210

experiments (Table 2). The decrease in cell recovery was expected as a previous study reported

211

that due to bacterial motility alone, cells may encounter the NAPL-aqueous interface more

212

frequently and subsequently be retained at/near the interface.22 In that study Law and Aitken22

213

observed greater accumulation of motile bacteria on the surface of an HMN drop by direct visual

214

inspection with microscopy. Another study documented that bacterial cells had a stronger

215

affinity for an oil-water interface compared to a glass-water interface and suggested that the

216

increased affinity was due to attractive energies and capillary forces.29

217

Differences in colloidal attachment and retention at immiscible liquid interfaces and solid-fluid

218

interfaces have been explained using interaction energy profiles.30 Thus, energy profiles of PpG7

219

interaction with HMN NAPL, and PpG7 interaction with sand were calculated based on the

220

Derjaguin-Landau-Verwey-Overbeek (DLVO) theory (see Figure S1 in the Supporting

221

Information). Repulsive energy barriers were compared given that this energy has to be

222

overcome for irreversible attachment of cells to occur.31 The repulsive energy barrier between

223

the cells and HMN NAPL was found to be 715 kT less than the repulsive energy barrier between

ACS Paragon Plus Environment

10

Page 11 of 23

Environmental Science & Technology

224

the cells and the sand. The reduced repulsive energy barrier for bacterial cell interaction with

225

NAPL implies that the cells have a greater tendency to be irreversibly retained near the NAPL

226

surface as compared to the sand surface. A greater attachment of cells to the NAPL surface

227

results in the decreased cell percent recovery observed in this study for NAPL experiments, even

228

when chemotaxis is not occurring.

229

Influence of Chemotaxis on Bacterial Transport in Sand Columns with Residual NAPL.

230

Chemotaxis-directed bacterial transport resulted in a distinct decrease in the recovery of PpG7

231

(10-13%, Figure 1a) as compared to both the NAPL control (18-25%, Figure 1c) conducted with

232

PpG7 and the chemotaxis control conducted with PpG7 Y1 (18-24%, Figure 1b). A one-tailed T-

233

test revealed a statistically significant decrease (p