Total Phosphate Influences the Rate of Hydrocarbon Degradation but

Apr 15, 2016 - The dominant mineralogical constraints on community composition were the relative amounts of adsorbed phosphate, brushite, and newberyi...
1 downloads 8 Views 2MB Size
Subscriber access provided by University of Pennsylvania Libraries

Article

Total phosphate influences the rate of hydrocarbon degradation but phosphate mineralogy shapes microbial community composition in cold-region calcareous soils. Steven D Siciliano, Tingting Chen, Courtney Leigh Phillips, Jordan Graeme Hamilton, David M Hilger, Blaine Chartrand, Jay Grosskleg, Kris Bradshaw, Trevor Carlson, and Derek Peak Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b05911 • Publication Date (Web): 15 Apr 2016 Downloaded from http://pubs.acs.org on April 15, 2016

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 33

Environmental Science & Technology

1

Total phosphate influences the rate of hydrocarbon degradation but phosphate mineralogy shapes

2

microbial community composition in cold-region calcareous soils.

3 4

Steven D Siciliano1, Tingting Chen1, Courtney Phillips1, Jordan Hamilton1, David Hilger1,

5

Blaine Chartrand2, Jay Grosskleg3, Kris Bradshaw3, Trevor Carlson3 and Derek Peak1*.

6 7

1

Department of Soil Science, University of Saskatchewan, Saskatoon, Canada 2

8 3

9

Saskatchewan Polytechnique, Saskatoon, Canada

Federated Cooperatives Limited, Saskatoon, Canada

10

Title Running Head

11

Phosphate mineralogy controls microbial hydrocarbon degradation *Corresponding Author

12

Contact Information

13

Derek Peak

14

51 Campus Drive, Department of Soil Science, Agriculture Bldg

15

University of Saskatchewan

16

Saskatoon, SK, Canada, S7N 5A8

17

Phone: 306-966-6806. Fax: 306-966-6881

18

Email: [email protected]

19

Key words: contaminated soil, bioremediation, mineralogy, hydrocarbons, phosphate,

20 21 22 23

XANES, catabolic genes Limit of 7000: Current = 4598 words: 2 small figures = 600, 3 large figures = 1500: One table = 300 Total = 6998 words

1 ACS Paragon Plus Environment

Environmental Science & Technology

24

Abstract

25

Managing phosphorus bioaccessibility is critical for the bioremediation of hydrocarbons in

26

calcareous soils. This paper explores how soil mineralogy interacts with a novel biostimulatory

27

solution to both control phosphorus bioavailability and influence bioremediation. Two large bore

28

infiltrators (1 m diameter) were installed at a PHC contaminated site and continuously supplied

29

with a solution containing nutrients and an electron acceptor. Soils from 8 contaminated sites

30

were prepared and pretreated, analyzed pre-trial, spiked with diesel, placed into nylon bags into

31

the infiltrators, and removed after 3 months. From XAS, we learned that 3 principal phosphate

32

phases had formed: adsorbed phosphate, brushite, and newberyite. All measures of

33

biodegradation in the samples (in-situ degradation estimates, mineralization assays, culturable

34

bacteria, catabolic genes) varied depending upon the soil’s phosphate speciation. Notably,

35

adsorbed phosphate increased anaerobic phenanthrene degradation and bzdN catabolic gene

36

prevalence. The dominant mineralogical constraints on community composition were the relative

37

amounts of adsorbed phosphate, brushite, and newberyite. Overall, this study finds that total

38

phosphate influences microbial community phenotypes whereas relative percentages of

39

phosphate minerals influences microbial community genotype composition.

40 41

2 ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33

Environmental Science & Technology

42

Introduction

43

Retail fuel distribution releases diesel and gasoline into soil. In North America’s calcareous

44

glacial till soils, natural attenuation of hydrocarbons can be limited because these soils are near

45

freezing for half the year and soil water moves through preferential flow paths1-4. One

46

remediation approach for these soils is to stimulate soil organisms to degrade hydrocarbons5 by

47

adding nitrogen, phosphorus, and electron acceptors (e.g., nitrate, ferric iron or ferric sulfate).

48

This approach works ex-situ6, 7 but often fails in-situ, perhaps due to limited microorganisms,

49

electron acceptors, or nutrients8, 9. In particular, bioavailability of phosphorus can limit

50

subsurface remediation efforts10.

51

Phosphate in solution is the most bio-accessible form of P. However, phosphate has a

52

strong tendency to adsorb to Fe, Al, and Ca minerals resulting in low solution concentrations 11,

53

12

54

Determining phosphate speciation in soils can be done using synchrotron-based X-ray absorption

55

near edge spectroscopy (XANES) to probe the average local bonding environment of phosphate

56

13-20

57

components using linear combination fitting and a reference library of phosphate standard

58

compounds21, 22, although this approach has difficulties distinguishing components with similar

59

spectral signatures17 and may need improved baseline correction algorithms23 to be quantitative.

60

. Therefore, identifying phosphate speciation is an important step in the remediation process.

. It is possible to semi-quantitatively deconvolute a heterogeneous sample into its

Carbonate mineralogy also plays an important role in the chemical fate of phosphate in

61

calcareous soils. The potential for adsorption and subsequent precipitation of calcium phosphates

62

onto carbonate minerals can limit the accessible fraction of soil phosphate16. However,

63

identifying minerals using laboratory X-ray diffraction (XRD) often requires sample pre-

64

treatments to remove organic matter, carbonates, and sand/silt coarse soil fractions24-26.

3 ACS Paragon Plus Environment

Environmental Science & Technology

65

Synchrotron powder XRD eliminates the need for sample pre-treatments27, 28, and enables clay

66

mineralogy to be determined without affecting soil carbonates.

67

Investigators usually use a range of approaches to demonstrate biodegradation, including

68

characterizing community composition29; conducting catabolic gene assays8, viability assays30,

69

and mineralization assays31; and applying environmental fate models of hydrocarbon

70

biodegradation32. Likewise, we used a range of approaches to evaluate how phosphorus

71

mineralogy is linked to hydrocarbon degradation. To characterize 16S rRNA composition, we

72

used the EMIRGE approach, which applies universal primers to amplify almost the entire 16S

73

rRNA gene and then fragments, sequences, and assembles these amplicons in silico33. EMIRGE

74

avoids some primer biases associated with amplifying a small target area and is best used to

75

assess dominant members of the microbial community34.

76

Measuring catabolic gene levels helps determine an increase in organisms containing the

77

catabolic pathway of interest. However, because primers do not amplify all organisms containing

78

the catabolic gene, PCR-based approaches are susceptible to false negatives. Using multiple gene

79

targets reduces this false negative risk. A key step in anaerobic degradation of alkanes is adding

80

fumarate to an alkane to form 1-methylalkyl-succinate, which is catalyzed by 1-methylalkyl-

81

succinate synthase35 and assessed by QPCR primers targeting the assA gene36. We assessed

82

anaerobic aromatic degradation using the following genes: (i) bssA, which encodes the α-subunit

83

of benzylsuccinate synthase that detects degraders using fumarate to activate the aromatic

84

substrate37 and form benzylsuccinate, which is further activated to produce benzylsuccinyl-CoA;

85

(ii) bzd or bcr, which encode class I benzoyl-CoA reductases to catalyze reduction of benzoyl-

86

CoA38, 39 by Azoarcus or Thauera strains, respectively40; and (iii) bamA, found in facultative and

87

obligate anaerobes41, which catalyzes the hydrolysis of 1-hydroxy-benzoyl-CoA.

4 ACS Paragon Plus Environment

Page 4 of 33

Page 5 of 33

88

Environmental Science & Technology

Gene presence is not necessarily linked to activity; therefore, assays demonstrating

89

mineralization and increases in culturable degraders are needed to confirm biodegradation30, 31.

90

Assays using 13C hydrocarbons demonstrate mineralization by detecting enriched 13CO2 in the

91

microcosm42 headspace. Most probable number (MPN) estimates of culturable organisms avoid

92

biases associated with DNA extraction and amplification and thus can confirm catabolic gene

93

prevalence estimates. Environmental fate models can estimate in-situ degradation rates32,

94

confirming mineralization assays and MPN estimates.

95

We hypothesized that biostimulatory solutions worked at some sites but not others due to

96

variation in phosphate speciation, whereas site-specific effects not linked to phosphate are

97

secondary. Soil samples were derived from the treated and untreated soils from 8 hydrocarbon

98

contaminated sites plus 2 common end-members of phosphorus mineralogy, hydroxyapatite and

99

phosphate absorbed on hematite. Treatment consisted of chemical extraction to remove Ca (and

100

add Mg) to force a greater range of phosphate solubility. These 18 samples were then duplicated

101

and incubated in duplicate large bore infiltrators at an active field site undergoing N and

102

phosphorus biostimulation with a solution we designed to stimulate bioremediation. At 0, 3, and

103

6 months, we removed the soils and assessed phosphate speciation (via laboratory methods and

104

XANES), solid-state mineralogy (via synchrotron powder XRD), in-situ degradation activity,

105

hydrocarbon mineralization potential, culturable degraders, catabolic gene prevalence, and

106

bacterial community composition. Our intent was to evaluate our hypothesis by using molecular

107

tools to link soil mineralogy, nutrient bioaccessibility, and hydrocarbon remediation.

108

Materials and Methods

109

Contaminated Sites

5 ACS Paragon Plus Environment

Environmental Science & Technology

110

Using a push core drill, soil from 8 (Colonsay, Young, 33rd Sand Lens, 33rd Street Clay

111

Lens, Broadway South Site, Broadway West Side, Meadow Lake, Winnipeg Street) retail

112

gasoline and diesel outlets in Alberta and Saskatchewan were collected from 2-4 m below

113

ground between May and August, assessed for hydrocarbon contamination, and stored at -20°C.

114

Prior to collection, all soils had been exposed to on-site gasoline and diesel leaks for a least a

115

year and had been frozen for 4-6 months per year.

116

Large Bore Injection System

117

In late 2013, two 5.5 m deep infiltrators (1 m diameter) were installed by a rotary bucket

118

caisson auger at an urban site and backfilled with coarse-grained silica sand. Each infiltrator had

119

1 infiltration well (5 cm), 1 overflow well connected to the infiltration well (30 cm), 3 wells (10

120

cm) in which we hung in situ reactors, and 1 standard monitoring well (5 cm). Submersible

121

pumps in the underground storage tank were connected to infiltrator lines, with floats and auto

122

shutoff switches installed in the tanks and infiltrators. Amendment delivery began August 14,

123

2014. The tank solution was pH 6.5 and contained 11 mM MgSO4, 1 mM H3PO4, and 0.08 mM

124

HNO3 in municipal potable water. These concentrations were selected to: (i) match local

125

groundwater chemistry for sulfate, (ii) be at the maximum allowed groundwater nitrogen based

126

on local guidelines, and (iii) be below solubility for any phosphate mineral formation in both the

127

storage tanks and the groundwater of the site based upon PHREEQC modeling of both tap and

128

groundwater. This experimental design provided a constant nutrient solution level to the in-situ

129

reactors over a relatively long time scale to satisfy our research needs, while also respecting the

130

industrial partner’s remediation requirements. It was pumped into the infiltrators 0.5 m below

131

ground at a flow rate of ~5000 L/day.

132

In-Situ Reactors

6 ACS Paragon Plus Environment

Page 6 of 33

Page 7 of 33

Environmental Science & Technology

133

Samples were air-dried, sieved to 850 µm, and divided into 2 treatment groups (in

134

duplicate): control and Ca-depleted. The control soil was not treated. The Ca-depleted soil was

135

prepared by adding deionized water to a 1:1 soil slurry, adjusting soil pH to 6.5 with 1.0 M

136

H2SO4, centrifuging, and washing with deionized water 3 times to remove excess ions as well as

137

CaSO4(s). The depleted, washed soil was added to 1.0 M MgCl2 in 1:4 ratio, stirred with a

138

magnetic bar overnight, centrifuged, and washed with deionized water 3 times. This entire

139

procedure was repeated 3 times to exchange bound Ca2+ for Mg2+. Hematite (#310050, Sigma-

140

Aldrich) and hydroxyapatite (CAS No.12167-74-7, Alfa Aesar) were used for the mineral control

141

groups.

142

Soils (15 g) were placed in an in-situ reactor (12.5 cm × 4 cm 60-micron mesh Nitex

143

membrane [Aqua Dynamic] bag) and spiked with retail diesel (100 µl) because the process of

144

altering mineralogy removed most volatile hydrocarbons. Reactors were zip-tied onto Tygon

145

tube (Fisher Scientific) and suspended 3 m underground (2.5 m below the water level) in the

146

research wells. The in-situ reactors were placed into the infiltrators on August 22, 2014. On

147

August 22 and December 9, 2014, and February 23, 2015, the reactors were removed and soil

148

was extracted from the in-situ reactors in an anaerobic chamber then immediately tested for

149

DNA, hydrocarbons, phosphate, and MPNs.

150

Chemical and Mineralogical Speciation

151

Synchrotron powder diffraction measurements were collected at the Canadian Light

152

Source (CLS) CMCF-BM beamline (08B1-1) using a wide-area detector. The data was

153

processed using GSAS-II43, with phase identification performed using X'pert Highscore Plus

154

software v. 3.0 (PANalytical, Inc.).

7 ACS Paragon Plus Environment

Environmental Science & Technology

155

Phosphorus speciation was performed in laboratory with a sequential chemical extraction

156

scheme modified from the Hedley extraction14. Extract analyses were performed via colorimetry

157

using Folio Instruments AA1 autoanalyzer. All P K-edge XANES were collected at the CLS

158

SXRMB beamline (08-B1-1) in fluorescence mode, using a 4-element solid-state detector

159

(Bruker). Multiple scans were collected to ensure adequate signal-to-noise for analysis; all data

160

was processed and analyzed using the Demeter software package44. A library of phosphate

161

standard compounds (pure and mixed Ca and Mg phosphates), organic phosphate compounds,

162

iron and aluminum phosphate minerals, and adsorption complexes of phosphate on mineral

163

surfaces were analyzed; standards were diluted to 1%wt phosphorus in boron nitride powder if

164

necessary.

165

DNA Extraction and Microbial Community Analysis

166

Soil DNA was extracted with the FastDNATM SPIN Kit for Soil (MP Biomedicals), and

167

DNA concentration and purity determined by Qubit® Fluorometer (Life Technologies). A 16S

168

rDNA PCR was performed with primers 27F and 1492R across a gradient of 8 temperatures (48-

169

54°C)33. Generation and sequencing of amplicons was as detailed by Handley et al. (2014),

170

except a Nextera DNA Sample Prep Kit was used to shear the DNA, and DNA was sequenced on

171

an Illumina Miseq. The supplemental information details library preparation and bioinformatics.

172

Catabolic Gene Prevalence

173

Quantitative PCR (qPCR) was performed on the 7500 Real-Time PCR System (Applied

174

Biosystems). The qPCR mixtures contained 10 µl SYBR Green master mix (Qiagen), 1 µl of 10

175

μM primer pairs, 6 µl water, and 2 µl template DNA, with a final volume of 20 µl. Primers were

176

ordered from Integrated DNA Technologies (Supplemental Table 1). For alkB, a conventional

177

qPCR was performed as follows: 95°C × 10’ for 1 cycle; 95°C × 1’, 50°C × 1’, and 72°C× 1’ for

8 ACS Paragon Plus Environment

Page 8 of 33

Page 9 of 33

Environmental Science & Technology

178

35 cycles; followed by a melt curve. For nah, rpoB, assA, bzdN, bcrC, and bamA, touchdown

179

qPCR reactions were performed as follows: 95°C × 10’ for 1 cycle; 95°C × 1’, 65°C × 1’

180

(decreasing 1°C per cycle), and 72°C × 1’ for 11 cycles; 95°C × 1’, 47-65°C (47°C for nah, 49°C

181

for rpoB, 56°C for assA and bamA, 58°C for bzdN, 65°C for bcrC) × 1’, and 72°C× 1’ for 35

182

cycles; followed by a melt curve. Standard strain for alkB was Pseudomonas putida ATCC 2934,

183

and for nah, Pseudomonas putida ATCC 17484. For other genes, clone libraries were created

184

with E. coli cells using the TOPO TA Cloning Kit (Invitrogen) from environmental DNA. DNA

185

was then extracted from clones and standardized to 109 copies.

186

Aerobic and Anaerobic Mineralization Assays

187

For aerobic assays, soil (2.5 g) from the in-situ bioreactor was transferred to 2 dram vials

188

and stored at 4°C for