Jiri Wald1
Blanka Vrchotova
Ondrej Uhlik- 1Department of Biochemistry and Microbiology, Faculty of Food and Biochemical Technology, University of Chemistry and Technology, Prague, Prague, Czech Republic
- 2Department of Genomics and Bioinformatics, Institute of Molecular Genetics, Czech Academy of Sciences, Prague, Czech Republic
- 3Laboratory of Fungal Biology, Institute of Microbiology, Czech Academy of Sciences, Prague, Czech Republic
Given that the degradation of aromatic pollutants in anaerobic environments such as sediment is generally very slow, aeration could be an efficient bioremediation option. Using stable isotope probing (SIP) coupled with pyrosequencing analysis of 16S rRNA genes, we identified naphthalene-utilizing populations in aerated polyaromatic hydrocarbon (PAH)-polluted sediment. The results showed that naphthalene was metabolized at both 10 and 20°C following oxygen delivery, with increased degradation at 20°C as compared to 10°C—a temperature more similar to that found in situ. Naphthalene-derived 13C was primarily assimilated by pseudomonads. Additionally, Stenotrophomonas, Acidovorax, Comamonas, and other minor taxa were determined to incorporate 13C throughout the measured time course. The majority of SIP-detected bacteria were also isolated in pure cultures, which facilitated more reliable identification of naphthalene-utilizing populations as well as proper differentiation between primary consumers and cross-feeders. The pseudomonads acquiring the majority of carbon were identified as Pseudomonas veronii and Pseudomonas gessardii. Stenotrophomonads and Acidovorax defluvii, however, were identified as cross-feeders unable to directly utilize naphthalene as a growth substrate. PAH degradation assays with the isolated bacteria revealed that all pseudomonads as well as Comamonas testosteroni degraded acenaphthene, fluorene, and phenanthrene in addition to naphthalene. Furthermore, P. veronii and C. testosteroni were capable of transforming anthracene, fluoranthene, and pyrene. Screening of isolates for naphthalene dioxygenase genes using a set of in-house designed primers for Gram-negative bacteria revealed the presence of such genes in pseudomonads and C. testosteroni. Overall, our results indicated an apparent dominance of pseudomonads in the sequestration of carbon from naphthalene and potential degradation of other PAHs upon aeration of the sediment at both 20 and 10°C.
Introduction
Polyaromatic hydrocarbons (PAHs) are ubiquitous pollutants of either anthropogenic or natural origin (Shuttleworth and Cerniglia, 1995). PAHs of anthropogenic origin are typically derived from fossil fuel combustion, refining of crude oil, tar production, waste incineration, and food technologies, whereas naturally derived PAHs originate from forest fires, natural oil seeps, and some microbial processes. Several microorganisms, including bacteria, yeasts, and other fungi, have been implicated in the metabolism of PAHs (Cerniglia, 1992; Pothuluri and Cerniglia, 1994; Cerniglia and Sutherland, 2010). Low molecular weight PAHs with up to three rings, such as naphthalene, phenanthrene, and anthracene, can serve as carbon sources for many bacteria whose complete degradation pathways have been well-described (van der Meer et al., 1992). However, it is less common for microorganisms to degrade higher molecular weight PAHs with four or more rings. With some exceptions (Shuttleworth and Cerniglia, 1996; Ho et al., 2000; Kim et al., 2005; Lai et al., 2012), high molecular weight PAHs are cometabolized with the low molecular PAHs, salicylate, and other plant-derived compounds being the primary substrates (van Herwijnen et al., 2003; Rentz et al., 2005, 2008; Peng et al., 2008).
In the initial steps of aerobic degradation, molecular oxygen is required to oxidize and cleave the aromatic ring (Cerniglia, 1992). In the absence of oxygen (anaerobic conditions), the degradation rates of PAHs are much lower (Eriksson et al., 2003). Due to the usually limited amount of air in many environments, oxygen is delivered through aeration when biodegradation needs to be accelerated (Lu et al., 2011). When assessing biodegradation potential, the presence of populations needs to be verified which are able to aerobically degrade, and ideally mineralize, the substrates. Stable isotope probing (SIP) is one of molecular methods for determining biodegradation potential (Chen and Murrell, 2010; Chen et al., 2010; Uhlik et al., 2013a). Using this method, many important findings regarding PAH-degrading organisms have been obtained, such as the discovery of a previously undescribed naphthalene-utilizing bacterium Polaromonas naphthalenivorans (Jeon et al., 2003, 2004), confirmation of salicylate as an inducer of naphthalene catabolism (Singleton et al., 2005), and detection of in situ PAH-degrading populations (Padmanabhan et al., 2003). Recently, SIP has been combined with high-throughput sequencing to gain a more in-depth insight into PAH-utilizing microbial populations (Singleton et al., 2011).
In addition to oxygen, temperature is an important factor in PAH biodegradation. The ambient temperatures typically used for the majority of cultivation-based studies often inadequately represent temperatures found in many environments. Although, various studies show that limited degradation is associated with lower temperatures (Eriksson et al., 2003; Brakstad et al., 2008; Okere et al., 2012; Leewis et al., 2013), psychrotolerant PAH-degrading organisms have also been found (Ahn et al., 2005; Sørensen et al., 2010).
In this study, we investigated microbial populations responsible for the biodegradation of PAHs in contaminated sediment after aeration. We hypothesized that (i) contaminated sediment contains populations with a potential to aerobically degrade naphthalene and other PAHs; (ii) aeration will result in metabolic activity in these populations; and (iii) aerobic biodegradation will occur at 10°C—a temperature more similar to that found in situ. SIP was coupled with the sequencing of 16S rRNA gene pyrotags to identify naphthalene-utilizing populations. The results of the study were confirmed by culture-based experiments and isolation of active naphthalene degraders in pure cultures.
Materials and Methods
Sediment Samples
Sediment samples were collected in November 2011 from a lake in Romania located near an oil refinery, where aeration is used to stimulate intrinsic bioremediation. Sediment dry weight (determined at 105°C) was 16% and pH was 8.0 ± 0.5, with PAH content being listed in Table 1. Immediately after transport to the laboratory, the samples were used for the construction of SIP microcosms and cultivation-based studies.
Table 1. The content of PAHs in the analyzed sediment.
SIP Microcosms
Duplicate SIP microcosms were constructed in 250 mL sterile flasks with 0.5 mg of uniformly 13C10-labeled naphthalene (>99% 13C, Sigma-Aldrich, USA), which was applied as a solution in acetone (concentration 0.05 g/ml). After the acetone was allowed to evaporate leaving behind naphthalene crystals, 2.5 g of settled sediment was added with basal mineral salt solution (Uhlik et al., 2009) to obtain a final volume of 25 mL. The flasks were sealed, incubated at 10 and 20°C for 4, 7, and 14 days with shaking at 130 rpm. Parallel incubation with unlabeled naphthalene (Sigma-Aldrich, USA) was carried out for 14 days. After the incubation time had passed, the microcosms were stored at −20°C prior to downstream processing.
Isotopic Composition of Headspace CO2
Microcosms were thawed at laboratory temperature and then kept in warm water for 20 min to release dissolved gas. Using a Hamilton GASTIGHT syringe, 250 μL of the headspace gas was transferred to a 12 mL W938 borosilicate vacutainer (IVA Analysentechnik e.K., Germany) that had been flush-filled with helium. The isotopic composition of headspace CO2 was analyzed using a Delta V Advantage IRMS (Thermo Fisher Scientific, USA) coupled with a Gasbench II device equipped with a single cryotrap.
Isolation of Labeled DNA
Sediment suspensions from SIP microcosms and from initial time zero (T0) soil samples were pelleted at 5000 × g for 10 min and frozen until further processing. The metagenomic DNA was isolated using PowerMax™ Soil DNA Isolation Kits (MoBio Laboratories Inc., USA) and concentrated by ethanol precipitation with glycogen (Roche, Switzerland). DNA concentration in all samples was normalized to 100 ng/μL by diluting with PCR grade water (Sigma-Aldrich, USA). Subsequently, 8 μL of DNA solution was mixed with cesium trifluoroacetate (Amersham, UK) pre-diluted to a concentration of 1.6 g/mL in a 2 mL ultracentrifugation cuvette. Labeled DNA was isolated by isopycnic centrifugation performed on a Discovery 90 Ultracentrifuge using a TFT-80.2 Fixed-Angle Ultraspeed Centrifuge Rotor (Sorvall, USA) under conditions of 145,000 × g for >70 h. Each gradient was fractionated into 60 μL fractions using a Beckman Fraction Recovery System (Beckman Coulter, USA) and a Harvard Pump 11 Plus Single Syringe (Harvard Apparatus, USA) by replacement with water at a flow rate of 200 μL/min.
Isopycnic centrifugation and gradient fractionation of the samples were run in parallel with two blanks (DNA replaced with water) which were used to ensure reproducibility and accuracy of both centrifugation and fractionation by measuring the buoyant density of each fraction with a Digital Handheld Refractometer (Reichert Analytical Instruments, USA). DNA was retrieved from cesium trifluoroacetate by isopropanol precipitation with glycogen (Uhlik et al., 2009).
The distribution of DNA as a function of buoyant density was determined by qPCR with primers 786f, 5′-GATTAGATACCCTGGTAG-3′ and 939r, 5′-CTTGTGCGGGC CCCCGTCAATTC-3′ targeting the 16S rRNA gene. DNA quantification was performed in 12 μL reactions using a DyNAmo Flash SYBR Green qPCR Kit (Thermo Scientific, USA), 5 pmol of each primer, and 2 μL of DNA. Cycling was performed in a CFX96 Real-Time System (Bio-Rad, USA) following the cycling program of a denaturation at 95°C for 7 min, followed by 35 cycles at 95°C for 20 s, 55°C for 30 s, and 72°C for 30 s, followed by the final extension at 72°C for 2 min. Fractions with 13C-enriched DNA were combined and further analyzed as 13C-DNA. Corresponding fractions from the gradients with unlabeled DNA and T0 DNA were combined and analyzed as controls.
Pyrosequencing
Amplicons of 16S rRNA genes were generated using a two-step PCR approach which was adapted and modified from Berry et al. (2011). The first PCR was performed with primers 515f, 5′-GTGCCAGCMGCNGCGG-3′ and 1406r, 5′-ACGGGCGGTGWGTRC-3′ in 20 μL reactions with 0.2 mM dNTPs (Thermo Scientific, USA), 0.5 μM primers, 0.1 mg/mL of BSA, and 0.4 U of Phusion Hot Start II DNA Polymerase with the corresponding HF buffer (Thermo Scientific, USA). Thermocycling conditions were set as follows: 98°C for 3 min, followed by 35 cycles of 98°C for 10 s, 61°C for 30 s, and 72°C for 30s, with the final elongation at 72°C for 10 min. The resulting PCR products were used as templates for a second PCR run, in which the primers 515f and 1406r were fused with sequencing adapters (454 Sequencing Application Brief No. 001-2009, Roche) at 5′-ends. The forward primers contained different tags (454 Sequencing Technical Bulletin No. 005-2009, Roche) which allowed for multiplexing. PCR was performed in 25 μL reactions with 0.5 μL of the template, 0.2 mM dNTPs (Thermo Scientific, USA), 1 μM primers, 0.1 mg/mL of BSA, and 0.5 U of Phusion Hot Start II DNA Polymerase with the corresponding HF buffer (Thermo Scientific, USA). The cycling program was set to 98°C for 3 min, followed by 10 cycles at 98°C for 10 s, 61°C for 30 s, and 72°C for 30 s, with the final elongation at 72°C for 10 min. Finally, the PCR products were pooled and purified using AMPure XP Beads (Agencourt, Beckman Coulter, USA) to remove residual primer-dimers according to the manufacturer's instructions. Amplicons were unidirectionally sequenced from the forward primer using GS FLX+ chemistry followed by standard analysis of signal processing (Roche).
Sequence Analysis
Pyrosequencing data were processed with the aid of the mothur software package, version 1.31.1 (Schloss et al., 2009). The flowgrams, trimmed to 650–800 flows, with the number of differences in barcode and primer being set to 0, were denoised until the change in flowgram correction reached 10−6. Barcode and primer sequences were removed from the fasta sequences generated, and unique sequences were aligned against merged SILVA bacterial and archaeal reference alignments. Sequences that did not align well or were shorter than 400 bp were removed, and unique sequences were pre-clustered using the pseudo-single linkage algorithm merging sequences with a difference of 1 bp per 100 bp of sequence length. Chimeric sequences were identified using the Perseus algorithm (Quince et al., 2011) and removed from the data set. Classification of valid sequences was performed against Ribosomal Database Project (RDP) (Cole et al., 2009; Cole and Tiedje, 2014) reference files (trainset 9). Singletons and contaminating sequences (mitochondria, chloroplasts, eukaryotic or unknown sequences) were removed from the data set. Final sequences were clustered by average linkage algorithm at 3% to create operational taxonomic units (OTUs). The most abundant sequence in each OTU was selected as a representative sequence and was submitted to RDP Seqmatch (Cole et al., 2009) in order to identify the closest type strains. The pyrosequencing error rate was determined by analyzing mock community sequences as described previously (Uhlik et al., 2013b).
Bacterial Extraction and Cultivation
The sediment samples were cultivated in liquid enrichment cultures (Kurzawova et al., 2012), with naphthalene as the sole carbon source. Sediment suspensions were prepared in liquid mineral salt solution (10%, w/v, Uhlik et al., 2009) with the addition of naphthalene crystals (~1 mg/mL) and transferred at 1 week intervals for 8 weeks at 10 and 20°C. Finally, 100 μL of each suspension were diluted to extinction in 0.85% (w/v) NaCl solution and plated on solid mineral media (Uhlik et al., 2011) with crystalline naphthalene in the lid (~10 mg). After colonies appeared, they were re-inoculated in a liquid mineral salt solution and incubated with naphthalene to confirm the growth on this substrate. At the same time, each of the colonies was re-inoculated on a Plate Count Agar (PCA), and after 48 h of growth at 28°C the cultures were identified using an Autoflex Speed MALDI-TOF/TOF mass spectrometer and MALDI Biotyper 3.1 software (Bruker Daltonik, Germany) equipped with 5627 MSP library (released in December 2013, Bruker Daltonics, Germany) in a manner similar to that described elsewhere (Uhlik et al., 2011). In order to match the sequences from pyrosequencing analysis with isolated cultured representatives, 16S rRNA genes were sequenced for isolates utilizing naphthalene and those that were not identified by MALDI TOF MS. Briefly, the genomic DNA was isolated using thermal lysis, in which a PCA-grown colony was transferred to a microtube, mixed with 50 μL of PCR grade water (Sigma-Aldrich, USA) and incubated at 95°C for 20 min. Afterwards, the cell debris was pelleted, and the supernatant containing the genomic DNA was used for PCR. The 16S rRNA genes were amplified with primers 8f (Suzuki and Giovannoni, 1996) and 1406r (described above) in 20 μL reactions with 0.2 mM dNTPs (Thermo Scientific, USA), 0.2 μM primers (Generi Biotech, Czech Republic), 0.1 mg/mL of BSA, and 0.02 U/μL of iProof™ DNA polymerase with the corresponding buffer (Bio-rad, USA). The cycling program was set to 98°C for 3 min, followed by 25 cycles at 98°C for 10 s, 56°C for 20 s, and 72°C for 40 s, with the final elongation at 72°C for 10 min. The PCR products were re-amplified in a reconditioning step (Thompson et al., 2002), in which 5 μL of PCR product was used as a template in another six-cycle PCR with the total volume of 50 μL with the same primers and other chemicals in the same concentrations as were used for the first PCR. The final PCR products were purified using AMPure XP Beads (Agencourt, Beckman Coulter, USA) following the manufacturer's instructions. Unidirectional sequencing from the forward primer was performed at the Sequencing Laboratory of the Faculty of Science, Charles University in Prague, Czech Republic. The 16S rRNA gene sequences in bacterial isolates and pyrosequencing libraries were compared by local BLAST using a BLAST+ Release 2.2.26 (Zhang et al., 2000).
Dioxygenase Gene Detection
Primers to detect naphthalene dioxygenase large subunit genes were designed on the basis of sequences obtained from a npah gene set in the Functional Gene Pipeline and Repository (Fish et al., 2013). Only sequences longer than 1200 bp and with a cutoff score >900 were used. Primers were designed manually in conserved regions of the gene after aligning the sequences using MEGA5 (Tamura et al., 2011). The sequences of the primers were as follows: 301f, 5′-TGCRRYTAYCAYGGCTGG-3′ and 1099r, 5′-CCATRTTSTCRKTRTCKTC-3′. The region to be amplified begins and ends at nucleotide positions 301 and 1099, respectively, of Pseudomonas putida NCIB 9816-4 nahAc, yielding an expected amplicon size of ~799 bp. The forward primer binding site flanks the Rieske center. In silico testing of the primer set using the Probe Match Search tool in FunGene (Fish et al., 2013) showed that the primers target: Pseudomonas-like nahAc, ndoAc, doxAc, pahAc genes (Kurkela et al., 1988; Denome et al., 1993; Simon et al., 1993; Takizawa et al., 1994; Chauhan et al., 2011); Comamonas-like pahAc genes (Moser and Stahl, 2001); Ralstonia or Polaromonas-like nagAc genes (Fuenmayor et al., 1998; Jeon et al., 2006); and Burkholderia, Acidovorax, or Alcaligenes-like phnAc genes (Laurie and Lloyd-Jones, 1999).
PCR amplification was performed in 20 μL reactions with 0.2 mM dNTPs (Thermo Scientific, USA), 0.8 μM primers, 0.1 mg/mL of BSA, and 0.02 U/μL of Phusion Hot Start II DNA Polymerase with the corresponding HF buffer (Thermo Scientific, USA). The cycling program was set to 98°C for 3 min, followed by 35 cycles at 98°C for 25 s, 56°C for 20 s, and 72°C for 30 s, with the final elongation at 72°C for 10 min. PCR products were processed in a way similar to that described for 16S rRNA genes of the isolated bacteria. The sequences obtained were translated into proteins and aligned with known sequences using MEGA5 (Tamura et al., 2011). The tree was constructed using the Maximum Likelihood method.
Degradation of PAHs
The growth curves of isolated strains were assessed in liquid mineral salt solution, with naphthalene as the sole carbon source at both 10 and 20°C. The cells were pre-grown in liquid mineral salt solution with naphthalene. After a late exponential phase had been achieved, the cultures were filtered through sterile cotton wool in order to remove remaining naphthalene crystals, centrifuged (10 min; 5000 × g) and washed twice with fresh mineral salt solution. Finally, the cells were resuspended in a small amount of mineral salt solution. The degradation assay was performed in triplicate in 8 mL glass vials with screw caps lined with polytetrafluoroethylene. Resuspended cells were diluted with mineral salt solution to the optical density OD600nm of 0.3 and final volume of 4 mL. The PAHs were spiked in the empty vial from PAH stock solutions in hexane which was allowed to evaporate in the fume hood before the cells were added. The final concentration of each PAH (naphthalene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, and pyrene, all from Sigma-Aldrich, USA) in the final mixture was 20 mg/L. After 4, 7, and 14 days of incubation (10 or 20°C; orbitally shaken at 130 rpm) in the dark, remaining PAHs were extracted twice with 1 mL of hexane and 150 mg/L benzene (internal standard) for 2 h (28°C; orbitally shaken at 130 rpm), eventually combining both extracts. Each assay was performed in triplicate. PAHs in the hexane fraction (10 μL injection) were quantified using HPLC (HP1100) analysis with a 150 × 2 mm × 5 μm Luna C18(2) column (Phenomenex) and isocratic elution with 35% acetonitrile:65% water with flow rate 0.8 mL/min. PAHs were detected spectrophotometrically (benzene 202 nm, naphthalene 220 nm, acenaphthene 226 nm with reference 276 nm for fluorene correction, pyrene 240 nm with reference 287 nm for fluoranthene correction, fluorene 263 with reference 308 nm for acenaphthene correction, fluoranthene 286 with reference 345 nm for pyrene correction, phenanthrene and anthracene 250 nm), the retention times were as follows: benzene 4.1 min, naphthalene 12.0 min, fluorene 28.6 min, acenaphthene 26.8 min, phenanthrene 36.1 min, anthracene 42.6 min, fluoranthene 58.6 min, and pyrene 61.7 min. ChemStation software was used for HPLC measurement, integration of the chromatograms as well as quantification of PAHs.
Degradation efficiency was evaluated using the Welch Two Sample t-test in R (R Development Core Team, 2009) at a 5% significance level.
Sequence Deposition
Pyrosequencing reads were deposited in NCBI Sequence Read Archive under accession number SRP063975. 16S rRNA gene and PAH dioxygenase gene sequences were deposited in the GenBank database under accession no. KJ643424–KJ643427 and KJ643430–KJ643433, respectively.
Results
SIP: Headspace 13CO2
The isotopic composition of headspace CO2 was monitored in order to assess the mineralization of 13C-naphthalene. The results showed that mineralization of 13C-naphthalene was faster at 20 than 10°C (Table 2). The amount of evolved 13CO2 did not increase significantly between day 7 and 14 at either temperature and indicated that between 50 and 55% of the original 13C-naphthalene was respired.
Table 2. Headspace CO2 analysis and the amount of 13CO2 evolved.
Quantitative PCR and “Heavy” DNA Isolation
Q-PCR analyses of gradient fractions were performed to examine the distribution of DNA throughout the gradient and to localize 13C-DNA-containing fractions (Figure 1). The dominant peak was recorded at buoyant densities (BDs) of 1.59–1.60 g/mL and represented unlabeled DNA. While a substantial amount of 13C-enriched (or “heavy”) DNA was recorded at BDs >1.62 g/mL, only minor quantities of DNA were detected at the same BDs in control samples (fractionated DNA from parallel incubations with unlabeled naphthalene or from initial (T0) sediment samples). The fractions with BDs >1.62 g/mL were combined and further analyzed as “heavy” DNA along with the corresponding fractions from T0 (control) samples.
Figure 1. Distribution of DNA in the gradient. CTRL and T0 refer to, respectively, DNA isolated from parallel incubations with unlabeled naphthalene and initial sediment sample. Symbols (a) and (b) refer to replicate experimental data.
16S rRNA Gene Profiling
Analysis of pyrosequencing reads of the mock community according to Schloss et al. (2011) showed that the sequence error was 10−3. The total number of valid sequence reads was 80,160 with an average number of 4715 reads per sample/replicate, including controls. Prior to further analyses the sizes of the libraries were normalized according to Schloss et al. (2011) to the number of reads in the least library, which was 1939.
16S rRNA gene profiling indicated that only members of Proteobacteria derived carbon from naphthalene (Table 3). In addition to the clearly predominant pseudomonads at both temperatures, naphthalene-derived 13C was detected in Stenotrophomonas, Acidovorax, and Comamonas, and sporadically in other Proteobacteria. The dominant pseudomonads comprised three OTUs, of which the most abundant OTU clustered with the Pseudomonas fluorescens group (with closest type strains being Pseudomonas veronii CIP 104663, Pseudomonas proteolytica CMS 64, Pseudomonas brenneri CFML 97-391T, and Pseudomonas gessardii CIP 105469) followed by the Pseudomonas syringae group (namely Pseudomonas cichorii LMG 2162T) and Pseudomonas peli R-20805. At 10°C, the number of Pseudomonas sequences in libraries from both 4- and 7-day incubation exceeded 99%, indicating a dominant role for this taxon in the primary consumption of naphthalene. Stenotrophomonas, Acidovorax, and Comamonas were detected only after 14 days of incubation. At 20°C, these sequences appeared earlier over the course of incubation. Other sequences also appeared in the 14-day incubation library at lower abundance rates (Table 3). Most of the sequences detected in SIP were not detected in the initial total community DNA, although more than 99% coverage was achieved according to Schloss et al. (2011). The composition of the initial community DNA is shown in Figure 2.
Table 3. Bacterial populations deriving carbon from naphthalene, either directly or via cross-feeding, as detected at different time points at 10 and 20°C.
Figure 2. Relative abundance of reads affiliated with phyla in the initial total community DNA isolated from the sediment prior to SIP incubations calculated from the total number of 4024 reads.
SIP vs. Cultivation-based Method
Aerobic naphthalene-utilizing bacteria were enriched in cultures, with naphthalene as the sole carbon source. Dilution to extinction and selective plating on mineral media with naphthalene resulted in the isolation of seven bacterial cultures, four of which were proved to utilize naphthalene as a sole carbon and energy source (Table 4). These naphthalene-utilizing cultures were identified as P. veronii, P. gessardii, Pseudomonas sp., and Comamonas testosteroni.
Table 4. Naphthalene-utilizing bacterial isolates extracted from the sediment and their identification based on MALDI-TOF MS with MALDI Biotyper and/or 16S rRNA gene analysis.
Furthermore, the culture-based study permitted a more reliable classification of the pyrosequencing reads to a species level. As shown by local BLAST alignment, all the 16S rRNA genes of the sequenced isolates shared >99% identity with the most abundant reads in 13C-DNA pyrosequencing libraries, with differences detected only in homopolymeric regions located toward the distal ends of the pyrosequencing reads. These were very likely to be technology-based errors as pyrosequencing is mainly affected by errors in homopolymeric regions (Kunin et al., 2010). Such errors tend to accumulate especially toward the distal end of the reads (Schloss, 2013). The local BLAST alignment-based comparison of 16S rRNA gene sequences from the pure cultures and the pyrosequencing libraries indicated that the most abundant OTU affiliated to Pseudomonas (Table 3) was associated with two taxa from the P. fluorescens group, namely P. veronii and P. gessardii (isolates 20a2 and 10a3, respectively; Table 4). Both of these isolates differed from the most abundant OTU's representative sequence only at one position in a homopolymeric region. The differences in the P. veronii and P. gessardii 16S rRNA genes were located outside the pyrosequenced region, resulting in their inclusion in a common OTU. The OTU described by the representative of P. cichorii (Table 3) seemed to be represented rather by P. taiwanensis from the P. putida group (isolate 20b3; Table 4); however, as the findings were not confirmed by mass spectrometry, species-level identification was inconclusive. The OTU identified as Comamonas was confirmed to be associated with C. testosteroni.
Naphthalene Dioxygenases
The detection of naphthalene dioxygenase large subunit genes (nahAc) was performed with newly designed primers reported in the Materials and Methods Section. Naphthalene dioxygenase genes were detected in all strains that utilized naphthalene in the culture-based experiments. Their sequences clustered with sequences previously detected in pseudomonads and comamonads (Figure 3). The genes in isolates 20a2 and 10a3 differed only in one nucleotide, resulting in one amino acid difference (corresponding to position 219 of Pseudomonas putida NCIB 9816-4 NahAc).
Figure 3. Naphthalene dioxygenase large subunit Maximum Likelihood tree. There were a total of 241 positions in the final dataset. Isolates from this study are highlighted with the black circle; public sequences are shown with their GenBank accession numbers.
Degradation Assay
The degradation assay focused on the seven most abundant PAHs in the sediment, namely naphthalene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, and pyrene. The results of the degradation assay (Figure 4) showed that the highest degradation capabilities could be ascribed to P. veronii 20a2 and C. testosteroni 20a3. These strains were able to degrade all the PAHs tested, including anthracene and the four-ring PAHs fluoranthene and pyrene. Degradation efficiency was higher at 20°C where statistically significant depletion of all PAHs was detected after just 4 days. Naphthalene and acenaphthene were efficiently degraded by all strains at both temperatures as early as on day 4.
Figure 4. Residual content of selected PAHs in medium after the degradation experiment with isolated bacteria. The degradation test was performed at 10 and 20°C. Initial concentration of each PAH was 20 mg/L (~100%, day 0), the experiment was terminated on day 4, 7, and 14. The results shown are mean values of three biological replicates, the colors used for different isolates (orange or red) are equal in meaning and indicate statistical significance of the degradation as determined by Welch Two Sample t-test (α = 0.05).
Discussion
Studying microbial diversity in contaminated environments is essential for predicting the bioremediation potential of autochthonous populations. In the context of bioremediation, the flow of carbon from contaminants through the microbial community is an important issue; detecting both primary degraders and cross-feeding populations enables all microorganisms essential to the completion of the degradation process to be identified (Uhlik et al., 2013a). SIP combined with time-course sampling is one such technique which allows for the detection of members of the community through which the labeled carbon flows.
Fate of Naphthalene-derived Carbon
When employing SIP, it is important not to exceed the environmentally relevant concentrations of the labeled substrate (Neufeld et al., 2007); in this study, the amount of 13C-naphthalene used was 0.5 mg per 2.5 g of the sediment, less than five times the naphthalene concentration measured in the sediment (Table 1). Unsurprisingly, the amounts of 13CO2 evolved by 13C-naphthalene respiration in this study indicated a faster metabolism at 20 vs. 10°C (Table 2). The diversity of populations incorporating the label from 13C-naphthalene did not differ significantly at the two temperatures tested, with only the carbon flow at 10°C being slower.
Pseudomonads dominated the acquisition of naphthalene-derived carbon (Table 3). They were the only taxon detected which contained the label on days 4 and 7 at 10°C and were also dominant at 20°C. The detection of pseudomonads as naphthalene-utilizing bacteria was not surprising, as many studies have described members of this genus to be involved in PAH utilization (Kuiper et al., 2001; Samanta et al., 2002; Singleton et al., 2005, 2011; Jones et al., 2011; Lu et al., 2011; Uhlik et al., 2012), some of which were members of the P. fluorescens group (Whitman et al., 1998; di Gennaro et al., 2006; Bučková et al., 2013). Our study reports an overwhelming dominance of one OTU (Table 3) affiliated with the P. fluorescens group during the degradation of naphthalene. The culture based experiments further determined that this OTU spanned two closely related species of pseudomonads, namely P. veronii and P. gessardii (Table 4). The degradation assay highlighted the exceptional degradation capabilities of the P. veronii culture (Figure 4); in addition to the virtual complete degradation of naphthalene, all other PAHs, including those with four-rings, were transformed by this culture at both 10 and 20°C. The naphthalene dioxygenase sequence detected in P. veronii was identical to that of P. putida NCIB 9816-4 (Figure 3). This enzyme has been reported to oxidize phenanthrene, anthracene (Jerina et al., 1976), acenaphthene (Selifonov et al., 1996) and fluorene (Resnick and Gibson, 1996) in addition to naphthalene.
Although the observed degradation capacity of P. veronii might thus have resulted from the activity of this enzyme, the presence of other dioxygenases in the isolate cannot be ruled out. It has previously been reported that Pseudomonas rhodesiae, which is closely related to P. veronii and P. gessardii, contains redundant dioxygenases responsible for PAH degradation (Kahng et al., 2002). Such a scenario is quite plausible given that the naphthalene dioxygenase genes of P. veronii and P. gessardii differed by only one amino acid although the degradation capabilities of the strains differed markedly (Figure 4). However, resolving these issues is beyond the scope of this study and remains to be elucidated in the future. Indeed, interactions among PAHs during their degradation have previously been reported to be highly complex, ranging from cometabolism (Cerniglia, 1992; Shuttleworth and Cerniglia, 1995; Habe and Omori, 2003; Rentz et al., 2005; Peng et al., 2008) via competition (Stringfellow and Aitken, 1995) to inhibition (Bouchez et al., 1995; Shuttleworth and Cerniglia, 1996).
C. testosteroni (Tables 3, 4), the other population deriving carbon directly from naphthalene, has also been previously described as a PAH-degrading species utilizing naphthalene, phenanthrene or anthracene as its carbon sources (Goyal and Zylstra, 1996, 1997). Additionally, the naphthalene dioxygenase large subunit from C. testosteroni clustered with that of previously described comamonads (Figure 3).
Primary Consumers vs. Cross-feeders
Through the combination of SIP and culture-based methods, we were able to observe that some populations directly involved in 13C-naphthalene degradation were likely to have utilized labeled intermediate metabolites of naphthalene rather than naphthalene itself. In fact, Stenotrophomonas sp., Stenotrophomonas maltophilia, and Acidovorax defluvii were among the bacteria isolated from liquid enrichment cultures which were unable to utilize naphthalene as the sole carbon source when their colonies were selected and reinoculated into liquid medium. It has previously been reported that such a secondary screening of colonies is required as some may grow on suspension carryover (Leigh et al., 2002). In addition, screening of Stenotrophomonas and Acidovorax isolates for the presence of naphthalene dioxygenase genes did not prove the presence of the dioxygenase genes studied. Stenotrophomonads are commonly associated with the utilization of aromatics, including those that are chlorinated (Padmanabhan et al., 2003; Gentry et al., 2004; Leigh et al., 2007; Koubek et al., 2013; Uhlik et al., 2013b), which are mineralized via catechol. Thus, cross-feeding by stenotrophomonads on catechol, or some other naphthalene degradation intermediate, appears to be likely. Although the results of SIP compared with those of cultivation experiments clearly demonstrate that SIP is an efficient method for describing microbial food webs involved in the complete degradation of a substrate, it is difficult to distinguish between primary utilizers of the substrate and cross-feeders. Cultivation could be an efficient supplementary method to SIP which aids in distinguishing between primary substrate utilizers and cross-feeding populations.
Conclusions
Although pseudomonads have been previously associated with the degradation of PAHs, our study reported an apparent dominance of pseudomonads in the sequestration of carbon from naphthalene. The culture-based component of this study demonstrated the extensive degradation capacities of P. veronii and C. testosteroni, spanning the seven most abundant PAHs found in the sediment studied. In conclusion, we demonstrated that natural microflora in the sediment had bioremediation potential, and upon delivery of oxygen initially very low abundant indigenous bacteria were capable of metabolizing naphthalene and potentially other PAHs.
Funding
Financial support was provided by grants awarded to OU from the Czech Science Foundation (no. 13-20414P) and the Czech Ministry of Education, Youth and Sports (no. LH14004). JJ also acknowledges the support provided by J. E. Purkyne's fellowship and the long-term development program RVO 61388971.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank Lucie Musilová and Antonio Pio Marcucci for their technical assistance. The authors also thank Dr. Monika Stavelová at AECOM CZ Ltd. for providing the sediment samples, Michal Strejcek for advice regarding bioinformatics, Jakub Rídl for pyrosequencing, and Dr. Nikolaus Gössweiner-Mohr for his comments on the manuscript and finally Dr. Mary-Cathrine Leewis and Michael O'shea for proofreading the text.
References
Ahn, T. S., Lee, G. H., and Song, H. G. (2005). Biodegradation of phenanthrene by psychrotrophic bacteria from Lake Baikal. J. Microbiol. Biotechnol. 15, 1135–1139.
Berry, D., Ben Mahfoudh, K., Wagner, M., and Loy, A. (2011). Barcoded primers used in multiplex amplicon pyrosequencing bias amplification. Appl. Environ. Microbiol. 77, 7846–7849. doi: 10.1128/AEM.05220-11
Bouchez, M., Blanchet, D., and Vandecasteele, J. P. (1995). Degradation of polycyclic aromatic hydrocarbons by pure strains and by defined strain associations: inhibition phenomena and cometabolism. Appl. Microbiol. Biotechnol. 43, 156–164. doi: 10.1007/BF00170638
Brakstad, O. G., Nonstad, I., Faksness, L. G., and Brandvik, P. J. (2008). Responses of microbial communities in Arctic sea ice after contamination by crude petroleum oil. Microb. Ecol. 55, 540–552. doi: 10.1007/s00248-007-9299-x
Bučková, M., Puškarová, A., Chovanová, K., Kraková, L., Ferianc, P., and Pangallo, D. (2013). A simple strategy for investigating the diversity and hydrocarbon degradation abilities of cultivable bacteria from contaminated soil. World J. Microbiol. Biotechnol. 29, 1085–1098. doi: 10.1007/s11274-013-1277-5
Cerniglia, C. (1992). Biodegradation of polycyclic aromatic hydrocarbons. Biodegradation 3, 351–368. doi: 10.1007/BF00129093
Cerniglia, C. E., and Sutherland, J. B. (2010). “Degradation of polycyclic aromatic hydrocarbons by fungi,” in Handbook of Hydrocarbon and Lipid Microbiology, ed K. N. Timmis (Berlin; Heidelberg: Springer), 2079–2110. doi: 10.1007/978-3-540-77587-4_151
Chauhan, A., Layton, A. C., Williams, D. E., Smartt, A. E., Ripp, S., Karpinets, T. V., et al. (2011). Draft genome sequence of the polycyclic aromatic hydrocarbon-degrading, genetically engineered bioluminescent bioreporter Pseudomonas fluorescens HK44. J. Bacteriol. 193, 5009–5010. doi: 10.1128/JB.05530-11
Chen, Y., and Murrell, J. C. (2010). When metagenomics meets stable-isotope probing: progress and perspectives. Trends Microbiol. 18, 157–163. doi: 10.1016/j.tim.2010.02.002
Chen, Y., Vohra, J., and Murrell, J. C. (2010). Applications of DNA-stable isotope probing in bioremediation studies. Methods Mol. Biol. 599, 129–139. doi: 10.1007/978-1-60761-439-5_9
Cole, J. R., and Tiedje, J. M. (2014). History and impact of RDP: a legacy from Carl Woese to microbiology. RNA Biol. 11, 239–243. doi: 10.4161/rna.28306
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–D145. doi: 10.1093/nar/gkn879
Denome, S. A., Stanley, D. C., Olson, E. S., and Young, K. D. (1993). Metabolism of dibenzothiophene and naphthalene in Pseudomonas strains: complete DNA sequence of an upper naphthalene catabolic pathway. J. Bacteriol. 175, 6890–6901.
di Gennaro, P., Conforti, P., Lasagni, M., Bestetti, G., Bernasconi, S., Orsini, F., et al. (2006). Dioxygenation of naphthalene by Pseudomonas fluorescens N3 dioxygenase: optimization of the process parameters. Biotechnol. Bioeng. 93, 511–518. doi: 10.1002/bit.20736
Eriksson, M., Sodersten, E., Yu, Z., Dalhammar, G., and Mohn, W. W. (2003). Degradation of polycyclic aromatic hydrocarbons at low temperature under aerobic and nitrate-reducing conditions in enrichment cultures from northern soils. Appl. Environ. Microbiol. 69, 275–284. doi: 10.1128/AEM.69.1.275-284.2003
Fish, J. A., Chai, B., Wang, Q., Sun, Y., Brown, C. T., Tiedje, J. M., et al. (2013). FunGene: the functional gene pipeline and repository. Front. Microbiol. 4:291. doi: 10.3389/fmicb.2013.00291
Fuenmayor, S. L., Wild, M., Boyes, A. L., and Williams, P. A. (1998). A gene cluster encoding steps in conversion of naphthalene to gentisate in Pseudomonas sp. strain U2. J. Bacteriol. 180, 2522–2530.
Gentry, T. J., Wang, G., Rensing, C., and Pepper, I. L. (2004). Chlorobenzoate-degrading bacteria in similar pristine soils exhibit different community structures and population dynamics in response to anthropogenic 2-, 3-, and 4-chlorobenzoate levels. Microb. Ecol. 48, 90–102. doi: 10.1007/s00248-003-1048-1
Goyal, A. K., and Zylstra, G. J. (1996). Molecular cloning of novel genes for polycyclic aromatic hydrocarbon degradation from Comamonas testosteroni GZ39. Appl. Environ. Microbiol. 62, 230–236.
Goyal, A. K., and Zylstra, G. J. (1997). Genetics of naphthalene and phenanthrene degradation by Comamonas testosteroni. J. Ind. Microbiol. Biotechnol. 19, 401–407. doi: 10.1038/sj.jim.2900476
Habe, H., and Omori, T. (2003). Genetics of polycyclic aromatic hydrocarbon metabolism in diverse aerobic bacteria. Biosci. Biotechnol. Biochem. 67, 225–243. doi: 10.1271/bbb.67.225
Ho, Y., Jackson, M., Yang, Y., Mueller, J. G., and Pritchard, P. H. (2000). Characterization of fluoranthene- and pyrene-degrading bacteria isolated from PAH-contaminated soils and sediments. J. Ind. Microbiol. Biotechnol. 24, 100–112. doi: 10.1038/sj.jim.2900774
Jeon, C. O., Park, M., Ro, H. S., Park, W., and Madsen, E. L. (2006). The naphthalene catabolic (nag) genes of Polaromonas naphthalenivorans CJ2: evolutionary implications for two gene clusters and novel regulatory control. Appl. Environ. Microbiol. 72, 1086–1095. doi: 10.1128/AEM.72.2.1086-1095.2006
Jeon, C. O., Park, W., Ghiorse, W. C., and Madsen, E. L. (2004). Polaromonas naphthalenivorans sp. nov., a naphthalene-degrading bacterium from naphthalene-contaminated sediment. Int. J. Syst. Evol. Microbiol. 54, 93–97. doi: 10.1099/ijs.0.02636-0
Jeon, C. O., Park, W., Padmanabhan, P., Derito, C., Snape, J. R., and Madsen, E. L. (2003). Discovery of a bacterium, with distinctive dioxygenase, that is responsible for in situ biodegradation in contaminated sediment. Proc. Natl. Acad. Sci. U.S.A. 100, 13591–13596. doi: 10.1073/pnas.1735529100
Jerina, D. M., Selander, H., Yagi, H., Wells, M. C., Davey, J. F., Mahadevan, V., et al. (1976). Dihydrodiols from anthracene and phenanthrene. J. Am. Chem. Soc. 98, 5988–5996. doi: 10.1021/ja00435a035
Jones, M. D., Crandell, D. W., Singleton, D. R., and Aitken, M. D. (2011). Stable-isotope probing of the polycyclic aromatic hydrocarbon-degrading bacterial guild in a contaminated soil. Environ. Microbiol. 13, 2623–2632. doi: 10.1111/j.1462-2920.2011.02501.x
Kahng, H. Y., Nam, K., Kukor, J., Yoon, B. J., Lee, D. H., Oh, D. C., et al. (2002). PAH utilization by Pseudomonas rhodesiae KK1 isolated from a former manufactured-gas plant site. Appl. Microbiol. Biotechnol. 60, 475–480. doi: 10.1007/s00253-002-1137-y
Kim, Y.-H., Engesser, K.-H., and Cerniglia, C. (2005). Numerical and genetic analysis of polycyclic aromatic hydrocarbon-degrading mycobacteria. Microb. Ecol. 50, 110–119. doi: 10.1007/s00248-004-0126-3
Koubek, J., Macková, M., Macek, T., and Uhlik, O. (2013). Diversity of chlorobiphenyl-metabolizing bacteria and their biphenyl dioxygenases in contaminated sediment. Chemosphere 93, 1548–1555. doi: 10.1016/j.chemosphere.2013.07.073
Kuiper, I., Bloemberg, G. V., and Lugtenberg, B. J. (2001). Selection of a plant-bacterium pair as a novel tool for rhizostimulation of polycyclic aromatic hydrocarbon-degrading bacteria. Mol. Plant Microbe Interact. 14, 1197–1205. doi: 10.1094/MPMI.2001.14.10.1197
Kunin, V., Engelbrektson, A., Ochman, H., and Hugenholtz, P. (2010). Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. Environ. Microbiol. 12, 118–123. doi: 10.1111/j.1462-2920.2009.02051.x
Kurkela, S., Lehväslaiho, H., Palva, E. T., and Teeri, T. H. (1988). Cloning, nucleotide sequence and characterization of genes encoding naphthalene dioxygenase of Pseudomonas putida strain NCIB9816. Gene 73, 355–362. doi: 10.1016/0378-1119(88)90500-8
Kurzawova, V., Stursa, P., Uhlik, O., Norkova, K., Strohalm, M., Lipov, J., et al. (2012). Plant-microorganism interactions in bioremediation of polychlorinated biphenyl-contaminated soil. N. Biotechnol. 30, 15–22. doi: 10.1016/j.nbt.2012.06.004
Lai, Q., Li, W., Wang, B., Yu, Z., and Shao, Z. (2012). Complete genome sequence of the pyrene-degrading bacterium Cycloclasticus sp. strain P1. J. Bacteriol. 194, 6677. doi: 10.1128/jb.01837-12
Laurie, A. D., and Lloyd-Jones, G. (1999). The phn genes of Burkholderia sp. strain RP007 constitute a divergent gene cluster for polycyclic aromatic hydrocarbon catabolism. J. Bacteriol. 181, 531–540.
Leewis, M. C., Reynolds, C. M., and Leigh, M. B. (2013). Long-term effects of nutrient addition and phytoremediation on diesel and crude oil contaminated soils in subarctic Alaska. Cold Reg. Sci. Technol. 96, 129–138. doi: 10.1016/j.coldregions.2013.08.011
Leigh, M. B., Fletcher, J. S., Fu, X., and Schmitz, F. J. (2002). Root turnover: an important source of microbial substrates in rhizosphere remediation of recalcitrant contaminants. Environ. Sci. Technol. 36, 1579–1583. doi: 10.1021/es015702i
Leigh, M. B., Pellizari, V. H., Uhlik, O., Sutka, R., Rodrigues, J., Ostrom, N. E., et al. (2007). Biphenyl-utilizing bacteria and their functional genes in a pine root zone contaminated with polychlorinated biphenyls (PCBs). ISME J. 1, 134–148. doi: 10.1038/ismej.2007.26
Lu, X.-Y., Zhang, T., and Fang, H.-P. (2011). Bacteria-mediated PAH degradation in soil and sediment. Appl. Microbiol. Biotechnol. 89, 1357–1371. doi: 10.1007/s00253-010-3072-7
Moser, R., and Stahl, U. (2001). Insights into the genetic diversity of initial dioxygenases from PAH-degrading bacteria. Appl. Microbiol. Biotechnol. 55, 609–618. doi: 10.1007/s002530000489
Neufeld, J. D., Wagner, M., and Murrell, J. C. (2007). Who eats what, where and when? Isotope-labelling experiments are coming of age. ISME J. 1, 103–110. doi: 10.1038/ismej.2007.30
Okere, U. V., Cabrerizo, A., Dachs, J., Jones, K. C., and Semple, K. T. (2012). Biodegradation of phenanthrene by indigenous microorganisms in soils from Livingstone Island, Antarctica. FEMS Microbiol. Lett. 329, 69–77. doi: 10.1111/j.1574-6968.2012.02501.x
Padmanabhan, P., Padmanabhan, S., Derito, C., Gray, A., Gannon, D., Snape, J. R., et al. (2003). Respiration of 13C-labeled substrates added to soil in the field and subsequent 16S rRNA gene analysis of 13C-labeled soil DNA. Appl. Environ. Microbiol. 69, 1614–1622. doi: 10.1128/AEM.69.3.1614-1622.2003
Peng, R.-H., Xiong, A.-S., Xue, Y., Fu, X.-Y., Gao, F., Zhao, W., et al. (2008). Microbial biodegradation of polyaromatic hydrocarbons. FEMS Microbiol. Rev. 32, 927–955. doi: 10.1111/j.1574-6976.2008.00127.x
Pothuluri, J. V., and Cerniglia, C. E. (1994). “Microbial metabolism of polycyclic aromatic hydrocarbons,” in Biological Degradation and Bioremediation Toxic Chemicals, ed G.R. Chaudry (London, UK: Chapman and Hall), 92–124.
Quince, C., Lanzen, A., Davenport, R. J., and Turnbaugh, P. J. (2011). Removing noise from pyrosequenced amplicons. BMC Bioinformatics 12:38. doi: 10.1186/1471-2105-12-38
R Development Core Team (2009). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rentz, J. A., Alvarez, P. J. J., and Schnoor, J. L. (2005). Benzo[a]pyrene co-metabolism in the presence of plant root extracts and exudates: implications for phytoremediation. Environ. Pollut. 136, 477–484. doi: 10.1016/j.envpol.2004.12.034
Rentz, J. A., Alvarez, P. J. J., and Schnoor, J. L. (2008). Benzo[a]pyrene degradation by Sphingomonas yanoikuyae JAR02. Environ. Pollut. 151, 669–677. doi: 10.1016/j.envpol.2007.02.018
Resnick, S. M., and Gibson, D. T. (1996). Regio- and stereospecific oxidation of fluorene, dibenzofuran, and dibenzothiophene by naphthalene dioxygenase from Pseudomonas sp. strain NCIB 9816-4. Appl. Environ. Microbiol. 62, 4073–4080.
Samanta, S. K., Singh, O. V., and Jain, R. K. (2002). Polycyclic aromatic hydrocarbons: environmental pollution and bioremediation. Trends Biotechnol. 20, 243–248. doi: 10.1016/S0167-7799(02)01943-1
Schloss, P. D. (2013). Secondary structure improves OTU assignments of 16S rRNA gene sequences. ISME J. 7, 457–460. doi: 10.1038/ismej.2012.102
Schloss, P. D., Gevers, D., and Westcott, S. L. (2011). Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE 6:e27310. doi: 10.1371/journal.pone.0027310
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Selifonov, S. A., Grifoll, M., Eaton, R. W., and Chapman, P. J. (1996). Oxidation of naphthenoaromatic and methyl-substituted aromatic compounds by naphthalene 1,2-dioxygenase. Appl. Environ. Microbiol. 62, 507–514.
Shuttleworth, K. L., and Cerniglia, C. E. (1995). Environmental aspects of PAH biodegradation. Appl. Biochem. Biotechnol. 54, 291–302. doi: 10.1007/BF02787927
Shuttleworth, K. L., and Cerniglia, C. E. (1996). Bacterial degradation of low concentrations of phenanthrene and inhibition by naphthalene. Microb. Ecol. 31, 305–317. doi: 10.1007/BF00171574
Simon, M. J., Osslund, T. D., Saunders, R., Ensley, B. D., Suggs, S., Harcourt, A., et al. (1993). Sequences of genes encoding naphthalene dioxygenase in Pseudomonas putida strains G7 and NCIB 9816-4. Gene 127, 31–37. doi: 10.1016/0378-1119(93)90613-8
Singleton, D. R., Powell, S. N., Sangaiah, R., Gold, A., Ball, L. M., and Aitken, M. D. (2005). Stable-isotope probing of bacteria capable of degrading salicylate, naphthalene, or phenanthrene in a bioreactor treating contaminated soil. Appl. Environ. Microbiol. 71, 1202–1209. doi: 10.1128/AEM.71.3.1202-1209.2005
Singleton, D. R., Richardson, S. D., and Aitken, M. D. (2011). Pyrosequence analysis of bacterial communities in aerobic bioreactors treating polycyclic aromatic hydrocarbon-contaminated soil. Biodegradation 22, 1061–1073. doi: 10.1007/s10532-011-9463-3
Sørensen, S. R., Johnsen, A. R., Jensen, A., and Jacobsen, C. S. (2010). Presence of psychrotolerant phenanthrene-mineralizing bacterial populations in contaminated soils from the Greenland High Arctic. FEMS Microbiol. Lett. 305, 148–154. doi: 10.1111/j.1574-6968.2010.01920.x
Stringfellow, W. T., and Aitken, M. D. (1995). Competitive metabolism of naphthalene, methylnaphthalenes, and fluorene by phenanthrene-degrading pseudomonads. Appl. Environ. Microbiol. 61, 357–362.
Suzuki, M. T., and Giovannoni, S. J. (1996). Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR. Appl. Environ. Microbiol. 62, 625–630.
Takizawa, N., Kaida, N., Torigoe, S., Moritani, T., Sawada, T., Satoh, S., et al. (1994). Identification and characterization of genes encoding polycyclic aromatic hydrocarbon dioxygenase and polycyclic aromatic hydrocarbon dihydrodiol dehydrogenase in Pseudomonas putida OUS82. J. Bacteriol. 176, 2444–2449.
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121
Thompson, J. R., Marcelino, L. A., and Polz, M. F. (2002). Heteroduplexes in mixed-template amplifications: formation, consequence and elimination by ‘reconditioning PCR’. Nucleic Acids Res. 30, 2083–2088. doi: 10.1093/nar/30.9.2083
Uhlik, O., Jecna, K., Mackova, M., Vlcek, C., Hroudova, M., Demnerova, K., et al. (2009). Biphenyl-metabolizing bacteria in the rhizosphere of horseradish and bulk soil contaminated by polychlorinated biphenyls as revealed by stable isotope probing. Appl. Environ. Microbiol. 75, 6471–6477. doi: 10.1128/AEM.00466-09
Uhlik, O., Leewis, M. C., Strejcek, M., Musilova, L., Mackova, M., Leigh, M. B., et al. (2013a). Stable isotope probing in the metagenomics era: a bridge towards improved bioremediation. Biotechnol. Adv. 31, 154–165. doi: 10.1016/j.biotechadv.2012.09.003
Uhlik, O., Musilova, L., Ridl, J., Hroudova, M., Vlcek, C., Koubek, J., et al. (2013b). Plant secondary metabolite-induced shifts in bacterial community structure and degradative ability in contaminated soil. Appl. Microbiol. Biotechnol. 97, 9245–9256. doi: 10.1007/s00253-012-4627-6
Uhlik, O., Strejcek, M., Junkova, P., Sanda, M., Hroudova, M., Vlcek, C., et al. (2011). Matrix-Assisted Laser Desorption Ionization (MALDI)-Time of Flight Mass Spectrometry- and MALDI Biotyper-based identification of cultured biphenyl-metabolizing bacteria from contaminated horseradish rhizosphere soil. Appl. Environ. Microbiol. 77, 6858–6866. doi: 10.1128/AEM.05465-11
Uhlik, O., Wald, J., Strejcek, M., Musilova, L., Ridl, J., Hroudova, M., et al. (2012). Identification of bacteria utilizing biphenyl, benzoate, and naphthalene in long-term contaminated soil. PLoS ONE 7:e40653. doi: 10.1371/journal.pone.0040653
van der Meer, J. R., de Vos, W. M., Harayama, S., and Zehnder, A. J. (1992). Molecular mechanisms of genetic adaptation to xenobiotic compounds. Microbiol. Rev. 56, 677–694.
van Herwijnen, R., Wattiau, P., Bastiaens, L., Daal, L., Jonker, L., Springael, D., et al. (2003). Elucidation of the metabolic pathway of fluorene and cometabolic pathways of phenanthrene, fluoranthene, anthracene and dibenzothiophene by Sphingomonas sp. LB126. Res. Microbiol. 154, 199–206. doi: 10.1016/S0923-2508(03)00039-1
Whitman, B. E., Lueking, D. R., and Mihelcic, J. R. (1998). Naphthalene uptake by a Pseudomonas fluorescens isolate. Can. J. Microbiol. 44, 1086–1093. doi: 10.1139/cjm-44-11-1086
Keywords: naphthalene, polyaromatic hydrocarbons, biodegradation, stable isotope probing, dioxygenase, Pseudomonas veronii, Pseudomonas gessardii, Comamonas testosteroni
Citation: Wald J, Hroudova M, Jansa J, Vrchotova B, Macek T and Uhlik O (2015) Pseudomonads Rule Degradation of Polyaromatic Hydrocarbons in Aerated Sediment. Front. Microbiol. 6:1268. doi: 10.3389/fmicb.2015.01268
Received: 21 March 2015; Accepted: 30 October 2015;
Published: 19 November 2015.
Edited by:
Pankaj Kumar Arora, Yeungnam University, South KoreaReviewed by:
Jun-Jie Zhang, Wuhan Institute of Virology, Chinese Academy of Sciences, ChinaWilliam James Hickey, University of Wisconsin-Madison, USA
Copyright © 2015 Wald, Hroudova, Jansa, Vrchotova, Macek and Uhlik. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ondrej Uhlik, ondrej.uhlik@vscht.cz