Synergistic Processing of Biphenyl and Benzoate: Carbon Flow Through the Bacterial Community in Polychlorinated-Biphenyl-Contaminated Soil.

Aerobic mineralization of PCBs, which are toxic and persistent organic pollutants, involves the upper (biphenyl, BP) and lower (benzoate, BZ) degradation pathways. The activity of different members of the soil microbial community in performing one or both pathways, and their synergistic interactions during PCB biodegradation, are not well understood. This study investigates BP and BZ biodegradation and subsequent carbon flow through the microbial community in PCB-contaminated soil. DNA stable isotope probing (SIP) was used to identify the bacterial guilds involved in utilizing (13)C-biphenyl (unchlorinated analogue of PCBs) and/or (13)C-benzoate (product/intermediate of BP degradation and analogue of chlorobenzoates). By performing SIP with two substrates in parallel, we reveal microbes performing the upper (BP) and/or lower (BZ) degradation pathways, and heterotrophic bacteria involved indirectly in processing carbon derived from these substrates (i.e. through crossfeeding). Substrate mineralization rates and shifts in relative abundance of labeled taxa suggest that BP and BZ biotransformations were performed by microorganisms with different growth strategies: BZ-associated bacteria were fast growing, potentially copiotrophic organisms, while microbes that transform BP were oligotrophic, slower growing, organisms. Our findings provide novel insight into the functional interactions of soil bacteria active in processing biphenyl and related aromatic compounds in soil, revealing how carbon flows through a bacterial community.

Polychlorinated biphenyls (PCBs) are a class of persistent organic pollutants that were widely used in such industrial applications as heat transfer fluids and electrical insulators until their production was banned globally in the 1980s. It is estimated that almost ten million kilograms still persist in the environment1. High-level PCB contamination has been reported at over 111 sites around the state of Alaska (http://dec.alaska.gov/) and at over 8770 sites across the U.S. (http://www.epa.gov/epawaste/hazard/tsd/pcbs/pdf/national.pdf). Ex situ physical-chemical remediation strategies or landfilling are commonly used for PCB-contaminated soils, while bioremediation using microorganisms and plants has potential to be applied either ex situ or in situ. Cleanup of sites through physical removal of soils and transport to a treatment facility can be extremely costly and result in much of the contaminant left behind in association with loose soil. Bioremediation is a potentially less expensive alternative for contaminated soil cleanup, yet remains challenging for several reasons such as differential recalcitrance of different congeners, extended times required to achieve regulatory cleanup limits, and concerns about the potential for degradative intermediates to remain234.

For more than three decades, the aerobic microbial degradation of PCBs, its analogue biphenyl (BP)5 and its degradation intermediate benzoate (BZ) has been the subject of a large body of research6. The aerobic biodegradation of (chloro)biphenyl to (chloro)benzoate and 2-hydroxypenta-2,4-dienoate, termed the upper PCB degradation pathway, is characterized by four separate reactions depicted in Fig. 1 678. (Chloro)BZ may be further metabolized via catechol and the β–ketoadipate pathway and funneled into central metabolism69. (Chloro)BP and (chloro)BZ can also be degraded via cometabolism, a process in which microbes metabolize chlorinated compounds fortuitously using enzymes synthesized to degrade the primary substrates, in this case BP or BZ101112. Some microorganisms, such as certain Arthrobacter, Pseudomonas, Rhodococcus and Burkholderia, are capable of performing complete mineralization of BP along with some (chloro)BP11213 and growing on a range of chlorinated substrates such as chlorinated BZ and BP51314. Other organisms are known to possess only the upper pathway and/or to generate dead-end metabolites of BP615. BZ catabolic pathways such as the β-ketoadipate pathway are relatively widespread in the soil environment9. In the context of the contaminated soil environment, it is not known if BP is mineralized largely by the same group of organisms possessing both the upper and lower pathways, or if it is processed synergistically by separate microbial guilds specializing in BP and BZ biodegradation. Knowing this could reveal novel insights into community interactions during processing of PCBs and other aromatics (such as plant-associated compounds) and aid in efforts to stimulate relevant organisms during bioremediation1617 or phytoremediation, since plant-associated compounds, such as phenols including flavonoids, have been shown to induce bioremediation of contaminated soils and microbial processing of PCBs71819.

Figure 1

The upper (chloro)biphenyl degradation pathway and lower transformation of (chloro)benzoic acid to catechol.

The enzymes involved in each reaction step are biphenyl 2,3 dioxygenase (i), 2,3-dihydro-2,3-dihydroxybiphenyl dehydrogenase (ii), 2,3-dihydroxybiphenyl 1,2-dioxygenase (iii), 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate hydrolase (iv), benzoate 1,2-dioxygenase (v), and cis-diol dehydrogenase (vi).

Microbial ecologists have long sought to understand the mechanisms underlying microbial community processing of complex carbon (C) substrates. Stable isotope probing (SIP) methods202122 have greatly advanced understanding of which organisms within a complex community are active in degrading certain substrates, whether directly or indirectly. To date, SIP studies have provided valuable insight into the processing of C from individual substrates or mixtures of compounds23242526. SIP also has great potential to reveal how members of the microbial community interact to degrade complex substrates and to show how C from a complex substrate flows through different microbial guilds over time2728.

Because bacteria mediate biodegradative processes in contaminated environments, it is valuable to identify the microorganisms involved and elucidate their interactions and ecological growth strategies. This study aims to reveal information concerning synergistic biodegradation and C flow in weathered PCB-contaminated soil from Kodiak, Alaska, USA. We employed DNA-SIP in parallel incubations with two different substrates to identify the microbes performing the upper (BP) and lower (BZ) PCB biodegradation pathways. We then used bacterial 16S rRNA gene sequencing to identify the individual bacterial taxa that derived C from the substrates over time, including those that actively assimilated C from both substrates. We hypothesized that since not all bacteria possess genes for the upper and lower pathways, (i) a portion of the bacterial community will be active only in either the upper or lower pathway of BP degradation, however (ii) a subset of the community will be able to complete the entire degradation pathway, and (iii) given the widespread nature of BZ degradation genes and potential for faster gain of C and energy for the organism, a larger proportion of the community will be involved in BZ degradation than BP degradation. By using non-halogenated compounds that are standard proxies for PCBs and chlorobenzoates12172930, our study is relevant not only to contaminant biodegradation but also to processing of naturally occurring aromatics that are important to soil carbon cycling.

Results

Within 1 day of incubation, excess 13CO2 was detected in microcosm headspaces containing labeled substrates and increased in abundance throughout the experiment (14 days). No change in 13CO2 was detected in control microcosms containing unlabeled BP or BZ. The amount of 13C-labeled substrate mineralized by the microbial community was calculated based on the amount of 13CO2 evolved (Table 1)31. In the 13C-BP incubated soils, approximately 1.81 μg (0.18% of total BP added) of 13C was mineralized within one day, and up to 5.91 μg (0.59%) within 14 days. In 13C-BZ-incubated microcosms, the amount of 13C mineralized was approximately 5.43 μg (0.54%) within one day and reached 127.8 μg (12.78%; 20-fold higher levels than for 13BP) within a 14-day incubation.

Table 1

Mineralization of 13C-BP and 13C -BZ as determined by 13CO2 evolution.

SIP incubation time (days)	13CO2 excess in headspace (nmol)		13C substrate mineralized (μg)		13C substrate mineralized (%)
	Benzoate	Biphenyl	Benzoate	Biphenyl	Benzoate	Biphenyl
1	311.3	141.1	5.431	1.814	0.543	0.181
4	4235	203.8	73.88	2.618	7.388	0.262
14	7327	458.8	127.8	5.895	12.78	0.590

Q-PCR analyses of gradient fractions were conducted to examine the relative abundance of 16S rRNA gene throughout the gradient and to identify fractions containing 13C-labeled DNA (Supplementary Fig. S1). There was a substantial peak of DNA at buoyant density (BD) of 1.576 g.ml−1, which corresponds to the BD of unlabeled (predominantly 12C) DNA. In 13C-incubated samples, a second peak at BD greater than 1.599 g/ml was detected indicating the presence of 13C -DNA. Based on the q-PCR results, we combined the fractions with BD of 1.584–1.628 g.ml−1 to create a compiled ‘heavy’ fraction for each sample prior to downstream molecular analyses as density gradients run in duplicate for each sample were similar in DNA quantity and distribution (Supplementary Fig. S1). Equivalent fractions were compiled from unlabeled controls and were analyzed identically to detect any unlabeled bacterial DNA contamination in heavy fractions. Mann-Whitney U tests on T-RFLP data indicated that the replicate samples were not significantly different, allowing for the combination of replicates for down-stream pyrosequencing analysis. Only a small subset of peaks in the T-RFLP profiles were found in the heavy fractions from control DNA, and fragments of the same size were not major contributors to T-RFLPs in the heavy fractions from 13C-incubated samples. To identify the bacterial populations in the total community and active in assimilating C from labeled substrates, we sequenced 16S rRNA gene amplicons. The biphenyl early (day 04) sample did not achieve sufficient labeling to distinguish the community from unlabeled controls according to ANOVA and post-hoc Tukey comparison of means (Supplementary Table S1) and was excluded from further analyses. The numbers of pyrosequencing reads of labeled DNA in the samples selected ranged from 5684 to 14848 per sample after sequence processing (Table 2). Despite an intense sequencing effort, the total diversity was not fully sampled as was indicated by rarefaction curves (Supplementary Fig. S2).

Table 2

Diversity estimates and number of sequences obtained after raw sequence data processing and normalization.

Sample ID		Sequences after processing	# OTUs (0.03)	Chao	1/D (Inverse Simpson)	Effective Number of Species
Biphenyl	Day 14	5684	938	1348	44.5	268.6
Benzoate	Day 01	13474	310	553	16.6	54.1
	Day 04	14848	515	856	4.6	10.9
Total Community		9745	2174	4876	159.5	717.3

The abundance of biphenyl dioxygenase genes (bphA), the gene that initiates the upper pathway of biphenyl degradation (Fig. 1), was quantified in total soil DNA and in heavy (labeled) fractions by q-PCR. The relative abundance of BphA gene, normalized by the total 13C-DNA concentration in each sample, was lowest in the total, unlabeled, soil community DNA (1.66 × 10−3 gene copies per ng DNA−1). Bacterial bphA abundance in heavy fractions obtained following 13C-BZ SIP was initially ten-fold higher than in the total soil DNA, and decreased in quantity over the time course (d01 BZ: 1.48 × 10−2 gene copies per ng 13C-DNA−1; d04 BZ: 8.15 × 10−3 gene copies per ng 13C-DNA−1). In 13C-BP labeled samples, bphA abundance was 5.64 × 10−3 gene copies per ng 13C-DNA−1 after 14 days, which was approximately five-fold higher than in total soil DNA and corresponded to d04 BZ bphA quantity.

Phylogenetic affiliations of pyrosequencing reads from 13C-DNA revealed diverse bacterial populations that were involved either directly or indirectly in acquiring C derived from BP or BZ (Supplementary Table S2). Actinobacteria predominated in the total soil bacterial community and in libraries obtained from incubations containing 13C-labeled BZ on day 4, and BP on day 14. Proteobacteria were present in all of the examined time points, but predominated in the BZ day 1 sample, where they accounted for 36.5% of total pyrosequencing reads. In addition, the phyla Acidobacteria, Actinobacteria, Gemmatimonadetes, and Nitrospira were found in all labeled DNA samples. In the early (day 1) 13C-BZ incubations, sequences of the order Rhizobiales predominated the labeled DNA, whereas Burkholderia and Actinomycetales dominated after 4 days of 13C-BZ. After 14 days, labeled sequences from 13C-BP samples were dominated by Actinobacteria. Some pyrosequencing reads could not be assigned by the RDP classifier, with 16.6% of the total community, 5.1-6.0% of the 13C-BZ, and 4.79% of 13C-BP remaining unclassified.

Bacteria labeled during incubation with 13C-BZ increased in richness over time (#OTUs and Chao1), but decreased in diversity index values over time, and bacterial sequences associated with the four-day incubation had the lowest diversity of all measured populations (Table 2). Species evenness shifted across the labeled populations over time, causing the reduction in diversity indices, which take into consideration both richness and evenness, but caused no reduction in richness (i.e. number of OTUs). Comparatively, the diversity of bacteria that derived C from 13C-BP, as measured by Chao1, Inverse Simpson and effective number of species, was higher than that measured in either of the 13C-BZ- incubated samples.

To investigate the overlap in bacterial taxa that utilized the two substrates, the percentage of unique and shared OTUs (defined at 97% sequence similarity) that derived C from the different substrates and time points was analyzed and is presented using Venn diagrams (Fig. 2). When all OTUs present in 13C-DNA from BP and BZ were considered, 22.2% of the OTUs were shared in labeled DNA from both substrates. In 13C-BZ incubated soils, 15.5% of sequences were shared among early (day 1) and late (day 4) time points, and the early harvested soils (day 1) had fewer (27.9%) unique OTUs than the late harvested (day 4) soils with 56.6% unique OTUs. To examine the potential transition between the upper and lower portions of the BP degradation pathway, we compared percentages of OTUs present in both the late (day 14) 13C-BP and early 13C-BZ (day 1) incubated soils: 12.6% of OTUs were shared between the two groups.

Figure 2

Overlap of OTUs across substrates.

Venn diagrams showing the overall overlap of OTUs that derived C from the two labeled substrates at different time points, and the overlap in total labeled OTUs from each substrate. The percentages represent the percent of OTUs unique to each of the individual time points and/or substrate. Sequences were grouped at 3% dissimilarity.

The relative abundance and overlap of OTUs that represent more than 0.2% of the total population in each sample were analyzed and visualized using a heat map and dendrogram (Fig. 3). The community that was labeled at the later 13C-BZ time point was the most unique, while that from the early incubated (day 1) 13C-BZ was more similar to both the late 13C-BP and total communities.

Figure 3

Clustering and heat map of OTUs.

Heat map and dendrogram generated using the relative abundance and overlap of OTUs which represent more than 0.2% of the total population in each sample. The color code indicates relative abundance, ranging from light yellow (low abundance) to dark red (high abundance). The dendrogram was generated using average linkage hierarchical clustering based on a Bray-Curtis distance measure.

To investigate the phylogenetic differences and overlaps between the most dominant utilizers of BZ and BP, we generated a table for the top ten most abundant OTUs from each incubation and time point (Table 3). In soils incubated with 13C-BP after 14 days of incubation the dominant taxa were most similar to Arthrobacter oryzae (OTU 2, 12%), Thermosipho geolei (OTU 7, 5%), Pseudomonas lutea (OTU 24, 4%), Acidothermus cellulolyticus (OTU 5, 3%), and Rhodanobacter soli (OTU 17, 2%). In soils incubated with 13C-BZ, a smaller set of OTUs became highly abundant over the course of the incubation than was observed in the 13C-BP sample. By day one (early), dominant labeled OTUs in the 13C-BZ incubation were Methylobacterium jeotgali (OTU 4), which comprised 19% of the total labeled sequences, and Rhodococcus qingshengii (OTU 8), which represented 11% of sequences. The late (four days) incubated 13C-BZ soils were highly dominated by sequences most similar to Burkholderia phytofirmans (OTU 1, 34%), Burkholderia sediminicola (OTU 3, 27%), and Cupriavidus basilensis (OTU 6, 16%).

Table 3

Top OTUs detected in 13C-DNA after incubation of soil with 13C-BP and 13C-BZ.

The percentage and number (in parenthesis) of sequences associated with each OTU is indicated under each investigated time point and substrate. Taxonomic assignments were generated using the Ribosomal Database Project’s classifier.

aThe numbers of pyrosequencing reads assigned using the RDP classifier (80% confidence threshold) are shown in parentheses.

bSequences represented <0.5% of total sequences present in the group.

cNearest RDP match; similarity score shown in parentheses.

Discussion

With this study, we sought to better understand how C from aromatic compounds flows through the soil microbial ecosystem. Using parallel SIP incubation studies, the experiment revealed taxa that participate in the assimilation of C from the polyaromatic compound BP (an unchlorinated analogue of PCBs) and BZ (product of partial BP degradation, analogue for chlorobenzoates and common monoaromatic in the environment) in a PCB-contaminated soil, and how they respond and work in synergy to process these substrates. An array of bacterial taxa derived C from biodegradation of 13C-BP and/or 13C-BZ, either directly or indirectly through cross feeding. Copy numbers of the BphA gene were enriched in labeled DNA from both BZ and BP incubations beyond that of total soil community DNA, supporting the conclusion that labeled DNA represented organisms truly active in aromatic biodegradation rather than exclusively contaminants or cross-feeders. The majority of taxa that assimilated the 13C-labeled aromatics appeared to perform either BP or BZ utilization, but not both. A subset of organisms derived C from both substrates, which could be explained by three possible scenarios: i) they possess both BP and BZ pathways, and thus can mineralize both substrates, ii) they are BZ utilizers that scavenge BZ released by BP-utilizers or iii) they scavenge other substrates or dead biomass generated by BP- and BZ-utilizing microbes. Through this time-course study, we observed a difference in ecological growth strategy for the two microbial guilds, with fast growing bacteria performing BZ degradation, and slower growing microbes degrading BP.

The rapid response to labeled BZ by the microbial community indicates that microbes involved in the lower pathway of BP degradation (BZ degradation) may be adapted to quickly react and utilize 13C-BZ when it is present in the environment, indicative of a copiotrophic population. Of the bacteria that assimilated BZ-carbon, all but one (early OTU 8; late OTU 2) of the top five most abundant OTUs from each time point have been described as copiotrophic microbes3233. The initially lower rates of 13CO2 evolution in 13C-BP-incubated microcosms might indicate that the BP-utilizing populations could be classified as “oligotrophic” and adapted for the optimal utilization of recalcitrant C, such as BP. Of the five most abundant OTUs from both BP early and late incubations, half have been previously described as oligotrophic populations (Table 3; Fierer et al. 2007). Also, approximately 49% of the BP incubated populations were Gram-positive, or mostly associated with oligotrophic microorganisms (Table 3, Supplementary Table S2).

Members of the genera Pseudomonas and Arthrobacter, implicated in BP metabolism in this study, have often been reported to be associated with the upper pathway of BP degradation, and less frequently the lower pathway, and are known to contain the bph operon and degrade PCBs1612172334. Pseudomonads have previously been reported to be the most prevalent group of bacteria in the biodegradation of complex organic compounds35. Both the Acidimicrobiales and Marmoricola have not been previously identified in the metabolism of BP or PCBs. The closest RDP match for the unclassified Acidimicrobiales was Terriglobus saanensis (Table 3), a bacterium isolated from tundra soils36. Although the ecological roles of many acidobacteria are poorly understood37, the Terriglobus species has wide potential physiological capabilities in soils including hydrolysis of different polysaccharides such as starch, pectin, laminarin and aesculin. Among organisms that derived C from BZ, the genera Rhodococcus, Burkholderia, and Cupriavidus have all previously been implicated in PCB degradation. Additionally, the genus Propionibacterium, which derived C from BZ, has previously been reported to use various aromatic hydrocarbons such as benzene, biphenyl, toluene, and styrene as sole carbon and energy sources38.

We investigated the diversity and overlap of organisms that derived C from each substrate over time to elucidate the flow of C through the community and to determine how the guilds may interact (Figs 2 and 3). While a relatively small number of labeled OTUs increased dramatically in relative abundance in 13C-BZ labeled populations during the incubation, increasing numbers of unique, but low-abundance, OTUs over time were also detected, suggesting that some members of the community were slower-growing organisms, or that C flowed out into the community via secondary feeding.

Most taxa that became labeled were involved in either BP or BZ degradation but not both, based on 16S rRNA analyses. Uhlik et al. (2012) found that, in legacy contaminated soil from the former Czechoslovakia, the majority of taxa involved in BP and BZ degradation overlapped. When all of the OTUs detected in this study in 13C-BP and 13C-BZ incubated soils were compared, 16% of OTUs were shared between the two substrates (Fig. 2), suggesting that the majority of the biodegradative microbial populations do not utilize both the upper and lower BP degradation pathways; rather the pathways occur primarily in separate populations. However, BphA genes were enriched in labeled DNA from both BP and BZ, further supporting the observation that BP-degradation capabilities do exist within a portion of the BZ-utilizing guild. Thirteen percent of the OTUs were shared between late 13C-BP and early 13C-BZ, with more unique OTUs present in late 13C-BP soils (72.1%; Fig. 2). At most, 12.6% of BP-degrading taxa may possess both pathways, but the number is likely lower given that cross-feeding on metabolites or dead biomass of BP-degraders may account for some portion of the overlap (Fig. 2; DeRito et al. 2005). The overlapping OTUs were mostly associated with the phylum Proteobacteria, which has been implicated in both BP and BZ degradation, but also contains known heterotrophs and potential cross feeders2326. Taxa similar to the OTU in the 13C-BP labeled community, Methylobacterium jeotgali (OTU 4), have previously been found as sequencing contamination, however because this OTU represents 19% of the labeled community in BP-incubated samples, and was not found in the other 13C-incubated samples, it does not appear to be a contaminant in this study39. Additionally, methane- and ammonia-oxidizing microbial communities, such as Methylobacterium jeotgali, were also previously hypothesized to co-metabolize PCBs in soils40.

The bacterial guilds that assimilated C from each of the substrates showed very different rates of substrate utilization and population shifts over the course of the incubations. Soil consortia that assimilated each of the two substrates differed in their patterns of substrate mineralization (13CO2 evolution); in which those mineralizing BP did so at a much slower rate than those mineralizing BZ, based on of 13CO2 evolution (Table 1). The BZ utilizing bacterial community became highly dominated by relatively few (1–2) OTUs throughout the course of the incubations, and BP-assimilating bacteria were relatively much more diverse and OTU-rich (Table 2). Although the observed differences in 13C processing among the guilds could be due to differences in bioavailability of BP versus BZ, the observed differences also suggest that communities possess diverse ecological patterns of growth.

Recent work recognizing broad, ecologically-based, classes of bacteria324142 could explain the differences in substrate processing rates and the diversity of organisms that dominate the assimilation of the more recalcitrant (BP) versus the more labile (BZ) substrates. These ecologically defined classes are similar to the r-K scheme which, in general, states that r-strategists are adapted to high rates of reproduction and K-strategists are adapted for optimal utilization of environmental resources43. When applied to microbial ecology, r-strategists are defined as ‘copiotrophic’ bacteria or organisms which show an increase in numbers in response to the addition of labile C sources to soil and K-strategists are defined as ‘oligotrophic’ bacteria which do not show a marked increase in numbers in response to the addition of a labile C source to the soil. It has been found through experimentation, observation, and meta-analytical approaches that Gram-negative bacteria, such as members of the phylum Bacteroidetes and subphyla Beta- and Alpha-proteobacteria, which rapidly colonize and grow under labile C-rich conditions, are copiotrophs. Gram-positive soil bacteria, such as some Actinobacteria, as well as Acidobacteria tend to be more successful in resource-limited situations and are considered oligotrophs3244.

The substrate utilization data, as measured by 13CO2 evolution (Table 1), provide one insight to the ecological growth strategies of the microorganisms associated with each of the substrates. The microcosms incubated with BZ showed a faster rate of 13C–labeled substrate mineralization than the BP-incubated microcosms (Table 1). This difference could be explained by the structural differences between BP and BZ, with BP having a more complex structure. The fact that the BP treatments contained approximately twice the molar quantity of 13C atoms underscores the notably higher rate of BZ substrate mineralization rate compared to BP. It has been found that compounds that require a large number of enzymatic steps to be degraded and are biochemically expensive to degrade, such as BP, may result in lowered substrate use efficiency45. For example, in some soil systems, a 50% loss of PCBs was reported after a five month incubation, whereas a similar loss was reported within 50 days for some chlorobenzoate enrichment cultures4647.

Another indication of differential growth strategies is evidenced in the shifts in populations associated with the two incubated substrates (Table 2). During the rapid utilization of BZ, a smaller number of BZ utilizers appeared to rapidly grow and dominate substrate assimilation over time than was observed during BP utilization. This is evident in the effective number of species in 13C-BZ incubated microcosms, which decreased over time (d01 to d04). Although Chao1 increased over time, diversity decreased when evenness was taken into account (effective number of species, Simpson Index), and the species that were detected increased in relative abundance. In contrast, the bacteria that assimilated BP-derived carbon were considerably more OTU-rich and mineralized the substrate at a much lower rate than for BZ degraders. These observations suggest that BP-metabolizing microorganisms are slower-growing populations than BZ utilizers, and slowly process recalcitrant substrates, both of which are characteristics of oligotrophic organisms.

Using a parallel substrate SIP study design, we revealed that the biotransformation of the aromatics BP and BZ in a PCB-contaminated system was mainly performed by distinct bacterial guilds with different ecological growth strategies. Our results suggest that the microbial populations associated with BZ are mostly fast-growing, copiotrophic organisms, while microbes that transform BP are oligotrophic, slower-growing, organisms. The different patterns of use of the available C sources could represent adaptations for the efficient use of monoaromatic versus polyaromatic substrates. A small but substantial proportion of the community assimilated C from both substrates, which may indicate that they possess both the upper and lower BP pathways, or that they are effective cross-feeders and scavengers of labeled metabolic products. Future studies should focus on the 13C-labeled metabolites produced by BP- and BZ-utilizing populations over time to examine the transfer of metabolites between microbial guilds during degradation of biphenyl-containing compounds. Because BP and BZ are analogues for PCBs and chlorobenzoates, respectively, these findings have implications for PCB biodegradation. The findings imply that guilds responsible for initial biodegradation of PCBs to chlorobenzoates are more of a limiting factor than chlorobenzoate degraders in terms of responsiveness to contaminants and biodegradation activity, and that biostimulation efforts targeting initial PCB-degraders may have the most benefit. Further studies should investigate synergistic biodegradation of these compounds with the added challenges of chloro substitutions.

Materials and Methods

SIP Microcosms

PCB-contaminated soils for this experiment were collected in clean sterile glass jars from the Coast Guard Station at Drury Gulch, Kodiak Island, Alaska, USA as a portion of cleanup actions in accordance with Federal Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA, http://dec.alaska.gov/spar/csp/ dod_sites.htm), and stored at 4 °C until analysis. The calculated total PCB concentration in investigated soils was 32.53 mg kg−1 as determined using Soxhlet extraction and run on an Agilent (Santa Clara, CA, USA) 6890N GC equipped with an ECD detector according to the method reported by Slater et al. (2011).

Microcosms for SIP incubations were established in sterile 125 mL glass serum bottles (Wheaton, Millville, NJ, USA). Uniformly labeled 13C-BP (13C12H10), 13C-BZ (13C7H5O2; both 99% labeled; Isotec, Miamisburg, OH, USA), unlabeled BP, or unlabeled BZ were dissolved in acetone and an appropriate volume was added to each microcosm by pipetting the solution onto the inside of the serum bottle and allowing the acetone to evaporate, leaving behind 1 mg BP or BZ crystals on the walls of the serum bottles23. After the BP or BZ solutions had dried fully, 7 g of homogenized contaminated soil and 400 μL sterile water was added to each serum bottle. Bottles were then sealed with a Teflon stopper and aluminum crimp top. Duplicate microcosms were set up for each substrate and each incubation time point (1, 4 and 14 days). These time points were selected based on times previously shown to reveal BP-degrading bacteria using SIP2326. Microcosms were incubated at room temperature (22 °C) until destructive harvesting and soils were stored at −80 °C until further analysis4849. An aliquot of soil without the addition of 13C labeled substrate was frozen at −20 °C to serve as the time 0 (T0) sample.

Headspace Gas Isotope Analysis

One milliliter of headspace gas was collected at microcosm initiation and at each sampling time point using an N2-purged syringe to quantify mineralization of labeled substrates. CO2 was analyzed for δ13C within seven days on a Trace Gas system interfaced to an IsoPrime mass spectrometer (GV instruments, Manchester, UK) at the Alaska Stable Isotope Facility (University of Alaska Fairbanks) as previously described2324. The mineralization of 13C-labeled BP or BZ to 13CO2 was calculated based on the assumption that in the aerobic incubation, oxygen was the terminal electron acceptor for microbial degradation of the supplied substrate.

13C-DNA Isolation

Total soil DNA was extracted using the PowerMax soil DNA isolation kit (MoBio Laboratories, Inc., USA) using the standard protocol. After the final elution, DNA was concentrated according to Uhlik et al. (2012)26 to increase yield. All solutions were diluted to obtain a concentration of 100 ng/μl. Isopycnic centrifugation and gradient fractionation was performed according to the method of Uhlik et al. (2012)26. Fractions were frozen at −20 °C until further analysis.

To locate and quantify bacterial DNA in SIP fractions, q-PCR targeting bacterial 16S rRNA was performed on every fraction for each time point in duplicate as described previously2450. Fractions that were determined to contain 13C-labeled bacterial DNA based on quantities of DNA in high density fractions that exceeded that of unlabeled controls were combined to create a combined “heavy DNA” sample for each gradient (buoyant density (BD) 1.584–1.628 g.ml−1). Equivalent heavy DNA samples were prepared for controls incubated with unlabeled substrates and they were subjected to the same downstream analyses as 13C-labeled DNA to identify and eliminate any background contamination present in controls.

T-RFLP Analysis

Bacterial community profiling using T-RFLP was performed on heavy fractions (13C-DNA), control heavy fractions (12C incubated soils), and total community DNA to provide an initial screening of the diversity of labeled DNA, to determine if replicates were statistically different prior to down-stream sequencing, and to assess any background contamination in heavy fractions from unlabeled DNA prior to pyrosequencing. PCR reactions to amplify the 16S rRNA gene were carried out as described previously5152. Terminal restriction fragments less than 50 bp and above 1000 bp in size were eliminated from all data sets.

Phylogenetic identification

Amplification and pyrosequencing were performed by a modification of a method described previously by Engelbrektson et al. (2010)26. A portion of the 16S rRNA gene (positions 8 to 357 [V1 and V2], Escherichia coli numbering) was amplified using the 8F primer containing a Roche 454 A pyrosequencing adaptor, a 10 base pair unique bar code and 357R primer also containing a Roche 454 B pyrosequencing adaptor.

The PCR mixture (final volume, 25 μl) contained 1 μl each primer (10 μM), 0.5 μl dNTP mix (10 mM), 2.5 μl FastStart 10 X Buffer #2, 0.25 μl FastStart HiFi Polymerase (5 U/μl), and 18.75 μl molecular biology grade water (Roche). The following thermal cycling scheme was used: initial denaturation at 95 °C for 3 min and 30 cycles of denaturation at 95 °C for 30 sec, annealing for 1 min at 55 °C, and extension at 72 °C for 1.5 min, followed by a final extension period at 72 °C for 10 min. Each PCR product was obtained in three parallel reactions, the resulting preparations were mixed, purified using Pure-Link PCR purification kit (Invitrogen, USA) and pooled for downstream sequencing. Roche 454 GS FLX Titanium sequencing (454 Life Sciences, USA) was performed on pooled reactions at the Institute of Molecular Genetics, Czech Academy of Sciences, Prague, Czech Republic.

Processing of Pyrosequencing Data

Raw pyrosequencing data (*.sff files) were processed using the mothur software package version 1.30.2 and following the standardized operating procedure5354 with minor modifications as described previously26. Sequence control and analyses were performed on 16S rRNA genes from 13C-DNA from the BP microcosms day 4 and 14, BZ day 1 and 4 according to He et al.2626. After sequencing, significant differences between labeled and unlabeled microbial communities present in heavy fractions were detected using one-way analysis of variance (ANOVA). Comparison of sample means was done using the Tukey post-hoc test at a significance level of P-value < 0.05. If the community from 13C-incubated samples were not significantly different from the Total Community or time 0 heavy fractions, these samples was excluded from further analyses. Pyrosequencing reads were deposited in the metagenomics (MG) - RAST server55 under the MG-RAST Project ID 8748.

Q-PCR of bphA

The abundance of portions of genes coding for the α subunit of biphenyl dioxygenase (BphA) were amplified from 13C-DNA after SIP with an ABI 7900 HT Fast Real-Time PCR System (Applied Biosciences, USA). Reactions were performed in a 10 μl final reaction volume containing 5 μl of Kapa SYBR FAST qPCR Master Mix (Kapa Biosystems, USA), 1 μM of each primer, approximately 5 ng of 13C-DNA. The bphA sequences were amplified using a primers 512f (tgr tbt tyg cva ayt ggg a) and 674r (tcb gcn gcr aay ttc cag tt). The reverse primer is a reverse complement sequence of a previously used primer BPHD-f356 and the forward primer was designed based on known sequences available in RDP FunGene database similarly as previously described57. The performance of the primers was verified using the genomic DNA of previously described PCB-degrading strains and the cycling conditions used were as follows: 3 min at 95 °C, 30 cycles of 5 s at 95 °C, 30 s at 56 °C, and 20 s at 72 °C, followed by a final extension at 72 °C for 5 min, and a melting curve.

Additional Information

How to cite this article: Leewis, M.-C. et al. Synergistic Processing of Biphenyl and Benzoate: Carbon Flow Through the Bacterial Community in Polychlorinated-Biphenyl-Contaminated Soil. Sci. Rep. 6, 22145; doi: 10.1038/srep22145 (2016).