Ya Fen Wang
Hao Wen Zhu2
Ying Wang
Xiang Ling Zhang
Nora Fung Yee Tam- 1Laboratory of Basin Hydrology and Wetland Eco-restoration, School of Environmental Studies, China University of Geosciences, Wuhan, China
- 2Department of Biology and Chemistry, City University of Hong Kong, Kowloon Tong, Hong Kong
- 3School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan, China
- 4State Key Laboratory in Marine Pollution, City University of Hong Kong, Kowloon Tong, Hong Kong
Little is known about the diversity and succession of indigenous microbial community during debromination of polybrominated diphenyl ethers (PBDEs). This study examined the diversity and dynamics of microbial community structure in eight saline (mangrove and marine) and freshwater sediment microcosms exhibiting different debrominating capabilities for hexa-BDE 153, a common congener in sediments, using terminal restriction fragment length polymorphism (T-RFLP) and clone library analyses. The results showed that microbial community structure greatly differed between the saline and freshwater microcosms, likely leading to distinct variations in their debrominating capabilities and pathways. Higher relative abundances of Chloroflexi and Deltaproteobacteria succeed by Alphaproteobacteria and Betaproteobacteria were detected in the two mangrove microcosms with the fastest debrominating capabilities mainly via para pathway, respectively; the dominance of Alphaproteobacteria resulted in less accumulation of tetra-BDEs and more complete debromination of lower brominated congeners (from di- to tetra-BDEs). Meanwhile, the shifts in both microbial community structure and PBDE profiles were relatively small in the less efficient freshwater microcosms, with relatively more ortho and meta brominated products of BDE-153 resulted. Coincidently, one of the freshwater microcosms showed sudden increases of Chloroflexi and Deltaproteobacteria by the end of incubation, which synchronized with the increase in the removal rate of BDE-153. The significant relationship between microbial community structure and PBDEs was confirmed by redundancy analysis (18.7% of total variance explained, P = 0.002). However, the relative abundance of the well-known dechlorinator Dehalococcoides showed no clear correlation with the debrominating capability across different microcosms. These findings shed light in the significance of microbial community network in different saline environments on enhancement of PBDE intrinsic debromination.
Introduction
Polybrominated diphenyl ethers (PBDEs) are highly diverse, toxic and persistent organohalides, prevalent in sediments and soils around the world (Oros et al., 2005; Binelli et al., 2007; Li et al., 2012; Zhu et al., 2014a; McGrath et al., 2017). They share the structural similarity and environmental fate with other highly halogenated compounds like polychlorinated biphenyls (PCBs) and have been included in Stockholm Convention (2009) on the list of persistent organic pollutants (POPs). Concerns about their bioaccumulation and potential endocrine, reproductive and behavioral toxicity effects to human and wildlife have posed urgent needs for their remediation (Lam et al., 2010; Chakraborty and Das, 2016).
Microbial debromination of PBDEs played a major role in PBDE dissipation in soil, compared to plant uptake and other physiochemical processes (Huang et al., 2010). It has been well documented that PBDEs could be transformed to less-brominated congers under both aerobic and anaerobic conditions by bacterial isolates or consortia (Robrock et al., 2008; Lee et al., 2011; Ding et al., 2013; Zhang et al., 2013). The only key biodegradation mechanism for certain recalcitrant and highly brominated BDE congeners is reductive dehalogenation, during which the organohalides could serve as electron acceptors for bacteria to derive energy for their growth under anaerobic conditions (Mohn and Tiedje, 1992). This specific group of bacteria is known as organohalide-respiring bacteria (OHRB), which have been identified from limited bacterial phyla of Chloroflexi, Firmicutes, and Proteobacteria (Hug et al., 2013). Different OHRB exhibited highly specialized capabilities for the debromination of different BDE congeners (He et al., 2006; Robrock et al., 2009; Yang et al., 2015), suggesting that complete or effective remediation of PBDEs required the cooperation of different OHRB, and with other bacterial functional groups such as those capable of cleavage of aromatic rings. Moreover, in situ microbial remediation process would largely depend on the degradation network of the indigenous microbial community present (Lovley, 2003). However, our current knowledge of the diversity and dynamics of indigenous microbial community involved in the whole PBDE debrominating process is very limited.
Previous studies of PBDE remediation were mainly conducted in electronic wastes contaminated soils, sewage sludge and river sediments (Yen et al., 2009; Shih et al., 2012; Song et al., 2015; Stiborova et al., 2015). Microbial debromination could be enhanced by addition of co-substrates or electron donors, priming with other halogenated compounds, and stimulation with electrochemical technology (Gerecke et al., 2005; Qiu et al., 2012; Yang et al., 2013; Huang et al., 2014). More recent studies followed the application of new functional materials such as tourmaline, nanoscale zero-valent iron and carbonaceous materials during PBDE remediation (Zhang et al., 2014; Cai et al., 2015; Ma et al., 2016; Zhu et al., 2016). But exploration of PBDE debromination under natural environments were rather limited, especially in the marine/saline environments (Zanaroli et al., 2015). Nevertheless, our previous study showed that different types of sediments varying in salinity had great impacts on debromination of hexa-BDE153. More rapid and complete reductive debromination was observed in mangrove and marine sediments than in freshwater sediments (Zhu et al., 2014b). Although microbial reductive dechlorination potential and activities have been reported in marine coastal and subseafloor sediments (Futagami et al., 2013; Matturro et al., 2016a), studies focusing on the role of indigenous microbial community in reductive debromination in marine environments is scarce.
Additionally, most of the known OHRB were isolated from freshwater environments, with only seven (out of 74) from marine and estuarine sediments, and the phylogeny and functional diversity of certain genus like Dehalococcoidia were found more likely regional specific, e.g., greatly differing between marine and freshwater environments (Atashgahi et al., 2016; Zinder, 2016). Thus, the distribution of the extensively studied Dehalococcoidia in marine environments and their role in intrinsic debromination are also of particular interest to be explored.
In the present study, the structure and dynamics of indigenous microbial community were examined across eight different sediment microcosms varying in reductive debrominating capabilities for BDE-153 over a course of 90 days. We aimed to test the hypothesis if the debromination potential of different sediments was determined by their indigenous microbial community structure. Besides, the phylogeny of the microbial community involved in PBDE debromination was revealed to infer the role of the known OHRB during intrinsic debromination.
Materials and Methods
Sample Description
Sediment slurry samples were taken from the batch microcosm experiment for evaluation of the intrinsic debromination potential of different sediments collected from Hong Kong, SAR. There were eight sediment microcosms, including five mangrove sediments from Mai Po Ramsar wetland (MP), Sha Tau Kok (STK), Ting Kok (TK), Ho Chung (HC), and Tai O (TO), two freshwater pond sediments from Mai Po (MPf) and Nam Sang Wai (NSW), and one marine sediment from Sai Kung (SK) (Supplementary Figure S1). The five intertidal mangrove sediments (MP, STK, HC, TK, and TO) were dominated by Kandelia obovata, with the salinity ranged from 29 to 35‰, pH from 7.55 to 7.92 and the total organic matter (TOM) from 2.07 to 8.98%, while the two freshwater sediments (MPf and NSW) were dominated by Phragmites australis with the salinity from 6 to 8‰, pH from 7.35 to 7.48 and the TOM from 6.59 to 7.22%. The marine sediment (SK) had a salinity of 34‰, pH of 8.08 and TOM of 3.54%.
The detailed set-up of the microcosms has been described by Zhu et al. (2014b). In brief, the minimal salt medium (MSM) containing (in g L-1): K2HPO4, 0.27; KH2PO4, 0.35; NH4Cl, 2.7; MgCl26H2O, 0.1; CaCl22H2O, 0.1; and the trace elements made up of FeCl24H2O, 0.009; MnCl24H2O, 0.004; ZnCl2, 0.0014; CoCl46H2O, 0.001 and (NH4)6 Mo7O244H2O, 0.001 was prepared as described previously (Li et al., 2009). The salinity was adjusted to be the same as that in the original sediment sample by adding appropriate amounts of NaCl into the MSM. The medium was autoclaved at 121°C for 30 min, then freshly prepared L-cysteine and sodium sulfide (Na2S) (0.2 mM each) were added aseptically to eliminate residual oxygen. For each microcosm, 100 mL of the sterilized MSM were dispersed into a 250-mL Quickfit conical flask, 20 g of fresh sediment were then aseptically added in an anaerobic glove box under the flux of nitrogen gas and mixed thoroughly with the medium to form the sediment slurry. The anaerobic condition was maintained using N2 refilling method, and all the microcosms were acclimated for 2 weeks. After acclimation, 200 μL stock solution of BDE-153 (Chem Service, AccuStandard, purity >99%) at 100 mg L-1 in acetone was spiked into each microcosm to obtain a nominal concentration of BDE-153 at 1.0 mg kg-1 (freeze-dried weight). The sediment slurry was shaken after spiking to let the slurry equilibrated with BDE-153 and let the small volume of acetone (0.2%) evaporated according to the previous study of Zhong et al. (2006). The actual spiked concentration of BDE-153 at day 0 was determined to be 1.08 ± 0.08 mg kg-1 (mean and standard deviation of three replicates). The solvent control with acetone was not included, as Fagervold et al. (2005) reported that the solvent control without the target pollutants, i.e., PCBs, showed little effects on the microbial dechlorinating activity and community profiles. All the microcosms in triplicates were incubated at 28°C on a horizontal shaker at 150 rpm for 90 days, to ensure the homogeneity of the sediment slurry. At regular time intervals, a 2 mL slurry was subsampled from each microcosm using sterile disposable syringes for PBDE analysis, in parallel, another 2 mL slurry was collected at days 1, 15, 30, 60, and 90, and stored at -70°C for later microbial community analysis.
DNA Extraction and PCR Amplification
Genomic DNA was extracted from 0.25 g of sediment using the FastDNATM SPIN Kit for Soil (MP Biomedicals, Irvine, CA, United States). The bacterial 16S rDNA genes were amplified using the universal primers 8F: 5′-AGA GTT TGA TCC TGG CTC AG-3′ and 1492R: 5′-GGC TAC CTT GTT ACG ACT T-3′ (Lane, 1991). The 8F primer was labeled at the 5′ end with 6-carboxy-fluorescein phosphoramidite (FAM). The PCR cocktail composed of 10 μL 5× colorless GoTaqTM Flexi Buffer, 3 μL of MgCl2 (25 mM), 1 μL of dNTP (10 mM), 0.5 μL of each primer (10 μM), 0.25 μL of GoTaqTM Hot Start Polymerase (5 U μL-1) (Promega, Madison, WI, United States), 1 μL template DNA (3–10 ng), and the volume was made up to 50 μL using nuclease-free H2O. The PCR condition was an initial denaturation at 94°C for 3 min; 30 cycles of 94°C, 45 s denaturation; 55°C, 45 s annealing; 72°C, 60 s extension; followed by a final extension at 72°C for 7 min. The PCR products (1.5 kb) were then verified by running a 1.2% (w/v) agarose gel electrophoresis in 1× Tris-acetate-EDTA (TAE) buffer (pH 7.4).
Terminal Restriction Fragment Length Polymorphism (T-RFLP) Analysis
The PCR products were concentrated and purified using the AxyPrepTM DNA Gel Extraction Kit (Axygen, United States), which resulted in sufficient abundance of T-RFs for detection after digestion with the restriction enzyme (the next step). Aliquots of the purified PCR products (500–800 ng) were digested with the restriction enzyme HaeIII at 37°C according to the manufacture’s protocol (New England Biolabs, Inc., United Kingdom). The total volume of the digestion system was 40 μL. Digested products were precipitated by two volumes of cold absolute ethanol at -20°C overnight, centrifuged at 14,000 rpm at 4°C for 20 min and washed with 100 μL of 70% cold ethanol. After centrifugation, the DNA pellets were dissolved in 20 μL of nuclease-free water. Purified products (10 μL) were mixed with 0.5 μL of the internal size standard (ET ROX-550, Amersham, United States). This mixture was denatured for 2 min at 95°C and immediately chilled on ice before electrophoresis on a MegaBASE genetic analyzer (Amersham, United States) operated in a Genotyping mode at the Coastal Marine Laboratory, Hong Kong University of Science and Technology.
The basic T-RF calling was performed using a Genetic Profiler in the MEGABACE software package (Amersham, United States). For each sample, peaks over a threshold of 50 fluorescence units were used, and T-RFs of <60 bp and >550 bp were excluded from the analysis to avoid detection of primers and uncertainties of size determination, respectively. The raw data file exported from MEGABACE was re-formatted and then processed by T-REX software developed by Culman et al. (2009), which helped to determine a baseline threshold for the identification of true peaks over noise, align T-RFs in all samples and construct data matrices based on the area of T-RF peaks.
A phylogenetic assignment tool, TRiFLe program for in silico T-RFLP analysis with user-defined sequence sets, was adopted (Junier et al., 2008). A total of 2,533 bacterial sequences tagged with “marine sediments” were selected from the RDP database and used as the reference data set for the simulation of the T-RF size distribution and identification of the candidate species in the sediment samples.
Cloning and Phylogenetic Analysis
16S rRNA gene fragments from the genomic DNA from microcosms were amplified with the universal bacterial primer pairs of 8F and 1492R as described above. The PCR amplicons were checked by agarose gel electrophoresis, purified and cloned to pGEM-T Easy Vector System I (Promega, Madison, WI, United States) according to the manufacturer’s instructions. Positive colonies were selected randomly and sent to a company (Beijing Genomics Institute, Shenzhen, China) for sequencing. Near-complete 16S rRNA gene sequences were obtained and compared to other sequences via the Basic Local Alignment Search Tool (BLAST) of the National Centre for Biotechnology Information (NCBI1). The 16S rRNA gene sequences obtained from the clone libraries in this study have been deposited in NCBI GenBank database under the following accession numbers: MH091074- MH091326. The phylogenetic tree was constructed by the neighbor-joining method in MEGA 5.0 package as described by Tamura et al. (2011).
Quantification Analysis of Dehalococcoides sp. by qPCR
Real-time PCR was performed using an ABI 7300 real-time PCR instrument (Applied Biosystems, Carlsbad, CA, United States) with genus-specific primers for Dehalococcoides sp. (Deh467F and Deh564R) and universal primers for total bacteria (341F and 534R) (Muyzer et al., 1993; Freeborn et al., 2005). The amplification was performed in 25-μL reaction mixtures composed of 1×ABI SYBR® Green Master Mix (Applied Biosystems, Carlsbad, CA, United States), 0.2 μM of each primer and 2 μL of DNA template, containing 1–10 ng μL-1 of DNA. The amplification program and efficiency has been described in our previous study (Wang et al., 2015). The relative abundance of Dehalococcoides sp. was determined by normalizing the quantitative results of the specific genes to the total amount of bacterial 16S rRNA gene copies within the same sample (van der Zaan et al., 2010).
Statistical Analysis
The peak area data from T-RFLP analysis were analyzed by clustering and multidimensional scaling (MDS) analyses, as suggested by Culman et al. (2008), for the T-RFLP datasets with high beta diversities (≥2). Clustering and MDS analyses for the matrix of peak area were based on [Log(x+1)] transformation of percentage values and Bray–Curtis similarity. Significant differences between neighboring samples were evaluated by similarity profile (SIMPROF) tests with a significance level at 0.05 in cluster analysis. Exploratory tool similarity percentage analysis (SIMPER) was also used to study the similarity between complex T-RFLP profiles. The above analyses were performed using the PRIMER 6 software package (Clarke and Gorley, 2006). The PAST program (PAlaeontological STatistics ver.1.64) was used to calculate the Shannon and Simpson diversity indices. The relationship between the profiles of T-RFLP and different PBDE congeners were investigated using redundancy analysis (RDA), after checking the length of gradient of T-RFLP dataset (to be 3.360) by detrended correspondence analysis (DCA) in the CANOCO program for Windows v4.5 from Microcomputer Power (Ithaca, NY, United States). Monte Carlo permutation test was performed to test the significance of the first and all canonical axes in RDA.
Results
Debrominating Pattern of BDE-153 in Different Sediments
Our previous study showed that reductive debromination of BDE-153 was detected more extensive in mangrove and marine sediments (the removal rate of BDE-153 up to approximately 98% after 90 days), than in the freshwater pond sediments. The half-lives of BDE-153 in the sediment microcosms varied from 7.6 days in STK to 165.0 days in MPf (Zhu et al., 2014b). The sequential debrominating products of hexa-BDE 153 ranging from penta- to di-BDEs were found in all the sediment microcosms after 90 days, however, the dominant BDE homolog groups varied greatly among sediments with different debrominating rates (Figure 1). Lower brominated congeners (di- and tri-) were detected within 15 days in STK and TO, the two most efficient sediment microcosms, and their proportions increased after 30 days of incubation, accounted for around 20% in STK and 14–34% in TO, while BDE-153 rapidly decreased during this period. In contrast, the proportions of these lower brominated congeners in the other sediment microcosms were very low (less than 3%). Penta- and tetra-BDEs were found as the main daughter BDE congeners in STK and TO, accounting for 53–62% of the total BDE congeners. However, tetra-BDEs were more dominant than penta-BDEs in TO after 30 days of incubation but the percentages of tetra- and penta-BDEs were more or less the same in STK, suggesting TO had a more rapid transformation of BDE-153 and penta-BDEs but slower transformation of tetra- to lower brominated products, which was quite different from that in STK. In the other mangrove and marine sediment microcosms, penta- and tetra-BDEs formed in a sequential manner, and their respective proportions ranged from 7 to 41% and from less than 3 to 24%. Only 4 to 21% of penta-BDEs were found in the freshwater microcosms (NSW and MPf) after 60 days onwards, and tetra-BDEs were even lower (less than 3%), suggesting much lower transformation abilities for hexa-BDE 153 in the freshwater microcosms.
FIGURE 1. Debromination pattern of BDE-153 to less-brominated congeners over time in the eight sediment microcosms. Data plotted are the mol percentage of the peak area of a BDE congener to the total peak areas of all BDE congeners identified in GC/MS chromatograms. Mean values of three replicates are shown. STK, TO, TK, HC, and MP are intertidal mangrove sediments collected from Sha Tau Kok, Tai O, Ting Kok, Ho Chung and Mai Po, respectively; SK is marine sediment collected from Sai Kung; MPf and NSW are fresh pond sediments collected from Mai Po nature reserve and Nam Sang Wai.
T-RFLP Analysis of Total Microbial Communities
The 16S rRNA fragment genes were fingerprinted with T-RFLP to monitor the dynamics of microbial community structure across 40 different sediment samples during debromination of BDE-153. Totally, there were 270 T-RFs generated by HaeIII digestion, with an average richness (number of T-RFs in each sample) of 36.1. The sample heterogeneity (also known as beta diversity) of T-RFLP dataset was very high (to be 6.5, greater than 5) as suggested by McCune and Grace (2002). The heterogeneity of the microbial community from replicated samples of the same treatment was small according to our previous study using fatty acid methyl ester (FAME) analysis and others with T-RFLP analysis (Sipilä et al., 2008; Wang Y.F. et al., 2014). Non-linear ordination method, MDS was thus employed to resolve the similarity of our T-RFLP data, which included 40 sediment samples at five sampling times from Days 1 to 90(Figure 2).
FIGURE 2. Sample scatter plot of multidimensional scaling (MDS) analysis of T-RFLP profiles of total bacterial community from the eight sediment microcosms at Days 1, 15, 30, 60, and 90 after spiked with BDE-153. Neighboring samples clustered at a similarity level of 20 and 30% are indicated.
All the sediment samples were firstly clustered along the salinity gradient, with mangrove and marine sediments on the left side and freshwater sediments on the right side of the MDS scatter plot. The mangrove sediments could be further separated with marine sediments from SK along the second axis, except for STK sediments collected after Days 60 and 90. The average similarity for each sediment site was calculated by SIMPER analysis, showing the temporal change was greatest in STK (21.94) and minimized in NSW (39.32). From the scatter plot, all the mangrove sites showed a clear temporal shift after Days 60 and 90 during the debrominating process, as well as for the freshwater sediment microcosm MPf after 90 days. By contrast, such temporal shifts were insignificant in marine microcosm SK and another freshwater microcosm NSW at a similarity level of 30%.
Dynamics of the dominant terminal restriction fragments (T-RFs) were examined to identify the key bacterial species involved in PBDE debromination (Figure 3). For the two most active mangrove microcosms STK and TO, they shared T-RFs 206, 271 and 208 at Days 1 and 15, while the number of their shared T-RFs decreased from 23 to 9 after 90 days, with T-RFs 216, 231, and 194 dominated at STK and T-RFs 266 and 316 dominated at TO after Day 60, respectively. T-RFs 271 and 208 were also detected as dominant T-RFs in TK and HC at Day 1, but the dominant T-RFs shifted to 254, 314, and 192 bp in TK and HC after 60 days. Accordingly, the T-RFLP patterns in TK and HC were the most similar (average dissimilarity = 74.20, SIMPER) compared with other pairs of sediment microcosms. The dominant T-RFs in marine sediment SK resembled those detected in mangrove sediments, except for T-RFs 238(at D1–D15), 261 (at D30–D90) and 411. The freshwater sediment NSW also held unique dominant T-RFs such as 230 (at D1–D15) and 337 (at D60–D90). However, it was interesting to find that both mangrove and freshwater sediments collected from MP shared certain dominant T-RFs, such as 193, 215, 221–225, and 306, but these T-RFs became dominant at different timing.
FIGURE 3. Relative abundance of HaeIII digested 16S rDNA T-RFs from total bacterial communities in the eight sediment microcosms from Days 1 to 90. The symbol color and size indicate the relative abundance of each T-RF as percentage of total peak area (log transformed).
Molecular identification of the dominant T-RFs with cumulative contribution to the differences in T-RFLP profiles over time greater than 50% in each sediment microcosm were summarized (Table 1). The dominant T-RFs in mangrove and marine microcosms from Days 1 to 30 were identified to contain 16S rRNA gene sequences similar to uncultured deltaproteobacterium clones (T-RFs 206, 208, 238) and two strains from the Firmicutes family of Clostridiaceae (T-RFs 257 and 271). The two T-RFs 225 and 220 were abundant in mangrove and freshwater microcosms from Mai Po at Day 1 and both belonged to the phylum of Chloroflexi. The partial 16S rRNA sequences of dominant T-RFs from Days 60 to 90 matched with more bacterial phyla, including Actinobacteria (T-RF 231, STK), Betaproteobacteria (T-RF 198, STK), Alphaproteobacteria (T-RFs 192, 253/254, mangrove microcosms), candidate division OP8 (T-RF 261, SK), and Bacteroidetes (T-RF 411, SK). T-RFs 266 and 316 enriched in TO after Day 60 were closely related to Clostridium algidicarnis and another uncultured Firmicutes, respectively. T-RF 215 shared by the two microcosms from Mai Po after Day 30 were 98% similar in the partial 16S rRNA sequence to one uncultured deltaproteobacterium clone, which was also detected from the anaerobic microbial community in oil-polluted subtidal sediments reported by Acosta-González et al. (2013).
TABLE 1. Molecular identification of the dominant terminal restriction fragments (T-RFs) with cumulative contribution to differences of T-RFLP profiles over time greater than 50% in each sediment microcosm by SIMER analysis.
Microbial Diversity and Phylogeny Involved in Debromination
Eight bacterial 16S rRNA clone libraries were constructed from four representative anaerobic microcosms of TO, STK, SK, and MPf at Days 1 and 90, respectively. Phylogenetic analysis revealed that 261 full-length 16S rRNA sequences were related to 20 bacterial phyla, including Gamma-, Alpha-, Delta-, and Beta-proteobacteria (n = 130, 50%), Bacteroidetes (n = 28, 11%), Firmicutes (n = 27, 10%), and Chloroflexi (n = 20, 8%). Other phyla with a few clones (less than 5% of the total), such as Verrucomicrobia, Actinobacteria, Planctomycetes, Acidobacteria, Gemmatimonadetes, and Chlorobi, as well as two candidate phyla WS3 and OP8, were also present (Figure 4). The phylum-level diversity distinctly decreased in the mangrove and marine microcosms (TO, STK, and SK) from 2.3 at Days 1 to 1.8 at Day 90 in average, while that in the freshwater microcosms (MPf) increased from 1.3 to 1.9. The distribution of dominant phyla also significantly differed between the saline and freshwater microcosms over time (Figure 4). The clones within Gamma-, Delta-proteobacteria, Bacteroidetes, Chloroflexi, and Firmicutes divisions from the saline microcosm libraries were more abundant than those from the freshwater ones, with distribution percentages accounting for 90, 90, 100, 55, and 74% of the total detected in the saline libraries, respectively. Furthermore, the percentages of the clones from the above first four phyla in the saline libraries detected at Days 1 and 90 dropped from 61 to 29%, 83 to 7%, 71 to 29%, and 50 to 5%, respectively. The respective values of the Alpha-, Beta-proteobacteria, and Firmicutes clones increased from 19 to 38%, 5 to 32%, and 22 to 52% in the saline libraries, but decreased from 24 to 19%, 50 to 14%, and 19 to 7% in the freshwater ones. The bloom of Chloroflexi (from 0 to 45%) was observed in MPf at Day 90, followed by that of Deltaproteobacteria (from 0 to 10%). Cluster analysis, in good accordance with the T-RFLP profiles, revealed a close clustering of the samples from MPf at Day 90 with those from STK and SK at Day 90, although the samples from STK, TO and SK microcosms differed greatly from MPf at Day 1 (Supplementary Figure S2).
FIGURE 4. Phylogenetic composition of the eight 16S rRNA gene clone libraries constructed from STK, TO, SK, and MPf sediment microcosms at Days 1 and 90, respectively.
A phylogenetic tree of all the 16S rRNA gene clones was constructed to clearly display the phylogenetic distribution of the whole microbial community across different microcosms and their phylogenetic association with 14 known dehalogenators (Figure 5). It was further explored on the major differences in the microbial community structure between STK and TO at Day 90, as these two most efficient microcosms also showed distinct PBDE profiles. Chromatiales, Rhodocyclales, and Clostridia were ranked as the three most abundant bacterial classes in TO, while Verrucomicrobia, Rhizobiales and unclassified Alphaproteobacteria were the three most dominant classes in STK (Table 2).
FIGURE 5. Neighbor-joining phylogenetic tree of the 261 16S rRNA nucleotide sequences obtained from the STK, TO, SK, and MPf sediment microcosms at Days 1 and 90 with reference to members of 14 known dehalogenators (1) Shewanella sediminis strainHAW-EB3 (NR_074819.1); (2) Ruegeria sp. TM1040 (CP000377.1:144962-146416); (3) Anaeromyxobacter dehalogenans strain 2CP-1 (AF382396.1); (4) Desulfomonile limimaris (AF282177.1); 5. Desulfovibrio sp. TBP-1 (AF090830.1); (6) Desulfovibrio dechloracetivorans strain SF3 (NR_025078.1); (7) Acetobacterium sp. AG (JQ627627.1); (8) Dehalobacter restrictus (U84497.2); (9) Dehalobium chlorocoercia DF-1 (AF393781.1); (10) Dehalococcoides sp. BAV1 (AY165308.1); (11) Dehalococcoides mccartyi GY50 (NC_022964.1:c841590-840096); (12) Dehalococcoides sp. MB (EU073964.1); (13) Dehalococcoides mccartyi strain 195 (NR_114415.1); and (14) Dehalospirillum multivorans (X82931.1).
TABLE 2. Relative abundance (%) of classes of Proteobacteria, Chloroflexi, and Firmicutes that known to contain organohalides-respiring bacteria (OHRB).
Of the 20 OTUs of Chloroflexi, only one OTU (D90-MPf-22), identified as Dehalogenimonas lykanthroporepellens (95% similarity), fell into the class of Dehalococcoidetes and closely clustered with the most extensively studied Dehalococcoides spp. (Bootstrap value = 87). Another OTU of two Firmicutes clones (D90-TO-14 and D90-TO-22) were identified as Desulfitobacterium sp. and affiliated with Dehalobacter restrictus (U84497.2) (Bootstrap value = 100). One OTU (D90-STK-4) identified as Ruegeria sp. were also clustered with the Alphaproteobacterial OHRB Ruegeria sp. TM1040 (CP000377.1:144962-146416) (Bootstrap value = 64).
Molecular Detection and Quantification of Dehalococcoides 16S rRNA Genes
All the sediment microcosms held 16S rRNA gene copies of the total bacteria and Dehalococcoides within the range of 106–109 and 102–106, respectively, suggesting the relative abundance of Dehalococcoides sp. within the whole microbial community was less than 1%. The relative abundances of Dehalococcoides sp. were above 10-4 at Day 1 in all the sediment microcosms, followed by “down-up” fluctuations from Days 1 to 30, and the fluctuations varied at different degrees among sediment microcosms in the later stage (Supplementary Figure S3). Notably, the relative abundances of Dehalococcoides sp. were found much higher in NSW at each sampling time (around 10-4 ∼ 10-3) but became the lowest in the marine sediment SK after 60 days (as low as 10-6).
Discussion
Relationship Between Microbial Community Structure and Their Debrominating Capability
To our knowledge, this study presents the first attempt to compare the microbial community structure from different types of sediments with varying reductive debrominating capabilities. This is a crucial step to explore the key debrominating microorganisms responsible for the in situ bioremediation processes (Lovley, 2003). The distribution of the well-known OHRB, i.e., Dehalococcoides, Dehalobacter, and Desulfitobacterium, has previously been examined across debromination microcosms with 28 different soils and sediments, but their presence or absence was not necessarily related to the debrominating activity of samples (Lee and He, 2010). In the present study, the distinct different debrominating capabilities of sediment microcosms between saline and freshwater environments have been shown to rightly correspond to the differences in their microbial community structures (Figure 2). The significant relationship between microbial community structure determined by T-RFLP profiles and PBDE debrominating pattern was further confirmed using RDA, as previously demonstrated by Xu et al. (2012)). The first and all canonical axes contributed to explain 6.4 and 18.7% of the total variance in microbial species data, with F-values of 2.321 (P = 0.005) and 1. 563 (P = 0.002), respectively. The biplot with PBDEs and samples displayed clearly that the large variance of T-RFLP profiles among different samples was mainly attributed to the debromination degree in terms of total PBDE concentration and types of BDE congeners (Figure 6). This result complemented our previous study (Zhu et al., 2014b), and provided valuable evidence into the role of indigenous microbial community in the fate of PBDEs during natural attenuation. In addition, it extended our knowledge regarding PBDE debrominating capability mainly derived from studies with isolated dehalogenating bacteria or enriched cultures (He et al., 2006; Robrock et al., 2009; Zeng et al., 2010), which supported the differences in microbial community structure could lead to different debrominating capabilities, owing to the substrate specificity and specific debrominating pathways adopted by different OHRB species present.
FIGURE 6. The biplot diagram with PBDEs and samples based on RDA analysis of the T-RFLP and PBDE profiles from the eight sediment microcosms from Days 1 to 90, with the size of the sample symbols corresponding to total PBDE concentration, and the pies of the sample symbols corresponding to the distribution of different groups of T-RFs that were defined in Table 1. (red: dominant T-RFs at Days 1 ∼ 30; green: dominant T-RFs at Days 60 ∼ 90; blank: the remaining T-RFs detected).
The clone library analysis revealed that the mangrove and marine microcosms exhibiting more rapid and extensive debromination of PBDEs harbored a more diverse microbial community than that in the freshwater microcosm at Day 1. The positive correlation between bacterial diversity and their degradation efficiency has also been reported by Dell’Anno et al. (2012). More importantly, the distinct distribution patterns of the dominant bacterial phyla were observed between the saline and freshwater microcosms. Deltaproteobacteria and Chloroflexi were more abundant in the saline microcosms, while Alpha- and Beta-proteobacteria dominated the freshwater microcosm at the 1st day. According to the phylogenetic distribution of the known OHRB, Deltaproteobacteria and Chloroflexi are the two major phyla that harbor non-obligate and obligate OHRB, respectively (Atashgahi et al., 2016). The Chloroflexi strains, including Dehalococcoides and o-17/DF-1 like bacteria, have frequently been identified as the key dehalogenators in the chlorinated organics-contaminated sediments, but the complete removal of the substrates usually took a few months or years (Zanaroli et al., 2012; Wang and He, 2013). However, one study reported that the co-culture consisting of Dehalococcoides and Desulfovibrio (Deltaproteobacteria) strains exhibited a rapid and complete debromination of tetra- and penta-BDEs within 14 days (Lee et al., 2011). These findings could partially explain the rapid and extensive debrominating capabilities of microbial communities from STK and TO microcosms in the present study.
Meanwhile, Alpha- and Beta-proteobacteria that dominated in MPf at Day 1 became one of the dominant phyla in STK and TO microcosms, respectively, at Day 90, when they exhibited distinctly different PBDE profiles (Figures 1, 4). The dominance of Alphaproteobacteria has also been documented in PCB dechlorinating microcosms of marine sediments by the later stage of incubation (Matturro et al., 2016c). The cooperation between Alphaproteobacterial members with the potential PBDE debrominating strains from Deltaproteobacteria and Chloroflexi was inferred by their succeed dominance in different stages of PBDE debromination, and suggested to contribute to more rapid and complete debromination in STK microcosms compared to in TO microcosms that dominated by Beta- and Gamma-proteobacteria at the same time (Table 2). The Alpha- and Beta-proteobacterial members were reported to be capable of degrading less recalcitrant and lower brominated BDE congeners (from mono- to hexa-BDEs), such as Sphingomonas sp. PH-07, Rhodococcus jostii RHA1 and Burkholderia xenovorans LB400, owing to their aromatic ring-cleavage activity (Kim et al., 2007; Robrock et al., 2009). However, information on the reductive debromination of Alpha- and Beta-proteobacteria was rather limited. It is thus essential to gain deep insights in their metabolic interactions with the known OHRB for the complete debromination of PBDEs in contaminated sediments.
Unique distribution of Planctomycetes, Gemmatimonadetes, and Verrucomicrobia in the saline sediments were detected. These phyla have also been revealed to be actively involved in organic carbon utilization in the deep-sea environments (Li et al., 2016). Besides, most of the key T-RFs enriched were identified to belong to the known OHRB phyla, except for T-RF 261 being assigned to candidate division OP8. In our previous study on constructed mangrove microcosms for treatment of a mixture of wastewater-borne PAHs and PBDEs, a candidate division OP8 clone showing a high affinity with the Dehalococcoides mccartyi 195 was also detected (Wang et al., 2015). It was also once reported in a hydrocarbon- and chlorinated-solvent contaminated aquifer during intrinsic remediation (Dojka et al., 1998) and as representative members in mangrove sediments (Ikenaga et al., 2010). Further investigation on the role of these phyla in reductive debromination is needed and may help to discover new OHRB from marine/saline environments.
Dynamics of Microbial Community Structure and Their Debrominating Capability
It was found that the more rapid and divergent succession of microbial communities in the saline microcosms largely coincided with the changes in PBDE profiles. In contrast, the shifts in both microbial community structure and PBDE profiles were relatively small in the freshwater microcosms (Figures 1, 2). Interestingly, a sudden change in the MPf microcosm till Day 90 was detected by T-RFLP analysis, with Chloroflexi and Deltaproteobacteria greatly enriched, and the removal rate of BDE-153 also greatly increased as observed in the mangrove and marine microcosms after Day 15. These results suggested that the microbial debrominating capability was not only determined by the microbial community structure, but greatly associated with the magnitude and dynamics of microbial community shifts over time. In fact, the succession of microbial community structure has been well documented in response to oil spill and could be used for prediction of the on-going remediation progress (Kimes et al., 2014). However, much less is known about microbial community shifts in response to PBDEs, further investigation is required to provide direct evidences, such as stable-isotope probing (SIP) and functional gene data, on the contribution of the enriched functional bacterial groups to PBDE debromination.
Relic DNA (referred to large amounts of extracellular DNA or in non-intact cells) has been demonstrated to significantly influence the structure of the microbial community, and its effects depended on the pH value of soil samples, which ranged from 3.5 to 8.0 as reported by Carini et al. (2016). Although the interference from relic DNA was not determined in the present study, the interference could be reflected by the dynamic analysis of the microbial community structure in each sediment microcosm, which just meant to resolve the succession of microbial communities. The pH values of the eight sediment samples in our study were all slight alkaline, from 7.4 to 8.0, so the effect of pH on relic DNA interference was considered insignificant as the pH gradient was relatively narrow. In addition, the presence of relic DNA did not obscure the detection of distinct shifts and sediment type-specific patterns in microbial communities from our sediment microcosms varying in the debromination capability.
The Role of Known OHRB During Intrinsic Debromination
So far only one isolate from the Chloroflexi, Dehalococcoides mccartyi GY50 could couple its growth with debromination of PBDEs and lead to more rapid and complete debromination compared to the co-metabolic pattern adopted by previously reported debrominating cultures (Ding et al., 2017). Detection and quantification of the well-known OHRB in dehalogenation microcosms have been widely applied to evaluate their role in intrinsic remediation processes (van der Zaan et al., 2010). Here, we only quantified 16S rRNA genes of Dehalococcoides because of the following reasons. Most of the known OHRB were isolated from the freshwater environments, with only seven originating from marine and estuarine sediments. Among the seven marine-originated OHRB isolates, two strains, i.e., o-17/DF-1 like bacteria, belong to the genus of Dehalococcoides, while the other five belonging to Delta- and Gamma-proteobacteria were all non-obligate OHRB, according to Atashgahi et al. (2016). We focused on monitoring the abundance of obligate OHRB, as many obligate degraders were identified as the key members responsible for the rapid bioremediation in marine environments, e.g., after oil spills (Yakimov et al., 2007). However, the obligate degraders were suggested as a minor seed population that would quickly withdraw after relief of the pollution stress (Yang et al., 2016). In the present study, Dehalococcoides were not dominant in the microbial community, and their relative abundances varied during the debromination process, with no clear pattern observed in relation to their debrominating efficiency across different microcosms (Supplementary Figure S1). For the role of OHRB in intrinsic debromination, Xu et al. (2012) found that that the known dechlorinating bacteria, such as Dehalobacter, showed low abundances in the enriched PBDE-degrading communities and some unclassified members could be involved in debromination by using barcoded pyrosequencing. In another recent study on tetrachloroethene (PCE) contaminated marine sediments, it has been found that the relative abundance of non-Dehalococcoides Chloroflexi retrieved from the original sediment took the predominance over Dehalococcoides sp. that added after 150 days, suggesting the important potential role of non-Dehalococcoides Chloroflexi during natural reductive dechlorination process (Matturro et al., 2016b). These results suggested that one or more unknown highly efficient debrominators rather than Dehalococcoides-like species would be present to perform extensive debromination in our microcosms.
Furthermore, different OHRB strains from even the same genus as Dehalococcoides showed quite different dehalogenation capabilities due to the different number and types of reductive dehalogenases (RDases) genes they possess (Hug et al., 2013). Recently, novel PCB dechlorinase genes pcbA1, pcbA4, and pcbA5 have been revealed in D. mccartyi strains by Wang S. et al. (2014), together with three PBDE RDases genes pbrA1, pbrA2, and pbrA3 in D. mccartyi strain GY50 by Ding et al. (2017). These important discoveries would provide more specific biomarkers targeting reductive dehalogenase genes for evaluation the role of the known OHRB in debromination processes in more depth.
Conclusion
This is the first study comparing the microbial community structure and dynamics across different sediments during intrinsic debromination of PBDEs. Both the microbial community structure and their succession largely attributed to variations in the debrominating capabilities and pathways among sediment types. Mangrove and marine sediment microcosms with relatively higher abundance of Chloroflexi and Delta-proteobacteria showed faster debromination compared with freshwater microcosms. Moreover, nearly complete debromination of BDE-153 could be achieved with Alpha-proteobacteria being enriched in the later stage. Freshwater microcosms exhibited relatively slow debromination but the regeneration of those functional debrominating microorganisms were detected by the end of incubation. Further investigation is needed to compare the functional diversity of microbial communities involved in debromination in saline and freshwater sediments, or under different salinity gradient, combined with culture-dependent approaches and powerful deep sequencing and environmental microbiome technologies. These findings would improve our understanding on the role of indigenous microbial community in the debromination of PBDEs in natural sediments, and indicate the great potential to enrich and isolate debrominating microorganisms from mangrove sediments for further bioremediation purposes.
Author Contributions
YW and NT designed the study. YW, HZ, and YW performed the microbial and PBDE analyses. All authors wrote the manuscript and the Supplementary Materials.
Funding
This present study was supported by the research funds from National Natural Science Foundation of China (Nos. 51509225 and 41630318), Laboratory of Basin Hydrology and Wetland Eco-restoration (No. BHWER201402), Hong Kong GRF Grant 9042542 and the State Key Laboratory in Marine Pollution, City University of Hong Kong.
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 thank Prof. Pei Yuan QIAN and Ms. Yin Ki TAM from the Coastal Marine Laboratory, the Hong Kong University of Science and Technology for assistance in T-RFLP sequencing work. We also express our great thanks for the reviewers for their thorough reviews and constructive comments, which greatly improved the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00952/full#supplementary-material
References
Acosta-González, A., Rosselló-Móra, R., and Marqués, S. (2013). Characterization of the anaerobic microbial community in oil-polluted subtidal sediments: aromatic biodegradation potential after the Prestige oil spill. Environ. Microbiol. 15, 77–92. doi: 10.1111/j.1462-2920.2012.02782.x
Atashgahi, S., Lu, Y. and Smidt, H. (2016). “Overview of known organohalide-respiring bacteria—Phylogenetic diversity and environmental distribution,” in Organohalide-Respiring Bacteria, eds L. Adrian, and F. E. Löffler (Berlin: Springer), 63–105.
Binelli, A., Sarkar, S. K., Chatterjee, M., Riva, C., Parolini, M., Bhattacharya, B. D., et al. (2007). Concentration of polybrominated diphenyl ethers (PBDEs) in sediment cores of Sundarban mangrove wetland, northeastern part of Bay of Bengal (India). Mar. Pollut. Bull. 54, 1220–1229. doi: 10.1016/j.marpolbul.2007.03.021
Cai, Y., Liang, B., Fang, Z., Xie, Y., and Tsang, E. P. (2015). Effect of humic acid and metal ions on the debromination of BDE209 by nZVM prepared from steel pickling waste liquor. Front. Environ. Sci. Eng. 9, 879–887. doi: 10.1007/s11783-014-0764-8
Carini, P., Marsden, P. J., Leff, J. W., Morgan, E. E., Strickland, M. S., and Fierer, N. (2016). Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat. Microbol. 2, 1–6. doi: 10.1038/nmicrobiol.2016.242
Chakraborty, J., and Das, S. (2016). Molecular perspectives and recent advances in microbial remediation of persistent organic pollutants. Environ. Sci. Pollut. Res. Int. 23, 16883–16903. doi: 10.1007/s11356-016-6887-7
Clarke, K. R., and Gorley, R. N. (eds). (2006). PRIMER v6: User Manual/Tutorial. Plymouth: PRIMER-E.
Culman, S. W., Bukowski, R., Gauch, H. G., Cadillo-Quiroz, H., and Buckley, D. H. (2009). T-REX: software for the processing and analysis of T-RFLP data. BMC Bioinformatics 10:171. doi: 10.1186/1471-2105-10-171
Culman, S. W., Gauch, H. G., Blackwood, C. B., and Thies, J. E. (2008). Analysis of T-RFLP data using analysis of variance and ordination methods: a comparative study. J. Microbiol. Methods 75, 55–63. doi: 10.1016/j.mimet.2008.04.011
Dell’Anno, A., Beolchini, F., Rocchetti, L., Luna, G. M., and Danovaro, R. (2012). High bacterial biodiversity increases degradation performance of hydrocarbons during bioremediation of contaminated harbor marine sediments. Environ. Pollut. 167, 85–92. doi: 10.1016/j.envpol.2012.03.043
Ding, C., Chow, W. L., and He, J. (2013). Isolation of Acetobacterium sp. Strain AG, which reductively debrominates octa- and pentabrominated diphenyl ether technical mixtures. Appl. Environ. Microbiol. 79, 1110–1117. doi: 10.1128/AEM.02919-12
Ding, C., Rogers, M. J., Yang, K. L., and He, J. (2017). Loss of the ssrA genome island led to partial debromination in the PBDE respiring Dehalococcoides mccartyi strain GY50. Environ. Microbiol. 19, 2906–2915. doi: 10.1111/1462-2920.13817
Dojka, M. A., Hugenholtz, P., Haack, S. K., and Pace, N. R. (1998). Microbial diversity in a hydrocarbon-and chlorinated-solvent-contaminated aquifer undergoing intrinsic bioremediation. Appl. Environ. Microbiol. 64, 3869–3877.
Fagervold, S. K., Watts, J. E., et al. (2005). Sequential reductive dechlorination of meta-chlorinated polychlorinated biphenyl congeners in sediment microcosms by two different Chloroflexi phylotypes. Appl. Environ. Microbiol. 71, 8085–8090. doi: 10.1128/AEM.71.12.8085
Freeborn, R., West, K., Bhupathiraju, V., Chauhan, S., Rahm, B., Richardson, R., et al. (2005). Phylogenetic analysis of TCE-dechlorinating consortia enriched on a variety of electron donors. Environ. Sci. Technol. 39, 8358–8368. doi: 10.1021/es048003p
Futagami, T., Morono, Y., Terada, T., Kaksonen, A. H., and Inagaki, F. (2013). Distribution of dehalogenation activity in subseafloor sediments of the Nankai Trough subduction zone. Philos. Trans. R. Soc. B Biol. Sci. 368:20120249. doi: 10.1098/rstb.2012.0249
Gerecke, A. C., Hartmann, P. C., Heeb, N. V., Kohler, H. P., Giger, W., Schmid, P., et al. (2005). Anaerobic degradation of decabromodiphenyl ether. Environ. Sci. Technol. 39, 1078–1083. doi: 10.1021/es048634j
He, J., Robrock, K. R., and Alvarez-Cohen, L. (2006). Microbial reductive debromination of polybrominated diphenyl ethers (PBDEs). Environ. Sci. Technol. 40, 4429–4434. doi: 10.1021/es052508d
Huang, H. L., Zhang, S. Z., Christie, P., Wang, S., and Xie, M. (2010). Behavior of decabromodiphenyl ether (BDE-209) in the soil-plant system: Uptake, translocation, and metabolism in plants and dissipation in soil. Environ. Sci. Technol. 44, 663–667. doi: 10.1021/es901860r
Huang, H. W., Chang, B. V., and Lee, C. C. (2014). Reductive debromination of decabromodiphenyl ether by anaerobic microbes from river sediment. Int. Biodeter. Biodegr. 87, 60–65. doi: 10.1016/j.ibiod.2013.10.011
Hug, L. A., Maphosa, F., Leys, D., Loffler, F. E., Smidt, H., Edwards, E. A., et al. (2013). Overview of organohalide-respiring bacteria and a proposal for a classification system for reductive dehalogenases. Philos. Trans. R. Soc. B Biol. Sci. 368:20120322. doi: 10.1098/rstb.2012.0322
Ikenaga, M., Guevara, R., Dean, A. L., Pisani, C., and Boyer, J. N. (2010). Changes in community structure of sediment bacteria along the florida coastal everglades marsh-mangrove-seagrass salinity gradient. Microb. Ecol. 59, 284–295. doi: 10.1007/s00248-009-9572-2
Junier, P., Junier, T., and Witzel, K. P. (2008). TRiFLe, a program for in silico terminal restriction fragment length polymorphism analysis with user-defined sequence sets. Appl. Environ. Microbiol. 74, 6452–6456. doi: 10.1128/aem.01394-08
Kim, Y. M., Nam, I. H., Murugesan, K., Schmidt, S., Crowley, D. E., and Chang, Y. S. (2007). Biodegradation of diphenyl ether and transformation of selected brominated congeners by Sphingomonas sp. ph-07. Appl. Microbiol. Biotechnol. 77, 187–194. doi: 10.1007/s00253-007-1129-z
Kimes, N. E., Callaghan, A. V., Suflita, J. M., and Morris, P. J. (2014). Microbial transformation of the deepwater horizon oil spill-past, present, and future perspectives. Front. Microbiol. 5:603. doi: 10.3389/fmicb.2014.00603
Lam, C., Neumann, R., Shin, P. K., Au, D. W., Qian, P. Y., and Wu, R. S. (2010). Polybrominated diphenylethers (PBDEs) alter larval settlement of marine benthic polychaetes. Environ. Sci. Technol. 44, 7130–7137. doi: 10.1021/es1012615
Lane, D. J. (1991). “16S/23S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt, and M. Goodfellow (New York, NY: John Wiley & Sons), 115–175.
Lee, L. K., Ding, C., Yang, K. L., and He, J. (2011). Complete debromination of tetra- and penta-brominated diphenyl ethers by a coculture consisting of Dehalococcoides and Desulfovibrio species. Environ. Sci. Technol. 45, 8475–8482. doi: 10.1021/es201559g
Lee, L. K., and He, J. (2010). Reductive debromination of polybrominated diphenyl ethers by anaerobic bacteria from soils and sediments. Appl. Environ. Microbiol. 76, 794–802. doi: 10.1128/AEM.01872-09
Li, C. H., Zhou, H. W., Wong, Y. S., and Tam, N. F.-Y. (2009). Vertical distribution and anaerobic biodegradation of polycyclic aromatic hydrocarbons in mangrove sediments in Hong Kong, South China. Sci. Total Environ. 407, 5772–5779. doi: 10.1016/j.scitotenv.2009.07.034
Li, M., Jain, S., and Dick, G. J. (2016). Genomic and transcriptomic resolution of organic matter utilization among deep-sea bacteria in guaymas basin hydrothermal plumes. Front. Microbiol. 7:1125. doi: 10.3389/fmicb.2016.01125
Li, Y., Lin, T., Chen, Y., Hu, L., Guo, Z., and Zhang, G. (2012). Polybrominated diphenyl ethers (PBDEs) in sediments of the coastal East China Sea: occurrence, distribution and mass inventory. Environ. Pollut. 171, 155–161. doi: 10.1016/j.envpol.2012.07.039
Lovley, D. R. (2003). Cleaning up with genomics: applying molecular biology to bioremediation. Nat. Rev. Microbiol. 1, 35–44. doi: 10.1038/nrmicro731
Ma, W., Yan, Y., Ma, M., Su, F., and Zhang, Y. (2016). Effect of black carbon on the migration and biodegradation of pentabromodiphenyl ether (BDE-99) during natural groundwater recharge with treated municipal wastewater. Int. Biodeterior. Biodegrad. 113, 177–186. doi: 10.1016/j.ibiod.2016.01.016
Matturro, B., Di Lenola, M., Ubaldi, C., and Rossetti, S. (2016a). First evidence on the occurrence and dynamics of Dehalococcoides mccartyi PCB-dechlorinase genes in marine sediment during Aroclor1254 reductive dechlorination. Mar. Pollut. Bull. 112, 189–194. doi: 10.1016/j.marpolbul.2016.08.021
Matturro, B., Presta, E., and Rossetti, S. (2016b). Reductive dechlorination of tetrachloroethene in marine sediments: Biodiversity and dehalorespiring capabilities of the indigenous microbes. Sci. Total Environ. 545–546, 445–452. doi: 10.1016/j.scitotenv.2015.12.098
Matturro, B., Ubaldi, C., Grenni, P., Caracciolo, A. B., and Rossetti, S. (2016c). Polychlorinated biphenyl (PCB) anaerobic degradation in marine sediments: microcosm study and role of autochthonous microbial communities. Environ. Sci. Pollut. Res. Int. 23, 12613–12623. doi: 10.1007/s11356-015-4960-2
McCune, B. and Grace, J. B. (eds). (2002). Analysis of Ecological Communities. Gleneden Beach, OR: MJM Software Design.
McGrath, T. J., Ball, A. S., and Clarke, B. O. (2017). Critical review of soil contamination by polybrominated diphenyl ethers (PBDEs) and novel brominated flame retardants (NBFRs); concentrations, sources and congener profiles. Environ. Pollut. 230, 741–757. doi: 10.1016/j.envpol.2017.07.009
Mohn, W. W., and Tiedje, J. M. (1992). Microbial reductive dehalogenation. Microbiol. Rev. 56, 482–507. doi: 10.1111/j.1574-6941.2000.tb00693.x
Muyzer, G., de Waal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 59, 695–700.
Oros, D. R., Hoover, D., Rodigari, F., Crane, D., and Sericano, J. (2005). Levels and distribution of polybrominated diphenyl ethers in water, surface sediments, and bivalves from the San Francisco Estuary. Environ. Sci. Technol. 39, 33–41. doi: 10.1021/es048905q
Qiu, M., Chen, X., Deng, D., Guo, J., Sun, G., Mai, B., et al. (2012). Effects of electron donors on anaerobic microbial debromination of polybrominated diphenyl ethers (PBDEs). Biodegradation 23, 351–361. doi: 10.1007/s10532-011-9514-9
Robrock, K. R., Coelhan, M., Sedlak, D. L., and Alvarez-Cohen, L. (2009). Aerobic biotransformation of polybrominated diphenyl ethers (PBDEs) by bacterial isolates. Environ. Sci. Technol. 43, 5705–5711. doi: 10.1021/es900411k
Robrock, K. R., Korytár, P., and Alvarez-Cohen, L. (2008). Pathways for the anaerobic microbial debromination of polybrominated diphenyl ethers. Environ. Sci. Technol. 42, 2845–2852. doi: 10.1021/es0720917
Shih, Y. H., Chou, H. L., and Peng, Y. H. (2012). Microbial degradation of 4-monobrominated diphenyl ether with anaerobic sludge. J. Hazard. Mater. 213-214, 341–346. doi: 10.1016/j.jhazmat.2012.02.009
Sipilä, T. P., Keskinen, A. K., Akerman, M. L., Fortelius, C., Haahtela, K., and Yrjälä, K. (2008). High aromatic ring-cleavage diversity in birch rhizosphere: PAH treatment-specific changes of I.E.3 group extradiol dioxygenases and 16S rRNA bacterial communities in soil. ISME J. 2, 968–981. doi: 10.1038/ismej.2008.50
Song, M., Luo, C., Li, F., Jiang, L., Wang, Y., Zhang, D., et al. (2015). Anaerobic degradation of polychlorinated biphenyls (PCBs) and polychlorinated biphenyls ethers (PBDEs), and microbial community dynamics of electronic waste-contaminated soil. Sci. Total Environ. 502, 426–433. doi: 10.1016/j.scitotenv.2014.09.045
Stiborova, H., Vrkoslavova, J., Pulkrabova, J., Poustka, J., Hajslova, J., and Demnerova, K. (2015). Dynamics of brominated flame retardants removal in contaminated wastewater sewage sludge under anaerobic conditions. Sci. Total Environ. 533, 439–445. doi: 10.1016/j.scitotenv.2015.06.131
Stockholm Convention (2009). Stockholm convention on persistent organic pollutants (POPs). Secretariat of the Stockholm Convention on Persistent Organic Pollutants: United Nations Environment Programme (UNEP). Available at: http://www.wipo.int/edocs/lexdocs/treaties/en/unep-pop/trt_unep_pop_2.pdf
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
van der Zaan, B., Hannes, F., Hoekstra, N., Rijnaarts, H., de Vos, W. M., Smidt, H., et al. (2010). Correlation of Dehalococcoides 16S rRNA and chloroethene-reductive dehalogenase genes with geochemical conditions in chloroethene-contaminated groundwater. Appl. Environ. Microbiol. 76, 843–850. doi: 10.1128/AEM.01482-09
Wang, S., Chng, K. R., Wilm, A., Zhao, S., Yang, K. L., Nagarajan, N., et al. (2014). Genomic characterization of three unique Dehalococcoides that respire on persistent polychlorinated biphenyls. Proc. Natl. Acad. Sci. U.S.A. 111, 12103–12108. doi: 10.1073/pnas.1404845111
Wang, S., and He, J. (2013). Dechlorination of commercial PCBs and other multiple halogenated compounds by a sediment-free culture containing Dehalococcoides and Dehalobacter. Environ. Sci. Technol. 47, 10526–10534. doi: 10.1021/es4017624
Wang, Y. F., Wu, Y., Pi, N., and Tam, N. F. Y. (2014). Investigation of microbial community structure in constructed mangrove microcosms receiving wastewater-borne polycyclic aromatic hydrocarbons (PAHs) and polybrominated diphenyl ethers (PBDEs). Environ. Pollut. 187, 136–44. doi: 10.1016/j.envpol.2014.01.003
Wang, Y. F., Wu, Y., Wu, Z. B., and Tam, N. F. Y. (2015). Genotypic responses of bacterial community structure to a mixture of wastewater-borne PAHs and PBDEs in constructed mangrove microcosms. J. Hazard. Mater. 298, 91–102. doi: 10.1016/j.jhazmat.2015.05.003
Xu, M., Chen, X., Qiu, M., Zeng, X., Xu, J., Deng, D., et al. (2012). Bar-coded pyrosequencing reveals the responses of PBDE-degrading microbial communities to electron donor amendments. PLoS One 7:e30439. doi: 10.1371/journal.pone.0030439
Yakimov, M. M., Timmis, K. N., and Golyshin, P. N. (2007). Obligate oil-degrading marine bacteria. Curr. Opin. Biotechnol. 18, 257–266. doi: 10.1016/j.copbio.2007.04.006
Yang, C., Kublik, A., Weidauer, C., Seiwert, B., and Adrian, L. (2015). Reductive dehalogenation of oligocyclic phenolic bromoaromatics by Dehalococcoides mccartyi Strain CBDB1. Environ. Sci. Technol. 49, 8497–8505. doi: 10.1021/acs.est.5b01401
Yang, T., Nigro, L. M., Gutierrez, T., D’Ambrosio, L., Joye, S. B., Highsmith, R., et al. (2016). Pulsed blooms and persistent oil-degrading bacterial populations in the water column during and after the Deepwater Horizon blowout. Deep Sea Res. Part 2 Top. Stud. Oceanogr. 129, 282–291. doi: 10.1016/j.dsr2.2014.01.014
Yang, Y., Xu, M., He, Z., Guo, J., Sun, G., and Zhou, J. (2013). Microbial electricity generation enhances decabromodiphenyl ether (BDE-209) degradation. PLoS One 8:e70686. doi: 10.1371/journal.pone.0070686
Yen, J. H., Liao, W. C., Chen, W. C., and Wang, Y. S. (2009). Interaction of polybrominated diphenyl ethers (PBDEs) with anaerobic mixed bacterial cultures isolated from river sediment. J. Hazard. Mater. 165, 518–524. doi: 10.1016/j.jhazmat.2008.10.007
Zanaroli, G., Balloi, A., Negroni, A., Borruso, L., Daffonchio, D., and Fava, F. (2012). A Chloroflexi bacterium dechlorinates polychlorinated biphenyls in marine sediments under in situ-like biogeochemical conditions. J. Hazard. Mater. 209-210, 449–457. doi: 10.1016/j.jhazmat.2012.01.042
Zanaroli, G., Negroni, A., Häggblom, M. M., and Fava, F. (2015). Microbial dehalogenation of organohalides in marine and estuarine environments. Curr. Opin. Biotechnol. 33, 287–295. doi: 10.1016/j.copbio.2015.03.013
Zeng, X., Massey Simonich, S. L., Robrock, K. R., Korytár, P., Alvarez-Cohen, L., and Barofsky, D. F. (2010). Application of a congener-specific debromination model to study photodebromination, anaerobic microbial debromination, and Fe0 reduction of polybrominated diphenyl ethers. Environ. Toxicol. Chem. 29, 770–778. doi: 10.1002/etc.119
Zhang, S., Xia, X., Xia, N., Wu, S., Gao, F., and Zhou, W. (2013). Identification and biodegradation efficiency of a newly isolated 2,2′,4,4′-tetrabromodiphenyl ether (BDE-47) aerobic degrading bacterial strain. Int. Biodeterior. Biodegrad. 76, 24–31. doi: 10.1016/j.ibiod.2012.06.020
Zhang, Z., Wang, C., Li, J., Wang, B., Wu, J., Jiang, Y., et al. (2014). Enhanced bioremediation of soil from Tianjin, China, contaminated with polybrominated diethyl ethers. Environ. Sci. Pollut. Res. Int. 21, 14037–14046. doi: 10.1007/s11356-014-3313-x
Zhong, Y., Luan, T., Zhou, H., Lan, C., and Tam, N. F. Y. (2006). Metabolite production in degradation of pyrene alone or in a mixture with another polycyclic aromatic hydrocarbon by Mycobacterium sp. Environ. Toxicol. Chem. 25, 2853–2859. doi: 10.1897/06-042R.1
Zhu, B., Wu, S., Xia, X., Lu, X., Zhang, X., Xia, N., et al. (2016). Effects of carbonaceous materials on microbial bioavailability of 2,2′,4,4′-tetrabromodiphenyl ether (BDE-47) in sediments. J. Hazard. Mater. 312, 216–223. doi: 10.1016/j.jhazmat.2016.03.065
Zhu, H. W., Wang, Y., Wang, X., Luan, T. G., and Tam, N. F. Y. (2014a). Distribution and accumulation of polybrominated diphenyl ethers (PBDEs) in Hong Kong mangrove sediments. Sci. Total Environ. 468-469, 130–139. doi: 10.1016/j.scitotenv.2013.08.021
Zhu, H. W., Wang, Y., Wang, X., Luan, T. G., and Tam, N. F. Y. (2014b). Intrinsic debromination potential of polybrominated diphenyl ethers in different sediment slurries. Environ. Sci. Technol. 48, 4724–4731. doi: 10.1021/es4053818
Keywords: PBDEs, mangrove sediment, microbial community, reductive debromination, T-RFLP, salinity
Citation: Wang YF, Zhu HW, Wang Y, Zhang XL and Tam NFY (2018) Diversity and Dynamics of Microbial Community Structure in Different Mangrove, Marine and Freshwater Sediments During Anaerobic Debromination of PBDEs. Front. Microbiol. 9:952. doi: 10.3389/fmicb.2018.00952
Received: 29 December 2017; Accepted: 24 April 2018;
Published: 15 May 2018.
Edited by:
Shanquan Wang, Sun Yat-sen University, ChinaReviewed by:
Meiying Xu, Guangdong Institute of Microbiology (CAS), ChinaYue Lu, Hunan University, China
Copyright © 2018 Wang, Zhu, Wang, Zhang and Tam. 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) and the copyright owner 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: Ya Fen Wang, d2FuZ3lmQGN1Zy5lZHUuY24=