Next Article in Journal
The Influence of Diet and Obesity on Gene Expression in SLE
Next Article in Special Issue
Engineering CatM, a LysR-Type Transcriptional Regulator, to Respond Synergistically to Two Effectors
Previous Article in Journal
Inter-Individual Variability in Xenobiotic-Metabolizing Enzymes: Implications for Human Aging and Longevity
Previous Article in Special Issue
Microbial Genes for a Circular and Sustainable Bio-PET Economy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 5
10.3390/genes10050404
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biphenyl/PCB Degrading bph Genes of Ten Bacterial Strains Isolated from Biphenyl-Contaminated Soil in Kitakyushu, Japan: Comparative and Dynamic Features as Integrative Conjugative Elements (ICEs)

by
Jun Hirose
1,*,
Hidehiko Fujihara
2,
Takahito Watanabe
3,
Nobutada Kimura
4,
Hikaru Suenaga
5,
Taiki Futagami
6,
Masatoshi Goto
7,
Akiko Suyama
2 and
Kensuke Furukawa
2
1
Department of Applied Chemistry, Faculty of Engineering, University of Miyazaki, Miyazaki 889-2192, Japan
2
Department of Food and Fermentation Sciences, Faculty of Food and Nutrition Sciences, Beppu University, Beppu 874-8501, Japan
3
Research Institute for Sustainable Humanosphere, Kyoto University, Uji 611-0011, Japan
4
Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba 305-8566, Japan
5
Biotechnology Research Institute for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo 135-0064, Japan
6
Education and Research Center for Fermentation Studies, Faculty of Agriculture, Kagoshima University, Kagoshima 890-0065, Japan
7
Faculty of Agriculture, Saga University, Saga 840-8502, Japan
*
Author to whom correspondence should be addressed.
Genes 2019, 10(5), 404; https://doi.org/10.3390/genes10050404
Submission received: 30 March 2019 / Revised: 18 May 2019 / Accepted: 20 May 2019 / Published: 27 May 2019
(This article belongs to the Special Issue Genetics of Biodegradation and Bioremediation)
Browse Figures
Versions Notes

Abstract

:
We sequenced the entire genomes of ten biphenyl/PCB degrading bacterial strains (KF strains) isolated from biphenyl-contaminated soil in Kitakyushu, Japan. All the strains were Gram-negative bacteria belonging to β- and γ-proteobacteria. Out of the ten strains, nine strains carried a biphenyl catabolic bph gene cluster as integrative conjugative elements (ICEs), and they were classified into four groups based on the structural features of the bph genes. Group I (five strains) possessed bph genes that were very similar to the ones in Pseudomonas furukawaii KF707 (formerly Pseudomonas pseudoalcaligenes KF707), which is one of the best characterized biphenyl-utilizing strains. This group of strains carried salicylate catabolic sal genes that were approximately 6-kb downstream of the bph genes. Group II (two strains) possessed bph and sal genes similar to the ones in KF707, but these strains lacked the bphX region between bphC and bphD, which is involved in the downstream catabolism of biphenyl. These bph-sal clusters in groups I and II were located on an integrative conjugative element that was larger than 110 kb, and they were named ICEbph-sal. Our previous study demonstrated that the ICEbph-sal of Pseudomonas putida KF715 in group II existed both in an integrated form in the chromosome (referred to as ICEbph-salKF715 (integrated)) and in a extrachromosomal circular form (referred to as ICEbph-sal (circular)) (previously called pKF715A, 483 kb) in the stationary culture. The ICEbph-sal was transferred from KF715 into P. putida AC30 and P. putida KT2440 with high frequency, and it was maintained stably as an extrachromosomal circular form. The ICEbph-salKF715 (circular) in these transconjugants was further transferred to P. putida F39/D and then integrated into the chromosome in one or two copies. Meanwhile, group III (one strain) possessed bph genes, but not sal genes. The nucleotide sequences of the bph genes in this group were less conserved compared to the genes of the strains belonging to groups I and II. Currently, there is no evidence to indicate that the bph genes in group III are carried by a mobile element. Group IV (two strains) carried bph genes as ICEs (59–61 kb) that were similar to the genes found in Tn4371 from Cupriavidus oxalacticus A5 and ICEKKS1024677 from the Acidovorax sp. strain KKS102. Our study found that bph gene islands have integrative functions, are transferred among soil bacteria, and are diversified through modification.

1. Introduction

To date, a number of biphenyl-utilizing bacteria have been isolated and characterized in terms of the degradation of polychlorinated biphenyls (PCBs), which are serious environmental contaminants that are prevalent worldwide [1,2,3]. These strains include both Gram-negative and Gram-positive bacteria. Biphenyl catabolic enzymes co-metabolize certain PCBs into chlorobenzoic acids. It is well documented that PCB degradation is highly dependent on chlorine substitutions, such as the number and positions of the substituted chlorine [4]. Degradation capabilities are also strain dependent. For the first time, biphenyl catabolic bph genes were cloned from Pseudomonas furukawaii KF707 [5]. Since then, bph genes were cloned from various strains, including both Gram-negative and Gram-positive bacteria, and then they were analyzed in detail [3]. These studies indicated that some strains possessed bph genes that were very similar to the ones in KF707 in terms of gene organization and nucleotide sequences, although some strains possessed bph gene clusters that were different from KF707 and diversified from each other [3]. Some bph genes are located on chromosomes, whereas others are present on plasmids. The bph genes of C. oxalacticus A5 (formerly Ralstonia sp. strain A5) [6] and Acidovorax sp. KKS102 [7] are located on the ICEs. Gram-positive Rhodococcus jostii RHA1 possesses multiple bph genes on large linear plasmids [8,9].
The typical bph gene cluster shown in KF707 is composed of bphRA1A2A3A4BCX0X1X2X3D (Figure 1) [3]. Briefly, the biphenyl dioxygenase is a multi-component enzyme encoded by bphA1A2A3A4, and it catalyzes the initial oxygenation of biphenyl, converting the biphenyl into dihydrodiol, where bphA1 and bphA2 encode a large and a small subunit of the terminal dioxygenase, respectively. bphA3 encodes ferredoxin, and bphA4 encodes ferredoxin reductase. The dihydrodiol compound is then converted to a dihydroxy-compound by the dehydrogenase encoded by the bphB. The dihydroxy-compound is then degraded into 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid by the ring-cleavage dioxygenase encoded by the bphC. Then, the ring meta-cleavage compound is degraded into benzoic acid and 2-hydroxypenta-2,4-dienoic acid by the hydrolase (encoded by the bphD). BphX1X2X3 is responsible for the further degradation of 2-hydroxypenta-2,4-dienoic acid into acetyl CoA. These structural bph genes are regulated by the bphR located on the bph gene cluster [10,11] Among these bph genes, bphA1 is critically important for substrate specificity, i.e., the biodegradation capability for various aromatic compounds, including PCBs [12,13,14].
Previously, we isolated more than ten biphenyl-utilizing bacterial strains (KF strains) from biphenyl-contaminated soil in Kitakyushu, Japan [15]. Among these KF strains, we determined the complete nucleotide sequence of the P. putida KF715 genome [16], which revealed five replicons: one circular chromosome and four plasmids. Southern blot analysis indicated that the majority of the KF715 cell population carried the bph-sal cluster on its chromosome. However, a small population of cells carried the cluster on a huge extrachromosomal circular element called pKF715A (483 kb). In addition, this element carried the oriT sequence, the repA gene involved in replication, the conjugal transfer gene (tra), and the partitioning gene (par). In this study, we were interested in how the KF strains isolated from the same location carried bph gene clusters along with other catabolic genes. We performed whole genome sequencing of these strains, anticipating that their genome information would shed light on the diversity and evolution of biphenyl-utilizing bacteria. Our results indicated that specific DNA blocks, including the bph gene cluster, were integrated within glycine tRNA (tRNA-Gly) genes and that some blocks contained an integrase gene, illustrating that certain bph gene islands had integrative functions.

2. Materials and Methods

2.1. Bacterial Strains and Cultivation

The bacterial strains (KF strains) used in this study are presented in Table 1. These biphenyl/PCB degrading strains were isolated at the biphenyl-manufacturing factory in Kitakyushu, Japan [15], and were deposited to the National Biological Resource Center (NBRC). Strains KF701, KF703, KF707, KF708, and KF712 were renamed based on the 16S rRNA sequence [17,18,19,20,21]. P. putida AC30Bph+ and P. putida KT2440Bph+ were obtained through conjugation with P. putida KF715, and these two transconjugants grew on biphenyl as a sole source of carbon and energy, as described in Reference [22]. P. putida F39/D is a mutant of toluene utilizing P. putida F1 [23], in which the todD gene is defective. This strain was used as a recipient for the conjugation experiments. The growth of the KF strains and the transconjugants was examined on a basal salt agar medium with biphenyl and various aromatic compounds as described in Reference [5].

2.2. Genome Sequencing and Computational Analysis

The whole genome sequences of the KF strains were determined by the National Institute of Technology and Evaluation (NITE), using a combination of shotgun sequencing on a 454 GS FLX+ system (Roche, Basel, Switzerland) and paired-end sequencing on a HiSeq sequencing system (Illumina, San Diego, CA, USA) as previously reported in Reference [17]. The reads obtained by the two systems were assembled using the Newbler version 2.8 (Roche). The draft sequence data of the P. furukawaii KF707 and P. putida KF715 were further completed using the GenoFinisher computer program (http://www.ige.tohoku.ac.jp/joho/gf_e/). Remaining gaps between the contigs were closed using polymerase chain reaction (PCR) amplification and DNA sequencing with standard Sanger technology. The genome sequences were annotated using the RAST (Rapid Annotation using Subsystem Technology) server [30]. The identification of the coding genes was checked using a BLAST search (http://www.ncb.gov/BLAST/). Sequence comparison was performed using EasyFig Ver. 2.1 [31], and the map was generated using drawGeneArrows3 (http://www.ige.tohoku.ac.jp/joho/labhome/tool.html). The whole genome sequences of the 10 strains were deposited to the DDBJ/EMBL/GenBank under the accession numbers presented in Table 1. The nucleotide sequences of the integrative conjugative elements in KF701, KF702, KF703, KF707, KF708, KF710, KF712, and KF716 were deposited separately. Their accession numbers were LC469607, LC469608, LC469609, LC469610, LC469611, LC469612, LC469613, and LC469614, respectively.

2.3. Phylogenetic Analysis and Gene Alignment

The nucleotide sequences were aligned computationally using the ClustalW algorithm as in Reference [32]. Phylogenetic trees were generated using the neighbor-joining method with the Mega 6.0 program [33]. The trees were evaluated through bootstrap resampling (500 replicates).

2.4. Conjugation Experiments

Transfer of the Bph+ phenotype (P. putida AC30Bph+ and P. putida KT2440Bph+) by conjugation into the recipient cells (P. putida F39/D) was carried out through filter mating as described in Reference [22].

2.5. DNA Manipulation

DNA isolation, Southern blot analysis, PCR, DNA sequencing, and other DNA manipulations were performed according to standard procedures as described in Reference [34]. Pulsed-field gel electrophoresis was performed in accordance with the manufacturer’s instructions (Bio-Rad Laboratories, Hercules, CA, USA).

3. Results

3.1. KF Strains and their Genomic Features

In this study, we found that various types of bph genes were present in the ten biphenyl/PCB degrading strains isolated from the biphenyl-contaminated soil. Figure 2 shows the 16S rRNA phylogenetic tree of the ten strains. Among the ten strains, seven belonged to the Pseudomonas genus, two strains to Cupriavidus spp., and one to Comamonas spp. These were all Gram-negative bacteria belonging to β-proteobacteria (Comamonas and Cupriavidus) and γ-proteobacteria (Pseudomonas). The complete genome sequences were determined for P. furukawaii KF707 and P. putida KF715. KF707 possessed one circular chromosome of 6,242,949 bp and one plasmid (pKF707) of 59,819 bp. On the other hand, P. putida KF715 possessed one circular chromosome of 6,583,376 bp and four plasmids as previously reported in Reference [16]. The total length of the contigs (total number of contigs, >500 bp) of the remaining strains were as follows: P. abietaniphila KF701, 6,886,250 bp (140) [17]; P. aeruginosa KF702, 7,167,540 bp (91) [24]; P. putida KF703, 6,434,897 bp (135) [18]; C. basilensis KF708, 7,826,077 bp (62) [20]; C. pauculus KF709, 6,826,799 bp (227) [24]; P. toyotomiensis KF710, 5,596,721 bp (29) [27]; C. testosteroni KF712, 5,890,323 bp (97) [21]; and P. stutzeri KF716, 4,188,013 bp (30) [29].
All ten strains possessed the bph genes, while seven strains possessed the salicylate catabolic sal genes. The bph-sal cluster was localized on the chromosome in KF707. The same cluster was located on the chromosome in the majority of the KF715 cells and also existed as an extrachromosomal circular form (483,376 bp) as well in the minor part of the cells in the stationary phase culture [16]. Based on the features of the bph genes, we classified the ten strains into four groups. Group I (five strains) possessed the bph-sal cluster that was almost identical to that of KF707, where the sal genes were localized approximately 6-kb downstream of the bph genes. This group included strains KF702, KF703, KF707, KF710, and KF716. The bph gene cluster of this group was composed of bphRA1A2A3A4BCX0X1X2X3D (11.2 kb), which encodes the catabolic enzymes that degrade biphenyl into benzoate and acetyl CoA via the meta-cleavage pathway. Group II (strains KF701 and KF715) possessed the bph-sal cluster that was similar to that of KF707, although the bphX region (3.5 kb) that is involved in the metabolism of 2-hydroxypenta-2,6-dienoic acid to acetyl CoA (Figure 1, lower pathway of biphenyl metabolism) was deleted as described in Reference [22]. The structural features of the bphRABCD (7.8 kb) and the flanking regions, including the sal genes (10.2 kb), were almost identical between KF701 and KF715, indicating that the bph-sal clusters of these two strains were horizontally transferred to one another. Group III (KF709) possessed the bph genes, but not the adjacent sal genes as in the case of Burkholderia xenovorans LB400 [35], another well-characterized PCB-degrader. Overall alignment of the bph gene cluster in KF709 was similar to that of group I, but each bph component of the bph genes was relatively low. Group IV (KF708 and KF712) possessed bph genes similar to the genes found in Tn4371 from C. oxalacticus A5 [6], and the ones in ICEKKS1024677 from the Acidovorax sp. strain KKS102 [7]. Some rearrangements of the bph gene cluster were observed compared to that in groups I–III as described below. This group of strains did not carry the sal genes.
All the strains possessed the catabolic genes of benzoate, the lower pathway intermediate of biphenyl catabolism. Strains KF701, KF702, KF703, KF707, KF710, and KF715 possessed multiple benzoate catabolic genes that encoded both the extradiol cleavage meta-pathways (bza) and intradiol cleavage ortho-pathways (ben). KF708, KF709, and KF712 possessed the box genes encoding the benzoate catabolic pathway via benzoyl-CoA. Besides the bph, sal, and benzoate catabolic (bza, ben, and box) genes, these ten strains possessed catabolic genes for various aromatic compounds. Eight strains (KF701, KF703, KF707, KF708, KF710, KF712, KF715, and KF716) possessed a phenol-catabolic dmp gene cluster [36]. Five strains (KF702, KF703, KF707, KF710 and KF715) had protocatechuate catabolic pca genes. Five strains (KF703, KF707, KF708, KF709 and KF715) possessed phenylacetate catabolic paa genes. Other than the aromatic catabolic genes, many strains carried heavy metal-resistant genes. Six strains (KF701, KF702, KF708, KF710, KF712, and KF715) possessed a cluster of mer genes that were responsible for the resistance of inorganic mercury [37]. All of the ten strains possessed putative czc genes involved in the resistance of cobalt, zinc, and cadmium [38]. In KF708 and KF709, more than ten putative czc gene clusters were present. Arsenate-resistant genes were present in all the ten strains. These results, except for the arsenate-resistant genes, are summarized in Table S1.

3.2. Comparison of the bph Genes

The identities (percent) of the nucleotide sequences of the bph genes in groups I–IV strains (named as types I–IV bph genes, respectively) are shown in Tables S2–S13. The phylogenetic trees of bphA1, bphB, bphC, and bphD belonging to types I and II are shown in Figure 3. A comparative analysis of the nucleotide sequences of each gene revealed that all the bph genes were almost identical in types I and II (96.7–100%), except for that of the bphX region (3.5 kb) between bphC and bphD, which was missing in the bph gene cluster of type II (Figure 4).
The bphA1, bphA2, bphA3, bphA4, bphB, and bphC genes of KF716 were almost identical to that of the type I bph gene. However, the 3’-terminus of the bphX3 gene (754 bp) and 5’-terminus of the bphD gene (260 bp) of KF716 were different from that of type I (Figure S1a,b), and the 5’-terminus of the bphD gene (260 bp) of KF716 was almost identical to that of the type II bph gene. These genetic features are reflected in the phylogenetic tree of the bphD (Figure 3). Thus, the bph genes of KF716 showed the structural features of both groups I and II.
The KF709 (group III strain) possessed a type III bph gene similar to the one in group I, but the respective bph genes were less conserved and rearranged (Figure 4). The nucleotide sequences of the KF709 bph genes had less than 72.7% identities compared to types I and II in the bphABC region (5.9 kb). However, the bphX1, X2, X3 genes and bphD gene were more identical (88.6–90.6% in bphX1 and 74.3–84.9% in bphX2, bphX2, bphX3, or bphD) (Tables S2–S13). The bphX0 encoding putative glutathione S-transferase was not present between the bphC and bphX1, but it was located downstream of the tRNA-Gly gene (data not shown). orf3 was present between the bphA2 and bphA3 of types I–III. Despite the conserved feature of orf3, its function remains unclear [3].
The genetic features of KF708 and KF712 in the group IV strain bph genes (Type IV bph gene) were very different to types I and II, yet they were similar to those of the Acidovorax sp. strain KKS102 and C. oxalacticus A5 (Figure 5). The bph gene cluster of this type was composed of bphSX1X2X3(V)A1A2A3BCD(W)A4. Thus, the bphX1X2X3 region was located upstream of the bphA1A2A3BCD region, and bphA4 was present just downstream of the bphD (Figure 5). The insertion sequence (1190 bp) and the transposase gene (IS/tnp) were present between the bphS and bphX1 of KF708 and KKS102 [39], but these were not present in KF712 and C. oxalacticus A5. The identity of the KF708 and KF712 bph genes of type IV was compared to types I and II (Tables S2–S13). The nucleotide sequences of bphA1, bphA2, bphA3, bphA4, bphB, bphC, and bphD were conserved between types I and IV; however, the identities were less than 77%. Two unidentified gene components, bphV and bphW, were present in all the type IV bph genes, but they were not present in types I–III. The type IV bph genes possessed a transcriptional regulator, bphS [39]. It was reported that bphS acts as a repressor, whereas the bphR of Type I acts as an activator for biphenyl catabolism. The functions of bphS and bphR oppose each other, although they belong to the same GntR family [10,39].

3.3. ICEbph-sal in KF Strains Belonging to Groups I and II

The genome sequence analysis of KF701, KF702, KF703, KF707, KF715, and KF716 revealed highly conserved large “genomic islands” that included the bph-sal genes adjacent to the tRNA-Gly(CCC) gene (Figure 6). There was an 18 bp direct repeat (5’-TTCCC(T/A)(T/C)(C/T)(G/A)CCCGCTCCA-3’) on the border of the conserved region and the non-conserved region (attL and attR, Figure 7). The attL included 18 bp of the 3ʹ end of the tRNA-Gly(CCC) gene. These 18 bp direct repeat sequences could be generated by the integration of genomic islands into the chromosome. In strain KF707, attL was followed by a phage-related integrase (int) gene. The bph cluster was located just downstream of the int gene, followed by the sal gene cluster approximately 6-kb downstream, and by the bza gene approximately 49-kb downstream. The attR site was located downstream of the bza gene. Thus, the genomic island was estimated to be 122.0 kb in size (Figure 6). Genes identified as the VirB4 component (ATPase) and VirD4 component (coupling protein) of the type IV secretory pathway [40] were located 40-kb and 80-kb downstream of the tRNA-Gly(CCC) gene, respectively. The parA and parB genes encoding the replication partition proteins were present near the right end, and were proposed to act as a stabilization system for the maintenance of mobile elements in the bacterial genomes [41]. These structures corresponded to the common backbone of many integrative conjugative elements (ICEs) [42], such as ICEclc from Pseudomonas knackmussii B13 [43], which carries the catabolic genes of chlorocatechols, the tyrosine integrase gene, and type IV secretory machinery. We designated it as ICEbph-salKF707. Within this element, many putative mobile protein genes were present surrounding the sal gene cluster. To note, five IS genes at the upstream region and two IS genes at the downstream region were identified. ICEbph-sal was observed in other KF strains of group I and KF701 of group II. An ICE including the bph and sal gene clusters of P. aeruginosa KF702, named ICEbph-salKF702, was calculated to be 126.7 kb (Figure 6). Likewise, the ICEbph-salKF703 of P. putida KF703, ICEbph-salKF710 of P. toyotomiensis KF710, and ICEbph-salKF716 of P. stutzeri KF716 were calculated to be 120.8 kb, 130.3 kb, and 117.3 kb, respectively. The ICEbph-salKF701 of P. abietaniphila KF701 in the group II strain was 117.4 kb in size. As in the case of ICEbph-salKF707, ICEbph-salKF701, ICEbph-salKF702, ICEbph-salKF703, and ICEbph-salKF710 contained three gene clusters of bph, sal, and bza. ICEbph-salKF716 contained bph and sal genes, but not the bza genes. Sequence comparison revealed an inversion in ICEbph-salKF702. The sal:bza and bza:sal fusion gene clusters were found, in which half parts of the sal genes and the bza genes were replaced with each other. These fusion genes were likely generated via homologous recombination between the sal and the bza genes. ICEbph-salKF710 and ICEbph-salKF716 contained the 5-kb region encoding putative multidrug efflux pumps, but ICEbph-sals from the remaining KF strains was deficient of the corresponding region. Thus, ICEbph-sal in groups I and II contained highly conserved nucleotide sequences larger than 110 kb, which were larger than many other ICEs found in bacterial strains to date [44].
We previously reported that KF715 harbors an approximately 90-kb conjugative bph-sal gene cluster in the chromosome [22], and that the cluster could be transferred to P. putida AC30 and P. putida KT2440 with very high frequency. Whole genome sequencing of the KF715 studied here revealed that the bph-sal cluster was located on a huge 483-kb extrachromosomal element, previously designated as pKF715A [16]. An 18 bp DNA sequence identical to the 3’- end portion of a tRNA-Gly(CCC) gene as part of its attachment site (attP) was present just upstream of the int gene in this element. It was also confirmed that the bacterial integration site (attB) was present within the 3ʹ-end portion of a tRNA-Gly gene in the KF715 chromosome. Furthermore, the two SpeI digested bands, approximately 300 kb and 200 kb, hybridized with a bphA1 probe were observed as reported [16]. Since there was only one SpeI site in the element, the large, faint band was matched to 314 kb of the SpeI digested extrachromosomal element. The second major band of 220 kb was matched to that of the SpeI digested chromosome of strain KF715. This suggested that the majority of the KF715 cells carried an integrated ICEbph-salKF715, whereas fewer cells carried it as an extrachromosomal circular form. Thus, it was likely that the circular form could be obtained by recombining the attL site (present at the 3ʹ-end of the tRNA-Gly gene) and attR site (present far downstream, 483 kb from the attL site), forming the attP site.

3.4. ICEbph in the Group IV Strains

The bph genes of C. basilensis KF708 and C. testosteroni KF712 were also located on different types of ICEs. ICEbphKF708 (61.8 kb), including bph genes of KF708, was almost identical to ICEKKS1024677 of the Acidovorax sp. strain KKS102 [7] (Figure 8). ICEbphKF708 had covalently bound ends of the conserved 5ʹ-GATTTTAAG-3’ sequence (attL1 and attR1). This att sequence was identical to that of ICEKKS1024677. Nine nucleotide substitutions, three nucleotide deletions, and one nucleotide insertion were found in ICEbphKF708 when compared to ICEKKS1024677. The ICEbphKF712 (59.4 kb) carrying the bph genes of KF712 was almost identical to the Tn4371 from C. oxalaticus A5, being the first ICE found carrying the bph gene cluster [6]. ICEbphKF712 had the covalently bound ends of the conserved attL2 and attR2 sequence (5’-TTTTCAT-3’). This att sequence was identical to that of Tn4371, but it was different from the attL1 and attR1 sequence (5’-GATTTTAAG-3’) of ICEKKS1024677 and KF708. The core part, including the bph and trb gene of ICEbphKF712 in length of 33 kb, was almost identical to that of Tn4371, whereas the remaining part flanked by the attL2 site was less conserved. Major parts of these ICEs shared common structures, including the bph gene cluster; genes encoding the replication and partition proteins; parA; VirD2 component (relaxase); as well as two conjugative transfer elements, the trb gene cluster, and the tra genes (Figure 8). ICEbphKF708 and ICEKKS1024677 possessed a putative arsenate-resistant gene cluster at 15.5–18.5 kb downstream of the attL site, which encoded the transcriptional regulator, arsenate reductase, and arsenite efflux transporter. The gene cluster corresponding to them was not found in ICEbphKF712 and Tn4371.

3.5. Conjugal Transfer of Extrachromosomal ICEbph-salKF715 Into P. putida F39/D

Since P. putida AC30 Bph+ and P. putida KT2440 Bph+ harbor an extrachromosomal ICEbph-salKF715, from strain KF715, we tried to transfer this element into P. putida F39/D through filter mating. The results are shown in Figure 9. Interestingly, the Bph+ transconjugant of F39/D exhibited two bands hybridized to KF715 bphA1 DNA when KT2440Bph+ was used as a donor strain of pKF715A. On the other hand, only a single band was detected when AC30Bph+ was used as a donor strain. These hybridized bands were all different in size to that of the 310-kb SpeI fragment of pKF715A, indicating that the extrachromosomal ICEbph-salKF715 of the two donor strains were transferred into the recipient strain F39/D and then integrated into the chromosomes at different loci.

4. Discussion

Bacteria evolve through a variety of genetic events, such as mutations, intergenomic shuffling, and horizontal gene transfer. ICEs are being identified in increasing numbers via bacterial genomic analysis. Genomic islands, including ICEs, are discrete DNA segments and play an important role in bacterial evolution [42,45,46]
In this study, we revealed that bph genes differentially existed in ten biphenyl/PCB degrading strains isolated from biphenyl-contaminated soil. The types I and II bph genes of KF strains belonging to groups I and II were very similar in terms of gene organization and nucleotide sequences, except that the bphX region was missing in the type II bph gene cluster. The ICE structures carrying the bph genes were identified in several strains. Typical ICEbph-sal of approximately 120 kb was observed in the group I strains and KF701 of group II, and they were accompanied by the sal genes, approximately 6-kb downstream of the bph genes. The KF715 strain of group II carried bph and sal gene clusters in the chromosome and also on a plasmid (pKF715A, 483 kb) [16]. ICEs are thought to reintegrate into the recipient’s chromosome immediately after transfer. However, a recent study suggested that certain ICEs, such as ICEBs1 from Bacillus subtilis [47] and SXT/R391 ICEs from Vibrio cholerae [48], are capable of autonomous replication. Transconjugants of P. putida AC30Bph+ and P. putida KT2440Bph+ carry the bph-sal cluster as an extrachromosomal circular form [16]. The mobile element carrying the bph-sal cluster replicates autonomously like plasmid, maintained stably, and consists of genes sharing homologies to components of the DNA replication and stabilization machinery.
The first example of ICE-harboring genes for the degradation of xenobiotic compounds is the clc element (ICEclc, 105 kb) from P. knackmussii B13. ICEclc encodes for the catabolic pathway involved in 3- and 4-chlorocatechol degradation [43,49,50,51]. Besides this catabolic gene cluster, ICEclc has the core genes, such as type IV secretion system-encoding genes, relaxase, and integrase. ICEclc can excise through recombination between short direct repeats at either end (attL and attR). The excised ICEclc can transfer to a new recipient cell through the conjugation apparatus and it integrates into the recipient’s chromosome between the 18-bp sequence at the 3’ end of the tRNA-Gly gene on the chromosome (attB) and the identical sequence on the excised ICEclc (attP), thereby restoring the tRNA-Gly gene. Both excision and integration are mediated by the IntB13 integrase. The ICEbph-sals of group I (strains KF702, KF703, KF707, KF710, and KF716) were related to ICEclc in terms of the conjugation apparatus and the core gene set. However, the genetic organizations of ICEbph-sals are different from those of ICEclc. The integrases in the group I ICEbph-sal were almost identical (~99% in amino acid sequences), but it was as low as 59% between that of ICEbph-sal and ICEclc. Such low identity of the integrases between ICEbph-sal and ICEclc reflected the differences of the insertion sites. The ICEbph-sals of group I were inserted at the 3′- end of the tRNA-Gly gene (76 bp in length), which carries the CCC anticodon. ICEclc inserts into a number of tRNA-Gly genes, but only the genes which carry the GCC anticodon [51]. It should also be noted that the ICEbph-sal in strains belonging to groups I and II are present in Pseudomonas spp., whereas other strains than Pseudomonas possess different ICEbph as seen in group IV strains, indicating that ICEbph-sal have restricted host ranges.
In this study, we found that various types of bph genes are present in ten different strains isolated from the same soil sample. Several lines of evidence suggested that many bph genes in these strains were present on the chromosome as an integrated form. However, it is also true that certain ICEbph-sal is present stably as a plasmid. ICEbph-salKF715 was stably maintained as a circular form in the two transconjugants, P. putida AC30Bph+ and P. putida KT2440Bph+. On the other hand, ICEbph-salKF715 (circular) from P. putida AC30Bph+ and P. putida KT2440Bph+ seemed to integrate differently in another recipient P. putida F39/D. When P. putida AC30Bph+ was used as a donor strain, the largest SpeI DNA fragment of the F39/DBph+ transconjugant was hybridized with the KF715 bphA1 probe (Figure 9). When P. putida KT2440Bph+ was used as a donor strain, two copies of the bphA1 DNA were detected in the F39/DBph+ transconjugant strain at different positions (Figure 9). Investigations are currently underway to reveal how the ICEbph-salKF715 is integrated in the genome of the F39/DBph+ transconjugants. The ICEclc of the P. knackmussii strain B13 was transferred by conjugation and integrated into two nonadjacent sites on the chromosome of toluene utilizing P. putida F1 [49]. Our repeated attempts to conjugally transfer the ICEbph-sal of other strains than KF715 have not been successful. This may be due to a lack of expression of the integrase genes or mutations in certain gene(s) involved in the excision or conjugal transfer.
ICEbphKF708 and ICEbphKF712 were found in Cupriavidus and Comamonas, respectively. ICEbphKF708 was almost identical to the ICEKKS1024677 from Acidovorax sp. KKS102, and they had several nucleotide differences, indicating that this type of ICEbph could be transferred between Cupriavidus and Acidovorax. The right wing corresponding to the core legion (33 kb) of ICEbphKF712 and Tn4371 is highly conserved, but the left wing is diversified (Figure 8). This indicated that these two ICEbph were rearranged in the left wing. ICEbph-sal (from the groups I and II strain) and ICEbph (from the group IV strain) are typical ICEs, which possess type IV secretion machinery. However, there are few relationships between ICEbph-sal and ICEbph in terms of gene organization, nucleotide sequence, and size (Figure 6, Figure 8). No significant identity was detected between the integrases of ICEbph and ICEbph-sal. The region encoding two conjugative transfer components, the tra/trb genes adjacent to the bph gene cluster, was found in ICEbph. The corresponding gene clusters were not found in ICEbph-sal. Gene components downstream of the sal gene cluster are likely to be involved in conjugative transfer in ICEbph-sal; however, they have not been identified because of the lack of reliable homologous genes that are identified as conjugative transfer components in public databases. Although ICEbph-sal and ICEbph possess common genetic components involved in biphenyl catabolism, their platforms are different.
The biphenyl degrading bacteria are considered to be responsible for lignin degradation at the final stage. Lignin is a complex compound based on the phenylpropane structure and contains a variety of biphenyl related molecules. Thus, the bph genes could be very ancient and distributed across a wide range of soil bacteria. Mobilization of the bph genes in soil bacteria can be achieved through various mobile genetic elements, including ICEs, transposons, and plasmids. It is highly conceivable that these genes can be modified and rearranged in different ways in new host cells. The results in this study provide a better understanding as to how soil bacteria exchange genetic islands involved in the catabolism of aromatic compounds, as well as how such genes are rearranged and modified in the natural environment.

Supplementary Materials

The following tables and figure are available online at https://www.mdpi.com/2073-4425/10/5/404/s1. Table S1: Distribution of the catabolic genes for aromatic compounds and the heavy metal resistance genes in the biphenyl/PCB degrading KF strains. Tables S2–S13: Identity (%) of nucleotide sequence of bphR, bphA1, bphA2, bphA3, bphA4, bphB, bphC, bphX0, bphX1, bphX2, bphX3, and bphD. Figure S1: Comparison of the bphX3 (a) and bphD (b) genes belonging to types I and II.

Author Contributions

K.F. planned the genome sequencing project and performed Southern blot hybridization analysis of the transconjugants. H.F. performed the sequence analysis of KF701 and KF702. H.S. performed the sequence analysis of KF703, KF708, and KF715. N.K. performed the sequence analysis of KF707. T.W. performed the sequence analysis of KF709 and KF710. J.H. participated in the sequence analysis of KF712 and KF716 and annotation of ICEbph-sal/ICEbph. A.S., T.F., and M.G. were involved in the whole genome annotation of the KF strains and completed the manuscript. All authors read and approved the final manuscript.

Funding

This work was partly supported by a research grant from the Institute for Fermentation, Osaka (G-2018-3-020 to T.W.). The APC was funded by a research grant from Beppu University (to H.F.).

Acknowledgments

We thank Atsushi Yamazoe and Akito Nishi for their useful discussion and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Furukawa, K. Microbial degradation of polychlorinated biphenyls. In Biodegradation and Detoxification of Environmental Pollutants; Chakrabarty, A.M., Ed.; CRC Press: Boca Raton, FL, USA, 1982; pp. 33–57. [Google Scholar]
  2. Abramowicz, D.A. Aerobic and anaerobic biodegradation of PCBs: A review. Crit. Rev. Biotechnol. 1990, 10, 241–251. [Google Scholar] [CrossRef]
  3. Furukawa, K.; Fujihara, H. Microbial degradation of polychlorinated biphenyls: Biochemical and molecular features. J. Biosci. Bioeng. 2008, 105, 433–449. [Google Scholar] [CrossRef] [PubMed]
  4. Furukawa, K.; Tonomura, K.; Kamibayashi, A. Effect of chlorine substitution on the biodegradability of polychlorinated biphenyls. Appl. Env. Microbiol. 1978, 35, 223–227. [Google Scholar]
  5. Furukawa, K.; Miyazaki, T. Cloning of a gene cluster encoding biphenyl and chlorobiphenyl degradation in Pseudomonas pseudoalcaligenes. J. Bacteriol. 1986, 166, 392–398. [Google Scholar] [CrossRef]
  6. Toussaint, A.; Merlin, C.; Monchy, S.; Benotmane, M.A.; Leplae, R.; Mergeay, M.; Springael, D. The biphenyl- and 4-chlorobiphenyl-catabolic transposon Tn4371, a member of a new family of genomic islands related to IncP and Ti plasmids. Appl. Env. Microbiol. 2003, 69, 4837–4845. [Google Scholar] [CrossRef] [PubMed]
  7. Ohtsubo, Y.; Ishibashi, Y.; Naganawa, H.; Hirokawa, S.; Atobe, S.; Nagata, Y.; Tsuda, M. Conjugal transfer of polychlorinated biphenyl/biphenyl degradation genes in Acidovorax sp. strain KKS102, which are located on an integrative and conjugative element. J. Bacteriol. 2012, 194, 4237–4248. [Google Scholar] [CrossRef]
  8. Masai, E.; Sugiyama, K.; Iwashita, N.; Shimizu, S.; Hauschild, J.E.; Hatta, T.; Kimbara, K.; Yano, K.; Fukuda, M. The bphDEF meta-cleavage pathway genes involved in biphenyl/polychlorinated biphenyl degradation are located on a linear plasmid and separated from the initial bphACB genes in Rhodococcus sp. strain RHA1. Gene 1997, 187, 141–149. [Google Scholar] [CrossRef]
  9. Shimizu, S.; Kobayashi, H.; Masai, E.; Fukuda, M. Characterization of the 450-kb linear plasmid in a polychlorinated biphenyl degrader, Rhodococcus sp. strain RHA1. Appl. Env. Microbiol. 2001, 67, 2021–2028. [Google Scholar] [CrossRef] [PubMed]
  10. Watanabe, T.; Inoue, R.; Kimura, N.; Furukawa, K. Versatile transcription of biphenyl catabolic bph operon in Pseudomonas pseudoalcaligenes KF707. J. Biol. Chem. 2000, 275, 31016–31023. [Google Scholar] [CrossRef]
  11. Fujihara, H.; Yoshida, H.; Matsunaga, T.; Goto, M.; Furukawa, K. Cross-regulation of biphenyl- and salicylate-catabolic genes by two regulatory systems in Pseudomonas pseudoalcaligenes KF707. J. Bacteriol. 2006, 188, 4690–4697. [Google Scholar] [CrossRef]
  12. Kumamaru, T.; Suenaga, H.; Mitsuoka, M.; Watanabe, T.; Furukawa, K. Enhanced degradation of polychlorinated biphenyls by directed evolution of biphenyl dioxygenase. Nat. Biotechnol. 1998, 16, 663. [Google Scholar] [CrossRef] [PubMed]
  13. Suenaga, H.; Goto, M.; Furukawa, K. Emergence of multifunctional oxygenase activities by random priming recombination. J. Biol. Chem. 2001, 276, 22500–22506. [Google Scholar] [CrossRef] [PubMed]
  14. Suenaga, H.; Mitsuoka, M.; Ura, Y.; Watanabe, T.; Furukawa, K. Directed evolution of biphenyl dioxygenase: Emergence of enhanced degradation capacity for benzene, toluene, and alkylbenzenes. J. Bacteriol. 2001, 183, 5441–5444. [Google Scholar] [CrossRef]
  15. Furukawa, K.; Hayase, N.; Taira, K.; Tomizuka, N. Molecular relationship of chromosomal genes encoding biphenyl/polychlorinated biphenyl catabolism: Some soil bacteria possess a highly conserved bph operon. J. Bacteriol. 1989, 171, 5467–5472. [Google Scholar] [CrossRef]
  16. Suenaga, H.; Fujihara, H.; Kimura, N.; Hirose, J.; Watanabe, T.; Futagami, T.; Goto, M.; Shimodaira, J.; Furukawa, K. Insights into the genomic plasticity of Pseudomonas putida KF715, a strain with unique biphenyl-utilizing activity and genome instability properties. Environ. Microbiol. Rep. 2017, 9, 589–598. [Google Scholar] [CrossRef]
  17. Fujihara, H.; Yamazoe, A.; Hosoyama, A.; Suenaga, H.; Kimura, N.; Hirose, J.; Watanabe, T.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of Pseudomonas abietaniphila KF701 (NBRC110664), a polychlorinated biphenyl-degrading bacterium isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00473-15. [Google Scholar] [CrossRef] [PubMed]
  18. Suenaga, H.; Yamazoe, A.; Hosoyama, A.; Kimura, N.; Hirose, J.; Watanabe, T.; Fujihara, H.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of the polychlorinated biphenyl-degrading bacterium Pseudomonas putida KF703 (NBRC 110666) isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00142-15. [Google Scholar] [CrossRef]
  19. Kimura, N.; Watanabe, T.; Suenaga, H.; Fujihara, H.; Futagami, T.; Goto, M.; Hanada, S.; Hirose, J. Pseudomonas furukawaii sp. nov., a polychlorinated biphenyl-degrading bacterium isolated from biphenyl-contaminated soil in Japan. Int. J. Syst. Evol. Microbiol. 2018, 68, 1429–1435. [Google Scholar] [CrossRef]
  20. Suenaga, H.; Yamazoe, A.; Hosoyama, A.; Kimura, N.; Hirose, J.; Watanabe, T.; Fujihara, H.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of the polychlorinated biphenyl-degrading bacterium Cupriavidus basilensis KF708 (NBRC 110671) isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00143-15. [Google Scholar] [CrossRef]
  21. Hirose, J.; Yamazoe, A.; Hosoyama, A.; Kimura, N.; Suenaga, H.; Watanabe, T.; Fujihara, H.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of the polychlorinated biphenyl-degrading bacterium Comamonas testosteroni KF712 (NBRC 110673). Genome Announc. 2015, 3, e01214-15. [Google Scholar] [CrossRef]
  22. Nishi, A.; Tominaga, K.; Furukawa, K. A 90-kilobase conjugative chromosomal element coding for biphenyl and salicylate catabolism in Pseudomonas putida KF715. J. Bacteriol. 2000, 182, 1949–1955. [Google Scholar] [CrossRef]
  23. Zylstra, G.J.; Gibson, D.T. Toluene degradation by Pseudomonas putida F1. Nucleotide sequence of the todC1C2BADE genes and their expression in Escherichia coli. J. Biol. Chem. 1989, 264, 14940–14946. [Google Scholar]
  24. Fujihara, H.; Yamazoe, A.; Hosoyama, A.; Suenaga, H.; Kimura, N.; Hirose, J.; Watanabe, T.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of Pseudomonas aeruginosa KF702 (NBRC 110665), a polychlorinated biphenyl-degrading bacterium isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00517-15. [Google Scholar] [CrossRef]
  25. Triscari-Barberi, T.; Simone, D.; Calabrese, F.M.; Attimonelli, M.; Hahn, K.R.; Amoako, K.K.; Turner, R.J.; Fedi, S.; Zannoni, D. Genome sequence of the polychlorinated-biphenyl degrader Pseudomonas pseudoalcaligenes KF707. J. Bacteriol. 2012, 194, 4426–4427. [Google Scholar] [CrossRef]
  26. Watanabe, T.; Yamazoe, A.; Hosoyama, A.; Fujihara, H.; Suenaga, H.; Hirose, J.; Futagami, T.; Goto, M.; Kimura, N.; Furukawa, K. Draft genome sequence of Cupriavidus pauculus strain KF709, a biphenyl-utilizing bacterium isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00222-15. [Google Scholar] [CrossRef]
  27. Watanabe, T.; Yamazoe, A.; Hosoyama, A.; Fujihara, H.; Suenaga, H.; Hirose, J.; Futagami, T.; Goto, M.; Kimura, N.; Furukawa, K. Draft genome sequence of Pseudomonas toyotomiensis KF710, a polychlorinated biphenyl-degrading bacterium isolated from biphenyl-contaminated soil. Genome Announc. 2015, 3, e00223-15. [Google Scholar] [CrossRef]
  28. Suenaga, H.; Yamazoe, A.; Hosoyama, A.; Kimura, N.; Hirose, J.; Watanabe, T.; Fujihara, H.; Futagami, T.; Goto, M.; Furukawa, K. Complete genome sequence of the polychlorinated biphenyl-degrading bacterium Pseudomonas putida KF715 (NBRC 110667) isolated from biphenyl-contaminated soil. Genome Announc. 2017, 5, e01624-16. [Google Scholar] [CrossRef]
  29. Hirose, J.; Yamazoe, A.; Hosoyama, A.; Kimura, N.; Suenaga, H.; Watanabe, T.; Fujihara, H.; Futagami, T.; Goto, M.; Furukawa, K. Draft genome sequence of the polychlorinated biphenyl-degrading bacterium Pseudomonas stutzeri KF716 (NBRC 110668). Genome Announc. 2015, 3, e01215-15. [Google Scholar] [CrossRef]
  30. Aziz, R.K.; Bartels, D.; Best, A.A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M.; et al. The RAST server: Rapid annotations using subsystems technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef]
  31. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef]
  32. Li, K.-B. ClustalW-MPI: ClustalW analysis using distributed and parallel computing. Bioinformatics 2003, 19, 1585–1586. [Google Scholar] [CrossRef] [Green Version]
  33. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  34. Sambrook, J.; Russell, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; CSHL Press: Cold Spring Harbor, NY, USA, 2001. [Google Scholar]
  35. Erickson, B.D.; Mondello, F.J. Nucleotide sequencing and transcriptional mapping of the genes encoding biphenyl dioxygenase, a multicomponent polychlorinated-biphenyl-degrading enzyme in Pseudomonas strain LB400. J. Bacteriol. 1992, 174, 2903–2912. [Google Scholar] [CrossRef] [PubMed]
  36. Shingler, V.; Powlowski, J.; Marklund, U. Nucleotide sequence and functional analysis of the complete phenol/3,4-dimethylphenol catabolic pathway of Pseudomonas sp. strain CF600. J. Bacteriol. 1992, 174, 711–724. [Google Scholar] [CrossRef]
  37. Brown, N.L.; Misra, T.K.; Winnie, J.N.; Schmidt, A.; Seiff, M.; Silver, S. The nucleotide sequence of the mercuric resistance operons of plasmid R100 and transposon Tn501: further evidence for mer genes which enhance the activity of the mercuric ion detoxification system. Mol. Gen. Genet. 1986, 202, 143–151. [Google Scholar] [CrossRef]
  38. Intorne, A.C.; de Oliveira, M.V.V.; de M Pereira, L.; de Souza Filho, G.A. Essential role of the czc determinant for cadmium, cobalt and zinc resistance in Gluconacetobacter diazotrophicus PAl 5. Int. Microbiol. Off. J. Span. Soc. Microbiol. 2012, 15, 69–78. [Google Scholar]
  39. Ohtsubo, Y.; Delawary, M.; Kimbara, K.; Takagi, M.; Ohta, A.; Nagata, Y. BphS, a key transcriptional regulator of bph genes involved in polychlorinated biphenyl/biphenyl degradation in Pseudomonas sp. KKS102. J. Biol. Chem. 2001, 276, 36146–36154. [Google Scholar] [CrossRef]
  40. Christie, P.J.; Whitaker, N.; González-Rivera, C. Mechanism and structure of the bacterial type IV secretion systems. Biochim. Biophys. Acta 2014, 1843, 1578–1591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Mavrodi, D.V.; Loper, J.E.; Paulsen, I.T.; Thomashow, L.S. Mobile genetic elements in the genome of the beneficial rhizobacterium Pseudomonas fluorescens Pf-5. BMC Microbiol. 2009, 9, 8. [Google Scholar] [CrossRef]
  42. Johnson, C.M.; Grossman, A.D. Integrative and conjugative elements (ICEs): What they do and how they work. Annu. Rev. Genet. 2015, 49, 577–601. [Google Scholar] [CrossRef]
  43. Gaillard, M.; Vallaeys, T.; Vorhölter, F.J.; Minoia, M.; Werlen, C.; Sentchilo, V.; Pühler, A.; van der Meer, J.R. The clc element of Pseudomonas sp. strain B13, a genomic island with various catabolic properties. J. Bacteriol. 2006, 188, 1999–2013. [Google Scholar] [CrossRef]
  44. Bellanger, X.; Payot, S.; Leblond-Bourget, N.; Guédon, G. Conjugative and mobilizable genomic islands in bacteria: Evolution and diversity. FEMS Microbiol. Rev. 2014, 38, 720–760. [Google Scholar] [CrossRef]
  45. Juhas, M.; van der Meer, J.R.; Gaillard, M.; Harding, R.M.; Hood, D.W.; Crook, D.W. Genomic islands: Tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol. Rev. 2009, 33, 376–393. [Google Scholar] [CrossRef]
  46. Wozniak, R.A.F.; Waldor, M.K. Integrative and conjugative elements: Mosaic mobile genetic elements enabling dynamic lateral gene flow. Nat. Rev. Microbiol. 2010, 8, 552–563. [Google Scholar] [CrossRef]
  47. Lee, C.A.; Babic, A.; Grossman, A.D. Autonomous plasmid-like replication of a conjugative transposon. Mol. Microbiol. 2010, 75, 268–279. [Google Scholar] [CrossRef] [Green Version]
  48. Carraro, N.; Poulin, D.; Burrus, V. Replication and active partition of integrative and conjugative elements (ICEs) of the SXT/R391 family: The line between ICEs and conjugative plasmids is getting thinner. PLoS Genet. 2015, 11, e1005298. [Google Scholar] [CrossRef]
  49. Ravatn, R.; Studer, S.; Springael, D.; Zehnder, A.J.B.; van der Meer, J.R. Chromosomal integration, tandem amplification, and deamplification in Pseudomonas putida F1 of a 105-kilobase genetic element containing the chlorocatechol degradative genes from Pseudomonas sp. strain B13. J. Bacteriol. 1998, 180, 4360–4369. [Google Scholar]
  50. Sentchilo, V.; Ravatn, R.; Werlen, C.; Zehnder, A.J.B.; van der Meer, J.R. Unusual integrase gene expression on the clc genomic island in Pseudomonas sp. strain B13. J. Bacteriol. 2003, 185, 4530–4538. [Google Scholar] [CrossRef]
  51. Sentchilo, V.; Czechowska, K.; Pradervand, N.; Minoia, M.; Miyazaki, R.; van der Meer, J.R. Intracellular excision and reintegration dynamics of the ICEclc genomic island of Pseudomonas knackmussii sp. strain B13. Mol. Microbiol. 2009, 72, 1293–1306. [Google Scholar] [CrossRef]
Genes 10 00404 g001 550
Figure 1. Catabolic pathway of the biphenyl degradation and organization of the bph gene cluster in P. furukawaii KF707. Compounds: I, biphenyl; II, 2,3-dihydroxy-4-phenylhexa-4,6-diene (dihydrodiol compound); III, 2,3-dihydroxybiphenyl; IV, 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (biphenyl meta-cleavage compound: HOPD); V, benzoic acid; VI, 2-hydroxypenta-2,4-dienoic acid; VII, 4-hydroxy-2-oxovalerate. Enzymes: BphA1-BphA2-BphA3-BphA4, biphenyl dioxygenase; BphB, dihydrodiol dehydrogenase; BphC, 2,3-dihydroxybiphenyl dioxygenase; BphX0, glutathione S-transferase; BphX1, 2-hydroxypenta-2,4-dienoate hydratase; BphX2, acetaldehyde dehydrogenase (acylating); BphX3, 4-hydoxy-2-oxovalerate aldolase; BphD, 2-hydroxy-6-oxo-6-phenylhexa-2,4-dieonic acid hydrolase. The BphR protein, belonging to the GntR family, is a transcriptional regulator involved in the expression of bphR and bphX0X1X2X3D. The function of orf3 remains unclear.
Figure 1. Catabolic pathway of the biphenyl degradation and organization of the bph gene cluster in P. furukawaii KF707. Compounds: I, biphenyl; II, 2,3-dihydroxy-4-phenylhexa-4,6-diene (dihydrodiol compound); III, 2,3-dihydroxybiphenyl; IV, 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (biphenyl meta-cleavage compound: HOPD); V, benzoic acid; VI, 2-hydroxypenta-2,4-dienoic acid; VII, 4-hydroxy-2-oxovalerate. Enzymes: BphA1-BphA2-BphA3-BphA4, biphenyl dioxygenase; BphB, dihydrodiol dehydrogenase; BphC, 2,3-dihydroxybiphenyl dioxygenase; BphX0, glutathione S-transferase; BphX1, 2-hydroxypenta-2,4-dienoate hydratase; BphX2, acetaldehyde dehydrogenase (acylating); BphX3, 4-hydoxy-2-oxovalerate aldolase; BphD, 2-hydroxy-6-oxo-6-phenylhexa-2,4-dieonic acid hydrolase. The BphR protein, belonging to the GntR family, is a transcriptional regulator involved in the expression of bphR and bphX0X1X2X3D. The function of orf3 remains unclear.
Genes 10 00404 g001
Genes 10 00404 g002 550
Figure 2. 16S rRNA phylogenetic trees of the ten KF strains. The multiple alignment outputs were used to generate neighbor-joining phylogenetic trees using MEGA 6.0 [33]. The bar indicates expected nucleotide substitutions per site. Numbers indicate the percentage occurrence of the branch in the bootstrapped trees on 500 replicates.
Figure 2. 16S rRNA phylogenetic trees of the ten KF strains. The multiple alignment outputs were used to generate neighbor-joining phylogenetic trees using MEGA 6.0 [33]. The bar indicates expected nucleotide substitutions per site. Numbers indicate the percentage occurrence of the branch in the bootstrapped trees on 500 replicates.
Genes 10 00404 g002
Genes 10 00404 g003 550
Figure 3. Phylogenetic tree of bphA1, bphB, bphC, and bphD of KF strains in groups I and II (types I and II bph genes). The multiple alignment outputs were used to generate neighbor-joining phylogenetic trees using MEGA 6.0 [33]. The bar indicates expected nucleotide substitutions per site.
Figure 3. Phylogenetic tree of bphA1, bphB, bphC, and bphD of KF strains in groups I and II (types I and II bph genes). The multiple alignment outputs were used to generate neighbor-joining phylogenetic trees using MEGA 6.0 [33]. The bar indicates expected nucleotide substitutions per site.
Genes 10 00404 g003
Genes 10 00404 g004 550
Figure 4. Comparison of the bph gene clusters of KF strains belonging to groups I–III (types I–III). The shading in pink to red shows the identity (65–100%) of the gene clusters as indicated at the bottom of the figure. R, bphR; A1, bphA1; A2, bphA2; A3, bphA3; A4, bphA4; B, bphB; C, bphC; D, bphD; X0, bphX0, X1, bphX1; X2, bphX2; and X3, bphX3.
Figure 4. Comparison of the bph gene clusters of KF strains belonging to groups I–III (types I–III). The shading in pink to red shows the identity (65–100%) of the gene clusters as indicated at the bottom of the figure. R, bphR; A1, bphA1; A2, bphA2; A3, bphA3; A4, bphA4; B, bphB; C, bphC; D, bphD; X0, bphX0, X1, bphX1; X2, bphX2; and X3, bphX3.
Genes 10 00404 g004
Genes 10 00404 g005 550
Figure 5. Comparison of the bph genes of Cupriavidus basilensis KF708 and Comamonas testosteroni KF712 belonging to group IV (type IV) to those of Acidovorax sp. KKS102 and Cupriavidus oxalacticus A5. The identities (88–100%) between the gene clusters are shown by shading in pink to red as indicated at the right bottom of the figure. E, bphE; F, bphF; G, bphG; S, bphS; V, bphV; W, bphW; IS, insertion sequence; tnp, transposase. The bphE, bphF, and bphG of KKS102/Tn4371 are homologous genes of bphX1, bphX3, and bphX2 of KF708/KF712, respectively. The other signs are defined in the legend in Figure 4.
Figure 5. Comparison of the bph genes of Cupriavidus basilensis KF708 and Comamonas testosteroni KF712 belonging to group IV (type IV) to those of Acidovorax sp. KKS102 and Cupriavidus oxalacticus A5. The identities (88–100%) between the gene clusters are shown by shading in pink to red as indicated at the right bottom of the figure. E, bphE; F, bphF; G, bphG; S, bphS; V, bphV; W, bphW; IS, insertion sequence; tnp, transposase. The bphE, bphF, and bphG of KKS102/Tn4371 are homologous genes of bphX1, bphX3, and bphX2 of KF708/KF712, respectively. The other signs are defined in the legend in Figure 4.
Genes 10 00404 g005
Genes 10 00404 g006 550
Figure 6. Organization of the ICEbph-sal in KF strains of group I and KF701 of group II. ICEbph-salKF701, ICEbph-salKF702, ICEbph-salKF703, ICEbph-salKF707, and ICEbph-salKF710 carry the int gene, bph genes, sal genes, and bza genes. ICEbph-salKF716 carries the int gene, bph, and sal genes, but not the bza genes. The sal genes and bza genes in ICEbph-salKF702 are recombined. 1, tRNA-Gly(CCC) genes (partial); 2, int genes; 3, VirB4 components of the type IV secretory pathway; 4, VirD4 component of the type IV secretory pathway; 5, parB genes; 6, parA genes; and 7, putative multidrug efflux pumps. The other, undefined gene components are not shown in the figure.
Figure 6. Organization of the ICEbph-sal in KF strains of group I and KF701 of group II. ICEbph-salKF701, ICEbph-salKF702, ICEbph-salKF703, ICEbph-salKF707, and ICEbph-salKF710 carry the int gene, bph genes, sal genes, and bza genes. ICEbph-salKF716 carries the int gene, bph, and sal genes, but not the bza genes. The sal genes and bza genes in ICEbph-salKF702 are recombined. 1, tRNA-Gly(CCC) genes (partial); 2, int genes; 3, VirB4 components of the type IV secretory pathway; 4, VirD4 component of the type IV secretory pathway; 5, parB genes; 6, parA genes; and 7, putative multidrug efflux pumps. The other, undefined gene components are not shown in the figure.
Genes 10 00404 g006
Genes 10 00404 g007 550
Figure 7. Integration sites of ICEbph-sal in the KF strains of group I and KF701 of group II. The attL (18 bp) site, including 18 bp of the 3’ end of tRNA-Gly gene, and attR (18 bp) sites are indicated.
Figure 7. Integration sites of ICEbph-sal in the KF strains of group I and KF701 of group II. The attL (18 bp) site, including 18 bp of the 3’ end of tRNA-Gly gene, and attR (18 bp) sites are indicated.
Genes 10 00404 g007
Genes 10 00404 g008 550
Figure 8. Organization of ICEbphKF708 and ICEbphKF712 in comparison to ICEKKS1024677 and Tn4371. attL1 and attR1 represent integration sites for the ICEbphKF708 and ICEKKS1024677. attL2 and attR2 represent integration sites for the Tn4371 or ICEbphKF712. 1, integrase (int) gene; 2, putative arsenate-resistant gene; 3, parA gene; 4, VirD2 component (relaxase); 5, transcriptional regulator bphS gene; 6, traR gene; and 7, traG gene. The shading in pink to red shows the identities (68–100%) of the gene clusters as indicated at the right bottom of the figure.
Figure 8. Organization of ICEbphKF708 and ICEbphKF712 in comparison to ICEKKS1024677 and Tn4371. attL1 and attR1 represent integration sites for the ICEbphKF708 and ICEKKS1024677. attL2 and attR2 represent integration sites for the Tn4371 or ICEbphKF712. 1, integrase (int) gene; 2, putative arsenate-resistant gene; 3, parA gene; 4, VirD2 component (relaxase); 5, transcriptional regulator bphS gene; 6, traR gene; and 7, traG gene. The shading in pink to red shows the identities (68–100%) of the gene clusters as indicated at the right bottom of the figure.
Genes 10 00404 g008
Genes 10 00404 g009 550
Figure 9. Conjugative transfer of the bph genes from P. putida AC30Bph+ and P. putida KT2440Bph+ into P. putida F39/D. Genomic DNA was digested with SpeI, applied to pulse field gel electrophoresis, and subjected to Southern blot analysis. The bphA1 DNA of KF715 was used as a probe. Lanes 1 and 6, F39/D; lanes 2 and 7, F39/DBph+ transconjugant from AC30Bph+; Lanes 3 and 8, AC30Bph+; Lanes 4 and 9, F39/DBph+ transconjugant from KT2440Bph+; and Lanes 5 and 10, KT2440Bph+.
Figure 9. Conjugative transfer of the bph genes from P. putida AC30Bph+ and P. putida KT2440Bph+ into P. putida F39/D. Genomic DNA was digested with SpeI, applied to pulse field gel electrophoresis, and subjected to Southern blot analysis. The bphA1 DNA of KF715 was used as a probe. Lanes 1 and 6, F39/D; lanes 2 and 7, F39/DBph+ transconjugant from AC30Bph+; Lanes 3 and 8, AC30Bph+; Lanes 4 and 9, F39/DBph+ transconjugant from KT2440Bph+; and Lanes 5 and 10, KT2440Bph+.
Genes 10 00404 g009
Table
Table 1. Biphenyl/PCB degrading KF strains used in this study.
Table 1. Biphenyl/PCB degrading KF strains used in this study.
StrainNBRC Number DDBJ/EMBL/GenBank Accession NumberReferences
Pseudomonas abietaniphila KF701110664B BQJ01000001-BBQJ01000140[15,17]
Pseudomonas aeruginosa KF702110665B BQK01000001-BBQK01000091[15,24]
Pseudomonas putida KF703110666BBQL01000001-BBQL01000135[15,18]
Pseudomonas furukawaii KF707110670AP014862[15,19,25]
Cupriavidus basilensis KF708110671B BQM01000001-BBQM01000062[15,20]
Cupriavidus pauculus KF709110672BBQN01000001-BBQN01000227[15,26]
Pseudomonas toyotomiensis KF710110674BBQO01000001-BBQO01000029[15,27]
Comamonas testosteroni KF712110673BBQP01000001-BBQP01000097[15,21]
Pseudomonas putida KF715110667AP015029, AP015030-AP015033 [15,16,22,28]
Pseudomonas stutzeri KF716110668BBQQ01000001-BBQQ01000030[29]

Share and Cite

MDPI and ACS Style

Hirose, J.; Fujihara, H.; Watanabe, T.; Kimura, N.; Suenaga, H.; Futagami, T.; Goto, M.; Suyama, A.; Furukawa, K. Biphenyl/PCB Degrading bph Genes of Ten Bacterial Strains Isolated from Biphenyl-Contaminated Soil in Kitakyushu, Japan: Comparative and Dynamic Features as Integrative Conjugative Elements (ICEs). Genes 2019, 10, 404. https://doi.org/10.3390/genes10050404

AMA Style

Hirose J, Fujihara H, Watanabe T, Kimura N, Suenaga H, Futagami T, Goto M, Suyama A, Furukawa K. Biphenyl/PCB Degrading bph Genes of Ten Bacterial Strains Isolated from Biphenyl-Contaminated Soil in Kitakyushu, Japan: Comparative and Dynamic Features as Integrative Conjugative Elements (ICEs). Genes. 2019; 10(5):404. https://doi.org/10.3390/genes10050404

Chicago/Turabian Style

Hirose, Jun, Hidehiko Fujihara, Takahito Watanabe, Nobutada Kimura, Hikaru Suenaga, Taiki Futagami, Masatoshi Goto, Akiko Suyama, and Kensuke Furukawa. 2019. "Biphenyl/PCB Degrading bph Genes of Ten Bacterial Strains Isolated from Biphenyl-Contaminated Soil in Kitakyushu, Japan: Comparative and Dynamic Features as Integrative Conjugative Elements (ICEs)" Genes 10, no. 5: 404. https://doi.org/10.3390/genes10050404

APA Style

Hirose, J., Fujihara, H., Watanabe, T., Kimura, N., Suenaga, H., Futagami, T., Goto, M., Suyama, A., & Furukawa, K. (2019). Biphenyl/PCB Degrading bph Genes of Ten Bacterial Strains Isolated from Biphenyl-Contaminated Soil in Kitakyushu, Japan: Comparative and Dynamic Features as Integrative Conjugative Elements (ICEs). Genes, 10(5), 404. https://doi.org/10.3390/genes10050404

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop