- Open Access
Identification of nocamycin biosynthetic gene cluster from Saccharothrix syringae NRRL B-16468 and generation of new nocamycin derivatives by manipulating gene cluster
Microbial Cell Factoriesvolume 16, Article number: 100 (2017)
Nocamycins I and II, produced by the rare actinomycete Saccharothrix syringae, belong to the tetramic acid family natural products. Nocamycins show potent antimicrobial activity and they hold great potential for antibacterial agent design. However, up to now, little is known about the exact biosynthetic mechanism of nocamycin.
In this report, we identified the gene cluster responsible for nocamycin biosynthesis from S. syringae and generated new nocamycin derivatives by manipulating its gene cluster. The biosynthetic gene cluster for nocamycin contains a 61 kb DNA locus, consisting of 21 open reading frames (ORFs). Five type I polyketide synthases (NcmAI, NcmAII, NcmAIII, NcmAIV, NcmAV) and a non-ribosomal peptide synthetase (NcmB) are proposed to be involved in synthesis of the backbone structure, a Dieckmann cyclase NcmC catalyze the releasing of linear chain and the formation of tetramic acid moiety, five enzymes (NcmEDGOP) are related to post-tailoring steps, and five enzymes (NcmNJKIM) function as regulators. Targeted inactivation of ncmB led to nocamycin production being completely abolished, which demonstrates that this gene cluster is involved in nocamycin biosynthesis. To generate new nocamycin derivatives, the gene ncmG, encoding for a cytochrome P450 oxidase, was inactivated. Two new nocamycin derivatives nocamycin III and nocamycin IV were isolated from the ncmG deletion mutant strain and their structures were elucidated by spectroscopic data analyses. Based on bioinformatics analysis and new derivatives isolated from gene inactivation mutant strains, a biosynthetic pathway of nocamycins was proposed.
These findings provide the basis for further understanding of nocamycin biosynthetic mechanism, and set the stage to rationally engineer new nocamycin derivatives via combinatorial biosynthesis strategy.
Nocamycins I and II (Fig. 1), isolated from the broth of Saccharothrix syringae NRRL B-16468 by Russian scientists in 1977, belong to the tetramic acid (2, 4-pyrrolidinedione) family natural products [1,2,3]. The original structural assignment of nocamycin I was incorrect and it was revised by a Japanese research group . The Japanese research group reported two compounds Bu-2313A and Bu-2313B from the strain Microtetraspora caesia ATCC 31295 nearly at the same time . Further structural elucidation showed that Bu-2313B was virtually identical to nocamycin I . Beyond the common tetramic acid structure, a tricyclic ketal structure is another interesting motif in nocamycins. In terms of structural viewpoint, streptolydigin, tirandamycins and tirandalydigin are closely related to nocamycins (Fig. 1). Among these compounds, nocamycins I and II are unique because they have a fused oxolane ring system other than an oxirane (spiro or fused) ring in streptolydigin, tirandamycin and tirandalydigin.
Nocamycin I (Bu-2313B) displays broad antimicrobial activity toward a panel of Gram-positive and Gram-negative anaerobic bacteria as well as some aerobic bacteria. Inhibitions of anaerobic bacteria Bacteroides fragilis, Clostridium sp., Fusobacterium sp., Sphaerophorus sp. by nocamycins are particularly potent, and the minimum inhibitory concentrations (MICs) are in the range of 0.1–0.4 μg/mL [5,6,7]. Further in vivo experiments conducted in mice showed that nocamycin I is effective in protecting mice against B. fragilis A20928-1 and Clostridium perfingens A9635 when administered by both oral and subcutaneous routes . In addition, nocamycins show antitumor effects . Up to now, the exact antibacterial mold of action of nocamycins has not been investigated. The closely related compounds tirandamycin and streptolydigin are validated to be inhibitors of bacterial RNA polymerase (RNAP), thus nocamycins are probably to be inhibitors of RNAP. In recent years, the molecular evidences for the structural basis of the RNAP interaction mechanism of this class of natural products have been disclosed by co-crystal complexes of streptolydigin with RNAPs from Escherichia coli and Thermus thermophilus [8, 9]. The key affinities of both bicyclic ketal and tetramic acid structures with RNAPs have been observed from the co-crystal complexes, indicating the substitution or modification in these two structural motifs is critical for the biological activity [8, 9]. Meanwhile, results of antibacterial activities of streptolydigin, tirandamycin and their congeners also demonstrated that the two featured motifs are closely related to the activity of this family of natural products [10,11,12].
The intriguing structure, action mold and biological activity of this small class of natural products attract more and more attentions from biochemists. So far, the gene clusters responsible for tirandamycin and streptolydigin biosynthesis have been identified from three different Streptomyces species by Sherman, Salas and Ju group, respectively [10, 13, 14]. Both tirandamycins and streptolydigin are assembled by hybrid iterative type I polyketide synthase (PKS) and non-ribosomal peptide synthase (NRPS). The functions of a number of genes related to post-tailoring, regulator and resistance involved in tirandamycin and streptolydigin biosynthetic pathway have been fully elucidated [12, 13, 15,16,17,18,19,20]. Some streptolydigin derivatives were generated by using combinatorial biosynthesis method . In streptolydigin and tirandamycins biosynthetic pathway, a uniform strategy is employed to catalyze the formation of tetramic acid moiety . The mechanism of bicyclic ketal structure formation remains unclear since no related gene candidates have been discovered in the two gene clusters. To fully understand the biosynthetic pathway of nocamycins, provide insights into the formation of bicyclic ketal structure and generate diversified nocamycin derivatives, we started to identify the nocamycin biosynthetic gene cluster from S. syringae NRRL B-16468. Here, we report the identification of nocamycin biosynthetic gene cluster and new nocamycin derivatives generated by manipulating the gene cluster.
Bacterial strains, plasmids, medium and culture conditions
The bacteria and plasmids used in this study are listed in Table 1. S. syringae was maintained on ISP4 agar medium. The medium used for fermentation of S. syringae and its mutant strains consists of 1% soybean, 3% glycerol, 0.5% mycose, 0.2% NaCl and 0.2% CaCO3. All cultures for S. syringae were incubated at 28 °C. For E. coli, Luria–Bertani (LB) liquid or agar media were used with appropriate antibiotics at a final concentration of: 100 μg/mL ampicillin (Amp), 50 μg/mL apramycin (Apr), 50 μg/mL kanamycin (Kan), 25 μg/mL chloramphenicol (Cml) and 50 μg/mL trimethoprim (TMP).
DNA sequencing, assembly and analysis
After growing in TSB medium for 48–72 h, the genomic DNA of S. syringae NRRL B-16468 was extracted according to standard protocols . Then, the genomic DNA was shotgun sequenced and annotated by Shanghai South Gene Technology Co. Ltd. (Shanghai, China). The gene cluster responsible for secondary metabolite biosynthesis was analyzed by antiSMASH online analysis tool (http://antismash.secondarymetabolites.org/). DNA and corresponding protein sequences in nocamycin gene cluster were analyzed by ORF finder program (http://www.ncbi.nlm.nih.gov/gorf/gorf.html), Frameplot 2.3.2 program (http://www.nih.go.jp/~jun/cgi-bin/frameplot.pl), and BLAST program (http://blast.ncbi.nlm.nih.gov/).
Construction and screening of S. syringae genomic library
Genomic library of S. syringae NRRL B-16468 was constructed using SuperCos1 Vector Kit according to manufacturer’s instruction (Stratagene). The library was packaged using phage extracts and transduced into the E. coli LE392. About 2600 resulting transductants were picked up and transferred to twenty-seven 96-well microtiter plates containing 150 μL LB medium supplemented with Kan (50 μg/mL). After overnight incubation at 37 °C, 30 μL E. coli broth in every microtiter pore was absorbed and mixed every 12 clones in a horizontal line and every 8 clones in a vertical line for each 96-well plate. Glycerol was added to the remaining broth of the clones (20% final concentration) for permanent stock. The DNA of mixed clones was extracted as templates for PCR screening.
The primer pairs targeted the cytochrome P450 oxidase gene (NcmG-SF and NcmG-SR), Dieckmann cyclase gene (NcmC-SF and NcmC-SR) and DH domain at NcmAII gene (DH-SF and DH-SR) were designed and they were used as PCR primers to screen S. syringae NRRL B-16468 genomic library (Table 2). The positive clones were selected from the genomic library. The cosmids were extracted and further analyzed by terminal-sequencing. The PCR reaction (20 mL volume) contained 2 μL 10 × PCR buffer, 1.6 μL dNTPs (2.5 mM), 0.4 μL forward primer (10 μM), 0.4 μL reverse primer (10 μM), 1 μL dimethylsulfoxide (DMSO), 1 μL DNA template, 0.1 μL rTaq (5 U/μL), and 13.5 μL ddH2O. The following PCR program was used: 94 °C, 4 min, 30 cycles of 94 °C, 45 s, 59 °C, 45 s, 72 °C, 1 min, and a final extension cycle at 72 °C, 10 min. Eventually, two cosmids p5-C-9 and p2-H-12 were chosen for further gene-inactivation experiments.
Generation of S. syringae mutant strains
λ-RED recombination technology was employed to inactivate the target gene ncmB, ncmL and ncmG according to literature previously reported . The primer pairs used for PCR-targeting are listed in Table 2. The fragment oriT/acc(3)IV cassette was used to replace partial gene region of ncmB or ncmL in p5-C-9 to generate cosmid pMoS1001 (ΔncmB) or pMoS1002 (ΔncmL). For ncmG, partial gene region was replaced by fragment oriT/acc(3)IV cassette in cosmid p2-H-12 and plasmid pMoS1003(ΔncmG) were generated. After verified by PCR and restriction enzyme digestion analysis, the correct mutated cosmids were introduced into E. coli ET12567/pUZ8002 and conjugated with wild type S. syringae spores. The wild type S. syringae spores were germinated in LB medium for 4–5 h at 30 °C, 200 rpm. The E. coli ET12567/pUZ8002 containing each mutated cosmid was grown in LB medium supplemented with Kan (50 μg/mL), Amp (100 μg/mL), Cml (25 μg/mL) and Apr (50 μg/mL) to OD600 = 0.6–0.8. Then the cells were harvested, washed twice with LB medium, mixed with germinated wild type spores and plated on ISP4 medium. The plates were incubated at 30 °C for 24 h. Then, each plate was covered by 800 μL sterile water supplemented with 30 μL TMP (50 mg/mL) and 30 μL Apr (50 mg/mL). The plates were continued incubated at 30 °C for 7–10 days until exconjugants appeared. Double cross-over mutants were first selected by the phenotype of Kan sensitive (KanS) and Apr resistant (AprR), and the genotype of the mutants were further confirmed by PCR. Finally, the mutant strains S. syringae MoS-1001 (ΔncmB), S. syringae MoS-1002 (ΔncmL) and S. syringae MoS-1003 (ΔncmG) were obtained.
Fermentation and analysis of S. syringae and mutant strains
Saccharothrix syringae wild type and mutant strains were inoculated in 250 mL flasks with 50 mL medium and incubated on a rotary shaker at 28 °C, 200 rpm. After 7 days fermentation, each of the 50 mL culture was added with 100 mL ethyl acetate and then vigorously mixed for 30 min. The ethyl acetate phase was evaporated to dryness to yield a residue. The residue was dissolved in 1 mL methanol and centrifuged, then, the supernatant was subjected to HPLC analysis. Analytical HPLC was performed on Agilent 1260 HPLC system (Agilent Technologies Inc., USA) equipped with a binary pump and a diode array detector using a Phenomenex Prodigy ODS column (150 × 4.60 mm, 5 μ) with UV detection at 355 nm. The mobile phase comprises solvent A and B. Solvent A consists of 15% CH3CN in water supplemented with 0.1% trifluoroacetic acid (TFA). Solvent B consists of 85% CH3CN in water supplemented with 0.1% TFA. Samples were eluted with a linear gradient from 5 to 90% solvent B in 20 min, followed by 90–100% solvent B for 5 min, then eluted with 100% solvent B for 3 min, at a flow rate of 1 mL/min and UV detection at 355 nm.
Isolation of new produced nocamycin derivatives from ΔncmG mutant strain
Two-step fermentation was used to culture ΔncmG mutant strain. 250 mL flask containing 50 mL medium was used as seed culture and 500 mL flask containing 100 mL medium was used as fermentation medium. Appropriate spores were inoculated to seed culture and grown at 28 °C, 200 rpm for 3 days. Then, 5 mL seed medium was inoculated to 100 mL fermentation medium and continued 7 days culture. 15 L liquid medium was used in total. After incubation, the culture broth was collected and centrifuged. The supernatant was extracted by ethyl acetate for three times and the mycelium was extracted by acetone for three times. Then, the entire organic phases were evaporated to dryness to yield crude extract. The crude extract was dissolved in a mixture of CH3OH: CHCl3 (1:1) and mixed with appropriate amount of silica gel (100–200 mesh, Qingdao Marine Chemical Corporation, China). The sample was applied on normal phase silica gel column chromatography and eluted with CHCl3-CH3OH (100:0–50:50) to give 10 fractions. All the fractions were analyzed by HPLC. Fraction 4 and 5 contained the major targeted compound nocamycin III and fractions 7 and 8 contained the major targeted compound nocamycin IV. The fractions 4–5 and fractions 7–8 were further purified on reverse phase C-18 silica gel (YMC, Japan) by using medium-pressure liquid chromatography (MPLC, Agela corporation, China) eluted by a linear gradient from 20 to 90% CH3CN in water, respectively. The sub-fractions contained targeted compounds were further purified by Sephadex LH-20 (GE healthcare, Sweden) gel filtration chromatography to afford the purified nocamycin III and nocamycin IV.
Spectroscopy analysis of new produced nocamycin derivatives
1H and 13C NMR spectra of nocamycin III and nocamycin IV were recorded at 25 °C on Bruker AV 500 instruments. HR–ESI–MS spectra data were acquired on a Waters micro MS Q-Tof spectrometer.
Sequencing and identification of nocamycin gene cluster
Saccharothrix syringae NRRL B-16468 genome was shotgun sequenced by Hiseq4000 technologies and the sequence reads were assembled into 10.8 Mb nucleotides. Then, S. syringae NRRL B-16468 genome data was analyzed by using online antiSMASH tool . AntiSMASH analysis results demonstrated that a hybrid PKS-NRPS gene cluster designated as Ncm seems to be the candidate responsible for nocamycin biosynthesis since it shows high similarity to tirandamycin biosynthetic gene cluster. In the Ncm gene cluster, some deduced gene products such as NcmC, NcmE, NcmF and NcmB show high similarity to TrdC, TrdE, TrdF and TrdB originated from tirandamycin biosynthetic pathway, respectively . Thus, we assumed that this cluster is probably involved in nocamycin biosynthesis. We then screened S. syringae genomic library by using PCR method with the primer pairs targeted at ncmG, ncmC and dehydratase (DH) domain at module 4. In total of eight positive cosmids were obtained. The eight cosmids were end-sequenced and two cosmids p2-H-12 and p5-H-9 were used for further gene inactivation experiments. To verify our hypothesis, a gene ncmB encoding a NRPS was inactivated to afford the strain S. syringae MoS-1001 (Additional file 1: Figure S1). HPLC analysis of the extract of S. syringae MoS-1001 fermentation broth revealed that S. syringae MoS-1001 failed to produce nocamycin I and II (Fig. 2I) completely, indicating ncmB’s involvement in nocamycin biosynthesis. This result also demonstrated this PKS-NRPS gene cluster is responsible for nocamycin biosynthesis. On basis of bioinformatics analysis, about 61 kb DNA locus consisted of 21 open reading frames (ORFs) whose deduced products are likely to be involved in nocamycin biosynthesis (Fig. 3; Table 3). Corresponding homologues and deduced function of each ncm gene are listed in Table 3. The sequence data of nocamycin biosynthesis in this study have been deposited in Genbank under accession number KY287782.
Linear chain assembly and releasing
Hybrid PKS-NRPS are employed to construct the backbone structure of nocamycin. Five type I PKS genes ncmAI, ncmAII, ncmAIII, ncmAIV and ncmAV transcribed in the same direction were identified in the gene cluster (Fig. 3). The deduced products of the five PKS genes were constituted by four, one, one, one and two modules respectively to assemble the polyketide backbone (Fig. 4). Each PKS module minimally contains ketosynthase (KS), acyltransferase (AT) and acyl carrier protein (ACP) domains. The conserved motifs from PKS modules in nocamycin gene cluster are listed in Additional file 1: Table S1. Except for loading module, M2 and M8, each module possess a ketoreductase (KR) domain with conserved active motif. KR domain in module M5 is the only KR domain contains the characteristic of A-type KRs, and all the other KR domains display the conserved motif characteristic for the B-type KRs . A characteristic KSQ domain of loading module indicated that a malonoyl–CoA might be used to provide acetate as starter unit, and this phenomenon was observed in tirandamycin and streptolydigin gene clusters [10, 14]. As shown in Table 4 and Fig. 4, the AT domains in extension modules M3, M7 and M8 display conserved active motif specific for malonate-CoA incorporation [29, 30], whereas AT domains in extension modules M1, M2, M4, M5 and M6 show conserved active motif specific for methylmalonate-CoA incorporation [29, 30], which is in good agreement with the polyketide carbon skeleton. There are three DH domains with conserved active motif HXXXGXXXXP distributed in module M4, M6 and M7 .
NcmB, a NRPS, shows most similarity to TrdD (56% identity/66% similarity) from Streptomyces sp. SCSIO1666 involved in tirandamycin biosynthetic pathway . Three domains condensation (C), adenylation (A), and peptidyl carrier protein (PCP) are found in NcmD. The amino acid binding pocket DILQLGVI located in A domain is predicted to activate glycine, which is accord to nocamycin structure.
NcmC shows most similarity to TrdC (45% identity/58% similarity) from Streptomyces sp. SCSIO1666 involved in tirandamycin biosynthetic pathway . TrdC and its analogues SlgC, KirHI have been determined as Dieckmann cyclases, and they catalyze the formation of tetramic acid or pyridone moiety . Bioinformatics analyses revealed that NcmC also possesses the characteristic catalytic traid Cys-Asp-His (Additional file 1: Figure S4). Thus, in nocamycin biosynthesis pathway, NcmC is proposed to be responsible for the PK-NRP chain releasing and catalyze the formation of tetramic acid moiety.
Genes involved in post-tailoring steps
After linear chain released from PKS-NRPS and formation of tetramic acid moiety, several post tailoring processes including oxolane ring system, C-10 hydroxyl/ketone group, C-14 methoxycarbonyl group are required to synthesis nocamycin I. Within the identified gene cluster, there are six genes encoding two cytochrome P450 monooxygenase (ncmO and ncmG), one monooxygenase (ncmL), one carboxylate O-methyltransferase (ncmP), one short chain dehydrogenase (ncmD) and one glycoside hydrolase (ncmE) are likely to be involved in these steps.
The glycoside hydrolase NcmE shows identity to TrdE (60% identity/76% similarity) involved in tirandamycin biosynthesis . In tirandamycin pathway, TrdE functions as a dehydratase and it is responsible for the formation of C11–C12 double bond . Thus, we propose that NcmE is a dehydratase and it catalyzes the formation of C11–C12 double bond.
Both cytochrome P450 monooxygenases NcmG and NcmO possess the highly conserved heme-binding domain (GXXXCXG), K-helix (EEXLL) and oxygen binding region (Additional file 1: Figure S5) [32, 33]. NcmG shows similarity to Mur7 (51% identity/64% similarity) involved in muraymycin biosynthesis . NcmO shows similarity to QnmO (48% identity/63% similarity) involved in quartromicin biosynthesis . Since at least three oxidative tailoring steps, including the formation of tetrahydrofuran fused in bicyclic ketal structure, C-10 hydroxyl and C-14 carboxyl group are required, NcmG or NcmO is proposed to be bifunctional. Sequence alignments analysis revealed that NcmO and NcmG are distinct from the cytochrome P450 oxidases TrdI/TamI, SlgO1 and SlgO2 involved in tirandamycin or streptolydigin biosynthetic pathway [10, 15] (Additional file 1: Figure S6). The reason for this phenomenon may attribute to the different oxidative modification in bicyclic ketal structure of nocamycin, streptolydigin and tirandamycin.
Within nocamycin gene cluster, only ncmP encodes for a SAM-dependent carboxylate O-methyltransferase and it shows identity to NokK (44% identity/55% similarity) and NivG (43% identity/53% similarity). Both NokK and NivG are proposed to catalyze methyl esterification of the carboxylate group in biosynthesis of K-252a and nivetetracyclates, respectively [36, 37]. Hence, it should be evident that NcmP serves as the best candidate responsible for methyl esterification during nocamycin biosynthetic pathway.
NcmL shows similarities to monooxygenase from a series of actinomycetes. BLAST analysis revealed that NcmL displays FAD-binding domain (pfam01494). Unlike bicovalent flavinylation protein TrdL/TamL involved in tirandamycin biosynthetic pathway, NcmL has no conserved His and Cys dual active site residues that distributed in 8α-histidyl and 6-S-cysteinyl FAD linked monooxygenase family (Additional file 1: Figure S7) [15, 16]. To investigate the function of NcmL in nocamycin biosynthesis, the gene ncmL was inactivated (Additional file 1: Figure S2). The fermentation broth of ΔncmL mutant strain was analyzed by HPLC (Fig. 2III). The results revealed that the titer of nocamycin I and nocamycin II in ΔncmL deletion strain is identical to that in wild type, indicating NcmL is not involved in nocamycin biosynthesis.
The putative product of ncmD shows identity to a series of short chain dehydrogenase (SDR) family oxidoreductase originated from various bacteria. NcmD shares the Rossmann fold NAD-binding motif and characteristic NAD-binding and catalytic sequence patterns . NcmD shows closet similarity to BatM (40% identity/56% similarity) which was proposed to catalyze the conversion from hydroxyl to ketone in C-17 position during kalimantacin/batumin-related polyketide antibiotic biosynthesis . Thus, NcmD is proposed to be the candidate to convert hydroxy to ketone in C-10 position.
Genes involved in regulation, resistance and unknown functions
Five genes related to regulation and resistance are easily discerned from the nocamycin biosynthetic gene cluster. NcmN encodes for a LuxR family regulator and it shows similarity to a series of regulators from different actinomycetes, including QmnRg4 (43% identity/55% similarity) from Amycolatopsis orientalis involved in quartromicin biosynthesis and TamH (39% identity/52% similarity) involved in tirandamycin biosynthesis [14, 35]. The characteristic C-terminal helix-turn-helix (HTH) DNA binding domain signature and a N-terminal ATP-binding domain represented by discernible Walker A (GxxGxGK) and Walker B (R/K-X(7-8)-H(4)-D) motifs present in all members of this family of regulatory proteins are found in NcmN . NcmJ is similar to AAA family ATPase from different actinomycetes. AAA family ATPases are present in all kingdoms and they are often involved in DNA replication, repair, recombination and transcription . NcmJ contains the Walker A and Walker B motifs, which is the hallmark of ATP-binding domain in these proteins . NcmI encodes a PadR family transcriptional regulator and it shows similarities to several PadR-like proteins of unknown function from different actinomycetes. PadR-like proteins is a quite recently identified family of regulatory proteins, named after the phenolic acid decarboxylation repressor of Bacillus subtilis [42, 43]. The hallmark of this family transcriptional regulator is a highly conserved N-terminal winged helix-turn-helix (HTH) domain with about 80–90 residues [44, 45], which is also found in NcmI. NcmK encodes for a TetR family transcriptional regulator and it shows identity to TrdK (49% identity/64% similarity) involved in tirandamycin biosynthesis . The characteristic N-terminal helix-turn-helix (HTH) DNA binding domain signature (pfam00440) presented in all members of this family of regulatory proteins has been found in NcmK.
NcmH, a major facilitator superfamily (MFS) transporter, shows identity to ChaT1 (42% identity/60% similarity) from Streptomyces chartreusis involved in antitumor agent chartreuse in biosynthesis pathway, is a candidate protein for resistance . NcmQ is similar to the proteins belong to glyoxalase/bleomycin resistance protein/dioxygenase superfamily. The exact role of NcmQ in nocamycin biosynthesis is unclear and we assume that NcmQ is likely involved in resistance.
The deduced product of ncmF shows similarity to a series of prenyltransferase, including TrdF (50% identity/61% similarity) involved in tirandamycin biosynthesis and SlgF (51% identity/60% similarity) involved in streptolydigin biosynthesis, respectively [10, 14]. Previous studies on TrdF and SlgF demonstrated that both the proteins show no relationship with tirandamycin or streptolydigin biosynthesis [10, 14]. Thus, we hypothesize that NcmF maybe not involved in nocamycin biosynthesis.
Inactivation of ncmG and isolation the new derivatives from the mutant strain
Cytochrome P450 oxidases are often play important roles in post-tailoring steps during antibiotic biosynthesis. Generally, oxygenation modification is a vital approach to improve bioactivity of parent molecule. To obtain more nocamycin derivatives, we inactivated ncmG by λ-RED/ET technology and generated ΔncmG mutant strain S. syringae MoS1003 (Additional file 1: Figure S3). HPLC analysis revealed that S. syringae MoS1003 abolished nocamycin I and nocamycin II production completely and two new peaks with similar characteristic UV absorption to these of nocamycin I and nocamycin II are detected (Fig. 2II). Then, A 15-L scale fermentation of ΔncmG mutant strain led the purification of nocamycin III and nocamycin IV. The structures of nocamycin III and nocamycin IV were determined by multiple spectroscopy data analyses. Both nocamycin III and IV are new nocamycin derivatives. Compared to nocamyin I and II, nocamycin III and IV show less oxidative modification, lacking of tetrahydrofuran ring, C-10 and C-21 modification.
The molecular formula of nocamycin III (4) is C25H35NO6 (m/z = 445.25), which was determined by HR–ESI–MS ([M − H]− m/z = 444.2422, [M + H]+ m/z = 446.2548, [M + Na]+ m/z = 468.2354) (Additional file 1: Figure S8). Comparisons of the 1H and 13C NMR spectroscopic data of nocamycin III to those of nocamycin I (Bu-2313B) suggested that they share a similar structure. Complete spectral data including COSY, HSQC, and HMBC spectra were also acquired (Additional file 1: Figures S10–S16), thereby allowing full assignments of the 1H and 13C signals (Table 4). Comparisons of the 1H and 13C NMR data for nocamycin I and nocamycin III revealed that the tetrahydrofuran ring is not closed and a Δ11,12 double bond is apparent in nocamycin III. HMBC correlations from H-20 to C-11, C-12, and C-13, and the COSY correlations of H-10/H-11 further substantiated these assignments. Additionally, H-15 was shifted to δH 4.29 due to the ring opening, relative to the same position of the cyclic form. A keto group in nocamycin I was replaced by a methylene group (δH, 1.98, 2.4; δC 23.9) at C-10 in 4, which was confirmed by the HMBC correlations from H-8, H-9, and H-11 to C-10 and from the COSY correlations of H-9/H-10α, and of H-10β/H-11. Another obvious difference observed from the 1H and 13C spectroscopic data was the absence of a –COOCH3 in 4 compared to that of nocamycin I. In turn, a methyl group (δH, 0.79, δC 11.8) was found to be attached at C-14. Cross peak of H-14/H-21 in the COSY spectrum and the HMBC correlations from H-21 to C-13, C-14 and C-15 further confirmed this assignment. Inspection of other NMR data for nocamycin III revealed other structural elements are identical to those of nocamycin I. Consequently, the structure of nocamycin III was elucidated as shown in Fig. 5.
Nocamycin IV (5) was isolated as a yellowish amorphous solid. Its molecular formula was determined as C25H35NO7 (m/z = 461.24) by HR–ESI–MS ([M − H]− m/z = 460.2355, [M + H]+ m/z = 462.2503, [M + Na]+ m/z = 484.2303) (Additional file 1: Figure S9), 16 mass units greater than that of nocamycin III, indicating one more oxygen atom than that of nocamycin III. Complete spectral data including COSY, HSQC, and HMBC spectra were also acquired (Additional file 1: Figures S17–S21), thereby allowing full assignments of the 1H and 13C signals (Table 4). It shared a similar structure to that of nocamycin III, except that a methyl signal at δH 1.62 was disappeared and an oxygen-bearing methylene signal at δH 3.96 and 4.07 occurred. Key HMBC correlations from H-20 to C-12 and from H-11 to C-12 further confirmed the location of the –CH2OH group at C-12. Thus, the structure of nocamycin IV was elucidated as 20-hydroxy-nocamycin III (Fig. 5).
In this study, the gene cluster responsible for nocamycin biosynthesis identified from S. syringae consists of 21 ORFs: 12 coding for structural proteins, seven involved in regulator and resistance and two with unknown function. Like the reported biosynthetic gene clusters of tirandamycin and streptolydigin, a hybrid PKS-NRPS mechanism is employed to assemble the chain PK-NPR backbone by co-linearity rule [10, 13, 14]. The core structure of nocamycin is bicyclic ketal unit and tetramic acid moiety. To date, tetramic acid structure has been identified in numerous natural products and four phylogenetically different family enzymes have been characterized to catalyze the tetramic acid formation through Dieckmann cyclisation reaction [21, 47,48,49,50]. In previous report, TrdC and its homologous protein SlgL have been characterized as Dieckmann cyclases to catalyze the formation of tetramic acid moiety in tirandamycin and streptolydigin biosynthetic pathway, respectively . Thus, it is plausible to assume that NcmC, the homologous protein to TrdC, is employed to generate tetramic acid moiety through Dieckmann cyclisation during nocamycins biosynthetic pathway .
Formation of bicyclic ketal ring represents the most intriguing issue of nocamycin family natural products, which is not fully understood. In our previous study, an abnormal DH at module 3 in tirandamycin PKS was proposed to be involved in spiroketal structure formation . Comparing with tirandamycin and streptolydigin gene clusters, it is important to notice that all the three gene clusters possess a similar unexpected DH domain with conserved active motif in the corresponding PKS. This abnormal DH domain at module 4 are likely not to catalyze the dehydration reaction to afford C-10 and C-11 double bond because the C-11 hydroxy group is absolutely required for the C-13 spiroketal group formation and no nocamycin derivatives possess C-10 and C-11 double bond have been identified. Recently, linear 7,13,9,13-diseco-tirandamycin derivative tirandamycin K, a shunt pathway product in tirandamycin pathway, was isolated from Streptomyces sp. 307-9 and its P450 monooxygenase disruption mutant strain . C-9 hydroxyl in tirandamycin K clearly indicates that DH3-catalyzed dehydration can be avoided, and it also provides evidence to support the mechanism that DH3 is involved in bicyclic ketal formation . Due to the high similarity in polyketide structure and domains organization of PKS between tirandamycin, nocamycin and streptolydigin gene clusters, the abnormal DH catalytic mechanisms are likely to be common spiroketalization mechanisms in these three pathways.
Based on bioinformatics and genetic engineering data, post tailoring steps of nocamycin can be predicted as follows (Fig. 4). Firstly, the earliest intermediate released from the PCP protein possesses a hydroxyl group in C-11 position, NcmE catalyze the dehydration process to afford nocamycin III. Next, nocamycin III undergoes several oxidative and one methyl esterification steps to produce nocamycin I. At last, NcmD catalyzes the dehydrogenation process to afford nocamycin II. Comparisons of gene clusters of tirandamycin and nocamycin revealed an interesting phenomenon that the post-tailoring enzymes involved in modification of similar structure are varied. In tirandamycin biosynthetic pathway, a FAD-dependent dehydrogenase TrdL/TamL is responsible for the conversion from C-10 hydroxyl to C-10 ketone [15, 16]. In our initial hypothesis, a TrdL/TamL homologous protein is predicted to be responsible for the same process, however, no TrdL/TamL homologous protein has been observed within the gene cluster. Although NcmL shows FAD-binding domain, it lacks the conserved bicovalent FAD linked active sites to that in TamL/TrdL [15, 16]. Meanwhile, the gene inactivation experiments revealed that NcmL shows no relationship to nocamycin biosynthesis, and this result also indicates that diversified modification mechanism occurred in this class of natural products. Overview the gene cluster, the short-chain dehydrogenase NcmD is the best candidate to catalyze the last C-10 dehydrogenation step in nocamycin biosynthetic pathway. The complex oxidative modifications including formation of fused oxolane ring system in bicyclic ketal moiety and the conversion from methyl group to carboxyl are intriguing issues, and the two cytochrome P450 oxidase NcmG and NcmO are expected to be involved in these steps. Two new derivatives nocamycin III and nocamycin IV lacking of closed tetrahydrofuran ring from ΔncmG mutant strain indicates NcmG’s involvement in the formation of the fused oxolane ring system. In terms of oxolane ring system formation, four different biosynthetic routes have been envisioned [52,53,54,55]. The mechanism of tetrahydrofuran ring in nocamycin is proposed to be similar to that in nonactin biosynthesis pathway . NcmG is likely to catalyze conjugate addition of C-15 hydroxyl groups to the adjacent C-11 and C-12 alkenyl moiety to form oxolane ring (Fig. 4). We notice that the C-20 hydroxyl in nocamycin IV is similar to C-18 hydroxyl in tirandamycin B. In tirandamycin biosynthetic pathway, a multifunctional cytochrome P450 TamI has been verified to be responsible for the formation of C-18 hydroxyl group . However, C-20 hydroxylation modification is not required in nocamycin biosynthetic pathway (Fig. 4). Thus, we hypothesize that nocamycin IV is probably a shunt product in nocamycins biosynthetic pathway and an oxidase located elsewhere of the genome can catalyze the hydroxylation process. Considerations of several oxidative modifications are required to afford nocamycin II, one of NcmG and NcmO is potentially responsible for more than one oxidative tailoring steps. Elucidation of the exact roles of NcmG and NcmO and the timing of modification in nocamycin biosynthesis is our ongoing project.
Up to now, the biosynthetic gene clusters responsible for streptolydigin, tirandamycin and nocamycin biosynthesis have been identified. Comparisons of the three gene clusters will help us deeply understand the biosynthetic mechanisms of this small class of natural products. The genetic insights and elucidations of enzyme function will facilitate us to rationally generate new derivatives with improved pharmacological property by manipulating biosynthetic pathway.
The nocamycin I and II, bearing a tricyclic ketal moiety, belong to acyl tetramic acid natural products and they display broad antimicrobial activity. In this report, we identify nocamycins biosynthetic gene cluster from rare actinomycete Saccharothrix syringae, which provides us the genetic insights into nocamycins biosynthesis and enzyme candidates for several intriguing biochemical transformations. Inactivation of cytochrome P450 monoxygenase NcmG led to isolation of two novel nocamycin derivatives from the mutant strain. Based on gene cluster data and new derivatives isolated from gene inactivation mutant strain, a putative biosynthetic pathway of nocamycin is proposed. These findings provide insights into further investigation of nocamycin biosynthetic mechanism, and also set the stage to rationally engineer new nocamycin derivatives via manipulating biosynthetic pathway.
open reading frame
nonribosomal peptide synthetases
- A domain:
peptidyl carrier protein
- C domain:
acyl carrier protein
- AprR :
apramycin resistant phenotype
- KanS :
kanamycin sensitive phenotype
high-performance liquid chromatography
nuclear magnetic resonance
Gauze G, Sveshnikova M, Ukholina R, Komarova G, Bazhanov V. Formation of a new antibiotic, nocamycin, by a culture of Nocardiopsis syringae sp. nov. Antibiotiki. 1977;22:483–6.
Brazhnikova M, Konstantinova N, Potapova N, Tolstykh I. Physicochmemical characteristics of the new antineoplastic antibiotic, nocamycin. Antibiotiki. 1977;22:486–9.
Horváth G, Brazhnikova MG, Konstantinova NV, Tolstykh IV, Potapova NP. The structure of nocamycin, a new antitumor antibiotic. J Antibiot (Tokyo). 1979;32:555–8.
Tsunakawa M, Toda S, Okita T, Hanada M, Nakagawa S, Tsukiura H, Naito T, Kawaguchi H. Bu-2313, a new antibiotic complex active against anaerobes II. Structure determination of Bu-2313 A and B. J Antibiot. 1980;33:166–72.
Tsukiura H, Tomita K, Hanada M, Kobaru S, Tsunakawa M, Fujisawa M, Fujisawa K, Kawaguchi H. Bu-2313, a new antibiotic complex active against anaerobes I. production, isolation and properties of Bu-2313 A and B. J Antibiot. 1980;33(2):157–64.
Bansal M, Dhawan V, Thadepalli H. In vitro activity of Bu-2313B against anaerobic bacteria. Chemotharapy. 1982;28:200–3.
Toda S, Nakagawa S, Naito T. Bu-2313, a new antibiotic complex active against anaerobes III. Semi-synthesis of Bu-2313 A and B. J Antibiot. 1980;33:166–71.
Tuske S, Sarafianos SG, Wang X, Hudson B, Sineva E, Mukhopadhyay J, Birktoft JJ, Leroy O, Ismail S, Clark AD Jr, Dharia C, Napoli A, Laptenko O, Lee J, Borukhov S, Ebright RH, Arnold E. Inhibition of bacterial RNA polymerase by streptolydigin: stabilization of a straight-bridge-helix active-center conformation. Cell. 2005;122:541–52.
Temiakov D, Zenkin N, Vassylyeva MN, Perederina A, Tahirov TH, Kashkina E, Savkina M, Zorov S, Nikiforov V, Igarashi N, Matsugaki N, Wakatsuki S, Severinov K, Vassylyev DG. Structural basis of transcription inhibition by antibiotic streptolydigin. Mol Cell. 2005;19:655–66.
Olano C, Gómez C, Pérez M, Palomino M, Pineda-Lucena A, Carbajo RJ, Braña AF, Méndez C, Salas JA. Deciphering biosynthesis of the RNA polymerase inhibitor streptolydigin and generation of glycosylated derivatives. Chem Biol. 2009;16:1031–44.
Carlson JC, Li S, Burr DA, Sherman DH. Isolation and characterization of tirandamycins from a marine-derived Streptomyces sp. J Nat Prod. 2009;72:2076–9.
Horna DH, Gómez C, Olano C, Palomino-Schätzlein M, Pineda-Lucena A, Carbajo RJ, Braña AF, Méndez C, Salas JA. Biosynthesis of the RNA polymerase inhibitor streptolydigin in Streptomyces lydicus: tailoring modification of 3-methyl-aspartate. J Bacteriol. 2011;193:2647–51.
Carlson JC, Fortman JL, Anzai Y, Li S, Burr DA, Sherman DH. Identification of the tirandamycin biosynthetic gene cluster from Streptomyces sp. 307-9. ChemBioChem. 2010;11:564–72.
Mo X, Wang Z, Wang B, Ma J, Huang H, Tian X, Zhang S, Zhang C, Ju J. Cloning and characterization of the biosynthetic gene cluster of the bacterial RNA polymerase inhibitor tirandamycin from marine-derived Streptomyces sp. SCSIO1666. Biochem Biophys Res Commun. 2011;406:341–7.
Carlson JC, Li S, Gunatilleke SS, Anzai Y, Burr DA, Podust LM, Sherman DH. Tirandamycin biosynthesis is mediated by co-dependent oxidative enzymes. Nat Chem. 2011;3:628–33.
Mo X, Huang H, Ma J, Wang Z, Wang B, Zhang S, Zhang C, Ju J. Characterization of TrdL as a 10-hydroxy dehydrogenase and generation of new analogues from a tirandamycin biosynthetic pathway. Org Lett. 2011;13:2212–5.
Mo X, Ma J, Huang H, Wang B, Song Y, Zhang S, Zhang C, Ju J. Δ(11,12) double bond formation in tirandamycin biosynthesis is atypically catalyzed by TrdE, a glycoside hydrolase family enzyme. J Am Chem Soc. 2002;134:2844–7.
Gómez C, Horna DH, Olano C, Palomino-Schätzlein M, Pineda-Lucena A, Carbajo RJ, Braña AF, Méndez C, Salas JA. Amino acid precursor supply in the biosynthesis of the RNA polymerase inhibitor streptolydigin by Streptomyces lydicus. J Bacteriol. 2011;193:4214–23.
Gómez C, Olano C, Méndez C, Salas JA. Three pathway-specific regulators control streptolydigin biosynthesis in Streptomyces lydicus. Microbiology. 2012;158:2504–14.
Gómez C, Horna DH, Olano C, Méndez C, Salas JA. Participation of putative glycoside hydrolases SlgC1 and SlgC2 in the biosynthesis of streptolydigin in Streptomyces lydicus. Microb Biotechnol. 2012;5:663–7.
Gui C, Li Q, Mo X, Qin X, Ma J, Ju J. Discovery of a new family of Dieckmann cyclases essential to tetramic acid and pyridone-based natural products biosynthesis. Org Lett. 2015;17:628–31.
Macneil DJ, Gewain KM, Ruby CL, Dezeny G, Gibbons PH, Macneil T. Analysis of Streptomyces-avermitilis genes required for avermectin biosynthesis utilizing a novel integration vector. Gene. 1992;111:61–8.
Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000;97:6640–5.
Gust B, Chandra G, Jakimowicz D, Yuqing T, Bruton CJ, Chater KF. Lambda red-mediated genetic manipulation of antibiotic-producing streptomyces. Adv Appl Microbiol. 2004;54:107–28.
Paget MSB, Chamberlin L, Atrih A, Foster SJ, Buttner MJ. Evidence that the extracytoplasmic function sigma factor sigma(e) is required for normal cell wall structure in Streptomyces coelicolor A3(2). J Bacteriol. 1999;181:204–11.
Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA. Practical streptomyces genetics. Norwich: John Innes foundation; 2000.
Weber T, Blin K, Duddela S, Krug D, Kim HU, Bruccoleri R, Lee SY, Fischbach MA, Müller R, Wohlleben W, Breitling R, Takano E, Medema MH. antiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 2015;43:W237–43.
Caffrey P. Conserved amino acid residues correlating with ketoreductase stereospecificity in modular polyketide synthases. ChemBioChem. 2003;4:654–7.
Haydock SF, Aparicio JF, Molnar I, Schwecke T, Khaw LE, Konig A, Marsden AFA, Galloway IS, Staunton J, Leadlay PF. Divergent sequence motifs correlated with the substrate-specificity of (methyl)malonyl-coa-acyl carrier protein transacylase domains in modular polyketide syntheses. FEBS Lett. 1995;374(2):246–8.
Reeves CD, Murli S, Ashley GW, Piagentini M, Hutchinson CR, McDaniel R. Alteration of the substrate specificity of a modular polyketide synthase acyltransferase domain through site-specific mutations. Biochemistry. 2001;40:15464–70.
Bevitt DJ, Cortes J, Haydock SF, Leadlay PF. 6-Deoxyerythronolide-B synthase from Saccharopolyspora erythraea: cloning of the structural gene, sequence-analysis and inferred domain-structure of the multifunctional enzyme. Eur J Biochem. 1992;204:39–49.
Nagano S, Cupp-Vickery JR, Poulos TL. Crystal structures of the ferrous dioxygen complex of wild-type cytochrome P450eryF and its mutants, A245S and A245T: investigation of the proton transfer system in P450eryF. J Biol Chem. 2005;280:22102–7.
Parajuli N, Basnet DB, Chan Lee H, Sohng JK, Liou K. Genome analyses of Streptomyces peucetius ATCC 27952 for the identification and comparison of cytochrome P450 complement with other Streptomyces. Arch Biochem Biophys. 2004;425:233–41.
Tang GL, Cheng YQ, Shen B. Leinamycin biosynthesis revealing unprecedented architectural complexity for a hybrid polyketide synthase and nonribosomal peptide synthetase. Chem Biol. 2004;11:33–45.
He HY, Pan HX, Wu LF, Zhang BB, Chai HB, Liu W, Tang GL. Quartromicin biosynthesis: two alternative polyketide chains produced by one polyketide synthase assembly line. Chem Biol. 2012;19:1313–23.
Chiu HT, Chen YL, Chen CY, Jin C, Lee MN, Lin YC. Molecular cloning, sequence analysis and functional characterization of the gene cluster for biosynthesis of K-252a and its analogs. Mol BioSyst. 2009;5:1180–91.
Chen C, Liu X, Abdel-Mageed WM, Guo H, Hou W, Jaspars M, Li L, Xie F, Ren B, Wang Q, Dai H, Song F, Zhang L. Nivetetracyclates A and B: novel compounds isolated from Streptomyces niveus. Org Lett. 2013;15:5762–5.
Kavanagh KL, Jörnvall H, Persson B, Oppermann U. Medium- and short-chain dehydrogenase/reductase gene and protein families: the SDR superfamily: functional and structural diversity within a family of metabolic and regulatory enzymes. Cell Mol Life Sci. 2008;65:3895–906.
Mattheus W, Masschelein J, Gao LJ, Herdewijn P, Landuyt B, Volckaert G, Lavigne R. The kalimantacin/batumin biosynthesis operon encodes a self-resistance isoform of the FabI bacterial target. Chem Biol. 2010;17:1067–71.
Walker JE, Saraste M, Runswick MJ, Gay NJ. Distantlyrelated sequences in the alpha- and beta-subunits of ATP synthase, myosin, kinases and other ATP-requiring enzymes and a common nucleotide binding fold. EMBO J. 1982;1:945–51.
Iyer LM, Leipe DD, Koonin EV, Aravind L. Evolutionary history and higher order classification of AAA + ATPases. J Struct Biol. 2004;146:11–31.
Barthelmebs L, Lecomte B, Diviès C, Cavin JF. Inducible metabolism of phenolic acids in Pediococcus pentosaceus is encoded by an autoregulated operon which involves a new class of negative transcriptional regulator. J Bacteriol. 2000;182:6724–31.
Gury J, Barthelmebs L, Tran NP, Diviès C, Cavin JF. Cloning, deletion, and characterization of PadR, the transcriptional repressor of the phenolic acid decarboxylase-encoding padA gene of Lactobacillus plantarum. Appl Environ Microbiol. 2004;70:2146–53.
De Silva RS, Kovacikova G, Lin W, Taylor RK, Skorupski K, Kull FJ. Crystal structure of the virulence gene activator AphA from Vibrio cholerae reveals it is a novel member of the winged helix transcription factor superfamily. J Biol Chem. 2005;280:13779–83.
Madoori PK, Agustiandari H, Driessen AJ, Thunnissen AM. Structure of the transcriptional regulator LmrR and its mechanism of multidrug recognition. EMBO J. 2009;28:156–66.
Xu Z, Jakobi K, Welzel K, Hertweck C. Biosynthesis of the antitumor agent chartreusin involves the oxidative rearrangement of an anthracyclic polyketide. Chem Biol. 2005;12:579–88.
Blodgett JA, Oh DC, Cao S, Currie CR, Kolter R, Clardy J. Common biosynthetic origins for polycyclic tetramate macrolactams from phylogenetically diverse bacteria. Proc Natl Acad Sci USA. 2010;107:11692–7.
Wu Q, Wu Z, Qu X, Liu W. Insights into pyrroindomycin biosynthesis reveal a uniform paradigm for tetramate/tetronate formation. J Am Chem Soc. 2012;134:17342–5.
Sims JW, Schmidt EW. Thioesterase-like role for fungal PKS-NRPS hybrid reductive domains. J Am Chem Soc. 2008;130:11149–55.
Liu X, Walsh CT. Cyclopiazonic acid biosynthesis in Aspergillus sp.: characterization of a reductase-like R* domain in cyclopiazonate synthetase that forms and releases cyclo-acetoacetyl-l-tryptophan. Biochemistry. 2009;48:8746–57.
Zhang X, Zhong L, Du L, Chlipala GE, Lopez PC, Zhang W, Sherman DH, Li S. Identification of an unexpected shunt pathway product provides new insightsinto tirandamycin biosynthesis. Tetrahedron Lett. 2016;57:5919–23.
Woo AJ, Strohl WR, Priestley ND. Nonactin biosynthesis: the product of nonS catalyzes the formation of the furan ring of nonactic acid. Antimicrob Agents Chemother. 1999;43:1662–8.
Demydchuk Y, Sun Y, Hong H, Staunton J, Spencer JB, Leadlay PF. Analysis of the tetronomycin gene cluster: insights into the biosynthesis of a polyether tetronate antibiotic. ChemBioChem. 2008;9:1136–45.
Bode HB, Zeeck A. Biosynthesis of kendomycin: origin of the oxygen atoms and further investigations. J Chem Soc Perkin Trans. 2000;1:2665–70.
Richter ME, Traitcheva N, Knüpfer U, Hertweck C. Sequential asymmetric polyketide heterocyclization catalyzed by a single cytochrome P450 monooxygenase (AurH). Angew Chem Int Ed. 2008;47:8872–5.
XM and QW conceived and designed the project. XM, CS, CG, YZ and JJ conducted and analyzed the data, XM, JJ and QW wrote the paper. All authors read and approved the final manuscript.
We are so grateful to Agricultural Research Service Culture Collection (NRRL) for providing S. syringae NRRL B-16468. Research in our laboratory was supported by Grants from the National Natural Science Foundation of China (Nos. 31300063, 31300005), Natural Science Foundation of Shandong Province (NO. ZR2013CL020), Research Foundation for Advanced Talents of Qingdao Agricultural University (No. 631301) and Open Research Foundation of CAS Key Laboratory of Tropical Marine Bio-resources and Ecology (No. LMB141007).
The authors declare that they have no competing interests.
Availability of data and materials
The dataset(s) supporting the conclusions of this article is (are) included within the article [and its additional file(s)].
This article does not contain any studies with human participants or animals performed by any of the authors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.