放线菌样品的振荡

靶向基因组挖掘-从化合物的发现到生物合成途径的阐明Targeted Genome Mining-From Compound Discovery to Biosynthetic Pathway Elucidation靶向基因组挖掘-从化合物的发现到生物合成途径的阐明 2 of 17Microorganisms 2020, 8, 2034 Article Targeted Genome Mining-From Compound Discovery to Biosynthetic Pathway Elucidation Nils Gummerlich 1, Yuriy Rebets 1, Constanze Paulus 1, Josef Zapp’and Andriy Luzhetskyy 1,3,* 1 Department of Pharmaceutical Biotechnology, Saarland University, Campus C2.3, 66123 Saarbrucken, Germany; nils.gummerl i ch@gmx.de (N.G.); yurko.rebets@gmail.com (Y.R.); constanzepaulus@gmai l .com (C.P.) 2 Department of Pharmaceutical Biology, Saarland University, Campus C2.3, 66123 Saarbrucken, Germany; j.zapp@mx.uni-saarland.de 3 Actinobacteria Metabolic Engineering Group, Helmholtz-Institute for Pharmaceutical Research Saarland (HIPS), Campus E8.1, 66123 Saarbrucken, Germany 米 Correspondence: a.luzhetskyy@mx.uni-saarland.de; Tel .: +49-681-302-70200 （HIPS）放线菌代谢工 Received: 29 November 2020; Accepted: 15 December 2020; Published: 19 December 2020 check for updates 程组 Abstract: Natural products are an important source of nove l investigational compounds in drug discovery. Especially i n the f ield of antibiotics, Actinobacteria have been proven to be a reliable source for lead structures. The discovery of these natural products with activity- and structure-guided screenings has been impeded by the constant rediscovery of previously identified compounds. Additionally, a large discrepancy between produced natural products and biosynthetic potential in Actinobacteria, including representatives of the order Pseudonocardiales, has been revealed using genome sequencing.To turn this genomic potential into novel natural products, we used an approach including the i n-silico pre-selection of unique biosynthetic gene clusters followed by their systematic heterologous expression. As a proof of concept, fifteen Saccharothrix espanaensis genomic library clones covering predicted biosynthetic gene clusters were chosen for expression i n two heterologous hosts, Streptomyces lividans and Streptomyces albus . As a result , two novel natural products, an unusual angucyclinone pentangumycin and a new type II polyketide synthase shunt product SEK90, were identified. After puri f ication and structure elucidation, the biosynthetic pathways leading to the formation of pentangumycin and SEK90 were deduced using mutational analysis of the biosynthetic gene cluster and f eeding experiments with 13C-labelled precursors. Keywords: angucyc l inone; heterologous expression; natural products; Actinobacteria; biosynthesis 1. Introduction Actinobacteria,especially Streptomyces, represent one of the most important sources of natural products [1]. During recent decades, their discovery in Actinobacteria was significantly slowed due to re-isolation of known compounds in activity- and structure-guided screenings [2,3]. As an early solution for the rediscovery i ssue, natural product databases have been created. Metabolic profiles recorded with high performance l iquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) can be compared with information stored in these databases. This process, known as dereplication, is problematic due to it being time consuming and error prone [4]. An additional problem of this “compound f irst" approach i s the need for laborious studies in order to identify the biosynthetic gene clusters (BGCs) responsible for the production of the isolated compounds. This information is essential to elucidate the biosynthetic pathway, which helps researchers to increase and modify the production of identified natura l products [5]. In contrast, the "cluster to compound" approach aims to obtain new natural products after cloning and expression of a BGC of i nterest in heterologous hosts. This approach became increasingly reasonable due to the low cost of next-generation sequencing [6], which has resulted in a flood of genomic information [7]. It has been shown that Streptomyces contain an average of 21.9 BGCs per genome, with 40-48% of these BGCs being unique [8,9]. The "cluster to compound " approach aims to tap into this genomic potential and offers decisive advantages over the “compound f irs t " approach. First , genome analysis and BGC prediction tools such as antiSMASH [10] offer automated cluster homology comparisons within hours. In turn, this allows cataloguing and Microorganisms 2020, 8, x FOR PEER REVIEW 2 of 18of a BGC of interest in heterologous hosts. This approach became increasingly reasonable due to the low cost of next-generation sequencing [6], which has resulted in a flood of genomic information [7]. It has been shown that Streptomyces contain an average of 21.9 BGCs per genome, with 40–48% of these BGCs being unique [8,9]. The “cluster to compound” approach aims to tap into this genomic potential and offers decisive advantages over the “compound first” approach. First, genome analysis and BGC prediction tools such as antiSMASH [10] offer automated cluster homology comparisons within hours. In turn, this allows cataloguing and prioritization of BGCs with unique features. Second, the utilization of optimized heterologous hosts and the expression of BGCs from genomic librariesofferadvantagesover use nativeproducingstrains[11–14].Furthermore, produced compounds can be identified by a simple comparison of the metabolic profiles from the heterologous host with and without the expressed BGCs [11]. The aim of the presented study was to assess the untapped potential for the production of secondary metabolites hidden in the genomes of Actinobacteria and to identify their biosynthetic pathways [8]. To achieve this goal, a widely applicable workflow was developed (Figure 1). This workflow includes the analysis of genome sequences utilizing algorithms such as antiSMASH [15] PRISM[16],theconstructionof sequencedgenomiclibrary,insilicopre-selection expression of the selected BGCs in optimized heterologous hosts followed by the identification, isolation and characterization of natural products and their biosynthetic pathways. To challenge the applicability of the designed approach, we aimed to express BGCs from a strain of Actinobacteria that is more distantly related to commonly used heterologous Streptomyces hosts. The large-scale analysis by Doroghazi and co-authors revealed that the order of Pseudonocardiales has an average of 19.8 BGCs per genome [8]. The strain Saccharothrix espanaensis, discovered by Labeda and Lechevalier in 1989 [17], belongs to the order of Pseudonocardiales, harboring all of the desired elements to prove the potential of our approach. This strain has a 9.3 Mbp genome containing at least 31 antiSMASH-predicted BGCs [18]. Due to the strain being poorly genetically tractable, any genetic manipulations by S. espanaensis a .realmostimpossible.Furthermore,onlyasinglegroupofpolysaccharidenaturalproducts, saccharomicins [19], is known to be produced by S. espanaensis. Figure 1. Workflow for the identification of novel natural products using genome mining. As proof of concept, we report the discovery of two novel natural products which were obtained by applying the proposed genome mining approach. Neither of the isolated compounds were detectable in the metabolic profile of S. espanaensis. Additionally, we were able to elucidate the biosynthesis of both compounds using the available information about the BGC architecture, mutational analysis and feeding experiments wi t h labelled precursors. 2. Materials and Methods 2.1. General Experimental Procedures All oligonucleotides, plasmids and strains used in this work can be found in Tables S1 and S2. Standard procedures for Escherichia coli (transformation and plasmid preparations) were performed as described by Sambrook [20]. All BACs of the genomic l ibrary of S. espanaensis were i solated with the BACMAXIM DNA purif i cation kit (Lucigen, Middleton, WI, USA). All E. coli cultures were cultivated at 37°C at 180 rpm in LB-Lennox (5 g yeast extract, 10 g tryptone,5 g NaCl and 1000 mL H2Odest;pH=7.2). For DNA fragment amplification, Dream Taq polymerase (Thermo Fisher Scientific, Wal t ham, MA, USA) was used. DNA fragments were purified from agarose gels using the QIAquick Gel Extraction Kit (Qiagen, Germany). Intergeneric conjugation between Streptomyces and E. coli was performed as described by Kieser [21] on MS agar (20 g soya flour, 20 g agar, 20 g mannitol and 1000 mL H2O with nalidixic acid 30 ug/mL for donor suppression) plates using E . coli PUB307 as donor strain. Standard DNA manipulations (ligation, polymerase chain reactions and endonuclease digestion) were performed according to the corresponding manufacturer’s protocol. All Streptomyces were cultivated in 15 mL TSB (17 g Tryptone, 3 g Peptone, 5 g NaCl,2.5g K2HPO4, 2.5 g glucose and 1000 mL H2Odesti pH=7.2) i n a 100-mL f lask (4 baffles; 5 g glass beads), inoculated from spores of the corresponding Streptomyces from an MS plate. After 24 h for S. albus and 72 h for S. lividans,50 mL of the corresponding production medium (SG medium (20 g glucose,5 g yeast extract, 10 g bactosoytone, 2 g CaCO3 and 1000 mL H2Odest;pH=7.2); DNPM medium (dextrin 40 g, soytone 7.5 g, f resh yeasts 5 g, MOPS 21 g and 1000 mL H2Odest; pH=6.8); TSB medium; SGG medium (10 g starch, 10 g glycerol , 2.5 g corn steep solids, 5 g Peptone, 2 g yeast extract,1 g NaCl, 3 g CaCO3 and 1000 mL H2Odest;pH=7.2)) was inoculated with 1 mL of the seed culture in a 500-mL flask (3 baffles, 10 g glass beads) and incubated for 72-164 h at 180 rpm and 29°C. Production optimization steps were carried out by adjustment of different parameters (growth time, pH, medium, extrac t ing solvent, addition of XAD -16 Amberlite (2 g/100 mL) after 72 h). The optimization steps were monitored using HPLC-MS. When needed, media were supplemented with the following antibiotics: apramycin 50 ug/mL, hygromycin 50 ug/mL, kanamycin 100 ug/mL and nalidixic acid 30 ug/mL. All HPLC-MS spectra were recorded on a Dionex Ultimate 3000 (Thermo Fisher Scientific, Waltham, MA, USA) coupled with an AmaZon ETD SL speed, Apollo II ESI (Bruker, Billerica, MA, USA) on a Waters BEH C18 column (100×2.1 mm, 1.7 um) with an 18-min linear gradient from 5%acetoni t rile (0.1% formic acid) to 95% acetonitrile (0.1% formic acid). The mass spectra were recorded in centroid mode (200 to 2000 m/z) at a scan rate of 2 Hz. Prior to analysis, all samples were dissolved in methanol and centrifuged for 10 min at 4 °C and 15,000 rpm. The data were analyzed using Bruker Compass Data Analysis 4.2. High-resolution masses were recorded on a Dionex Ultimate 3000RSLC HPLC (Thermo F i sher Scientific) coupled with an LTQ Orbitrap (Thermo Fisher Scientific) on a BEH C18 100×2.1 mm, 1.7-um column (Waters) with an 18-min linear gradient from 5% acetonitrile (0.1% formic acid) to 95% acetonitrile (0.1% formic acid). Data analysis was carried out using Xcalibur 3.0. Preparative HPLC was carried out using a Waters Autopurification System equipped with a Waters 2545 Binary Gradient module, a Waters System Fluidics Organizer (SFO), a Waters 2998 Photodiode Array Detector (PAD), a Waters 2767 Sample Manager and a Waters SQ-Detector-2 on Nucleodur C18Htec 250/21-C18 5 um column (Macherey-Nagel,Duren, Germany). Data analysis was carried out using MassLynx. Structural elucidation was carried out using nuclear magnetic resonance (NMR) technology. The purified compounds were solved in 300 uL deuterated solvent (DMSO-d6, MeOD-d4, CDCl3) and measured in a corresponding 5-mm Shigemi tube (DEUTERO GMBH; Kastellaun, Germany). NMR data (H, HH-COSY; TOCSY; HMBC, HSQC,13C) were acquired either on a Bruker Ascend 700spectrometer equipped with a 5-mm TXA Cryoprobe or a Bruker Avance 500 spectrometer equipped with a 5-mm BBO Probe at 300 K (Bruker, BioSpin, GmbH, Rheinstetten, Germany). The data were analyzed using Bruker’s TopSpin 3.5a software. 2.2. In-Silico Analysis of Gene Clusters The genome of Saccharothrix espanaensis was analyzed with antiSMASH [15] and the Geneious software. Identified BGCs were mapped to the BAC l ibrary. A detailed analysis of 1E5's BGC was performed with BLAST [22]. 2.3. Extrac t ion, Dereplication and Isolation The biomass was separated from the production medium (SG medium (20 g glucose, 5 g yeast extrac t , 10 g bactosoytone,2g CaCO3 and 1000 mL H2Odest;pH=7.2), 172 h, 29°℃, 180 rpm,pH=7.2,1% XAD-16 after 24 h) by centrifugation. Then, 20 mL of the supernatant was extracted with 20 mL ethyl acetate followed by 20 mL butanol for 20 min using a Laboshake LS500 (C. Gerhardt GmbH, Konigswinter, Germany) at 160 rpm. The biomass was extracted with a 1:1 mixture of acetone and methanol for 60 min. The solvents were evaporated unti l dry with a rotary evaporator (150 rpm, 60℃,240 mbar for ethyl acetate or 25 mbar for butanol ) or under a nitrogen stream at 40°C. The extracts were analyzed using LC-MS and LC-HRMS. Singular peaks were compared to the databases Dictionary of Natural Products, version 27.1 (CRC Press, Cleveland, OH,USA), Sci-Finder (CAS, Columbus, OH, USA) and Chemspider (Royal Society of Chemistry, Raleigh,NC,USA). For production, 15 L of SG medium inoculated with the corresponding recombinant strain was cultivated. The extracts were evaporated to dryness and dissolved in 10 mL methanol. To purify single compounds, size exclusion chromatography (SEC) was performed. The column was packed with 600 mL Sephadex LH-20 resin (GE Healthcare Europe GmbH, 79,111 Freiburg, Germany) solved i n methanol. At a flow rate of 10 mL per 15 min, a fraction collector was used to separate the extract within 24 h. To identify fractions containing the novel substances, every 4th fraction was analyzed by HPLC-MS. Fractions with the corresponding mass were combined and evaporated. Prior to NMR analysis, a preparative HPLC was performed . Preparative HPLC information can be found in Tables S3 and S4. NMR data of pentangumycin (1) and SEK90 (2) can be found in Table 1. Table 1. NMR spec t roscopic data for pentangumycin (1)(DMSO-d6)and SEK90 (2) (DMSO-d6). Pentangumycin (1) Pos. 8c 8H (J in Hz) HMBC SEK90 (2) Pos. 8c 8H (J in Hz) HMBC 1 156.7C 1 166.5C 一 一 2 108.7 CH 7.10s 1,4,12b,19 2 102.7C 一 一 3 143.2C - 3 165.5C 一 4 117.5 CH 7.44s 1,2,4a,5,12b,19 4 101.9 CH 5.68 (s br) 2,5,6 4a 136.6C 一 一 5 160.6C 一 一 5 135.1 CH 7.88s 1, 4,4a,6a,7,12a, 12b,13 6 36.1 CH2 3.55 (s br) 4,5,7,8,12 6 139.1C 一 7 132.5C 6a 131.2C 8 120.1CH 6.74 dd (7.7,N/A) 9,10,11,12,13 7 188.0C 一 9 130.9 CH 7.2dd (7.9,N/A) 7,8,10,11,12 7a 115.8C 一 10 114.4 CH 6.77 dd (8.1,N/A) 8,9, 11, 12, 13 8 160.5C 一 11 153.3C 一 9 123.0 CH 7.26m 7a,11 12 130.6C 10 136.7 CH 7.73t(7.5) 8,11a 13 199.9C 一 一 11 117.2 CH 7.50dd (7.5,0.93) 7a,9,10,8,12 14 115.3C 一 11a 135.5C 一 15 164.9 C 12 185.5C 一 16 100.5 CH 6.11d(2.4) 14,15,17 12a 140.0C 一 17 166.3 C 一 12b 119.5C 一 18 111.4 CH 6.04 dd (0.64,2.4) 13,14,15,16,20 13 132.2C - 19 142.7C - 一 14/18 129.9 CH 7.26m 6,14/18,16 20 21.2 CH3 1.8s 14,18,19 15/17 114.8 CH 6.78 dt (8.5,2) 13,15/17,16 21 29.5 CH 4.36t(8.2) 1,2,3,22,23 16 156.6 C 22 31.5 CH。 1.75 dt (7.3,8.2) 2,21,23,24 19 42.4 CH 4.43d(6) 2,3,4,20 23 20.6 CH 1.05dq (7.3) 21,22,24 20 169.6C 一 24 13.8 CH3 0.89t(7.3) 22,23 21 22.8 CH3 1.93 s 20 一 一 22 56.1 CH3 3.88s 1 一 一 一 一 NH 8.53 t (br) 19,20 一 一 OH1 一 一 一 一 12.67 s 14,15,16 OH2 11.49 s (br) OH3 一 10.38 s 16,18 OH4 一 一 一 一 9.78 s 10, 11,12 2.4. Modification of BACs All constructed single knockout BACs were obtained by homologous recombination, as previously described by Gust [14] and Myronovskyi [23], and confirmed by PCR. The recombinant BACs were transformed into E. coli PUB307 and conjugated into S. lividans AYA6. All recombinant strains were cultivated , and their metabolome was analyzed as described above. 2.5. Construction of Plasmids For the construction of the plasmids pTOS _A3_R1 and pTOS _A3_R2, we amplified penR1and penR2 using PCR and used the ordered oligonucleotides to introduce the sequence of the A3promotor [24] and a KPNI restriction site in the 3'oligonucleotide and a HindII I restric t ion site in the 5' oligonucleotide. The amplified PCR product and the i ntegrative pTOS plasmid were cut using HindIII and KPNI and ligated by the T4 ligase. The obtained DNA fragment was transformed into E. coli PUB307 and conjugated i nto S. lividans △YA6. The backbone of pTOS was removed from the genome by expression of pUWL _DRE containing the Dre recombinase. Synthesis of 1E5_CMP1 and 1E5_CMP2 For the synthesis of 1E5_CMP1 and 1E5_CMP2, 6-benzyl-4-Hydroxy-2-pyrone (Wuxi AppTec, Saint Paul , MN, USA) was used under the conditions described by de March. The pyrone was mixed (2:1) with either butyraldehyde or formaldehyde in 3 mL ethanol containing 10 uL piperidine and 10 uL glacial acetic acid. The reaction was performed at room temperature for 24 h. All compounds were purified using the Waters Autopurification System and confirmed by NMR. Feeding experiments.To prove the incorporation from L-Phenylalanine, L-Tyrosine and L-Tryptophan, a culture of either S. lividans AYA6_1E5 or S. lividans △YA6_1E5_A3_R1 was prepared as described above. Thereafter, 10 mg of the l abelled amino acid (either L-Phenylalanin, L-Tyrosine or L-Tryptophan) was supplemented in 5 steps: 2 mg after 24, 48, 72, 96 and 120 h. After 164 h of incubation, the supernatant of the production cultures was extracted and the incorporation was quantified with LC-MS. 3. Results and Discussion 3.1. Selection of Putative Biosynthetic Gene Clusters Initially, the biosynthetic potential of S. espanaensis (txid1179773) was assessed. An antiSMASH analysis of its genome predicted the presence of 31 BGCs, while the only known products are saccharomicins (Table S5) [19]. The BGC of saccharomicins was identified as cluster number 28 [18]. It was originally cloned as a cosmid c l one by Berner et al., while attempts at its complete expression in heterologous hosts had failed. Nevertheless, its aglycon was successfully expressed and its biosynthesis was elucidated [25,26]. A bacterial artificial chromosome (BAC) library of S. espanaensis was constructed using the pSMART-BAC-S vector as a backbone with an average inserted fragments size of 100 kb. The library was end-sequenced, and the sequences were mapped to the genome of S. espanaensis. This allowed the alignment of predicted BGCs to the clones of the genomic l ibrary. To minimize the laboratory effort and reduce the putative rediscovery rate of natural products in the described approach, BGCs were prioritized by their class and predicted products. Eight BGCs were not fully covered by the genomic library and, therefore, could not be included into the l aboratory approach. The remaining 23 BGCs of the BAC library were composed of three BGCs that were predicted to produce geosmin, melanin and bacteriocin and were, therefore, excluded. Three BGCs shared more than 75%similarity to clusters of assigned products and were, therefore, not used for heterologous expression. Finally, a predicted terpene cluster lacking i ts synthase and a l antipeptide BGC were excluded as well . Thereafter, 15 remaining BACs covering 17 BGCs were selected for heterologous expression in S. lividans AYA6 and S. albus J1074 (Table S5). Among the 17 BGCs, two were terpene clusters, five were either type I or type II polyketide synthase (PKS) clusters, five were non-ribosomal peptide synthetase (NRPS) clusters, two were type I PKS/NRPS hybrid clusters, one was a lantipeptide cluster, one was an aminoglycoside cluster and one was t he polysaccharide cluster encoding saccharomicins BGC. 3.2. Success Rate of the Expression System The chosen BACs were conjugated into S. lividans AYA6 and S. albus J1074. Every further-used transconjugant was tested with colony polymerase chain reaction (PCR) for the presence of the conjugated BAC. Verified transconjugants were cultivated in different production media and metabolites were extracted using different solvents. Their metabolic profiles were analyzed with LC-MS and compared to the metabolic profiles of the empty heterologous hosts. Singular peaks were identified and assumed to be the products, shunt products or intermediates of the expressed BGC. Extracts containing singular peaks were further analyzed by LC-HRMS, and the obtained exact masses were compared to masses in common natural product databases (Dictionary of Natural Products (http://dnp.chemnetbase.com CRC Press, Cleveland, OH, USA), Supernatural [http://bioinf-applied .charite. de/supernatural _new/index.php Charite University of Medicine, Berlin, Germany], StreptomeDB [http://132.230.56.4/streptomedb/University of Freiburg, Freiburg, Germany) [27,28]. The recombinant strain S. lividans △YA6_3C18, containing a 115-kb insert of the S. espanaensis genome, revealed a singular peak (Figure S1) with a mass of 308.1987 Da (m/z 309.2060 [M+H]+, Rt = 11.5 min). BAC clone 3C18 carries a type I PKS cluster similar to lavendiols BGC [29] and an unknown NRPS BGC. The low production rate of the produced compound rendered isolation, structure-and biosynthetic elucidation unfeasible. The expression of 3C18 was unsuccessful in S. albus J1074. Microorganisms 2020, 8, x FOR PEER REVIEW 6 of 18S. espanaensis, produced two singular peaks with masses of 467.1368 (m/z 468.1438[M+H]*,Rt = 9.1 min, further designated as pentangumycin (1)) and 790.2261 Da (m/z 791.2339 [M+H]+, Rt = 10.6 min; further designated as SEK90 (2)), as shown in Figure 2. S. albus J1074_1E5 produced compound 1(Figure S2), but its low production titer rendered further work infeasible. Both compounds 1 and 2have not been observed in extracts obtained from S. espanaensis (Figure S3). Figure 2. Structural formulas of pentangumycin (1) and SEK90 (2). In the case of the remaining 13 BAC clones, no production was observed under the tested conditions described in Section 2.1. An overall success rate of 11% (corresponding to 2 out of 17 BGCs, BAC_3C18andBAC_1E5)wasobservedwhenexpressingBGCsfromPseudonocardialesin Streptomyces hosts. In similar experiments with BGCs from a Streptomyces strain, a success rate of 35%(6 out of 17 BGCs) was observed (unpublished data). The generally low expression rate can be explained by a variety of factors: Firstly, differences in the regulatory network of strains can cause major shifts of natural product production. Prominent examples include the mutations in the genes rpoB (C1298T, S443L) and rpsL (A262G, K88E) in the strains S. coelicolor and S. lividans, which activated and increased the production of natural products in both strains [12]. Secondly, codon bias of different strains can be restrictive for the expression of BGCs in heterologous hosts. In fact, it has been shown to be one of the most important factors for the expression of prokaryotic genes [30,31]. In Escherichia coli, the heterologous expression of genes was increased by expanding the intracellular concentration of rare tRNA [32]. In S. lividans, the expression of a transglutaminase was increased by tRNA [32]. In S. lividans , the expression of a transglutaminase was increased by 73.6%, when using a codon optimized gene for expression [33]. Thirdly, promotors and ribosomal binding sites not derived from the native strain can be unfitting for the polymerase of the host strain and, therefore, disrupt the expression of BGCs. The utilization of strong promotors can activate the expression of BGCs in native and host strains. Salas et al. successfully activated the expression of an NRPS and PKS-NRPS hybrid BGC by introducing the strong and constitutive ermE*p promotor in front of the NRPS (sshg _00313) and PKS-NRPS (sshg _05713) genes. The activation of the NRPS BGC was confirmed by the production of a blue pigment. Furthermore, they identified the natural products 6-epi-alteramid A and B, which were produced by the PKS-NRPS BGC [34]. Fourthly, the toxicity of either expressed proteins or produced natural products is an additional issue to keep in mind when expressing BGCs in heterologous hosts. For example, a mutated version of an outer membrane protein (OmpA) of E. coli K12 caused a toxic lysis of the cell when expressed [35]. Furthermore, the toxicity of avermectin was shown to be a restrictive factor in i ts own production in S. avermitilis. After the introduction of multiple copies of the ABC transporter AvtAB, the avermectin production increased from 3.3 to 4.8 g/L [36]. Finally, environmental factors can drastically influence the production of natural products. A well-studied example is the produc t ion of jadomycin, since it was only produced after the strain was either treated with a heat shock or the production medium was supplemented with ethanol [37]. Furthermore, i t has been proven that the addition of sugars [38] and sub-toxic concentrations of antibiotics can activate and strongly alter the expression of BGCs in bacteria [39,40]. 3.3. Isolation of Pentangumycin and SEK90 Ultimately, 3.5 mg of pentangumycin 1 and 14 mg of SEK90 2 were isolated f rom 15 L of a S. lividans △YA6_1E5 production culture. Both compounds were purified using size exclusion chromatography and preparative HPLC. Through analysis of LC-HRMS data and NMR data ( H-NMR,HSQC,HMBC,1H-1H-COSY, 13C-NMR), the structures of both compounds were elucidated (Figure 2). NMR data and the detailed data analysis are presented in Table 1 and Figures S4-S50. All chemical shifts observed for 1 were within expected ranges as calculated for its structure. The connection between rings E and B was observed through HMBC correlations between H-5 and C-13 and between H-14/H-18and C-6 (Figure 2, Figure S4). Compound 1 i s a member of the angucyclinone group containing a unique fifth phenol ring attached through a direct carbon-carbon bond at position C-6. Additional ring systems, typically derived from amino acids, have previously been observed in the angucycline family of aromatic polyketides. While such structures have been found in several jadomycins and urdamycins, their position and type of connection are different [41-43]. Additionally, formicamycins and f asamycins, which are examples of anthracycline antibiotics, carry a similar structural motif with the fifth aromatic ring at position C-7 of ring B [44]. Furthermore, 1 has a methylated hydroxy group at position C-1 in ring A, a modification that can be found only in a few angucyclines, such as the chlorocyclinones [45]. In addition, the drastic struc t ural feature of 1 i s t he presence of an aminated and, subsequently, the acetylated methyl group at position C-19, which has not been previously reported. The biological activity of 1 and 2 was tested. Both compounds showed no significant antimicrobial and anticancer activity. In a CAS assay, compound 2 showed low iron-binding siderophore activity [46]. The NMR spec t ra of 2 contained less carbon and proton signals than anticipated regarding its observed exact mass. The lH-NMR signal of H-21 (CH) i ndicated the presence of a mirror plane, since all other integrated proton signals were a multiple of one (2-6). Another strong indicator was the low number of carbon peaks observed in the 13C-NMR (24 peaks in 13c-NMR), which only sums up to the corresponding exact mass of 2 with an unreasonable number of incorporated heteroatoms. After the identification of the mirror plane, the monomer SEK43 (Figure S51) could easily be elucidated [47]. The protons H-16, H-18 and H-20 as well as H-6, H-8, H-9 and H-10 showed all necessary HMBC and 1H-1H-COSY correlations f or the elucidation of t heir corresponding aromatic system. The connecting ketone moiety was elucidated through the correlations of H-18 and H-10 to C-13 and the carbon shift of C-13 (199.9 ppm). The pyrone moiety of SEK43 was elucidated through correlations between H-6 to C-5 and C-4, H-4 to C-1, C-2, C-4 and C-6, H-21 to C-1 and C-3 and H-22 to C-2. The connecting butyl group showed all expected H-H-COSY and HMBC correlations (Figure S28). The putative axial stereochemistry of 2 can be assumed to be a racemic mixture due to its spontaneous formation as described below. 3.4. Both 1 and 2 Are the Result of Type II PKS Gene Cluster Expression The antiSMASH analysis of the sequence of BAC 1E5 revealed the presence of three putative BGCs: a t erpene cluster, a type-I I PKS and a PKS-NRPS hybrid . From a structural perspective, i t is obvious that neither 1 nor 2 are products of the terpene BGC and that 1 is a product of the type II polyketide synthase. However, 2 cannot easily be associated with the type II PKS or the PKS-NRPS hybrid. Even though the SEK43 moiety of 2 i s a known shunt product of type II polyketide synthases, the connection through a butyl group has never been reported. To verify the biosynthetic origins of 1 and 2, we aimed to identify the core genes of the corresponding BGCs and deleted them from BAC 1E5. A BLAST analysis of the type II PKS BGC showed that gene BN6_54860 is coding for a B-ketoacyl-synthase (Figure 3) and gene BN6_55110 is encoding a type I polyketide synthase in the PKS-NRPS BGC. Both genes were deleted separately i n frame from BAC 1E5. The resulting BACs,1E5△54860 and 1E5▲55110, were conjugated into S. lividans AYA6, and the metabolic profiles of both recombinant strains were analyzed. The production in the strain carrying mutant BAC 1E5^55110remained unchanged, whereas the deletion of gene BN6_54860 resulted in the loss of production of 1and 2 (Figure S52). This clearly proves the origin of both compounds to be the type II PKS gene cluster. Microorganisms 2020, 8, x FOR PEER REVIEW 9 of 18 (a) C D R1 H M N Q S T R2 Cytochrome P450 Monooxygenase >PKS Core Red u c t ase / Oxygenase M e thyltransferase >Cyclase Gene# Name Homology to: Annotation Ident.[%] Cover.[%] BN6_54730 penA WP_124773694 cytochrome P450 71.1 99 BN6_54740 penB WP_124773693 ferredoxin 64.6 98 BN6_54750 penC WP_124773692 aminotransferase class II-fold, pyridoxal phosphate-dependent 59.5 99 BN6_54760 penD WP_124773691 methyltransferase 68.8 99 BN6_54770 penE WP_015037167(jadG)* antibiotic biosynthesis monooxygenase 58.8 90 BN6_54780 penF WP_007904366 ester cyclase 46.3 96 BN6_54790 penG WP_005164912 FAD-dependent oxidoreductase 41.7 90 BN6_54800 penR1 WP_034626020 XRE family transcriptional regulator 57.4 95 BN6_54810 penH WP_124773686 MFS transporter 49.2 89 BN6_54820 penl WP_124776097 aldo/keto reductase 61.1 76 BN6_54830 penJ WP_124776095 nuclear transport factor 2 family protein 75.5 77 BN6_54840 penK WP_015035778 acetyl-/propionyl-CoA carboxylase subunit alpha 69.8 99 BN6_54850 penL AAD13535.1 (lanF) Polyketide cyclase 74.1 98 BN6_54860 penM AAD13536.1(lanA) Beta-ketoacyl-ACP synthase homolog 73.3 99 BN6_54870 penN AAD13537.1(lanB) Polyketide chain length factor 68.2 99 BN6_54880 penO AAD13538.1 (lanC) acyl carrier protein 69.4 97 BN6_54890 penP AAD13539.1(lanD) ketoreductase 77.0 99 BN6_54900 penQ AAD13540.1 (lanL) aromatase/cyclase 65.6 96 BN6_54910 penS WP_015037166(jadF)* oxidoreductase 69.7 72 BN6_54920 penT AAV52248(jadH)* oxidoreductase 66.8 99 BN6_54930 penU WP_046087459 acyl-CoA carboxylase subunit beta 81.8 99 BN6_54940 penV WP_062205789 acyl-CoA carboxylase subunit epsilon 38.3 73 BN6 54950 penR2 AAB36584(jadR1) Phosphate regulon transcriptional regulatory protein PhoB 57.9 95 * putat i vl e y involved i n r in g o pening of Ri n g B Figure 3. Schematic representation of the biosynthetic gene cluster of pentangumycin (a) and a BLAST analysis of the pen genes (b). 3.5. Detailed Analysis and Border Prediction of the Biosynthetic Gene Cluster of Pentangumycin The putative borders ofpentangumycin’s BGC (MIBiG accession number: BGC0002070) are defined by the genes BN6_54730 (further as penA) on the left edge and BN6_54950 (further as penR2) on the right (Figure 3a). To predict their functions, individual pen genes were analyzed using BLASTx. the right (Figure 3a). To predict their functions, individual pen genes were analyzed using BLASTx. Detailed results are illustrated in Figure 3b. The core of the pen cluster is formed by genes penM, penN and penO, coding for the minimal polyketide synthase. Together with two cyclase genes, penL and penF, the cyclase/aromatase penQ and the ketoreductase penP are forming a minimal set of genes required for the biosynthesis of the angucyclinone backbone. The polyketide core region is surrounded by seven genes putatively involved in oxidative or reductive reactions. The gene penA encodes a cytochrome P450 enzyme and is adjacent to penB, which encodes a ferredoxin protein. Both genes have previously been described as a functiona l unit [48-50]. The genes penE, penG, penI, penS and penT encode a monooxygenase (penE), three oxidoreductases (penG, penI and penS) and an aldo-/ketoreductase (penT ). Except for PenE, PenS and PenT, the function of these enzymes cannot be predicted from the BLAST analysis. penE,penS and penT encode orthologues of JadG, JadF and JanH, respectively, which are three enzymes putatively i nvolved in the oxidative opening of ring B of j adomycin [51]. Other genes such as penC (encoding putative aminotransferase) and penD (encoding putative methyltransferase) are also supposed to be involved in post-PKS modifications of 1. penK, penU and penV encode c-, B-and e-subunits of a carboxyltransferase, respectively. The function of these enzymes in the biosynthesis of 1 is unclear. Two regulatory genes-penR1, encoding a putative transcriptional regulator with the predicted helix-turn-helix motif at the N-terminus, and penR2, coding for a SARP protein—are present in the cluster [52]. PenR1 most probably controls the transcription of the outward-oriented penH (Figure 3), similar to the landomycins lanK, lanJ regulatory system [53-55]. The gene penH encodes a putative major facilitator superfamily transporter [56]. In turn, PenR2 shows a high degree of similarity compared to many well-studied SARP proteins from angucycline biosynthetic gene clusters, including JadR1 from S. venezuelae [52]. Due to the similar structures of 1 and landomycin and jadomycin [41,57], we aimed to identify homologues to genes involved in the biosynthe t ic pathways of both aforementioned compounds within the pen gene cluster. The PKS core genes penL, penN, penO, penP and penQ showed a significant similarity towards the landomycin biosynthetic genes l anL, lanA, lanB, l anC and lanD, respectively (Figure 3b). The oxygenase penE and the oxidoreductases penS and penT showed a remarkable homology to the oxidase jadG and the oxidoreductases jadF and jadH, respectively [58]. Additionally, a similarity between individual genes including penA, penB, penC, penD and penH including their organization and the putative biosynthetic gene cluster with unknown product from Micromonospora sp. LB39 was observed (Figure 3b and Figure S53). We aimed to confirm the predicted borders (penA and penR2) of pentangumycins biosynthetic gene cluster. Therefore, two BACs (1E5_DEL _LEFT and 1E5_DEL _RIGHT) were designed. BAC 1E5_DEL _LEFT was constructed through deletion of 44.1 kb covering the left flanking region of the pen cluster upstream from the gene penA (Figure 3a), including the entire terpene cluster. The construction of 1E5_DEL_RIGHT was carried out i n a similar way by deleting a 36.3-kb large region downstream to the penR2 gene, including the entire predicted PKS-NRPS gene cluster. The constructed BACs were conjugated into S . lividans AYA6, and the production of 1 was analyzed. The recombinant strains carrying the mutated BACs still produced 1. Deletions of the terminal genes penA and penR2 caused complete cessation of pentangumycin biosynthesis in the recombinant strains carrying corresponding BAC clones 1E5ApenA and 1E5ApenR2. The metabolic profiles of S. lividans transconjugants showed that both genes are required for production of 1. These findings confirm that al l genes necessary for the production of 1 are within the predicted borders of the pen gene cluster. 3.6. Biosynthesis of Pentangumycins Core Structure Based on the available bioinformatic information described above (Figure 3b) and our experimental results, a biosynthetic route towards 1 was proposed (Figure 4) [59,60]. PenM, PenN and PenO utilize acetyl-CoA and malonyl-CoA to form the initial 20 carbon polyketide chain, which is reduced at position C-9 by PenP and cyclized by PenL and PenQ. Like in the case of all other angucycline type aromatic polyketides, the biosynthesis of 1 seems to proceed though the common intermediate UWM6 [57]. UWM6 is oxidized by PenE, Pens [60] and PenT, which leads to the formation of I [58]. The cleavage of ring B undergoes a Baeyer-Villiger oxidation through the putative intermediates II and III. The latter one is proposed to have an aldehyde functional group at position C-5 and a carboxyl group at position C-6. Similar to dauD in daunomycin biosynthesis, the ester cyclase PenF could perform a Knoevenagel condensation, which results in the formation of the C-C bond between the B-position of 4-hydroxyl-phenylpyruvic acid, an intermediate of the tyrosine degradation catalyzed by an aromatic amino acid transaminase [61], and the aldehyde of III, resulting in intermediate IV. An intramolecular nucleophilic reac t ion between the y-hydroxyl and the c-keto groups cleaves oxalic acid and thus creates a double bond between the former B and y positions and a conjugated system. The subsequent reformation of ring B is catalyzed by the decarboxylation at position C-5 and directed through the six-membered transition state towards the former B-position of the 4-hydroxyl-phenylpyruvic acid, leading to intermediate V. The aromatase/cyclase PenQ is proposed to be responsible for the aromatization of ring B, resulting i n i ntermediate VI. To prove t hat t he f ifth ring of 1 derives from the incorporation of 4-hydroxyl-phenylpyruvic acid, a culture of S. lividans △YA6_1E5 was supplemented with fully labeled 13C-9_15N L-tyrosine. As a result, a mass shift to m/z 475 [M+H]+ was observed for the peak that corresponds to 1 (Figures S72-S77). This shift clearly shows that t he fifth ring of 1 derives from i ncorporation of L-tyrosine , as it i s i n the case of j adomyc i n, and not f rom acetate, as it was shown for formicamycin [44]. Furthermore, the lack of two tyrosine-derived carbons in the final structure shows t hat the insertion of the f ifth ring most probably occurs by the proposed mechanism with the loss of a-and B-carbons of tyrosine and the exchange of carbon six of IV with one of L-tyrosine’s in Microorganisms 2020, 8, x FOR PEER REVIEW 12 of 18 Figure 4. Biosynthetic pathway of pentangumycin; Roman numerals: proposed intermediates; Arabic Figure 4. numerals: intermediates with identified fitting masses after single gene knockouts and overexpression of the regulatory genes; Black pathway: methylation as first biosynthetic step after aromatization; Grey pathway: methylation as last biosynthetic step to complete the formation of pentangumycin (1).3.8. Analysis of Regulatory Gene Functions To elucidate the functions of the identified regulatory genes on the biosynthesis of 1, we 3.7. Tailoring Steps of Pentangumycin’s Biosynthesis To elucidate t he role of the remaining genes encoding the putative tailoring enzymes involved in post PKS modifications of 1, a set of single knockout mutant BACs (1E5△penA; 1E5ApenC; 1E5ApenD;1E5ApenE, 1E5ApenG; 1E5ApenI; 1E5ApenJ; 1E5ApenU; 1E5ApenV) was created and expressed in S. lividans AYA6. The produc t ion of 1 was not affected by deletions of penG, penI, penJ, penU and penV genes from BAC 1E5. At the same t i me, recombinant strains carrying 1E5 with deleted penA, penC, penD and penE genes lack the production of 1, but were found (except penE mutant) to accumulate related angucyclinone compounds based on optical absorption and mass spectrometry data (Figure 4) The inactivation of penA resulted in the production of a new compound with m/z 411.1218[M+H]+that is very close to the calculated exact mass of 410.1154 Da of proposed intermediate 3 (A2.190 ppm)(Figure S78). The deletion of penC led to accumulation of the proposed intermediate 4 (detected m/z 427.1174[M+H]+, calculated mass of 426.1103 Da, A0.409ppm) (Figure S79). The BAC lacking penD gene facilitated the production of another compound with m/z 454.1278[M+H]+ that corresponds to the mass of the proposed intermediate 5 (calculated exact mass of 453.1212,△1.506 ppm) (Figure S80). With the information gained from the analysis of all metabolic profiles from the recombinant strains carrying 1E5with deletion of the aforementioned genes, it can be hypothesized that PenA catalyzes a hydroxylation of the methyl group at position C-19 to form a primary alcohol. The alcohol is subsequently oxidized by one of the oxidoreductases to an aldehyde, resulting in predicted intermediate VII. A similar biosynthetic reaction was proposed for the borrelidin biosynthetic pathway [62]. PenC performs an amination of the aldehyde at this position, resulting in structure VIII. However, in the extract of S. lividans AYA6_1E5ApenC, only proposed intermediate 4 can be found, most probably because the aldehyde of VI I is not stable and is spontaneously reduced to an alcohol. Lastly, the deletion of penD leads to accumulation of 5, the last intermediate before formation of 1 lacking methylation of the hydroxyl group at position C-1 (Figure 4). The presence of t his intermediate, as well as methylated 3and 4, makes us believe that PenD can act on both substrates 5 and VI as well as that the presence of a methyl group at this posi t ion is not crucial for amination and acylation events. The deletion of penE leads to a complete loss of production of 1. Interestingly, the shunt product rabelomycin [63] can be ident i fied in all of our extracts obtained from knockout mutants (identi f ied with an external standard), but its production is drastically reduced in the extracts of S. lividans △YA6 1E5ApenE (Figure S81). This result indicates an involvement of penE in the early steps of the biosynthesis of 1. 3.8. Analysis of Regulatory Gene Functions To elucidate the functions of the identified regulatory genes on the biosynthesis of 1, we constructed recombinant BACs with deletions of penR1 and penR2. Both BACs were conjugated into S. lividans AYA6. The constructed recombinant strains S. lividans AYA6 1E5_ApenR1 and S. lividans AYA6 1E5_ApenR2were analyzed for production of 1 (Figure S82). While 1 was still detectable in the ApenR1 mutant with a 50-fold decrease in yield, i t was no longer present i n the S. lividans AYA6 1E5_ApenR2 extracts. From this observation and high similarity of penR2 to other SARP-encoding genes, it becomes obvious that this gene is encoding a pathway-specific regulator controlling the expression of structural pen genes. Furthermore, the overexpression of either penR1 or penR2 might lead to the increase in the production of 1. To test this hypothesis, both genes were cloned under the control of the moderately strong A3promoter [24] and introduced into S. lividans ▲YA6 using a marker-free vector system [64]. The resulting strains S. lividans AYA6_A3penR1 and S. l ividans AYA6_A3_penR2 were tested for the production of 1 and 2. Surprisingly, the production level of 1 remained unchanged i n both strains (Figure S83). S. lividans AYA6_A3_penR1 1E5 still produced 2, but 2 was absent in extracts of S. lividans AYA6_A3_penR21E5 (Figure S84). These results indicate an interdependence between the two regulatory elements, similar to the promoter effects observed in the production of bottromycin, and lead to the conclusion that the unbalanced expression of pentangumycin’s BGC is indeed the reason for accumulation of SEK90 and i ts derivatives by the heterologous strain carrying BAC 1E5 [65]. Furthermore, two novel peaks, 6 and 7, were detected in the extract of S. lividans AYA6_A3_penR1 1E5 (Figure S85). Based on LC-HRMS data, 6 could be proposed to be an analogue of 1 with tryptophan incorporated into ring B instead of tyrosine (detected m/z 491.15934 [M+H]+, calculated exact mass of 490.1529 Da,△1.646), while 7 seems to carry phenylalanine-derived moiety at the same position (detected m/z 452.14903[M+H]+, calculated exact mass 451.1419 Da, A0.485 ppm) (Figure S85). To prove t he possible nature of the identified derivatives 6 and 7, the culture of S. lividans AYA6_A3_penR11E5 was supplemented with either L-tryptophan (13C11H12N202) or L-phenylalanine (ring-Ds). As result, in the correspondingly fed cultures, the mass of 6 shifted to m/z 500.19 [M+H]+ (Figure S86), and the mass of 7 shifted to m/z 457.20 [M+H]+ (Figure S87), which corresponds to the incorporation of labelled L-tryptophan and labelled L-phenylalanine. 3.9. Origin of SEK90 From a structural perspective, 2 is a dimer of SEK43 l inked by a butyl group. It was shown by McDaniel that SEK43 is formed by spontaneous cyclization of the polyketide chain in mutants lacking the cyc l ase [47]. Since then, SEK43 has been reported as a shunt product in several studies of aromatic polyketide cyclases/aromatases [66-68]. At the same time, SEK43 was found i n the extract of natural strains only in the case of the aranciamycin-producing Streptomyces sp. Tu6384 [69]. On the other hand, its derivative, named SEK43F, resulting from fusion of SEK43 with the pyrole-like moiety, was isolated from the recombinant S. albus expressing fluostatin biosynthetic gene cluster [70]. Nevertheless, a dimerization of t his compound has not been described t hus far. SEK43 has the 4-hydroxypyrone structural motif. The reactiv i ty of 4-hydroxypyrones towards saturated aldehydes was reported previously [71]. It can be hypothesized that dimerization towards SEK90 can occur spontaneously inside the cell after the biosynthetic formation of SEK43 (Figure 5). To verify this theory, an experimental setup was designed; 6-Benzyl-4-hydroxy-2-pyrone that possesses the same structural motif as SEK43was incubated with butanal or formaldehyde i n distilled water, production medium i noculated with S. lividans AYA6 and production medium inoculated with S. lividans △YA6_1E5. If the dimerization of the pyrone with the corresponding aldehyde is performed by an enzymatic reaction, the expected dimers (1E5_CMP1, further designated as 8, and 1E5_CMP2, further designated as 9; see Figure S51) should only be present in the extracts of S. lividans AYA6 or S. l ividans △YA6_1E5, depending on the presence of the necessary enzymes in the heterologous host or the introduced BAC. If the reaction towards 2 is spontaneous, the dimerization of pyrone and aldehydes should occur in the distilled water mixture. After incubation, all reactions were extracted and analyzed with LC-HRMS. Peaks that correspond to dimerized pyrones were present in all mixtures, including the distilled water-based one, confirming our hypothesis about the spontaneous nature of the 4-hydroxypyrone dimerization. Furthermore, we have synthesized 8 and 9 to verify their structure by NMR (Tables S6 and S7; Figures S54-S70) [71]. The peaks observed in the distilled water mixtures had identical exact masses and HPLC retention times as the synthesized standards (Figures S71 and S72). Detailed analysis of the S. lividans △YA6 1E5extract revealed the presence of SEK43 (detected m/z 396.09679 [M+H]+, A0.242 ppm to SEK43 with the calculated exact mass of 368.0860 Da) (Figures S51 and S73) and a derivative of SEK90 connected by a methyl-instead of a butyl group named as SEK87 (detected m/z 749.1859 [M+H]+, A0.777ppm to SEK87 with the calculated exac t mass of 748.1792 Da) (Figures S51 and S74). None of the three compounds can be found in the extract of S. espanaensis. Thus, it is obvious that SEK90 as well as SEK87 derives from the combination of an unbalanced performance of pentangumycins minimal PKS, which seems to be uncoordinated with the activity of the Pen cyclases, leading to accumulation of SEK43 and the primary metabolism of S . lividans. A similar situation was proposed for the assembly of SEK43F, which resulted from interplay between the primary metabolism of the host strain S. albus and the heterologously expressed f luostatins type I I PKS gene cluster [70]. which resulted from interplay between the primary metabolism of the host strain S. albus and the heterologously expressed fluostatins type II PKS gene cluster [70]. Figure 5. Origin of SEK90. 4. Conclusions In conclusion, we have shown that the described “cluster to compound” approach based on a systematic analysis of the genome of an Actinobacteria strain for unique BGCs combined with heterologousexpressioncanbeasourceofnewnaturalproductswithuniquestructural characteristics. As result, two novel compounds were isolated in sufficient amounts for structure elucidationandtestingoftheirbiologicalactivity.Throughgeneticmanipulationandfeeding experiments,thebiosynthesisof1waselucidated.SEK90representsanewshuntproductin biosynthesis of aromatic polyketides that arises from interplay between expressed BGC and primary metabolism of the host strain. At the same time, pentangumycin, the actual final product of the pen gene cluster is a new member of the angucyclinone family with a heavily modified ring system and a new, unusual tailoring modification which was not observed in other members of this group of natural products before. Importantly,wehavedemonstratedthatStreptomycesspeciescanbeusedashostsfor heterologous expression of biosynthetic gene clusters derived from distantly related Actinobacteria species, even if their expression is silent in the wild-type strain. However, a still relatively low success rateshouldtriggerthedevelopmentofhoststrainsotherthanStreptomycesforactinobacterial secondary metabolism gene cluster expression. The Supporting Information is available and contains supplementary figures and tables including HPLC-MS SupplementaryMaterials:Thefollowingareavailableonlineatwww.mdpi.com/xxx/s1,TheSupporting Information is available and contains supplementary figures and tables including HPLC-MS chromatograms, NMR spectra and additional structures. AuthorContributions:Laboratorywork,N.G.;analytics,N.G.;structureanalysis,N.G.,C.P.andJ.Z.; manuscript preparation, N.G.; review and editing, Y.R. and A.L.; supervision, Y.R. and A.L.; funding acquisition A.L. All authors have read and agreed to the published version of the manuscript. Funding: This research received no external funding. Conflicts of Interest: The authors declare no conflict of interest. Refer W e a n l c sh e ,s C.T. A chemocentric view of the natural product inventory. Nat. Chem. Biol. 2015, 11, 620-624. 1. Berdy, J. Bioactive microbial metabolites. J. Antibiot. (Tokyo) 2005, 58, 1–26. 3. Kong, D.X.; Guo, M.Y.; Xiao, Z.H.; Chen, L.L.; Zhang, H.Y. Historical variation of structural novelty in a natural product library. Chem. Biodivers 2011, 8, 1968-1977. [CrossRef ] 4 Cech, N.B.; Yu, K. Mass spectrometry for natural products research: Challenges, pitfalls, and opportunities. LCGC N. Am. 2013,31,938-947. 5. Zhang, Y.; Pan, G.; Zou, Z .; Fan, K.; Yang, K.; Tan, H. Jadr-mediated feed-forward regulation of cofactor supply in jadomycin biosynthesis . Mol. Microbiol. 2013, 90,884-897. [CrossRef] [PubMed ] 6. National-Human-Genome-Research-Institute; Wetterstrand,K.A. DNA Sequencing Costs: Data from the Nhgri Genome Sequencing Program (gsp). 2019. Available online: www.ge nome.g o v /sequencing co stsdata (accessed on 16 April 2019). 7. Ahmed, N. A flood of microbial genomes-do we need more? PLoS ONE 2009, 4, e5831. [CrossRef ][PubMed ] 8. Doroghazi, J .R.; Albright, J.C.; Goering, A.W.; Ju, K.S.; Haines, R.R.; Tchalukov, K.A.; Labeda, D.P.; Kelleher, N.L .; Metcalf, W.W. A roadmap for natural product discovery based on large-scale genomics and metabolomics. Nat. Chem. Biol.2014,10,963-968. [CrossRef] [Pu bMed ] 9 Tokovenko, B.; Rebets, Y.; Luzhetskyy, A. Automating assessment of the undiscovered biosynthetic potential of actinobacteria. bioRxiv 2016,036087. [CrossRef ] 10. Weber, T.; Blin, K.;Duddela, S.; Krug, D.; Kim, H.U.; Bruccoleri, R.; Lee, S.Y.; Fischbach, M.A.;Muller, R.; Wohlleben, W.;et al . Ant i smash 3.0-a comprehensive resource for the genome mining of biosynthet i c gene clusters. Nucleic Acids Res. 2015,43,W237-W243. [CrossRef ] 11 Myronovskyi, M.; Rosenkranzer, B.; Nadmid, s.; Pujic, P.; Normand, P.; Luzhetskyy, A. Generation of a cluster-free streptomyces albus chassis strains for improved heterologous expression of secondary metabol i te clusters. Metab. Eng. 2018, 49, 316-324. [CrossRef ] 12. Gomez-Escribano, J.P.; Bibb, M.J . Engineering streptomyces coelicolor for heterologous expression of secondary metabolite gene clusters. Microb. Biotechnol. 2011, 4, 207-215. [CrossRef] [PubMed ] 13. Mazodier, P; Petter, R.; Thompson, C. Intergeneric conjugation between escherichia coli and streptomyces species. J. Bacteriol. 1989,171,3583-3585. [CrossRef] [Pu bMed 14. Gust, B.; Chandra, G.; Jakimowicz, D.; Yuqing, T; Bruton, C.J .; Chater, K.F. Lambda red-mediated genetic manipulation of antibiotic-producing streptomyces. Adv. Appl . Microbiol.2004,54,107-128.[PubMed ] 15. Blin, K.; Wolf, T.; Chevrette, M.G.; Lu, X.; Schwalen, C.J.; Kautsar, S.A.; Suarez Duran, H.G.; de Los Santos, E.L .C .; Kim, H.U.;Nave, M.; et al. Antismash 4.0-improvements in chemistry prediction and gene cluster boundary identification. Nucleic Acids Res. 2017, 45, W36-W41. [CrossRef ] [Pu bMed ] 16. Skinnider, M.A.; Merwin, N.J.; Johnston, C.W.; Magarvey, N.A. Prism 3: Expanded prediction of natural product chemica l structures from microbial genomes. Nucleic Acids Res. 2017,45, W49-W54. [CrossRef ][PubMed l 17. Labeda, D.P.; Lechevalier, M.P. Amendment of the genus saccharothrix labeda et al. 1984 and desc r iptions of saccharothrix espanaensis sp. Nov ., saccharothrix cryophilis sp. Nov., and saccharothrix mutabilis comb. Nov. Int . Syst . Evol . Microbiol. 1989,39,420-423. [CrossRef ] 18. Strobel, T.; Al-Dilaimi, A.; Blom, J.; Gessner, A.; Kalinowski,J.; Luzhetska, M.; Puhler, A.; Szczepanowski , R.; Bechthold, A.; Ruckert, C. Complete genome sequence of saccharothrix espanaensis d s m 44229(t) and comparison to the other completely sequenced pseudonocardiaceae. BMC Genom. 2012,13,465. [CrossRef ] 19. Singh, M.P.; Petersen, P.J .; Weiss, W.J.; Kong, F.;Greenstein, M. Saccharomicins, novel heptadecaglycoside antibiot i cs produced by saccharothrix espanaensis: Ant i bacterial and mechanistic activities. Antimicrob . Agents Chemother. 2000, 44,2154-2159. [CrossRef] 20. Green, M.R.; Sambrook, J .; Sambrook, J . Molecular Cloning: A Laboratory Manual; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2012. 21. Kieser, T.; Bibb,M.J .; Buttner, M.J.; Chater, K.F.; Hopwood, D.A. Practical Streptomyces Genetics; John Innes Founda t ion: Norwich, UK, 2010; p. 613. 22 Altschul , S.F.; Gish, W.; Miller, W.;Myers, E.W.; Lipman, D.J. Basic loca l alignment search tool. J . Mol. Biol .1990,215,403-410. [CrossRef ] 24. Siegl , T.; Tokovenko, B.; Myronovskyi , M.;Luzhetskyy, A. Design, construc t ion and characterisation of a synthetic promoter library for fine-tuned gene expression in ac t inomycetes. Metab. Eng. 2013,19,98-106.[CrossRef ][PubMed] 25. Brenner, M. Molekularbiologische Untersuchungen zur Biosynthese der Oligosaccharid-Antibiotika Saccharomicin a und b in Saccharothrix Espanaensis. Ph.D. Thesis, Albert-Ludwig-Universitat Freiburg, Freiburg im Breisgau , Germany, 2016. 26. Berner, M.; Krug, D.; Bihlmaier, C.; Vente, A.; Muller, R.; Bechthold, A. Genes and enzymes i nvolved in caf f eic acid biosynthesi s in the actinomycete saccharothrix espanaensis. J. Bacteriol. 2006, 188,2666-2673.[CrossRef ] [PubMed 27. Banerjee, P.; Erehman,J .; Gohlke, B.O.; Wilhelm, T.; Preissner, R.;Dunkel, M. Super natural ii一A database of natural products. Nucleic Acids Res. 2015, 43,D935-D939. [CrossRef] [PubMed ] 28. Klementz, D.; Doring, K.; Lucas, X.; Telukunta, K.K.; Erxleben, A.; Deubel, D.; Erber, A.; Santillana, I .; Thomas, O.S.; Bechthold, A.; et al. Streptomedb 2.0-An extended resource of natural products produced by streptomycetes. Nucleic Acids Res. 2016, 44, D509-D514. [CrossRef ] 29. Pait, I.G.U.; Kitani, S.; Roslan, FW.; Ulanova, D.; Arai, M.; Ikeda, H.; Nihira , T. Discovery of a new diol-containing polyketide by heterologous expression of a silent biosynthetic gene cluster from streptomyces lavendulae f ri-5. J . Ind. Microbiol. Biotechnol. 2018,45,77-87. [CrossRef ] [Pu bMed ] 32. Kane, J.F. Effects of rare codon clusters on high-level expression of heterologous proteins in escherichia coli. Curr. Opin. Chem. Biol . 1995, 6,494-500. [CrossRef ] 33. Liu, S.; Wang , M.; Du, G.; Chen, J . Improving the active expression of t ransglutaminase in streptomyces lividans by promoter engineering and codon optimization. BMC Biotechnol . 2016,16,75. [CrossRef ] 34. Olano, C.; Garcia, I.; Gonzalez, A.; Rodriguez, M.; Rozas, D.; Rubio, J .; Sanchez-Hidalgo, M.; Brana, A.F.; Mendez, C.; Salas, J.A. Activation and identification of five c l usters for secondary metabolites i n streptomyces albus j1074. Microb. Biotechnol. 2014,7,242-256. [CrossRef ] 35. Freudl, R.; Braun, G.; Hindennach, I; Henning, U. Lethal mutations in the structural gene of an outer membrane protein (ompa) of escherichia coli k12. Mol. Gen. Genet . 1985, 201,76-81. [CrossRef ] [PubMed] 36. Qiu, J.;Zhuo, Y.; Zhu, D.; Zhou, X.; Zhang, L.; Bai, L.; Deng, Z. Overexpression of the abc t ransporter avtab increases avermectin production in streptomyces avermit i lis. Appl. Microbiol. Biotechnol. 2011, 92,337-345.[CrossRef ][PubMed1] 37. Doull, J.L.; Singh, A.K.; Hoare, M.; Ayer, S.W. Conditions for t he produc t ion of jadomycin b by streptomyces venezuelae isp5230: Effects of heat shock, ethanol treatment and phage infection. J . Ind. Microbiol. Biotechnol.1994,13,120-125. 38. Pishchany, G.; Mevers, E.; Ndousse-Fetter, S.; Horvath, D.J., Jr.; Paludo, C.R.; Silva-Junior, E .A.; Koren, S.; Skaar, E.P.; Clardy, J .; Kolter, R. Amycomicin is a potent and specific antibiotic discovered with a targeted interaction screen. Proc . Natl. Acad. Sci. USA 2018, 115,10124-10129.[CrossRef ] [PubMed ] 39. Moon, K.; Xu, F; Zhang, C.; Seyedsayamdost , M.R. Bioactivity-hites unveils cryptic antibiotics encoded in actinomycete bacteria. ACS Chem. Biol. 2019, 14,767-774.[CrossRef] [PubMed ] 40. Seyedsayamdost, M.R. Toward a global picture of bacterial secondary metabolism. J . Ind. Microbiol. Biotechnol .2019,46,301-311. [CrossRef 41. Doul l, J .L.; Ayer, S.W.; Singh, A.K.; Thibault, P. Production of a novel polyketide antibiotic, j adomycin b, by streptomyces venezuelae following hea t shock. J. Antibiot. (Tokyo) 1993, 46,869-871.[CrossRef ] 42. Drautz, H.; Zahner, H.; Rohr,J.; Zeeck, A. Metabolic products of microorganisms. 234. Urdamycins, new angucycline antibiotics from streptomyces fradiae. I . Isolation, characterization and biological properties J . Antibiot. (Tokyo) 1986, 39, 1657-1669. [CrossRef] 43. Robertson, A.W.;MacLeod,J .M.;MacIntyre,L.W.;Forget, S.M.; Hall, S.R.; Bennett, L.G.; Correa, H.;Kerr,R.G.; Goralski, K.B.; Jakeman, D.L. Post polyketide synthase carbon-carbon bond formation in t ype-ii pks-derived natural products from streptomyces venezuelae. J. Org. Chem. 2018, 83, 1876-1890. [CrossRef ] 44. Qin,Z.;Munnoch,J.T; Devine,R.; Holmes,N.A.; Seipke,R.F .; Wi l kinson, K.A.; Wilkinson, B.; Hutchings,M.I . Formicamycins, antibacterial polyketides produced by streptomyces formicae isolated from african tetraponera plant-ants. Chem. Sci. 2017,8,3218-3227.[CrossRef ] 45. Potterat , O.;Puder, C.; Wagner, K.; Bolek, W.; Vettermann,R.; Kauschke, S.G.Chlorocyclinones a-d, chlorinated angucyclinones f rom streptomyces sp. Strongly antagonizing rosigli t azone-induced ppar-gamma activation. J. Nat . Prod. 2007,70,1934-1938. [CrossRef ] [Pu bMed] 46. Louden, B.C.; Haarmann, D.; Lynne, A.M. Use of blue agar cas assay for siderophore detection. J. Microbiol. Biol . Educ. 2011,12,51-53. [CrossRef ][PubMed ] 47 McDaniel , R.; Ebert-Khosla, S.; Hopwood, D.A.; Khosla, C. Rational design of aromatic polyketide natural products by recombinant assembly of enzymatic subunits. Nature 1995, 375,549-554. [CrossRef ][PubMed ] 48. Girvan, H.M.; Munro, A.W. Applications of microbial cytochrome p450 enzymes in biotechnology and synthetic biology. Curr. Opin. Chem. Biol. 2016, 31, 136-145. [CrossRef ] 49. Waskell, L.; Kim, J .-J .P. Elec t ron transfer partners of cytochrome p450. In Cytochrome p450: Structure, Mechanism, and Biochemistry; de Ortiz Montellano, P.R., Ed.; Springer: Cham, Switzer l and, 2015; pp. 33-68 511..(Chen, Y.H.; Wang, C.C.; Greenwel l , L.; Rix, U .; Hoffmeister, D.; Vining, L.C.; Rohr, J; Yang, K.Q. Functional analyses of oxygenases in jadomycin biosynthesis and identification of jadh as a bifunctional oxygenase/dehydrase. J . Biol. Chem. 2005, 280, 22508-22514. [CrossRef ] 52. Wang, L.; Tian, X.; Wang, J; Yang, H.; Fan, K.; Xu, G.; Yang, K.; Tan, H. Autoregulation of antibiotic biosynthesis by binding of t he end product to an atypical response regulator. Proc . Natl. Acad. Sci. USA 2009,106,8617-8622. [CrossRef ] 53. Novakova, R.; Homerova, D.; Feckova, L.; Kormanec, J. Characterization of a regulatory gene essential for the production of the angucycline-like polyketide antibiotic auricin in streptomyces aureofaciens ccm 3239. Microbiology 2005,151,2693-2706.[CrossRef ] 54 Ostash, I .; Ostash , B.; Walker, S.; Fedorenko, V. Proton-dependent transporter gene lndj confers resistance to landomycin e in streptomyces globisporus 1912. Russ. J. Genet. 2007,43,852-857. [CrossRef ] 55. Rebets, Y.; Dutko, L.; Ostash, B.; Luzhetskyy, A.; Kulachkovskyy, O.; Yamaguchi, T.; Nakamura, T.; Bechthold, A.; Fedorenko, V. Function of lani in regulation of landomycin a biosynthesis in streptomyces cyanogenus s136 and cross-complementation studies with streptomyces antibiotic regulatory proteins encoding genes. Arch. Microbiol. 2008,189,111-120.[CrossRef ] Pickens, L.B.; Tang, Y. Oxytetracycline biosynthesis. J . Biol . Chem. 2010,285,27509-27515. [CrossRef ]26. Westrich, L.;Domann, S.; Faust, B.; Bedford, D.;Hopwood, D.A.; Bechthold , A. Cloning and characterization of a gene cluster from streptomyces cyanogenus s136 probably involved in landomycin biosynthesis. FEMS Microbiol. Lett . 1999, 170,381-387. [CrossRef ] [P ubMed ] 58. Rix, U.; Zheng, J;Remsing Rix, L .L.; Greenwell , L .; Yang, K.;Rohr, J . The dynamic structure of jadomycin b and the amino acid incorporation step of its biosynthesis. J. Am. Chem. Soc. 2004, 126,4496-4497. [CrossRef ][PubMed l 59. Sharif, E.U.; O'Doherty, G.A. Biosynthesis and total synthesis studies on the jadomycin family of natural products. Eur. J. Org. Chem. 2012, 2012. [CrossRef] [PubMed ] 60. Yushchuk, O.; Kharel , M.;Ostash, I.; Ostash, B. Landomycin biosynthesis and i ts regulation i n streptomyces. Appl. Microbiol. Biotechnol. 2019, 103,1659-1665. [CrossRef ] 61. Dickens, M.L.; Ye,J.;Strohl, W.R. Analysis of c l ustered genes encoding both early and late steps in daunomycin biosynthesis by streptomyces sp. Strain c5. J. Bacteriol . 1995, 177,536-543. [CrossRef ] [PubMed ] 622..Olano, C.; Moss, S.J .; Brana, A.F.; Sheridan, R.M.; Math, V.; Weston, A.J .; Mendez, C.; Leadlay, P.F; Wilkinson, B.; Salas, J .A. Biosynthesis of t he angiogenesis inhibitor borrelidin by streptomyces parvulus tu4055: Insights into nitrile formation. Mol. Microbiol. 2004,52,1745-1756. [CrossRef ] 63. Liu, W.C.; Parker, L.; Slusarchyk, S.; Greenwood, G.L.; Grahm, S.F.; Meyers, E. Isolation, characterization, and structure of rabelomycin, a new antibiotic. J. Ant i biot. (Tokyo) 1970, 23,437-441. [CrossRef ] 65. Horbal , L.; Marques,F;Nadmid, S.; Mendes, M.V; Luzhetskyy, A. Secondary metabolites overproduction through transcriptional gene cluster refactoring. Metab. Eng. 2018, 49,299-315. [CrossRef ] 66. Meurer, G.; Gerlitz, M.; Wendt-Pienkowski, E.; Vining, L.C.; Rohr, J.; Hutchinson, C.R. I terative type i i polyketide synthases, cyclases and ketoreductases exhibit context-dependent behavior i n the biosynthesis of linear and angular decapolyketides. Chem. Biol. 1997,4, 433-443. [CrossRef ] 67. Zawada, R.J.; Khosla, C. Domain analysis of the molecular recognition features of aromatic polyketide synthase subunits. J . Biol . Chem. 1997,272, 16184-16188. [CrossRef ][PubMed ] 68. Bao, W.; Sheldon, P.J; Wendt-Pienkowski , E.; Hutchinson, C.R. The streptomyces peucetius dpsc gene determines the choice of starter unit in biosynthesis of the daunorubicin polyketide. J. Bacteriol. 1999,181,4690-4695. [CrossRef ][PubMed ] 69. Nachtigall, J. Strukturaufklarung von Naturstoffen aus Bakterien der Gattung Streptomyces. Ph.D. Thesis, Technischen Universitat Berlin, Fakultat II-Mathematik und Naturwissenschaften der Technischen, Berlin, Germany, 2012. 70. Huang, C.; Yang, C.; Zhu, Y.; Zhang, W.; Yuan, C.; Zhang, C. Marine bacterial aromatic polyketides from host-dependent heterologous expression and fungal mode of cyclization. Front. Chem. 2018, 6,528.[CrossRef ][PubMed ] 71.de March, P; Moreno-Manas, M.; Casado, J; Pleixats, R.; Roca, J .L.; Trius, A. The reactivity of 4-hydroxy-6-methyl-2-pyrone t t o o w wards 紡 aliphatic :saturatedand a,B-unsaturatedaldehydes. J. Heterocycl . Chem. 1984, 21,85-89. [CrossRe f Publisher's Note: MDPI stays neutral with regard to jurisdictional claims in published maps and i nstitutional affiliations. @ 2020 by the authors. Licensee MDPI , Basel, Switzerland. This article is an open access article distributed under the terms and condi t ions of the Creative Commons Attribution (CC BY) l icense (http://creativecommons.org/licenses/by/4.0/).