Autochthonous Bioaugmentation-Modified ... - ACS Publications

Jan 29, 2018 - and carcinogenicity.3 Phenanthrene (PHE) is employed as a model PAH compound due to its ..... a mass spectrometric detector (model 5975...
1 downloads 16 Views 504KB Size
Subscriber access provided by READING UNIV

Article

Autochthonous bioaugmentation modified bacterial diversity of phenanthrene degraders in PAH-contaminated wastewater as revealed by DNA-stable isotope probing Jibing Li, Chunling Luo, Dayi Zhang, mengke song, Xixi Cai, Longfei Jiang, and Gan Zhang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b05646 • Publication Date (Web): 29 Jan 2018 Downloaded from http://pubs.acs.org on January 31, 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 36

Environmental Science & Technology

1

Autochthonous bioaugmentation modified bacterial diversity of phenanthrene

2

degraders in PAH-contaminated wastewater as revealed by DNA-stable isotope

3

probing

4 5

Jibing Li,†, Ѱ Chunling Luo,†, ‡* Dayi Zhang,§ Mengke Song, ‡ Xixi Cai,# Longfei Jiang,† Gan

6

Zhang†

7 8



9



Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China College of Natural Resources and Environment, South China Agricultural University, Guangzhou,

10

510642, China

11

§

12

#

13

350002, China

14

School of Environment, Tsinghua University, Beijing 100084, China College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou,

Ѱ

University of Chinese Academy of Sciences, Beijing, 100049, China

15 16

*Corresponding author: Dr. Chunling Luo

17

E-mail: [email protected]

18

Tel.: +86-20-85290290; Fax: +86-20-85290706

1

ACS Paragon Plus Environment

Environmental Science & Technology

19 20

TOC

21 22

2

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

Environmental Science & Technology

23

ABSTRACT

24

To reveal the mechanisms of autochthonous bioaugmentation (ABA) in wastewater contaminated

25

with polycyclic aromatic hydrocarbons (PAHs), DNA-stable-isotope-probing (SIP) was used in

26

the present study with the addition of an autochthonous microorganism Acinetobacter tandoii LJ-5.

27

We found LJ-5 inoculum produced a significant increase in phenanthrene (PHE) mineralization,

28

but LJ-5 surprisingly did not participate in indigenous PHE degradation from the SIP results. The

29

improvement of PHE biodegradation was not explained by the engagement of LJ-5 but attributed

30

to the remarkably altered diversity of PHE degraders. Of the major PHE degraders present in

31

ambient wastewater (Rhodoplanes sp., Mycobacterium sp., Xanthomonadaceae sp. and

32

Enterobacteriaceae sp.), only Mycobacterium sp. and Enterobacteriaceae sp. remained functional

33

in the presence of strain LJ-5, but five new taxa Bacillus, Paenibacillus, Ammoniphilus,

34

Sporosarcina and Hyphomicrobium were favoured. Rhodoplanes, Ammoniphilus, Sporosarcina

35

and Hyphomicrobium were directly linked to, for the first time, indigenous PHE biodegradation.

36

Sequences of functional PAH-RHDα genes from heavy fractions further proved the change in PHE

37

degraders by identifying distinct PAH-ring hydroxylating dioxygenases between ambient

38

degradation and ABA. Our findings indicate a new mechanism of ABA, provide new insights into

39

the diversity of PHE-degrading communities, and suggest ABA as a promising in situ

40

bioremediation strategy for PAH-contaminated wastewater.

3

ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 36

41

1. INTRODUCTION

42

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental pollutants consisting of

43

two or more fused aromatic rings.1 These compounds have the potential to pose serious health

44

risks to all organisms.2 They are of great environmental concern due to their high toxicity,

45

mutagenicity and carcinogenicity.3 Phenanthrene (PHE) is employed as a model PAH compound

46

due to its ubiquity in nature and fused-ring angular structure.4 Bioremediation is considered as a

47

cost-effective and eco-friendly approach for elimination of PAHs from natural environments.5

48

PAHs in a polluted environment can be partially degraded by the autochthonous microbial

49

population via natural attenuation. When indigenous microbes lack the capacity to degrade PAHs,

50

the introduction of PAH degraders isolated from other PAH-contaminated sites,6,

51

bioaugmentation (BA), is a promising technology to encourage PAH degradation.8, 9 However, the

52

efficacy of BA is variable, since the survival and the PAH-degrading ability of introduced

53

microorganisms are highly dependent on environmental conditions.10,

54

difficulties, autochthonous bioaugmentation (ABA) has been proposed as an alternative method.12

55

ABA is defined as a BA technology that uses indigenous degraders in the polluted sites (soil,

56

sediment and water) with the aim of accelerating and enhancing the biodegradation potential.13

57

11

7

known as

To overcome these

Many microorganisms capable of degrading PAHs have been isolated using cultivation-based

58

approaches and have been evaluated for their potential application in ABA strategies.7,

59

Cultivation-based methods are used to determine the phenotypes and metabolic characteristics of

60

PAH degraders and the functional genes associated with PAH metabolism such as PAH-ring

61

hydroxylating dioxygenase.16 Particularly, the gene encoding the alpha subunit of the PAH-ring

62

hydroxylating dioxygenase (PAH-RHDα) is usually used as the biomarker to quantify the bacterial 4

ACS Paragon Plus Environment

14, 15

Page 5 of 36

Environmental Science & Technology

63

populations capable of degrading PAH.16 However, these approaches share an critical drawback in

64

that the majority of microbes are uncultivable,17 and many autochthonous PAH-degrading

65

microorganisms remain resistant to traditional cultivation approaches.4,

66

cultivation greatly underestimates the microbial diversity in ecosystems and fails to explain the

67

complex interactions between microorganisms and the natural environment.19

18

In addition, direct

68

In recent years, cultivation-independent methods have been used to evaluate the prokaryotic

69

diversity of complex systems and estimate the PAH degradation potential of indigenous

70

consortia.20-22 High-throughput sequencing has revolutionised our ability to investigate the

71

microbial communities in environmental samples by providing higher resolution of microbial taxa

72

compared with conventional cloning techniques.23 Nevertheless, the metabolic features of

73

organisms are difficult to infer accurately using this approach.23 Stable-isotope probing (SIP) is a

74

cultivation-independent technique that circumvents the requirement of isolating an organism to

75

assess metabolic responses and link its identity to function.24 This technique relies on the

76

incorporation of a stable isotope-labelled substrate and the identification of active microorganisms

77

by analysing isotope-enriched cellular components, such as DNA, RNA and protein.4, 25 We can

78

therefore characterise the target microorganisms resistant to cultivation.18 To date, SIP has been

79

used to identify a large number of indigenous bacteria capable of degrading PAHs.20, 26-28

80

In our previous work, we isolated a novel PHE degrading strain, Acinetobacter tandoii LJ-5,

81

from PAH-contaminated wastewater (37°49′N, 118°25′E; altitude, 37.49 m) collected in Shandong

82

Province, China in September 2015.18 A. tandoii LJ-5 is a major autochthonous PHE degrader due

83

to its PHE metabolism confirmed by its enrichment in the heavy DNA fraction in the presence of

84

13

C-PHE.18 A. tandoii LJ-5 is a Gram negative, rod-shaped, obligate aerobe lacking flagella. This 5

ACS Paragon Plus Environment

Environmental Science & Technology

85

strain grew well in minimal medium supplemented with 100–1,000 mg/L PHE under optimal

86

conditions (pH 7.0 and 30°C), and a satisfactory PHE degradation efficiency (> 80%) was

87

achieved within 3 days for low PHE concentrations (100 and 200 mg/L).18

88

Although the ABA strategy can significantly promote PHE mineralisation by introducing

89

autochthonous strains, altering the overall microbial activity and population,15, 29 no studies have

90

confirmed the degrading capacity of the reintroduced degraders, and little is known about changes

91

in the diversity of indigenous PHE degraders during ABA. In the present study, wastewater

92

microcosms were established to evaluate the bioremediation potential of ABA in

93

PAH-contaminated wastewater with A. tandoii LJ-5 inoculum. Additionally, DNA-SIP was used to

94

investigate the changes in the diversity of indigenous PHE degraders by comparing treatments in

95

the presence versus absence of LJ-5, and to link the indigenous bacterial taxa with their PHE

96

biodegradation phenotypes. The PHE-degrading bacteria in different wastewater treatments were

97

successfully characterised using DNA-SIP and high-throughput sequencing. Furthermore, the

98

functional PAH-RHD gene was investigated by analysing relevant sequences amplified from the

99

heavy fractions of

13

C-PHE treatments. We have provided useful information regarding the

100

mechanisms of ABA in aiding bioremediation of PAH-contaminated wastewater and other

101

potentially polluted environments.

6

ACS Paragon Plus Environment

Page 6 of 36

Page 7 of 36

Environmental Science & Technology

102

2. MATERIALS AND METHODS

103

2.1. Sample collection. Water samples (~10 L) were collected in March 2016 from industrial

104

wastewater at the same oil refinery (37°49′N, 118°25′E; altitude, 37.49 m) in Shandong Province

105

(China) as described in our previous work.18 After transport to the laboratory, the wastewater was

106

stored immediately at 4°C. Portions of the samples were stored at −80°C for subsequent DNA

107

extraction, and the remainder was promptly used in PHE degradation and SIP experiments. PAHs

108

identified in the wastewater are listed in the Supporting Information Table S1 (determined using

109

gas chromatography–mass spectrometry as described below).

110

2.2. Cultivation of bacteria. A. tandoii LJ-5 was previously isolated from wastewater in

111

September 2016 at the same oil refinery, using PHE as the carbon source,18 and deposited in the

112

Marine Culture Collection of China (MCCC; 1K02142). A. tandoii LJ-5 is a promising candidate

113

for remediation of PAH-contaminated water, as it has been shown to be an indigenous PHE

114

degrader by both DNA-SIP and direct cultivation. A. tandoii LJ-5 was stored in glycerol solution

115

(15% v/v) at −80°C, thawed and re-cultivated as described previously.18 Briefly, it was cultivated

116

in minimal medium supplemented with 1,000 mg/L PHE in the dark for 18–24 h with shaking at

117

180 rpm at 30°C. Enriched cultures were transferred into a sterile centrifuge tube. Pelleted cells

118

were obtained after centrifugation for 10 min (3000 × g; Eppendorf, Centrifuge 5810R, Hamburg,

119

Germany) and resuspended in phosphate-buffered saline (PBS). Finally, the population of A.

120

tandoii LJ-5 was determined to be approximately 2 × 108 colony forming units/mL (CFU/mL),

121

using the dilution plate counting method.30 Details of the surface tension measurement was

122

provided in the Supporting Information.

123

2.3. SIP microcosms. Microcosms were established in 150 mL serum bottles with 50 mL 7

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 36

124

wastewater with and without inoculation of strain LJ-5. For treatment with strain LJ-5, 100 µL of

125

the cell/PBS suspension were added to the wastewater, thereby establishing an initial inoculum

126

population of approximately 4 × 105 CFU/mL at the beginning of each treatment.10 The unlabelled

127

PHE (99%) or

128

Tewksbury, MA, USA) was added to the above bottles at a final concentration of 10 mg/L. The

129

biotic treatments, which allowed comparisons of ambient degradation (AD) and ABA, are referred

130

to as

131

(13C-PHE without strain LJ-5), and

132

PHE (non-PHE control) and treatments with unlabelled PHE in filter-sterilised wastewater (sterile

133

control) were also established. Each treatment was performed in triplicate. The microcosms were

134

incubated using the method described in our previous work.18 The residual PHE on days 3 and 6

135

after incubation was extracted and detected. Almost all of the PHE was degraded on day 6 (Figure

136

S3), indicating that the majority of PHE was mineralised and incorporated by the microbes during

137

the first 3 days. Thus, samples from each treatment were collected on day 3 for PHE analysis and

138

DNA extraction to avoid cross-feeding.

13

C-labeled PHE (13C14-PHE, 99%, Cambridge Isotope Laboratories, Inc.,

12

C_AD (12C-PHE without strain LJ-5), 13

12

C_ABA (12C-PHE with strain LJ-5),

13

C_AD

C_ABA (13C-PHE with strain LJ-5). Microcosms without

139

2.4. Nucleic acid extraction and ultracentrifugation. DNA extraction from triplicate

140

microcosms (25 mL) of each treatment was performed and quantified as described previously.18

141

Ultracentrifugation of DNA was conducted as described in our previous report.31 Briefly,

142

approximately 3 µg DNA from each microcosm were added to Quick-Seal polyallomer tubes (13

143

× 51 mm, 5.1 mL, Beckman Coulter, Pasadena, CA, USA) and mixed with Tris-EDTA (pH

144

8.0)/caesium chloride solution at a final buoyant density (BD) of ~1.77 g/mL. After balancing and

145

heat-sealing, density gradient centrifugation was performed in an ultracentrifuge (Optima 8

ACS Paragon Plus Environment

Page 9 of 36

Environmental Science & Technology

146

L-100XP, Beckman Coulter, USA) at 45,000 rpm (20°C) for 48 h. Centrifuged gradients were

147

fractionated into different fractions of 400 µL each. The DNA fractions were purified after the BD

148

of each fraction was measured. The relationships between BD and the fraction number or DNA

149

concentration for all four biotic treatments are listed in Figures S1 and S2, respectively.

150

2.5. Quantitative polymerase chain reaction (qPCR). Abundance of the bacterial 16S

151

ribosomal RNA (rRNA) gene in each fraction was determined by qPCR using the universal

152

bacterial primer pair Bac519F/Bac907R (Table 1). The 20 µL PCR mixture contained 10 µL

153

SYBR green PCR Premix Ex Taq II (Takara Bio Inc., Japan), 0.2 µL each primer (10 µM), and 1

154

µL DNA template. PCR was performed using the ABI 7500 real-time PCR system (Applied

155

Biosciences, USA). Tenfold serial dilutions of known copy numbers of the plasmid DNA extracted

156

from Escherichia coli were generated to produce a standard curve. The reactions were conducted

157

as follows: initial denaturation at 94°C for 10 min, followed by 40 cycles at 94°C for 30 s, 55°C

158

for 30 s, and 72°C for 15 s. The SYBR green signal was measured after a 20 s step at 72°C in each

159

cycle. The melt curve was obtained by heating from 60 to 95°C at a rate of 0.2°C/s.

160

2.6. High-throughput sequencing and computational analyses. The hypervariable V4

161

region of bacterial 16S rRNA gene fragments from the DNA fractions derived from the 12C_AD,

162

12

163

(Table S2), as described previously.32 Unique heptad-nucleotide sequences (12 bases) were added

164

to the reverse primers as barcodes to assign sequences to the different fractions. PCR was

165

performed using the method described by Song et al.33 Sequencing was conducted using 2 × 250

166

base pair paired-end technology on the Illumina MiSeq sequencer in a standard pipeline. The

167

qualified sequences were analysed after read filtering34,

C_ABA,

13

C_AD and

13

C_ABA microcosms was amplified using the 515f/806r primer set

35

and assigned using an operational

9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 36

168

taxonomic unit (OTU)-based method to generate microbiome profiles, as described previously.36-38

169

OTU assignment was based on a 97% cut-off in the present study.

170

The relative abundance of each OTU was determined, and the top 100 most relatively

171

abundant were selected for analysis according to previous studies.18, 39 The PHE degraders were

172

identified by OTUs enriched in the heavy fractions from the

173

compared with the 12C_AD and 12C_ABA samples. In the present study, we identified four OTUs

174

(OTU_56, OTU_65, OTU_106 and OTU_111) from the

175

(OTU_8, OTU_15, OTU_23, OTU_43, OTU_83, OTU_106 and OTU_111) from the

176

treatment. These sequences were further trimmed using the Greengenes database and aligned to

177

Bacillus spp. (OTU_8, accession number: MF037428), Ammoniphilus spp. (OTU_15, accession

178

number: MF037425), Hyphomicrobium spp. (OTU_23, accession number: MF037426),

179

Paenibacillus spp. (OTU_43, accession number: MF037429), Rhodoplanes spp. (OTU_56,

180

accession number: MF037430), Xanthomonadaceae spp. (OTU_65, accession number:

181

MF037431), Sporosarcina spp. (OTU_83, accession number: MF037427), Mycobacterium spp.

182

(OTU_106, accession number: MF037432), and Enterobacteriaceae spp. (OTU_111, accession

183

number: MF037433), respectively. Phylogenetic analysis of these sequences was performed as

184

described previously.4

13

13

C_AD and

13

C_ABA samples

C_AD treatment and seven OTUs 13

C_ABA

2.7. Detection of PAH-RHD genes. The PAH-RHDα genes in the heavy DNA fractions of

185 186

13

187

642f/933r) and Gram negative (610f/911r) degraders,33 respectively (Table S2). Gradient PCR and

188

the amplification reactions were performed as described previously.16 However, in the present

189

study, only the PAH-RHDα GP primer set produced a strong and specific amplicon and was

C-PHE treatment samples were amplified using two primer sets for Gram positive (GP;

10

ACS Paragon Plus Environment

Page 11 of 36

Environmental Science & Technology

190

selected for analysis. We identified two PAH-RHDα GP genes (PAH-RHDα C1 and C15) from the

191

13

C_AD treatment, and three PAH-RHDα GP genes (PAH-RHDα L1, L6 and C10) from the

192

13

C_ABA treatment. These sequences are available in GenBank with the following accession

193

numbers: PAH-RHDα C1: MF037434, PAH-RHDα C15: MF037435, PAH-RHDα L1: MF037436,

194

PAH-RHDα L6: MF037437 and PAH-RHDα L10: MF037438.

195

2.8. Chemical analysis. The concentrations of PHE and other PAHs in each microcosm were

196

analysed by gas chromatography (model 7890, Agilent, Santa Clara, CA, USA) using a capillary

197

column (DB-5MS, 30 m, 0.25 mm, 0.25 µm) and a mass spectrometric detector (model 5975,

198

Agilent) on days 0, 3 and 6, as described previously.18 Briefly, the water sample was spiked with

199

1,000 ng deuterated PAHs and extracted twice with dichloromethane. The extracted organic phase

200

was concentrated to approximately 0.5 mL and purified using a silica gel/alumina column (8 mm

201

i.d.). The eluent was concentrated to approximately 50 µL using a gentle stream of N2, and 1,000

202

ng hexamethylbenzene were added as an internal standard to all samples before instrumental

203

analysis. The components and concentrations of the deuterated PAHs, standards and the internal

204

standard are listed in Table S3.

205

11

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 36

206

3. RESULTS

207

3.1. PHE degradation. The PHE biodegradation curves in the 12C_AD, 12C_ABA, 13C_AD

208

and 13C_ABA treatments are shown in Figure S3. The recovery rates of PHE during the extraction

209

procedure were 79-85% in this study. The PHE concentration in the sterile control exhibited less

210

decrease than those in the biotic treatments, consistent with our previous observations. On day 3,

211

there was 36.7%, 38.2%, 22.6% and 23.8% residual PHE in the 12C_AD, 13C_AD, 12C_ABA and

212

13

213

treatments. No significant difference (p > 0.05) was observed between either the

214

13

215

in the

216

than those in the

217

degradation efficiency increased by 13.8% when strain LJ-5 was added to the wastewater.

218

Statistical analysis showed a significant difference between the AD and ABA treatments after 3

219

days of incubation, suggesting that ABA strategy by A. tandoii LJ-5 inoculum could significantly

220

improve the PHE biodegradation efficiency in PAH-contaminated wastewater. In addition, there

221

was little loss of other PAHs in the wastewater over the 3-day period (data not shown).

C_ABA microcosms, respectively, indicating that PHE biodegradation occurred in the biotic 12

C_AD and

C_AD treatments, or the 12C_ABA and 13C_ABA treatments. Moreover, the PHE concentrations 12

C_ABA (22.6%) and 12

13

C_ABA (23.8%) microcosms were significantly lower (p < 0.05)

C_AD (36.7%) and

13

C_AD (38.2%) microcosms, suggesting that the

222

3.2. PHE degraders revealed by DNA-SIP. DNA extracted from all four biotic treatments,

223

after 3 days of incubation, were separated by isopycnic caesium chloride density gradient

224

centrifugation, followed by high-throughput sequencing of each fraction. The relative abundance

225

of the total 16S rRNA defined by the genus indicated a slight difference in indigenous microbial

226

communities between the samples from either the

227

12

12

C_AD and

13

C_AD treatments, or the

C_ABA and 13C_ABA treatments (Figure S4). Moreover, the addition of strain LJ-5 induced a 12

ACS Paragon Plus Environment

Page 13 of 36

Environmental Science & Technology

228

significant change in the composition and structure of microbial communities in the ABA

229

treatments, compared with the AD microcosms. Obviously, the relative abundances of the genus

230

Acinetobacter (2.08%) increased significantly after addition of strain LJ-5 in the ABA treatments,

231

and was much higher than that in the AD microcosms (0.05%). The majority of the abundant

232

bacteria (> 5% total abundance), including members of the genera Kaistobacter, Rhodanobacter,

233

Burkholderia, and unclassified Xanthomonadaceae, were enriched in the AD microcosms (Figure

234

S4). The relative abundance of the genus Kaistobacter (38.4%) increased and was significantly

235

higher (p < 0.05) in ABA microcosms than those in AD microcosms (22.9%), whereas the relative

236

abundances of the genera Rhodanobacter, Burkholderia, and unclassified Xanthomonadaceae

237

were significantly lower (p < 0.05) in ABA microcosms (15.8%, 1.9% and 3.3%, respectively)

238

than those in AD microcosms (24.3%, 7.2% and 10.5%, respectively). As for rare bacteria, the

239

ABA strategy stimulated growth of the genera Conexibacter and Rhodoplanes, owing to their

240

higher relative abundances (2.09% and 2.01%) compared with the AD microcosms (0.76% and

241

1.23%, p < 0.05). No significant changes were detected in the non-PHE controls (Figure S5).

242

The bacterial 16S rRNA gene abundance was quantified by qPCR using DNA recovered

243

from each fraction of all samples as the template. As shown in Figure S6, in both AD and ABA

244

microcosms after the 3-day incubation, the bacterial 16S rRNA in fractions with higher BDs

245

(1.7382 or 1.7556 g/mL; see asterisks in Figures 1 and 2) was significantly higher in

246

treatments than those in the 12C-PHE control (marked in grey). The abundance of total 16S rRNA

247

of additive bioaugmented strain LJ-5 was 4 × 105 copies/mL, accounting for 2% of that in the

248

ABA treatment (1.96 × 107 copies/mL; Figure S7). Accordingly, the relative abundance of the

249

genus Acinetobacter after ABA treatment increased to 2.08%, and the total 16S rRNA gene copy 13

ACS Paragon Plus Environment

13

C-PHE

Environmental Science & Technology

Page 14 of 36

250

number of ABA microcosms was slightly higher than that of raw wastewater (Figure S7),

251

indicating that the bioaugmented strain LJ-5 persisted in these treatments. The indigenous microorganisms responsible for

252

13

C-PHE assimilation were detected by

253

comparing the relative abundances of specific OTUs in the 12C-PHE and 13C-PHE treatments from

254

each fraction. As shown in Figure 1, OTU_56, OTU_65, OTU_106 and OTU_111 were enriched

255

at higher BDs (1.7382 or 1.7556 g/mL) in the

256

treatments. Comparing to the relative abundances of OTU_56, OTU_65, OTU_106 and OTU_111

257

in the same fractions of the

258

the higher abundance in the heavy fractions from the 13C_AD sample (1.49%, 1.68%, 0.94% and

259

4.45%, respectively) indicated that the microorganisms represented by the above OTUs played a

260

primary role in PHE degradation. However, the inoculation of strain LJ-5 produced a significant

261

change in the diversity of the indigenous PHE-degrading communities. Seven main types of

262

bacteria represented by OTU_8, OTU_15, OTU_23, OTU_43, OTU_83, OTU_106 and OTU_111

263

at a higher BD (1.7382 or 1.7556 g/mL) were enriched in the

264

12

265

OTU_43, OTU_83, OTU_106 and OTU_111 were significantly higher in the heavy fractions of

266

the 13C_ABA treatment (0.93, 1.71, 1.26, 1.08, 0.11, 0.31 and 7.78%, respectively) than those in

267

the

268

Nevertheless, OTU_4, representing strain LJ-5, was not enriched in the heavy fractions of the

269

13

270

BDs, 1.7382 and 1.7556 g/mL) of the 13C_ABA microcosm (0.02 and 0.02%, respectively) than in

271

the 12C_ABA microcosm (0.13 and 0.05%, respectively). Although PHE-degrading communities

12

13

C_AD microcosms, but not in the

12

C_AD

C_AD treatment (0.18%, 0.15%, 0.11% and 0.82%, respectively),

13

C_ABA sample, but not in the

C_ABA sample (Figure 2). Similarly, the relative abundances of OTU_8, OTU_15, OTU_23,

12

C_ABA microcosm (0.05, 0.08, 0.01, 0.01, 0.007, 0.13 and 1.55%, respectively).

C_ABA treatment samples because of its lower relative abundance in the heavy fractions (higher

14

ACS Paragon Plus Environment

Page 15 of 36

Environmental Science & Technology

272

differed between the AD and ABA treatments, OTU_106 and OTU_111 played a vital role in PHE

273

degradation under both conditions.

274

Figure 3 shows phylogenetic information for the PHE degraders represented by the above

275

OTUs. OTU_8, OTU_15, OTU_43 and OTU_83 belong to the genus Bacillus (family

276

Bacillaceae), Ammoniphilus (family Aneurinibacillus), Paenibacillus (family Paenibacillaceae)

277

and Sporosarcina (family Planococcaceae), respectively, within the same order Bacillales

278

(phylum Firmicutes, class Bacilli). OTU_15 shares 100% similarity with Ammoniphilus resinae

279

CC-RT-E (HM193518) and forms a subclade with a high bootstrap value of 99. OTU_43 exhibits

280

98.02% similarity with the partial 16S rRNA gene sequence of strain Paenibacillus xylanisolvens

281

X11-1 (AB495094). OTU_83 has 100% similarity with strains Sporosarcina ureae DSM 2281

282

(AF202057), Sporosarcina aquimarina SW28 (AF202056), Sporosarcina saromensis HG645

283

(AB243859) and Sporosarcina luteola Y1 (AB473560) and forms a subclade with a bootstrap

284

value of 84. OTU_23 and OTU_56 are assigned to the genus Hyphomicrobium (family

285

Hyphomicrobiaceae) and Rhodoplanes (family Bradyrhizobiaceae), respectively, within the same

286

order Rhizobiales (phylum Proteobacteria, class Alphaproteobacteria). OTU_23 exhibits 98.81%

287

similarity with the partial 16S rRNA gene sequence of strain Hyphomicrobium denitrificans ATCC

288

51888 (ACVL01000012), and forms a subclade with a high bootstrap value of 96. OTU_56 shares

289

100% similarity to the partial 16S rRNA gene sequence of uncultured Rhodoplanes sp. clone

290

Leob163 (KF226093.1), and forms a subclade with a high bootstrap value of 98. OTU_65 and

291

OTU_111 are classified in the family Xanthomonadaceae (order Xanthomonadales) and

292

Enterobacteriaceae

293

Betaproteobacteria (phylum Proteobacteria). Additionally, OTU_106 is characterised as genus

(order

Enterobacteriales),

respectively,

15

ACS Paragon Plus Environment

within

the

same

class

Environmental Science & Technology

Page 16 of 36

294

Mycobacterium (phylum Actinobacteria, class Actinobacteria, order Corynebacteriales, family

295

Mycobacteriaceae) and shares 100% similarity with many strains in this genus, including

296

Mycobacterium brisbanense ATCC 49938 (AY012577), Mycobacterium elephantis 484

297

(AJ010747), Mycobacterium moriokaense DSM 44221 (AJ429044) and Mycobacterium pulveris

298

DSM 44222 (AJ429046), and forms a subclade with a high bootstrap value of 100.

299

3.3. Presence of PAH-RHDα genes involved in PHE metabolism. In the present study, the 13

13

300

PAH-RHDα GP genes were analyzed in the heavy fractions of the

301

treatments (marked with an asterisk in Figures 1 & 2). PAH-RHDα C1 and C15 were detected in

302

the heavy fractions of

303

conspicuous change in the types of PAH-RHDα GP genes, compared with the AD microcosms,

304

and PAH-RHDα L1, L6 and L10 were detected (Figure 4). Thereinto, PAH-RHDα C1 and L1 share

305

100% similarity and are the same PAH-RHDα GP gene, indicating that they exist in the heavy

306

fractions of both the

307

bootstrap value of 96 in the phylogenetic tree of amplified PAH-RHDα GP genes, and show 89%

308

similarity with the PAH-RHD genes of Mycobacterium sp. S23 (ALS30442.1). The other

309

PAH-RHD gene in the heavy fractions of the 13C_AD treatment, PAH-RHDα C15, exhibits 99%

310

similarity with the PAH-RHD genes of uncultured bacterium (AMM73080.1). In the heavy

311

fractions of the

312

genes of Mycobacterium vanbaalenii PYR-1 (AAY85176.1), Mycobacterium novocastrense

313

(GAT12202.1) and Terrabacter sp. FLO (ABA87073.1), and PAH-RHDα L10 shows 100%

314

similarity with the PAH-RHD genes of Bacillus sp. CL1-1 (AND66078.1), Enterobacter sp.

315

CL1-2 (AND66079.1) and uncultured Comamonas sp. (CAN85224.1), and forms a subclade with

13

13

13

C_AD treatment (Figure 4). The

C_AD and

13

13

C_AD and

C_ABA

C_ABA treatment produced a

C_ABA microcosms. They form a subclade with a high

C_ABA treatment, PAH-RHDα L6 reveals 99% similarity with the PAH-RHD

16

ACS Paragon Plus Environment

Page 17 of 36

Environmental Science & Technology

316

a bootstrap value of 61.

17

ACS Paragon Plus Environment

Environmental Science & Technology

317

Page 18 of 36

4. DISCUSSION

318

DNA-SIP has been used successfully to demonstrate that native microorganisms collected at

319

field sites are involved in PHE biodegradation.20, 27 Gutierrez et al. used DNA-SIP with fully

320

13

321

PHE in surface slicks and plume waters.27 Jones et al. performed DNA-SIP with 13C-labelled PHE

322

as part of a larger project investigating strategies for bioremediation of PAH-contaminated soil

323

from a former manufactured-gas plant site.20 Our work applied

324

DNA-SIP and revealed the mechanisms of ABA strategy in PAH-contaminated wastewater with

325

the addition of the autochthonous microorganism A. tandoii LJ-5. Although ABA strategies have

326

been successfully applied for remediation of PAH-contaminated sites,14, 15 limited studies have

327

addressed the changes that occur in the microbial community during the ABA process.14, 40 In the

328

present study, ABA with strain LJ-5 inoculum produced a significant increase in PHE

329

biodegradation efficiency in PAH-contaminated wastewater. This suggests that strain LJ-5 is

330

potentially an ABA agent encouraging the remediation of PAH-contaminated sites. Moreover, the

331

addition of LJ-5 markedly modified the bacterial community structure. By comparing the bacterial

332

community structures of microbial communities between AD and ABA, we identified three

333

bacterial genera, affiliated with Kaistobacter, Conexibacter and Rhodoplanes, that were

334

stimulated by LJ-5 inoculation. Kaistobacter demonstrated relatively high abundance in the

335

indigenous microbial communities of a tailings dump contaminated with antimony41 and the

336

ability to biodegrade both S-ethyldipropylthiocarbamate and atrazine in soils.42 Members of the

337

genus Conexibacter within the class Actinobacteria are affiliated with aromatic hydrocarbon

338

degraders and contain bph genes encoding the biphenyl degradative pathway.43 Rhodoplanes has

C-labelled PHE to link the phylogenetic identity of bacterial taxa with their ability to mineralize

18

ACS Paragon Plus Environment

13

C-PHE as the substrate in

Page 19 of 36

Environmental Science & Technology

339

been proposed to accommodate taxa that are primarily phototrophic and dominant in rhizospheric

340

soils and roots.44 Furthermore, the addition of LJ-5 also affected the composition of abundant

341

bacteria, including members of the genera Kaistobacter, Rhodanobacter, Burkholderia, and

342

unclassified Xanthomonadaceae. Among them, Rhodanobacter spp. are dominant in crude oil plus

343

dispersant and rhizospheres of maize cultivars,45, 46 and Burkholderia and unclassified phylotypes

344

within the family Xanthomonadaceae such as Stenotrophomonas, are the dominant genera of PAH

345

degraders with high metabolic activities in PAH-contaminated soils and sediments.47, 48 However,

346

there was no previous evidence directly linking the above activated or abundant bacteria (except

347

Xanthomonadaceae and Burkholderia) to PHE degradation. They might play important roles in

348

stabilizing the microbial community or degrading other organic pollutants in PAH-contaminated

349

wastewater.18

350

Besieds the structure and dynamics of the microbial community, ABA might also shape the

351

composition of functional PHE-degraders community, although this has not yet been reported. Our

352

study is the first to show that ABA also influences the abundance and diversity of PAH-degrading

353

bacteria in PAH-polluted wastewater, suggesting that different PHE degradation activities can be

354

achieved by distinct microbial communities within the same environment. In the AD treatment,

355

the indigenous microorganisms responsible for PHE degradation were affiliated with

356

Rhodoplanes, Mycobacterium, Xanthomonadaceae (genus unclassified) and Enterobacteriaceae

357

(genus unclassified). The genus Rhodoplanes classified under the family Hyphomicrobiaceae of

358

the order Rhizobiales in the class Alphaproteobacteria, was first described in 1994.49 Members of

359

this genus are characterized as Gram negative, phototrophic, rod-shaped, motile and are

360

widespread in aquatic habitats.50 Rhodoplanes is the dominant genus in wastewater treatment 19

ACS Paragon Plus Environment

Environmental Science & Technology

361

plants,51 and some strains in this genus can grow on organic carbon sources such as acetate,

362

pyruvate, glucose and fructose.50 However, there is still no proof of its ability to degrade PAHs,

363

including PHE, and their roles in PAH-contaminated sites remain unclear. Our results provide

364

unequivocal evidence that some microbes in this genus are primarily responsible for PHE

365

degradation in PAH-polluted wastewater. Mycobacterium is well known to degrade various

366

environmental contaminants such as chlorinated compounds, polychlorobiphenyls and PAHs.52

367

Degradation of PHE by Mycobacterium has been reported previously,53 and several strains in this

368

genus can effectively degrade high molecular weight PAHs.54 However, prior to this study, PHE

369

degradation by indigenous Mycobacterium has not been identified using DNA-SIP. The families

370

Xanthomonadaceae and Enterobacteriaceae are members of the class Betaproteobacteria.

371

Previous studies have reported that Stenotrophomonas sp. and Pseudoxanthomonas sp. in the

372

family Xanthomonadaceae possess the functions of metabolizing a wide range of PAHs, e.g.

373

naphthalene, PHE, anthracene, fluorene, pyrene and benzo[a]pyrene.55-57 Members of the family

374

Enterobacteriaceae, such as Raoultella sp. and Klebsiella sp., can also degrade PAHs efficiently,

375

including acenaphthene, fluorene, pyrene, benzo[a]pyrene, PHE and fluoranthene.58-60 Here, our

376

results provide strong evidence that some microbes in these families are active PHE degraders in

377

PAH-contaminated wastewater.

378

It should be noted that the bacterial community structure of the wastewater in this study was

379

significantly different from that of our previously reported wastewater,18 thus resulting in the

380

change in diversity of microorganisms responsible for in situ PHE degradation. The majority of

381

the abundant bacteria in wastewater collected in September 2015 were affiliated with the genera

382

Pseudomonas (23.5 ± 1.2%), unclassified Chitinophagaceae (4.7 ± 0.3%) and Comamonadaceae 20

ACS Paragon Plus Environment

Page 20 of 36

Page 21 of 36

Environmental Science & Technology

383

(6.4 ± 0.3%).18 We identified four OTUs directly responsible for indigenous PHE biodegradation,

384

including the phylotypes affiliated with Acinetobacter, Sphingobium, Kouleothrix and

385

Sandaracinobacter. In the present study, among the PHE degraders identified in the AD

386

microcosms, only two taxa Mycobacterium and Enterobacteriaceae (genus unclassified) were also

387

identified as PHE degraders in the ABA treatments with strain LJ-5. Five new indigenous

388

PHE-degrading organisms became active, namely phylotypes affiliated with Bacillus,

389

Paenibacillus, Ammoniphilus, Sporosarcina and Hyphomicrobium. The genera Bacillus and

390

Paenibacillus are known to be metabolically versatile, degrading aromatic and hydroxylated

391

aromatic compounds.55, 61, 62 Moreover, these two genera can degrade PHE,55 and appear to be

392

highly competitive by the addition of root exudates.63 However, no studies have considered the

393

roles of Bacillus and Paenibacillus in the remediation of PAH-contaminated sites by ABA

394

strategy. The genus Ammoniphilus is a member of the family Aneurinibacillus in the phylum

395

Firmicutes. Members of this genus are typically Gram variable, oxidase- and catalase-positive,

396

obligately oxalotrophic, and motile by peritrichous flagella.64 Strains of the genus Ammoniphilus

397

are able to utilize high concentrations of ammonium ions.65 Similar to Ammoniphilus,

398

Sporosarcina also belongs to the phylum Firmicutes. It is usually detected in clinical specimens

399

and raw cow’s milk.66 Sporosarcina has been shown to principally degrade aflatoxin B-1, which is

400

the most potent naturally occurring carcinogen known.67 Hyphomicrobium, belonging to the

401

Alphaproteobacteria and Hyphomicrobiaceae family, has versatile metabolic capabilities such as

402

the ability to degrade dichloromethane,68 methamidophos,69 dimethyl sulphide and methanol.70

403

However, the taxa Ammoniphilus, Sporosarcina and Hyphomicrobium have not been linked

404

previously with PHE degradation; thus, their mechanisms of degradation in PAH-contaminated 21

ACS Paragon Plus Environment

Environmental Science & Technology

405

sites remain unknown. Our results provide strong evidence that members of the above genera are

406

the primary indigenous PHE degraders induced by ABA strategy in wastewater, expanding our

407

knowledge of the diversity of indigenous PHE-degrading communities.

408

Despite the satisfactory performance of ABA strategies in bioremediation by the addition of

409

autochthonous degraders, no studies have confirmed the activities or degrading capacity of the

410

reintroduced strains. It is worth mentioning that, although the ABA strategy with LJ-5 inoculum

411

encouraged PHE biodegradation, the participation of A. tandoii LJ-5 in PHE degradation in situ

412

was questioned due to their limited enrichment in the heavy DNA fraction according to the

413

DNA-SIP results. Previous researchers have focused primarily on the effects of reintroduced

414

degraders on degradation efficiency and the total microbial activities or population during ABA.15,

415

29

416

strain to the polluted sites.40,

417

autochthonous strain LJ-5 in ABA strategy, modifying the diversity of indigenous PHE degraders

418

instead of participating in in situ PAH degradation. Such a mechanism has never been addressed

419

before and remains unclear in PAH-contaminated wastewater. One possible reason might be

420

attributed to the change in wastewater composition, consequently leaving strain LJ-5 from

421

functional PHE-degraders in previous wastewater to inactive microbes in the present study.

422

Another reason might be the contribution of Acinetobacter biosurfactants to the bioremediation of

423

PAHs and the shift in community composition, which has been confirmed by other researchers.72

424

For instance, a wide range of Acinetobacter sp. can produce biosurfactants (amphiphilic molecules

425

consisting of hydrophilic and hydrophobic domains)73, which might increase the bioavailability of

426

hydrophobic compounds by solubilisation and/or emulsification or alter the cell surface properties

The diversity and functions of the whole microbial community can be enhanced by adding a 71

Our work is the first study to uncover the roles of the

22

ACS Paragon Plus Environment

Page 22 of 36

Page 23 of 36

Environmental Science & Technology

427

of other microorganisms.72 In the present study, A. tandoii LJ-5 could reduce the surface tension

428

of water from 70.0 to 35.2 mN/m. Our results indicated that strain LJ-5 might be capable of

429

producing biosurfactants, simultaneously improving PHE bioavailability and activating the

430

functions of other PHE-degrading candidates. Moreover, the relative abundance of the

431

bioaugmented strain LJ-5 was 4×105 copies/mL, accounting for 2% of the total 16S rRNA in ABA

432

treatments and approximately 2 times as that of the initial inoculum abundance in surface tension

433

measurements. It suggested sufficient amount of strain LJ-5 in all the bioaugmentation treatments

434

to enhance PHE degradation via the production of biosurfactants.

435

More evidence of a change in the PHE-degrading community was provided by the sequences

436

of PAH-RHDα genes between AD and ABA treatments. Only one novel PAH-RHDα GP gene,

437

PAH-RHDα C1 and L1, was detected in the heavy DNA fraction of both the

438

13

439

(genus unclassified, OTU_111) which were also both detected as functional PHE-degrading

440

bacteria in the two microcosms. In addition, the distinctive PAH-RHDα L15 gene is potentially the

441

functional gene from Rhodoplanes or Enterobacteriaceae (genus unclassified) as they were all

442

detected in the heavy DNA fraction of the

443

microcosm, PAH-RHDα L6 and L10 genes might be linked to the PHE degraders of Bacillus,

444

Paenibacillus, Ammoniphilus, Sporosarcina or Hyphomicrobium. Here, the active PHE degraders

445

might possess other functional genes that were not targeted by the primers used in this study, and

446

we could not accurately attribute one functional gene to one bacterium for the lack of available

447

database information. Indeed, little information is available regarding the functional genes of

448

certain dominant PHE-degrading genera, such as Rhodoplanes, Ammoniphilus, Sporosarcina and

13

C_AD and

C_ABA treatments, possibly associated with Mycobacterium (OTU_106) or Xanthomonadaceae

13

C_AD treatment. Accordingly, in the

23

ACS Paragon Plus Environment

13

C_ABA

Environmental Science & Technology

449

Hyphomicrobium, since this is first report of PHE biodegradation for these genera.

450

ABA by A. tandoii LJ-5 was first applied as a potential strategy to enhance the remediation of

451

PAH-contaminated wastewater in this study. Besides resulting in a significant increase in the PHE

452

biodegradation efficiency, ABA remarkably modified the functional PHE degrading community of

453

the wastewater. The indigenous microorganisms responsible for PHE degradation were affiliated

454

with Rhodoplanes, Mycobacterium, Xanthomonadaceae and Enterobacteriaceae in the AD

455

treatments, whereas five new taxa (Bacillus, Paenibacillus, Ammoniphilus, Sporosarcina and

456

Hyphomicrobium) were acitvated in in situ PHE biodegradation in ABA microcosms. Of all the

457

above PHE degraders, Rhodoplanes, Ammoniphilus, Sporosarcina and Hyphomicrobium were

458

linked to indigenous PHE biodegradation for the first time. Nevertheless, LJ-5 did not participate

459

in indigenous PHE degradation. The change in PAH-RHDα gene diversity further confirmed our

460

findings by different PAH-RHDɑ genes involved in PHE metabolism in the AD and ABA

461

treatments. Collectively, our findings raise a new mechanism of ABA, provide new insights into

462

the diversity of PHE-degrading communities, and suggest ABA as a promising in situ strategy for

463

PAH-contaminated wastewater.

24

ACS Paragon Plus Environment

Page 24 of 36

Page 25 of 36

Environmental Science & Technology

464

ASSOCIATED CONTENT

465

Supporting Information

466

(Table S1) Concentrations of PAHs in wastewater; (Table S2) primers used for the PCR of

467

16S rRNA and PAH-RHD gene; (Table S3) the components of deuterated PAHs, standards

468

and internal standard; (Figure S1) correlation between the fraction number and buoyant

469

density (BD; g/mL) of DNA extracted from (a)

470

12

471

(ng/µL) and buoyant density (g/mL) from DNA extracted from (a)

472

treatments, and (b) 12C-ABA and 13C-ABA treatments; (Figure S3) residual PHE percentage

473

in wastewater after 3 days of incubation; (Figure S4) relative abundance of 16S rRNA

474

defined taxa by genus in 12C_AD, 13C_AD, 12C_ABA, and 13C_ABA microcosms; (Figure S5)

475

relative abundance of 16S rRNA defined taxa by genus in non-PHE microcosms; (Figure S6)

476

the quantitative distribution of density-resolved bacterial 16S rDNA obtained from

477

wastewater samples in (a)

478

treatments; (Figure S7) the 16S rRNA gene copies in

479

13

480

Information

C-ABA and

13

12

C-AD and

13

C-AD treatments, and (b)

C-ABA treatments; (Figure S2) correlation between DNA concentration

12

C-AD and

13

C-AD treatments, and (b) 12

C_AD,

12

13

12

C-AD and

C-ABA and

C_AD,

12

13

C-AD

13

C-ABA

C_ABA, and

C_ABA microcosms. There are 12 pages, 3 tables and 7 figures in the Supporting

481 482

ACKNOWLEDGMENTS

483

Financial support was provided by the Scientific and Technological Planning Project of

484

Guangzhou, China (Nos. 201707020034), the National Natural Science Foundation of China (No.

485

41673111), and the Department of Science and Technology of Guangdong province

486

(2016TQ03Z938).

487

25

ACS Paragon Plus Environment

Environmental Science & Technology

488

REFERENCES

489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530

1. Hamdi, H.; Benzarti, S.; Manusadžianas, L.; Aoyama, I.; Jedidi, N., Bioaugmentation and biostimulation effects on PAH dissipation and soil ecotoxicity under controlled conditions. Soil Biol Biochem 2007, 39, (8), 1926-1935. 2. Xue, W.; Warshawsky, D., Metabolic activation of polycyclic and heterocyclic aromatic hydrocarbons and DNA damage: a review. Toxicol. Appl. Pharmacol. 2005, 206, (1), 73-93. 3. Keith, L.; Telliard, W., Priority pollutants Ia perspective view. Environ Sci Technol 1979 13 (4), 416-423. 4. Jiang, L.; Song, M.; Luo, C.; Zhang, D.; Zhang, G., Novel Phenanthrene-Degrading Bacteria Identified by DNA-Stable Isotope Probing. PloS one 2015, 10, (6), e0130846. 5. Bento, F. M.; Camargo, F. A.; Okeke, B. C.; Frankenberger, W. T., Comparative bioremediation of soils contaminated with diesel oil by natural attenuation, biostimulation and bioaugmentation. Bioresour Technol 2005, 96, (9), 1049-1055. 6. Edgehill, R. U., Bioremediation by inoculation with microorganisms. Bioremediation of Contaminated Soils 1999, 289-313. 7. Hosokawa, R.; Nagai, M.; Morikawa, M.; Okuyama, H., Autochthonous bioaugmentation and its possible application to oil spills. World Journal of Microbiology and Biotechnology 2009, 25, (9), 1519-1528. 8. Krüger, U. S.; Johnsen, A. R.; Burmølle, M.; Aamand, J.; Sørensen, S. R., The potential for bioaugmentation of sand filter materials from waterworks using bacterial cultures degrading 4‐chloro‐2‐methylphenoxyacetic acid. Pest Manage. Sci. 2015, 71, (2), 257-265. 9. Yao, L.; Teng, Y.; Luo, Y.; Christie, P.; Ma, W.; Liu, F.; Wu, Y.; Luo, Y.; Li, Z., Biodegradation of Polycyclic Aromatic Hydrocarbons (PAHs) byTrichoderma reeseiFS10-C and Effect of Bioaugmentation on an Aged PAH-Contaminated Soil. Biorem. J. 2015, 19, (1), 9-17. 10. Owsianiak, M.; Dechesne, A.; Binning, P. J.; Chambon, J. C.; Sørensen, S. R.; Smets, B. F., Evaluation of bioaugmentation with entrapped degrading cells as a soil remediation technology. Environ Sci Technol 2010, 44, (19), 7622-7627. 11. Schultz-Jensen, N.; Aamand, J.; Sorensen, S. R., Bioaugmentation potential of free and formulated 2,6-dichlorobenzamide (BAM) degrading Aminobacter sp. MSH1 in soil, sand and water. AMB Express 2016, 6, (1), 33. 12. Ueno, A.; Ito, Y.; Yumoto, I.; Okuyama, H., Isolation and characterization of bacteria from soil contaminated with diesel oil and the possible use of these in autochthonous bioaugmentation. World Journal of Microbiology and Biotechnology 2007, 23, (12), 1739-1745. 13. Nikolopoulou, M.; Pasadakis, N.; Kalogerakis, N., Evaluation of autochthonous bioaugmentation and biostimulation during microcosm-simulated oil spills. Marine Pollution Bulletin 2013, 72, (1), 165-173. 14. Madueño, L.; Alvarez, H. M.; Morelli, I. S., Autochthonous bioaugmentation to enhance phenanthrene degradation in soil microcosms under arid conditions. 26

ACS Paragon Plus Environment

Page 26 of 36

Page 27 of 36

Environmental Science & Technology

531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574

International Journal of Environmental Science and Technology 2014, 12, (7), 2317-2326. 15. Madueño, L.; Coppotelli, B. M.; Alvarez, H. M.; Morelli, I. S., Isolation and characterization of indigenous soil bacteria for bioaugmentation of PAH contaminated soil of semiarid Patagonia, Argentina. Int. Biodeterior. Biodegrad. 2011, 65, (2), 345-351. 16. Cebron, A.; Norini, M. P.; Beguiristain, T.; Leyval, C., Real-Time PCR quantification of PAH-ring hydroxylating dioxygenase (PAH-RHDalpha) genes from Gram positive and Gram negative bacteria in soil and sediment samples. J. Microbiol. Methods 2008, 73, (2), 148-159. 17. Amann, R. I.; Ludwig, W.; Schleifer, K. H., Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiological Reviews 1995, 59, (59), 143-169. 18. Li, J.; Luo, C.; Song, M.; Dai, Q.; Jiang, L.; Zhang, D.; Zhang, G., Biodegradation of Phenanthrene in Polycyclic Aromatic Hydrocarbon-Contaminated Wastewater Revealed by Coupling Cultivation-Dependent and -Independent Approaches. Environmental Science and Technology 2017, 51, 3391-3401. 19. Oren, A., Prokaryote diversity and taxonomy: current status and future challenges. Philosophical transactions of the Royal Society of London. Series B, Biological sciences 2004, 359, (1444), 623-638. 20. Jones, M. D.; Crandell, D. W.; Singleton, D. R.; Aitken, M. D., Stable-isotope probing of the polycyclic aromatic hydrocarbon-degrading bacterial guild in a contaminated soil. Environmental Microbiology 2011, 13, (10), 2623-2632. 21. Huang, W. E.; Ferguson, A.; Singer, A. C.; Lawson, K.; Thompson, I. P.; Kalin, R. M.; Larkin, M. J.; Bailey, M. J.; Whiteley, A. S., Resolving genetic functions within microbial populations: in situ analyses using rRNA and mRNA stable isotope probing coupled with single-cell raman-fluorescence in situ hybridization. Appl. Environ. Microbiol. 2009, 75, (1), 234-241. 22. Singleton, D. R.; Hunt, M.; Powell, S. N.; Frontera-Suau, R.; Aitken, M. D., Stable-isotope probing with multiple growth substrates to determine substrate specificity of uncultivated bacteria. J. Microbiol. Methods 2007, 69, (1), 180-187. 23. Gutierrez, T., Identifying polycyclic aromatic hydrocarbon-degrading bacteria in oil-contaminated surface waters at Deepwater Horizon by cultivation, stable isotope probing and pyrosequencing. Reviews in Environmental Science and Bio/Technology 2011, 10, (4), 301-305. 24. Dumont, M. G.; Murrell, J. C., Stable isotope probing - linking microbial identity to function. Nature Reviews Microbiology 2005, 3, (6), 499-504. 25. Jiang, B.; Jin, N.; Xing, Y.; Su, Y.; Zhang, D., Critical review on unraveling uncultivable pesticide degraders via stable isotope probing (SIP). Critical Reviews in Biotechnology 2018, DOI:10.1080/07388551.07382018.01427697. 26. Jeon, C. O.; Park, W.; Padmanabhan, P.; DeRito, C.; Snape, J. R.; Madsen, E. L., Discovery of a bacterium, with distinctive dioxygenase, that is responsible for in situ biodegradation in contaminated sediment. Proceedings of the National Academy of Sciences of the United States of America 2003, 100, (23), 13591-13596. 27

ACS Paragon Plus Environment

Environmental Science & Technology

575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618

27. Gutierrez, T.; Singleton, D. R.; Berry, D.; Yang, T.; Aitken, M. D.; Teske, A., Hydrocarbon-degrading bacteria enriched by the Deepwater Horizon oil spill identified by cultivation and DNA-SIP. The ISME journal 2013, 7, (11), 2091-2104. 28. Jones, M. D.; Singleton, D. R.; Carstensen, D. P.; Powell, S. N.; Swanson, J. S.; Pfaender, F. K.; Aitken, M. D., Effect of incubation conditions on the enrichment of pyrene-degrading bacteria identified by stable-isotope probing in an aged, PAH-contaminated soil. Microb. Ecol. 2008, 56, (2), 341-349. 29. Kuppusamy, S.; Thavamani, P.; Megharaj, M.; Naidu, R., Bioaugmentation with novel microbial formula vs. natural attenuation of a long-term mixed contaminated soil—treatability studies in solid-and slurry-phase microcosms. Water, Air, Soil Pollut. 2016, 227, (1), 25. 30. Li, J.; Lin, T.; Lu, Q.; Wang, J. J.; Liao, C.; Pan, Y.; Zhao, Y., Changes in physicochemical properties and bactericidal efficiency of acidic electrolyzed water ice and available chlorine decay kinetics during storage. LWT-Food Science and Technology 2014, 59, (1), 43-48. 31. Song, M.; Jiang, L.; Zhang, D.; Luo, C.; Wang, Y.; Yu, Z.; Yin, H.; Zhang, G., Bacteria capable of degrading anthracene, phenanthrene, and fluoranthene as revealed by DNA based stable-isotope probing in a forest soil. J Hazard Mater 2016, 308, 50-57. 32. Bates, S. T.; Berglyons, D.; Caporaso, J. G.; Walters, W. A.; Knight, R.; Fierer, N., Examining the global distribution of dominant archaeal populations in soil. ISME J 2011, 5, (5), 908. 33. Song, M.; Luo, C.; Jiang, L.; Zhang, D.; Wang, Y.; Zhang, G., Identification of Benzo[a]pyrene-Metabolizing Bacteria in Forest Soils by Using DNA-Based Stable-Isotope Probing. Applied & Environmental Microbiology 2015, 81, (21), 7368-7376. 34. Caporaso, J. G.; Kuczynski, J.; Stombaugh, J., QIIME allows analysis of high-throughput community sequencing data. Nature Methods 2010, 7, (5), 335–336. 35. Schloss, P. D.; Westcott, S. L.; Ryabin, T.; Hall, J. R.; Hartmann, M.; Hollister, E. B.; Lesniewski, R. A.; Oakley, B. B.; Parks, D. H.; Robinson, C. J., Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Applied & Environmental Microbiology 2009, 75, (23), 7537-7541. 36. Edgar, R. C., Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, (19), 2460-2461. 37. Mcdonald, D.; Price, M. N.; Goodrich, J.; Nawrocki, E. P.; Desantis, T. Z.; Probst, A.; Andersen, G. L.; Knight, R.; Hugenholtz, P., An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. Isme Journal 2012, 6, (3), 610-618. 38. Werner, J. J.; Koren, O.; Hugenholtz, P.; Desantis, T. Z.; Walters, W. A.; Caporaso, J. G.; Angenent, L. T.; Knight, R.; Ley, R. E., Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. Isme Journal 2012, 6, (1), 94-103. 39. Sun, W.; Xie, S.; Luo, C.; Cupples, A. M., Direct link between toluene 28

ACS Paragon Plus Environment

Page 28 of 36

Page 29 of 36

Environmental Science & Technology

619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662

degradation in contaminated-site microcosms and a Polaromonas strain. Appl. Environ. Microbiol. 2010, 76, (3), 956-959. 40. Teng, Y.; Luo, Y.; Sun, M.; Liu, Z.; Li, Z.; Christie, P., Effect of bioaugmentation by Paracoccus sp. strain HPD-2 on the soil microbial community and removal of polycyclic aromatic hydrocarbons from an aged contaminated soil. Bioresour. Technol. 2010, 101, (10), 3437-3443. 41. Xiao, E. Z.; Krumins, V.; Dong, Y. R.; Xiao, T. F.; Ning, Z. P.; Xiao, Q. X.; Sun, W. M., Microbial diversity and community structure in an antimony-rich tailings dump. Appl. Microbiol. Biotechnol. 2016, 100, (17), 7751-7763. 42. Liu, X.; Zhang, S.; Jiang, Q.; Bai, Y.; Shen, G.; Li, S.; Ding, W., Using community analysis to explore bacterial indicators for disease suppression of tobacco bacterial wilt. Scientific reports 2016, 6, 36773. 43. Brossi, M.; Mendes, L. W.; Germano, M. G.; Lima, A. B.; Tsai, S. M., Assessment of Bacterial bph Gene in Amazonian Dark Earth and Their Adjacent Soils. PloS one 2014, 9, (6), e99597. 44. Tan, Y.; Cui, Y.; Li, H.; Kuang, A.; Li, X.; Wei, Y.; Ji, X., Diversity and composition of rhizospheric soil and root endogenous bacteria in Panax notoginseng during continuous cropping practices. Journal of Basic Microbiology 2017, 57, (4), 337-344. 45. Al-Jawasim, M.; Yu, K. W.; Park, J. W., Synergistic effect of crude oil plus dispersant on bacterial community in a louisiana salt marsh sediment. FEMS Microbiol. Lett. 2015, 362, (17), fnv144. 46. Wen, X.; Wang, M.; Ti, J. S.; Wu, Y.; Chen, F., Bacterial community composition in the rhizosphere of maize cultivars widely grown in different decades. Biol. Fertility Soils 2017, 53, (2), 221-229. 47. Kuppusamy, S.; Thavamani, P.; Megharaj, M.; Venkateswarlu, K.; Lee, Y. B.; Naidu, R., Pyrosequencing analysis of bacterial diversity in soils contaminated long-term with PAHs and heavy metals: Implications to bioremediation. J. Hazard. Mater. 2016, 317, 169-179. 48. Muangchinda, C.; Chavanich, S.; Viyakarn, V.; Watanabe, K.; Imura, S.; Vangnai, A. S.; Pinyakong, O., Abundance and diversity of functional genes involved in the degradation of aromatic hydrocarbons in Antarctic soils and sediments around Syowa Station. Environmental Science and Pollution Research 2015, 22, (6), 4725-4735. 49. Hiraishi, A.; Ueda, Y., Rhodoplanes gen. nov., a New Genus of Phototrophic Bacteria Including Rhodopseudomonas rosea as Rhodoplanes roseus comb. nov. and Rhodoplanes elegans sp. nov. International Journal of Systematic and Evolutionary Microbiology 1994, 44, (4), 665-673. 50. Srinivas, A.; Sasikala, C.; Ramana Ch, V., Rhodoplanes oryzae sp. nov., a phototrophic alphaproteobacterium isolated from the rhizosphere soil of paddy. Int. J. Syst. Evol. Microbiol. 2014, 64, (7), 2198-2203. 51. Ma, Q.; Qu, Y. Y.; Shen, W. L.; Zhang, Z. J.; Wang, J. W.; Liu, Z. Y.; Li, D. X.; Li, H. J.; Zhou, J. T., Bacterial community compositions of coking wastewater treatment plants in steel industry revealed by Illumina high-throughput sequencing. Bioresour. Technol. 2015, 179, 436-443. 29

ACS Paragon Plus Environment

Environmental Science & Technology

663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706

52. Hormisch, D.; Hormisch, D.; Brost, I.; Kohring, G. W.; Giffhorn, F.; Kroppenstedt, R. M.; Stackebradt, E.; Färber, P.; Holzapfel, W. H., Mycobacterium fluoranthenivorans sp. nov., a Fluoranthene and Aflatoxin B1 Degrading Bacterium from Contaminated Soil of a Former Coal Gas Plant. Systematic and Applied Microbiology 2004, 27, (6), 653-660. 53. Hennessee, C. T.; Li, Q. X., Effects of Polycyclic Aromatic Hydrocarbon Mixtures on Degradation, Gene Expression, and Metabolite Production in Four Mycobacterium Species. Appl. Environ. Microbiol. 2016, 82, (11), 3357-3369. 54. Kwak, Y.; Li, Q. X.; Shin, J. H., Draft genome sequence of Mycobacterium rufum JS14(T), a polycyclic-aromatic-hydrocarbon-degrading bacterium from petroleum-contaminated soil in Hawaii. Standards in genomic sciences 2016, 11, 47. 55. Cebron, A.; Louvel, B.; Faure, P.; France-Lanord, C.; Chen, Y.; Murrell, J. C.; Leyval, C., Root exudates modify bacterial diversity of phenanthrene degraders in PAH-polluted soil but not phenanthrene degradation rates. Environ. Microbiol. 2011, 13, (3), 722-736. 56. Patel, V.; Cheturvedula, S.; Madamwar, D., Phenanthrene degradation by Pseudoxanthomonas sp. DMVP2 isolated from hydrocarbon contaminated sediment of Amlakhadi canal, Gujarat, India. J. Hazard. Mater. 2012, 201, (1), 43-51. 57. Arulazhagan, P.; Al-Shekri, K.; Huda, Q.; Godon, J.; Basahi, J.; Jeyakumar, D., Biodegradation of polycyclic aromatic hydrocarbons by an acidophilic Stenotrophomonas maltophilia strain AJH1 isolated from a mineral mining site in Saudi Arabia. Extremophiles 2017, 21, (1), 163-174. 58. Alegbeleye, O. O.; Opeolu, B. O.; Jackson, V., Bioremediation of polycyclic aromatic hydrocarbon (PAH) compounds: (acenaphthene and fluorene) in water using indigenous bacterial species isolated from the Diep and Plankenburg rivers, Western Cape, South Africa. Braz. J. Microbiol. 2016, 48, (2), 314-325. 59. Ping, L. F.; Zhang, C. R.; Zhang, C. P.; Zhu, Y. H.; He, H. M.; Wu, M.; Tang, T.; Li, Z.; Zhao, H., Isolation and characterization of pyrene and benzo a pyrene-degrading Klebsiella pneumonia PL1 and its potential use in bioremediation. Appl. Microbiol. Biotechnol. 2014, 98, (8), 3819-3828. 60. Xu, X. Y.; Chen, X.; Su, P.; Fang, F.; Hu, B. B., Biodegradation potential of polycyclic aromatic hydrocarbons by bacteria strains enriched from Yangtze River sediments. Environ. Technol. 2016, 37, (5), 513-520. 61. Raj, A.; Reddy, M. M. K.; Chandra, R., Identification of low molecular weight aromatic compounds by gas chromatography-mass spectrometry (GC-MS) from kraft lignin degradation by three Bacillus sp. Int. Biodeterior. Biodegrad. 2007, 59, (4), 292-296. 62. Daane, L. L.; Harjono, I.; Barns, S. M.; Launen, L. A.; Palleroni, N. J.; Haggblom, M. M., PAH-degradation by Paenibacillus spp. and description of Paenibacillus naphthalenovorans sp nov., a naphthalene-degrading bacterium from the rhizosphere of salt marsh plants. Int. J. Syst. Evol. Microbiol. 2002, 52, 131-139. 63. Haichar, F. Z.; Marol, C.; Berge, O.; Rangelcastro, J. I.; Prosser, J. I.; Balesdent, J.; Heulin, T.; Achouak, W., Plant host habitat and root exudates shape soil bacterial community structure. ISME J 2008, 2, (12), 1221-1230. 30

ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36

Environmental Science & Technology

707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742

64. Zaitsev, G. M.; Tsitko, I. V.; Rainey, F. A.; Trotsenko, Y. A.; Uotila, J. S.; Stackebrandt, E.; Salkinoja-Salonen, M. S., New aerobic ammonium-dependent obligately oxalotrophic bacteria: description of Ammoniphilus oxalaticus gen. nov., sp. nov. and Ammoniphilus oxalivorans gen. nov., sp. nov. International Journal of Systematic Bacteriology 1998, 48, (1), 151. 65. Lin, S.-Y.; Wu, Y.-H.; Hameed, A.; Liu, Y.-C.; Young, C.-C., Ammoniphilus resinae sp. nov., an endospore-forming bacterium isolated from resin fragments. International Journal of Systematic and Evolutionary Microbiology 2016, 66, (8), 3010-3016. 66. Wolfgang, W. J.; Coorevits, A.; Cole, J. A.; De Vos, P.; Dickinson, M. C.; Hannett, G. E.; Jose, R.; Nazarian, E. J.; Schumann, P.; Van Landschoot, A.; Wirth, S. E.; Musser, K. A., Sporosarcina newyorkensis sp. nov. from clinical specimens and raw cow's milk. Int. J. Syst. Evol. Microbiol. 2012, 62, (2), 322-329. 67. Adebo, O. A.; Njobeh, P. B.; Mavumengwana, V., Degradation and detoxification of AFB1 by Staphylocococcus warneri, Sporosarcina sp. and Lysinibacillus fusiformis. Food Control 2016, 68, 92-96. 68. Nikolausz, M.; Kappelmeyer, U.; Nijenhuis, I.; Ziller, K.; Kastner, M., Molecular characterization of dichloromethane-degrading Hyphomicrobium strains using 16S rDNA and DCM dehalogenase gene sequences. Syst. Appl. Microbiol. 2005, 28, (7), 582-587. 69. Wang, L.; Wen, Y.; Guo, X. Q.; Wang, G. L.; Li, S. P.; Jiang, J. D., Degradation of methamidophos by Hyphomicrobium species MAP-1 and the biochemical degradation pathway. Biodegradation 2010, 21, (4), 513-523. 70. Hayes, A. C.; Liss, S. N.; Allen, D. G., Growth Kinetics of Hyphomicrobium and Thiobacillus spp. in Mixed Cultures Degrading Dimethyl Sulfide and Methanol. Appl. Environ. Microbiol. 2010, 76, (16), 5423-5431. 71. Dejonghe, W.; Boon, N.; Seghers, D.; Top, E. M.; Verstraete, W., Bioaugmentation of soils by increasing microbial richness: missing links. Environ. Microbiol. 2001, 3, (10), 649-657. 72. Lawniczak, L.; Marecik, R.; Chrzanowski, L., Contributions of biosurfactants to natural or induced bioremediation. Applied Microbiology and Biotechnology 2013, 97, (6), 2327-2339. 73. Zhao, Z.; Wong, J. W., Biosurfactants from Acinetobacter calcoaceticus BU03 enhance the solubility and biodegradation of phenanthrene. Environ. Technol. 2009, 30, (3), 291-299.

743

31

ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 36

744

Captions

745

Figure 1 The shift tendency of OTU_56, OTU_65, OTU_111 and OTU_106 fragments. The

746

relative abundances of the OTU_56, OTU_65, OTU_111 and OTU_106 fragments over a range of

747

buoyant density (BD) of DNA extracted from the

748

standard deviation (SD) relative abundance from triplicate microcosms are shown.

749

Figure 2 The shift tendency of OTU_8, OTU_15, OTU_43, OTU_83, OTU_23, OTU_106 and

750

OTU_111 fragments. The relative abundance of the OTU_8, OTU_15, OTU_43, OTU_83,

751

OTU_23, OTU_106 and OTU_111 fragments over a range of buoyant density (BD) of DNA

752

extracted from the 12C_ABA and

753

abundance from triplicate microcosms are shown.

754

Figure 3 Phylogenetic tree of identified OTUs responsible for PHE degradation.

755

Neighbour-joining tree based on 16S rRNA gene sequences showing the phylogenetic position of

756

the bacteria corresponding to OTU_8, OTU_15, OTU_23, OTU_43, OTU_83, OTU_56, OTU_65,

757

OTU_111, OTU_106 and their representatives of other related taxa. Bootstrap values (expressed

758

as percentages of 1000 replications) > 50% are shown at the branch points. Bar 0.05 substitutions

759

per nucleotide position.

760

Figure 4 Phylogenetic tree of amplified PAH-RHDα GP genes from the heavy fractions of

761

13

762

amplified from heavy fractions of 13C_AD treatment, and PAH-RHD L1, L6 and L10 represent the

763

PAH-RHDα genes from the 13C_ABA treatment.

C_AD and

13

13

12

C_AD and

13

C_AD treatments. Mean ±

C_ABA treatments. Mean ± standard deviation (SD) relative

C_ABA microcosms. PAH-RHD C1 and C15 represent the PAH-RHDα genes

32

ACS Paragon Plus Environment

Page 33 of 36

Environmental Science & Technology

764

Figure 1

3

3

Relative Abundance (%)

Rhodoplanes (OTU_56-bacteria)

Xanthomonadaceae (OTU_65-bacteria)

12C 13C

2

2

1

1

0 1.68

1.70

1.72

1.74

1.76

0 1.68

1.78

1.5

1.70

1.72

1.74

1.76

1.78

1.74

1.76

1.78

8

Mycobacterium (OTU_106-bacteria)

6

Enterobacteriaceae (OTU_111-bacteria)

1.0

4

0.5 2

0.0 1.68

1.70

1.72

1.74

1.76

1.78

0 1.68

1.70

1.72

Buoyant Density ( g/mL)

765 766

33

ACS Paragon Plus Environment

Environmental Science & Technology

767

Page 34 of 36

Figure 2 1.6

Relative Abundance (%)

1.2

1.6

Bacillus (OTU_8-bacteria)

0.8

0.8

0.4

0.4

0.0 1.70

1.72

1.74

1.76

1.78

3.0

2.4

Paenibacillus (OTU_43-bacteria)

1.2

0.0 1.68

1.70

1.72

1.74

12C 13C

1.76

1.78 0.5

0.30

Ammoniphilus (OTU_15-bacteria)

Sporosarcina (OTU_83-bacteria)

0.24

0.4

Mycobacterium (OTU_106-bacteria)

0.3 1.8

0.18

1.2

0.12

0.6

0.06

0.2

0.1

0.0 0.0 1.68

1.70

1.72

1.74

1.76

1.78

2.0

1.6

0.00 1.68

1.70

1.72

1.74

1.76

10

Hyphomicrobium (OTU_23-bacteria)

8

1.2

0.8

1.68

1.78

1.70

1.72

1.74

1.76

1.78

1.76

1.78

1.0

Enterobacteriaceae (OTU_111-bacteria)

0.8

6

0.6

4

0.4

2

0.2

Acinetobacter tandoii (OTU_4-bacteria)

0.4

0.0 1.68

768 769

1.70

1.72

1.74

1.76

1.78

0 1.68

1.70

1.72

1.74

1.76

1.78

0.0 1.68

Buoyant Density (g/ mL)

34

ACS Paragon Plus Environment

1.70

1.72

1.74

Page 35 of 36

Environmental Science & Technology

770

Figure 3 Salmonella enterica subsp. houtenae DSM 9221 (U92195) Escherichia coli EP5 (KU904492.1) Escherichia fergusonii ATCC 35469 (CU928158) 72

OTU_111 bacteria Shigella dysenteriae ATCC 13313 (X96966)

62 100

Kluyvera georgiana ATCC 51603 (AF047186) Raoultella planticola ATCC 33531 (JMPP01000074)

88

75

Shigella sonnei GTC 781 (AB273732)

Gammaproteobacteria

Enterobacteriales

57

Xanthomonadales

Rosenbergiella nectarea 8N4 (HQ284827) Luteibacter anthropi CCUG 25036 (FM212561) 100

OTU_65 bacteria Frateuria aurantia DSM 6220 (AGRV01000006)

82 95

Dyella soli JS12-10 (EU604272) OTU_23 bacteria

Alphaproteobacteria

96 90

Hyphomicrobium denitrificans ATCC 51888 (ACVL01000012)

100

Hyphomicrobium chloromethanicum CM2 (AF198623) Hyphomicrobium facile ATCC 27489 (Y14311)

100

98 100

OTU_56 bacteria Uncultured Rhodoplanes sp. clone Leob163 (KF226093.1)

Rhodoplanes elegans AS130 (D25311) 63 72

Rhodoplanes piscinae JA266 (AM712913) Rhodoplanes roseus 941 (D25313)

Actinobacteria

Mycobacterium brisbanense ATCC 49938 (AY012577) Mycobacterium elephantis 484 (AJ010747) 100

Mycobacterium moriokaense DSM 44221 (AJ429044) Mycobacterium pulveris DSM 44222 (AJ429046)

OTU_15 bacteria 99 Ammoniphilus Resinae CC-RT-E (HM193518) 96 Ammoniphilus oxalaticus RAOx1 (Y14578) Ammoniphilus oxalivorans RAOx-FS (Y14580) OTU_43 bacteria 94

Paenibacillus xylanisolvens X11-1 (AB495094) Paenibacillus naphthalenovorans PR-N1 (AF353681)

86

Sporosarcina aquimarina SW28 (AF202056) 84 Sporosarcina saromensis HG645 (AB243859) OTU_83 bacteria Sporosarcina luteola Y1 (AB473560) Uncultured Bacillus sp. Jel2 8B (LT601541.1) Bacillus sp. IS 1a (KU244571.1) Bacillus sp. ADMK38 (KU850994.1) 95

OTU_8 bacteria Bacillus sp. MJAU D0051 (KT350461.1)

0.02

771

35

ACS Paragon Plus Environment

Bacillaceae

99

Planococcaceaee

Paenibacillus soli DCY03 (DQ309072) Sporosarcina ureae DSM 2281 (AF202057)

Firmicutes

59 100

Paenibacillaceae

96

Aneurinibacillus

OTU_106 bacteria

Environmental Science & Technology

772

Figure 4

Mycobacterium vanbaalenii PYR-1 PAH-RHD gene (AAY85176.1) PAH-RHDα L6 Mycobacterium sp. py143 NidA3 gene (ADH94645.1) Mycobacterium novocastrense PAH-RHD gene (GAT12202.1) Terrabacter sp. FLO PAH-RHD gene (ABA87073.1) Marine bacterium UBF-Py PAH-RHD gene (AFT63035.1) Mycobacterium sp. AP1 PAH-RHD gene (AFT63034.1) Uncultured bacterium PAH-RHD gene (AEW70580.1) Uncultured bacterium PAH-RHD gene (AMM73050.1) Uncultured bacterium NidA3 gene (AIC78910.1) Uncultured bacterium PAH-RHD gene (AMM73086.1) 99

Uncultured bacterium PAH-RHD gene (AEW70585.1) Uncultured bacterium PAH-RHD gene (AMM73080.1) Uncultured bacterium PAH-RHD gene (AEW70614.1) 63 PAH-RHDα C15 Uncultured bacterium PAH-RHD gene (AEW70574.1) Uncultured bacterium PAH-RHD gene (AEW70635.1)

97 PAH-RHDα L1 PAH-RHDα C1

76 87

Mycobacterium sp. S23 PAH-RHD gene (ALS30442.1) Uncultured bacterium RHD gene (AGO66328.1)

57

96 Mycobacterium sp. CH-1 PAH-RHD gene (ABD97978.1) 76 Mycobacterium sp. SNP11 PAH-RHD gene (ABK27720.1) Uncultured Comamonas sp. PAH-RHD gene (CAN85239.1) 99

Aeromonas sp. XF3-2 PAH-RHD gene (AND66083.1)

PAH-RHDα L10 70 Bacillus sp. CL1-1 PAH-RHD gene (AND66078.1) 61 Enterobacter sp. CL1-2 PAH-RHD gene (AND66079.1) Uncultured Comamonas sp. PAH-RHD gene (CAN85224.1) Uncultured bacterium NagA gene (AHZ97886.1)

773

0.1

36

ACS Paragon Plus Environment

Page 36 of 36