Isolation and identification of the DBP-degrading bacterial strain
Strain 0426 was isolated from a soil sample by enrichment culture using DBP as the sole source of carbon and energy. After sequential subculturing on CFMM agar plates supplemented with DBP, a pure culture was yielded. 16S rRNA gene sequence analysis revealed that this new DBP-degrading isolate belonged to the genus Glutamicibacter, a reclassified novel genus which previously belonged to Arthrobacter (Busse, 2016). Until now, although there have been no reports of a Glutamicibacter sp. capable of degrading PAEs, some Arthrobacter strains have been shown to degrade PAEs (Chatterjee and Dutta, 2008; Wu et al., 2010; Wen et al., 2014).
Optimization of culture conditions for DBP degradation by strain 0426
Because the two-factor interaction effects (X1X2) were not significant (P > 0.05) and not contained in the model (Table 1) based on the ANOVA results, the corrected quadratic polynomial equation (Eq. (4)) was adopted to describe the effects of temperature and pH on the DBP degradation rate:
Y = -712.50 + 13.08 X1 + 173.43 X2 − 0.2064 X12 − 12.560X22 (4)
Where Y is the predicted DBP degradation rate (%), X1 and X2 are the coded values for the pH, temperature, respectively.
The results of residual error diagnosis for the regression model were satisfactory (Fig. S2 A). The value of the regression coefficient (R2 = 99.24%) is large enough to prove the regression model equation that can be used to explain most of the variation in the response. The adjusted R2 (98.74%) and predicted coefficient R2 (96.99%) both indicated the accuracy of the model. As shown in Figure. 1, the effects of temperature and pH on DBP degradation rates can be observed visually in the three dimensional response surface and contour plots. The optimal values of temperature and pH predicted by the regression model Eq. (4) were 31.7 and 6.9, respectively, and the theoretical maximum DBP degradation rate was 93.27% within 10 hours (Fig. S2 B).
|
Source
|
DFa
|
SSb
|
MSc
|
F-value
|
P-value*
|
Table 1
Analysis of variance (ANOVA) for the fitted quadratic polynomial model for DBP degradation.
|
Model
|
5
|
704.720
|
140.944
|
166.87
|
0.000
|
|
Linear
|
2
|
106.444
|
53.222
|
63.01
|
0.000
|
|
X1
|
1
|
71.691
|
71.691
|
84.88
|
0.000
|
|
X2
|
1
|
34.752
|
34.752
|
41.14
|
0.001
|
|
Square
|
2
|
597.121
|
298.560
|
353.47
|
0.000
|
|
X1 X1
|
1
|
67.452
|
67.452
|
79.86
|
0.000
|
|
X2 X2
|
1
|
399.642
|
399.642
|
473.15
|
0.000
|
|
Two-factor interaction
|
1
|
1.156
|
1.156
|
1.37
|
0.295
|
|
X1 X2
|
1
|
1.156
|
1.156
|
1.37
|
0.295
|
|
Error
|
5
|
4.223
|
0.845
|
|
|
|
Lack of fit
|
3
|
3.768
|
1.256
|
5.51
|
0.157
|
|
Pure error
|
2
|
0.455
|
0.228
|
|
|
|
Total
|
10
|
708.943
|
|
|
|
a refers to degrees of freedom, b refers to sum of sequences, c refers to mean square.
R2 = 99.40% (adjusted R2 = 98.81%, predicted R2 = 95.60%)
* P<0.05 indicates that the model terms are significant.
Biodegradation of DBP and kinetic analysis
Complete degradation of 300 mg/L DBP was achieved within 12 h companied with an increase in the cell density of strain 0426 (Fig. 2A). By comparison, the non-inoculated control showed an insignificant change in the DBP amount. Similarly, DBP at a higher initial concentration ranging from 500 to 1000 mg/L was completely removed after 20–40 h (Fig. 2B), indicating the high catabolic activity and strong virulence tolerance of strain 0426. Different bacterial strains have shown divergent catabolic abilities to transform DBP. It was reported that most strains degrade 100–1000 mg/L DBP after more than 48 h (Gao and Wen, 2016). Among these strains, Gordonia sp. strain QH-11 showed strong degradation ability removing 300–750 mg/L DBP within 45 h (Jin et al., 2012). Compared with previously reported strains, strain 0426 displayed even stronger catabolic ability for DBP degradation.
During the degradation of DBP at different concentrations by strain 0426, there was no initial lag phase detected. The exponential model represents a vigorous reaction without an initial lag. In the present study, DBP degradation by strain 0426 therefore fit well with the exponential model. First-order kinetic models were constructed for the degradation of DBP at varying initial concentrations by strain 0426. The DBP degradation rate constant and the degradation half-life are shown in Table 2. The results showed that the degradation rate decreased and the degradation half-life was prolonged by an increase in the initial DBP concentration. It has been demonstrated that microbial degradation activity for DBP decreased with an increase in its concentration, and this is related to the toxicity caused by highly concentrated DBP (Ahuactzinpérez et al., 2014).
Table 2
First order kinetic analysis of DBP degradation by Glutamicibacter sp. strain 0426.
|
Initial concentration (mg/L)
|
Degradation rate constant
(h− 1)
|
Half-life
(h)
|
R2
|
|
300
|
0.2351
|
2.95
|
0.9148
|
|
500
|
0.1006
|
6.89
|
0.9769
|
|
700
|
0.0719
|
9.64
|
0.9571
|
|
1000
|
0.0552
|
12.56
|
0.9384
|
DBP biodegradation in the soil bioaugmented with strain 0426
To investigate the application potential of strain 0426, bioaugmentation of DBP-supplemented soil samples with strain 0426 was conducted. On the one hand, in the sterile soil sample bioaugmented with strain 0426, 96.6% of DBP was removed after 3 d, and 99.5% was degraded after 4 d, compared with the non-bioaugmented sterile soil sample (Fig. 3A). On the other hand, in the non-sterile soil sample, bioaugmentation with strain 0426 led to the complete removal of DBP at the 6th day (Fig. 3B). The non-sterile soil sample without bioaugmentation also showed DBP depletion; however, only 45.6% of DBP was removed after 6 d and DBP remained after this time point. Although the bioaugmented sterile-soil group removed DBP faster than the bioaugmented nonsterile-soil group in the whole experimental process, it should be noticed that the latter achieved a smaller half-time (1.10 d) than the former (1.99 d) within the first two days (Fig. 3).
Soil is thought to be a major reservoir for PAEs (Jin et al., 2015). Previous studies have shown that biodegradation of PAEs under natural environmental conditions was not as efficient as under laboratory conditions (Gao and Wen, 2016). In the present study, 99.5% of DBP (1 mg/g soil) was degraded after 4 d treatment by strain 0426, indicating the potential applications of this strain in contaminated soil. However, the degradation activity of strain 0426 was decreased in soil compared with aqueous phase tests, which might be ascribed to the adsorption of the substrate in soil organic matrices (Gao and Wen, 2016). Bioaugmentation in contaminated environments is recommended when indigenous microorganisms are not competent at cleaning up the pollutants. As shown in the present study, although indigenous bacteria can remove some of the DBP in soil, more than half of the contaminant residues remained in the soil. Bioaugmentation with strain 0426 enhanced the biodegradation process, and complete removal of DBP was achieved. During bioremediation in contaminated environments, the interaction of the inoculated microorganism with biotic and abiotic factors can influence the effect of bioaugmentation (Herrero and Stuckey, 2015), which could explain the slow-down of DBP removal after 2 days in the bioaugmented nonsterile-soil group. Therefore, the influence of biotic and abiotic factors on DBP degradation by strain 0426 in contaminated environments requires further study.
Bioinformatic analysis of the strain 0426 genome
In order to explore the genetic mechanism of PAEs catabolism by strain 0426, PAEs degrading-related genes in genome of strain 0426 were identified and compared with that of other degraders. The genome size of strain 0426 was 3.55 Mb with G + C content of 62.00%, and the annotated genome contains 3194 protein-coding sequences (CDSs), 72 RNA genes, 3 ncRNA genes, 64 tRNA genes and 47 pseudo genes (Table S4). A graphical circular genome map, which displays the general characteristic of the strain genome, was constructed (Fig. S1). The function of identified genes annotated by the NR, Swiss-Prot, GO, KEGG, and COG databases accounted for 94.2%, 59.6%, 59.4%, 46.5%, and 45.9% of the CDSs, respectively.
Based on the genome sequence, genes involved almost completed DBP degradation were identified, including many hydrolases, two gene clusters encoding for degradation of benzoic acid and protocatechuic acid. Strain 0426 contains 105 CDSs annotated as esterases or hydrolases, which may account for the strain with great potential for hydrolyzing various xenobiotics (Vamsee-Krishna and Phale, 2008). However, there were no homology of these hydrolases with dialkyl or monoalkyl PAEs hydrolases as reported, which were responsible for ester bond hydrolysis of PAEs. Protein sequences alignment among an alpha/beta fold hydrolase (WP_083586847.1) of strain 0426, two dialkyl PAEs hydrolases (AJO67804.1 and AGY55960.1), two monoalkyl PAEs hydrolases (AMJ52171.1 and BAE78500.1), and two ester hydrolases (YP_352784.1 and AAN66920.1) from different bacteria was performed. Regarding the hydrolases to be compared, the two dialkyl PAEs hydrolases can catalyze the hydrolysis of dialkyl PAEs to monoalkyl PAEs (Whangsuk et al., 2015, Jiao et al., 2013), and the two monoalkyl PAEs hydrolases catalyze monoalkyl PAEs to phthalic acid (Nahurira et al., 2018, Nishioka et al. 2006). The two ester hydrolases can concurrently catalyze the cleavage and formation of ester bond with remarkable activity and broad substrate spectrum (Liu et al. 2010). The results showed that a catalytic triad and pentapeptide motif (GX1SX2G) of WP_083586847.1 are well conserved (Fig. 4). The activities of a wide variety of ester hydrolases and typical lipases rely mainly on the catalytic triad usually formed by Ser, His and Asp residues, which is functionally (but not structurally) identical with that of trypsin and subtilisin (Arpigny et al. 1999). It was indicated that the function of WP_083586847.1 in strain 0426 was similar to that of typical lipases and other ester hydrolases, which could hydrolyze peptide linkage and ester linkage (Ahmad et al. 1995, Arpigny et al. 1999, Mahajan et al., 2019). Additionally, lipolytic enzymes are included into alpha/beta hydrolases, and the true lipases possess peculiar catalytic properties that they can markedly catalyze water-insoluble substrates in heterogeneous system and they can adsorb on the oil/water interface before hydrolysis by the changes of the enzyme's architecture during catalysis (Arpigny et al. 1999). Coincidentally, this hydrolase was encoded by a gene within pca cluster and located between pcaB and pcaL (Fig. 5). All the results may account for the marked performance of DBP hydrolysis and the corresponding intermediates by hydrolysis.
The pca cluster is involved in the degradation of protocatechuic acid, which is a central intermediate during aerobic degradation of PAEs and polycyclic aromatic hydrocarbons by bacteria (Fan et al., 2018). Similar to pca clusters of other strains, pcaGH-pcaB-pcaL-pcaIJ was present in the pca cluster of strain 0426 (Fig. 5). However, pcaF was absent but a benzoate transporter gene located on the 899 bp upstream of pcaL and the gene encoding 4-hydroxybenzoate 3-monooxygenase followed behind pcaH (Fig. 5). 4-hydroxybenzoate 3-monooxygenase is the enzyme that catalyzes 4-hydroxybenzoate to protocatechuic acid. This may indicate a pathway of DBP degradation from phthalic acid to protocatechuic acid different from the general pathway, which depends on the pht cluster as reported in Mycobacterium vanbaalenii PYR-1(Stingley et al., 2004), Terrabacter sp. DBF63 (Habe et al., 2003), Arthrobacter keyseri 12B (Eaton, 2001), and Gordonia sp. HS-NH1 (Li et al., 2016). The genome of strain 0426 contains all genes responsible for complete degradation of PAEs, which may expand its environments bioremediation potential, and are conducive to reveal the mechanism of PAEs metabolism.
Metabolites, putative enzymes, and the deduced pathway of DBP degradation
During DBP degradation by strain 0426, the intermediates generated were monitored by GC/MS at different time points. The mass spectra for the main metabolites are shown in Fig. S3, and correspond to phthalic acid (Fig. S3 A), butyl methyl phthalate (BMP, Fig. S3 B), ethyl methyl phthalate (EMP, Fig. S3 C), dimethyl phthalate (DMP, Fig. S3 D), methyl phthalate (Fig. S3 E), and butanoic acid (Fig. S3 F).
Based on the identified metabolites and the putative enzymes involved in degradation revealed by genome sequence analysis, a complete DBP degradation pathway was deduced in strain 0426 (Fig. 6). In the proposed pathway, DBP was converted to phthalic acid via hydrolysis, besides phthalic acid, BMP, EMP and DMP were formed. This process was different from the PAEs initial hydrolysis occurred in other bacteria that DBP was firstly hydrolyzed into mono-n-butyl phthalate and further hydrolyzed into phthalic acid (Jin et al., 2012; Huang et al., 2018; Mahajan et al., 2019; Wu et al., 2013; Whangsuk et al., 2015). It was indicated that there may be a special hydrolysis mechanism of DBP in strain 0426, which related to the peculiar hydrolyases harbored by strain 0426 as described above. Particularly, the hydrolase WP_083586847.1 may catalyze a reverse reaction of hydrolysis and esterification of mono-n-butyl phthalate to BMP, and phthalic acid to EMP/DMP. Phthalic acid was converted to benzoate by decarboxylation and then formed 4-hydroxybenzoate under the role of benzoate 4-monooxygenase. Protocatechuic acid was produced from 4-hydroxybenzoate by the catalysis of 4-hydroxybenzoate 3-monooxygenase (WP_073707416.1). Under the role of pca cluster, protocatechuic acid underwent successive oxygenolytic ring-cleavage, rearomatization, decarboxylation, and hydrolysis, and transformed into 3-ketoadipayl-CoA, which finally transformed into succinyl-CoA and acetyl-CoA, and entered into the TCA cycle. Although protocatechuate, 3-carboxy-cis, cis-muconate, 4-carboxy muconolactone, 3-ketoadipate enol lactone, 3-ketoadipate, and 3-ketoadipayl-CoA were not detected by GC-MS, the presence of the coding genes of the enzymes responsible for their transformation was revealed by analysis of the genome sequence of strain 0426. As to the existence of benzoic acid degradation gene cluster, benzoic acid in the proposed pathway might branch into cis-1,2-dihydroxycyclohexa-3,5-diene-1-carboxylate pathway, followed by the formation of catechol, cis, cis-muconate, and (+)-muconolactone, then joined into the 3-ketoadipate enol-lactone degradation.
In Gram-positive bacteria, phthalic acid is oxygenated to form 3,4-dihydro-3,4-dihydroxyphthalate, then dehydrogenized to 3,4-dihydroxyphthalate and finally decarboxylated to protocatechuate (Eaton, 2001; Habe et al., 2003; Stingley et al 2004; Li et al 2016). In Gram-negative bacteria, the process will undergo oxygenation and dehydrogenation at 4 and 5 position of carbons to form 4,5-dihydroxyphthalate, followed by decarboxylation to form protocatechuate (Chang and Zylstra, 1998; Nomura et al., 1992). In Rhodococcus sp. 2G, phthalic acid transformation is divided into two pathways, one is phthalic acid progresses through decarboxylation to form benzoic acid, and the other is phthalic acid conversion into protocatechuic acid like common Gram-negative bacteria. However, phthalic acid in strain 0426 was decarboxylated to form benzoic acid firstly, then twice hydroxylated to form 4-hydroxybenzoate and protocatechuic acid, which is a novel pathway for DBP degradation.