- Open Access
Comparative transcriptomic analysis reveals the significant pleiotropic regulatory effects of LmbU on lincomycin biosynthesis
Microbial Cell Factories volume 19, Article number: 30 (2020)
Lincomycin, produced by Streptomyces lincolnensis, is a lincosamide antibiotic and widely used for the treatment of the infective diseases caused by Gram-positive bacteria. The mechanisms of lincomycin biosynthesis have been deeply explored in recent years. However, the regulatory effects of LmbU that is a transcriptional regulator in lincomycin biosynthetic (lmb) gene cluster have not been fully addressed.
LmbU was used to search for homologous LmbU (LmbU-like) proteins in the genomes of actinobacteria, and the results showed that LmbU-like proteins are highly distributed regulators in the biosynthetic gene clusters (BGCs) of secondary metabolites or/and out of the BGCs in actinomycetes. The overexpression, inactivation and complementation of the lmbU gene indicated that LmbU positively controls lincomycin biosynthesis in S. lincolnensis. Comparative transcriptomic analysis further revealed that LmbU activates the 28 lmb genes at whole lmb cluster manner. Furthermore, LmbU represses the transcription of the non-lmb gene hpdA in the biosynthesis of l-tyrosine, the precursor of lincomycin. LmbU up-regulates nineteen non-lmb genes, which would be involved in multi-drug flux to self-resistance, nitrate and sugar transmembrane transport and utilization, and redox metabolisms.
LmbU is a significant pleiotropic transcriptional regulator in lincomycin biosynthesis by entirely activating the lmb cluster and regulating the non-lmb genes in Streptomyces lincolnensis. Our results first revealed the pleiotropic regulatory function of LmbU, and shed new light on the transcriptional effects of LmbU-like family proteins on antibiotic biosynthesis in actinomycetes.
Streptomycetes are famous for the ability to produce a great variety of valuable secondary metabolites, which have pharmacological activities such as antibacteria, antifungi, antiparasites, anticancer, and immunosuppression . Genes encoding for biosynthesis of secondary metabolites are usually located together into the biosynthetic gene clusters (BGCs) in streptomycetes. For timely and coordinated expression, BGCs are under strict controls at the transcription level of intertwined regulation by pleiotropic regulators and pathway specific regulators . The pathway specific regulators encoded by genes in BGCs of secondary metabolites activate or repress the expression of the BGCs for the biosynthesis of corresponding products. And the pleiotropic regulators associate various aspects involved in primary and secondary metabolisms, and the morphology of sporulation and pigmentation.
Lincomycin, produced by Streptomyces lincolnensis, is a lincosamide antibiotic used for the treatment of the infective diseases caused by Gram-positive bacteria . The lincomycin biosynthetic (lmb) gene cluster has been partially characterized in S. lincolnensis [4, 5]. And the genomes of S. lincolnensis NRRL 2936 and S. lincolnensis LC-G have been sequenced up to date [6, 7]. The lmb cluster harbors three resistance genes (lmrA, lmbrB and lmrC) and twenty-six open reading frames (ORFs) with putative biosynthetic and regulatory functions. The biosynthetic pathway of lincomycin (Fig. 1) has been deeply explored in recent years. Lincomycin consists of a methylthiolincosaminide (MTL) moiety and an N-methylated 4-propyl-l-proline (PPL) moiety. The moiety MTL originates from a GDP-activated C8 sugar lincosamide (LSM), which is synthesized from the condensation of C5 donor α-d-ribose-5-phosphate and C3 donor fructose-6-phosphate (or d-sedoheptulose-7-phosphate) by LmbR and LmbN . The biosynthesis of PPL begins with hydroxylation of l-tyrosine catalyzed by LmbB2 . The ergothioneine (EGT) S-conjugated LSM and the ATP-activated PPL are assembled together by the condensation enzymes [10,11,12,13,14]. Finally, lincomycin is synthesized after modifications of thiols exchange, N-methylation, C–S bond cleavage and S-methylation [15, 16].
For transcriptional regulation of lincomycin biosynthesis, the pleiotropic regulators have been reported in previous studies. The orphan response regulator GlnR showed to enhance lincomycin biosynthesis by up-regulating the transcription of nitrate-specific ABC transporter and nitrate assimilation genes, and lincomycin exporter gene lmrA . The TetR family regulator SLCG_2919 was identified to negatively regulate the transcription of the lmb cluster and its adjacent gene SLCG_2920, which encodes an ATP/GTP-binding protein associated to resistance against lincomycin . The bldA gene encoding the tRNA with rare leucine codon UUA regulated morphological differentiation and lincomycin biosynthesis via controlling the translation of TTA-containing genes including lmbB2 and lmbU in S. lincolnensis . The developmental regulator BldD played a positive regulation role in both the sporulation of S. lincolnensis and the transcription of the partial lmb genes . The A-factor AdpA was involved in the cascade regulation of lincomycin biosynthesis by activating some lmb genes and the bldA gene .
When the lmb cluster was firstly cloned, the lmbU gene was presumed to be a pathway specific regulatory gene . Until recently, LmbU was partially characterized as a positive transcriptional regulator involved in lincomycin biosynthesis by its inactivation and overexpression experiments . LmbU could activate several lmb genes by binding to DNA sites of the promoter regions through its helix-turn-helix motif . However, the in vivo regulatory effects of LmbU on lincomycin biosynthesis have not been fully understood. On the other hand, several homologous LmbU-like proteins have been reported as pathway specific regulators in BGCs of novobiocin [23, 24], himastatin , hormaomycin  and actinomycin D . SACE_5599 of Saccharopolyspora erythraea, a LmbU-like protein out of the erythromycin BGC, cross-regulates erythromycin production and the morphological differentiation . Furthermore, SACE_5599 could regulate lincomycin biosynthesis by directly binding to promoters of lmbA and lmbW in vitro acting similarly to LmbU .
In this study, to comprehensively understand the regulatory effects of LmbU, we first conducted the bioinformatic analysis which showed that LmbU homologous proteins were universally distributed in actinomycetes and belong to a novel regulator family. Then, we constructed the lmbU deletion strain to perform transcriptome sequencing. Comparative transcriptomic analysis showed LmbU is a pleiotropic regulator. It positively regulates the transcription of the 28 lmb genes at whole lmb cluster manner. Moreover, LmbU positively regulates nineteen non-lmb (out of lmb) genes and negatively regulates one non-lmb gene in S. lincolnensis genome.
In silico analysis of the putative pleiotropic regulatory function of LmbU
LmbU-like proteins with similarity over 49% were collected by BLASTP search against NCBI database using LmbU protein as query sequence. Then 93 LmbU-like proteins were selected based on their complete flanking sequences, which were used to locate lmbU-like genes either in the BGCs or not. The 16S rRNA genes from strains hosting these LmbU-like proteins were used to construct a phylogenetic tree (Fig. 2). LmbU-like proteins distribute in 72 species of 10 families of actinomycetes. More than half species (38 of 72) belong to genus Streptomyces, and others (34 of 72) belong to 21 genera from families Actinosynnemataceae, Pseudonocardiaceae, Micromonosporaceae, Frankiaceae, Streptosporangiaceae, Thermomonosporaceae, Nocardiopsaceae, Gaiellaceae and Solirubrobacteraceae. Sequence analysis of upstream and downstream regions around the lmbU-like genes revealed that nearly two-thirds of lmbU-like genes are located in the BGCs (mainly nonribosomal peptide BGCs, or polyketide-nonribosomal peptide BGCs) with single copy, similar to lmbU in the lmb cluster of S. lincolnensis. Furthermore, several BGCs have two lmbU-like genes. For example, the actinomycin C BGC of Streptomyces anulatus contains acmO and acmJ . One-third of lmbU-like genes are out of the BGCs, such as SACE_5599 is not in the erythromycin or any predicted BGC of Saccharopolyspora erythraea. It is interesting that some actinobacteria possess lmbU-like genes both in and out of the BGCs. For instance, Streptomyces lunaelactis strain MM109, a ferroverdin A (anticholesterol) producer, contains two lmbU-like genes, one in a predicted nonribosomal peptide BGC and another out of BGC . Thus, LmbU-like proteins are highly conserved and spread widely in actinobacteria, and they might be a new regulator protein family with pleiotropic regulatory roles in the biosynthesis of secondary metabolites and other biological process in actinomycetes.
Functionality of lmbU regulating lincomycin biosynthesis
To characterize the lmbU gene, the in-frame lmbU deletion strain SyBE2904 was constructed (Fig. 3a), which was confirmed by PCR (Fig. 3b), and further verified by DNA sequencing. The fermentation was conducted and lincomycin production in broth was measured by HPLC (Fig. 4A). As shown in Fig. 4B, compared to the original strain SyBE2901, the lmbU deletion strain SyBE2904 completely lost the ability to produce lincomycin during the fermentation, suggesting that inactivation of lmbU blocked the biosynthesis of lincomycin.
To complement the lmbU gene, three vectors (pHZ1358-U, pANT1201-U, and pLCY010-U) for the lmbU expression with different promoters and resistance genes (Table 1) were constructed and then introduced into the lmbU deletion strain SyBE2904, respectively. However, they all failed to resume lincomycin production. Analyzing the sequences of the lmbU gene and its downstream genes (lmbY and lmbX), the intergenic region between lmbU and lmbY is only 145-bp long, and the initiation codon of lmbX is overlapped with the termination codon of lmbY, which resulted in tetranucleotide sequence AUGA as stop-start codon. Thus, in-flame deletion of the lmbU gene might cause the polar effect on the expression of the lmbY and lmbX genes. The vector with fragment contains the lmbU-lmbY-lmbX genes under PermE* promoter was constructed and introduced into the lmbU deletion strain SyBE2904, and the resulting strain SyBE2913 restored the production of lincomycin (Fig. 4A). Moreover, the productivity of the lmbU overexpression strain SyBE2915 increased to 118.2% of the original titer of lincomycin (Fig. 4B). All results suggested that the lmbU gene encodes a positive regulator for activation of lincomycin biosynthesis.
Organization of transcriptional units in the lmb cluster
The transcription patterns of the 29 ORFs in the lmb cluster were complicated. Of which, four genes (lmbC, lmbD, lmbK, and lmbW) were considered to be transcribed as the monocistronic units, as deduced from their chromosomal arrangements (Fig. 5a). The rest genes of the lmb cluster could be transcribed as polycistronic units. The contiguous genes overlapped coding sequences (usually GUGA or AUGA as stop–start codon) were considered to be transcribed and translated together, including lmbB1-B2 (overlapping 4 bp), lmbF-G (overlapping 4 bp), lmbIH-J (overlapping 8 bp), lmbL-M–N-Z-P-O (overlapping 4 bp, respectively), lmbS-R (overlapping 4 bp), lmbQ-T-V (overlapping 4 bp, respectively), and lmbY-X (overlapping 4 bp) (Fig. 5a). Thus, there were at least 12 transcriptional units (TUs) for the 29 ORFs in the lmb cluster.
In order to identify TUs in the lmb cluster, strand specific RNA sequencing (ssRNA-Seq) was conducted from the original strain SyBE2901 at 3 day fermentation. The number of FPKM, fragments per kilobase of transcript sequence per millions base pairs sequenced, was used to estimate the expression level of each gene, and Rockhopper soft  was used to predict Transcription Start Site (TSS) paired with Transcription Termination Site (TTS) in the lmb cluster. The distribution of TTSs in the lmb cluster and the FPKM value of each gene were shown in Fig. 5a. According to the existence of contiguous overlapping reads, the distribution of TTSs and the degressive expression levels of following genes in one transcription unit, 17 TUs in the lmb cluster were preliminarily deduced.
For the subcluster lmbAB1B2, two TTSs were predicted at the upstreams of the lmbA and lmbB1B2 genes, respectively, and FPKM of the lmbB1 gene was much larger than that of the lmbA gene, demonstrating the lmbA gene and lmbB1B2 gene were transcribed separately. The distributions of TTSs and TSSs showed three TUs (lmbE, lmbG-F and lmbJ-IH) in the subcluster lmbJIHGFE. However, the FPKMs of lmbE and lmbG were not significantly different, indicating that lmbG-F-E might be co-transcribed as one TU. In the case of the compact subcluster lmbRSOPZNML, FPKMs were degressive from lmbV to lmbQ, lmbR to lmbZ, and lmbN to lmbL, respectively. Although TSS was predicated in the upstream of lmbO instead of lmbN, the FPKM values of lmbNML were much larger than those of lmbRSOPZ genes in their upstreams, indicating the transcription of lmbNML was separated from those of lmbRSOPZ. Thus, the subcluster lmbVTQRSOPZNML was predicted to be transcribed into three TUs (lmbV-T-Q, lmbR-S-O-P-Z and lmbN-M-L). For the subcluster lmbUYX, although the FPKMs of lmbU, lmbY, and lmbX were degressive following transcriptional direction, one TSS was predicted in the upstreams of lmbU and lmbY, respectively, suggesting that lmbU and lmbYX might be transcribed as two TUs.
In order to eliminate the non-determinacy of TSS prediction and verify the co-transcription of the subclusters lmbJIHGFE, lmbVTQRSOPZNML and lmbXYU, reverse transcription polymerase chain reaction (RT-PCR) was further performed. The gene specific antisense primers listed in Additional file 1: Table S2 were used to synthesize the first strand DNA (1st cDNA), and then the upstream genes were amplified by PCR to detect whether the genes were cotranscribed or not. As shown in Fig. 5b, in the subcluster lmbJIHGFE, using the 1st DNAlmbE (synthesized by antisense primer of lmbE) as template, the ds-cDNAs of the upstream genes lmbF and lmbG were detected, and the ds-cDNA of lmbIH was undetected. Furthermore, the ds-cDNA of lmbJ was detected using 1st DNAlmbIH as template, indicating the transcription of the subcluster lmbJIHGFE was divided into two TUs (lmbG-F-E and lmbJ-IH). In the case of the subcluster lmbRSOPZNML, the ds-cDNA of lmbN was detected while that of lmbZ was undetected using the 1st DNAlmbL as template, indicating the individual transcription of lmbNML. Meanwhile, the ds-cDNAs of lmbR and lmbO were detected using the 1st DNAlmbZ as template, indicating the co-transcription from lmbR to lmbZ genes. These results confirmed the subcluster lmbRSOPZNML was transcribed as two TUs (lmbN-M-L and lmbR-S-O-P-Z). A similar result was observed between lmbX and lmbU, indicating the co-transcription of lmbX, lmbY, and lmbU, which was consistent with the complementation to the lmbU deletion by expression of the lmbUYX. Taking together, 15 transcription units of the lmb cluster were identified, including eight monocistrons (TU1, TU2, TU4, TU5, TU8, TU12, TU13, TU15), two bicistrons (TU3, TU7), four tricistron (TU6, TU9, TU11, TU14) and one pentacistron (TU10) (Fig. 5c).
Comparison of transcriptomes between the original strain and the lmbU deletion strain
To obtain insight into the changes in gene expression levels between the original strain SyBE2901 and the lmbU deletion strain SyBE2904, RNA was isolated and subjected to whole-transcriptome sequencing via ssRNA-seq. Approximately 36.59 million 150 bp paired-end clean reads per sample were obtained after cleaning and checking the reads quality. Approximately 98% of clean reads were aligned uniquely to the S. lincolnensis genome. The expression level of each gene was calculated using FPKM. The correlation clustering among the two biological replicates of each sample was conducted based on the expression level of all genes. All biological replicates (two RNA samples for SyBE2901 and three RNA samples for SyBE2904) showed correlation coefficients over 0.9, indicating good reproducibility between biological replicates. To investigate transcriptional regulation effects of LmbU, the whole transcriptomes were compared.
Volcano plot was presented to identify genes with both high fold change and significance between the original strain SyBE2901 and the lmbU deletion strain SyBE2904 (Fig. 6). For high change, the foldchange of a gene between SyBE2904 and SyBE2901 should be > 1/2 and < 2, that means − 1 < log2(foldchange) < 1. For significance, the p value of a gene should be < 0.005, that means − log10(p value) > 2.3. Compared to the original strain SyBE2901, the transcription expressions of 49 genes were statistically different, of which 48 genes displayed decreased transcription levels, and only one gene displayed increased transcription level in the lmbU deletion strain SyBE2904. The results demonstrated that LmbU acted as a pleiotropic regulator for the positive or negative control expression of the genes including the lmb cluster and the non-lmb genes for lincomycin biosynthesis.
LmbU upregulates the transcription of all genes in the lmb cluster
Of 49 differential expression genes observed from ssRNA-Seq results, 29 genes are the lmb genes encompassing the entire lmb cluster. Figure 7a showed that the reads from strain SyBE2901 mapping to the lmb cluster were extremely more than those from SyBE2904, which demonstrated that the significant change of the transcriptional levels of all 29 lmb genes between these two strains. No reads mapped to the lmbU gene showing the successful deletion of the lmbU in strain SyBE2904. When the lmbU gene was inactive in strain SyBE2904, the transcriptional levels of the rest 28 lmb genes were declined by 93.71% to 99.98% compared to SyBE2901. These variation multiples could be also expressed by the − log2 foldchange value (Fig. 7b). The genes in entire lmb cluster were barely transcribed in SyBE2904, which were consistent with the results of semi-quantitative PCR analysis (Fig. 7c). PCR products of TU1, TU5, and TU8 were detected for 30 thermocycles instead of 26 thermocycles in strain SyBE2901, indicating that those four TUs were transcribed at low levels in strain SyBE2901. In the case of the other twelve TUs (TU2, TU3, TU4, TU6, TU7, TU9, TU10, TU11, TU12, TU13, TU14, and TU15), PCR products were all strongly detected in strain SyBE2901 for both 26 and 30 thermocycles, indicating that those TUs were transcribed at high levels in strain SyBE2901. However, PCR products of 15 TUs were not detected in strain SyBE2904 for 26 and 30 thermocycles, respectively, indicating that the transcriptions of all 15 TUs were completely blocked by the deletion of the lmbU gene. Taking together, LmbU positively and strongly regulates the transcription of the 28 lmb genes in the lmb cluster.
LmbU regulates transcription of non-lmb genes in S. lincolnensis genome
In addition to positive regulation on the expression of the whole lmb cluster, the results of RNA-Seq analysis showed that twenty non-lmb genes are regulated by LmbU (Table 2). Of which, eight genes were selected to perform semi-quantitative PCR analysis (Additional file 1: Fig. S1), and the results were consistent with those of ssRNA-Seq.
SLCG_3526 encodes a 4-hydroxyphenylpyruvate dioxygenase (HpdA), which is an Fe(II)-containing non-heme enzyme, and catalyzes the oxidation and decarboxylation of 4-hydroxyphenylpyruvate (precursor of l-tyrosine) to homogentisic acid . The transcript of hpdA in the original strain was minor, while significantly increased after the lmbU deletion in strain SyBE2904, indicating that LmbU negatively regulated the hpdA expression and might be involved in the biosynthetic pathway of l-tyrosine via blocking the degradation of precursor 4-hydroxyphenylpyruvate.
In contrast to the hpdA gene, the nineteen non-lmb genes were transcriptionally decreased in strain SyBE2904 after the lmbU deletion, compared with the original strain, which indicated that LmbU positively regulates those genes. Based on biochemical and cellular functions, the nineteen non-lmb genes could be classified into four groups. Group 1 genes encode the transporters , including three major facilitator superfamily (MFS) transporters (SLCG_3493, SLCG_6997, and SLCG_7561), four ATP-binding cassette (ABC) superfamily transporters (SLCG_0089, SLCG_7139, SLCG_7140, and SLCG_7316), and a drug/metabolite transporter (DMT) superfamily protein (SLCG_1725). Group 2 genes encode the five redox proteins and enzymes, including a flavoprotein (SLCG_3218), a succinate dehydrogenase (SLCG_7376), an enoyl reductase (SLCG_2976), an antibiotic biosynthetic monooxygenase (SLCG_6211), and a SgcJ/EcaC family oxidoreductase (SLCG_7612). Group 3 genes encode cellular metabolic enzymes, including an epimerase (SLCG_1534), an acetyl-CoA synthetase (SLCG_6998), an acetolactate synthase (SLCG_7000), and a carbonic anhydrase (SLCG_0221). The last group 4 genes encode two hypothetical proteins (SLCG_0224 and SLCG_1709) without definite properties deduced from the NCBI database. To sum up, the results indicated LmbU would regulate the degradation, transport, oxidoreduction, and anabolism for lincomycin biosynthesis.
Homologous lmbU-like regulators are highly conserved proteins, and their genes widely co-exist in the BGCs of secondary metabolites or are distributed out of the BGCs, or both of them in actinomycetes which suggests that LmbU could be a new family of the transcription regulators, targeting the specific pathways or/and the other cellular metabolisms. Here, our comparative transcriptomics revealed that LmbU is a pleiotropic regulator via activating the transcription of the whole lmb cluster and regulating the transcription of non-lmb genes that scattered in the chromosome of S. lincolnensis, both involved in lincomycin biosynthesis.
The positive regulation of LmbU on the lincomycin biosynthesis was firstly confirmed by the lmbU deletion which abolished lincomycin production, and validated by the lmbU complementation and overexpression which restored and increased the titers of lincomycin, while no obvious morphology changes were observed after genetic manipulation in S. lincolnensis. It was consistent with the previous reports on the lmbU deletion and overexpression strains . To obtain the insight on the regulatory properties of LmbU on lincomycin biosynthesis, analysis of comparative transcriptomes between the original strain and the lmbU deletion strain was performed in this study. Unexpectively, LmbU positively regulates all of the 28 lmb genes as the regulatory targets, including the structural genes encoding enzymes for catalyzing biosynthesis of lincomycin, and resistance genes (lmrA, lmrB, and lmrC) in the lmb cluster. Several LmbU-like proteins in the BGCs also have been reported with positively regulatory effects. NovE activated the pathway-specific regulator gene novG for novobiocin biosynthesis [23, 24, 34]. And HmtD up-regulated the transcription of NRPS gene hmtI for himastatin biosynthesis at a fermentation time-dependent manner . Different from those reported LmbU-like proteins, LmbU in S. lincolnensis acts as a positive cluster-scale regulator for lincomycin biosynthesis.
In addition to the lmb cluster, LmbU positively or negatively regulates the transcription of non-lmb genes in S. lincolnensis. 4-Hydroxyphenylpyruvate is precursor substrate for biosynthesis of l-tyrosine, which is hydroxylated by LmbB2 into l-dopa and initiates the biosynthesis of PPL moiety of lincomycin . However, 4-hydroxyphenylpyruvate could be degraded by catalysis of HpdA  that was down-regulated by LmbU, thus it would lead to maintenance of l-tyrosine availability for PPL biosynthesis. In silico genomic sequence analysis of S. lincolnensis showed that the hpdA and hpaR (SLCG_3525) genes are clustered as the hpd operon. HpdR is an AsnC family transcriptional regulator and activates the hpdA transcription in the presence of 4-hydroxyphenylpyruvate and l-tyrosine . In this study, the transcriptional level of hpdR was low and no obviously different between the original strain SyBE2901 and the lmbU deletion strain SyBE2904 (log2(FoldChange) = 0.55136), suggesting that the inhibitive effect of LmbU on the hpdA transcription might be HpdR-independent. It was interesting that LmbU not only represses the hpdA expression to increase l-tyrosine pool through inhibiting degradation of 4-hydroxyphenylpyruvate, but also promotes the expression of lmbB1B2. Both of them are benefit for PPL biosynthesis for l-tyrosine. Thus, LmbU plays at least a dual regulatory role at metabolic node of l-tyrosine as a substrate for lincomycin biosynthesis.
Absorptive ability of mycelia to the carbon and nitrogen sources from the fermentation medium is crucial for cell growth and lincomycin biosynthesis. LmbU facilitated the expression of the nitrate/nitrite transporter SLCG_3493 and sugar-binding protein SLCG_0089, which would force nitrate and sugar uptake and utilization. The similar function were reported on the nitrate transporters regulated by GlnR , suggesting that SLCG_3493 may be involved in lincomycin biosynthesis. SLCG_6997 similar to an oxalate:formate antiporter  and SLCG_1725 similar to threonine and homoserine efflux transporter  might be helpful for exchange of primary metabolites including organic acids and amino acids. Four proteins SLCG_7139, SLCG_7140, SLCG_7316, and SLCG_7561 could possess the multidrug efflux transport function. Especially, SLCG_7561 is a homolog of LmrA (48% of similarity in amino acid sequence), which also widely exists in lincomycin-resistant bacteria . In addition to the resistance conferred by lmrA, lmrB and lmrC in the lmb cluster, the elevated expression of the non-lmb multidrug efflux transporter genes by up-regulation of LmbU might contribute higher self-resistance of S. lincolnensis to lincomycin and lead to the great amount of lincomycin production during the fermentation.
The steady-state redox potential is essential for primary and secondary metabolisms. Dehydrogenase SLCG_7376 might serve as a vital link between the tricarboxylic acid cycle and oxidative phosphorylation , and SLCG_2976, SLCG_6211 and SLCG_7612 would be involved in the reductive bioprocess . SLCG_3218 functions as NADPH-dependent FMN reductase, and provides reducing power NADPH. LmbU up-regulates those oxidoreductase genes, and participates in cellular redox metabolisms.
Although majority of regulator genes in clusters are identified as pathway specific regulators of their cognate BGCs, we demonstrated that LmbU plays regulatory roles both in the lmb cluster and in non-lmb genes as aforementioned, exhibiting pleiotropic regulation. A few in-cluster regulators have been reported to have such regulatory effects on the BGCs and/or genes out of the BGCs. The PAS-LuxR family regulator FscRI within the candicidin BGC, not only activated genes in its cluster, but also provided a direct and essential regulation for expression of the antimycin BGC . Moreover, apart from the cross-regulation of filipin and oligomycin biosynthesis, PteF in the pentaene filipin BGC could regulate DNA replication and repair, carbohydrate, lipid and energy metabolisms and morphological differentiation . These results showed that the pleiotropic regulatory effects of the in-BGC regulators maybe not rare, and will be more frequently discovered following high throughput RNA-Seq technologies.
The results of this study demonstrated that LmbU is a pleiotropic transcriptional regulator on lincomycin biosynthesis. In addition to direct activation of the whole lmb cluster, the regulatory effects of LmbU would be involved in l-tyrosine biosynthesis, nitrate and sugar transmembrane transport and utilization, multi-drug resistance, and oxidoreduction metabolisms by regulating the genes out of the lmb cluster. Our results provided evidence to elucidate the regulatory functions of LmbU-like proteins in actinomycetes.
Phylogenetic analyses of actinomycetes with LmbU-like proteins
LmbU-like proteins as were collected by BLASTP search against NCBI database. The flanking sequences of homologous lmbU genes (from upstream 50 kb to downstream 50 kb) were used to predict secondary metabolite BGCs by BLASTP and antiSMASH 5.0 . The 16S rRNA genes of the host strains containing homologous LmbU proteins were used to construct phylogenetic tree alignments and phylogenetic analysis was performed using MEGA7  by the neighbor-joining method (Kimura 2-parameter model + G)  and 500 bootstrap replications. The GenBank accession numbers of 16S sRNA genes and homologous LmbU proteins were listed in Additional file 1: Table S1.
Strains, plasmids and culture conditions
The bacterial strains and plasmids used in this study are described in Table 1. Primers used in this study are listed in Additional file 1: Table S2. The sporulation and fermentation of S. lincolnensis strains were carried out as previously described . Briefly, the spores of S. lincolnensis were routinely cultivated on the modified Gauze’s Medium No.1 (20 g l−1 soluble starch, 5 g l−1 soybean flour, 1 g l−1 KNO3, 0.5 g l−1 NaCl, 0.5 g l−1 MgSO4, 0.5 g l−1 K2HPO4, 0.01 g l−1 FeSO4, 18 g l−1 agar, pH 7.2 to 7.4) for 7 days at 30 °C. For shake flask fermentation, the spores inoculated into 25 ml of seed medium (20 g l−1 soluble starch, 10 g l−1 glucose, 10 g l−1 soybean, 30 g l−1 cream corn, 1.5 g l−1 (NH4)2SO4, 4 g l−1 CaCO3, pH 7.1) and grown for 2 days at 30 °C, 250 rpm. Then 2 ml of seed culture was added into a 250-ml flask containing 25 ml of fermentation medium (100 g l−1 glucose, 20 g l−1 soybean, 1.5 g l−1 cream corn, 8 g l−1 NaNO3, 5 g l−1 NaCl, 6 g l−1 (NH4)2SO4, 0.3 g l−1 K2HPO4, 8 g l−1 CaCO3, pH 7.1) and cultivated for 7 days at 30 °C and 250 rpm. Appropriate antibiotics were added in the medium when necessary.
Inactivation and complementation of lmbU
For the lmbU in-frame deletion, the vector pAP07 was first constructed by the homologous recombination method. The 3.06-kb PCR fragment from the upstream (1105 bp) to the downstream (1275 bp) of lmbU was amplified from the genomic DNA of S. lincolnensis (isolated by Kirby mix procedure ) using primers DU-F and DU-R (Additional file 1: Table S2), digested by EcoR I, and cloned into SuperCos1 to yield vector pCosU. Using pIJ773 as template, the fragment acc(3)IV-oriT with 39-bp homologous sequences of lmbU at each terminus was amplified by primers QU-F and QU-R, and transformed into E. coli BW25113/pIJ790/pCosU. The lmbU gene on pCosU was replaced with acc(3)IV-oriT through the homologous recombination, resulting vector pAP07. To delete lmbU, the vector pAP07 was introduced into the original strain SyBE2901 by intergeneric conjugation from E. coli ET12567/pUZ8002/pAP07, following the procedure described previously . The single-crossover mutants were identified with apramycin resistance, and then cultivated for several rounds without antibiotics. The double-crossover mutants for the lmbU deletion with kanamycin sensitivity and apramycin resistance were selected and confirmed by PCR using primers QU-F2 and QU-R2. Further sequencing verified that the 656-bp coding sequence of the lmbU gene on the genome was replaced with acc(3)IV-oriT, resulting the strain SyBE2904 with the in-frame deletion of lmbU.
For complementation of lmbU, the 1454-bp fragment containing the lmbU gene was amplified using 10U-F and 10U-R as primers, and ligated to pLCY010, pHZ1358, and pANT1201, resulting pLCY0010-U, pHZ1358-U and pANT1201-U, respectively. The 3245-bp fragment containing lmbUYX genes was amplified and ligated to pLCY010, resulting pLCY010-UYX. The vectors pLCY0010-U, pHZ1358-U pANT1201-U and pLCY010-UYX were introduced into the lmbU deletion strain SyBE2904, generating the complementation strain SyBE2929, SyBE2930, SyBE2931 and SyBE2913, respectively. The vector pLCY0010-U was introduced into the original strain SyBE2901, generating the lmbU overexpression strain SyBE2915.
Analysis of lincomycin in fermentation broth
After cultivation, the supernatant was collected from fermentation broth (8000 rpm for 10 min), then mixed with 1.5 volumes of methanol and centrifuged again to remove the residue. The sample filtered through 0.22 μm of nylon membrane was subjected to HPLC analysis (Agilent 1200, USA) on C18 column (4.6 × 250 mm, 5 μm, Agela Technology, China) at UV 214 nm, 25 °C. The mobile phase was 50% ammonium acetate solution (0.005 mol l−1, pH 9.0) and 50% methanol. The concentration of lincomycin was estimated by HPLC using the lincomycin standard (National Institute for the Control of Pharmaceutical and Biological Products, China).
RNA-Seq transcriptomic analysis
Because of the complexity of bacterial transcriptome, we used strand specific RNA sequencing (ssRNA-Seq) to investigate transcription of S. lincolnensis. The total RNA was extracted from S. lincolnensis after growth for 3 days. To enrich mRNA, rRNA was removed by Ribo-zero kit. The strand-specific libraries for high throughput RNA sequencing were generated using NEBNext Ultra Directional RNA Library Prep Kit for Illumia (NEB). Then fragmentation buffers were added to break the mRNA into short fragments. The 1st cDNA was synthesized with random primers using the mRNA as template. The buffer solution, dNTPs (dTTP was replaced by dUTP), DNA polymerase I and RNase H were added to synthesize the double-stranded cDNA. After been purified by AMPure XP beads, the double-stranded cDNAs were repaired at the end, attached with the A-tail, ligated with Y-shape adapter and selected with the fragment size by AMPure XP beads. And then the “U” containing strand of the adapter ligated DNA was degraded by USER-enzyme. Finally, the strand-specific libraries were enriched by PCR amplification, and then deep sequenced on an Illumina Hiseq 2500 platform. The raw sequence data were filtered by removing reads containing adapter, reads containing more than 10% N (undetermined base), and low-quality reads (reads containing more than 50% bases with mass value Qpred < 5). The clean reads were aligned with the genome of S. lincolnensis LC-G (GenBank ID: 1435096411) by Bowtie2-2.2.3 . The number of FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) was used to quantify gene expression. Genes with |log2(FoldChange)| > 1 and p value < 0.005 found by DEGSeq were assigned as differentially expressed.
Semi-quantitative PCR analysis
For RNA extraction, mycelia of S. lincolnensis original strain SyBE2901 and the lmbU deletion strain SyBE2904 were collected from 25-ml culture (3 days of growth in fermentation medium) by centrifugation, and ground into powder in liquid nitrogen. Total RNA was extracted using RNApure kit (Bio Teke, China) according to the manufacturer’s protocol.
Reverse transcription was conducted using TransScript II One-Step gDNA Removal and cDNA Synthesis SuperMix kit (TransGene Biotech, China). Using gene specific primers (reverse primer used in semi-quantitative PCR) and 500 ng of total RNA as template, the first strand was generated. For analyzing transcription, the 1st cDNA reaction mixture was used as the template to amplify ds-cDNA in the following semi-quantitative PCR. In each case the experiments were conducted in triplicate. The PCR products were detected using agarose gel electrophoresis, and then exposed under UV to analyze relative intensity. Primers used in semi-quantitative PCR experiments were listed in Additional file 1: Table S2.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Biosynthetic gene cluster
- lmb cluster:
Lincomycin biosynthetic gene cluster
Open reading frames
Strand specific RNA sequencing
Transcription start site
Transcription termination site
Reverse transcription polymerase chain reaction
- 1st cDNA:
First strand DNA
- MFS transporter:
Major facilitator superfamily transporter
Barka EA, Vatsa P, Sanchez L, Gaveau-Vaillant N, Jacquard C, Meier-Kolthoff JP, et al. Taxonomy, physiology, and natural products of actinobacteria. Microbiol Mol Biol Rev. 2015;80(1):1–43.
van der Heul HU, Bilyk BL, McDowall KJ, Seipke RF, van Wezel GP. Regulation of antibiotic production in Actinobacteria: new perspectives from the post-genomic era. Nat Prod Rep. 2018;35(6):575–604.
Spížek J, Řezanka T. Lincomycin, cultivation of producing strains and biosynthesis. Appl Microbiol Biotechnol. 2004;63(5):510–9.
Peschke U, Schmidt H, Zhang HZ, Piepersberg W. Molecular characterization of the lincomycin-production gene cluster of Streptomyces lincolnensis 78-11. Mol Microbiol. 1995;16(6):1137–56.
Koběrská M, Kopecký J, Olšovská J, Jelinkova M, Ulanova D, Man P, et al. Sequence analysis and heterologous expression of the lincomycin biosynthetic cluster of the type strain Streptomyces lincolnensis ATCC 25466. Folia Microbiol. 2008;53(5):395–401.
Hou B, Lin Y, Wu H, Guo M, Petkovic H, Tao L, et al. The novel transcriptional regulator LmbU promotes lincomycin biosynthesis through regulating expression of its target genes in Streptomyces lincolnensis. J Bacteriol. 2017;200(2):e00447-17.
Xu Y, Tan G, Ke M, Li J, Tang Y, Meng S, et al. Enhanced lincomycin production by co-overexpression of metK1 and metK2 in Streptomyces lincolnensis. J Ind Microbiol Biotechnol. 2018;45(5):345–55.
Sasaki E, Lin CI, Lin KY, Liu HW. Construction of the octose 8-phosphate intermediate in lincomycin A biosynthesis: characterization of the reactions catalyzed by LmbR and LmbN. J Am Chem Soc. 2012;134(42):17432–5.
Neusser D, Schmidt H, Spizèk J, Novotnà J, Peschke U, Kaschabeck S, et al. The genes lmbB1 and lmbB2 of Streptomyces lincolnensis encode enzymes involved in the conversion of l-tyrosine to propylproline during the biosynthesis of the antibiotic lincomycin A. Arch Microbiol. 1998;169(4):322–32.
Kadlčík S, Kučera T, Chalupská D, Gažák R, Koběrská M, Ulanová D, et al. Adaptation of an l-proline adenylation domain to use 4-propyl-l-proline in the evolution of lincosamide biosynthesis. PLoS ONE. 2013;8(12):e84902.
Janata J, Kadlcik S, Koberska M, Ulanova D, Kamenik Z, Novak P, et al. Lincosamide synthetase—a unique condensation system combining elements of nonribosomal peptide synthetase and mycothiol metabolism. PLoS ONE. 2015;10(3):e0118850.
Zhao Q, Wang M, Xu D, Zhang Q, Liu W. Metabolic coupling of two small-molecule thiols programs the biosynthesis of lincomycin A. Nature. 2015;518(7537):115.
Wang M, Zhao Q, Zhang Q, Liu W. Differences in PLP-dependent cysteinyl processing lead to diverse S-functionalization of lincosamide antibiotics. J Am Chem Soc. 2016;138(20):6348–51.
Ushimaru R, Lin CI, Sasaki E, Liu HW. Characterization of enzymes catalyzing transformations of cysteine S-conjugated intermediates in the lincosamide biosynthetic pathway. ChemBioChem. 2016;17(17):1606–11.
Najmanová L, Kutejová E, Kadlec J, Polan M, Olšovská J, Benada O, et al. Characterization of N-demethyllincosamide methyltransferases LmbJ and CcbJ. ChemBioChem. 2013;14(17):2259–62.
Kamenik Z, Kadlcik S, Radojevic B, Jiraskova P, Kuzma M, Gazak R, et al. Deacetylation of mycothiol-derived ‘waste product’ triggers the last biosynthetic steps of lincosamide antibiotics. Chem Sci. 2016;7(1):430–5.
Meng S, Wu H, Wang L, Zhang B, Bai L. Enhancement of antibiotic productions by engineered nitrate utilization in actinomycetes. Appl Microbiol Biotechnol. 2017;101(13):5341–52.
Xu Y, Ke M, Li J, Tang Y, Wang N, Tan G, et al. TetR-type regulator SLCG_2919 is a negative regulator of lincomycin biosynthesis in Streptomyces lincolnensis. Appl Environ Microbiol. 2018;85(1):e02091-18.
Hou B, Tao L, Zhu X, Wu W, Guo M, Ye J, et al. Global regulator BldA regulates morphological differentiation and lincomycin production in Streptomyces lincolnensis. Appl Microbiol Biotechnol. 2018;102(9):4101–15.
Li J, Wang N, Tang Y, Cai X, Xu Y, Liu R, et al. Developmental regulator BldD directly regulates lincomycin biosynthesis in Streptomyces lincolnensis. Biochem Biophys Res Commun. 2019;518(3):548–53.
Kang Y, Wang Y, Hou B, Wang R, Ye J, Zhu X, et al. AdpAlin, a pleiotropic transcriptional regulator, is involved in the cascade regulation of lincomycin biosynthesis in Streptomyces lincolnensis. Front Microbiol. 2019;10:2428.
Hou B, Zhu X, Kang Y, Wang R, Wu H, Ye J, et al. LmbU, a Cluster-situated regulator for lincomycin, consists of a DNA-binding domain, an auto-inhibitory domain, and forms homodimer. Front Microbiol. 2019;10:989.
Eustáquio AS, Luft T, Wang ZX, Gust B, Chater KF, Li SM, et al. Novobiocin biosynthesis: inactivation of the putative regulatory gene novE and heterologous expression of genes involved in aminocoumarin ring formation. Arch Microbiol. 2003;180(1):25–32.
Dangel V, Härle J, Goerke C, Wolz C, Gust B, Pernodet JL, et al. Transcriptional regulation of the novobiocin biosynthetic gene cluster. Microbiology. 2009;155(12):4025–35.
Cai X, Teta R, Kohlhaas C, Crüsemann M, Ueoka R, Mangoni A, et al. Manipulation of regulatory genes reveals complexity and fidelity in hormaomycin biosynthesis. Chem Biol. 2013;20(6):839–46.
Xie Y, Li Q, Qin X, Ju J, Ma J. Enhancement of himastatin bioproduction via inactivation of atypical repressors in Streptomyces hygroscopicus. Metab Eng Commun. 2019;8:e00084.
Liu M, Jia Y, Xie Y, Zhang C, Ma J, Sun C, et al. Identification of the actinomycin D biosynthetic pathway from marine-derived Streptomyces costaricanus SCSIO ZS0073. Mar Drugs. 2019;17(4):240.
Kirm B, Magdevska V, Tome M, Karničar K, Petek M, Vidmar R, et al. SACE_5599, a putative regulatory protein, is involved in morphological differentiation and erythromycin production in Saccharopolyspora erythraea. Microb Cell Fact. 2013;12:126.
Keller U, Lang M, Crnovcic I, Pfennig F, Schauwecker F. The actinomycin biosynthetic gene cluster of Streptomyces chrysomallus: a genetic hall of mirrors for synthesis of a molecule with mirror symmetry. J Bacteriol. 2010;192(10):2583–95.
Maciejewska M, Pessi IS, Arguelles-Arias A, Noirfalise P, Luis G, Ongena M, et al. Streptomyces lunaelactis sp. nov., a novel ferroverdin A-producing Streptomyces species isolated from a moonmilk speleothem. Antonie van Leeuwenhoek. 2015;107(2):519–31.
McClure R, Balasubramanian D, Sun Y, Bobrovskyy M, Sumby P, Genco CA, et al. Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res. 2013;41(14):e140.
Crouch NP, Lee MH, Iturriagagoitia-Bueno T, MacKinnon CH. Cloning, expression, and purification of mammalian 4-hydroxyphenylpyruvate dioxygenase/alpha-ketoisocaproate dioxygenase. Methods Enzymol. 2000;324:342–55.
Saier MH, Reddy VS, Tsu BV, Ahmed MS, Li C, Moreno-Hagelsieb G. The Transporter Classification Database (TCDB): recent advances. Nucleic Acids Res. 2016;44:D372–9.
Dangel V, Eustáquio AS, Gust B, Heide L. novE and novG act as positive regulators of novobiocin biosynthesis. Arch Microbiol. 2008;190(5):509–19.
Loprasert S, Whangsuk W, Dubbs JM, Sallabhan R, Somsongkul K, Mongkolsuk S. HpdR is a transcriptional activator of Sinorhizobium meliloti hpdA, which encodes a herbicide-targeted 4-hydroxyphenylpyruvate dioxygenase. J Bacteriol. 2007;189(9):3660–4.
Heymann JA, Sarker R, Hirai T, Milne JL, Maloney PC, Subramaniam S. Projection structure and molecular architecture of OxlT, a bacterial membrane transporter. EMBO J. 2001;20(16):4408–13.
Livshits VA, Zakataeva NP, Aleshin VV, Vitushkina MV. Identification and characterization of the new gene rhtA involved in threonine and homoserine efflux in Escherichia coli. Res Microbiol. 2003;154(2):123–35.
van Veen HW, Margolles A, Müller M, Higgins CF, Konings WN. The homodimeric ATP-binding cassette transporter LmrA mediates multidrug transport by an alternating two-site (two-cylinder engine) mechanism. EMBO J. 2000;19(11):2503–14.
Cook GM, Greening C, Hards K, Berney M. Energetics of pathogenic bacteria and opportunities for drug development. Adv Microb Physiol. 2014;65:1–62.
Huang T, Chang CY, Lohman JR, Rudolf JD, Kim Y, Chang C, et al. Crystal structure of SgcJ, an NTF2-like superfamily protein involved in biosynthesis of the nine-membered enediyne antitumor antibiotic C-1027. J Antibiot. 2016;69(10):731–40.
McLean TC, Hoskisson PA, Seipke RF. Coordinate regulation of antimycin and candicidin biosynthesis. mSphere. 2016;1(6):e00305–16.
Vicente CM, Payero TD, Santos-Aberturas J, Barreales EG, de Pedro A, Aparicio JF. Pathway-specific regulation revisited: cross-regulation of multiple disparate gene clusters by PAS-LuxR transcriptional regulators. Appl Microbiol Biotechnol. 2015;99(12):5123–35.
Blin K, Shaw S, Steinke K, Villebro R, Ziemert N, Lee SY, et al. antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019;47(W1):W81–7.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.
Pang A-P, Du L, Lin C-Y, Qiao J, Zhao G-R. Co-overexpression of lmbW and metK led to increased lincomycin A production and decreased byproduct lincomycin B content in an industrial strain of Streptomyces lincolnensis. J Appl Microbiol. 2015;119(4):1064–74.
Gust B, Chandra G, Jakimowicz D, Tian Y, Bruton CJ, Chater KF. Red-mediated genetic manipulation of antibiotic-producing Streptomyces. Adv Appl Microbiol. 2004;54:107–28.
Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA. Practical streptomyces genetics. Norwich: John Innes Foundation; 2000.
Du L, Liu R-H, Ying L, Zhao G-R. An efficient intergeneric conjugation of DNA from Escherichia coli to mycelia of the lincomycin-producer Streptomyces lincolnensis. Int J Mol Sci. 2012;13(4):4797–806.
Sun Y, Zhou X, Liu J, Bao K, Zhang G, Tu G, et al. ‘Streptomyces nanchangensis’, a producer of the insecticidal polyether antibiotic nanchangmycin and the antiparasitic macrolide meilingmycin, contains multiple polyketide gene clusters. Microbiology. 2002;148(2):361–71.
DeSanti CL, Strohl WR. Characterization of the Streptomyces sp. strain C5 snp locus and development of snp-derived expression vectors. Appl Environ Microbiol. 2003;69(3):1647–54.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Semi-quantitative PCR of the non-lmb genes in strains SyBE2901 and SyBE2904 for 30 and 34 thermocycles, respectively. In each case the experiments were repeated three times. Table S1. The actinomycetes with 16S rRNA genes and LmbU-like proteins used for constructing the phylogenetic trees in Fig. 1. Table S2. Primers used in this study.
About this article
Cite this article
Lin, C., Pang, A., Zhang, Y. et al. Comparative transcriptomic analysis reveals the significant pleiotropic regulatory effects of LmbU on lincomycin biosynthesis. Microb Cell Fact 19, 30 (2020). https://doi.org/10.1186/s12934-020-01298-0
- Streptomyces lincolnensis
- Transcriptional regulator
- Pleiotropic regulation