INTRODUCTION
Squalene is an acyclic triterpenoid and a precursor for the biosynthesis of sterols, steroids, ubiquinones, and hopanoids (
1). Because of its unique properties, squalene is widely used in the cosmetic industry as an antioxidant. Squalene is also used as an emulsion for the delivery of vaccines and other therapeutic substances and as a protective and preventive agent in cancer treatment (
2,
3). Currently, squalene is mainly isolated from the liver oil of deep-sea sharks and plant oils (
2). However, owing to the increasing global demand for squalene and concerns about the protection of marine animals and the fate of food crops, alternative production of squalene by microbial fermentation has drawn much attention as a sustainable source (
4). So far, squalene production has been improved in many microorganisms. Common strategies in these works were engineering the methylerythritol-4-phosphate (MEP) or mevalonate (MVA) pathway to increase precursor pools and metabolic flux toward squalene and optimization of fermentation processes for increased accumulation of squalene. For example, the squalene content and volumetric titer of
Aurantiochytrium sp. strain 18 W-13a increased to 171 mg/g dry cell weight (DCW) and 0.9 g/liter via optimization of fermentation conditions (
5). In another study,
Pseudozyma sp. strain JCC207 was identified and shown to produce up to 340.52 mg/liter squalene (
6), and
Escherichia coli accumulated squalene up to 230 mg/liter or 55 mg/g DCW by coexpressing a chimeric mevalonate pathway with human or
Thermosynechococcus squalene synthase (
1,
6–8). The squalene production by engineered
Saccharomyces cerevisiae increased to 350 mg/g when the entire squalene biosynthesis pathway was overexpressed (
9). The squalene accumulation in yeast mainly depended on the flux of the MVA pathway, downregulation of squalene degradation, and increasing the capacity of yeast to store the accumulated squalene in lipid droplets (
10–12). Recently, we demonstrated that the increased lipid content of
S. cerevisiae, providing storage for squalene accumulation, led to a 250-fold increase in squalene production compared to a parental strain (
12). However, the limited storage capacity might still be a bottleneck to accommodate squalene in
S. cerevisiae due to the restricted intracellular lipophilic storage space.
The oleaginous yeast
Yarrowia lipolytica is an attractive host for the industrial production of hydrophobic compounds such as terpenoids due to its available genetic tools, large intercellular pool size of acetyl-coenzyme A (CoA), and rapid-growth phenotypes to high densities with high lipid content. Indeed, efficient production of a wide range of hydrophobic products, including lycopene, β-carotene, farnesene, and limonene, has been demonstrated in
Y. lipolytica (
13–17). Prior studies reported that
Y. lipolytica is capable of producing intracellular lipid bodies where large amounts of hydrophobic compounds can be accumulated. For instance, lycopene production could be improved through the manipulation of lipid body formation (
15). Additionally, Gao et al. and Larroude et al. reported a positive correlation between β-carotene and lipid contents and resulting high-level production of β-carotene in
Y. lipolytica (
18,
19).
As wild-type
Y. lipolytica can synthesize squalene, we sought a metabolic engineering strategy for the efficient production of squalene by engineered
Y. lipolytica. Specifically, genetic perturbations resulting in enhancement of the endogenous upstream MVA pathway and improved lipid accumulation were employed (
Fig. 1). Notably, we first characterized the biological function of
URE2 in
Y. lipolytica, which encodes a bifunctional protein that is involved in both nitrogen catabolite repression and oxidative stress response. Deletion of
URE2 led to improved squalene production, and double-gene deletion of
URE2 and
PEX10 resulted in improved squalene productivity.
DISCUSSION
Many metabolic engineering strategies have been employed to generate high terpenoid production in microbial hosts as a sustainable alternative to current production methods. Successful terpenoid production has been reported in the model microorganisms
E. coli and
S. cerevisiae through a multimodular strategy by engineering the MEP or MVA pathway and the redox cofactor supply modules and precursor production (e.g., acetyl-CoA) and restraining competing pathways and optimizing fermentation processes. In addition to the work on terpenoid production in model microorganisms, the potential of several other hosts, including the oleaginous yeast
Y. lipolytica, has been demonstrated in recent years. In
Y. lipolytica, successful efforts have been made to overexpress the structural genes in the MVA pathway for the improvement of terpenoid production. An effective strategy is the overexpression of HMG-CoA reductase, which is a rate-limiting enzyme in the MVA pathway. In
S. cerevisiae, the overexpression of
tHMG1 lacking a regulatory domain significantly improves terpenoid production (
44,
45). When it came to
Y. lipolytica, full-length
HMG1 was overexpressed to boost lycopene, limonene, and linalool production (
15,
17,
46). Kildegaard et al. compared the performance of complete and truncated
HMG1 variants and demonstrated that the expression of the full-length
HMG1 was more effective than truncated
HMG1 in enhancing β-carotene production (
47). In our previous study, the overexpression of endogenous
HMG1 led to a significant increase in squalene production up to 3-fold higher than that of the wild-type strain. However, the overexpression of heterologous
tHMG1 from
S. cerevisiae impaired squalene accumulation in
Y. lipolytica (
23). For high-level squalene production, extensive engineering of the MVA pathway flux is then bolstered by simultaneous overexpression of
HMG1 with other structural genes in the MVA pathway. In our study, the optimal combination was obtained from coexpression of
ERG9 and
ERG20 with
HMG1, resulting in a 4.2-fold increase of squalene content compared to the
HMG1-overexpressing strain and reaching 4.2 mg/g dry cell weight (DCW). To further increase the squalene productivity, enhancement of the
HMG1 gene dosage may be considered, because it was reported that the β-carotene content of
Y. lipolytica could be improved by optimizing the expression level of
HMG1, which subsequently enhanced HMG-CoA catalytic activity (
18). With the recently established CRISPR/Cas9 system, further optimization of the genetic tools for multicopy integration could be easily transferred to our metabolic engineering strategy for
Y. lipolytica (
48,
49).
Our previous studies on squalene production in
S. cerevisiae have suggested the introduction of genetic perturbations known to increase lipid contents to enhance squalene accumulation, as the lipid body can store squalene. Another study found overexpressed fatty acid synthases FAS1/FAS2 led to a significant increase in squalene production (
50). Other groups working on
Y. lipolytica reported a positive correlation between carotenoid synthesis and lipid content (
18,
19). Therefore, we speculated that the aggregation of lipid in
Y. lipolytica could also contribute to squalene overproduction. To this end, we turn our attention to lipid accumulation. While there have been significant efforts to increase lipid production in
Y. lipolytica and other yeasts, fewer efforts have been focused on engineering
Y. lipolytica to improve lipid metabolism and enable increase squalene accumulation. In the present study, peroxisomes were impaired by deletion of the peroxisomal biogenesis gene
PEX10 to increase lipid content and synergistically boosted squalene production in engineered
Y. lipolytica. However, several research groups have targeted peroxisomes as dynamic depots for the storage of terpenoids by introducing a complete MVA pathway in peroxisomes and achieved high-level productivity of terpenoids in
S. cerevisiae (
9,
51). Future work will be dedicated to examining the compartmentalization strategy in
Y. lipolytica for an enhanced production of acetyl-CoA derived chemicals.
In general, induction of lipid body accumulation in microorganisms, such as the oleaginous yeast, microalgae, and bacteria, is linked to nitrogen limitation (
52,
53). Thus, it is tempting to hypothesize that nitrogen catabolite repression (NCR) (i.e., the genes involved in the utilization and transport of poor nitrogen sources are repressed in the presence of preferred nitrogen sources) is involved in coordinating lipid metabolism with growth. Several NCR transcription factors have been identified, including the GATA-family transcription activators Gln3p and Gat1p, the repressors Dal80p and Gzf3p, and the zinc finger transcription factors Dal81p and Dal82p (
54,
55). Among these, Gln3p activity is inhibited by Ure2p, a negative nitrogen regulator that plays a central role in NCR in
S. cerevisiae under lush, nitrogen-replete conditions. As nitrogen becomes limiting or rapamycin is added, TOR is inactivated; consequently, Gln3p relocates to the nucleus and depresses NCR-sensitive transcription (
56–60). It has been assumed by several authors that inhibition of TOR by rapamycin mimics nitrogen starvation (
61,
62). A more recent study conducted on
S. cerevisiae showed that TOR inhibition results in fast lipid droplet accumulation with a lower growth rate, which, on the contrary, is inhibited by rapamycin in mammals (
63,
64). Their result further demonstrated that Gln3p and Gat1p could regulate the full induction of lipid droplets by rapamycin (
63). In this study, we investigated the regulatory roles of Ure2p in a naturally minimized genomic background (
Fig. 5). The data from qRT-PCR showed that the
GLN3 and
GAT1 transcription levels were higher in Δ
URE2 than WT strains in nitrogen-rich medium, which is consistent with previous findings in
S. cerevisiae (
31).
Nutrient restriction can induce lipogenesis with a significant reduction in cell growth in microorganisms. For economically feasible biobased lipid production, the metabolic uncoupling of cell growth from intracellular lipid accumulation is a prerequisite. Previous studies showed that higher rapamycin concentrations resulted in growth inhibition in
S. cerevisiae. In our study, with the increase in lipid accumulation, no obvious differences in optical density (OD) between the set of Δ
PEX10 Δ
URE2 mutant strains were observed during cultivation, indicating that cell growth was not affected by the abolishment of the
URE2 gene in the nitrogen-replete medium. This phenomenon is in striking contrast to the inhibition of TOR by supplementation with rapamycin on algae and
S. cerevisiae growth, where lipid is aggregated with decreased biomass formation. A recent study demonstrated, for the first time, that lipogenesis could be induced with the addition of rapamycin in an oleaginous yeast,
Trichosporon oleaginosus, without compromising growth (
65). Although the structure and functions of TOR complexes, such as TORC1 activation of translation and transcription, are conserved from yeast to mammals, the downstream processes are vastly different depending on the organism (
66). In
Y. lipolytica, such regulation of lipid metabolism by TOR is not yet well established. We proposed a model of the URE2-related regulatory network in
Y. lipolytica by combining the literature with our work (
Fig. 6). However, to gain deeper knowledge of the TOR signaling pathway associated with lipogenesis, additional experiments are needed.
Taken together, our findings demonstrate the potential of the combination of the MVA pathway optimization and lipid metabolic engineering for the overproduction of squalene in oleaginous organisms. The discovery of the latent role of URE2 in squalene accumulation enabled a novel strategy to increase squalene production that can be combined with other strategies to engineer cell factories for the efficient production of squalene in the future. More research into how the TOR signaling pathway is associated with lipogenesis will yield additional insights into the control and induction of lipid accumulation in oil yeasts.