Identification and engineering of 32 membered antifungal macrolactone notonesomycins
Microbial Cell Factories volume 19, Article number: 71 (2020)
Notonesomycin A is a 32-membered bioactive glycosylated macrolactone known to be produced by Streptomyces aminophilus subsp. notonesogenes 647-AV1 and S. aminophilus DSM 40186. In a high throughput antifungal screening campaign, we identified an alternative notonesomycin A producing strain, Streptomyces sp. A793, and its biosynthetic gene cluster. From this strain, we further characterized a new more potent antifungal non-sulfated analogue, named notonesomycin B. Through CRISPR–Cas9 engineering of the biosynthetic gene cluster, we were able to increase the production yield of notonesomycin B by up to 18-fold as well as generate a strain that exclusively produces this analogue.
Natural products (NPs) represent a diverse source of bioactive compounds with clinical relevance, as evidenced by the fact that half of current drugs are NP or NP-derived [1, 2]. They can be further engineered to best suit various clinical, flavor and agrochemical applications [3, 4]. The A*STAR Natural Organism Library (NOL)  is a Singapore-based collection of over 120,000 microbial strains, which include 58,000 actinomycetes collected from diverse habitats. Previous screening of this collection has uncovered a number of bioactive compounds, including antibacterial compounds microsphaerins A–D , gargantulide A , and anthracimycin , and antagonists of the human CCR5 receptor; fuscinarin, fuscin, and cochlioquinones .
Candida species are common pathogens associated with invasive fungal diseases contributing to high disease burden [10,11,12]. Increasing number of immunocompromised patients , emerging new fungal infections  and drug-resistant fungal pathogens [15, 16] have rekindled interest in new antifungal therapeutics. Bioactivity-guided profiling of NOL extracts against Candida albicans uncovered the 32-membered macrolactone notonesomycin A (Fig. 1). Known to have antifungal property, notonesomycin A was first reported in Streptomyces aminophilus subsp. notonesogenes 647-AV1 in 1986 , from an assay against rice sheath blight disease. Characterization of similar compounds, such as brasilinolides [18, 19], ibomycin , langkolide , novonestmycins A and B , PM100117 and PM100118 [23, 24], demonstrates the bio-versatility of these 32–36 membered macrolactones. For this family of compounds, immunosuppressive, anti-tumor, and antifungal bioactivities have been reported. A major difference between these macrolactones and notonesomycin A lies in the addition of a sulfate group at position C-17 by O-sulfation, which is a modification seen in a limited number of microbial natural products. Interestingly, we also identified for the first time, a non-sulfated bioactive analogue of notonesomycin A, henceforth named notonesomycin B (Fig. 1).
Advances in DNA sequencing and sequence analysis have driven deeper genomic insights into streptomycetes [25, 26]. In addition, genetic engineering of previously challenging strains is now possible with CRISPR–Cas9 [27,28,29,30,31]. Combining the understanding of biosynthetic gene clusters (BGCs) and improved ability to edit streptomycete genomes, combinatorial engineering of analogue libraries for the purpose of producing a final therapeutically relevant product can now be accomplished at a faster pace.
In this study, we present the characterization of notonesomycin BGC from Streptomyces sp. A793 as well as the isolation of a new non-sulfated notonesomycin analogue. We also demonstrated the application of CRISPR–Cas9 in Streptomyces sp. A793 to increase the production of the more bioactive non-sulfated notonesomycin B.
Isolation and identification of notonesomycins A (1) and a new analogue notonesomycin B (2) from Streptomyces sp. A793
From bioactivity based profiling of NOL strain library, extracts generated from Streptomyces sp. A793 were active against Candida albicans. Further bioassay-guided isolation and identification of compounds from Streptomyces sp. A793 led to the isolation and characterization of 2 compounds (1, 2).
Compound 1 (Fig. 1, Additional file 1: Figure S1) was isolated as a brownish amorphous solid with [α]D − 25 (c 0.9, MeOH). ESI-MS showed an [M+H]+ peak at m/z 1454 and the HR-ESI-MS measurement established the molecular formula as C68H111NO30S. The 13C NMR data (Additional file 1: Table S1 and Figure S2) showed a total of 68 carbon resonances, comprising 4 carbonyl (δ 167.9, 168.4, 170.5, 172.3), eight olefinic (δ 109.1, 116.9, 117.2, 124.4, 127.5, 144.8, 146.9, 148.5), three carbons bonded to 2 heteroatoms (δ 98.2, 99.1, 102.6), 21 oxymethine (δ 60.3–85.6), five methine (δ 36.0, 38.3, 38.6, 40.2, 40.5), 17 methylene (δ 22.7–44.4), one methoxy (δ 61.1), one N-methyl (δ 29.8), and eight methyl (δ 4.8, 9.5, 10.4, 11.1, 14.4, 15.3, 17.5, 18.4) carbon atoms. The 1H NMR spectrum (Additional file 1: Table S1 and Figure S3) showed the presence of three aromatic signals [δ 7.63 (d, J = 2 Hz), 6.58 (d, J = 9 Hz), 7.68 (dd, J = 9, 2 Hz)], indicating a tri-substituted aromatic ring system. The remaining two olefinic proton signals resonated at δ 5.95 and 7.04 with a large coupling constant of 15 Hz suggesting a E configuration and were assigned to an α, β-unsaturated ester group based on the observed δH 5.95 (H-2) and 7.04 (H-3) to δC 168.4 (C-1) HMBC correlations (Additional file 1: Figure S1). Two downfield methyl singlets observed at δ 3.50 and 2.87 (3 protons) were assigned to the O-methyl and N-methyl, respectively. Eight methyl signals were observed at δ 0.76 (d, J = 7 Hz), 0.84 (d, J = 7 Hz), 0.90 (d, J = 7 Hz), 0.96 (d, J = 7 Hz), 0.97 (t, J = 7 Hz), and 1.06 (d, J = 7 Hz), 1.27 (d, J = 7 Hz), and 1.32 (d, J = 6 Hz) and were assigned to H-44, H-42, H-43, H-40, H-39, H-41, H-6ʹʹʹ and H-6ʹ, respectively. The two oxymethines observed at δH 2.79 (δC 63.8) and 2.73 (δC − 60.3) were characteristic of an epoxide moiety with cis-orientation based on the observed small coupling constant of 2 Hz [32, 33]. The anomeric protons of the two deoxysugar moieties were deduced from the observed doublet and doublet of doublets signals at δH 4.64 and 5.11 based on COSY, HSQC and HMBC NMR data (Additional file 1: Figures S4–S6). The two deoxysugar units were linked together by the trisubstituted benzene ring B via ester group to deoxysugar ring A from H-3ʹ to C-1ʹʹ and ether linkage of ring B to deoxysugar ring C from H-1ʹʹʹ to C-4ʹʹ 3J HMBC correlations (Additional file 1: Figures S1, S6). The attachment of deoxysugar ring A to C-37 was evident from anomeric proton (H-1ʹ) to C-37 HMBC correlation. Presence of malonate ester was confirmed by two carbonyl signals at δC 170.5 and 172.3 but the isolated methylene signals (δH 3.17 and δC 43.0) were observed only in DMSO-d6 and not in methanol-d4 due to deuterium exchange [17, 34]. The placements of malonate ester group and the remaining sulfate group at C-21 [34, 35] and C-17 , were assigned based on the observed carbon values of δC 71.0 and 78.2, respectively. From 1D and 2D NMR data, compound 1 was a 32-membered macrocyclic lactone (HMBC: from H-31 to lactone carbonyl) consisting of a sulfate group, a malonate ester, an epoxide, a trisubstituted benzene ring and 2 deoxysugar moieties similar to that of notonesomycin A. In addition, the 1H and 13C NMR data of 1 were consistent with those reported for notonesomycin A [17, 34]. Hence, 1 was identified as notonesomycin A.
Compound 2 (Fig. 1) was isolated as a white amorphous solid with [α]D − 56 (c 1.3, MeOH). ESI-MS showed an [M+H]+ peak at m/z 1374 and the HR-ESI-MS measurement established the molecular formula as C68H111NO27. The difference in 80 mass units as compared to 1 suggested a loss of a [SO3] group. The 1H, 13C, COSY, HSQC and HMBC NMR spectra (Additional file 1: Figures S7–S11) indicated the presence of two deoxysugars [δ 4.64 (brd, J = 8 Hz), δ 5.11 (dd, J = 9, 2 Hz)], an epoxide (δ 2.72, 2.78), a trisubstituted benzene ring [δ 7.62 (d, J = 2 Hz), δ 6.58 (dd, J = 9, 2 Hz), δ 5.11 (dd, J = 9, 2 Hz)], a macrolide incorporating a hemiketal ring, and an α, β-unsaturated lactone moiety [δ 5.95 (d, J = 15 Hz), δ 7.03 (ddd, J = 15, 8, 8 Hz)] similar to those of 1. The 1D and 2D NMR (Additional file 1: Figures S7–S11) of 2 indicated that 1 and 2 have similar core structures, except for the signals of C-16 (δH 3.63, δC 75.9), C-17 (δH 3.87, δC 70.0) and C-18 (δH 1.89, 1.25, δC 40.7) of the 6-membered hemiketal moiety embedded within the macrocyclic ring. Since there was no change in the total number of carbon and protons in 2 as compared to 1, the observed 80 mass unit loss in HRMS data indicated the presence of a hydroxyl group instead of the sulfate group at C-17. This was further supported by the observed 8.2 ppm upfield shift for C-17 in 2 (2: δC 70.0; 1: δC78.2)  and HMBC correlation from H-16 to C-17 (Additional file 1: Figures S1, S11). As in 1, the presence of malonate ester group was evidenced by the methylene singlet at δH 2.91, δC 45.9 in DMSO-d6 NMR data. Hence, compound 2 was identified as a desulfated derivative of 1 and the gross structure of 2 is shown in Fig. 1. The relative configurations at most stereogenic centers were determined by comparison with the 13C NMR data of novonestmycins A and B  and brasilinolides A and C  except for those at C-5, -7, -12, -13, and -14, which remain unassigned. The relative configurations of 1 were assumed to be the same as 2 based on high similarity of their NMR data (Fig. 1). 2 is therefore a desulfated derivative of 1 and named notonesomycin B.
Biological activities of notonesomycins A and B
Antifungal activities of notonesomycins A and B were determined against C. albicans ATCC® 90028™ and Aspergillus fumigatus ATCC® 46645™ (Table 1). Compared to notonesomycin A, the non-sulfated notonesomycin B was four times more active against A. fumigatus and five times more active against C. albicans. Both compounds were further profiled against commonly used laboratory human cell lines and they displayed IC50 values of ≤ 1 μM (Table 1). Between the two compounds, notonesomycin B was three times less active against the tested human cancer cell lines. Both notonesomycins A and B also displayed activity against Staphylococcus aureus Rosenbach ATCC® 25923™ with a minimum inhibition concentration of around 20 μM (Table 1) but were inactive against tested Gram-negative bacteria.
In silico identification and sequence architecture of a biosynthetic gene cluster indicative for notonesomycin biosynthesis
Whole-genome sequencing by Pacific Biosciences RS II of Streptomyces sp. A793 yielded a ~ 8.0 Mb contig and an 80 kb plasmid. The 8.0 Mb contig had a GC content of 72.4% and within this, 7027 coding sequences and 73 RNAs were identified by RAST (Rapid Annotation using Subsystem Technology, ). Based on 16S rRNA sequence alignment, the strain is shown to be closely associated with Streptomyces armeniacus. Analysis by antiSMASH 3.0  identified a 153 kb BGC with 59 open reading frames indicative for notonesomycin biosynthesis (Fig. 2a). BLASTp against the NCBI non-redundant sequence database  was used to annotate the 59 orfs (Additional file 1: Table S2). The polyketide scaffold is encoded by 18 polyketides synthase (PKS) modules within 8 proteins (nbc20, 21 and nbc38-43), which are consistent with the polyketide backbone of notonesomycin. The putative biosynthetic pathway (order of reactions, extender unit incorporation and tailoring reactions in the polyketide synthases) align mostly with the macrolactone (Fig. 2), with some exceptions. AT (module 17, Nbc43) is similar to AT (module 19, GonP7) from PM100117 and PM100118 BGC (62.4% similarity). Both acyltransferases are predicted to be specific for methylmalonyl CoA, according to homology predictions (Additional file 1: Figure S12) . However, in their putative biosynthetic pathways, a malonyl CoA extender unit is predicted to be incorporated. A redundant dehydratase was also observed in module 12, Nbc40. In notonesomycins A and B, a moiety consisting of 4-amino, 3-hydroxy benzoic acid and two deoxysugars is appended onto C-37 (Fig. 2b). The 4-amino, 3-hydroxy benzoic acid is predicted to be produced by aminobenzoate synthase, a chorismate lyase and a 4-hydroxybenozate 3-monooxygenase, which are encoded by nbc57-59 respectively. Thymidine diphosphate (TDP)-4-keto-6-deoxy-d-glucose is synthesized by glucose-1-phosphate thymidylyltransferase and dTDP-glucose-4,6-dehydratase that are encoded by nbc24 and nbc25 respectively . A putative set of genes (nbc15-17) encoding for an epimerase, TDP-hexose-3-ketoreductase and TDP-hexose 2, 3-dehydratase are proposed for the biosynthesis of the final two deoxysugars from TDP-4-keto-6-deoxy-d-glucose (Additional file 1: Figure S13). No additional ketoreductase and dehydratase are found for sugar biosynthesis within the cluster, leading us to suspect that nbc15-17 may be promiscuous in their substrate and stereochemistry specificities . O-methylation on the deoxysugar is predicted to occur via a methyltransferase encoded by nbc56. Within the cluster, two putative glycosyltransferases were observed (Nbc18 and Nbc22) (Additional file 1: Figure S14). Nbc18 is a predicted YjiC type glycosyltransferase. Based on the demonstrated flexibility of YjiC , we proposed that the two deoxysugar moieties are loaded onto the 4-amino-3-hydroxybenzoic acid by Nbc18. Subsequently, glycosylation of the C–OH of the polyketide is most likely catalyzed by the Nbc22 (Additional file 1: Figure S14).
Genes nbc33, nbc36, nbc37 encode for three P450 homologs, which are predicted to participate in the hydroxylation of C-16, C-14 and the epoxidation of C-28 (Additional file 1: Figure S15). Esterification is also predicted to proceed via PapA5-like alpha/beta-hydrolase acyltransferase, Nbc13, using malonyl-CoA as a substrate. Nbc1 is predicted to be an S-adenosylmethionine-dependent methyltransferase responsible for the C-methylation of C-38 in the polyketide chain (Additional file 1: Figure S16).
Experimental validation that the described BGC is responsible for the production of notonesomycins A and B
In order to validate that the BGC described above is related to the production of notonesomycins A and B, the CRISPR–Cas9 mediated genomic deletion of nbc20 and nbc21 was carried out. As expected, deletion of this cluster segment ablated the production of notonesomycins A and B (Fig. 2c, Additional file 1: Figure S17).
Increase in yield of production by SARP overexpression
The presence of a predicted Streptomyces antibiotic regulatory protein (SARP) family regulator (nbc14) encoded in the notonesomycin BGC presented an opportunity to improve metabolite production. Using CRISPR–Cas9-mediated genetic editing, we inserted a strong constitutive kasO* promoter  in front of nbc14. Since the insertion is located in the intergenic region between nbc13 and nbc14, a P8 promoter  was also inserted in front of nbc13 to minimize disruption to nbc13 expression. Consistent with nbc14 being a transcriptional activator, increased transcription of nbc14 and selected BGC gene was observed in the engineered strain (nbc14 overexpression, OE, mutant) (Fig. 3). In the wild type Streptomyces sp. A793, notonesomycin production was limited in CA10LB media with soluble starch as the sole carbon source, while in SV2 media containing glucose and glycerol, notonesomycin production level was at least 100-fold higher than in CA10LB. The engineered strain was grown in SV2 and CA10LB to assess changes in notonesomycin production. Corresponding to the increase in BGC gene expression, 18- and 11-fold increase in notonesomycin A production in SV2 and CA10LB media was observed in the engineered strain (Fig. 3). Similarly, notonesomycin B yields were increased 3- and 12-fold, respectively in SV2 and CA10LB media with nbc14 overexpression.
Sulfonate incorporation is predicted to occur via a common sulfate donor, 3′-phosphoadenosine-5′-phosphosulfate (PAPS). Biosynthesis of PAPS is predicted to be directed by an operon encoding for a heterodimeric sulfate adenylytransferase (nbc10 and nbc11), and adenylyl-sulfate kinase (nbc9), which are similar to CysC, D and N in E. coli  and AziH1-3 of the azinomycin BGC , respectively. Sulfonate transfer onto the C-17 accepting hydroxyl group is predicted to be catalyzed by sulfotransferase Stf3 enzyme (Nbc48). Phylogenetic analysis of the sulfotransferases showed that Nbc48 and other sulfotransferases from available Streptomyces genomes formed a clade with Stf3 from Mycobacterium tuberculosis (Additional file 1: Figure S18). Similar to Nbc48, Stf3 is a sulfotransferase found in the sulfomenaquinone BGC, which also includes the PAPS operon [46, 47].
Notably, ratios of sulfated notonesomycin A to the more bioactive notonesomycin B varied significantly in different media conditions. For wild type Streptomyces sp. A793 grown in SV2, notonesomycin B is the major product with a notonesomycin A:B ratio of 1:8 (Fig. 3a). Whilst in NotA media, notonesomycin A is the dominant product with a 3:1 notonesomycin A:B ratio (Fig. 4c). Sources of sulfate in microorganisms come from the trans-sulfuration of media amino acids such as methionine, serine, and cysteine , as suggested by the presence of a cystathionine gamma-synthase (nbc3) that catalyzes the biosynthesis of cysteine. In general, the preferred production of notonesomycin A in NotA media by Streptomyces sp. A793 correlated to the higher levels of sulfur-containing amino acids, which is at least fourfold higher in NotA media than SV2 media (Additional file 1: Table S3). To selectively produce notonesomycin B, we deleted nbc48 in Streptomyces sp. A793 and as expected, notonesomycin B was exclusively produced by nbc48 KO strains (Fig. 4b, c).
Within the 32–36 membered macrolactones, notonesomycin is notable for its unique sugar moieties and also sulfated functional group. In notonesomycin, its deoxysugars are combined with 4-amino, 3-hydroxy benzoic acid, an exception among the 32–36 membered ring compounds with repertoires of solely deoxysugars combinations. In this study, we not only characterized the biosynthetic pathway of notonesomycin A but we also identified a new non-sulfated analogue, notonesomycin B. We have shown that O-sulfation of the C-17 hydroxyl group reduces its fungicidal activity while increasing its cytotoxicity. In natural product biosynthesis, sulfation is an unusual, poorly understood phenomenon and there are only a handful of naturally occurring sulfated compounds, such as azinomycin , glycopeptide A47934 , linear polyenes clethramycin and mediomycins , liponucleoside antibiotics Caprazamycins , β-lactam sulfated carbapenem MM4550  and deplelides A and B . Enzymatic sulfation has been deployed to increase structural diversity in several classes of natural products . In particular, engineered sulfation of glycopeptide Teicoplanin A47934  yielded bioactive derivatives. With genes encoding both substrates and sulfation reactions elucidated within the notonesomycin BGC, it is possible to engineer sulfation. Here, we used this knowledge to selectively produce the more bioactive notonesomycin analogue. The deletion of Nbc48 sulfotransferase is sufficient to abate notonesomycin A production in favor of notonesomycin B.
Polyketide biosynthesis pathways are modulated by carbon and nitrogen sources through complex regulatory networks involving gene expression, enzyme activities, distribution and abundance of precursors [55,56,57,58]. For notonesomycin biosynthesis, the flow of precursors depends on available carbon sources and nutrient utilization, thereby limiting the production yield and relative abundance of the analogues. Yields, analogue distributions and gene expression were vastly different in glycerol and glucose based SV2 culture compared to a starch only CA10LB culture. In addition, through rational engineering, we demonstrated that Streptomyces sp. A793 has not yet reached its full production potential. Overexpression of a pathway specific SARP regulator further improved the yields of notonesomycin A and B in both media.
In a screen against Candida albicans, Streptomyces sp. A793, which produces notonesomycins A and B, was identified. Compared to notonesomycin A, the non-sulfated notonesomycin B demonstrated higher fungicidal and lower cytotoxicity activity. With the structure and BGC elucidation, along with CRISPR–Cas9-mediated editing, we were able to selectively produce notonesomycin B in Streptomyces sp. A793 and increase the production yields.
Bacterial strains, media and cultivation conditions
Notonesomycin-producing strain Streptomyces sp. A793 was obtained from the strain library of A*STAR . The strain was sub-cultured on Bennett’s agar (1% glucose, 0.1% yeast extract, 0.1% lab lemco, 5% glycerol, 0.2% casitone, 1.5% agar) for 7 days at 28 °C. Seed culture was grown in SV2 for 7 days at 28 °C with shaking at 200 rpm. Shake flasks fermentation was carried out in CA12LB, CA10LB, SV2 or NotA media  (Additional file 1: Table S3). Fermentation was carried out for 9 days at 28 °C with shaking at 200 rpm. Large scale fermentation of 6 L for (1) isolation was carried out in NotA media and 4 L fermentation for (2) isolation was carried out in SV2 media. Chemicals and reagents for media preparation were purchased from Sigma or Oxoid. Cell enumeration was determined using drop plate technique and cell viability was carried out using CellTiter-Blue® cell viability assays kits (Promega) according to manufacturer’s instructions and fluorescence recorded with a plate reader.
Metabolites extraction and isolation
For notonesomycin B (2): The culture broths of Streptomyces sp. A793 were combined and centrifuged to separate the supernatant and the mycelia. The combined supernatant was freeze-dried, partitioned with DCM:MeOH:H2O (1:1:1). The organic layer was then evaporated to dryness. The dried crude extract (1.5 g) was re-dissolved in DMSO and separated by C18 reversed-phase preparative HPLC (solvent A: H2O+ 0.1% HCOOH, solvent B: acetonitrile+ 0.1% HCOOH; flow rate: 30 mL/min, gradient conditions: 70:30 isocratic for 5 min; 30% to 70% of solvent B over 85 min, followed by 70% to 100% of solvent B over 10 min, and finally isocratic at 100% of solvent B for 15 min). Fractions eluting at 41–45 min (Fraction A) and 51 min (Fraction B) from HPLC were concentrated and dried under reduce pressure. The dried fraction A from HPLC was then subjected to C18 reversed-phase preparative HPLC for further purification and yielded 1.75 mg/L of notonesomycin B (2). Fraction B (notonesomycin A detected) was not sufficient for further purification. For notonesomycin A (1): The culture broths of Streptomyces sp. A793 were combined and centrifuged to separate the supernatant and the mycelia. The combined supernatant was charged to 10 cm × 8 cm Sepra C18-E (50 μm, 65 Å, Phenomenex) column. The column was eluted by isocratic gradient of 20%, 40%, 60%, 80%, and 100% aqueous methanol. The 80% aqueous methanol fraction was concentrated and dried under reduced pressure to yield 1.6 g of partially enriched fraction of notonesomycin A (1). The dried extract was then re-dissolved in 2 mL of MeOH and further separated by Sephadex LH-20 column (mobile phase: MeOH) to give 1.67 mg/L of notonesomycin A (1).
Chemical structural data
The NMR spectra of notonesomycins A and B are provided (Additional file 1: Figures S2–S11).
Notonesomycin A (1)
Brownish amorphous solid; [α]D − 25 (c 0.9, MeOH); UV (MeOH) λmax (log ε) 228 sh (4.27) and 307 (3.91) nm; HR-ESI-MS m/z 1454.6985 [M+H]+ (calcd for C68H111NO30S + H, 1454.6984); 1H and 13C NMR data, see Additional file 1: Table S1.
Notonesomycin B (2)
White amorphous solid; [α]D − 56 (c 1.3, MeOH); UV (MeOH) λmax (log ε) 229 sh (3.94) and 309 (3.89) nm; HR-ESI-MS m/z 1374.7414 [M+H]+ (calcd for C68H111NO27 + H, 1374.7416);1H and 13C NMR data, see Additional file 1: Table S1.
General experimental procedures
Optical rotations were recorded on a JASCO P-2000 digital polarimeter. UV spectra were obtained on a GE Healthcare Ultrospec 9000 spectrophotometer. NMR spectra were collected on a Bruker DRX-400 NMR spectrometer with Cryoprobe, using 5-mm BBI (1H, G-COSY, multiplicity-edited G-HSQC, and G-HMBC spectra) or BBO (13C spectra) probe heads equipped with z-gradients. Spectra were calibrated to residual protonated solvent signals (CD3OD δH 3.30 and CD3OD δC 49.0; DMSO-d6 δH 2.49 and DMSO-d6 δC 39.5). Preparative HPLC analysis was performed on the Agilent 1260 Infinity Preparative-Scale LC/MS Purification System, completed with Agilent 6130B single quadrupole mass spectrometer for LC and LC/MS Systems. The samples were separated on an Agilent Prep C18 column (100 × 30 mm) by gradient elution with a mixture of 0.1% formic acid in water (solvent A) and 0.1% formic acid in acetonitrile (solvent B). The HR-ESI-MS spectra were acquired on Agilent UHPLC 1290 Infinity coupled to Agilent 6540 accurate-mass quadrupole time-of-flight (QTOF) mass spectrometer equipped with a splitter and an ESI source. The analysis was performed with a C18 4.6 × 75 mm, 2.7 µm column at flowrate of 2 mL/min, under standard gradient condition of 0.1% formic acid in water and 0.1% formic acid in acetonitrile over 15 min.
Minimum inhibition concentration (MIC) was determined using the microbroth dilution method according to guidelines of the Clinical Laboratory Standards Institute (CLSI) with slight modifications. Assays were carried out with C. albicans ATCC® 90028™ at 2.5 × 105 cells mL−1, A. fumigatus ATCC® 46645™ at 2.5 × 104 spores mL−1 and both S. aureus Rosenbach ATCC® 25923™ and Escherichia coli (ATCC® 25922™) at 5 × 105 cells mL−1. The cells were incubated with the compounds at 35 °C for 24 h (C. albicans) and 48 h (A. fumigatus); and 37 °C for 24 h for both bacteria. The effect of the compounds on bacterial growth was evaluated by measuring the optical density at 600 nm using a microplate reader. Mammalian cell cytotoxicity was assessed with A549 human lung carcinoma cells (ATCC® CCL-185™) seeded at 1500 cells well−1, while both PANC-1 human pancreatic carcinoma cells (ATCC® CRL-1469™) and HepG2 human hepatocellular carcinoma cells (ATCC® HB-8065™) which were both seeded at 2500 cells well−1. These cells were treated with the compounds for 72 h at 37 °C in the presence of 5% CO2. PrestoBlue™ cell viability reagent (Life Technologies) was used to assess the cytotoxic effect of the compounds. Microplates were incubated with this dye for 2 h before measuring the fluorescence at excitation 560 nm and emission 590 nm. All assays were performed in triplicate on two different test runs.
Isolation of genomic DNA from Streptomyces sp. A793
High-quality and molecular weight genomic DNA (> 20 kb) was extracted from Streptomyces sp. A793 with modified CTAB method . Cells were suspended in TE buffer (10 mM Tris; 1 mM EDTA, pH 8.0) containing 20 µL of lysozyme (100 mg mL−1) and 15 µL of RNase (10 mg mL−1) at 37 °C overnight. Subsequent procedures of the CTAB method were carried out, followed by phenol–chloroform purification. An additional purification step was carried out using Mag-Bind® RxnPure Plus (OMEGA bio-tek) as per manufacturer’s protocol.
Genome sequencing, putative biosynthetic gene cluster annotation and analysis
Genome sequencing was performed with the Pacific Biosciences (PacBio) RSII. Purified genomic DNA was sheared to approximately 20 kb using a g-Tube (Covaris). A SMRTbell library was prepared according to manufacturer’s instructions, loaded with a MagBead bound library protocol onto two SMRTCells at 0.125 and 0.3 nM, and sequenced using the P5-C3 chemistry on the PacBio RSII instrument (Pacific Biosciences) with a 180 min movie time. De novo assembly was performed with the Hierarchical Genome Assembly Process 3 (HGAP3)  in the SMRT Analysis suite (version 2.3) using all default parameters. Putative BGCs were identified using antiSMASH 3.0 server .
A nucleotide BLAST search of the 16S rRNA obtained from the A793 genome sequence against the NCBI 16S ribosomal RNA database  revealed that their copies shared 99.7% to 99.9% sequence identity (E-value = 0.0 for all hits) to 16S rRNAs of S. armeniacus strain JCM 3070 (accession numbers NR_112046 and NR_112047) covering a length of 1530 nucleotides. In addition, we found a 100% identity match of a 1330 nucleotide sequence stretch of the A793 strains’ 16S rRNAs to the 16S rRNA partial sequence available in the draft genome of S. armeniacus strain ATCC 15676 (accession CP031320; ). This evidence is further supported by the finding that the nucleotide sequences of four housekeeping genes, rpoB, dnaK, atpD, and recA, found in both genomes are 100% identical. Consequently, strain A793 was identified to be closely associated to S. armeniacus.
Protein annotation and analysis
The identification of open reading frames (ORFs) was obtained with an in-house pipeline that takes into account the gene predictions by Prodigal 2.6.2 , Glimmer 3.02 [63, 64] and GeneMark 3.26 . BLASTp searches against the NCBI non-redundant sequence database  were carried out for the Nbc1-59 functional annotation based on sequence homology (Additional file 1: Table S2). Proteins without a significant hit were further analysed with ANNOTATOR  indicating their small size with considerable amount of low complexity and predicted disordered regions. The sequence boundaries for the modules and domains of the polyketides were based on the antiSMASH predictions. Protein and domain sequences were aligned with different MAFFT algorithms  and phylogenetic trees were generated with Maximum Likelihood method based on the JTT matrix-based model  with MEGA7 .
Construction of genome editing plasmids
All DNA manipulations were carried out in E. coli DH5α or OmniMAX™ (Thermo Fisher). Primers used in this study are listed in Additional file 1: Table S4. Restriction enzymes were obtained from New England Biolabs. Protospacers were first inserted into pCRISPomyces-2 plasmids  via BbsI-mediated Golden Gate assembly before the respective homology flanks were inserted via Gibson assembly. The protocol for plasmid construction was previously described .
Promoter knock-in constructs are transformed into conjugating E. coli strains and colonies were picked into LB with 50 mg/L apramycin (Sigma). A diaminopimelic acid (DAP) negative ET12567/pUZ8002 donor strain was used . Overnight cultures were diluted 1:100 into fresh LB with apramycin and grown to an OD600 of 0.4–0.6. 400 µL of the culture was pelleted, washed twice and resuspended in LB without apramycin. The washed E. coli cells were then mixed with spores at 1:5 volume ratio and spotted on R2 without sucrose plates. After incubation for 16–20 h at 30 °C, the plates were flooded with nalidixic acid and apramycin and incubated until exconjugants appear. Exconjugants were streaked onto MGY plates containing apramycin at 30 °C followed by re-streaking to MGY plates at 37 °C to cure the CRISPR–Cas9 plasmid containing a temperature-sensitive origin of replication. Apramycin-sensitive clones growing at 37 °C were then subjected to validation of promoter knock-in and genome editing as described below.
Validation of engineered strains
Genomic DNA from wild type and exconjugants from the indicated strains were isolated from liquid cultures using the Blood and Tissue DNeasy kit (Qiagen) after pretreating the cells with 20 mg/mL lysozyme for 30 min at 37 °C. PCR was performed using control primers beyond the homology regions or knock-in specific primers (Additional file 1: Table S4) with KODXtreme Taq polymerase (Millipore). Where indicated, PCR products were subjected to digest with specific restriction enzymes to differentiate between PCR products of wild type genomic sequences and successful genome editing by knock-ins. Positive samples were validated by Sanger sequencing. Oligos used for PCR and sequencing are listed in Additional file 1: Table S4.
Quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) analysis of notonesomycin biosynthetic gene cluster
Total RNA from wild type and engineered Streptomyces sp. A793 were isolated using RNeasy Midi Kit (Qiagen) after 9 days growth and treated with TURBO™ DNase (Ambion). The quality of purified RNA was determined on an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). DNA free RNA samples were used for first strand cDNA synthesis with random hexamers using iScript™ (Biorad) reverse transcription supermix. RT-qPCR was performed on Applied Biosystems™ QuantStudio™ 6 Flex Real-Time PCR System using SYBR FAST qPCR master mix (KAPA). Analysis was carried out with QuantStudio™ Real-Time PCR Software, Version 1.3. The relative abundance of transcripts was calculated by delta-delta CT relative quantitation method/comparative threshold cycle (CT) method . The average CT of housekeeping rpoB, hrdB, recA and atpB genes of Streptomyces sp. A793 was used as endogenous reference to normalize gene expression of each target gene (Additional file 1: Table S5). Triplicates of three biological repeats were performed per condition.
Mass spectroscopy analysis of gene deletion and insertion mutants
A 10 mL culture was freeze dried and extracted with equal volume of methanol and sonicated for 10 min. A wet equivalent of the extract was centrifuged and injected into LC-QTOF for analysis. The HRESIMS spectra were acquired on Agilent UHPLC 1290 Infinity coupled to Agilent 6540 accurate-mass quadrupole time-of-flight (QTOF) mass spectrometer equipped with a splitter and an ESI source. The analysis was performed with a C18 4.6 × 75 mm, 2.7 µm column at flowrate of 2 mL/min, under standard gradient condition of 100% water with 0.1% formic acid to 100% acetonitrile with 0.1% formic acid over 15 min. The typical QTOF operating parameters were as follows: positive ionization mode; sheath gas nitrogen, 12 L/min at 295 °C; drying gas nitrogen flow, 8 L/min at 275 °C; nebulizer pressure, 30 psig; nozzle voltage, 1.5 kV; capillary voltage, 4 kV. Lock masses in positive ion mode: purine ion at m/z 121.0509 and HP-921 ion at m/z 922.0098.
Pairwise comparisons between means abundances of the compounds detected from the mutant and wild type strains were carried out using a Student t-test (two tailed, unpaired) where the null hypothesis states that the means are equal. The difference between the mutant and wild type strains is considered statistically significant if the Student t-test gives a P value < 0.05.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files]. Cluster sequence has been deposited in NCBI GenBank under accession number MN586816.
Biosynthetic gene clusters
Minimum inhibition concentration
- IC50 :
Half maximal inhibitory concentration
- IC90 :
90% maximal inhibitory concentration
Extracted ion chromatograph
Streptomyces antibiotic regulatory protein
Natural Organism Library
Open reading frame
Newman DJ, Cragg GM. Natural products as sources of new drugs from 1981 to 2014. J Nat Prod. 2016;79:629–61.
Bérdy J. Thoughts and facts about antibiotics: where we are now and where we are heading. J Antibiot (Tokyo). 2012;65:385–95.
Parkinson EI, Bair JS, Nakamura BA, Lee HY, Kuttab HI, Southgate EH, et al. Deoxynybomycins inhibit mutant DNA gyrase and rescue mice infected with fluoroquinolone-resistant bacteria. Nat Commun. 2015;6:6947.
Rodrigues T, Reker D, Schneider P, Schneider G. Counting on natural products for drug design. Nat Chem. 2016;8:531–41.
Ng SB, Kanagasundaram Y, Fan H, Arumugam P, Eisenhaber B, Eisenhaber F. The 160 K Natural Organism Library, a unique resource for natural products research. Nat Biotechnol. 2018;36:570–3.
Yoganathan K, Cao S, Crasta SC, Aitipamula S, Whitton SR, Ng S, et al. Microsphaerins A–D, four novel benzophenone dimers with activity against MRSA from the fungus Microsphaeropsis sp. Tetrahedron. 2008;64:10181–7.
Rho J-R, Subramaniam G, Choi H, Kim E-H, Ng SP, Yoganathan K, et al. Gargantulide A, a complex 52-membered macrolactone showing antibacterial activity from Streptomyces sp. Org Lett. 2015;17:1377–80.
Sirota FL, Goh F, Low K-N, Yang L-K, Crasta SC, Eisenhaber B, et al. Isolation and Identification of an anthracimycin analogue from Nocardiopsis kunsanensis, a halophile from a saltern, by genomic mining strategy. J Genomics. 2018;6:63–73.
Yoganathan K, Yang L, Rossant C, Huang Y, Ng S, Butler MS, et al. Cochlioquinones and epi-cochlioquinones: antagonists of the human chemokine receptor CCR5 from Bipolaris brizae and Stachybotrys chartarum. J Antibiot (Tokyo). 2004;57:59–63.
Vallabhaneni S, Mody RK, Walker T, Chiller T. The global burden of fungal diseases. Infect Dis Clin North Am. 2016;30:1–11.
Guinea J. Global trends in the distribution of Candida species causing candidemia. Clin Microbiol Infect. 2014;20(Suppl 6):5–10.
Paulussen C, Hallsworth JE, Álvarez-Pérez S, Nierman WC, Hamill PG, Blain D, et al. Ecology of aspergillosis: insights into the pathogenic potency of Aspergillus fumigatus and some other Aspergillus species. Microb Biotechnol. 2016;10:296–322.
Brown GD, Denning DW, Gow NAR, Levitz SM, Netea MG, White TC. Hidden killers: human fungal infections. Sci Transl Med. 2012;4:165rv13.
Sun S, Lui Q, Han L, Ma Q, He S, Li X, et al. Identification and characterization of Fusarium proliferatum, a new species of fungi that cause fungal keratitis. Sci Rep. 2018;8:4859.
Smith KD, Achan B, Hullsiek KH, McDonald TR, Okagaki LH, Alhadab AA, et al. Increased antifungal drug resistance in clinical isolates of Cryptococcus neoformans in Uganda. Antimicrob Agents Chemother. 2015;59:7197–204.
Sharma C, Kumar R, Kumar N, Masih A, Gupta D, Chowdhary A. Investigation of multiple resistance mechanisms in voriconazole-resistant Aspergillus flavus clinical isolates from a Chest Hospital Surveillance in Delhi, India. Antimicrob Agents Chemother. 2018;62:e01928-17.
Sasaki T, Furihata K, Shimazu A, Seto H, Iwata M, Watanabe T, et al. A novel macrolide antibiotic, notonesomycin A. J Antibiot (Tokyo). 1986;39:502–9.
Chiu H-T, Weng C-P, Lin Y-C, Chen K-H. Target-specific identification and characterization of the putative gene cluster for brasilinolide biosynthesis revealing the mechanistic insights and combinatorial synthetic utility of 2-deoxy-l-fucose biosynthetic enzymes. Org Biomol Chem. 2016;14:1988–2006.
Tanaka Y, Komaki H, Yazawa K, Mikami Y, Nemoto A, Tojyo T, et al. Brasilinolide A, a new macrolide antibiotic produced by Nocardia brasiliensis: producing strain, isolation and biological activity. J Antibiot (Tokyo). 1997;50:1036–41.
Robbins N, Spitzer M, Wang W, Waglechner N, Patel DJ, O’Brien JS, et al. Discovery of ibomycin, a complex macrolactone that exerts antifungal activity by impeding endocytic trafficking and membrane function. Cell Chem Biol. 2016;23:1383–94.
Helaly SE, Kulik A, Zinecker H, Ramachandaran K, Tan GYA, Imhoff JF, et al. Langkolide, a 32-membered macrolactone antibiotic produced by Streptomyces sp. Acta 3062. J Nat Prod. 2012;75:1018–24.
Wan Z, Fang W, Shi L, Wang K, Zhang Y, Zhang Z, et al. Novonestmycins A and B, two new 32-membered bioactive macrolides from Streptomyces phytohabitans HBERC-20821. J Antibiot (Tokyo). 2015;68:185–90.
Salcedo RG, Olano C, Gómez C, Fernández R, Braña AF, Méndez C, et al. Characterization and engineering of the biosynthesis gene cluster for antitumor macrolides PM100117 and PM100118 from a marine actinobacteria: generation of a novel improved derivative. Microb Cell Factories. 2016;15:44.
Salcedo RG, Olano C, Fernández R, Braña AF, Méndez C, de la Calle F, et al. Elucidation of the glycosylation steps during biosynthesis of antitumor macrolides PM100117 and PM100118 and engineering for novel derivatives. Microb Cell Factories. 2016;15:187.
Blin K, Medema MH, Kottmann R, Lee SY, Weber T. The antiSMASH database, a comprehensive database of microbial secondary metabolite biosynthetic gene clusters. Nucleic Acids Res. 2017;45:D555–9.
Gomez-Escribano JP, Alt S, Bibb MJ. Next generation sequencing of actinobacteria for the discovery of novel natural products. Mar Drugs. 2016;14:78.
Zhang MM, Wong FT, Wang Y, Luo S, Lim YH, Heng E, et al. CRISPR–Cas9 strategy for activation of silent Streptomyces biosynthetic gene clusters. Nat Chem Biol. 2017;13:607–9.
Tong Y, Charusanti P, Zhang L, Weber T, Lee SY. CRISPR–Cas9 based engineering of actinomycetal genomes. ACS Synth Biol. 2015;4:1020–9.
Jia H, Zhang L, Wang T, Han J, Tang H, Zhang L. Development of a CRISPR/Cas9-mediated gene-editing tool in Streptomyces rimosus. Microbiol Read Engl. 2017;163:1148–55.
Cobb RE, Wang Y, Zhao H. High-efficiency multiplex genome editing of Streptomyces species using an engineered CRISPR/Cas system. ACS Synth Biol. 2015;4:723–8.
Huang H, Zheng G, Jiang W, Hu H, Lu Y. One-step high-efficiency CRISPR/Cas9-mediated genome editing in Streptomyces. Acta Biochim Biophys Sin. 2015;47:231–43.
Takeuchi T, Hatano M, Umekita M, Hayashi C, Wada S-I, Nagayoshi M, et al. ATP depletion assay led to the isolation of new 36-membered polyol macrolides deplelides A and B from Streptomyces sp. MM581-NF15. Org Lett. 2017;19:4207–10.
Komatsu K, Tsuda M, Tanaka Y, Mikami Y, Kobayashi J. Absolute stereochemistry of immunosuppressive macrolide brasilinolide A and its new congener brasilinolide C. J Org Chem. 2004;69:1535–41.
Sasaki T, Furihata K, Nakayama H, Seto H, Otake N. The structure of a novel macrolide antibiotic, notonesomycin A. Tetrahedron Lett. 1986;27:1603–6.
Hideyuki S, Yasushi T, Katsukiyo Y, Yuzuru M, Jun’ichi K. Brasilinolide A, new immunosuppressive macrolide from actinomycete Nocardia brasiliensis. Tetrahedron. 1996;52:9031–4.
Wan Z, Fang W, Shi L, Wang K, Zhang Y, Zhang Z, et al. Novonestmycins A and B, two new 32-membered bioactive macrolides from Streptomyces phytohabitans HBERC-20821. J Antibiot. 2015;68:185–90.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42:D206–14.
Weber T, Blin K, Duddela S, Krug D, Kim HU, Bruccoleri R, et al. AntiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 2015;43:W237–43.
NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2018;46:D8–13.
Salcedo RG, Olano C, Gómez C, Fernández R, Braña AF, Méndez C, de la Calle F, Salas JA. Characterization and engineering of the biosynthesis gene cluster for antitumor macrolides PM100117 and PM100118 from a marine actinobacteria: generation of a novel improved derivative. Microb Cell Factories. 2016;15–62:44.
Wang L, White RL, Vining LC. Biosynthesis of the dideoxysugar component of jadomycin B: genes in the jad cluster of Streptomyces venezuelae ISP5230 for L-digitoxose assembly and transfer to the angucycline aglycone. Microbiol Read Engl. 2002;148:1091–103.
Pandey RP, Gurung RB, Parajuli P, Koirala N, Tuoi LT, Sohng JK. Assessing acceptor substrate promiscuity of YjiC-mediated glycosylation toward flavonoids. Carbohydr Res. 2014;393:26–31.
Luo Y, Zhang L, Barton KW, Zhao H. Systematic identification of a panel of strong constitutive promoters from Streptomyces albus. ACS Synth Biol. 2015;4:1001–10.
Leyh TS, Taylor JC, Markham GD. The sulfate activation locus of Escherichia coli K12: cloning, genetic, and enzymatic characterization. J Biol Chem. 1988;263:2409–16.
Zhao Q, He Q, Ding W, Tang M, Kang Q, Yu Y, et al. Characterization of the azinomycin B biosynthetic gene cluster revealing a different iterative type I polyketide synthase for naphthoate biosynthesis. Chem Biol. 2008;15:693–705.
Mougous JD, Senaratne RH, Petzold CJ, Jain M, Lee DH, Schelle MW, et al. A sulfated metabolite produced by stf3 negatively regulates the virulence of Mycobacterium tuberculosis. Proc Natl Acad Sci USA. 2006;103:4258–63.
Sogi KM, Holsclaw CM, Fragiadakis GK, Nomura DK, Leary JA, Bertozzi CR. Biosynthesis and regulation of sulfomenaquinone, a metabolite associated with virulence in Mycobacterium tuberculosis. ACS Infect Dis. 2016;2:800–6.
Sekowska A, Kung HF, Danchin A. Sulfur metabolism in Escherichia coli and related bacteria: facts and fiction. J Mol Microbiol Biotechnol. 2000;2:145–77.
Pootoolal J, Thomas MG, Marshall CG, Neu JM, Hubbard BK, Walsh CT, et al. Assembling the glycopeptide antibiotic scaffold: the biosynthesis of A47934 from Streptomyces toyocaensis NRRL15009. Proc Natl Acad Sci USA. 2002;99:8962–7.
Hong H, Samborskyy M, Usachova K, Schnatz K, Leadlay PF. Sulfation and amidinohydrolysis in the biosynthesis of giant linear polyenes. Beilstein J Org Chem. 2017;13:2408–15.
Tang X, Eitel K, Kaysser L, Kulik A, Grond S, Gust B. A two-step sulfation in antibiotic biosynthesis requires a type III polyketide synthase. Nat Chem Biol. 2013;9:610–5.
Li R, Lloyd EP, Moshos KA, Townsend CA. Identification and characterization of the carbapenem MM 4550 and its gene cluster in Streptomyces argenteolus ATCC 11009. Chembiochem Eur J Chem Biol. 2014;15:320–31.
van der Horst MA, Hartog AF, El Morabet R, Marais A, Kircz M, Wever R. Enzymatic sulfation of phenolic hydroxy groups of various plant metabolites by an arylsulfotransferase. Eur J Org Chem. 2015;2015:534–41.
Banik JJ, Craig JW, Calle PY, Brady SF. Tailoring enzyme-rich environmental DNA clones: a source of enzymes for generating libraries of unnatural natural products. J Am Chem Soc. 2010;132:15661–70.
Schniete JK, Cruz-Morales P, Selem-Mojica N, Fernández-Martínez LT, Hunter IS, Barona-Gómez F, et al. Expanding primary metabolism helps generate the metabolic robustness to facilitate antibiotic biosynthesis in streptomyces. mBio. 2018;9:e02283-17.
Romero-Rodríguez A, Rocha D, Ruiz-Villafan B, Tierrafría V, Rodríguez-Sanoja R, Segura-González D, et al. Transcriptomic analysis of a classical model of carbon catabolite regulation in Streptomyces coelicolor. BMC Microbiol. 2016;16:77.
Jin X-M, Chang Y-K, Lee JH, Hong S-K. Effects of Increased NADPH concentration by metabolic engineering of the pentose phosphate pathway on antibiotic production and sporulation in Streptomyces lividans TK24. J Microbiol Biotechnol. 2017;27:1867–76.
Wentzel A, Bruheim P, Øverby A, Jakobsen ØM, Sletta H, Omara WAM, et al. Optimized submerged batch fermentation strategy for systems scale studies of metabolic switching in Streptomyces coelicolor A3(2). BMC Syst Biol. 2012;6:59.
Wilson K. Preparation of genomic DNA from bacteria. Curr Protoc Mol Biol. 2001;Chapter 2:Unit 2.4.
Chin C-S, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563–9.
Labana P, Gosse JT, Boddy CN. Draft genome sequence of the type strain Streptomyces armeniacus ATCC 15676. Microbiol Resour Announc. 2018;7:e01107-18.
Hyatt D, Chen G-L, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.
Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27:4636–41.
Salzberg SL, Delcher AL, Kasif S, White O. Microbial gene identification using interpolated Markov models. Nucleic Acids Res. 1998;26:544–8.
Lukashin AV, Borodovsky M. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.
Eisenhaber B, Kuchibhatla D, Sherman W, Sirota FL, Berezovsky IN, Wong W-C, et al. The recipe for protein sequence-based function prediction and its implementation in the ANNOTATOR software environment. Methods Mol Biol Clifton NJ. 2016;1415:477–506.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci CABIOS. 1992;8:275–82.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Lim YH, Wong FT, Yeo WL, Ching KC, Lim YW, Heng E, et al. Auroramycin: a potent antibiotic from Streptomyces roseosporus by CRISP–Cas9 activation. ChemBioChem. 2018;19:1716–9.
Allard N, Garneau D, Poulin-Laprade D, Burrus V, Brzezinski R, Roy S. A diaminopimelic acid auxotrophic Escherichia coli donor provides improved counterselection following intergeneric conjugation with actinomycetes. Can J Microbiol. 2015;61:565–74.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods San Diego Calif. 2001;25:402–8.
The authors gratefully acknowledge the late Sydney Brenner for his discussion and support. We also thank Raquel Peh for her technical assistance.
This work was supported by the Agency for Science, Technology and Research (A*STAR), Singapore. WLY, EH, MMZ, FTW are supported by National Research Foundation, Singapore [NRF2013-THE001-094].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
NMR spectral dataa of notonesomycin A (1) and notonesomycin B (2). Table S2. Annotation of 59 ORFs of the notonesomycin BGC by sequence homology searches using BLASTp against the NCBI nr database. Table S3. Media compositions used in shake flask fermentation. Table S4. Oligos used in this study. Table S5. List of primers used for RT-qPCR. Figure S1. Structure and gross structures of notonesomycins. Selected HMBC correlations of 1 and 2. Figure S2.13C NMR spectrum (methanol-d4, 100 MHz) of Notonesomycin A (1). Figure S3.1H NMR spectrum (methanol-d4, 400 MHz) of Notonesomycin A (1). Figure S4. COSY spectrum (methanol-d4, 400 MHz) of Notonesomycin A (1). Figure S5. HSQC spectrum (methanol-d4, 400 MHz) of Notonesomycin A (1). Figure S6. HMBC spectrum (methanol-d4, 400 MHz) of Notonesomycin A (1). Figure S7.13C NMR spectrum (methanol-d4, 100 MHz) of Notonesomycin B (2). Figure S8.1H NMR spectrum (methanol-d4, 400 MHz) of Notonesomycin B (2). Figure S9. COSY spectrum (methanol-d4, 400 MHz) of Notonesomycin B (2). Figure S10. HSQC spectrum (methanol-d4, 400 MHz) of Notonesomycin B (2). Figure S11. HMBC spectrum (methanol-d4, 400 MHz) of Notonesomycin B (2). Figure S12. Phylogenetic analysis of 59 acyltransferases (AT) amino acid sequences from notonesomycin (Ntc AT), brasilinolides (Nbr AT), PM100117 and PM100118 (GonP AT) BGCs. Figure S13. Proposed biosynthesis pathways for A) 4-amino 3-hydroxybenzoic acid and B) deoxysugars. Figure S14. Sequence analysis of glycosyltransferases Nbc18 and Nbc22 from the notonesomycin BGC with known glycosyltransferases. Figure S15. Sequence analysis of cytochrome P450 enzymes (Nbc33, Nbc36 and Nbc37) from the notonesomycin BGC with other cytochrome P450 with either hydroxylation or epoxidation function. Figure S16. Sequence analysis of methyltransferases Nbc1 and Nbc56 from Streptomyces sp. A793 present in the notonesomycin BGC with other known O-, N- and C-methyltransferases. Figure S17. LC–MS spectra (base peak chromatogram, BPC) for wild type and nbc20 and nbc21 double KO mutant. Figure S18. Sequence analysis of sulfotransferase Nbc48 from the notonesomycin BGC with sulfotransferases that have been associated with other BGCs.
About this article
Cite this article
Goh, F., Zhang, M.M., Lim, T.R. et al. Identification and engineering of 32 membered antifungal macrolactone notonesomycins. Microb Cell Fact 19, 71 (2020). https://doi.org/10.1186/s12934-020-01328-x