Plackett-Burman design to improve lincomycin A production in S. lincolnensis fermentation
PBD was developed using Design-Expert software (Wang et al. 2019). The design matrix selected for the screening of significant variables for lincomycin A and lincomycin B production (Table 2). The confidence level was set at 5%; therefore, the variables that scored a probability (p) value less than 0.05 were considered significantly influential factors that could affect lincomycin production. The statistical analysis of the model was used an analysis of variance (ANOVA) to evaluate the significance and effectiveness of the design. The ANOVA for lincomycin A content showed that soybean powder, corn steep liquor, glucose, and CaCO3 were components of the culture medium that significantly affected lincomycin A production, as their p-values were less than 0.05 (Table S2). ANOVA analysis of lincomycin B content showed that soybean powder, corn steep liquor, glucose, and NaNO3 were components that had a significant effect on lincomycin A production, as their p-values were also less than 0.05 (Table S2). Lincomycin A and lincomycin B content was significantly affected in media that contained soybean powder, corn steep liquor, and glucose. Thus, the important factors identified by the initial Plackett-Burman screening method, which influenced lincomycin A and lincomycin B content, were soybean powder, corn steep liquor, and glucose.
Table 2
Results of Plackett-Burman experiment
Run | X1 | X2 | X3 | X4 | X5 | X6 | X7 | X8 | Lincomycin A (mg/L) | Lincomycin B (%) |
0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 3555 | 4.23 |
1 | + 1 | -1 | + 1 | + 1 | -1 | -1 | -1 | + 1 | 2884 | 1.54 |
2 | -1 | + 1 | + 1 | + 1 | + 1 | + 1 | -1 | -1 | 2625 | 1.89 |
3 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | 3706 | 5.43 |
4 | + 1 | -1 | -1 | + 1 | + 1 | + 1 | + 1 | + 1 | 3937 | 4.36 |
5 | + 1 | + 1 | + 1 | + 1 | -1 | -1 | + 1 | -1 | 3083 | 1.08 |
6 | -1 | + 1 | -1 | + 1 | + 1 | -1 | -1 | + 1 | 3124 | 2.15 |
7 | + 1 | + 1 | -1 | -1 | -1 | + 1 | -1 | + 1 | 1837 | 1.66 |
8 | + 1 | -1 | + 1 | -1 | + 1 | + 1 | -1 | -1 | 1568 | 2.62 |
9 | -1 | -1 | -1 | -1 | -1 | + 1 | + 1 | -1 | 4364 | 4.69 |
10 | + 1 | + 1 | -1 | -1 | + 1 | -1 | + 1 | -1 | 2624 | 4.38 |
11 | -1 | -1 | + 1 | -1 | + 1 | -1 | + 1 | + 1 | 4542 | 4.83 |
12 | -1 | + 1 | + 1 | -1 | -1 | + 1 | + 1 | + 1 | 4439 | 4.62 |
The values were obtained from three independent experiments. |
Steepest ascent improved culture conditions for optimizing lincomycin A production
Because soybean powder and corn steep liquor exerted negative effects on lincomycin A and glucose exerted a positive effect (Table S2), the direction of steepest ascent indicated that the concentration of glucose should increase and the concentrations of soybean powder and corn steep liquor should decrease for optimal experimental conditions that maximize lincomycin A production. Five sets of experiments using the steepest ascent and corresponding experimental results showed that the optimum value area was located in group Run3 (Table 3).
Table 3
Three-factors steepest ascent experiment
Run | X1 (g/L) | X2 (g/L) | X7 (g/L) | Lincomycin A (mg/L) | Lincomycin B (%) |
0 | 25 | 2 | 100 | 3536 ± 82 | 4.4 ± 0.2 |
1 | 15 | 1.2 | 140 | 3421 ± 102 | 5.3 ± 0.2 |
2 | 20 | 1.3 | 133 | 3875 ± 48 | 2.8 ± 0.1 |
3 | 25 | 1.4 | 126 | 4606 ± 101 | 0.8 ± 0.1 |
4 | 30 | 1.5 | 119 | 3671 ± 152 | 2.4 ± 0.1 |
5 | 35 | 1.6 | 112 | 3415 ± 72 | 3.7 ± 0.2 |
6 | 40 | 1.7 | 105 | 2299 ± 96 | 3.6 ± 0.1 |
7 | 45 | 1.8 | 98 | 1029 ± 22 | 7.2 ± 0.3 |
The values were obtained from three independent experiments. |
These results helped define the optimal concentrations as 25 g/L soybean powder, 1.4 g/L corn steep liquor, and 126 g/L glucose; further suggesting that lincomycin A and lincomycin B were proximal to the region of maximum production. Accordingly, the concentrations of the three medium ingredients in the the third out of five experiments were considered the center point of BBD.
Box-Behnken design experimental optimization and response surface analysis of growth conditions that affect lincomycin production
BBD is a well-known optimization method based on the establishment of a mathematical model that assesses the statistical significance of the effects of different factors on the final response (Annadurai and Sheeja 1998). The BBD method was used to assess the effect of soybean powder, corn steep liquor, glucose concentrations on lincomycin A and B production. A total of 17 experiments with different combinations of nutrient concentrations were performed and their effects on lincomycin A and B production are presented in Table 4. Multiple regression analysis of the experimental data used a second-order polynomial equation derived for the lincomycin A yield and lincomycin B content using the significant terms:
Table 4
Box-Behnken Design experimental design
Run | X1 (g/L) | X2 (g/L) | X7 (g/L) | Lincomycin A (mg/L) | Lincomycin B (%) |
0 | 25 | 2 | 100 | 3519 ± 71 | 4.8 ± 0.2 |
1 | 25 | 1.68 | 151.2 | 4078 ± 39 | 3.0 ± 0.1 |
2 | 20 | 1.40 | 151.2 | 3464 ± 64 | 3.6 ± 0.1 |
3 | 25 | 1.68 | 100.8 | 4511 ± 36 | 3.5 ± 0.2 |
4 | 30 | 1.40 | 100.8 | 3752 ± 75 | 4.3 ± 0.2 |
5 | 20 | 1.12 | 126 | 3712 ± 61 | 2.8 ± 0.2 |
6 | 30 | 1.68 | 126 | 3452 ± 54 | 2.9 ± 0.1 |
7 | 20 | 1.68 | 126 | 4375 ± 122 | 2.6 ± 0.2 |
8 | 30 | 1.40 | 151.2 | 2558 ± 102 | 4.3 ± 0.2 |
9 | 25 | 1.12 | 100.8 | 4711 ± 58 | 3.9 ± 0.1 |
10 | 30 | 1.12 | 126 | 3734 ± 193 | 2.5 ± 0.2 |
11 | 25 | 1.12 | 151.2 | 3820 ± 19 | 2.7 ± 0.1 |
12 | 20 | 1.40 | 100.8 | 3712 ± 15 | 3.4 ± 0.1 |
13 | 25 | 1.40 | 126 | 4500 ± 41 | 1.0 ± 0.1 |
14 | 25 | 1.40 | 126 | 4580 ± 46 | 0.8 ± 0.1 |
15 | 25 | 1.40 | 126 | 4608 ± 78 | 1.3 ± 0.1 |
16 | 25 | 1.40 | 126 | 4612 ± 43 | 1.2 ± 0.1 |
17 | 25 | 1.40 | 126 | 4548 ± 57 | 0.8 ± 0.1 |
The values were obtained from three independent experiments. |
Y1 = 4569.6-220.88×X1 + 54.87×X2-345.75×X7-236.25×X1-236.50×X1X7 + 114.50×X2X7-829.93×X12 + 78.57×X22-368.17×X72 (2)
Y2 = 0.96 + 0.20×X1 + 0.013×X2-0.19×X7 + 0.15×X1X2-0.05×10− 3×X1X7 + 0.17×X2X7+1.18×X12 + 0.56×X22 + 1.76×X72 (3)
where Y1 (mg/L) is the predicted lincomycin A yield; Y2 (%) is the predicted lincomycin B content; and X1, X2, and X7 are the coded values of soybean powder, corn steep liquor, and glucose, respectively. Statistical significance of the response surface model and all the coefficient estimates were assessed with ANOVA (Tables S3 and S4). The high F-value (180.64) and a very low p-value (< 0.0001) suggested that the model was highly significant, while an insignificant lack of fit (p = 0.2201 > 0.05) revealed the effectiveness of the regression analysis, suggesting that the regression model could fit the effect of the three culture factors on lincomycin A content. The ratio of the explained and total variation indicated that the coefficient of determination (R2) could be used to assess the goodness of the model. The value of R2 was 99.57%, which indicated that only 0.43% of the lincomycin A content variability could not be explained by the predicted equation of the model. The Adj-R2 value of 99.02% further validated the significance of this model. A low coefficient of variation (CV) (CV = 1.44%) value revealed that the deviations between the predicted and experimental values were low, and it displayed a high degree of precision and reliability in the conducted experiments. "Adeq Precision" provides the signal-to-noise ratio and a ratio greater than four is desirable. In this study, a ratio (47.512) greater than 4 indicated the use of this model in future studies will be supported.
The ANOVA showed that the model for lincomycin B content had a p-value = 0.0006, suggested that the model was highly significant (Table S3). The lack of fit (p = 0.1252 > 0.05) was insignificant, suggesting that the regression model could fit the effect of the three culture factors on lincomycin B content.
The R2 was 0.9628, Adj-R2 was 0.9149, and the ratio was 14.541; these results indicated a high degree of precision and high reliability for this model, which supported its use in future studies. The content of lincomycin A and lincomycin B changed with changes in soybean powder, corn steep liquor and glucose concentrations, and their corresponding 3D response surfaces were generated to better determine the interaction of variables with the corresponding variables (Fig. 1). The maximum lincomycin A and lincomycin B contents were 4569.6 mg/L and 0.96%, respectively, with 25 g/L soybean powder (X1), 1.40 g/L corn steep liquor (X2), and 126 g/L glucose (X7) .
Fermentation was conducted and lincomycin production in broth was evaluated by HPLC (Fig. 2). The results for the validation experiment showed that the values for the two responses of lincomycin A and lincomycin B contents were in close agreement with the predicted values. When 25 g/L soybean powder, 1.40 g/L corn steep liquor, and 126 g/L glucose were used, the content of lincomycin A increased from 3585 ± 110 mg/L to 4600 ± 134 mg/L, and that of lincomycin B decreased from 4.5%±0.2% to 0.8 ± 0.1%. These concentrations verified the accuracy of the statistical model.
Assay of lincomycin A, lincomycin B, and fermentation parameters
Glucose is not only a carbon source but can also affect the osmotic stress of a culture medium at certain concentrations, while corn steep liquor provides various amino acids, vitamins, and metal ions, etc. To explore possible causes of the increased lincomycin A and decreased lincomycin B after medium optimization, the fermentation parameters for S. lincolnensis 24 were measured (Fig. 3) using samples collected every 24 h.
The optimized medium not only influenced the nutrient content of the culture, such as sugar (Fig. 3A), it also affected the growth environment, such as osmotic stress and pH parameters (Fig. 3A and B). As shown in Fig. 3A and 3B, reducing the sugar and osmotic stress in the optimized medium were always higher than the initial medium before 144 h. At the end of 168 h of fermentation, the reducing sugars in and osmotic stress of the optimized medium were the same as the initial medium, with the reducing sugar consumed completely. These results indicated that changes in glucose affected the reducing sugar and osmotic stress of the medium. We hypothesized that these variations in the environment and nutrients affected the primary metabolism of S. lincolnensis and caused PMV changes (Fig. 3C), as well as changes in the secondary metabolism of S. lincolnensis, which improved lincomycin A quantity and quality (Fig. 3D). As shown in Fig. 3D, the optimized medium produced lincomycin A before 48 h, while the initial medium produced lincomycin A after 48 h. For lincomycin B, only the initial medium produced lincomycin B before 48 h. For both the initial and optimized media, both quickly increased lincomycin A and decreased lincomycin B in the 72–96 h stage.
Influence of corn steep liquor, glucose, and osmotic stress on the optimized medium conditions
In the present study, the increased glucose concentration affected the osmotic stress environment of the fermentation, which benefited S. lincolnensis growth and the reduction of lincomycin B accumulation (Fig. 3). We hypothesized that the change in osmotic stress had a significant influence on cell growth and lincomycin synthesis. Therefore, to explore the role of osmotic stress in lincomycin fermentation, NaCl was used to regulate osmotic stress in the initial medium. To determine the effect of salt-stress sensitivity on lincomycin production, S. lincolnensis was grown in the initial medium (control, containing 5 g/L NaCl), optimized medium (OP, initial medium with 1.4 g/L corn steep liquor and 126 g/L glucose), initial medium with 10 g/L NaCl (M1), initial medium with 126 g/L glucose (M2), or initial medium with 1.4 g/L corn steep (M3). The osmotic pressure of the OP and M1-M2 media was similar. As shown in Fig. 4, the OP medium produced the highest lincomycin A and lowest lincomycin B concentrations. Increasing the salt stress with 10 g/L NaCl in the initial medium (M1) resulted in a dramatic decrease in lincomycin B production. Similarly, under glucose (M2) and corn steep liquor (M3) salt-stress conditions, a knockdown of about 31.8–63.6% was achieved for lincomycin B production when compared with the control medium. These observations indicated that the osmotic conditions decreased lincomycin B concentrations during shake flask fermentation.
Transcriptional analysis of lincomycin biosynthesis-related genes
To further understand the observed variations between lincomycin production in initial and optimized media, quantitative real-time PCR (qRT-PCR) gene expression analysis was performed. In the biosynthesis of lincomycin, lmbW is an important methylase gene that is responsible for the C-methylation of 4-n-propyl-L-proline and determines whether 4- n-propyl-L-proline synthesizes lincomycin A or lincomycin B (Pang et al. 2015). The regulatory genes in lincomycin synthesis are divided into the lincomycin biosynthesis gene cluster (lmb cluster) and non-lmb genes, which are lmbU and SLCG_Lrp, respectively. These genes play positive regulatory roles in lincomycin biosynthesis (Hou et al. 2018; Xu et al. 2020). The lmbB1 gene encodes for a 2,3-extradiol cleavage enzyme that breaks down the L-3,4-dihydroxyphenylalanine (L-DOPA) aromatic ring (Novotná et al. 2004). Replacing the tyrosinase involved in melanin synthesis, LmbB1 is assumed to be a dioxygenase that catalyzes the cleavage of 2,3-extradiol in the aromatic ring of L-3,4-dihydroxyphenyl alanine (Zhong et al. 2017). LmbB1 is a key gene that determines whether L-DOPA follows the lincomycin or the melanin biosynthesis pathway. Overexpression of lmbB1 can increase lincomycin A and decrease lincomycin B and melanin (Yang et al. 2020). These genes were selected based on their roles in the synthesis of lincomycin A and lincomycin B. Furthermore, osmotic stress in the optimized medium was higher than that of the initial medium (Fig. 3A); therefore, the transcript level of the osmotic stress related gene (mscl; GenBank: ANS65846.1) was also determined.
RNA was isolated from both initial and optimized media samples collected after 48 h of fermentation. qRT-PCR was conducted in triplicate for each sample, using 16S rRNA as an internal control (RT primers are listed in Table 5). The qRT-PCR data showed that lmbU, SLCG_Lrp, lmbW and lmbB1 expression was significantly higher in the optimized medium than in the initial medium (2- to 4-fold increase; Fig. 5). One study recently reported that LmbU is a significant pleiotropic transcriptional regulator in lincomycin biosynthesis that activates the lmb cluster, including lmbW and lmbB1 in S. lincolnensis (Lin et al. 2020). Another study showed that SLCG_Lrp is a positive regulator for lincomycin biosynthesis through the direct induction of lmb cluster genes, such as lmbU (Xu et al. 2020). Our results indicated that increased lincomycin A production in the optimized medium was likely to be associated with the overexpression of the key genes involved in the lincomycin biosynthesis pathway, as well as some of the key regulators of the pathway.
Table 5
name | primer sequences 5’→3’ |
SLCG_Lrp | (F)-TCGTCGTACAGCCGCTGGTAG (R)-GATCGCGGAAGTGGTGGATGC |
lmbU | (F)-GCGTAGTTGCGGATCGTCTGG (R)-ACTCATCGGCCTGGTGTCTGG |
mscl | (F)-CCTCATGGTCCTGCCGATGT (R)-AGCTCGCTCACCTCGATGAC |
lmbW | (F)-A G C T G C T G G C C G A G G G C G T A (R)-G C C G C C G G A C T T G G A C G A C A |
lmbB1 | (F)-AGTAAAGTCAATGCCGCCCGTATC (R)-GAATGTGTCGAGGGTCCAGGAAAC |
16S rRNA | (F)-GCATCTGTGGTGGTTGAAAG (R)-CGTGTCTCAGTCCCAGTGTG |
Microorganisms cope with environmental stress factors during industrial fermentation processes, such as osmotic, temperature, and oxidative stress stresses, which can significantly impact their primary and secondary metabolism (Lee et al. 2005; Li et al. 2009). Osmotic stress has an important effect on the secondary metabolism of Streptomyces. For example, disruption of either osaBSa, which encodes a response regulator protein, led to increased production of oligomycin up to 200%, and avermectin up to 37% (Godinez et al. 2015). Regulation of the osmotic stress response can impact both development and antibiotic production in the model streptomycete, S. coelicolor (Bishop et al. 2004; Martínez et al. 2010). Bacteria have two families of mechanosensitive (MS) channels: small conductance (MscS) and large conductance (MscL) channels (Blount and Moe 2005; Blount and Iscla 2012). The majority of bacteria contain a single copy of mscL, which is highly conserved between species (Wray et al. 2019). In bacteria, MS channels act as emergency release valves; when bacteria are exposed to high osmotic pressure, they will transport (K+, glutamate, betaine, and proline) and synthesize (glutamate, trehalose, proline, and betaine) solutes to balance the increase in external osmotic pressure to maintain high cell turgor, which is a requisite for cell growth and division.
In this study, the optimized medium enhanced the osmotic stress of the fermentation condition and led to a 6.2-fold increase in the transcription of mscl when compared with the initial medium (Fig. 5). The expression levels of lmbU were 3.9-fold up-regulated in the optimized medium after 48 h compared to those of initial medium (Fig. 5). Recent studies have suggested LmbU can bind to the regions upstream of the lmbA and lmbW genes at the consensus and palindromic sequence 5’-CGCCGGCG-3’ (Hou et al. 2018). By sequence alignment, this palindromic sequence was located 184 bp upstream of the mscl gene (Fig. 6). LmbU binds directly to the regions upstream of lmbW and activates transcription. Thus, LmbU may also actively inhibit the production of lincomycin B (Hou et al. 2018; Pang et al. 2015). Higher osmotic stress might result in the down-regulation of lincomycin B biosynthesis (Fig. 4) due to LmbU targeting mscl and affecting the expression of lmbW during lincomycin B biosynthesis.
As shown in Fig. 3D, under the optimized conditions, the emergence time of lincomycin A was enhanced, which effectively reduced the fermentation period. Previous studies showed that the pneumocandin B0 fermentation period was significantly shorter when Glarea lozoyensis were cultured in a high osmolarity medium (Song et al. 2018); thus, osmotic stress could shorten the fermentation period in other bacteria. Moreover, simultaneously with the enhanced production time of lincomycin A, that of lincomycin B was delayed by approximately 24 h (Fig. 3D). These results suggest a temporal order and interconversion of lincomycin A and lincomycin B biosynthesis that are related to osmotic stress. The mechanism stated for lincomycin A and lincomycin B interconversion has not been reported. This study provides a foundation for interconversion of lincomycin A and lincomycin B in S. lincolnensis, although further research is needed to develop its application potential.