Comparative genomics and transcriptomics analyses provide insights into the high yield and regulatory mechanism of Norvancomycin biosynthesis in Amycolatopsis orientalis NCPC 2-48
Microbial Cell Factories volume 20, Article number: 28 (2021)
Norvancomycin has been widely used in clinic to treat against MRSA (Methicillin-resistant Staphylococcus aureus) and MRSE (Methicillin-resistant Staphylococcus epidermidis) infections in China. Amycolatopsis orientalis NCPC 2-48, a high yield strain derived from A. orientalis CPCC 200066, has been applied in industrial large-scale production of norvancomycin by North China Pharmaceutical Group. However, the potential high-yield and regulatory mechanism involved in norvancomycin biosynthetic pathway has not yet been addressed.
Here we sequenced and compared the genomes and transcriptomes of A. orientalis CPCC 200066 and NCPC 2-48. These two genomes are extremely similar with an identity of more than 99.9%, and no duplication and structural variation was found in the norvancomycin biosynthetic gene cluster. Comparative transcriptomic analysis indicated that biosynthetic genes of norvancomycin, as well as some primary metabolite pathways for the biosynthetic precursors of norvancomycin were generally upregulated. AoStrR1 and AoLuxR1, two cluster-situated regulatory genes in norvancomycin cluster, were 23.3-fold and 5.8-fold upregulated in the high yield strain at 48 h, respectively. Over-expression of AoStrR1 and AoLuxR1 in CPCC 200066 resulted in an increase of norvancomycin production, indicating their positive roles in norvancomycin biosynthesis. Furthermore, AoStrR1 can regulate the production of norvancomycin by directly interacting with at least 8 promoters of norvancomycin biosynthetic genes or operons.
Our results suggested that the high yield of NCPC 2-48 can be ascribed to increased expression level of norvancomycin biosynthetic genes in its cluster as well as the genes responsible for the supply of its precursors. The norvancomycin biosynthetic genes are presumably regulated by AoStrR1 and AoLuxR1, of them AoStrR1 is possibly the ultimate pathway-specific regulator for the norvancomycin production. These results are helpful for further clarification of the holistic and pathway-specific regulatory mechanism of norvancomycin biosynthesis in the industrial production strain.
Glycopeptide antibiotics, exhibiting outstanding activity against Gram-positive pathogens, are a class of widely known natural compounds produced by Actinomycetes with typical representatives of vancomycin, balhimycin, teicoplanin, A40926 etc. . Some of glycopeptide antibiotics have been approved for clinical use to treat persistent infections by Gram-positive multi-resistant pathogens since vancomycin was first approved in 1958 . Norvancomycin is an important glycopeptide antibiotic of vancomycin group, and has been used as a first-line empiric antibiotic therapy to prevent and treat the intracranial infections of MRSA (Methicillin-resistant Staphylococcus aureus) and MRSE (Methicillin-resistant Staphylococcus epidermidis) in China since 1969, although its structure was finally determined in the 1980s . Norvancomycin (Fig. 1a) has a closely similar chemical structure to vancomycin and shows a comparable antibacterial spectrum and activity to those of vancomycin. The fermentation potency of norvancomycin in the original strain Amycolatopsis orientalis CPCC 200066 was about 200 μg/ml (unpublished data). After series of physical and chemical mutagenesis, a high yield norvancomycin strain (A. orientalis NCPC 2-48) was obtained from A. orientalis CPCC 200066 by North China Pharmaceutical Group . Using the patented fermentation medium, the fermentation potency of high-yield strain could be increased up to more than 6000 μg/ml. Although A. orientalis NCPC 2–48 has been successfully applied to industrial large-scale production, the potential high-yield mechanism and its biosynthetic regulatory mechanism remain obscure.
The biosynthetic gene clusters (BGCs) for some known glycopeptide antibiotics were reported, including vancomycin (vcm) , balhimycin (bal) , A40926 (dbv) , teicoplanin (tcp) [8, 9], pekiskomycin (pek)  and kistamicin (kis)  etc. In our previous studies, we reported the genome sequence of A. orientalis CPCC 200066 and localized the biosynthetic gene cluster of norvancomycin (nvcm)  which displays a high resemblance with reported gene clusters of vcm in A. orientalis sub sp. orientalis KCTC 9412 T . More than forty genes were identified in the nvcm cluster (Fig. 1b), including four resistance proteins, three large NRPSs, several post-assembly tailoring enzymes, and a series of biosynthetic proteins for the supply of non-proteinaceous amino-acid precursors. Combining bioinformatics and biochemical studies, most biosynthetic steps of glycopeptide antibiotics are gradually deciphered [13, 14]. Since chemical structure of norvancomycin is closely similar to vancomycin (Fig. 1a), with only an N-methyl group absence at the N-terminus of the polypeptide, we speculate that the biosynthesis of norvancomycin also includes three steps by analogy with vancomycin [15, 16].
The regulatory mechanisms of glycopeptide antibiotics biosynthetic clusters are still not completely characterized. The VanR-VanS pair, a two-component regulatory system which controls the transcription of glycopeptide antibiotics resistance genes , have been found in the A40926 cluster (Dbv6 and Dbv22) , balhimycin cluster (VnlRS)  and teicoplanin cluster (Tei2 and Tei3, also named as Tcp7 and Tcp6) . Nevertheless, the vancomycin and norvancomycin cluster lack VanRS although components of self-resistance system, e.g., vanHAXY, are present in the vcm and nvcm cluster. StrR-family regulators are most commonly present in glycopeptide antibiotics biosynthetic clusters as specific pathway regulators. Bbr in balhimycin cluster binds to five promoter regions within bal cluster in vitro, and these genes are responsible for dehydrovancosamine biosynthesis (dvaA), oxidative cross-link (oxyA), transportation (tba), sodium-proton antiporter (orf7) and self-regulation (bbr) . Likewise, Dbv4 and Tcp28 (also named as Tei15*) were proved to be positive regulators in the A40926 cluster and teicoplanin cluster, respectively [21,22,23]. Moreover, the proposed palindromic sequence (GTCCAaN17TtGGAC) for Bbr-homologous regulator-binding site is also conserved within bal, dbv and tcp BGCs [20, 21, 23]. In addition, LuxR-family regulators have been characterized in A40926 and teicoplanin cluster. Dbv3 positively regulates production of A40926 and hierarchically control the transcription of dbv4 (StrR-family) . However, unlike StrR-family regulator, the direct target of LuxR-family regulator in the glycopeptide antibiotic clusters remain to be elucidated, since there is no direct binding evidence currently. The bioinformatic analysis indicates that there are four putative regulators (LuxR-family, StrR-family, TetR-family and AraC-family) present in or nearby the vcm or nvcm cluster. The regulatory mechanism of gene expression for the vancomycin or norvancomycin BGCs has not been reported yet.
The purpose of this study is to uncover the high-yield and biosynthetic regulatory mechanism of industrial strain NCPC 2-48 through exploring the genomic and transcriptomic features compared with the original strain CPCC 200066. Here, we present comprehensive comparative genomics and transcriptomics analyses of these two strains and show that the whole biosynthetic gene cluster of norvancomycin, as well as some primary metabolite pathways for the amino acid and other precursors of norvancomycin were comprehensively upregulated. The cluster situated regulators, AoStrR1 and AoLuxR1, play a positive role in norvancomycin production, and almost all of the norvancomycin biosynthetic genes are directly controlled by AoStrR1. These results have provided insights into the mechanism of high yield production in the industrial strain, as well as the role of cluster-specific regulators in the norvancomycin biosynthesis.
Results and discussion
Comparative genomics analysis of norvancomycin producing strains
The complete genome of industrial producing strain A. orientalis NCPC 2-48 and that of original strain A. orientalis CPCC 200066  are both circular chromosomes of 9.5 Mb with a G+C content of 68.84% (Fig. 1c). The genome analysis of industrial strain showed that it contained 8,705 genes, and the total length of genes was 8,579,274 bp, which makes up to 90.41% of genome. There are 404 tandem repeat sequence (35,253 bp), which makes up to 0.3715% of genome, 325 minisatellite DNAs, 11 microsatellite DNAs, 50 tRNAs, and 12 rRNAs.
We compared the genomes of NCPC 2-48 to CPCC 200066 and found that two genomic sequences are extremely similar with a high identity of 99.97%. There was no large fragment duplication or deletion in the entire genome of NCPC 2-48 compared with CPCC 200066. Moreover, the internal structure of the chromosome and gene order were largely conserved without rearrangements appeared in the genome of NCPC 2-48. Both strains have extremely similar codon usage of 50 tRNA genes. For the norvancomycin biosynthesis, further comparative analysis of secondary metabolism gene clusters showed that no duplication, SVs (structural variations) or InDels (insertions and deletions of small fragments (≤ 50 bp)) were occurred in nvcm biosynthetic gene cluster in the industrial strain. It suggested that the high-yield of the industrial strain is not caused by the increased copy number of the nvcm biosynthetic gene cluster or gene mutations (SVs or InDels) within the cluster.
The genomic difference between these two strains is mainly accounted for three SVs, including two deletions of 965 bp (SV1) and 12,212 bp (SV2) fragments and one insertion of 12,076 bp (SV3) fragment in strain NCPC 2-48 (Fig. 1c, Additional file 1: Table S3). All these SVs appear at coding DNA sequence (CDS) regions but far away from the nvcm biosynthetic gene cluster. There are 34 protein coding genes involved in these fragments of deletion and insertion, including 3 regulators, 1 transporter, 3 transposases, 9 other enzymes and 18 unknown proteins (Additional file 1: Table S3). The coding gene information of the three SVs is shown in Additional file 1: Table S4.
The first structural variation, SV1, occurred in gene B37_4355 (Additional file 1: Table S4). The 965 bp of B37_4355 sequence is deleted in genome of NCPC 2-48. SV2 is a deleted fragment of 12,212 bp in genome of NCPC 2-48, corresponding to the original strain’s chromosome position from gene B37_6566 to B37_6583. There are 17 possible CDSs in this region (Additional file 1: Table S4). The further analysis of the sequence flanking SV2 showed that one lacI family transcriptional regulatory gene is located upstream. Usually, lacI family regulators control the expression of some key enzymes involved in carbon metabolism, and regulate the transcription of a series of downstream genes including some transcriptional factors . The deletion of genes downstream lacI might block the regulation of lacI and enhance glucose catabolism of primary metabolism, which is conducive to the precursor sugar synthesis of norvancomycin. SV3 is a 12,076 bp of insertion, which contains 16 possible CDSs (Additional file 1: Table S4). Analysis of function of flanking genes revealed that there are some key enzymes of primary metabolism such as acetyl-CoA dehydrogenase and acetyl-CoA synthetase at downstream of SV3.
In addition to three fragments of deletion and insertion, there are also 216 InDels present in the genome of high-yield strain. The coding regions of some important enzymes related to primary and secondary metabolism together with transcription factors, such as LacI, LysR, TetR, MerR, YebC/PmpR, SARP family proteins and two-component regulators were present in these InDels (Additional file 1: Tables. S5, S6). The InDels in these transcriptional factors may affect the expression of norvancomycin through some unknown regulatory mechanism. However, the function of these transcriptional factors remains unexplored. In these SVs and InDels, there are no functionally known genes directly related to the biosynthetic pathway of norvancomycin and its precursors based on their functional annotation. Thus, the genomic variations of industrial strain would be difficult to give a simple explanation on its high-yield production of norvancomycin. Since some key enzymes of primary metabolism were found in these SVs and InDels, or in flanking segments, we speculated that these genomic mutations may change the metabolic flow by affecting expression of the important enzymes of primary metabolism. Meanwhile, more than ten regulatory genes were detected in SVs and InDels. Along with the loss, insertion or mutation of these regulatory genes, the holistic regulation of primary or secondary metabolism and cell growth may be changed in the high-yield producer strain, which ultimately makes the overall metabolic flow more favorable for the biosynthesis of norvancomycin.
Transcriptomics profiling at norvancomycin producing strains
In order to further investigate the mechanism of high yield production of norvancomycin, the transcriptomic analysis of norvancomycin original strain and industrial strain were carried out at three different time points (12 h, 24 h and 48 h). RNA from NCPC 2-48 and CPCC 200066 were extracted and sequenced, and an average of 23,578,714 raw reads were generated. After removing low-quality reads, the average number of remaining clean reads was 23,523,435, and the average comparison rates of clean reads to the reference gene and reference genome were 80.54% and 96.69%, respectively. The statistics on the sequencing data for each sample is shown in Additional file 1: Table S7. Pair-wise differentially expressed gene (DEG) analyses revealed that more than 2000 DEGs, about one-fourth of the total genes in the genome, had significantly lower or higher transcript abundance (fold change (FC) > 2 and false discovery rate (FDR) ≤ 0.001) at each time points (12 h, 24 h and 48 h) in NCPC 2–48 relative to CPCC 200066, as shown in Fig. 2a. The transcriptional levels of some genes within the nvcm cluster were verified by RT-qPCR (Additional file 1: Fig. S1). The results of DEGs hierarchical clustering analysis showed that differential gene expression pattern was similar at 24 h and 48 h (Fig. 2b), more than two-thirds of the differential genes were up-regulated in the industrial strain. Interestingly, more genes were downregulated at 12 h compared to other two time points (Fig. 2b), and there are 105 genes downregulated at 12 h but upregulated at 24 h and 48 h time points. KEGG pathway search showed that 35 out of these 105 genes are located in nvcm cluster. The biosynthetic pathways of secondary metabolites are usually activated in a growth phase-dependent manner, so that the genes responsible for secondary metabolism coincide with the onset of stationary phase in liquid fermentation in microorganisms. In the case of high-yield strain, transcriptional levels of genes responsible for biosynthesis of norvancomycin are lower in the early stage of growth (12 h), and then upregulated abruptly from 24 to 48 h, showing that the norvancomycin biosynthesis is more strictly controlled during the different growth stages. Due to the similarity of gene expression patterns between 24 h and 48 h, and most of DEGs included at 24 h (2,039 up-regulated and 628 down-regulated), we then analyzed the functional pathway enrichment of DEGs at 24 h based on KEGG database. Enrichment analysis of the functional categories of the transcriptome indicated that a total of 1,764 differential genes were annotated into 150 metabolic pathways, most of them related to the primary metabolism such as nitrogen metabolism and arginine biosynthesis, as well as the biosynthesis of secondary metabolites such as norvancomycin, tetracyclines and other type II polyketides, degradation of naphthalene and aromatic compounds, tyrosine and inositol phosphate metabolism processes. Top 20 of most specific KEGG enrichment results as shown in Fig. 2c.
In particular, the visualization of the transcriptome (Fig. 2d) showed that in the nvcm cluster, all key enzymes related to the biosynthesis of respective unit substrates such as Bht (VcmD, OxyD, Vhp), Hpg (Pdh, HmaS, HmO, HpgT), Dpg (DpgA/B/C/D, HpgT) and vancosamine (VasA/B/C/D/E), heptapeptide assemblage (VcmA/B/C) and post-modifications (OxyA/B/C, Vhal, GtfD/E) were transcribed significantly higher in high yield strain than that of original strain at 24 h and 48 h (Fig. 3, Additional file 1: Table S8). It suggested that the increased transcriptional level of norvancomycin's biosynthetic genes (1.9–8.6-fold upregulated at 24 h, and 3.0–18.3-fold upregulated at 48 h) directly promote high yield of norvancomycin in NCPC 2-48. In addition, primary pathways for amino acid (Leu, Asn, Tyr) and glucose production were also upregulated (Fig. 3, Additional file 1: Table S9), including genes B37_4517 and B37_7337 for prephenate (the precursor of Hpg); B37_6997, B37_3479 and B37_6785 for Tyr (the precursor of Bht); B37_7779 for Malonyl-CoA (Dpg's precursor); B37_4701, B37_8154 and B37_2225 for Leu; B37_7110 and B37_3479 for Asn and its precursor Asp synthesis; B37_7117 (RfbA) for TDP-D-Glucose (vancosamine precursor). The results suggested that every step of the whole pathway of the biosynthesis of norvancomycin was significantly upregulated in the industrial production strain, from the abundant supply of amino acid and glucose precursors to the NRPS assembling and the post modification of the glycopeptide antibiotic.
AoStrR1 and AoLuxR1 positively regulate the biosynthesis of norvancomycin
The expression of antibiotic biosynthetic genes is usually regulated by cluster-specific regulators within the gene cluster. There are four putative regulatory genes (AoLuxR1, AoStrR1, AoTetR1, AoAraC1) located within or adjacent to the nvcm cluster in the norvancomycin producing strains. The homologues of the four genes also present in a vancomycin producing strain A. orientalis KCTC 9412T (Additional file 1: Fig. S2). But in another vancomycin producing strain A. keratiniphila HCCB 10007, only homologues of AoLuxR1, AoStrR1 and AoTetR1 are present (Additional file 1: Fig. S2). Differential expression analysis revealed that the regulatory genes AoLuxR1 and AoStrR1 were significantly upregulated in the high-yield strain, with the same trend as the structural genes of nvcm clusters (Fig. 4a, Additional file 1: Table S8). The transcriptional levels of these two regulatory genes increased since 24 h and were 23.3-fold and 5.8-fold at 48 h, respectively (Fig. 4a, Additional file 1: Table S8). AoTetR1 and AoAraC1 located near the nvcm cluster showed 2 ~ 3-fold higher transcription level at 12 h, 24 h and 48 h in the industrial strain than that in the original strain, but in a different trend to the structural genes of nvcm clusters, and the transcription level (Fragments Per Kilobase of exon model per Million mapped fragments, FPKM value) were much lower than that of other genes in nvcm cluster (Fig. 4a, Additional file 1: Table S8).
The AoStrR1 gene is located in the nvcm cluster (7,787,412–7,788,377 nt), with a total length of 966 bp, encoding the AoStrR1 protein of 321 amino acids. Aligning the amino acid sequences of 22 BGCs of known glycopeptide antibiotics from GenBank (Additional file 1: Table S10) showed that the homologues of AoStrR1 appeared in almost all searched clusters (Additional file 1: Fig. S3A, Table S11), with only exception of corbomycin, the newly discovered compound . These amino acid sequences were used to build phylogenetic tree by MEGA-X program  based on Neighbor-Joining method. Compared with the clades of Tcp28 (from teicoplanin BGC), AoStrR1 tended to group with the regulators encoded in BGCs for vancomycin, balhimycin, chloroeremomycin and A40926 (Additional file 1: Fig. S3A). Among them, Bbr and Dbv4, which were confirmed as positive regulators in balhimycin  and A40926  cluster respectively, are highly homologous to the AoStrR1 protein with an identity of 84% and 80%, respectively (Additional file 1: Table S11). Thus, we speculated that AoStrR1 might be a pathway-specific regulator of norvancomycin biosynthetic genes, and the increased transcriptional level of AoStrR1 may be an important factor to trigger production of norvancomycin.
The AoLuxR1 gene is close to the right border of nvcm cluster (7,724,319–7,724,996 nt), with a total length of 678 bp, encoding the AoLuxR1 protein of 225 amino acids. There are 14 possible LuxR-like regulators identified in or adjacent to the 22 glycopeptide BGCs (Additional file 1: Fig. S3B, Table S12). The LuxR phylogenetic tree revealed two main clades (Additional file 1: Fig. S3B). AoLuxR1 appeared to be related to the regulators encoded in BGCs for vancomycin, decaplanin, keratinimicin and nogabecin, with an identity of 83% ~ 98% (Additional file 1: Fig. S3B, Table S12). The other clade included Dbv3 (from A40926 BGC) and Tcp29 (from teicoplanin BGC), which have larger size (more than 500 amino acids) and share low consistency with AoLuxR1 (Additional file 1: Fig. S3B, Table S12). These results suggested that the function or regulatory target of AoLuxR1 may be different from well-characterized Dbv3 or Tcp29. Although the AoLuxR1 homologous gene was not found in the well-known balhimycin BGC in A. balhimycina DSM 5908 (there is no adjacent ORF sequences available in GenBank), it is highly conserved in the reported two vancomycin producing strains (Additional file 1: Fig. S2, Table S12).
In order to determine the regulatory function of AoStrR1 and AoLuxR1, we constructed AoStrR1 and AoLuxR1 over-expression plasmids based on pULVK2A vector , under its native promoter or ermE*p, a strong constitutive promoter, respectively, and then conjugated into A. orientalis CPCC 200066. Norvancomycin yield in the fermentation broths of the AoStrR1/AoLuxR1 over-expressing strains were detected by HPLC and LC–MS. Fermentation results showed both AoStrR1 and AoLuxR1 genes could increase norvancomycin production, especially under the ermE*p promoter, which led to a 2 ~ 5 times higher norvancomycin yield (Fig. 4b, d). To confirm the role of AoStrR1 and AoLuxR1 in transcriptional regulation of nvcm cluster, the gene expression analysis was conducted by RT-qPCR analysis in over-expression strains. As expected, transcripts of the seven biosynthetic enzyme genes for norvancomycin biosynthesis, vcmA, oxyA, oxyB, vph, hmaS, vasA and ald, significantly increased in both over-expression strains (Fig. 4c, e). These results indicated that AoStrR1 and AoLuxR1 both acted as activators for norvancomycin biosynthesis.
AoStrR1 binds to eight promoter regions in the nvcm cluster
Cluster-specific regulatory proteins generally activate the transcription by binding to the promoter regions of one or more structural genes within a biosynthetic gene cluster, thereby promoting the biosynthesis of secondary metabolites and ultimately increasing their yield. In order to determine the potential target genes of AoStrR1 and AoLuxR1, we tried to express and purify his-tagged AoStrR1 and AoLuxR1 in E. coli BL21(DE3), and then perform electrophoretic mobility shift analysis (EMSA). Unfortunately, although different expression conditions were conducted, no soluble His-tagged AoLuxR1 was detected in the supernatant of recombinant E. coli BL21(DE3). The soluble AoStrR1 was expressed and purified as a fusion protein with the N-terminal His10-tag in E. coli BL21(DE3) (Fig. 5a). The expressed protein was checked by western blot with anti-His-tag antibody. As shown in Additional file 1: Fig. S4, a smaller His-tagged protein band was found in addition to the intact His-tagged AoStrR1 and it is likely that a fraction of AoStrR1 is degraded at the C-terminus of the protein. Since the DNA binding domain was predicted present in the C-terminus of AoStrR1, the degraded protein is supposed not be able to bind promoter regions and thus not affect the formation of binding bands between the intact protein and the DNA fragments in subsequent EMSA experiments.
To identify AoStrR1 target genes, 20 intergenic regions in the norvancomycin cluster responsible for the transcription of almost all structural genes and possible regulator genes were amplified and labeled with biotin as probes (200–500 bp) for EMSA experiments. The results showed that the recombinant His10-AoStrR1 could form a stable complex with the promoter regions upstream of vanY, AoStrR1, oxyA, oxyB, vhp, hmaS, vasA and ald (Fig. 5b, d). Addition of 100-fold unlabeled specific competitive DNA attenuated the shift band, indicating that the bindings of His10-AoStrR1 to the above probes are specific (Fig. 5b, d). The sequences of above eight AoStrR1 binding sites were input into the GLAM2  and MEME  software and the consensus binding sequence of AoStrR1 was identified as GTCCAN18TTGGAC containing an incomplete palindromic sequence (Fig. 5c). The AoStrR1 binding motif sequence is highly similar to the reported consensus sequence of its homologues Bbr, Dbv4 and Tcp28 (Tei15*) [20, 21, 23]. This indicated that the StrR-like regulators may have a conservative regulatory mechanism in different glycopeptide BGCs.
Based on RT-qPCR and the EMSA results, AoStrR1 might be the ultimate overall positive regulator responsible for norvancomycin production by binding directly to 8 promoters in nvcm cluster, including those of the genes and operons responsible for Bht biosynthesis (Vhp), Hpg biosynthesis (HmaS), vancosamine biosynthesis (VasA), linear heptapeptide cyclization (OxyAB) and self-resistance of norvancomycin (VanY). In addition, AoStrR1 did not bind to the promoters of the other 3 regulators (AoLuxR1, AoTetR1, AoAraC1, Additional file 1: Fig. S5) described here but its own promoter. When comparing the regulatory targets of those reported StrR-like regulators in different glycopeptide BGCs, it is interesting that except for Tcp28 (Tei15*), other three regulators (AoStrR1, Bbr and Dbv4) control one common biosynthetic step, heptapeptide cyclization (oxyA), and both AoStrR1 and Bbr could bind to promoter region of genes responsible for specialized amino sugar biosynthesis (vasA / dvaA) in their respective gene clusters. Different from Dbv4 and Tcp28 (Tei15*), Bbr and AoStrR1 are able to bind to its own upstream region. It suggested that the positive feedback mechanism of AoStrR1 might be responsible for the significant upregulation of norvancomycin production in the industrial strain. Moreover, among the four StrR-family regulators discussed above, only AoStrR1 could bind to the upstream of vanY, part of the putative self-resistance genes. Given that self-resistance regulatory system VanRS is lack in nvcm cluster, it is tempting to speculate that AoStrR1 is somehow involved in the regulation of self-resistance.
In this study, we compared the genomes of original norvancomycin producing strain CPCC 200066 and industrial strain A. orientalis NCPC 2–48. Three SVs and some InDels were found in the genome of high-yield strain compared with the original strain, and it is difficult to give a clear clue on its high-yield mechanism of norvancomycin. We further compared the transcriptomes of CPCC 200066 and NCPC 2-48 in three time points to identify differential expression genes involved in the production of norvancomycin. The results showed that the transcriptional upregulation of most nvcm biosynthetic genes and genes for precursor supply led to the high yield of norvancomycin in the industrial strain. Furthermore, two positive pathway-specific regulators of norvancomycin production were confirmed by overexpression of AoStrR1 and AoLuxR1. AoStrR1 positively regulates most nvcm genes by directly binding to multiple promoter regions in nvcm cluster. Thus, this study provides insights into the high yield mechanism and regulatory mechanism of norvancomycin in A. orientalis NCPC 2–48, which sets a foundation for future strain improvement.
Materials and methods
Strains, plasmids and growth conditions
The norvancomycin-producing strains A. orientalis CPCC 200066 and NCPC 2-48 were grown at 28 ℃ on solid Bennet medium (peptone 0.2%, glycerol 1%, beef extract 0.1%, glucose 1%, malt extract 0.3%, yeast extract 0.1%, agar 2%, pH 7.2.) for sporulation and liquid Bennet medium for fermentation. The Mannitol soya flour (MS) medium agar  was used for conjugation between E. coli and A. orientalis. Tryptic soy broth liquid medium  was used to grow strains for isolation of genomic DNA. E. coli ET12567/pUZ8002 was used for conjugal transfer according to the established protocol . All the strains and plasmids used in this study are listed in Additional file 1: Tables S1 and S2.
Genomic DNA extraction and sequencing
Genomic DNA were extracted from A. orientalis CPCC 200066 and NCPC 2-48 strains using the DNA extraction kit (TANBead, China) according to the manufacturer’s instructions. Genomic DNA library preparation, Illumina sequencing, chromosome assemble and annotation were carried out at Beijing Genomics Institute (Shenzhen, China), and performed as described . The second-generation sequencing platform Illumina Hiseq 2000 and a third-generation sequencing platform Pacbio RSII were used for high-throughput sequencing.
RNA extraction and sequencing
RNA extraction. Strains were collected at 12 h, 24 h and 48 h at the beginning of the fermentation, and total RNA were extracted using the TRIzol reagent (Invitrogen, USA) and chloroform followed by a PureLink™ RNA Mini Kit (Invitrogen, USA) according to the kit’s instructions. The RNA samples from the three different time points of the original strain were named B12, B24, and B48, respectively. The RNA samples of high-yield strains were named N12, N24 and N48, respectively. One biological replicate of each sample was used to perform RNA-seq analysis.
cDNA library construction and RNA-sequencing. The above 6 RNA samples were sent to Beijing Genomics Institute (Shenzhen, China) to build cDNA library. The second-generation sequencing platform BGISEQ-500 was used for high-throughput sequencing. After removing low-quality reads, the filtered data was compared with the reference sequence using HISAT  and Bowtie2  tools.
Differentially expressed genes (DEGs). The RSEM tool  was used to quantify gene expression. To eliminate the influence of gene length and sequencing quantity, the results of gene expression quantification were mapped to FPKM (fragments per kilobase of exon per million fragments mapped). We calculated the differential expression of the gene between different samples based on the FPKM value. In our analysis, differentially expressed genes (DEGs) were defined as genes with a fold change more than two times and FDR ≤ 0.001. Hierarchical clustering analysis on DEGs were performed by Cluster [34, 35] and show on javaTreeView . The functional enrichment of DEGs was analyzed by KEGG . The software Integrated Genome Browser  was used to create and view the whole genome expression on transcriptome level for each sample.
Reverse transcription quantitative PCR (RT-qPCR) analysis
Total RNAs were isolated as described above. 2 μg of each of the total RNA was used as a template for reverse transcription (RT), which was performed with the TransScript® One-Step gDNA Removal and cDNA Synthesis SuperMix (Transgen), using random primers following the manufacturer’s instructions. qPCR reaction was detected using CFX96 Touch Real-Time PCR Detection System (Bio-Rad). Each reaction (20 μl) contained 2.5 μl cDNA, 12.5 μl Fast Start Universal SYBR Green Master ROX (Roche) and 0.25 μM of forward and reverse primers (Additional file 1: Table S2). The relative cDNA level of target genes was normalized to the level of hrdB according to Pfaffl’s method .
AoStrR1 and AoLuxR1 over-expression in A. orientalis CPCC 200066
The vector pULVK2A  and pULVK2A-ermE*p which containing an ermE*p promoter, were used for gene overexpression. The 1310 bp fragment containing the coding region of AoStrR1 and its native promoter or 1005 bp fragment containing the coding region of AoLuxR1 and its native promoter were amplified from A. orientalis CPCC 200066 genomic DNA by PCR with specific primers (Additional file 1: Table S2), and were cloned into HindIII-BamHI sites of pULVK2A to obtain the p2A-AoStrR1 and p2A-AoLuxR1. The 981 bp fragment containing coding region of AoStrR1 or 678 bp fragment containing coding region of AoLuxR1 were amplified from A. orientalis CPCC 200066 genomic DNA by PCR with specific primers (Additional file 1: Table S2) and were cloned into NdeI-BamHI sites of pULVK2A-ermE*p to obtain the p2A-ermE*p-AoStrR1 and p2A-ermE*p-AoLuxR1. The constructed vector were introduced into A. orientalis CPCC 200066 by intergeneric conjugation from E. coli ET12567/pUZ8002 according to the established protocol .
Strain fermentation and norvancomycin yield detection
Amycolatopsis orientalis CPCC 200066 and its derivatives were cultured in 100 ml Bennet liquid medium at 28 ℃ for 2 days, and then 10% seed culture was transferred into 100 ml fresh Bennet for continuous fermentation at 28 ℃ for 4 days. Each fermented supernatant was collected by centrifugation at 12,000 rpm for 10 min, and then filtered with microporous membrane of 0.22 μm and analyzed by HPLC–MS (Agilent 1290-Agilent 1956 single quadrupole MS coupled system). The HPLC conditions were as follows: Agilent Eclipase plus C18 column (250 mm × 4.6 mm, 5 μm), mobile phase Solvent A was 100% MeOH, Solvent B was water with 0.1% formic acid. The HPLC program included column elution with a linear gradient of 5 to 30% solvent A over 30 min at 25 ℃. The flow rate was set at 0.8 ml/min. MS spectra data were collected in the positive-ion mode in which a mass range of m/z 150 to 2,000 covered. The NVCM peak was extracted ion chromatogram (EIC) of m/z 717.9 [M + 2H]2+.
AoStrR1 and AoLuxR1 protein expression and purification
The AoStrR1 and AoLuxR1 coding sequence were amplified by PCR from A. orientalis CPCC 200066 genomic DNA using specific primers (Additional file 1: Table S2). The amplified DNA fragment were cloned into NdeI-BamHI sites of pET-16b (Novagen, Madison, USA) to give pET16b-AoStrR1 or pET16b-AoLuxR1, and then expressed as fusion proteins with the N-terminal His10-tag in E. coli BL21(DE3). The transformed strains were grown in LB medium until an OD600 of 0.8 and induced by adding 0.1 mM IPTG for a further 7 h at 28 ℃. The bacteria were harvested by centrifugation (5000 rpm, 10 min, 4 ℃), and resuspended in 25 ml binding buffer (20 mM NaH2PO4, 500 mM NaCl, and 20 mM imidazole, pH 7.4), subsequently, lysed by high pressure continuous flow cell cracker (Constant systems, TS 0.75 kW). Cellular debris was removed by centrifugation (12,000 rpm, 20 min, 4 ℃). A HisTrap™ FF crude kit was used to purify recombinant protein, as described by the manufacturer (GE Healthcare). Fractions eluted from the column with 500 mM imidazole were passed through the PD-10 Desalting Columns (GE Healthcare) and eluted with 1 × TGEK buffer (250 mM Tris, 50% Glycerol, 5 mM EDTA, 500 mM KCl, pH 8.0) at 4 ℃, and then stored at − 80 ℃. The concentration of His10-tagged AoStrR1 was determined using BCA™ Protein Assay Kit (Pierce Biotechnology, Rockfold, USA) and its purity assessed using SDS-PAGE analysis.
Electrophoretic mobility shift analysis (EMSA)
DNA fragments containing the promoter regions of the nvcm genes were obtained by PCR using primers labeled at their 5′-ends with Biotin (Additional file 1: Table S2) and used as probes in EMSAs. Each 20 μl binding reaction consisted of 2 μl 10× binding buffer (100 mM Tris–HCl, 500 mM KCl, 10 mM DTT, pH 7.5), 2 μg salmon sperm DNA, 20 fmol labeled probe and 1 μM purified His10-tagged protein. The specific competitive reactions were carried out by adding 2 pmol competitors (unlabeled probes) apart from ingredients above. Reaction mixtures were incubated at room temperature for 20 min and then analyzed using a native 5% TBE-acrylamide gel and run at 4 °C, 100 V for 80 min. After electrophoretic transfer (15 V for 40 min) to nylon membrane (Hybond-H+, GE Amersham) by Trans-BLOT SD Semi-dry transfer cell (Bio-Rad), the probes were visualized by Lightshift Chemiluminescent EMSA Kit (Pierce Biotechnology), according to the manufacturer’s instructions.
Amino acid sequences of proteins homologous to AoStrR1 and AoLuxR1 were retrieved using Blastp tools and were further manually curated to ensure its location within or near the glycopeptide biosynthetic gene cluster. The phylogenetic trees were calculated using the MEGA-X program based on Neighbor-Joining method . The consensus motif of AoStrR1 binding sequence was represented as logos, obtained at the WebLogo website, using GLAM2  and MEME  algorithm.
The RNA-sequencing data are publicly available at NCBI’s GenBank’s repository under NCBI BioProject ID: PRJNA624813 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA624813), BioSample ID: SAMN03273949 & SAMN14585755, SRA ID: SRR11529125 ~ 127 & SRR11529348 ~ 350.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files.
Kahne D, Leimkuhler C, Lu W, Walsh C. Glycopeptide and lipoglycopeptide antibiotics. Chem Rev. 2005;105:425–48.
Van Bambeke F. Lipoglycopeptide antibacterial agents in gram-positive infections: a comparative review. Drugs. 2015;75:2073–95.
Liu YF, Zhou Y, Ling DK. Separation and identification of N-demethylvancomycin. Yao Xue Xue Bao. 1985;20:931–2.
Zhao Y, Dai M, Zhang L, Zhang XX. The establishment of genetic manipulation system for norvancomycin producing strain Amycolatopsis orientalis. J Microbiol. 2016;36:35–9.
Jeong H, Sim YM, Kim HJ, Lee DW, Lim SK, Lee SJ. Genome sequence of the vancomycin-producing Amycolatopsis orientalis subsp. orientalis strain KCTC 9412T. Genome Announc. 2013; 1:25.
Pelzer S, Sussmuth R, Heckmann D, Recktenwald J, Huber P, Jung G, Wohlleben W. Identification and analysis of the balhimycin biosynthetic gene cluster and its use for manipulating glycopeptide biosynthesis in Amycolatopsis mediterranei DSM5908. Antimicrob Agents Chemother. 1999;43:1565–73.
Sosio M, Stinchi S, Beltrametti F, Lazzarini A, Donadio S. The gene cluster for the biosynthesis of the glycopeptide antibiotic A40926 by nonomuraea species. Chem Biol. 2003;10:541–9.
Li TL, Huang F, Haydock SF, Mironenko T, Leadlay PF, Spencer JB. Biosynthetic gene cluster of the glycopeptide antibiotic teicoplanin: characterization of two glycosyltransferases and the key acyltransferase. Chem Biol. 2004;11:107–19.
Sosio M, Kloosterman H, Bianchi A, de Vreugd P, Dijkhuizen L, Donadio S. Organization of the teicoplanin gene cluster in Actinoplanes teichomyceticus. Microbiology. 2004;150:95–102.
Thaker MN, Wang W, Spanogiannopoulos P, Waglechner N, King AM, Medina R, Wright GD. Identifying producers of antibacterial compounds by screening for antibiotic resistance. Nat Biotechnol. 2013;31:922–7.
Nazari B, Forneris CC, Gibson MI, Moon K, Schramma KR, Seyedsayamdost MR. Nonomuraea sp. ATCC 55076 harbours the largest actinomycete chromosome to date and the kistamicin biosynthetic gene cluster. Medchemcomm. 2017;8:780–8.
Lei X, Zhang C, Jiang ZB, Li XX, Shi YY, Liu M, Xie YY, Wang LF, Hong B. Complete genome sequence of Amycolatopsis orientalis CPCC200066, the producer of norvancomycin. J Biotechnol. 2017;247:6–10.
Yim G, Thaker MN, Koteva K, Wright G. Glycopeptide antibiotic biosynthesis. J Antibiot (Tokyo). 2014;67:31–41.
Kaniusaite M, Tailhades J, Kittila T, Fage CD, Goode RJA, Schittenhelm RB, Cryle MJ: Understanding the early stages of peptide formation during the biosynthesis of teicoplanin and related glycopeptide antibiotics. FEBS J. 2020.
Xu L, Huang H, Wei W, Zhong Y, Tang B, Yuan H, Zhu L, Huang W, Ge M, Yang S, et al. Complete genome sequence and comparative genomic analyses of the vancomycin-producing Amycolatopsis orientalis. BMC Genomics. 2014;15:363.
Hubbard BK, Walsh CT. Vancomycin assembly: nature’s way. Angew Chem Int Ed Engl. 2003;42:730–65.
Alduina R, Tocchetti A, Costa S, Ferraro C, Cancemi P, Sosio M, Donadio S. A two-component regulatory system with opposite effects on glycopeptide antibiotic biosynthesis and resistance. Sci Rep. 2020;10:6200.
Kilian R, Frasch HJ, Kulik A, Wohlleben W, Stegmann E. The VanRS homologous two-component system VnlRSAb of the glycopeptide producer Amycolatopsis balhimycina activates transcription of the vanHAXSc genes in Streptomyces coelicolor, but not in A. balhimycina. Microb Drug Resist. 2016;22:499–509.
Beltrametti F, Consolandi A, Carrano L, Bagatin F, Rossi R, Leoni L, Zennaro E, Selva E, Marinelli F. Resistance to glycopeptide antibiotics in the teicoplanin producer is mediated by van gene homologue expression directing the synthesis of a modified cell wall peptidoglycan. Antimicrob Agents Chemother. 2007;51:1135–41.
Shawky RM, Puk O, Wietzorrek A, Pelzer S, Takano E, Wohlleben W, Stegmann E. The border sequence of the balhimycin biosynthesis gene cluster from Amycolatopsis balhimycina contains bbr, encoding a StrR-like pathway-specific regulator. J Mol Microbiol Biotechnol. 2007;13:76–88.
Alduina R, Lo Piccolo L, D’Alia D, Ferraro C, Gunnarsson N, Donadio S, Puglia AM. Phosphate-controlled regulator for the biosynthesis of the dalbavancin precursor A40926. J Bacteriol. 2007;189:8120–9.
Lo Grasso L, Maffioli S, Sosio M, Bibb M, Puglia AM, Alduina R. Two master switch regulators trigger A40926 biosynthesis in Nonomuraea sp. strain ATCC 39727. J Bacteriol. 2015;197:2536–44.
Horbal L, Kobylyanskyy A, Truman AW, Zaburranyi N, Ostash B, Luzhetskyy A, Marinelli F, Fedorenko V. The pathway-specific regulatory genes, tei15* and tei16*, are the master switches of teicoplanin production in Actinoplanes teichomyceticus. Appl Microbiol Biotechnol. 2014;98:9295–309.
Swint-Kruse L, Matthews KS. Allostery in the LacI/GalR family: variations on a theme. Curr Opin Microbiol. 2009;12:129–37.
Culp EJ, Waglechner N, Wang W, Fiebig-Comyn AA, Hsu YP, Koteva K, Sychantha D, Coombes BK, Van Nieuwenhze MS, Brun YV, Wright GD. Evolution-guided discovery of antibiotics that inhibit peptidoglycan remodelling. Nature. 2020;578:582–7.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Kumar CV, Coque JJ, Martin JF. Efficient transformation of the Cephamycin C producer Nocardia lactamdurans and development of shuttle and promoter-probe cloning vectors. Appl Environ Microbiol. 1994;60:4086–93.
Frith MC, Saunders NF, Kobe B, Bailey TL. Discovering sequence motifs with arbitrary insertions and deletions. PLoS Comput Biol. 2008;4:e1000071.
Bailey TL, Williams N, Misleh C, Li WW. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006;34:W369–73.
Kieser T, Bibb M, Buttner M, Chater K, Hopwood D. Pratical Streptomyces Gnetics. Norwich: The John Innes Foundation; 2000.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–8.
de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20:1453–4.
Saldanha AJ. Java Treeview–extensible visualization of microarray data. Bioinformatics. 2004;20:3246–8.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Nicol JW, Helt GA, Blanchard SG Jr, Raja A, Loraine AE. The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics. 2009;25:2730–1.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Li XX, Yu TF, He Q, McDowall KJ, Jiang BY, Jiang ZB, Wu LZ, Li GW, Li QL, Wang SM, et al. Binding of a biosynthetic intermediate to AtrA modulates the production of lidamycin by Streptomyces globisporus. Mol Microbiol. 2015;96:1257–71.
We thank CPCC (China Pharmaceutical Culture Collection) for providing the original strain CPCC 200066. We also thank Beijing Genomics Institute (Shenzhen, China) for Illumina Hiseq sequencing, assembly and annotation of the genome and transcriptome.
This work was supported by the National Natural Science Foundation of China (81703398, 81872780 and 81630089), the Drug Innovation Major Project of China (2018ZX09711001-007-001 and 2014ZX09201001-004-001), CAMS Innovation Fund for Medical Sciences (2016-I2M-3-012 and 2018-I2M-3-005).
Ethics approval and consent to participate
Consent for publication
All authors agreed to publish this article.
The authors declare no financial or commercial conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, X., Zhang, C., Zhao, Y. et al. Comparative genomics and transcriptomics analyses provide insights into the high yield and regulatory mechanism of Norvancomycin biosynthesis in Amycolatopsis orientalis NCPC 2-48. Microb Cell Fact 20, 28 (2021). https://doi.org/10.1186/s12934-021-01521-6