- Open Access
Engineering global transcription to tune lipophilic properties in Yarrowia lipolytica
Biotechnology for Biofuels volume 11, Article number: 115 (2018)
Evolution of complex phenotypes in cells requires simultaneously tuning expression of large amounts of genes, which can be achieved by reprograming global transcription. Lipophilicity is an important complex trait in oleaginous yeast Yarrowia lipolytica. It is necessary to explore the changes of which genes’ expression levels will tune cellular lipophilic properties via the strategy of global transcription engineering.
We achieved a strategy of global transcription engineering in Y. lipolytica by modifying the sequences of a key transcriptional factor (TF), SPT15-like (Yl-SPT15). The combinatorial mutagenesis of this gene was achieved by DNA assembly of up to five expression cassettes of its error-prone PCR libraries. A heterologous beta-carotene biosynthetic pathway was constructed to research the effects of combined Yl-SPT15 mutants on carotene and lipid production. As a result, we obtained both an “enhanced” strain with 4.7-fold carotene production and a “weakened” strain with 0.13-fold carotene production relative to the initial strain, nearly 40-fold changing range. Genotype verification, comparative transcriptome analysis, and detection of the amounts of total and free fatty acids were made for the selected strains, indicating effective tuning of cells’ lipophilic properties. We exploited the key pathways including RNA polymerase, ketone body metabolism, fatty acid synthesis, and degradation that drastically determined cells’ variable lipophilicity.
We have examined the effects of combinatorial mutagenesis of Yl-SPT15 on cells’ capacity of producing beta-carotene and lipids. The lipophilic properties in Y. lipolytica could be effectively tuned by simultaneously regulating genome-wide multi-gene expression levels. The exploited gene targets and pathways could guide design and reconstruction of yeast cells for tunable and optimal production of other lipophilic products.
Metabolic engineering and synthetic biology has developed versatile tools to engineer certain metabolic pathways to get optimal product synthesis and hyperburst of target functions [1,2,3,4,5]. However, evolution of complex traits in cells always demands arranging the expression levels of larger range of genes. If we can construct a wide change of gene expression, we can also understand why and how cells choose and accept special changes to meet target functions. An accessible strategy is global transcriptional machinery engineering (gTME), which was firstly reported by Stephanopoulos’ group [6, 7]. Gene sequences of transcription factors are subjected to mutation via error-prone polymerase chain reaction (PCR) or site-directed mutagenesis methods, leading to changeable promotor preferences and efficiency of RNA polymerase, which further modulates genome-wide transcription levels to varied extents. The diversity of changed transcriptomes confers diverse probable phenotypes on cells [6,7,8,9,10,11,12]. This strategy has already been used in several organisms, such as Escherichia coli, Saccharomyces cerevisiae, and Zymomonas mobilis, offering effective methods to generate novel mutants with enhanced environmental tolerance, metabolite production, and substrate utilization [6,7,8, 11, 12]. Besides, several transcription factors have been successfully engineered to optimize performance of strains, such as sigma factor, CRP, Spt15, H-NS, and Hha [10, 13,14,15]. In oleaginous yeast Yarrowia lipolytica, a similar strategy of global transcription engineering is in demand to tune cellular complex phenotypes such as lipophilicity which usually requires participation of large amounts of genes.
In recent years, Y. lipolytica attracts more attention as this yeast produces high-level lipids and its gene manipulation is available [16,17,18,19]. Metabolic pathway engineering has been used to improve production of several valuable products such as biofuels, polyunsaturated fatty acids, and carotenoids [20,21,22,23,24,25]. The experience obtained from rational engineering of model organism E. coli and S. cerevisiae has also been proved effective in Y. lipolytica, for example, for carotenoid synthesis. Wu et al.  improved the production of beta-carotene in E. coli to 44.2 mg/g DCW in flasks by combining strategies of modular pathway engineering and membrane engineering. Xie et al.  constructed a S. cerevisiae strain producing 1.156 g/L (20.79 mg/g DCW) of carotenoids using sequential control strategy in fermentation. Chen et al.  got 1.65 g/L (55.56 mg/g DCW) of lycopene from fermentation of S. cerevisiae by combining chassis engineering and heterologous pathway engineering. As for Y. lipolytica, the native tHMGR and ERG series genes were overexpressed in Gao’s work to improve fermentative production of beta-carotene to 4 g/L . A recent work reported that the strain highly accumulating lipids could also highly produce carotenoids . They achieved a fermentative production of beta-carotene of 6.5 g/L (90 mg/g DCW) with a concomitant production of 42.6 g/L of lipids. However, the deep tuning mode of competitive synthesis of fatty acids, lipids, and other heterologous lipophilic products like carotenoids was not clearly elucidated. Some distant pathway enzymes such as transporters and coenzymes could also tune lipid product synthesis [21, 30]. It was necessary to research how Y. lipolytica tunes its gene expression to affect internal lipophilicity to gain appropriate levels of native lipids and heterologous lipophilic products.
Here we report a strategy of engineering global transcription in Y. lipolytica to tune the complex lipophilicity and detailed lipophilic properties such as carotenoids, fatty acids and lipid bodies. In eukaryote’s general polymerase II transcription factor D (TFIID), a central Spt15 factor is TATA-binding protein (TBP), which is able to globally control the transcription of various metabolic and regulatory genes . This factor has four conservative functional regions, namely “repeat element 1,” “helix 2,” “repeat element 2,” and “helix 2’.” Some mutations of Spt15 resulted in pleiotropic changes of gene transcription levels [32, 33]. The SPT15 gene mutants were demonstrated feasibility to phenotypic evolution in S. cerevisiae in many ways, such as enhancing ethanol production , improving xylose fermentation , and enhancing adaptation to corn cob acid hydrolysate . For modifying Yl-Spt15 (a potential TF) in Y. lipolytica, we constructed combinatorial mutagenesis of its coding gene and screened the optimal mutant combinations for enhanced production of beta-carotene, which was synthesized as a non-native product by a heterologous pathway (Fig. 1). A 40-fold changing range of carotene production was obtained and the further analysis of the selected strains indicated effective tuning of lipophilic properties. We exploited the key changed pathways and offered guidance for tuning yeast lipophilicity to produce other high value-added lipophilic products.
Expression of Yl-SPT15 mutant libraries in assembled cassettes
According to the amino acid sequence of Spt15 in S. cerevisiae, we searched the conserved domains in Y. lipolytica by BLAST in NCBI database and got a Spt15-like protein named as Yl-Spt15 coded by YALIOB23056g on chromosome B (Fig. 1a and Additional file 1: Figure S1). This Yl-Spt15 owned the same conserved regions of “repeat element 1,” “helix 2,” “repeat element 2,” and “helix 2'” with Spt15 (Additional file 1: Figure S1, Additional file 2). Since the mutants of Spt15 in S. cerevisiae were proved effective to promote phenotypic evolution in many ways, we decided to construct the mutants of its potential counterpart Yl-Spt15 in Y. lipolytica to test its effects on product synthesis such as lipophilic products [10, 11, 34]. The manipulated Yl-SPT15 mutant libraries (see “Methods” for detailed construction process) were inserted in designed five expression cassettes and assembled meanwhile integrated in yeast chromosomal GUT2 site by in vivo homologous recombination (HR) (Fig. 1b, Additional file 1: Figure S2a, b and “Methods”). The correct rate of three-cassette assembly and integration at GUT2 site was about 7.3% and could be further improved to 18.0% in the strain with ku70 knockout (Additional file 1: Figure S3). Although the efficiency was not so high, this site primarily could be employed for integration of multiple Yl-SPT15 mutants. This kind of combinatorial mutagenesis of Yl-SPT15 permitted increased frequencies of both wild-type copies and mutant-type copies, offering expanded tuning range to identify dominant mutations that cause new functions in the presence of other unaltered genes.
To test whether the design of combinatorial Yl-SPT15 mutation was effective to tune metabolism of lipophilic products, we constructed a heterologous hybrid pathway to synthesize beta-carotene, a typical lipophilic product. Beta-carotene was used here for its lipophilic trait and visible color convenient for available yeast colony selection . The heterologous carotene would introduce a change in the metabolism and composition of cell’s lipophilic properties, and through selection of drastically changed carotene production in yeast colonies, we obtained a chance to learn how cells tuned both global gene transcription and intrinsic lipophilicity. The carotene pathway containing four modules of H0-EXP1p-crtE-XPR2t-H1, H1-TEFp-crtB-LIP2t-H2, H2-GPDp-crtI-OCTt-H3, H3-GPATp-crtY-PEX16t-H4 and a left module “rDNAL-Ura3-H0” and a right module “H4-rDNAR” were assembled and integrated in chromosomal rDNA site (Additional file 1: Figure S2c, the strain was named as Yl_ini). Both “rDNAL” and “rDNAR” were 700 bp referenced to Gao’s design . We got a similar phenomenon that the integration in rDNA was actually not exactly single copy. Quantitative PCR (Q-PCR) indicated that Yl_ini slightly underwent unexpected additional integrations of partial cassettes containing TEFp, GPDp, and GPATp (Additional file 1: Figure S4, Additional file 3: Table S1, “Methods”). However, we considered that these already obtained integrations were stable based on the nature of rDNA site, and the initial carotene production was also stable as detected. Based on this Yl_ini strain, we assembled sequential five cassettes of mutant libraries in yeast cells and successfully selected obvious color-changed yeast colonies. A most “enhanced” colony and a most “weakened” colony was selected by eye judgement for further analysis. It was proved that only cassettes of wild-type Yl-SPT15 did not affect colony’s color and carotene production (Additional file 1: Figure S5a).
Effects of combinatorial mutagenesis of Yl-SPT15 on beta-carotene production
The control strain containing cassettes of wild-type Yl-SPT15 was named as Yl_5_1, the visible “enhanced” strain was named as Yl_5_2, and the visible “weakened” strain was named as Yl_5_3 as follows (Fig. 2a). Firstly, the genotypes of selected strains were tested. According to the results of normal PCR and sequencing with special primers (“Methods” and Additional file 3: Table S2), the left integration site was exactly located upstream of GUT2’s open read frame (ORF), but the right integration site was not completely explicit. The designed homologous arms (H8) were 70% in accordance with the sequenced PCR products, but GUT2R could not be verified in alignment. This unexpected integration might lead to partial gene knockdown instead of knockout. The results of Q-PCR proved this speculation, as the expression of GUT2 was only partially suppressed to nearly half of previous level (Additional file 1: Figure S6 and Additional file 3: Table S1, “Methods”). By contrast, the expression of Yl-SPT15 was almost at the same level across all the tested strains including Yl_ini, Yl_5_1, Yl_5_2, and Yl_5_3, meaning the cells kept the Yl-SPT15’s transcription at a conservative level, although extra five cassettes were integrated.
The sequencing results showed that Yl_5_2 did receive more than one mutation in different cassettes, one (Glu208Ala) localized in the first module and three (Pro25Thr, Ala213Thr, Ala216Val) localized in the fifth module (Fig. 2b, Additional file 3: Table S3). The Yl_5_3 contained only one mutation (Phe194Leu) localized in the second module. All of the other modules at the expected locations in strain Yl_5_1, Yl_5_2, and Yl_5_3 were the cassettes of wild-type Yl-SPT15. As seen, four of all five mutations existed in the conserved transient region across “repeat element 2” (beta sheets) to “helix 2’” (alpha helix) [32, 36], indicating that this region was necessary to determine TF’s activity, similar within S. cerevisiae .
The production of beta-carotene was tested for these three strains (Fig. 2c). We chose defined media instead of rich media to remove interference caused by rich nutritional conditions. The carotene production in Yl_5_1 was 3.03 mg/L (0.87 mg/g DCW). Yl_5_2 got enhanced production as 12.34 mg/L (4.09 mg/g DCW, 4.7-fold of the Yl_5_1’s value), and Yl_5_3 got decreased production as 0.40 mg/L (0.11 mg/g DCW, 0.13-fold of the Yl_5_1’s value). The distinguished productions were directly reflected by the varied depth of strains’ yellow color (Fig. 2a). The growth rate was slightly different between strains (Additional file 1: Figure S5b). The Yl_5_2 with highest carotene production got lowest OD600, but the Yl_5_3 with lowest carotene production got highest OD600. An evaluation of the effects of each mutation on carotene production was also done. The exact single or combinatorial point mutants were introduced to natural Yl-SPT15 gene by PCR instead of directly cloning the sequencing-verified mutants from genome. Their individual expression cassette on plasmid pLD-EcYl was transformed into Yl_5_1. The carotene detection results suggested that the mutants in module 5 in Yl_5_2 and the mutant in module 2 in Yl_5_3 played significant role of impacting carotene production in each strain. It was also proved that the enhanced phenotype in Yl_5_2 was a result of combined efforts from different mutants (Additional file 1: Figure S7).
To test the influences of different levels of glucose, the strains were cultivated under 20, 30, and 40 g/L glucose keeping other conditions constant (Fig. 2d). Generally, the carotene yields increased as the glucose concentration increased, but there was slight impact on milligram per gram DCW values. When the glucose concentration was 40 g/L, we got the most enhanced production of carotene as 15.29 mg/L (4.29 mg/g DCW) in strain Yl_5_2, a limited 23.96% improvement. It was supposed that Y. lipolytica strains under SC-Ura-Leu fermentation could not present highest milligram per liter yield, and the SC other than YPD medium was used here for clearly detecting and analyzing cells’ internal change .
Exploration of key changed genes and pathways by transcriptome analysis
In a recent work, the transcriptome analysis was done for Y. lipolytica strains producing heterologous lycopene or not . Our work focused on the changes in different strains producing varied levels of beta-carotene. To reveal the effects of combinatorial mutation of Yl-SPT15 on global transcription, the transcriptome analysis was done for the three strains of Yl_5_1 (as control), Yl_5_2, and Yl_5_3. Primarily, according to the Pearson Correlation analysis, the reasonableness of samples and reliability of the detection were confirmed (Additional file 1: Figure S8). The transcriptome data was analyzed under two conditions, a normal condition of p < 0.05, and a stringent condition of p < 0.001 and change-fold < 0.5 or > 2.0. According to the normal condition (p < 0.05) results, both Yl_5_2 and Yl_5_3 exhibited changed expressions of almost 3000 genes, nearly half of the total genes throughout genome (Additional file 1: Figure S9; Additional file 4). This ratio was coincident with the ratio of the genes under the control of Polymerase II transcriptional system in S. cerevisiae. This highly proved that Yl-SPT15 did own the function of global TF. When the condition was stringent (p < 0.001 and fold < 0.5 or > 2.0), Yl_5_2 and Yl_5_3 exhibited altered expression of hundreds of genes throughout genome compared with the control (Fig. 3). Only one-tenth of these transcription-modulated genes (39 out of total 416) in two strains were in accordance, quite differentiated with that of normal condition (1943 out of total 4687). Most of the genes (164 out of total 196) with differential expression in Yl_5_2 were up-regulated, while the majority in Yl_5_3 (201 out of total 259) were down-regulated, indicating a whole “enhanced” transcription in Yl_5_2 but a whole “weakened” transcription in Yl_5_3. The highest two up-regulated targets in Yl_5_2 were YALI0D17556g(RER2) that expressed a major enzyme participating in polyprenol synthesis in both endoplasmic reticulum (ER) and lipid droplets, and YALI0F19184g(RGT2) that was plasma membrane high glucose sensor for regulating glucose transport.
The KEGG pathway enrichment was done for each pair of Yl_5_2 relative to Yl_5_1 and Yl_5_3 relative to Yl_5_1 (p < 0.05). The “enhanced” strain Yl_5_2 obtained some prominent enrichment of up-regulated genes in special pathways such as synthesis and degradation of ketone bodies, alpha-linolenic acid metabolism (an actual non-existing pathway in Y. lipolytica but related to three disperse real reactions later explained), mismatch repair, and homologous recombination (rich factors (RF) > 0.5) (Additional file 1: Figure S10 and Additional file 3: Table S4). The “weakened” strain Yl_5_3 gained some prominent enrichment of down-regulated genes in special pathways such as 10 amino acid metabolisms (valine, leucine, isoleucine, tryptophan, phenylalanine, tyrosine, lysine, histidine, cysteine, and methionine) and fatty acid degradation (RF > 0.6) (Additional file 1: Figure S11). The two strains shared an essential commonly changed pathway, RNA polymerase. The related genes in Y_5_2 were obviously down-regulated (RF = 0.62), while that in Yl_5_3 were notably up-regulated (RF = 0.55) (Additional file 1: Figures S12, S13). As mentioned above, another work concluded that the strain highly producing lipid bodies also highly produced carotene. However, in our work, a somewhat paradoxical phenomenon was observed. The genes associated with fatty acid biosynthesis in Yl_5_2 were obviously down-regulated (RF = 0.80), including the genes of YALI0C11407g(ACC1), YALI0B15059g(FAS1), YALI0B19382g(FAS2), and YALI0D17864g(FAA) (Additional file 1: Figure S13, Additional file 3: Table S5). By contrast, the fatty acid degradation was down-regulated in Yl_5_3 (RF = 0.78), including the genes of YALI0D17864g(FAA), YALI0D24750(ACX), YALI0F10857g(ACX), YALI0B10406g(ECH), YALI0E18568g(POT1), YALI0E11099g(ERG10), YALI0B08536g(ERG10), and YALI0F23749g(ACDH) (Additional file 1: Figure S11, Additional file 3: Table S6). There were two basic possibilities. First, the total carbon distribution to fatty acid synthesis and mevalonic acid (MVA) synthesis was constant and the increment of carotene caused down-regulated synthesis of fatty acids and vice versa. Second, the quantitative changes of carotene and fatty acids might be in accordance. The contradictory transcription regulation might indicate negative feedback regulation which commonly existed in fatty acid metabolism. The latter possibility was verified by later detection of the amounts of cellular fatty acids.
Modulation of the amounts of total fatty acids, free fatty acids, and lipid bodies
As transcriptome analysis revealed that fatty acid metabolism was regulated in both “enhanced” strain Yl_5_2 and “weakened” strain Yl_5_3, it was necessary to detect the amounts of fatty acids and lipids in these strains. The concentration of total fatty acids in Yl_5_2 was enhanced to 174.76 mg/L (58.03 mg/g DCW), 1.14-fold (1.32-fold) of the production contained in Yl_5_1 as 153.34 mg/L (44.11 mg/g DCW) (Fig. 4a, Additional file 1: Figure S14). As for free fatty acids, the concentration in Yl_5_2 was 29.73 mg/L (10.14 mg/g DCW), 1.15-fold (1.19-fold) of that in Yl_5_1 as 25.75 mg/L (8.55 mg/g DCW) (Fig. 4b, Additional file 1: Figure S15). By contrast, the concentration of total fatty acids in Yl_5_3 was reduced to 131.02 mg/L (36.73 mg/g DCW), 0.85-fold (0.83-fold) of that in Yl_5_1 (Fig. 4a, Additional file 1: Figure S14), and the free fatty acids was 23.54 mg/L (6.88 mg/g DCW) in Yl_5_3, 0.91-fold (0.80-fold) of that in Yl_5_1 (Fig. 4b, Additional file 1: Figure S15). The composition of FFA and total fatty acids was consistent among different strains (Additional file 1: Figures S14, S15). The increased production of both total and free fatty acids in Yl_5_2 might explain why its fatty acid synthesis pathway was down-regulated because this pathway owned typical mode of drastic feedback regulation by end products. On the other hand, the decreased production of fatty acids in Yl_5_3 explained why the fatty acid degradation pathway was down-regulated as cells had to down-regulate fatty acid degradation to protect its inner lipophilic environment. The microscopic observation showed that Yl_5_2 got more droplets of lipids but not a larger one compared with Yl_5_1 and Yl_5_3, as the increased lipophilic compartments were beneficial to carotene preservation and was created by whole cell’s complex tuning forces (Fig. 4c and Additional file 1: Figure S12).
Among the transcriptional regulated genes, ACC1, FAS1, and FAS2 played essential roles in the biosynthesis of fatty acids. The carboxylation of acetyl-CoA by ACC1 was considered as rate-limiting step, and this gene was usually overexpressed for lipid accumulation in S. cerevisiae and Y. lipolytica [21, 24, 25, 39, 40]. Besides, the coordinate expression of FAS1 and FAS2 was necessary for heteromultimeric fatty acid synthase complex . In Y. lipolytica, acyl-CoA was a key intermediate closely correlated with both free fatty acids and lipids (triacylglycerol, TAG), and its changing tendency was probably same as FFA and lipids (Fig. 4, Additional file 1: Figures S14, S15). As increased acyl-CoA plays the key role of negative feedback regulation of ACC1, it was understandable that the transcriptions of YALI0C11407g(ACC1), YALI0B15059g(FAS1), and YALI0B19382g(FAS2) were down-regulated in Yl_5_2 at the cell’s late growth stage. In addition, YALI0E11099g(ERG10) and YALI0F10857g(ACX) might play key roles as these two genes were regulated by obviously reverse trends in Yl_5_2 and Yl_5_3 (Figs. 4c, 5a, b). The two enzymes, acetyl-CoA C-acetyltransferase and acyl-CoA oxidase, also covered the crossroad locations for regulating distribution of acetyl-CoA flux in Y. lipolytica.
Enhancement of carotene production by overexpressing newly revealed intrinsic genes
As previously mentioned, in this work we revealed new up-regulated pathways in Yl_5_2 that might be correlated with the phenotype of higher carotene production. To determine these obviously up-regulated pathways in Yl_5_2 contributed to the improved beta-carotene production, we chose seven typical up-regulated genes from the two pathways and transformed their individual expression cassette on plasmids into Yl_5_2 (Fig. 5a, b, Additional file 1: Figure S2c). The results showed that all the seven transformed strains with recombinant plasmids (from Yl_5_2_1 to Yl_5_2_7) gained higher beta-carotene production than the initial strain Yl_5_2 with blank plasmid (Yl_5_2_0). These genes did increase the flux of acetyl-CoA towards terpene synthesis. The increment was up to about 50–60%, although all the values were slightly reduced perhaps owing to the influences of cell’s growth on the medium with hygromycin B (Fig. 5c). The collaboration of four enzymes expressed individually from YALI0E11099g(EGR10), YALI0F30481g(ERG13), YALI0B22550g(HMGL), and YALI0F26587g(OXCT) generated a special recycling of acetyl-CoA in Y. lipolytica and its flux was obviously up-regulated in Yl_5_2 for enhancing the pull of acetyl-CoA towards HMG-CoA and offering more precursors to MVA pathway (Additional file 3: Table S4). The single-gene overexpression proved that the flux of this cycling could be further enhanced for improving carotene synthesis and therein the YALI0E11099g(EGR10) might be the bottleneck. Another up-regulated pathway alpha-linolenic acid metabolism was actually a mistake of KEGG enrichment analysis as it did not exist in Y. lipolytica. The actual functions of the coupling of YALI0F10857g(ACX) and YALI0E18568(POT1) shall be reacting repeated beta-oxidation of acyl-CoA to produce final acetyl-CoA (Fig. 5b, Additional file 3: Table S4). Another gene YALI0F10010g(TGL4) actually charged degradation of multiple compounds to recycle free fatty acids (Fig. 5b, Additional file 3: Table S4). Single-gene especially YALI0E18568(POT1) overexpression showed that it was beneficial to further enhance this pathway to gain more carotene products.
We have designed a strategy of global transcription engineering, like gTME, in Y. lipolytica for tuning lipophilic properties. The gTME promotes identification of genes requiring perturbation for optimizing cell’s complex phenotypes, which has been proved previously [6, 12, 34, 42]. Even recently, the heterologous transcriptional sigma factor has been modified to make E. coli accept amazing direct transcription of promoters from other prokaryotic resources. Here we constructed the coexistence of different ratios of wild-type and mutated Yl-SPT15 genes in Y. lipolytica, permitting exploration of dominant mutations in the condition that other chromosomal genes could be unaltered. Fortunately, we obtained the assured changes of both genotype and phenotype in the selected strains with enhanced carotene production or weakened carotene production.
The results of transcriptional profiling revealed that the “enhanced” strain Yl_5_2 and the “weakened” strain Yl_5_3 exhibited differential expression of large amounts of genes relative to Yl_5_1 as control in the SC medium. Except for the P25T in Yl_5_2-module 5, the other four mutants all localized in the region from 194 to 216, which was corresponding to the region from 207 to 229 in S. cerevisiae (Additional file 1: Figure S1). Three out of total four mutants in Yl_5_2 (E208A in module 1 and A213T, A216V in module 5) were in “helix 2’” region. The only mutant F194L in Yl_5_3-module 2 was in the “repeat element 2” region. This phenomenon was coincident with previous work as this region was charged in the correlation with DNA sequence and other functional proteins [6, 35, 36].
The transcriptional changes in Yl_5_2 and Yl_5_3 were quite broad, exhibiting respective enrichment of certain pathways and functional gene groups (Fig. 3 and Additional file 1: Figures S10–S13). Yl_5_2 gained prominent enrichment of up-regulated pathways especially the metabolism of ketone bodies and three disperse reactions enriched in the fictitious alpha-linolenic pathway, but also a down-regulated pathway of fatty acid synthesis relative to Yl_5_1 (control). Yl_5_3 got prominent enrichment of down-regulated pathways especially up to 10 amino acid metabolisms and fatty acid degradation relative to Yl_5_1 (control). RNA polymerase was quite differentiated between the two strains as was down-regulated in Yl_5_2 but up-regulated in Yl_5_3. Further detection of total and free fatty acid contents in each selected strain proved that the total lipophilic products were tuned along with the tuned global transcription.
Seven genes in the up-regulated pathways with highest rich factors in Yl_5_2 were investigated as overexpression targets in this strain again. It was certified that the pathway of ketone body metabolism and three other disperse reactions were contributable to the improvement of beta-carotene production (Fig. 5). The results here showed difference from the S. cerevisiae gTME for screening strains with high production and tolerance of ethanol. In Alper’s work, no particular pathways or genetic networks were discovered as predominately responsible targets . A reasonable explanation was that no stressed conditions were used in our work, although somewhat weakening the selection strength, providing chances to reveal those predominately and progressively tunable pathways under normal condition . Our work offered some new targets to tune cell’s lipophilic products, as well as the previous approaches of molecular construction [21, 29]. It was a fact that the yields of beta-carotene in this work were still far lower than those in Gao’s work and Larroude’s work [23, 29]. Here only SC medium with minimal nutrients was used to cultivate cells, resulting not in optimal carotene production but offering a chance to research modes of tuning cellular lipophilic products.
A speculated model was that cells tuned special pathways to cope with competitive relations between lipophilic carotene and fatty products (Fig. 6). It was assured that there was a close correlation around the crossroad of acetyl-CoA between synthesis of carotenoids and fatty acid products, both as lipophilic properties sharing the common lipophilic space in cells. Primarily, the strain overproducing lipids also produced more beta-carotene after the complex tuning of global gene transcription. The proved highest effective resources of acetyl-CoA for synthesizing beta-carotene in Yl_5_2 were ketone body metabolism (pathway 1 in Fig. 6) and fatty acid beta-oxidation (pathway 2 in Fig. 6). This made the cells pull more fatty acids from lipid bodies. Meanwhile, the whole tuning tendency led to increment of lipid bodies generating more membrane structure versus volume, which might be beneficial to the storage of newly produced beta-carotene (Fig. 4c and Additional file 1: Figure S16). The driving force from pathway 1 and 2 also increased the amounts of free fatty acids and conversely down-regulated the synthesis pathway by negative feedback regulation . The increment of fatty acids and the number of lipid droplets in the “enhanced” strain potentially indicated the phenotype of co-evolution of native and heterologous lipophilic product synthesis. Next, we would research on how the cell tunes itself between “multiple lipid droplets” stage and “single lipid droplet” stage. Our research offered another angle to study how to firstly construct a robust oil-rich chassis and then transfer it into a producer of target lipophilic products, as well as other strategies [30, 44].
In this work a global transcription engineering strategy is designed and operated in Y. lipolytica, an oleaginous yeast. These results demonstrate the availability of gTME in this yeast for altering complex phenotypes such as lipophilicity. The synthesis of beta-carotene and fatty acid products can be meanwhile tuned by a close correlation model. The strategy is proved a promising tool for the gene manipulation in Y. lipolytica to optimize other desirable high value-added lipophilic metabolites.
Strains and media
The Y. lipolytica starting strain was ATCC 201249 (MATA ura3-302 leu2-270 lys8-11 PEX17-HA) from American Type Culture Collection (ATCC). All the strains and plasmids newly constructed in this work were given at Additional file 3: Table S7. Y. lipolytica cells were grown at 28 °C on petri dish (pH 6.5) or in liquid SC medium (pH 6.0) consisting of 20 g glucose, 6.7 g yeast nitrogen base (YNB), 2 g mixed amino acid powders only lacking both uracil and leucine (per liter). The fermentation was done as follows. First day, the cells were cultivated overnight. Next day, the cells were inoculated in new medium until 12 h and the harvested cells were inoculated again as initial OD = 0.1 in 50 mL fermentation medium in flask and cultivated for 96 h. E. coli competent cells (TransGen Biotech Company) for cloning DNA were grown at 37 °C in normal LB medium. Phosphatic buffer solution consisted of 137 mM NaCl, 2.7 mM KCl, 4.3 mM Na2HPO4, and 1.4 mM KH2PO4.
Construction of Yl-SPT15 mutant libraries
The mutation of Yl-SPT15 was done by error-prone PCR. The wild-type YL-SPT15 was primarily ligated in the plasmid pEASY-Blunt (TransGen Biotech Company). Then an error-prone PCR was operated using this constructed plasmid as template. The 100 µL PCR system contained 10 µL 10× FastTaq Buffer, 10 µL EP dNTP mixture (containing 5 mM dCTP, 5 mM dTTP, 2 mM dATG, and 2 mM dGTP), 0.05 mM M13F primer, 0.05 mM M13R primer, 2 µL FastTaq enzyme, 0.15 mM MnCl2, 1 µL template DNA, and ddH2O. The reaction was operated as follows: initial template denaturation at 94 °C for 3 min, followed by 30 cycles of denaturation at 94 °C for 30 s, primer annealing at 50 °C for 30 s, and elongation at 72 °C for 30 s, following that, a last amplification at 72 °C for 10 min was run. After that, the PCR products of Yl-SPT15 mutant libraries were purified after gel electrophoresis and digested by BsmBI enzyme. This digestion made all the DNA segments own standard left sticky end of 5′-AATG-3′ and right sticky end of 5′-TTTA-3′.
Assembly of expression cassettes of Yl-SPT15 libraries and integration into chromosome
After digested by BsmBI, the mixed segments were ligated into five cassettes (Additional file 1: Figure S2) to get sequential expression modules. All the cassettes on plasmid pLD2 (modified from pEASY-Blunt with mutation of BsaI and BsmBI sites) owned standard left 5′-CATT-3′ and right 5′-TAAA-3′ sticky ends for ligation with Yl-SPT15 libraries. The constructed expression cassettes were H3-GPATp-Yl-SPT15 library-PEX16t-H4, H4-YAT1p-Yl-SPT15 library-LIP1t-H5, H5-XPR2p-Yl-SPT15 library-PEX20t-H6, H6-FBAp-Yl-SPT15 library-CYC1t-H7, H7-LEUp-Yl-SPT15 library-ACOt-H8. Most promoters and terminators were cloned according to Blazeck’s work except for some little modification . All “Hx” were homologous arms at the length of 200 bp with low G/C contents (20–30%), which were all cloned directly from S. cerevisiae genome (primers were listed in Additional file 3: Table S2). The transformed E. coli colonies were scraped off the plates and plasmids were isolated. Then the cassettes were digested by NotI from plasmid backbones and the total DNA mixture was directly purified without running gel electrophoresis. 20 μL of each segment and an extra left module of “GUT2L-Leu2-H3” and a right module of “H8-GUT2R” were mixed to form a total 74 μL for further transformation into the Y. lipolytica strain already containing beta-carotene pathway (Additional file 1: Figure S2, the modules were constructed by a similar way and then assembled in yeast).
The “GUT2L” was 1000 bp upstream of GUT2’s ORF, and “GUT2R” was 1000 bp downstream of GUT2’s ORF. The Y. lipolytica transformation method was generally similar as the “LiAc/SS-DNA/PEG procedure” used for S. cerevisiae . The total system was 360 μL including 74 μL DNA fragment, 10 μL boiled ss-carrier DNA (10 mg/mL), 240 μL 50% PEG, and 36 μl 1 M LiAc. The transformed cells were plated on SC–Ura–Leu plates and incubated for 3–4 days at 28 °C until the colonies were present.
Detection of Yl-SPT15 mutant genotypes in selected yeast strains
The yeast colonies with desired colors could be picked best after a week until the colors completely deposited. The genome of selected yeast cells was extracted using normal methods. The Yl-SPT15 mutants in all five sequential cassettes were cloned from genome by a normal PCR method (Primers were shown in Additional file 3: Table S3). Also, the integration site GUT2 was verified using the upstream primer pair of “GUT2L” and “HR-F,” and downstream primer pair of “GUT2R” and “HR-R” (Additional file 3: Table S3). The PCR enzyme was FastPfu purchased from TransGen Biotech Company. The amplificated DNA products were purified after gel electrophoresis and sent to GENEWIZ Company for sequencing.
The concentration of the extracted whole-cell RNA was measured with Qubit® RNA Assay Kit in Qubit® 2.0 Fluorometer (Life Technologies, CA, USA). Then an amount of 3 μg RNA per sample was used as input. Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA). Finally, the PCR products were purified using AMPure XP system and the library quality was assessed by Agilent Bioanalyzer 2100 system. Then the libraries of PCR products were sequenced on an Illumina Hiseq 4000 platform and 150 bp paired-end reads were generated. The data analysis was all done in Novogene Company. Especially, the enrichment of differential expression genes in KEGG pathways (http://www.genome/jp/kegg/) was tested by KOBAS software.
Extraction of cellular beta-carotene from Y. lipolytica cells
At the end of fermentation, the cells were collected and the supernatant was discarded after centrifugation. One milliliter of 3 M hydrochloric acid was added to the centrifuge tubes containing cells, shocked to blending. The centrifuge was put in boiling water for 5 min, then immediately put in ice bath for 5 min. The supernatant was discarded after centrifugation. Then the cells were washed with sterile water. One milliliter of acetone containing 1 g BHT (butylated hydroxytoluene) per liter was added to the cells and the tubes were shaken for 10 min. The supernatant was filtered with 0.2 μm hydrophobic membranes for collecting carotene products.
Extraction of cellular total fatty acid from Y. lipolytica cells
Cells were harvested after centrifugation. One milliliter of methanol containing 3 M HCl and 10 mg/L heptadecanoic acid and 100 μL chloroform were added to the centrifuge tube containing cells, shocked to blending. Then the tubes were incubated at 70 °C for 3 h, reversed several times per 40 min. After that, the tubes were cooled down to room temperature. Then sodium chloride particle was added to take saturation, shaken for 1 min. After added with 500 μL hexyl hydride, the tube was shaken for a last minute. The supernatant was filtered with 0.2 μm hydrophobic membranes for collecting fatty ester products.
Extraction of cellular free fatty acid from Y. lipolytica cells
800 μL of cell culture was taken into centrifuge tube from 96 h incubated cultures, then 400 μL dichloromethane containing 200 mM methyl iodide and 30 mg/L pentadecanoic acid as an internal standard was added into the tube and 20 μL 40% tetrabutylammonium hydroxide (base catalyst) was added immediately. The mixtures were shaken for 30 min at 1400–1600 rpm using a vortex mixer, and then centrifuged at 6000 rpm to promote phase separation. A 300 μL dichloromethane layer was transferred into a new centrifuge tube, and evaporated 1 h to dryness. The extracted methyl esters were resuspended in 150 μL hexane and then filtered with 0.2 μm hydrophobic membranes for collecting fatty ester products [46, 47].
Measurement of beta-carotene
The quantitative determination of beta-carotene was achieved by HPLC. 10 μL sample was injected into a series e2695 HPLC, combined with a 2489 mass spectrometer (Waters, UV, USA), on a Ascentis® Express C18 column (5 cm × 2.1 mm × 2.7 μm). The mobile phase consisted of 80% methanol, 18% acetonitrile, and 2% dichloromethane. HPLC conditions were mobile phase flow rate at 0.3 mL per min, sample room temperature programed at 22 °C, column temperature at 25 °C. The detection wavelength was at 450 nm. Every tested strain had three parallel samples. The beta-carotene standard sample was provided by SIGMA.
Measurement of fatty acids
The quantitative detection of fatty acids was available by GC–MS method. 10 μL sample was injected into a series 6890N GC (Agilent Corp., USA), combined with a GCT Premier™ MICROMASS mass spectrometer (Waters Corp., USA), on a DB-5 ms column with a split injection of 5:1. GC conditions were carrier gas (helium) at 91 kPa per min in constant pressure mode, oven temperature programed from 70 °C (3 min) to 280 °C at the rate of 3 °C per min, source temperature at 250 °C, and interface temperature at 250 °C. Electron impact (EI) spectra were obtained at − 70 eV. GC–MS raw data were analyzed using the software package Masslynx4.1 (Waters Corp., USA), avoiding detector overload and isotope fractionation as described.
Q-PCR was carried out on a CFX96 Cycler-Real Time PCR Detection System (Bio-Rad Laboratories, Inc., Hercules, CA, USA), in white-walled PCR plates (96 wells). A ready-to-use master mix contained a fast proof-reading Polymerase, dNTPs, stabilizers, MgCl2, and SYBR® Green dye and were used according to the manufacturer’s instructions (Bio-Rad). Reactions were prepared in a total volume of 18 μL containing 400 nM each primer (MWG), 2× SsoAdvanced™ SYBR® Green Supermix (Bio-Rad), and 2 μL cDNA. The cycle conditions were set as follows: initial template denaturation at 95 °C for 3 min, followed by 40 cycles of denaturation at 95 °C for 10 s, and combined primer annealing/elongation at 57 °C for 20 s. The amount of fluorescence for each sample, given by the incorporation of SYBR® Green into dsDNA, was measured at the end of each cycle and analyzed via CFX Manager™ software 2.1 (Bio-Rad Laboratories, Inc.). With Yl_ini as control, the copies of gene GUT2 and Yl-SPT15 were represented by the expression level of target gene divided by that of ACT gene which performed as the internal reference gene . In another Q-PCR for quantification of carotene pathway gene integrations, the blank strain ATCC 201249 was chosen as control and the ACT was still used as a reference gene, in contrast, the Q-PCR was done on only genome extraction sample.
- S. cerevisiae :
- Y. lipolytica :
global transcription machinery engineering
polymerase chain reaction
Luo Y, Li BZ, Liu D, Zhang L, Chen Y, Jia B, et al. Engineered biosynthesis of natural products in heterologous hosts. Chem Soc Rev. 2015;44(15):5265–90.
Song MC, Kim EJ, Kim E, Rathwell K, Nam SJ, Yoon YJ. Microbial biosynthesis of medicinally important plant secondary metabolites. Nat Prod Rep. 2014;31(11):1497–509.
Ahmadi MK, Pfeifer BA. Recent progress in therapeutic natural product biosynthesis using Escherichia coli. Curr Opin Biotechnol. 2016;42:7–12.
d’Espaux L, Mendez-Perez D, Li R, Keasling JD. Synthetic biology for microbial production of lipid-based biofuels. Curr Opin Chem Biol. 2015;29:58–65.
Rogers JK, Taylor ND, Church GM. Biosensor-based engineering of biosynthetic pathways. Curr Opin Biotechnol. 2016;42:84–91.
Alper H, Moxley J, Nevoigt E, Fink GR, Stephanopoulos G. Engineering yeast transcription machinery for improved ethanol tolerance and production. Science. 2006;314(5805):1565–8.
Alper H, Stephanopoulos G. Global transcription machinery engineering: a new approach for improving cellular phenotype. Metab Eng. 2007;9(3):258–67.
Lam FH, Hartner FS, Fink GR, Stephanopoulos G. Enhancing stress resistance and production phenotypes through transcriptome engineering. Method Enzymol. 2010;470:509–32.
Liu HM, Xu L, Yan M, Lai CG, Ouyang PK. gTME for construction of recombinant yeast co-fermenting xylose and glucose. Chin J Biotechnol. 2008;24(6):1010–5.
Liu HM, Yan M, Lai CG, Xu L, Ouyang PK. gTME for improved xylose fermentation of Saccharomyces cerevisiae. Appl Biochem Biotechnol. 2010;160(2):574–82.
Liu HM, Liu K, Yan M, Xu L, Ouyang PK. gTME for improved adaptation of Saccharomyces cerevisiae to corn cob acid hydrolysate. Appl Biochem Biotechnol. 2011;164(164):1150–9.
Tan F, Wu B, Dai L, Qin H, Shui Z, Wang J, et al. Using global transcription machinery engineering (gTME) to improve ethanol tolerance of Zymomonas mobilis. Microb Cell Fact. 2016;15(1):1–9.
Klein-Marcuschamer D, Stephanopoulos G. Assessing the potential of mutational strategies to elicit new phenotypes in industrial strains. Proc Natl Acad Sci USA. 2008;105(7):2319–24.
Zhang H, Chong H, Ching CB, Song H, Jiang R. Engineering global transcription factor cyclic AMP receptor protein of Escherichia coli for improved 1-butanol tolerance. Appl Microbiol Biotechnol. 2012;94(4):1107–17.
Hong SH, Lee J, Wood TK. Engineering global regulator Hha of Escherichia coli to control biofilm dispersal. Microb Biotechnol. 2010;3(6):717–28.
Blazeck J, Liu L, Redden H, Alper H. Tuning gene expression in Yarrowia lipolytica by a hybrid promoter approach. Appl Environ Microbiol. 2011;77(22):7905–14.
Gao S, Han L, Zhu L, Ge M, Yang S, Jiang Y, et al. One-step integration of multiple genes into the oleaginous yeast Yarrowia lipolytica. Biotechnol Lett. 2014;36(12):2523–8.
Schwartz CM, Hussain MS, Blenner M, Wheeldon I. Synthetic RNA polymerase III promoters facilitate high-efficiency CRISPR–Cas9 mediated genome editing in Yarrowia lipolytica. ACS Synth Biol. 2016;5(4):356–9.
Celińska E, Ledesma-Amaro R, Larroude M, Rossignol T, Pauthenier C, Nicaud JM. Golden Gate assembly system dedicated to complex pathway manipulation in Yarrowia lipolytica. Microb Biotechnol. 2017;10(2):450–5.
Qiao K, Abidi SHI, Liu H, Zhang H, Chakraborty S, Watson N, et al. Engineering lipid overproduction in the oleaginous yeast Yarrowia lipolytica. Metab Eng. 2015;29:56–65.
Qiao K, Wasylenko TM, Zhou K, Xu P, Stephanopoulos G. Lipid production in Yarrowia lipolytica is maximized by engineering cytosolic redox metabolism. Nat Biotechnol. 2017;35(2):173–7.
Xue Z, Sharpe PL, Hong SP, Yadav NS, Xie D, Short DR, et al. Production of omega-3 eicosapentaenoic acid by metabolic engineering of Yarrowia lipolytica. Nat Biotechnol. 2013;31(8):734–40.
Gao S, Tong Y, Zhu L, Ge M, Zhang Y, Chen D, et al. Iterative integration of multiple-copy pathway genes in Yarrowia lipolytica for heterologous β-carotene production. Metab Eng. 2017;41:192–201.
Ledesma-Amaro R, Dulermo R, Niehus X, Nicaud JM. Combining metabolic engineering and process optimization to improve production and secretion of fatty acids. Metab Eng. 2016;38:38–46.
Tai M, Stephanopoulos G. Engineering the push and pull of lipid biosynthesis in oleaginous yeast Yarrowia lipolytica for biofuel production. Metab Eng. 2013;15(1):1–9.
Wu T, Ye LJ, Zhao DD, Li SW, Li QY, Zhang BL, et al. Membrane engineering—a novel strategy to enhance the production and accumulation of β-carotene in Escherichia coli. Metab Eng. 2017;43:85–91.
Xie WP, Ye LD, Lv XM, Xu HM, Yu HW. Sequential control of biosynthetic pathways for balanced utilization of metabolic intermediates in Saccharomyces cerevisiae. Metab Eng. 2015;28:8–18.
Chen Y, Xiao WH, Wang Y, Liu H, Li X, Yuan YJ. Lycopene overproduction in Saccharomyces cerevisiae through combining pathway engineering with host engineering. Microb Cell Fact. 2016;15(1):113–25.
Larroude M, Celinska E, Back A, Thomas S, Nicaud JM, Ledesma-Amaro R. A synthetic biology approach to transform Yarrowia lipolytica into a competitive biotechnological producer of β-carotene. Biotechnol Bioeng. 2018;115(2):464–72.
Wasylenko TM, Ahn WS, Stephanopoulos G. The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica. Metab Eng. 2015;30:27–39.
Seong YJ, Park H, Yang J, Kin SJ, Choi W, Kim KH, et al. Expression of a mutated SPT15 gene in Saccharomyces cerevisiae, enhances both cell growth and ethanol production in microaerobic batch, fed-batch, and simultaneous saccharification and fermentations. Appl Microbiol Biotechnol. 2017;101(9):3567–75.
Cang Y, Auble DT, Prelich G. A new regulatory domain on the TATA-binding protein. EMBO J. 1999;18(23):6662–71.
Eisenmann DM, Dollard C, Winston F. SPT15, the gene encoding the yeast TATA binding factor TFIID, is required for normal transcription initiation in vivo. Cell. 1989;58(6):1183–91.
El-Rotail AAMM, Liang Z, Li Y, Liu SP, Shi GY. A novel constructed SPT15 mutagenesis library of Saccharomyces cerevisiae by using gTME technique for enhanced ethanol production. AMB Express. 2017;7(1):111–22.
Özaydın B, Burd H, Lee TS, Keasling JD. Carotenoid-based phenotypic screen of the yeast deletion collection reveals new genes with roles in isoprenoid production. Metab Eng. 2013;15(1):174–83.
Chasman DI, Flaherty KM, Sharp PA, Kornberg RD. Crystal structure of yeast TATA-binding protein and model for interaction with DNA. Proc Natl Acad Sci USA. 1993;90(17):8174–8.
Sestric R, Munch G, Cicek N, Sparling R, Levin DB. Growth and neutral lipid synthesis by Yarrowia lipolytica on various carbon substrates under nutrient-sufficient and nutrient-limited conditions. Bioresour Technol. 2014;164(7):41–6.
Zhao C, Gao Q, Chen J, Wei L, Imanaka T, Hua Q. Metabolomic changes and metabolic responses to expression of heterologous biosynthetic genes for lycopene production in Yarrowia lipolytica. J Biotechnol. 2017;251:174–85.
Liu YB, Yap SA, Koh CM, Ji LH. Developing a set of strong intronic promoters for robust metabolic engineering in oleaginous Rhodotorula (Rhodosporidium) yeast species. Microb Cell Fact. 2016;15(1):200–8.
Shi SB, Chen Y, Siewers V, Nielsen J. Improving production of malonyl coenzyme A-derived metabolites by abolishing Snf1-dependent regulation of Acc1. MBio. 2014;5(3):e01130-14.
Kajiwara S, Oura T, Shishido K. Cloning of a fatty acid synthase component FAS1 gene from Saccharomyces kluyveri and its functional complementation of S. cerevisiae fas1 mutant. Yeast. 2001;18(14):1339–45.
Si HM, Zhang F, Wu AN, Han RZ, Xu GC, Ni Y. DNA microarray of global transcription factor mutant reveals membrane-related proteins involved in n-butanol tolerance in Escherichia coli. Biotechnol Biofuels. 2016;9(1):114–27.
Beopoulos A, Mrozova Z, Thevenieau F, Le Dall MT, Hapala I, Papanikolaou S, et al. Control of lipid accumulation in the yeast Yarrowia lipolytica. Appl Environ Microbiol. 2008;74(24):7779–89.
Xu P, Qiao K, Stephanopoulos G. Engineering oxidative stress defense pathways to build a robust lipid production platform in Yarrowia lipolytica. Biotechnol Bioeng. 2017;114(7):1521–30.
Gietz RD, Schiestl RH, Willems AR, Woods RA. Studies on the transformation of intact yeast cells by the LiAc/SS-DNA/PEG procedure. Yeast. 1995;11(4):355–60.
Zhou YJ, Buijs NA, Zhu ZW, Qin JF, Siewers V, Nielsen J. Production of fatty acid-derived oleochemicals and biofuels by synthetic yeast cell factories. Nat Commun. 2016;7:11709–17.
Steen EJ, Kang Y, Bokinsky G, Hu ZH, Schirmer A, McClure A, et al. Microbial production of fatty-acid-derived fuels and chemicals from plant biomass. Nature. 2010;463(7280):559–62.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(−Delta Delta C) method. Methods. 2001;25(4):402–8.
MW, DL, GNL, HL, and LZ designed and performed the experiments. XL, BZL, DL, and YJY conceived the project and guided writing the manuscript. All authors read and approved the final manuscript.
The authors are grateful for the financial support from the Ministry of Science and Technology of China (“973”Program: 2014CB745100), and the National Natural Science Foundation of China (21390203, 81502976).
The authors declare that they have no competing interests.
Availability of supporting data
The transcriptome raw data were uploaded to GEO and the Accession Number was GSE106929. Other data will be made available from the corresponding author on reasonable request.
Consent for publication
All authors read and approved the final manuscript.
Ethics approval and consent to participate
Ministry of Science and Technology of the People’s Republic of China (“973”Program: 2014CB745100), National Natural Science Foundation of China (21390203, 81502976).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.