Heterologous Expression of Two Malate Transporters From an Oleaginous Fungus Mucor circinelloides Improved the Lipid Accumulation in Mucor lusitanicus

The fungus, Mucor lusitanicus, is of great interest for microbial lipids, because of its ability to accumulate intracellular lipid using various carbon sources. The biosynthesis of fatty acid requires the reducing power NADPH, and acetyl-CoA, which is produced by the cleavage of citrate in cytosol. In this study, we employed different strategies to increase lipid accumulation in the low lipid-producing fungi via metabolic engineering technology. Hence, we constructed the engineered strain of M. lusitanicus CBS 277.49 by using malate transporter (mt) and 2-oxoglutarate: malate antiporter (sodit) from M. circinelloides WJ11. In comparison with the control strain, the lipid content of the overexpressed strains of mt and sodit genes were increased by 24.6 and 33.8%, respectively. These results showed that mt and sodit can affect the distribution of malate in mitochondria and cytosol, provide the substrates for the synthesis of citrate in the mitochondria, and accelerate the transfer of citrate from mitochondria to cytosol, which could play a significant regulatory role in fatty acid synthesis leading to lipids over accumulation.


INTRODUCTION
Oleaginous microorganisms, generally producing cellular lipids up to 20-80% of cell dry weight (CDW), can accumulate lipids when grown in the media with high C/N ratio, since the excess carbon sources can be converted into fatty acids, and stored in the form of triacylglycerols (TAGs) in the cells (Ratledge and Cohen, 2008). Nowadays, microbial lipids have gained an increasing interest as the alternative polyunsaturated fatty acids sources to fish oil, due to the pressure of the marine fish reserve depletion and ocean environmental problems (Li et al., 2008;Meng et al., 2009;Shi et al., 2011;Liang and Jiang, 2013). Mucor lusitanicus, previously classified as M. circinelloides f. lusitanicus (Wagner et al., 2020), is a model organism to study the mechanism of lipid accumulation (Ratledge and Wynn, 2002), gene silence (Nicolas et al., 2007), and light response (Silva et al., 2006). In addition, it is high in γ-linolenic acid (GLA, C18:3) production, which has many beneficial effects for human (Ratledge, 2013).
Gene manipulations and metabolic engineering have been successfully applied to enhance the production of GLA in this organism . During nitrogen starvation, the exceeding glucose is used to generate pyruvate, which then, enters mitochondria and participates in the biosynthesis of citrate under the catalyzing of citrate synthase (Ratledge and Wynn, 2002). Then the mitochondrial citrate is transported to cytosol and cleaved by ATP: citrate lyase (ACL) to generate acetyl-CoA, one of indispensable substrates for fatty acid synthesis. Whereas, the other product oxaloacetate is reduced to malate by malic dehydrogenase, which is transferred back to mitochondria in exchange for citrate (Evans et al., 1983). Meanwhile, the cytosolic malate can be converted to pyruvate by malic enzyme and produce NADPH, the essential supply of reducing power for fatty acid biosynthesis (Vongsangnak et al., 2012;Rodriguez-Frometa et al., 2013), which was once considered as a rate-limiting step in fatty acid synthesis of oleaginous microorganisms (Wynn et al., 1999;Song et al., 2001;Zhang et al., 2007). Hence, the malate pool regulates the activity of malic enzyme and the transportation of citrate in mitochondria. Thus, malate has an important role as the intermediate in fatty acid biosynthesis. Some studies have also proposed that malate is the key metabolite in a wide range of central metabolic processes, such as the maintenance of cellular osmolarity, pH, and stomatal regulation (Emmerlich et al., 2003;Fernie and Martinoia, 2009).
Although malate is an important intermediate for microbial cell, it could only go through the impermeable membrane equipped with malate transporters (MT) responded for the translocation of malate. Previously, MT have been studied in mammals (Gnoni et al., 2009), plants (Palmieri et al., 2008;Dolce et al., 2014), yeasts (Evans et al., 1983), and fungi (Thammarongtham et al., 2018). Recently, the malate transporter (named MT) located in the mitochondrial membrane was found for the first time by the genome-scale analysis of the metabolic networks in M. lusitanicus CBS 277.49 with a low lipid content (13% of CDW; Vongsangnak et al., 2013). Genetic modification work has shown that homologous overexpression of MT increased lipid content by 70%, and on the contrary, when it was knockout, the lipid content decreased by 27% (Zhao et al., 2016). Subsequently, 13 C-labeled metabolic flux analysis was used to demonstrate the mechanism of MT and suggested that it can lead to the increase of malate influx into mitochondria and promote lipid synthesis (Wang et al., 2019). Further investigation of the genome of M. lusitanicus CBS 277.49 found that one 2oxoglutarate: malate antiporter (named SODIT) could contribute for cytosolic malate and citrate exporting, which may mediate the distribution of malate and citrate used for fatty acid biosynthesis in cells (Yang et al., 2020). Very recently, these two MT with a high sequence identity [90.6% of MT, belonging to tellurite resistance/dicarboxylate transporter family and 66% of SODIT, belonging to 2.A.47 family (the divalent anion: Na + symporter (DASS) family or SLC13 family)] (Yang et al., 2020) were found in a high-lipid producing strain M. circinelloides WJ11, which could accumulate as much as 36% (w/w) lipid of CDW. However, the potential involvement of MT and SODIT of M. circinelloides WJ11 in promoting lipid synthesis in M. lusitanicus CBS 277.49 still remains unanswered.
In this study, the mt and sodit genes from M. circinelloides WJ11 were cloned and overexpressed in M. lusitanicus CBS 277.49. The growth characteristics, gene expressing levels and lipid accumulation of transformants were analyzed. This work is the first time to explore the roles of these two MT from high-lipid producing strain in a low-lipid accumulation strain and reveals new insights into the mechanism of fatty acid accumulation in M. lusitanicus and will provide more contribution on the development of M. lusitanicus as a model fungus.

Strains, Cultivation, and Transformation Conditions
The bacterium Escherichia coli Top10 was used for plasmids construction, preservation, and propagation. LB medium (10 g/L tryptone, 5 g/L yeast extract, 10 g/L NaCl, and 20 g/L agar) were used to grow E. coli and supplemented with ampicillin (100 µg/mL) when necessary. The uridine auxotroph pleu-MU402 derived from wild-type strain CBS 277.49 (Yang et al., 2020) was used as the parent strain in genetic modification experiments (Yang et al., 2019). The cDNA library of WJ11 (CCTCC No. M 2014424) was used as the source of the genes mt and sodit. The fungus was maintained on complete media YPG (20 g/L glucose, 10 g/L peptone, 3 g/L yeast extract, and 20 g/L agar; Bartnicki-Garcia and Nickerson, 1962) or uridineless media MMC (20 g/L glucose, 10 g/L casaminoacid, 0.5 g/L yeast nitrogen base without amino acids and ammonium sulfate, and 20 g/L agar; Nicolas et al., 2007), which were adjusted to pH 4.5 or 3.2 for the mycelia and colonial growth, respectively, and the culture temperature was set at 28 • C. Transformation was carried out as previously described (Gutierrez et al., 2011;Garre et al., 2015), and the transformants were grown at 28 • C in MMC medium supplemented with uracil at 200 mg/mL when required.

Plasmid Construction
Plasmid pMAT1552 (Zhao et al., 2015), contained the pyrG of CBS 277.49 was used as a selectable maker and the 1 kb fragments of up-and downstream of CarRP sequences, was adopted to construct the mt-overexpressing and soditoverexpressing plasmids. Mt and sodit genes were amplified from the cDNA of WJ11 with corresponding primers mt-F/R, sodit-F/R (Supplementary Table 1). The PCR fragments were inserted into plasmid pMAT1552 digested by XhoI after a strong promoter pzrt1 to generate gene overexpressing plasmids pMAT1552-mt and pMAT1552-sodit ( Figure 1A) using the One Step cloning kit (Takara).

Genomic DNA Preparation, RNA Isolation and Quantitative RT-qPCR
For DNA extraction and RNA isolation, M. lusitanicus transformants were grown for 4 days in fermenter with the modified K&R medium at 28 • C. After fermentation period, the samples were harvested by filtration, and then washed three times with distilled water. Genomic DNA was extracted from mycelium after crumbling under liquid nitrogen using the DNA Quick Plant System kit (Tiangen Biotech Co., Ltd). Mycelia, harvested at 24 h, were fully grounded with Trizol in liquid nitrogen to obtain the total RNA of the transformants. Then the total RNA was used to reverse transcribed to cDNA via ReverTra Ace qPCR RT Kit (Roche) as described in the instructions. RT-qPCR was carried out in LightCycler 96 (Roche) using the SYBR Green Realtime PCR Master Mix in accordance with the manufacturer's instructions, which was based on the method of 2 − Ct . Primers for RT-qPCR are listed on the Supplementary Table 1.

Determination of Biomass and Lipid Content
The biomass was collected on a weighed and dried filter paper by filtered via a Buchner funnel, purged three or four times with distilled water, frozen overnight at −80 • C and freezedried to constant weight, the weight of which was determined gravimetrically. Cellular lipids were extracted from ∼10 mg freeze-dried samples by chloroform/methanol (2:1, v/v; Folch et al., 1956), pentadecanoic acid (C15:0) was used as the internal standard, and then 10% HCl/methanol (w/w) was added for transesterification. The fatty acid methyl esters were separated with n-hexane (HPLC grade) and analyzed by gas chromatography (GC) which is equipped with a 30 m × 0.32 mm DB-Waxetr column. The GC program was set at 120 • C for 3 min and ramped up to 200 • C at the speed of 5 • C/min, then increased to 220 • C and hold for 2 min.

Determination of the Concentration of Glucose, Ammonium Ion and Citrate in the Media
In the case of M. lusitanicus transformants cultivations, the concentration of glucose and ammonium ion in media were determined by glucose oxidase Perid-test kit (Shanghai Rongsheng Biotech Co., Ltd.) and indophenol method as described (Chaney and Marbach, 1962), respectively. Meanwhile, the concentration of extracellular citrate was measured by using a citrate kit (Suoqiao Biotec.).

Determination of Extracellular Malate Concentration
The extracellular malate concentration in the medium for 72 h was measured by HPLC equipped with a 4.6 mm × 250 mm Diamonsil C18 column (Dikma, China; . The mixture of methanol/phosphate solution (Na 2 HPO 4 7.01 g/L; 2:98, v/v) served as mobile phase mobile phase. The initial column temperature was set at 25 • C, flow rate was 0.5 mL/min and the running time was 26 min.

Statistical Analysis
All the experiments were performed in three replicates and the results were presented as mean ± SD. Student's t test of SPSS 16.0 was used for statistical analysis of the data and p < 0.05 was considered as significant different.

Generation of mt-and sodit-Overexpressing Strains of M. lusitanicus by Genetic Engineering
The gene sequences of mt and sodit were found based on the genomic data of WJ11 and sequence alignment analysis in the National Center for Biotechnology Information (Yang et al., 2020). To investigate the involvement of mt and sodit in lipid accumulation in M. lusitanicus CBS 277.49, overexpression strains of these two genes were generated. The expression vector, pMAT 1552, contains the pyrG gene as a selectable maker, a strong promoter pzrt1 for regulating the expression of the target genes, and the up-and down-sequences of carRP. The control plasmid pMAT1552 and the target gene overexpressing plasmids pMAT1552-mt and pMAT1552sodit were transformed into defective strain MU402 to obtain the transformants Mc-1552, Mc-mt and Mc-sodit, respectively ( Figure 1A). The implementation of transformants selection was performed as described in a previous study (Rodriguez-Frometa et al., 2013). For each overexpression plasmid, three independent transformants, named Mc-mt-1, Mc-mt-2 and Mcmt-3 for pMAT1552-mt, Mc-sodit-1, Mc-sodit-2 and Mc-sodit-3 for pMAT1552-sodit, and one Mc-1552 as the control strain were selected, and the obtained recombinant transformants were verified by PCR analysis. The target genes, mt and sodit, and the 585 bp sequences of plasmid pMAT1552 were amplified by a primer pair 1552-F/R (Supplementary Table 1 , respectively, which were consistent with the control strain and the corresponding transformants. The PCR amplification results verified that the target genes, mt and sodit, were transferred into the fungus. The transformants were cultured in baffled flask with K&R medium for 4 days. Since there was little difference in lipid content among transformants (data not show), the highest lipid content of the obtained transformants, was selected for further experiments. Notably, in our preliminary findings of homologous overexpression of target strains, the lipid content was not changed significantly if compared with the heterologous overexpression. Thus, the heterologous overexpression was performed in further experiments, and the compared data in their lipid content between the homologous and heterologous overexpression were shown in Figure 2 (homologous overexpression related details were showed in Supplementary Figures 1, 2).  Table 1) were designed. As shown in Figure 3, the expression levels of the genes from WJ11, mt and sodit, in the recombinant transformants were significantly higher than that of the endogenous genes, which indicates that these genes were successfully overexpressed. However, the expression levels of endogenous genes were not affected.

Cell Growth and Lipid Content in mt-and sodit-Overexpressing Strains
The biomass and lipid content of mt-and sodit-overexpressing strains were measured when cultured in 2 L fermenters for 96 h (Figure 4). The growth patterns of the transformants were similar with the control strain. The biomass of the mt-overexpressing strain was lower than that of control strain throughout the fermentation, whereas the sodit-overexpressing strain had a higher biomass than the control strain in the early stages of fermentation but decreased after 72 h ( Figure 4A). Therefore, the lipid content of mt-and sodit-overexpressing strains was significantly increased (Figure 4B). The lipid accumulation patterns of three strains were almost similar, lipid accumulation increased rapidly during in the first 48 h, afterward entered a stable period, and then ceased after 72 h. The amounts of lipid accumulated in the mt-overexpressing strain reached 16.2%, an increase of 24.6% compared to the control strain, which was coincides with the results of the previous study (Zhao et al., 2016). While the lipid content of sodit-overexpressing strain was increased by 33.8% compared to Mc-1552 (from 13.0% in control to 17.4% in Mc-sodit), nevertheless, the sodit gene deletion in WJ11 experiments, showed lipid accumulation was increased, which indicated this gene had a different function in regulating malate transportation in oleaginous and non-oleaginous fungi . The fatty acid profiles of these transformants illustrated that GLA content in total fatty acids (TFAs) of mtoverexpressing strain (18.41% at 72 h) and sodit-overexpressing strain (18.66% at 72 h) were significantly lower than that of Mc-1552 (27.35% at 72 h; Figure 5). However, the content of hexadecanoic acid (C16:0) in TFAs of mt-overexpressing strain (23.80% at 72 h) and sodit-overexpressing (23.23% at 72 h) strain were significantly higher than that of strain Mc-1552 (17.10% at 72 h).

Glucose and Nitrogen Consumption and Extracellular Citrate Accumulation in mtand sodit-Overexpressing Strains
The consumption rate of glucose and nitrogen sources was nearly similar in all studied strains (Figures 6A,B). Notably, when The total RNA of the transformants at 24 h was extracted and the mRNA expressing levels was quantified by RT-qPCR. Column 1 and 2, the relative expressing level of mt gene located in the genome was qualified and amplified by cbs-mt-F/R primers. Column 3, the relative expressing level of mt gene in the plasmid was qualified and amplified by WJ11-mt-F/R primers. Column 4 and 5, the relative expressing level of sodit gene located in the genome was qualified and amplified by cbs-sodit-F/R primers. Column 6, the relative expressing level of sodit gene in the plasmid was qualified and amplified by WJ11-mt-F/R primers. The values were mean of three independent fermentation experiments. Error bars represent the standard error of mean. Values which showed different superscripts were significantly different to each other (p < 0.05).
the nitrogen was depleted at 12 h, the consumption glucose rate, was rapidly increased in transformants and the control strain ( Figure 6A). Interestingly, the amounts of extracellular citrate in the transformants were greatly increased during the whole cultivation period compared to the control. There were no significant differences between the amount of extracellular citrate in gene overexpressed transformants, Mc-sodit and Mcmt during the stage of cell growth and rapid lipid accumulation ( Figure 6C). Whereas, the concentration of extracellular citrate of Mc-sodit was much higher than that of Mc-mt after 48 h during the period of slow lipid accumulation, while the biomass decreased relatively as shown in Figure 6C. These results indicated that the intracellular glycolysis, lipid biosynthesis and citrate metabolism had been affected by the heterologous expression of mt and sodit genes and led to the increased accumulation of lipid in M. lusitanicus.

The Concentration of Extracellular Malate in mt-and sodit-Overexpressing Strains
To investigate whether overexpression of two genes (mt and sodit) from WJ11 could affect the malate metabolism, the extracellular malate concentration of samples cultured for 72 h in 2 L fermenter with K&R medium was analyzed by HPLC. As shown in Figure 7, the concentration of extracellular malate in sodit-overexpressing strain was higher compared to the control, whereas, the malate secretion in mt-overexpressing strain was the lowest, which illustrated that these two MT could have different regulation to malate transportation.

The Transcription Level of Key Genes for Fatty Acid Biosynthesis in mt-and sodit-Overexpressing Strains
The expression levels of key genes related to fatty acid synthesis at 24 h revealed the underlying mechanisms of lipid accumulation in mt-and sodit-overexpressing strains (Figure 8). There are five genes encoding malic enzymes which could convert malate to pyruvate and produce reducing power NADPH for fatty acid synthesis, three localized in mitochondria (mme1, gene ID: 78524; mme2, gene ID: 166127; and mme3, gene ID: 11639) and two localized in cytosol (cme1, gene ID:182779; cme2, gene ID: 186772; Vongsangnak et al., 2012). All of these genes were significantly up-regulated in both gene overexpressing strains, the transcription levels of mme1, mme2, and mme3 were higher in mt-overexpressing  strain, while the transcription levels of cme1 and cme2 were higher in sodit-overexpressing strain. There are two genes (fas1, gene ID: 72770; fas2, gene ID: 183529) encoding fatty acid synthase which could catalyze de novo fatty acid synthesis. As shown in Figures 8F,G, the transcription levels of fas1 were slightly down-regulated whereas the fas2 was significantly up-regulated in both strains. In addition, ACL (encoded by acl, gene ID: 110808), which could catalyze the cleavage of citrate to acetyl-CoA for fatty acid synthesis, was also analyzed. The transcription levels of acl were significantly up-regulated in mt-and sodit-overexpressing strains. The mitochondrial citrate transporter (encoded by ct gene, gene ID: 155787) and tricarboxylate transporter (encoded by tct gene, gene ID: 141738) expression levels are investigated, since they are involved in mitochondrial citrate transportation which could be affected by the malate metabolism. The transcription levels of ct and tct were significantly up-regulated in both strains, however, the expression levels of these two genes were higher in sodit-overexpressing strain, which was in agreement with the extracellular citrate concentration.

DISCUSSION
Lipid accumulation in oleaginous microorganisms is a result of a combination of multiple metabolic events. When the cells are exhausting a kind of essential nutrients (usually nitrogen source), the growth is inhibited by harsh conditions, thus the redundant carbon substrates continue to be taken up and converted to fatty acids and stored in the form of TAG (Ratledge and Wynn, 2002), therefore the lipid production increases significantly. Previous studies showed that nitrogen restriction leads to increase in AMP deaminase activity, which in turn reduces AMP concentration and decreases the activity of NAD + -isocitrate dehydrogenase, so that the tricarboxylic acid (TCA) cycle is retarded then the citrate was accumulated in mitochondria Ratledge and Wynn, 2002;Lazar et al., 2018). The excess mitochondrial citrate is exported to cytosol with malate as a counter-substrate (Huypens et al., 2011;Yang W. et al., 2021), then cleaved to oxaloacetate and acetyl-CoA which is the supplement of essential precursor for fatty acid synthesis (Ratledge, 2014). Consequently, the oxaloacetate is reduced by malate dehydrogenase to malate which is transported back into mitochondria by a MT to supply substrate for citrate synthesis, meanwhile the cytosolic malate also can be converted to pyruvate by malic enzymes, accompanied by the reducing force NADPH for fatty acid synthesis (Rodriguez-Frometa et al., 2013). In addition, intracellular malate can also be secreted into the environment by a plasma membrane MT when the cytosolic concentration is more than enough (Zhao et al., 2016). Therefore, the MTs directly affect the metabolisms in mitochondria and lipid biosynthesis.
Mucor lusitanicus was widely recognized as a model microorganism for investigating lipid synthesis because of its known complete sequenced genome and established genetic manipulation tools. Recently, the genome, proteome and biochemical activities of CBS 277.49 and WJ11 have been studied and compared (Tang et al., 2015a(Tang et al., ,b, 2017. By analyzing their genomic data, two different kinds of MTs were found in both strains, one is the MT located in mitochondria inner membrane and the other is the plasma membrane malate transporter (SODIT; Yang et al., 2020). The gene modification of the native mt in CBS 277.49 elucidated its vital role for lipid accumulation (Zhao et al., 2016). Considering the previous work, heterologous expression of the citrate transporters of WJ11 led to a lipid synthesis enhancement in the low-lipid producing fungus M. lusitanicus (Yang et al., 2019). Hence, heterologous expression of the mt and sodit of WJ11 might affect the cytosolic malate and regulate the lipid accumulation in CBS 277.49.
Based on our preliminary results of homologous overexpression, the results showed no significant differences in lipid content of mt-and sodit-overexpressing strains (20.2% in SD0081 and 21.2% in SD0082, respectively) compared to the control strain SD0080 (22% of CDW; Figure 2). Whereas, the heterologous overexpression of these strains was performed and had a significant influence on their lipid content in CBS 277.49. The lipid content was significantly increased by 24.6 and 33.8% in the Mc-mt and Mc-sodit, respectively (Figure 4). Furthermore, FIGURE 8 | (A-J) Expression of key genes in mt-and sodit-overexpressing strains. Relative mRNA levels were determined after 24 h of growth by RT-qPCR. mme, mitochondrial malic enzyme; cme, cytosolic malic enzyme; fas, fatty acid synthase; acl, ATP: citrate lyase; ct, citrate transporter; tct, tricarboxylate transporter. The values are the means of three biological replicates and significantly different from each other (p < 0.05) when they do not share common superscripts. the two overexpressed strains secreted more citrate out of the cell compared to the control strain, which demonstrated that both of the MT regulated the lipid synthesis by adjusting the synthesis and transport of citrate (Mizuarai et al., 2005;Gnoni et al., 2009). The excessive cytosolic citrate that cannot be utilized was exported outside the cells and resulted in the cell growth restriction and thereby affecting the final biomass (Figure 6). This is consistent with previous findings that citrate is secreted into medium when the citrate concentration is above a critical level, resulting in diminished cell growth (Makri et al., 2010;Bellou et al., 2016). There was no significant difference between the amount of extracellular citrate in these two transformants, in the stages of cell growth and rapid lipid accumulation (Figure 6). Nevertheless, when lipid accumulation became slow, the amount of citrate required was reduced leading to the accumulation of citrate in the cytosol and then secreted to the environment, since the excessed cytosolic citrate could disrupt the balance osmotic pressure of cells, resulting in impaired cell growth and biomass reduction in both transformants.
In order to explore the enzymes associated with malate secretion, citrate metabolism and fatty acid synthesis in mtand sodit-overexpressing strains, transcriptional analysis of genes encoding malic enzymes, fatty acid synthase, ATP citrate lyase, citrate transporter and TCA transporter was carried out (Figure 8). The overexpression of MT in mitochondria (mt) could promote more malate influx to mitochondria as a countersubstrate for the mitochondrial citrate transporters, thus upregulated the expression levels of these two genes (Figures 8I,J). Increasing mitochondrial malate needs more malic enzyme(s) to produce pyruvate for citrate synthesis. Therefore, the expression levels of mmes in Mc-mt are the highest in three strains. Whereas, sodit overexpression could promote malate transport out of the cell, which may drive more cytosolic citrate to produce malate and could upregulate the expression levels of genes of malic enzyme and citrate transporters. The fas1 expression was slightly decreased, however, the expression levels of acl and fas2 were increased in both strains and more acetyl-CoA was produced for fatty acid synthesis. Meanwhile, due to the expression levels of cytosolic malic enzymes (cme1 and cme2) were increased slightly, the reducing power NADPH was generated mainly for the saturated fatty acid synthesis rather than unsaturated fatty acids, which might be the explanation for the decreasing of GLA content in both transformants. Furthermore, the expression levels of ct and tct were higher in Mc-mt and Mc-sodit than that in the control, which could be the evidence for the higher citrate accumulation in cells, subsequently leading to the increased lipid production.
It is worth noting that the activities of malic enzymes decreased rapidly after lipid accumulation commences , which may cause the high cytosolic malate concentration and accelerate the amount of malate transferred into the mitochondria which used for the mitochondrial citrate synthesis, and the excessive citrate was moved to the cytosol by citrate transporters and transferred malate as the countersubstrate at the same time, thus, promoting the synthesis of lipid in Mc-mt (Huypens et al., 2011). Nevertheless, the overexpression of sodit could secrete the cytosolic malate outside of the cell during the initiation state of lipid accumulation, leads to promote the metabolism of citrate flux toward to malate biosynthesis, which provided more acetyl-CoA supply for fatty acid synthesis and promoted lipid accumulation. But, after the lipid accumulation slowed down, the continued malate transportation to the media led a lower malate concentration in cytoplasm which may decrease the citrate transportation carried by citrate transporters due to the lack of exchange substrate (Yang W. et al., 2021). Therefore, the accumulated citrate was ejected into the environment of Mc-sodit which could explain the low utilization rates of glucose and biomass in soditoverexpressing strain.

CONCLUSION
The mt and sodit genes of the high-lipid producing strain WJ11 were heterologous expressed in CBS 277.49. The levels of mRNA for mt and sodit were increased significantly during the lipid accumulation, resulted in the improvement of the lipid production in mt-and sodit-overexpressing strains. These results indicated that these two MTs have a prominent effect on malate transportation and played a vital role in lipid accumulation in CBS 277.49.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.