Integrating multi-omics analyses of Nonomuraea dietziae to reveal the role of soybean oil in [(4′-OH)MeLeu]4-CsA overproduction

Background Nonomuraea dietziae is a promising microorganism to mediate the region-specific monooxygenation reaction of cyclosporine A (CsA). The main product [(4′-OH)MeLeu]4-CsA possesses high anti-HIV/HCV and hair growth-stimulating activities while avoiding the immunosuppressive effect of CsA. However, the low conversion efficiency restricts the clinical application. In this study, the production of [(4′-OH)MeLeu]4-CsA was greatly improved by 55.6% from 182.8 to 284.4 mg/L when supplementing soybean oil into the production medium, which represented the highest production of [(4′-OH)MeLeu]4-CsA so far. Results To investigate the effect of soybean oil on CsA conversion, some other plant oils (corn oil and peanut oil) and the major hydrolysates of soybean oil were fed into the production medium, respectively. The results demonstrated that the plant oils, rather than the hydrolysates, could significantly improve the [(4′-OH)MeLeu]4-CsA production, suggesting that soybean oil might not play its role in the lipid metabolic pathway. To further unveil the mechanism of [(4′-OH)MeLeu]4-CsA overproduction under the soybean oil condition, a proteomic analysis based on the two-dimensional gel electrophoresis coupled with MALDI TOF/TOF mass spectrometry was implemented. The results showed that central carbon metabolism, genetic information processing and energy metabolism were significantly up-regulated under the soybean oil condition. Moreover, the gas chromatography-mass spectrometry-based metabolomic analysis indicated that soybean oil had a great effect on amino acid metabolism and tricarboxylic acid cycle. In addition, the transcription levels of cytochrome P450 hydroxylase (CYP) genes for CsA conversion were determined by RT-qPCR and the results showed that most of the CYP genes were up-regulated under the soybean oil condition. Conclusions These findings indicate that soybean oil could strengthen the primary metabolism and the CYP system to enhance the mycelium growth and the monooxygenation reaction, respectively, and it will be a guidance for the further metabolic engineering of this strain. Electronic supplementary material The online version of this article (doi:10.1186/s12934-017-0739-0) contains supplementary material, which is available to authorized users.


Background
Cyclosporine A (CsA), produced by soil fungus Tolypocladium niveum, is a natural cyclic undecapeptide possessing immunosuppressive activity, and is the active ingredient of Sandimmune ® and Neoral ® for preventing organ transplant rejection [1]. Besides, CsA is reported to substantially inhibit the virus replication of HCV (hepatitis C virus) [2,3] and HIV-1 (human immunodeficiency virus type 1) [4]. However, developing a CsA-based anti-HIV/HCV drug should exclude the immunosuppressive activity of CsA since it will antagonize the host immune
Nonomuraea dietziae (Sebekia benihana), a rare soil actinomycete, is one of the best microorganisms that can efficiently convert CsA to the region-specific hydroxylation product [(4′-OH)MeLeu] 4 -CsA [13] by cytochrome P450 hydroxylases (CYPs) [14]. To improve the bioconversion rate of CsA, various strategies have been applied so far, such as genetic manipulation [14,15], medium optimization and traditional mutation. Lee et al. overexpressed CYP-sb21 gene in wild-type N. dietziae and increased the conversion rate of CsA by twofold (reaching 29%) [16]. By the same token, a 54% conversion rate was obtained when molybdenum salt was added into the optimized medium [17]. However, the product titers and yields of these strains are still at a low level, and impede the industrial application of [(4′-OH)MeLeu] 4 -CsA.
Recently, we have successfully obtained a mutant strain with the high [(4′-OH)MeLeu] 4 -CsA production by UV-LiCl complex mutation. Interestingly, addition of 0.1% (w/v) soybean oil could further increase the conversion rate by 55.6%. Although soybean oil has been proved to accelerate the strain's growth and improve the antibiotic production [18][19][20], the specific effects on [(4′-OH) MeLeu] 4 -CsA production may be complex and comprehensive because the conversion process from CsA to [(4′-OH)MeLeu] 4 -CsA is a monooxygenation reaction and seems irrelevant to soybean oil. Hence, to unveil the potential mechanism, in this study, the two-dimensional gel electrophoresis (2-DE) coupled with matrix-assisted laser-desorption/ionization time-of-flight/time-of-flight mass spectrometry (MALDI-TOF/TOF-MS) and gas chromatography-mass spectrometry (GC-MS), were employed for proteomic and metabolomic analyses, respectively. Moreover, the transcription levels of CYP genes were analyzed by quantitative real-time PCR (qRT-PCR) to investigate the role of the soybean oil on the overproduction of [(4′-OH)MeLeu] 4 -CsA.

Effects of the exogenous soybean oil on the fermentation properties of N. dietziae
As shown in Fig. 1, the major fermentation features of N. dietziae in soybean oil medium (assigned as the medium MO) were distinctly different from the control group (medium MC). The initial 12 h was the lag phase and the both biomass changed little, while the following log phase (12-48 h) displayed significant discrepancies (Fig. 1a). During the log phase, biomass in the medium MO increased sharply from 3.63 to 6.4 g/L (improved by 76.3%, Fig. 1a), compared with a relatively gentle change in the medium MC (improved by 43.5%, Fig. 1a), indicating a more appropriate growth environment under the MO condition. Throughout the whole fermentation process, the highest biomass (7.58 g/L) and [(4′-OH)MeLeu] 4 -CsA titer (233.4 μmol/L or 284.4 mg/L) were achieved at 120 h in the medium MO, 1.16-and 1.56-folds compared to that in the medium MC, respectively ( Fig. 1a, b). After 120 h, the [(4′-OH)MeLeu] 4 -CsA production declined gradually due to the product degradation.
In this study, the intrinsic conversion capacity of N. dietziae per unit of biomass was characterized by the specific production rate (μ p ), an important kinetic parameter that removed the impact of biomass on the [(4′-OH) MeLeu] 4 -CsA production. As shown in Fig. 1c, the maximum μ p was achieved in the medium MO (0.67 μmol/g/h, 36 h), which was 1.3-folds and appeared 12 h earlier than that in the MC medium (0.51 μmol g/h, 48 h). Although the μ p of MO was lower than MC during the late fermentation period (from 84 h to 120 h, Fig. 1c), yet the most production of [(4′-OH)MeLeu] 4 -CsA achieved at 84 h in both media (Fig. 1b), indicating that the catalytic capacity of medium MO was superior to MC.
In addition, the concept "total mole concentration" (TMC) was put forward to describe the total mole concentration of the residual CsA and the produced [(4′-OH)MeLeu] 4 -CsA in the fermentation broth media based on the fact that the hydroxylation process is a monooxygenation reaction, in which the conversion coefficient of CsA to [(4′-OH)MeLeu] 4 -CsA is 1 mol:1 mol. In Fig. 1d, the TMC reached the minimum value at 48 h, but significantly increased during the late fermentation period in both medium. This phenomenon indicated that part of CsA was probably converted into other intermediates, such as the CsA-CYPs complex, as mentioned previously [16,21]. Accordingly, the lower value of TMC in the MO medium represented a stronger intermediate processing capacity (or catalytic ability) than that in the control.

Effects of the major hydrolysates of soybean oil and some other plant oils (corn oil and peanut oil) on [(4′-OH) MeLeu] 4 -CsA production
Since soybean oil is a natural lipid mixture and can be hydrolyzed into fatty acids and glycerol, both of which will subsequently participate in the metabolic system [19,22]. Considering the potential association between the lipid metabolism and the soybean oil, series of feeding experiments were designed and carried out, including the main hydrolysates of soybean oil, such as oleic acid, linoleic acid and glycerol, as well as the soybean oil-like compounds, such as corn oil (0.1%, w/v), peanut oil (0.1%, w/v) and Tween 80 (polyoxyethylene sorbitan monooleate, 0.05%, w/v). All the compounds were fed into the medium MC at the beginning of fermentation.
As shown in Fig. 2, oleic acid, linoleic acid and Tween 80 all failed to enhance the conversion of CsA to [(4′-OH) MeLeu] 4 -CsA. Glycerol was an exception since the [(4′-OH)MeLeu] 4 -CsA production had a 7% increase. Considering that glycerol was probably used as a complementary carbon source, and that increasing the concentration of glycerol could not further enhance the CsA conversion (data not shown), a 7% increase, nevertheless, was still much smaller compared with the impact of soybean oil (55.6%). Thus, it could be concluded that the main hydrolysates of soybean oil contributed little or even negatively to the improvement of CsA conversion.
Addition of corn oil and peanut oil, by contrast, led to the 20.7 and 25.8% increases in [(4′-OH)MeLeu] 4 -CsA production, respectively, suggesting that plant oils (soybean oil, corn oil and peanut oil) could efficiently improve the conversion of CsA. Therefore, although N. dietziae cannot directly utilize the hydrolysates of plant oil to strengthen mycelium growth and [(4′-OH)MeLeu] 4 -CsA  4 -CsA indicated that plant oils worked in another way rather than being the substrates. Taken the above results, soybean oil played a significant role in improving biomass and specific production rate. However, soybean oil might not play its role in the lipid metabolic pathway and more intercellular metabolic details still remained to be investigated. To this end, a combined proteomic and metabolic analysis was implemented to further reveal the role of soybean oil in [(4′-OH)MeLeu] 4 -CsA overproduction.

Comparative proteomic and metabolomic analyses in response to the soybean oil addition Comparative proteomic profile analysis
Protein extracts from both conditions were sampled at 48 h and 96 h and then subjected to 2-DE. The results of gel electrophoresis were presented in Additional file 1: Figure S4 and 2-DE could effectively separate most of the proteins, indicating the feasibility of our experimental methods. In this study, a total of 95 protein spots in the 2-DE gels with significantly differential expressions (fold change >1.5) were identified (Table 1). They were classified into eight functional groups by their cellular roles, mainly including central carbon metabolism, energy metabolism, genetic information processing, amino acid metabolism, regulatory proteins, nucleotide metabolism, proteins of unknown function and hypothetical proteins. Figure 3 presented the protein distribution at the sampling time with different abundances, and most of the proteins in each group were present at higher levels under the MO condition, suggesting a comprehensive strengthening effect of soybean oil on the whole cell metabolism.

Comparative metabolic profile analysis
Intercellular metabolites of N. dietziae (48, 72, 96 and 120 h) in the medium MO and MC with different [(4′-OH)MeLeu] 4 -CsA producing capabilities were analyzed by GC-MS. As a result, a total of 50 intracellular metabolites were identified, including amino acids, organic acids, sugars and fatty acids. Abundance changes of metabolites are depicted by the heat map in Fig. 4.
The metabolomic profiling showed obvious differences in tricarboxylic acid (TCA) cycle and amino acid metabolism. More specifically, succinic acid increased to 2.8-and 5.5-fold of that in the medium MC at 48 and 96 h, respectively. Similar change patterns also exhibited in fumarate (1.4-and 1.2-fold vs. the control) and α-ketoglutaric acid (1.2-and 2.2-fold vs. the control). Although the abundances of these metabolites decreased gradually along the fermentation process, they were still much higher than the control, indicating the enhancement of TCA cycle in medium MO. On the contrary, the abundance of amino acids such as threonine, glutamate, lysine, valine and leucine were relatively lower in medium MO. Particularly, both glutamate and lysine decreased by 78% at 72 h. A reasonable explanation was that large amount of amino acids was effectively consumed to synthesize cellular building blocks for cell growth.
Based on the proteomic and metabolomic results, we therefore summarized the metabolic profile in Fig. 5, which visually described the intracellular metabolic pathway in response to the exogenous soybean oil. In this model, most proteins with high levels (in red font) were mainly involved in the central carbon metabolism, amino acid metabolism, genetic information processing and oxidation-reduction process. As regarding to the hydroxylation of CsA to [(4′-OH)MeLeu] 4 -CsA, it was directly mediated by CYPs and P450 oxidoreductase (POR) along with the formation of water and NAD(P) + , which was seemingly uncorrelated to other intracellular metabolic activities. Subsequently, pathway analysis was implemented to dissect the effects of soybean oil on the metabolism and regulation of N. dietziae in detail.

Central carbon metabolism and amino acid metabolism
Central carbon metabolism represents the backbone of the cellular metabolism and provides the precursors required for the cell growth and the synthesis of target products. As shown in Fig. 5, it was worth noting that some key enzymes involved in glycolytic pathway, such as 6-phosphofructokinase (Spot 6, PFK) and pyruvate kinase (Spot 8, PYK) were both present at higher levels.   Different from many other actinomycetes, PFK of Nonomuraea is not allosterically regulated by ATP, AMP, ADP, or phosphoenolpyruvate and pyruvate, but controlled by the availability of pyrophosphate (PPi) produced in nucleic acid and protein biosynthesis and in the cycling between glycogen and glucose-1-phosphate [23]. Additionally, the activity of PPi-dependent PFK is reversible, suggesting that N. dietziae had a flexible glycolytic pathway as the regulatory node. Apart from PFK, dihydrolipoamide dehydrogenase (Spot 5, DLD), pyruvate dehydrogenase (Spot 7, PDH) and enolase (Spot 3, ENO) also showed higher levels under the MO condition, suggesting a strengthened Embden-Meyerhof-Parnas (EMP) pathway. In regard to TCA cycle, citrate synthase (Spot 11, CS), which converts acetyl-CoA and oxaloacetate to citrate in the initial step and controls flux into the TCA cycle [24], showed a higher level (3.23-folds) under the MO condition while the increase in the level of malate dehydrogenase (Spot 9, Mdh) could supply more oxaloacetate that reacts with acetyl-CoA [25]. In addition, the higher levels of alpha-ketoglutarate decarboxylase (Spot 13, SucA) and dihydrolipoyl dehydrogenase (Spot 51, DLDH) further improved the metabolic rate of TCA cycle [26]. As a consequence, TCA-related metabolites, such as fumarate, α-ketoglutaric acid and succinic acid, were all present at higher levels (Fig. 4), indicating that the TCA cycle was active under the MO condition from both proteomic and metabolomic insights.
Additionally, glucose-6-phosphate dehydrogenase (Spot 1, G6PD), the first key enzyme of pentose phosphate (PP) pathway, had a higher level (1.67-fold) at 48 h under the MO oil condition, but at 96 h the abundance of G6PD decreased to 56% of the control. G6PD-catalyzed reaction dictates and limits the fluxes between Embden-Meyerhof-Parnas (EMP) and PP pathways [27], and the PP pathway oxidizes glucose to generate NADPH for reductive biosynthesis reactions and ribose-5-phosphate for the synthesis of the nucleotides [28]. The decreased level of G6PD under the MO condition restricted the specific cell growth (Fig. 1, 96 h) and the specific production rate of [(4′-OH)MeLeu] 4 -CsA (Fig. 1, 96 h) since the cofactor NADPH also participated in the CsA hydroxylation process in addition to biomass accumulation.
Another significantly changed protein, phosphogluconate dehydratase (Spot 2, Edd), was only detected at 96 h and had a 64% increase after the soybean oil addition. Gunnarsson et al. had identified the ED pathway in Nonomuraea sp. by the 13 C labelling-based method [29]. ED pathway catabolizes glucose to pyruvate by using the 6-phosphogluconate dehydratase (Edd) and 2-keto-3-deoxyphosphogluconate aldolase (KDPG aldolase), which connects the EMP, PP pathways and TCA cycle, and tunes the carbon flux distribution, the production of ATP, reducing equivalent for cell growth. As the first critical enzyme of ED pathway, the higher level of Edd under the MO condition probably indicated a more robust sugar metabolism at the end of fermentation than the control.

Fig. 3
Protein distribution with significant changes in each functional category at 48 and 96 h respectively. 2-DE profiles of N. dietziae were presented in Additional file 1: Figure S4 and the significantly differential proteins and their characteristics were listed in Table 1 Amino acid metabolism is essential for cell growth by supplying the building blocks and metabolic intermediates along with EMP, TCA and PPP. In this study, three proteins involved in glutamate synthesis, i.e. Spot 42 (urocanate hydratase, HutU), Spot 43 (histidine ammonia-lyase, HAL) and Spot 35 (glutamate synthase, GS), as well as some other amino acid metabolism-related proteins (Spot 16, 2-isopropylmalate synthase, LeuA, Spot 40, cysteine synthase, CysO) had significantly higher levels in the medium MO. Although abundances of glutamate, leucine and alanine were much lower in the medium MO than that in the medium MC (Fig. 4), it could probably be resulted from the effective utilization for protein synthesis and cell growth. In addition, two proteins in S-adenosyl-l-methionine (SAM) cycle, i.e. Spot 33 (SAM synthetase, MetK) and Spot 34 (S-adenosyl-l-homocysteine hydrolase, SahH) were detected at higher abundance in the medium MO, which were involved in the synthesis of Fe-S protein clusters and methionine [30].
Therefore, it can be concluded that addition of soybean oil could significantly upregulate the intracellular central carbon metabolism and amino acid metabolism, which supplies a more suitable intracellular environment for maintaining the cell robustness and promoting the [(4′-OH)MeLeu] 4 -CsA overproduction.

Genetic information processing and nucleotide metabolism
Bacterial growth is directly correlated to the synthesis of protein and DNA. In this study, the higher abundances of phenylalanine-tRNA ligase (PheRS, Spot 36), aspartyl/ glutamyl-tRNA amidotransferase subunit B (GatB, Spot 37), proline-tRNA ligase (ProRS, Spot 38) and tryptophanyl-tRNA synthetase (TrpRS, Spot 39) were observed under the MO condition, indicating that aminoacyl-tRNA biosynthesis was activated. Additionally, elongation factor Ts (EF-Ts, Spot 46) and 30S ribosomal protein S2 (Spot 49) were also present at higher levels compared with the control. EF-Ts mediates the regeneration of EF-Tu-GDP complex, which catalyzes the addition of aminoacyl-tRNA into ribosome [31,32]. Meanwhile, DNA-directed RNA polymerase subunit alpha (Spot 47), DNA polymerase III subunit beta (Spot 48), transcription termination factor NusA (Spot 28) and transcription termination factor Rho (Spot 50), were present at higher levels. They are involved in the transcription and replication processes and inosine-5-monophosphate dehydrogenase (Spot 44) provides precursors (guanine) for RNA and DNA synthesis [33]. Higher levels of the above proteins could efficiently enhance the cell growth, leading to a higher specific growth rate under the MO condition as shown in Fig. 1.

Cytochrome P450 (CYP) and energy metabolism
CYPs belong to a family of terminal monooxygenases that transfer one oxygen atom to X-H bonds (X:-C,-N, S) of a substrate with the concomitant reduction of the other oxygen atom to water [14]. CYPs directly participate in the bioconversion of CsA to [(4′-OH)MeLeu] 4 -CsA. In this study, five CYP isoforms ( Spot 15,[17][18][19][20] were identified and most of them showed higher abundance compared with the control (Table 1). Proteins involved in energy metabolism, especially the redox reaction, displayed a significant improvement (Table 1) under the MO condition. 15 of 16 proteins involved in energy metabolism had higher levels at 48 h, and 10 proteins demonstrated the similar expression pattern at 96 h (Table 1). In addition, proteins related to the electron transfer chain showed significantly higher levels in the medium MO. In particular, NADH dehydrogenase (Spot 21), the first proton pump in oxidative phosphorylation, showed a 2.4-fold higher abundance than the control at 48 h. Throughout the fermentation process, the abundances of NADH dehydrogenase (Spot 21), FADlinked oxidase (Spot 26) and flavoprotein oxidoreductase (Spot 27) were significantly increased (>3.1-folds) in the medium MO. It is worth noting that NADH dehydrogenase showed the maximum change (20.9-folds) among all the identified proteins. Since the growth of Nonomuraea sp. is strictly dependent on the aerobic metabolism and oxidative phosphorylation [34], as the key members of electron transfer chain, these proteins with higher levels would enhance the accessibility and the turnover rate of NADH and FADH to CYPs so as to improve the hydroxylation activity [15].
The sampling time points 36 h (12 h after soybean oil addition) and 48 h (24 h after soybean oil addition) were selected since the greatest specific production rate appeared before 48 h under both conditions. As shown in Table 2, only 11 CYP genes were detected under both conditions while the other 10 genes were detected only in certain conditions. At 36 h, not only the number of the genes but also the transcription abundance of the CYPs in the soybean oil group were significantly increased compared with the control group (Table 2), which could explain why the specific production rate of the soybean oil group was much higher than that of the control group at 36 h.
At 48 h, the number of the detected CYPs in the control group was larger than that of the soybean oil group, but the transcription levels of the crucial CYP gene sb21, sb13 and sb3-2 were still much lower than the soybean oil group, which could explain the phenomenon that in Fig. 1c the specific production rate of the control group at 48 h was increased, but was still lower than that in the MO condition. In addition, the upregulation of the transcription level of the CYPs under the MO condition also implied the formation of the CYP-CsA complex since the TMC (total mole concentration) in the fermentation broth reached bottom at 48 h and the TMC in the MO medium was lower than that in the MC medium as shown in Fig. 1d. Interestingly, the transcription level of CYP-sb21 in the control group greatly reduced at 48 h ( Table 2) but the catalytic efficiency reached the maximum (Fig. 1c). Although sb21 played a leading role in the conversion of CsA to [(4′-OH)MeLeu] 4 -CsA [14], yet the role of other CYPs should not be neglected because the expression of sb3-1, sb4, sb7, sb10, sb16, sb17, sb23 and sb24, and the Table 2   upregulation of sb3-2, sb8, sb9, sb11 and sb15 counteracted the effect of the down-regulated sb21.

Proposed metabolic mechanism of [(4′-OH)MeLeu] 4 -CsA overproduction under the MO condition by N. dietziae
Soybean oil could improve the production of antibiotics, such as tetracycline, cephamycin C and tacrolimus [18][19][20], but the specific mechanism was not always the same. Soybean oil was found to enhance tetracycline production by extracting antibiotic into the oil phase of the culture broth, thereby relieving product inhibition and decreasing damage to microbial cells by foam formation [18]. While in cephamycin C production, soybean oil mainly functioned as a carbon source for Streptomyces sp. p6621 fermentation [19]. Additionally, soybean oil could induce the expression of lipase to produce CoAesters, the precursors of FK506 by S. tsukubaensis [20,22]. These three products have a similarity that they are produced by the de novo synthesis pathway. However, [(4′-OH)MeLeu] 4 -CsA is converted from CsA by the monooxygenation reaction and is directly associated with CYPs, POR (cytochrome P450 reductase), oxygen and reducing equivalent. Thus, the role of soybean oil herein is probably not one of the above mentioned. In this study, the apparent catalytic efficiency (or conversion rate) is closely associated with two factors, i.e., the conversion capacity per unit of the biomass (μ p ) and the quantity of biomass (DCW). μ p is the intrinsic property which is dependent on the abundance of CYPs, PORs and the accessibility of reducing equivalent supplied by the cell metabolism. The synthesis of biomass relies on the abundance of intracellular building blocks, the rate of genetic information processing and the energy availability (ATP and reducing equivalent). On one hand, the proteomic analysis showed that the significantly changed proteins were involved in EMP, TCA cycle, amino acid metabolism and redox process, thereby enhancing the flux into central carbon metabolism and supplying sufficient cellular building blocks, ATP and reducing equivalents for maintaining the cell robustness and promoting the hydroxylation of CsA. Meanwhile, the improvement of transcription and translation process helps to accelerate the protein synthesis and mycelium growth. Moreover, metabolomic analysis indicated that soybean oil had a great effect on amino acid metabolism and tricarboxylic acid cycle.
On the other hand, the transcriptional analysis of all the CYPs under both conditions confirm that soybean oil can strengthen the CYP system for the conversion of CsA to [(4′-OH)MeLeu] 4 -CsA. Additionally, the enhanced hydroxylation ability is not only dependent on the elevated expression of CYP-sb21, but also on other CYPs, such as CYP-sb13 and CYP-sb8, although they work with different degrees. These findings demonstrate that the CsA conversion process is under a sophisticated and systematic regulation, although it's a simple in vivo monooxygenation reaction.
Another phenomenon is that besides the soybean oil, some other plant oils, such as corn oil and peanut oil, could also exert a positive influence on the [(4′-OH) MeLeu] 4 -CsA improvement, implying a common feature of these plant oils. In the proteomic analysis, it was noteworthy that no differentially expressed proteins were observed in lipid metabolism ( Table 1), implying that the major contribution of soybean oil may not be related to lipid metabolism. As regarding to the metabolomic data, oleic acid, stearic acid, and myristic acid were more abundant in MC samples, especially in the initial fermentation stage (48 and 72 h) (Fig. 4). These results demonstrated that lipid metabolism had not been activated under the MO condition. Moreover, the main hydrolysates of soybean oil did not strengthen the CsA conversion ability (Fig. 2). All these results indicated that plant oils could influence the fermentation characteristics. A reasonable explanation was that plant oils increased the oxygen transfer efficiency due to their lower polarity and the stronger oxygen-carrying capacity compared with water [35], thus strengthening the oxygen supply for intracellular redox metabolism and CsA conversion process, just as shown in metabolic pathway analysis in response to soybean oil.

Conclusions
In summary, a systematic comparative proteomic and metabolomic analysis was successfully conducted to gain insights into the role of soybean oil in improving [(4′-OH)MeLeu] 4 -CsA production by N. dietziae. Soybean oil could strengthen the primary pathways and the CYP system of N. dietziae, thereby improving the rate of biomass synthesis and the hydroxylation efficiency. The omics-based analysis provides a number of intracellular biomarkers, and is a starting point for exploring the regulatory mechanism between cell growth and [(4′-OH) MeLeu] 4 -CsA production, which will be a guidance for the further metabolic engineering of this strain.

Microorganism and cultivations
Nonomuraea dietziae used in this study was stocked in our laboratory and cultivated on ISP-2 agar slant [36]. Seed medium was prepared as Zhang's work [37]. Production medium: 20 g/L glucose, 3 g/L yeast extract, 10 g/L peptone, 12 g/L dextrin, 15 g/L corn steep liquor and pH 6.5. N. dietziae spores were washed from the fresh agar slant by 5 mL sterile 0.9% NaCl solution and transferred into 100 mL seed medium of a 500 mL Erlenmeyer flask and then incubated at 28 °C, 220 rpm for 72 h. The [(4′-OH)MeLeu] 4 -CsA production was carried out in 50 mL medium of a 250 mL flask at 220 rpm for 120 h at 30 °C after inoculating 10 mL seed culture. 0.1% (w/v) soybean oil was fed at the beginning of cultivation in the experimental group. Here, CsA was predissolved in 95% ethanol solution and added into the cultivation medium at 24 h with an initial concentration of 600 mg/L (499.2 μmol/L).

Analytical methods
The biomass yield was determined by dry cell weight (DCW). For the determination of the concentration of CsA and [(4′-OH)MeLeu] 4 -CsA, 5 mL culture fluid was immediately mixed with 5 mL 95% ethanol and shaken intermittently for 1 h. After centrifugation, the supernatant was subjected to HPLC (Agilent 1200, USA) equipped with an Eclipse XDB-C18 column (5 μm; 150 mm × 4.6 mm; Agilent Technologies) and a UV detector at 210 nm. The mobile phase was acetonitrile-0.1% phosphoric acid water solution (70:30, v/v) with a flow rate of 1.0 mL/min, and the column temperature was 60 °C.
The specific production rate (μp) is defined as μp = dp/ dt/M, and the specific consumption rate is defined as μc = dc/dt/M, where p is the [(4′-OH)MeLeu] 4 -CsA production and c is the residual CsA in the fermentation broth, t is fermentation time and M is the dry weight of mycelium. dp/dt and dc/dt was calculated by the slope calculating plug-in "Tangent.opk" in Origin 8.1 (Origin-Lab, USA).
The staining gels were scanned at 300 dpi resolution by Umax Powerlook 2100XL Flatbed Scanner (UMAX Technologies Inc., Dallas, TX, USA) [41]. Subsequently, the image was analyzed with the Bio-Rad PDQuest Basic 2-D image processing software (version 8.0.1). By using PDQuest software, the average ratios and t test values for each spot and a value below 0.05 for t test was regarded as being significant. Protein spots with an average abundance change of greater than 1.5-fold and present in all biological replicates were subjected to MS analysis.
Protein spots were detained and digested as previously described [41]. More details were present in Additional file 1 (Protein extraction and proteomics analysis). The digested peptides were analyzed using a 4700 Proteomics Analyzer (Applied Biosystems, USA). The instrument was performed at a maximum accelerating potential of 20 kV and an m/z range from 700 to 4000. Six standards (Applied Biosystems, USA) were used as the internals to calibrate each spectrum to a mass accuracy within 0.1 Da. Protein candidate spots were analyzed using MALDI-TOF/TOF-MS in positive ion mode. Because the database of N. dietziae was uncompleted, proteins were identified by automated peptide mass fingerprinting using the Global Proteome Server Explorer software 3.0 (Applied Biosystems, USA) against a self-built protein sequence database of N. candida, N. coxensis DSM 45129 and Nonomuraea sp. SBT364 and some cyclosporinespecific P450 hydroxylases. For the algorithm, parameter settings were described by Wang et al. [42]. All the proteins identified were presented by MASCOT report protein scores for MS or total ion scores for MS/MS with greater than 95% confidence intervals.

Sample preparation of intracellular metabolites for GC-MS
The experimental data was obtained from four replicates of each treatment. The samples of mycelium at 48, 72, 96 and 120 h were harvested for quenching and extraction of intracellular metabolites at a low temperature. Subsequently, the extracts were derived with a two-step method. The methods of sample preparation for GC-MS had been previously described [43] in the Additional file 1 (sample preparation of intracellular metabolites).

Data acquisition and processing of GC-MS
GC-MS was performed by Agilent 6890N-5975C MSD system (Agilent Technologies, USA) equipped with a DB-5MS capillary column (30 m × 0.25 mm, 0.25 μm film thickness; Agilent Technologies, USA) and an autosampler. Parameter settings of GC-MS system were consistent with Wang' work [44]. Metabolomic data was processed with the Agilent MSD ChemStation (Agilent Technologies, USA) for spectrum deconvolution, denoising, retention time aligning, peak area integration, compound identification combined with NIST mass spectrum database (http://webbook.nist.gov/chemistry/) [44].