INTRODUCTION
Carbofuran (2,3-dihydro-2,2-dimethyl-7-benzofuranoyl-
N-methylcarbamate), a representative broad-spectrum
N-methylcarbamate pesticide (
Fig. 1A), was first introduced in 1967 by the FMC Corporation (Princeton, NJ, USA). Carbofuran was widely used to control insect and nematode pests in crops (
1). Carbofuran is a potent inhibitor of cholinesterase in mammals (
2,
3). Moreover, the pesticide also acts as an endocrine disruptor (
4). Therefore, carbofuran was prohibited in many countries, but it is still used in some developing countries due to its insecticidal effects. Although carbofuran is not stable chemically because of hydrolysis in the environment, its relatively good water solubility (700 mg/liter at 25°C) and low adsorption (mean organic carbon-water partition coefficient [
Koc] of 30) (
https://www3.epa.gov/pesticides/endanger/litstatus/effects/carbofuran/riskanalysis.pdf) result in the contamination of surface and groundwater. Therefore, great concern and interest have been raised regarding the environmental behavior and degradation mechanisms of carbofuran.
Microbial degradation plays an important role in the elimination of carbofuran in the environment. Various carbofuran-degrading microbes, mostly bacteria and some fungi, have been isolated from diverse geographical origins. The reported carbofuran-degrading bacteria are from the genera
Pseudomonas,
Flavobacterium,
Achromobacter,
Sphingomonas,
Novosphingobium, and
Paracoccus (
5–14). Among these bacteria, only sphingomonads (bacteria of the genus
Sphingomonas and the closely related genera
Novosphingobium,
Sphingopyxis, and
Sphingobium, commonly referred to as sphingomonads [
15,
16]) can covert carbofuran to carbofuran phenol and further transform carbofuran phenol (
8,
10–13), while others just perform the hydrolysis of carbofuran to carbofuran phenol. In the
Sphingomonas sp. strain SB5, it was proposed that following the production of carbofuran phenol, a putative hydrolase has the ability to catalyze the cleavage of the furanyl ring (
10). Recently, the complete catabolic pathway of carbofuran mineralization by strain KN65.2 was proposed according to the metabolic profiling of wild-type KN65.2 and its plasposon mutants (
Fig. 1A) (
13); in strain KN65.2, carbofuran phenol is likely hydroxylated by a monooxygenase described below, leading to the cleavage of the furanyl ring. Here, the catabolic pathways before and after the cleavage of the benzene ring are called the upstream pathway and the downstream pathway, respectively.
The genetic basis of the carbofuran catabolic pathway is far from clear. Two carbofuran hydrolases, Mcd and CfdJ, which catalyze the conversion of carbofuran to carbofuran phenol, have been reported. Mcd was purified from
Achromobacter sp. strain WM111, and the corresponding gene
mcd was cloned (
17). The homologues of
mcd have been detected in phylogenetically diverse carbofuran-degrading isolates (
18), but no
mcd homologue has been found in carbofuran-degrading sphingomonads.
cfdJ was found in the genome of
Novosphingobium sp. strain KN65.2 (
19), which is a close homologue of the carbaryl hydrolase-encoding gene
cehA cloned from the
Rhizobium sp. strain AC100. CehA
AC100 catalyzes the conversion of carbaryl to 1-naphthol (
20). Although CehA
KN65.2 (CfdJ is the same as CehA
KN65.2) and CehA
AC100 differ by only four amino acids, CehA
KN65.2 but not CehA
AC100 can recognize carbofuran as a substrate (
20,
21). Hashimoto et al. found that CehA
AC100, purified from the
Rhizobium sp. strain AC100, showed no detectable activity against carbofuran (
20). Öztürk et al. revealed that none of the CehA
KN65.2 amino acid residues that differ from CehA
AC100 were silent regarding carbofuran hydrolytic activity, but the substitution of Phe152 to Leu152 proved crucial (
21). Therefore, it is likely that the close homologues of
cehA are responsible for the conversion of carbofuran to carbofuran phenol in carbofuran-degrading sphingomonads, although this must be further confirmed by genetic evidence.
To identify the genes involved in the mineralization of carbofuran in strain KN65.2, Nguyen et al. generated 27 plasposon mutants of strain KN65.2, which abolished or diminished its capability to degrade and mineralize carbofuran (
13). On the basis of the results of metabolite analysis of the mutants and sequence analysis of the genes disrupted by mini-Tn
5, it was proposed that the gene cluster
cfdABCDEFGH is involved in the mineralization of carbofuran phenol. CfdC is predicted to be a flavin-dependent monooxygenase, and CfdE is a putative dioxygenase; CfdC and CfdE may be responsible for the hydroxylation and the cleavage of the benzene ring of carbofuran phenol, respectively. However, more genetic and biochemical data are needed to validate the functional and physiological roles of CfdC and CfdE.
The
linABCDEF genes for the catabolism of hexachlorocyclohexane (HCH) are found in the HCH-mineralizing sphingomonads isolated all over the world (
22). Similarly, the genes involved in the isoproturon-catabolic pathway are highly conserved among the isoproturon-mineralizing sphingomonads (
23). Considering these facts, we hypothesized that the genes involved in carbofuran catabolism might also be conserved (likely with >95% similarity) among carbofuran-degrading sphingomonads. Both
Sphingomonas sp. strain CDS-1 and
Novosphingobium sp. strain KN65.2 can convert carbofuran and carbofuran phenol (
11,
13). However, unlike strain KN65.2, strain CDS-1 cannot utilize carbofuran or carbofuran phenol as a sole carbon source for growth. Strain CDS-1 was isolated from activated sludge in the Jiangsu province of China, and strain KN65.2 was isolated from soil sampled from a vegetable field with a long history of carbofuran treatment in the Soc Trang province of Vietnam. The draft genome sequence of strain KN65.2 was released recently (
19). In this work, the draft genome of strain CDS-1 was sequenced; the genes involved in the carbofuran catabolism were first predicted from the open reading frames (ORFs) that share ≥95% nucleic acid similarities between strains CDS-1 and KN65.2 and then were subjected to experimental validation.
DISCUSSION
Carbofuran has been widely used as an insecticide since 1967. To date, eight carbofuran-degrading sphingomonads have been described; all of these strains can transform carbofuran to carbofuran phenol and further convert carbofuran phenol (
5–14). Two of these strains (CF06 and KN65.2) can grow on carbofuran (
8,
13); only strain KN65.2 was reported to be able to mineralize the benzene ring of carbofuran (
13). In this work, both genetic and biochemical results show that in strain CDS-1, CehA
CDS-1 hydrolyzes carbofuran to carbofuran phenol and CfdC
CDS-1 hydroxylates carbofuran phenol at the benzene ring. We assume that the close homologues of
cehA and
cfdC are present in all carbofuran-degrading sphingomonads and are responsible for the initial two steps of carbofuran catabolism.
The IS element plays an important role in the evolution of the catabolic pathway of xenobiotic compounds, and the involved genes are usually associated with an IS element (
31).
cehACDS-1 and
cfdCCDS-1 are linked to IS
6100, while
cehAAC100 and
cehAKN65.2 are flanked by IS
Rsp3 (
Fig. 1). The IS
6100 element is a member of the IS
6 family and has an extremely broad host range (
32). IS
6100 has been found to be associated with the genes (
lin,
mpd, and
adoQTA1A2B) involved in the catabolism of hexachlorocyclohexane, methyl parathion, and aniline (
23,
33,
34). IS
Rsp3 containing a transposase gene pair (
istA-istB) belongs to the IS
21 family (
20). IS
Rsp3 is also adjacent to the chloroacetanilide
N-dealkylase encoding gene
cndA in
Sphingomonas sp. strain DC-6 (
35). The G+C contents of
cehA (57.9%) and
cfdC (60.5%) are obviously lower than those of the genomes of their host strains KN65.2 (63.1%) and CDS-1 (62.7%), showing the features of horizontal gene transfer. IS
6100 and IS
Rsp3 might contribute to the transposition of
cehA and
cfdC, facilitating the horizontal transfer of
cehA and
cfdC.
The unexpected finding of this study is that CehA
AC100 can hydrolyze carbofuran to carbofuran phenol, which is inconsistent with two previous reports (
20,
21). Hashimoto et al. purified CehA
AC100 from the carbaryl-degrading strain AC100 and found that the purified CehA
AC100 showed no measurable activity against carbofuran (
20). Öztürk et al. expressed CehA
AC100 and its four variants in
E. coli; it should be noted that Öztürk et al. did not purify these five proteins, and the carbofuran hydrolytic activity assay was performed using cell extracts. The results indicate that CehA
AC100 could not transform carbofuran, and only the variants containing the Phe152Leu substitution can convert carbofuran to carbofuran phenol (
21). Why the carbofuran-hydrolytic activity of CehA
AC100 was not detected in these two previous works may have two explanations. First, the HPLC retention times of carbofuran and carbofuran phenol are very close; neither the positive control nor the authentic chemical of carbofuran phenol was mentioned in the assay of CehA
AC100 purified from strain AC100 (
20). Thus, if the HPLC peaks of carbofuran and carbofuran phenol could not be separated well, the generation of carbofuran phenol may be neglected. Second, CehA
AC100 is poorly expressed in
E. coli; moreover, CehA
AC100 exhibits weak activity against carbofuran. If the soluble CehA
AC100 was present at only a very low concentration in the cell extracts, it would be difficult to detect its activity toward carbofuran. In this work, recombinant CehA
AC100 was purified to homogeneity. The result of the carbofuran-hydrolytic activity assay using the purified CehA
AC100 confirms that CehA
AC100 can convert carbofuran to carbofuran phenol. Additionally, in line with previous reports (
21), we found that the Phe152Leu substitution significantly elevated the affinity and catalytic efficiency of CehA
AC100 against carbofuran.
Among the five putative flavin reductases tested, two supported the activity of CfdC. Furthermore, CfdC activity was also supported by the flavin reductases in
S. wittichii RW1 and
P. putida KT2440. However, the
E. coli strains (DH5α or BL21) expressing CfdC failed to transform carbofuran phenol, indicating the flavin reductases of
E. coli did not support the activity of CfdC. In strain DH5α, there is no flavin reductase gene present that is related to
cfdX or CBW64_RS08340 (data not shown). These results show that CfdC has some degree of specificity in the recognition of its partner flavin reductase. This phenomenon has been described in other flavin-dependent monooxygenase systems (
36–39). For example, the
E. coli expressing only the oxygenase PheA1 of phenol hydroxylase PheA1A2 was unable to convert phenol, while coexpressing reductase PheA2 with PheA1 enabled the
E. coli to convert phenol (
36). Similarly, the
E. coli cells synthesizing 2,5-diketocamphane monooxygenase (2,5-DKCMO) or 3,6-diketocamphane monooxygenase (3,6-DKCMO) without their cognate flavin reductase showed very little activity against camphor (
37). In the absence of its partner reductase, MeaX that is involved in the catabolism of chloroacetanilide herbicides can be functionally expressed in strains
S. wittichii RW1 and
P. putida KT2440 but not in
E. coli strains (DH5α or BL21) (
38). Also, DszA that catalyzes the desulfurization of dibenzothiophene cannot be functionally expressed in
E. coli (
39). We postulate that some form of interaction between CfdC and its partner flavin reductase is required for the transfer of reduced FAD or FMN. The optimal ratio of 6His-CfdC and CfdX-6His is approximately 25:1, which is far less than that of phenol hydroxylase PheA1A2 (200:1) (
40).
pheA1 and
pheA2 are organized in one operon in the genome of
Bacillus thermoglucosidasius A7 (
36), indicating that PheA1A2 is a well-evolved, two-component flavin-dependent monooxygenase system. In contrast, there was no evidence of a gene for a flavin reductase located in the immediate vicinity of
cfdC in the genome of strain CDS-1. Although CfdX supported the activity of CfdC, CfdX is not the cognate reductase of CfdC; therefore, CfdCX is not a mature two-component monooxygenase system. Thus, PheA1 and PheA2 cooperate more efficiently than CfdC and CfdX.
The MS/MS data indicate that CfdCX hydroxylates the benzene ring of carbofuran phenol, likely at the
para position. Notably, the MS/MS analyses show that this hydroxylation event does not result in the cleavage of the furanyl ring (
Fig. 4C), which is different from previous speculations in strains KN65.2 and SB5 (
10,
13). In strain KN65.2, the authors suggested that the
ortho hydroxylation of the benzene ring of carbofuran phenol by CfdC
KN65.2 led to the cleavage of the furanyl ring (
Fig. 1A). In strain SB5, it was proposed that a putative hydrolase catalyzed the hydrolysis of the furanyl ring of carbofuran phenol. To explain these differences, genetic and biochemical evidence are required in both strains KN65.2 and SB5. In addition, since strain CDS-1 failed to grow on carbofuran or carbofuran phenol and no benzene ring cleavage product of carbofuran phenol was found during the conversion of carbofuran phenol by strain CDS-1 (data not shown), this strain lacks the genes for the downstream pathway of carbofuran catabolism, and these genes should be investigated in the sphingomonads that can mineralize carbofuran, such as strain KN65.2.
MATERIALS AND METHODS
Chemicals, bacterial strains, plasmids, and culture conditions.
Carbofuran (97% purity) was kindly supplied by Jiangsu Kuaida Agrochemical Co., Ltd. Carbofuran phenol (99% purity) was purchased from Ehrenstorfer-Schafers (Augsburg, Germany). The stock solutions of carbofuran and carbofuran phenol were prepared in methanol and sterilized by membrane filtration (pore size, 0.22 μm). The bacterial strains and plasmids used in this study are listed in
Table 2. All the bacterial strains were grown in Luria-Bertani (LB) medium. The
E. coli strains were grown at 37°C, and other bacterial strains were grown at 30°C. Mineral salt medium (MSM; pH 7.0) contained the following (g/liter): NaCl, 1.0; NH
4NO
3, 1.0; K
2HPO
4, 1.5; KH
2PO
4, 0.5; MgSO
4, 0.1; and CaCl
2·2H
2O, 0.1. The concentrations of various carbon sources supplemented in MSM are indicated in the experimental procedures section below. The following antibiotics were used for selections: streptomycin (Str), 100 μg/ml; ampicillin (Ap), 100 μg/ml; gentamicin (Gm), 30 μg/ml; chloramphenicol (Cm), 20 μg/ml; and kanamycin (Km), 50 μg/ml.
DNA manipulation, sequencing, and analysis.
The isolation and manipulation of recombinant DNA were performed using standard techniques. All the enzymes were commercial preparations and were used as specified by the supplier (TaKaRa, Dalian, China). The oligonucleotides used in this study are listed in
Table 3. The standard PCR mixture contained 10 μl PrimeSTAR buffer (5×), 4 μl deoxynucleoside triphosphate (dNTP) mix (2.5 mM each), 0.5 μM forward primer, 0.5 μM reverse primer, 10 ng template DNA, and 5 U PrimeSTAR HS DNA polymerase in a total volume of 50 μl; the cycling parameters were an initial denaturation at 98°C for 2 min and subsequent steps of 98°C for 10 s, annealing at 55 to 60°C for 10 s, and extension at 72°C for 1 min per kb, for 30 cycles total.
The total genomic DNA was extracted as described by Pitcher et al. (
41). Genome sequencing was performed by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The draft genome sequence of strain CDS-1 was produced using the Illumina MiSeq platform. The annotation was carried out using Glimmer 3.02 (
42), tRNAscan-SE version 1.3.1 (
43), and Barrnap 0.4.2 (
44). BLASTN and BLASTP were used for the nucleotide sequence and deduced amino acid identity searches, respectively. To obtain the highly conserved genes between strains CDS-1 and KN65.2, the genome sequences of strain CDS-1 were searched against the genome of strain KN65.2 using BLASTN with default parameters; the BLASTN output was further filtered by an in-house PERL script using the strict criteria of alignment length of ≥100 bp, E value of <1e−5, and similarity of ≥95%. Then, the highly conserved sequences were retrieved with a filtered BLASTN output. For the phylogenetic analysis, Clustal X 2.1 (
45) was used to align all the protein sequences. The multiple-sequence alignment was then imported into MEGA version 5.0 (
46), and the phylogenetic tree was constructed by the neighbor-joining method.
Gene inactivation and complementation.
To inactivate
cehACDS-1 (GenBank gene locus tag
CBW64_RS22415) through a single-crossover event, an 893-bp DNA fragment (in the middle of
cehACDS-1) was generated by PCR using strain CDS-1 genomic DNA as the template and the primers cehA-KO-F and cehA-KO-R. The resulting product was then cloned between the XbaI and XhoI sites of the suicide plasmid pJQ200SK (
47) to give pJQKA. pJQKA was delivered to strain CDS-1 via triparental mating (
48), and the transconjugants were selected on LB plates supplemented with Str and Gm. The mutant strain CDS-1-cehA was confirmed by PCR and DNA sequencing. pBBR-A for gene complementation was constructed by ligating the fragment containing
cehACDS-1 and its putative promoter region to the KpnI and EcoRI sites of the broad-host-range plasmid pBBR1MCS-2 (
49). The plasmid was then transferred into mutant strain CDS-1-cehA via triparental mating. Similarly,
cfdC (GenBank gene locus tag
CBW64_RS25315) was inactivated using plasmid pJQKC, and the resulting mutant strain CDS-1-cfdC was complemented by plasmid pBBR-C.
Degradation test.
The tested cells growing in LB broth overnight were washed twice with MSM and suspended in MSM containing 0.2 mM substrate to an optical density at 600 nm (OD600) of approximately 1.0. The substrate and its metabolites were extracted by dichloromethane, and anhydrous sodium sulfate was used to dehydrate the organic phase. The solvent dichloromethane was removed by a stream of nitrogen gas, and the resulting samples were redissolved in methanol. After filtering through a 0.22-μm-pore Millipore membrane, the samples were analyzed using HPLC as described below.
Expression and purification of recombinant CehAAC100 and CehACDS-1.
The expression and purification of recombinant
cehAAC100 and
cehACDS-1 were carried out according to a method described previously (
21). To express
cehA in
E. coli BL21(DE3) using the pET-28a system,
cehACDS-1 was amplified using primer pair cehA-Ex-F/cehA-Ex-R; since there is only single nucleotide transversion between
cehAAC100 and
cehACDS-1,
cehAAC100 was generated by overlap PCR using primer pairs cehA-Ex-F/TTGC-R and TTGC-F/cehA-Ex-R with
cehACDS-1 as the template. The N-terminal His-tagged CehA
AC100 and CehA
CDS-1 were purified on a Ni-nitrilotriacetic acid (NTA) agarose resin matrix (Sangon Biotech, Shanghai, China). The protein concentration was determined by the Bradford method (
50) with bovine serum albumin as the standard.
Identification of a proper reductase component for CfdCCDS-1.
Five putative FMN/FAD reductase-encoding genes were predicted from the genome of strain CDS-1. The locus tags of these genes in the GenBank database are
CBW64_RS07175,
CBW64_RS07655,
CBW64_RS08340,
CBW64_RS14135, and
CBW64_RS17905. These five genes were ligated to plasmid pET-28a(+), producing pET07175, pET07655, pET08340, pET14135, and pET17905. Then, each of these five plasmids was introduced into
E. coli BL21(DE3)/pBBR-C. The recombinant strains were grown to an OD
600 of 0.5 before 0.2 mM IPTG (isopropyl-β-
d-thiogalactopyranoside) was added. After 12 h of incubation at 16°C, the ability of the cells to convert carbofuran phenol was assayed as described above.
Synthesis and purification of recombinant CfdC and CfdX.
The sequence of cfdC was PCR amplified using primers cfdC-F and cfdC-R and cloned into the expression vector pET-28a(+) to give pET-C. The plasmids pET-C and pET17905 were transferred into E. coli BL21(DE3) for protein expression and purification. When the OD600 reached approximately 0.5, 0.2 mM IPTG was added, and further cultivation was continued for 12 h at 16°C. The cells were harvested by centrifugation, suspended in binding buffer (20 mM Tris-HCl, 300 mM NaCl, 10 mM imidazole [pH 7.5]), and lysed by sonication. After centrifugation at 12,000 × g for 30 min at 4°C, the supernatant was charged onto a 1-cm3 Ni-NTA gravity column (Sangon Biotech, Shanghai, China). After washing with washing buffer (20 mM Tris-HCl, 300 mM NaCl, 50 mM imidazole [pH 7.5]), the target proteins were eluted with 5 ml elution buffer (20 mM Tris-HCl, 300 mM NaCl, 500 mM imidazole [pH 7.5]). The purified protein was dialyzed against 20 mM Tris-HCl buffer (pH 7.5) and concentrated by ultrafiltration using the 10,000 molecular weight cutoff (MWCO) centrifugal filter (Merck Millipore, Germany). The purified proteins were detected by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE).
Enzyme assays.
All the enzyme reactions were performed in Tris-HCl buffer (20 mM, pH 7.5) with a final volume of 1 ml. In the kinetic assays of CehA
AC100 and CehA
CDS-1 against carbofuran, seven concentrations of carbofuran (0.05, 0.1, 0.25, 0.5, 0.7, 1, and 1.25 mM) were added to the reaction mixtures. Each reaction mixture contained 28 nM purified protein. The reaction was performed at 30°C for 0, 10, 20, 30, 40, and 50 min and was terminated by adding 25 μl of 2 mM HgCl
2 (
20). The consumption of carbofuran and the production of carbofuran phenol were analyzed by HPLC as described below.
NAD(P)H consumption at 340 nm was used for the determination of oxidoreductase activity of CfdX-6His. At this wavelength, the molar extinction coefficient for NAD(P)H was 6,220 M−1 · cm−1. In the kinetic assays of CfdX-6His, specific activities of CfdX for NADH were measured at nine concentrations of NADH (1, 3, 5, 10, 20, 40, 60, 80, and 100 μM) in the presence of 10 μM FMN/FAD (for the determination of the kinetic parameters for NADH) or at nine concentrations of FMN/FAD (0.1, 0.3, 0.5, 1, 3, 5, 7, 10, and 15 μM) in the presence of 200 μM NADH (for the determination of the kinetic parameters for FMN/FAD). In each of these assays, 9 nM CfdX-6His was used; the reactions were started by the addition of NAD(P)H, and the reaction mixtures were incubated at 30°C for 5, 10, 30, and 60 s.
To evaluate the relationship between the molar ratio of 6His-CfdC:CfdX-6His and carbofuran phenol hydroxylation activity, purified 6His-CfdC was kept constant at 1 μM, while the concentration of CfdX-6His increased from 0 to 0.06 μM. The concentration of carbofuran phenol was 250 μM. The reaction mixture was incubated at 30°C for 10 min and stopped by adding an equal volume of dichloromethane. To test the activity of recombinant CfdCX against carbofuran phenol, seven concentrations of substrate (25, 50, 70, 100, 150, 250, and 500 μM) were used in the presence of 200 μM NADH, 10 μM FMN/FAD, 1 μM 6His-CfdC, and 0.04 μM CfdX-6His. The reaction mixture was incubated at 30°C for 5, 10, and 15 min, and the reaction was stopped by adding an equal volume of dichloromethane. Carbofuran phenol in the organic phase was analyzed by HPLC. The products in the water phase and organic phase were detected by HPLC and MS/MS, as described below.
The kinetic data were evaluated by a nonlinear regression analysis with the Michaelis-Menten equation v = Vmax × [S]/(Km + [S]), where v is the reaction rate and [S] is the concentration of substrate; the kcat was calculated using the equation Vmax = kcat × [total enzyme concentration]. All the data were collected from three independent determinations. One unit of enzyme activity was defined as the amount of enzyme required to catalyze 1 μmol substrate per min at 30°C.
Analytical methods.
The prepared samples were filtered through a 0.22-μm-pore Millipore membrane before HPLC analysis. The analysis was carried out as follows: the separation column was Shim-pack VP-ODS C18 (250 mm by 4.6 mm), the mobile phase was methanol/water at an 80:20 (vol/vol) ratio with a flow rate of 1.0 ml/min, and an SPD-20A detector was used to measure the UV absorption at 280 nm. The mass spectrum data were collected using a TripleTOF 5600 (AB SCIEX) mass spectrometer. The metabolites were ionized by electrospray with positive polarity, and characteristic fragment ions were detected using MS/MS.
Accession number(s).
The draft genome sequence of
Sphingomonas sp. CDS-1 has been deposited in the GenBank database under accession number
NHRH00000000.