Molecular Genetic Characterization of an Anthrabenzoxocinones Gene Cluster in Streptomyces Sp. FJS31-2 for the Biosynthesis of BE-24566B and Zunyimycin Ale

Genome mining is an effective tool used to discover novel natural products from actinomycetes. Genome sequence analysis of Streptomyces sp. FJS31-2 revealed the presence of one putative type II polyketide gene cluster (ABX), which may correspond to type II polyketide products including BE-24566B and its chloro-derivatives. The addition of natural humus acid successfully activated the biosynthsis of the abx gene cluster. BE-24566B and its chloro-derivatives, named zunyimycin A, were also detected. The targeted deletion of the polyketide skeleton synthesis genes such as abxp, abxk, and abxs was performed in the wild strain to identify the gene cluster for BE-24566B biosynthesis.


Comparative Analysis of ABX Gene Clusters
The draft genome sequence of Streptomyces sp. FJS31-2 contains 11595310 bp, with a GC content of 70.73%. A total of 10877 putative protein coding sequences (CDSs) were identified in the draft genome covering 8840283 bp, accounting for 76.2% of the genome. Prediction of biosynthetic gene clusters on the whole genome level deduced by antiSMASH online reveals that the genome contains 23 secondary metabolism gene clusters such as polyketide synthase systems of every formally classified families (types I and II PKS), nonribosomal peptide synthetases (NRPS), PKS-NRPS hybrids, terpenes, siderophores, lantibiotics, aminocyclitol, melanin, ectoine, bacteriocin, and so on. As shown in Figure 2, genome sequencing and comparative analysis revealed that Streptomyces sp. FJS31-2 contains only one putative type II polyketide gene cluster (GenBank accession No. KU243130). This gene cluster shares a high similarity to its homologous gene cluster from Streptomyces sp. FXJ1.264 at the amino acid sequence (from 86.3% to 94.5%), gene content, and gene order levels (Table 1). Thirteen open reading frames are involved in the ABX polyketide backbone biosynthesis (abxp, abxk, and abxs), regulation (abxR2, abxR1, and abxR), drug resistance (abxA), cyclization (abxC and abxD), methylation (abxM), and halogenation (abxO and abxH) within the gene cluster of Streptomyces sp. 31-2 ( Figure 2A). Ten homologous genes were found in Streptomyces sp. 1.264 ( Figure 2B).

Comparative Analysis of ABX Gene Clusters
The draft genome sequence of Streptomyces sp. FJS31-2 contains 11595310 bp, with a GC content of 70.73%. A total of 10877 putative protein coding sequences (CDSs) were identified in the draft genome covering 8840283 bp, accounting for 76.2% of the genome. Prediction of biosynthetic gene clusters on the whole genome level deduced by antiSMASH online reveals that the genome contains 23 secondary metabolism gene clusters such as polyketide synthase systems of every formally classified families (types I and II PKS), nonribosomal peptide synthetases (NRPS), PKS-NRPS hybrids, terpenes, siderophores, lantibiotics, aminocyclitol, melanin, ectoine, bacteriocin, and so on. As shown in Figure 2, genome sequencing and comparative analysis revealed that Streptomyces sp. FJS31-2 contains only one putative type II polyketide gene cluster (GenBank accession No. KU243130). This gene cluster shares a high similarity to its homologous gene cluster from Streptomyces sp. FXJ1.264 at the amino acid sequence (from 86.3% to 94.5%), gene content, and gene order levels (Table 1). Thirteen open reading frames are involved in the ABX polyketide backbone biosynthesis (abxp, abxk, and abxs), regulation (abxR2, abxR1, and abxR), drug resistance (abxA), cyclization (abxC and abxD), methylation (abxM), and halogenation (abxO and abxH) within the gene cluster of Streptomyces sp. 31-2 ( Figure 2A). Ten homologous genes were found in Streptomyces sp. 1.264 ( Figure 2B). gene cluster shares a high similarity to its homologous gene cluster from Streptomyces sp. FXJ1.264 at the amino acid sequence (from 86.3% to 94.5%), gene content, and gene order levels ( Table 1). Thirteen open reading frames are involved in the ABX polyketide backbone biosynthesis (abxp, abxk, and abxs), regulation (abxR2, abxR1, and abxR), drug resistance (abxA), cyclization (abxC and abxD), methylation (abxM), and halogenation (abxO and abxH) within the gene cluster of Streptomyces sp. 31-2 ( Figure 2A). Ten homologous genes were found in Streptomyces sp. 1.264 ( Figure 2B).

Characterization of the ABX Gene Cluster
Results of PCR products sequencing of the three open reading fragments, namely, abxp, abxk and abxs, from the knockout mutant strain ( abxpks) indicated that the 2.7 kp DNA fragment was replaced by the kanamycin resistance gene ( Figure 3).

Characterization of the ABX Gene Cluster
Results of PCR products sequencing of the three open reading fragments, namely, abxp, abxk and abxs, from the knockout mutant strain (△abxpks) indicated that the 2.7 kp DNA fragment was replaced by the kanamycin resistance gene ( Figure 3). Wild-type Streptomyces sp. FJS 31-2 and its abxpks-deleted mutant strain (△abxpks) were fermented, isolated, and chemically prepared under the same conditions, followed by high-pressure liquid chromatography (HPLC) detection. The result of the metabolite profiling analysis showed that the mutant ceased to produce BE-24566B. However, BE-24566B was still detected in the wild-type strain; this finding confirms that abxp, abxk, and abxs are involved in the biosynthesis of BE-2456B ( Figure 4).  Wild-type Streptomyces sp. FJS 31-2 and its abxpks-deleted mutant strain ( abxpks) were fermented, isolated, and chemically prepared under the same conditions, followed by high-pressure liquid chromatography (HPLC) detection. The result of the metabolite profiling analysis showed that the mutant ceased to produce BE-24566B. However, BE-24566B was still detected in the wild-type strain; this finding confirms that abxp, abxk, and abxs are involved in the biosynthesis of BE-2456B ( Figure 4).

Characterization of the ABX Gene Cluster
Results of PCR products sequencing of the three open reading fragments, namely, abxp, abxk and abxs, from the knockout mutant strain (△abxpks) indicated that the 2.7 kp DNA fragment was replaced by the kanamycin resistance gene ( Figure 3). Wild-type Streptomyces sp. FJS 31-2 and its abxpks-deleted mutant strain (△abxpks) were fermented, isolated, and chemically prepared under the same conditions, followed by high-pressure liquid chromatography (HPLC) detection. The result of the metabolite profiling analysis showed that the mutant ceased to produce BE-24566B. However, BE-24566B was still detected in the wild-type strain; this finding confirms that abxp, abxk, and abxs are involved in the biosynthesis of BE-2456B ( Figure 4).

Chemical Identification of BE-24566B and Zunyimycin A
Chromatography and preparative high-pressure liquid chromatography (HPLC) yielded two compounds from strain Streptomyces sp. FJS 31-2. The NMR data of the two compounds were obtained in MeOD solvent. Compound 1 was obtained as a pale yellow powder. The 13  The HMBC experiment showed a correlation between the proton at δH 7.01 and carbons at δC 39.6, 39.5, 109.9, and 122.1, which indicated that one hydroxyl atom connected with the C-8 position. The signal at δH 6.45 (1H, s), the carbons at δC 164.5, 112.6, and 164.5, and, combined with the ROSEY spectra, the δH 6.45 (1H, s) had no correlation, which implied that the signal at δH 6.45 connected with a benzene ring at position C-12. The proton signal at δH 6.27 (1H, s) correlated with δC 152.6, 150.6, and 115.4, which indicated that the signal at δH 6.27 connected with a benzene ring at position C-4. The proton signal at δH 6.24 (1H, s) had a key correlation with the signal at δC 98.2, which indicated that the signal at δH 6.24 connected with a benzene ring at position C-16. According to the above information, the relative structure of 2 was substituted by chlorine atoms at the C-2 and C-10 positions of BE-24566B; therefore, Compound 2 was named zunyimycin A ( Figure 5).

Chemical Identification of BE-24566B and Zunyimycin A
Chromatography and preparative high-pressure liquid chromatography (HPLC) yielded two compounds from strain Streptomyces sp. FJS 31-2. The NMR data of the two compounds were obtained in MeOD solvent. Compound 1 was obtained as a pale yellow powder. The 13  The HMBC experiment showed a correlation between the proton at δH 7.01 and carbons at δC 39.6, 39.5, 109.9, and 122.1, which indicated that one hydroxyl atom connected with the C-8 position. The signal at δH 6.45 (1H, s), the carbons at δC 164.5, 112.6, and 164.5, and, combined with the ROSEY spectra, the δH 6.45 (1H, s) had no correlation, which implied that the signal at δH 6.45 connected with a benzene ring at position C-12. The proton signal at δH 6.27 (1H, s) correlated with δC 152.6, 150.6, and 115.4, which indicated that the signal at δH 6.27 connected with a benzene ring at position C-4. The proton signal at δH 6.24 (1H, s) had a key correlation with the signal at δC 98.2, which indicated that the signal at δH 6.24 connected with a benzene ring at position C-16. According to the above information, the relative structure of 2 was substituted by chlorine atoms at the C-2 and C-10 positions of BE-24566B; therefore, Compound 2 was named zunyimycin A ( Figure 5).

Discussion
Cluster analysis results show that AbxH is clustered with the identified halogenases involved in the halogenated modification of pyrrole or its derivatives. Furthermore, AbxH exhibits the highest similarity (i.e., 49% at the amino acid sequence level) to NapH2 (Figure 6), which is involved in the halogenation reaction of napyradiomycin [29]. These results suggest that AbxH may possess a capacity for the halogenation reaction of polyketide. No abxH homologous gene was predicted in the genome, and no halogenated natural product was detected from the fermented products of

Discussion
Cluster analysis results show that AbxH is clustered with the identified halogenases involved in the halogenated modification of pyrrole or its derivatives. Furthermore, AbxH exhibits the highest similarity (i.e., 49% at the amino acid sequence level) to NapH2 (Figure 6), which is involved in the halogenation reaction of napyradiomycin [29]. These results suggest that AbxH may possess a capacity for the halogenation reaction of polyketide. No abxH homologous gene was predicted in the genome, Streptomyces sp. 1.264; hence, abxH (halogenase) may play a vital role in zunyimycin biosynthesis (Scheme 1).

Strains, Plasmids, and Reagents
Streptomyces sp. FJS31-2 was isolated from a soil sample collected from the Fanjing Mountain of the Guizhou Province at an altitude of 800 m in October 2012. The strain was deposited in the China General Microbiological Culture Collection Center under accession number CGMCC 4.7321; the strain exhibits biological activities against Candida albicans, Bacillus subtilis, Staphylococcus aureus, and Micrococcus luteus. The nucleotide sequencing was performed by BGI (Shenzhen, China) using an Illumina/Hiseq 2000 sequencer. Assembly was performed using SOAP de novo software [30]. Genomic DNA was directly extracted from the expanded culture of the ISP 2 liquid medium and used for genome sequencing and genomic library construction [31]. The bacterial strain EPI300-T1 R and the plasmid pCC2FOS were obtained from Epicentre (Madison, WI, USA) and used to construct a genomic library by using CopyControl™ HTP Fosmid Library Production Kit (Cat. No. CCFOS059; Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Cloning plasmids, namely, pMD18-T and pSIMPLE-18 EcoR V/BAP, were purchased from Takara Biotechnology Company (Dalian, China) and used to construct a sub-clone library. Restriction enzymes, T4 DNA ligase, DNA end blunting Streptomyces sp. 1.264; hence, abxH (halogenase) may play a vital role in zunyimycin biosynthesis (Scheme 1).

Strains, Plasmids, and Reagents
Streptomyces sp. FJS31-2 was isolated from a soil sample collected from the Fanjing Mountain of the Guizhou Province at an altitude of 800 m in October 2012. The strain was deposited in the China General Microbiological Culture Collection Center under accession number CGMCC 4.7321; the strain exhibits biological activities against Candida albicans, Bacillus subtilis, Staphylococcus aureus, and Micrococcus luteus. The nucleotide sequencing was performed by BGI (Shenzhen, China) using an Illumina/Hiseq 2000 sequencer. Assembly was performed using SOAP de novo software [30]. Genomic DNA was directly extracted from the expanded culture of the ISP 2 liquid medium and used for genome sequencing and genomic library construction [31]. The bacterial strain EPI300-T1 R and the plasmid pCC2FOS were obtained from Epicentre (Madison, WI, USA) and used to construct a genomic library by using CopyControl™ HTP Fosmid Library Production Kit (Cat. No. CCFOS059; Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Cloning plasmids, namely, pMD18-T and pSIMPLE-18 EcoR V/BAP, were purchased from Takara Biotechnology Company (Dalian, China) and used to construct a sub-clone library. Restriction enzymes, T4 DNA ligase, DNA end blunting Scheme 1. Proposed biosynthetic pathway of zunyimycin A.

Strains, Plasmids, and Reagents
Streptomyces sp. FJS31-2 was isolated from a soil sample collected from the Fanjing Mountain of the Guizhou Province at an altitude of 800 m in October 2012. The strain was deposited in the China General Microbiological Culture Collection Center under accession number CGMCC 4.7321; the strain exhibits biological activities against Candida albicans, Bacillus subtilis, Staphylococcus aureus, and Micrococcus luteus. The nucleotide sequencing was performed by BGI (Shenzhen, China) using an Illumina/Hiseq 2000 sequencer. Assembly was performed using SOAP de novo software [30]. Genomic DNA was directly extracted from the expanded culture of the ISP 2 liquid medium and used for genome sequencing and genomic library construction [31]. The bacterial strain EPI300-T1 R and the plasmid pCC2FOS were obtained from Epicentre (Madison, WI, USA) and used to construct a genomic library by using CopyControl™ HTP Fosmid Library Production Kit (Cat. No. CCFOS059; Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Cloning plasmids, namely, pMD18-T and pSIMPLE-18 EcoR V/BAP, were purchased from Takara Biotechnology Company (Dalian, China) and used to construct a sub-clone library. Restriction enzymes, T4 DNA ligase, DNA end blunting kination ligation, and PCR kits were also purchased from the same company. The primers used for gene cloning and DNA sequencing were synthesized by Beijing Invitrogen Biotechnology Co., Ltd. (Beijing, China). The DNA gel extraction kit was obtained from Axygen Biosciences, Inc. (Hangzhou, China).

DNA Isolation, Manipulation, and Sequencing
Genomic DNA was isolated from Streptomyces sp. FJS31-2 by using the standard method described by Hopwood [32]. Molecular manipulation was conducted according to the standard methods described by Sambrook and Russel [33]. Plasmid and fosmid sequencing was performed on the 3730xl DNA analyzer at Beijing Invitrogen Biotechnology Co., Ltd. by using the shotgun cloning strategy. The pCC2/pEpiFOS forward and reverse primers (sequences in the CopyControl protocol) were used to sequence the fosmid ends. Moreover, the M13 primers were used to sequence the subclone ends. Draft genome sequencing of Streptomyces sp. FJS31-2 was performed on Illumina HiSeq2000 (BGI Biotechnology Co., Ltd. Shenzhen, China) by using a proprietary reversible terminator-based strategy.

The Genomic Library and Subclone Library Construction and Screening
The full-length sequence of the predicted lobophorin biosynthetic gene cluster was obtained by constructing the genomic library of Streptomyces sp. FJS31-2 in CopyControl pCC2FOS according to the manufacturer 1 s instructions. The library contains more than 6000 clones with approximately 40 kb inserted genomic DNA fragments. The genomic library was screened through colony PCR with primer pairs designed based on the nucleic acid sequences of abxH according to the results of the predictive analysis of the draft genome sequencing. The end of the positive clone was sequenced to confirm that the clone was located in the abx cluster. Next-round screening was performed with primer pairs designed based on the end sequences of the last round of screening results until three positive clones overlapping one another were obtained. Accordingly, 1 µg of fosmid DNA was fragmented using an ultrasonic cell disruptor under 200 W for 1 s on ice. DNA fragments within 2-4 kb size were recovered using the AxyPrep DNA gel extraction kit. The fragments were cloned into pSIMPLE-18 EcoR V/BAP vector to construct subclone libraries after being end blunted and 5 1 -phosphorylated. PCR screening was performed using sequence-specific primers designed based on the end sequences of each gap. The positive clones were then sequenced with M13-47F/M13-48R and specific primers.

Analysis of the Anthrabenzoxocinone Biosynthetic Gene Cluster
A 16.5-kb contiguous DNA segment was obtained by subcloned library sequencing and deposited in GenBank under accession number KU243130. The sequence was uploaded to the antiSMASH website [34] and used to determine gene content and order. The open reading frames proposed for the biosynthesis pathway of the anthrabenzoxocinone skeleton, post modification, and regulation were deduced based on functional gene analysis. The analysis was performed by comparing the sequence homology with its homologous clusters by using BlastP. These clusters were found by antiSMASH. The comparative analysis and the global rearrangement structure of the two anthrabenzoxocinone biosynthetic gene clusters were conducted using Mauve software [35].

Inactivation of the Anthrabenzoxocinone Biosynthetic Gene Cluster
A 2.7-kb DNA fragment (from 6165 nt to 8872 nt of KU243130) encoding the ABX polyketide backbone biosynthesis (abxp, abxk, and abxs) was replaced by the kanamycin resistance protein-coding region (AAF85969.1). This region was constructed by PCR, inserted into the pKC1139 plasmid, and transformed into E. coli ET12567. The constructs for gene inactivation were introduced into Streptomyces sp. FJS31-2 by conjugal transfer following the procedure of Hopwood 2 . Disruption was confirmed by PCR analysis. HPLC (LC20AT, Shimadzu, Beijing, China) was used to confirm the abolishment of anthrabenzoxocinone production.

Fermentation, Isolation, and Chemical Identification of BE-24566B and Zunyimycin A
Streptomyces sp. FJS31-2 was cultured using 140ˆ500 mL shake flasks containing 100 mL of ISP 2 agar medium with 10% natural humus acid water extracts and then incubated for 15 d at 28˝C. The solid culture was mashed and extracted three times with 140 L of ethanol after cultivation. The organic portion was then concentrated in vacuo to remove the solvent. The crude extract was applied to silica gel column chromatography using the CHCl 3 /MeOH gradient to obtain the crude products. Further purification was conducted using Sephadex LH-20 (GE Healthcare, Tokyo, Japan) (MeOH) column and RP-HPLC (Shimadzu SPD-M20A with Xbridge ODS 10 mmˆ150 mm column). Compounds were identified with an HRESI-MS (Waters Xevo G2 QTOF mass spectrometer (Waters corportion, Milford, MA, USA) and NMR (Bruker AV 600 MHz) (Bruker Corporation, Karlsruhe, Germany) analyses.

Conclusions
In summary, a cryptic anthrabenzoxocinone gene cluster was identified by genome mining. This cluster was successfully activated by adding natural humus acid. Novel halogenated natural products were isolated, characterized, and named as zunyimycin A. Furthermore, the results indicated that halogenase AbxH is responsible for the halogenation reaction of zunyimycin A and may provide potential for the halogenation of natural products.