Skip to main content

Low dissolved oxygen supply functions as a global regulator of the growth and metabolism of Aurantiochytrium sp. PKU#Mn16 in the early stages of docosahexaenoic acid fermentation



Thraustochytrids accumulate lipids with a high content of docosahexaenoic acid (DHA). Although their growth and DHA content are significantly affected by the dissolved oxygen (DO) supply, the role of DO on the transcriptional regulation of metabolism and accumulation of intracellular metabolites remains poorly understood. Here we investigate the effects of three different DO supply conditions (10%, 30%, and 50%) on the fed-batch culture of the Aurantiochytrium PKU#Mn16 strain to mainly reveal the differential gene expressions and metabolite profiles.


While the supply of 10% DO significantly reduced the rates of biomass and DHA production in the early stages of fermentation, it achieved the highest amounts of biomass (56.7 g/L) and DHA (6.0 g/L) on prolonged fermentation. The transcriptome analyses of the early stage (24 h) of fermentation revealed several genes involved in the central carbon, amino acid, and fatty acid metabolism, which were significantly downregulated at a 10% DO level. The comparative metabolomics results revealed the accumulation of several long-chain fatty acids, amino acids, and other metabolites, supporting the transcriptional regulation under the influence of a low oxygen supply condition. In addition, certain genes involved in antioxidative systems were downregulated under 10% DO level, suggesting a lesser generation of reactive oxygen species that lead to oxidative damage and fatty acid oxidation.


The findings of this study suggest that despite the slow growth and metabolism in the early stage of fermentation of Aurantiochytrium sp. PKU#Mn16, a constant supply of low dissolved oxygen can yield biomass and DHA content better than that with high oxygen supply conditions. The critical information gained in this study will help to further improve DHA production through bioprocess engineering strategies.


Docosahexaenoic acid (DHA, C22:6) is an omega-3 long-chain polyunsaturated fatty acid (LC-PUFA), which has been shown to provide diverse human health benefits and is a well-recognized nutraceutical [1]. To meet the increasing demand for DHA and overcome the negative impacts of traditional resources, it is essential to establish a sustainable way for DHA production [2]. Like many other oleaginous microorganisms [3, 4], thraustochytrids—a group of unicellular marine heterotrophic protists—have attracted biotechnological interests because of their potential to accumulate high levels of DHA [1, 5, 6]. Interestingly, thraustochytrids have become one of the most potential microbial sources for DHA production at the industrial scale [6, 7]. However, the industrial production of DHA requires the identification and application of optimal bioprocess conditions to achieve the highest output [1, 8].

Dissolved oxygen (DO) is an important environmental condition influencing the growth and production of thraustochytrids [9,10,11]. Previous studies on thraustochytrids reported that a medium level of oxygen supply can yield the highest DHA productivity [12, 13]. For instance, an aeration rate of 1.0 vvm to Schizochytrium sp. HX-308 culture yielded the highest DHA productivity (206.7 mgL−1 h−1) compared with 0.5 vvm and 1.5 vvm [12]. With medium agitation speed (450 rpm) the Schizochytrium sp. 31 culture showed the highest concentration (21.26 g/L) of DHA, compared to that with 300 or 600 rpm [13]. Conversely, the highest volumetric oxygen transfer coefficient (kLa) of 1802 ± 105 h−1 achieved a DHA concentration of 28.93 g/L and productivity of 301 mgL−1 h−1, which were the best among those achieved with lower KLa values [14]. Furthermore, with an intermittent oxygen feeding method to maintain a 50% DO level for Aurantiochytrium limacinum SR21, a DHA concentration of 20.3 g/L was achieved [15]. However, these previous studies used different methodologies to control the oxygen levels, therefore it is somewhat difficult to compare their outputs and the effects of DO on DHA fermentation. In addition, most of the earlier investigations focused on the effects of various fermentation strategies for controlling the DO levels to maximize DHA production [10, 11, 15,16,17,18] and minimize the oxidative stress [9, 19, 20]. To determine more intuitively the effects of oxygen on the growth and metabolism of thraustochytrids, it is imperative to understand the dynamics of biomass and lipid accumulation under different DO levels. Besides, the transcriptional and metabolic changes within thraustochytrid cells in response to different DO levels remain poorly understood.

Two distinct pathways, the fatty acid synthase (FAS) pathway and PUFA synthase (PKS), are known to involve in the de novo biosynthesis of DHA in thraustochytrids [21, 22]. The former is an aerobic pathway and carries out a series of desaturation and elongation steps to introduce double bonds and extend the acyl chain. The other is an anaerobic pathway that synthesizes DHA from malonyl-CoA through the PKS-like machinery. Whether and how the DO levels affect the expression of biosynthetic genes involved in these two pathways, particularly, the ‘anaerobic’ pathway remain largely unknown. Comparative transcriptomics analysis of Schizochytrium sp. HX-308 culture between agitation of 300 rpm and 400 rpm suggested upregulation of genes associated with central carbon and fatty acid metabolism under 400 rpm agitation [23]. Particularly, a higher oxygen supply resulted in more NADPH and acetyl‑CoA production for cell growth and lipid synthesis. Furthermore, comparative metabolomics analysis of Schizochytrium sp. HX-308 culture indicated that a low oxygen level can yield significantly higher levels of succinate and several amino acids, suggesting their important roles in resisting oxygen deficiency and regulating DHA synthesis [24]. Interestingly, enhanced cell growth and PUFA production has been found to closely associate with the amino acid levels in engineered Schizochytrium sp. HX-RS [25]. Overall, these previous studies focused on the responses of thraustochytrids either at the transcriptional or metabolic level and did not provide a holistic view of cellular metabolism under the influence of different DO supply conditions.

In this study, the growth and DHA production of a previously isolated thraustochytrid strain (Aurantiochytrium sp. PKU#Mn16) were investigated in variable volume fed-batch fermentation at three different DO levels. The cellular responses to the different DO levels were further evaluated using comparative transcriptomics and metabolomics approaches. This study aimed to provide a comprehensive report on the impact of DO levels on growth, DHA production, and relevant metabolic pathways.


Growth and lipid production under different DO levels

To study the effects of DO on growth, TFA, and DHA production, variable-volume fed-batch fermentation experiments were conducted under three different DO levels (10%, 30%, and 50%). The specific growth rates (µ) of PKU#Mn16 culture cultivated under 10%, 30%, and 50% DO were 0.103 h−1, 0.111 h−1, and 0.116 h−1, respectively. At 10% oxygen saturation, the culture exhibited a comparatively delayed initiation of the exponential growth phase than those cultivated at 30% and 50% DO levels. Before the first feeding at 36 h, the biomass reached 13.2 ± 2.4 g/L, 23.9 ± 0.76 g/L, and 26.6 ± 0.34 g/L, respectively (Additional file 1: Table S1). These results indicated that 30% and 50% saturation levels can significantly increase biomass production in the early exponential phase. Furthermore, the cultures at 30% and 50% DO levels were fed twice, at 36 h and 48 h, whereas the culture under 10% level was fed once, at 48 h, during the exponential phase (Additional file 1: Fig. S1). This suggests that PKU#Mn16 culture under 30% and 50% DO levels consumed more amount of the carbon source than that under 10% DO level during exponential growth, resulting in higher biomass production. However, the culture under 50% DO level consumed a significantly lower amount of the carbon source during the 48 h–60 h period and thereafter reached the stationary growth phase (Additional file 1: Fig. S1). Interestingly, the cultures under 10% and 30% saturation reached their stationary phases much later at 96 h and meanwhile consumed more amounts of the carbon source (Fig. 1a). However, in the stationary phase, the biomass of PKU#Mn16 culture was significantly higher at 10% DO level than that at 30% and 50% DO levels. Overall, the biomass at 10%, 30%, and 50% DO levels ranged between 0 g/L–56.7 g/L, 0 g/L–47.5 g/L, and 0 g/L–38.5 g/L, respectively, during the fermentation period of 144 h (Additional file 1: Table S1). These findings suggest that although a 10% DO level can produce a higher amount of biomass at the end of fermentation, it significantly retards the growth in the exponential phase.

Fig. 1
figure 1

Dynamics of a biomass, b TFA, c DHA, and d DHA yield during fed-batch cultivation of Aurantiochytrium sp. PKU#Mn16 at 10% and 30% DO levels

The TFA productivity by PKU#Mn16 culture at 10% DO level exhibited a much slower rate (0.092 g/L.h−1) at 24 h of fermentation than those cultivated at 30% (0.275 g/L.h−1) and 50% (0.253 g/L.h−1) DO level (Additional file 1: Fig. S2). Before the first feeding at 36 h, the TFA productions at 10%, 30%, and 50% DO levels were 3.37 ± 1.05 g/L, 5.27 ± 0.78 g/L, and 7.62 ± 0.56 g/L, respectively (Additional file 1: Table S1). However, after the first feeding, the production rate of culture at 10% DO level increased gradually and remained greater than those at 30% and 50% (Additional file 1: Fig. S2). The highest rates under 10%, 30%, and 50% DO levels were observed at 60 h (0.232 g/L.h−1), 24 h (0.275 g/L.h−1), and 36 h (0.352 g/L.h−1), respectively. The maximum TFA concentrations were obtained at 144 h (19.0 ± 0.77 g/L), 120 h (13.1 ± 0.56 g/L), and 120 h (15.2 ± 0.66 g/L), respectively (Additional file 1: Table S1, Fig. 1b). However, the best TFA yields (mg/g-DCW) at 10%, 30%, and 50% saturation were achieved at 12 h (360 ± 80 mg/g-DCW), 24 h (260 ± 10 mg/g-DCW), and 72 h (380 ± 30 mg/g-DCW), respectively.

As the major PUFA produced by PKU#Mn16 (Additional file 1: Fig. S3), DHA occupied more than 30% of TFA produced under 10%, 30%, and 50% DO levels during the exponential and stationary phases. The maximum DHA concentrations were 6.0 ± 0.19 g/L (132 h), 4.7 ± 0.23 (120 h), and 5.3 ± 0.23 g/L (108 h), respectively (Additional file 1: Table S1, Fig. 1c). The DHA yields at 10% (112.5 ± 14.08 mg/g-DCW) and 30% (92.8 ± 4.20 mg/g-DCW) DO levels were the highest at 24 h and were significantly different (Fig. 2d, Additional file 1: Table S1). Whereas at 50% DO level, the highest DHA yield was at 72 h (135.3 ± 1.40 mg/g-DCW). After the first feeding at 36 h, the DHA yields at 10% and 30% DO levels did not change drastically (Fig. 1d). These findings indicate that prolonged fermentation at low DO levels can produce more amount of DHA. However, a low DO level can significantly reduce DHA concentrations in the early stages of fermentation.

Fig. 2
figure 2

Differential gene expressions in central carbon metabolism between 10 and 30% DO levels at 24 h of fermentation. Genes in red or blue are supported by at least one transcript with significant regulation (|log2FC|≥ 1 and FDR < 0.05). Red: upregulated; Blue: downregulated. pgm phosphoglucomutase, glgB 1,4-alpha-glucan branching enzyme, GPI glucose-6-phosphate isomerase, PFK 6-phosphofructokinase, ALDO fructose-bisphosphate aldolase, GAPDH glyceraldehyde 3-phosphate dehydrogenase, PGK phosphoglycerate kinase, PGAM 2,3-bisphosphoglycerate-dependent phosphoglycerate mutase, ENO enolase, PK pyruvate kinase, aceE pyruvate dehydrogenase E1 component, aceF pyruvate dehydrogenase E2 component, G6PD glucose-6-phosphate dehydrogenase, GPD glycerol-3-phosphate dehydrogenase (NAD +), SCOAS succinyl-CoA synthetase, TPI triosephosphate isomerase, HMGCS hydroxymethylglutaryl-CoA synthase, HMGCR hydroxymethylglutaryl-CoA reductase (NADPH), GPAT glycerol-3-phosphate O-acyltransferase, PLPP phosphatidate phosphatase, suhB myo-inositol monophosphatase, iolA methylmalonate-semialdehyde dehydrogenase, AGPAT 1-acyl-sn-glycerol-3-phosphate acyltransferase, DGK diacylglycerol kinase. Glu-6P Glucose 6-phosphate, PG 6-Phospho-D-glucono-1,5-lactone, G3P Glyceraldehyde 3-phosphate, PEP Phosphoenolpyruvate, Glycerol-3P Glycerol 3-phosphate, Glycerone-P Glycerone phosphate, HMG-CoA Hydroxymethylglutaroyl coenzyme A, DAG 1,2 Diacylglycerol

Overall, these findings suggest that a low DO level can markedly influence the rates of biomass and DHA production of PKU#Mn16 strain in the early stages of fermentation. Interestingly, at 10% DO level, although a contrastingly slower growth dynamics was evident, the biomass and DHA concentrations at the end of fermentation were the highest compared to that achieved at higher DO levels.

De Novo assembly and functional annotation of PKU#Mn16 transcriptome

To investigate the mechanism underlying the differential dynamics of DHA production under different DO levels, the PKU#Mn16 transcriptome was analyzed between cultures cultivated under 10% and 30% DO levels. Sequencing of cDNA libraries generated an average of 46,255,072 and 43,535,297 raw reads for 10% and 30% DO level samples, respectively. After filtering, an average of 44,774,899 and 41,958,984 clean reads, respectively, were obtained. A total of 34,723 transcripts and 29,907 unigenes were further generated after assembly (Additional file 1: Table S2). The sequence length of transcripts/unigenes varied from 201 to 13,157 bp (Additional file 1: Fig. S4). The N50 values of transcripts and unigenes were 1963 bp and 2,024 bp, respectively. Functional annotation of all the unigenes revealed that 28.35%, 20.02%, 29.18%, 25.73%, 32.95%, and 25.77% of them match with the known proteins in GO, KEGG, Pfam, SwissProt, eggNOG, and NR databases (Additional file 1: Table S2). Furthermore, a total of 7,053 unigenes were identified as DEGs, among which, 3,107 were upregulated and 3,946 were downregulated (Additional file 1: Table S2).

Transcriptional regulation of central carbon and amino acid metabolism

The comparative transcriptomics results showed that certain genes involved in the central carbon metabolism, including glycolysis, glycogen metabolism, oxidative phosphorylation, and TCA cycle were differentially regulated between 10 and 30% DO levels (Fig. 2). The DEGs encoding glycolytic enzymes included glucose-6-phosphate isomerase, phosphofructokinase, aldolase, glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate kinase, phosphoglycerate mutase, enolase, and pyruvate kinase. These DEGs, except phosphoglycerate mutase, were downregulated at a 10% DO level, suggesting that low DO can reduce glycolytic flux which may limit pyruvate formation. In addition, the genes encoding pyruvate dehydrogenase, methylmalonate semialdehyde dehydrogenase, hydroxymethylglutaryl CoA synthase (HMGCS), and hydroxymethylglutaryl CoA reductase (HMGCR) were downregulated under 10% DO level. The lower expression of these genes suggests that a low DO level can suppress the conversion of pyruvate to other precursor molecules such as mevalonate and succinyl-CoA.

The genes encoding phosphoglucomutase, 1,4 α-glucan branching enzyme, and glucose-6-phosphate dehydrogenase exhibited lower expression levels under 10% DO. Their lower expression further suggests a reduced glycolytic flux because glucose-6-phosphate (glycolytic intermediate) is one of the substrates in their reactions. The lower gene expression level of 1,4 α-glucan branching enzyme also indicated limited glycogen biosynthesis. Furthermore, the expression of the gene encoding myo-inositol monophosphatase, in the pathway of my-inositol formation from glucose-6-phosphate, was found to be lower at a 10% DO level. The genes encoding triosephosphate isomerase, glycerol-3-phosphate dehydrogenase, and glycerol-3-phosphate O-acyltransferase, involved in glycerol metabolism also exhibited lower expressions at 10% DO level.

Conversely, certain genes involved in the glycerolipid metabolism, including 1-acyl-sn-glycerol-3-phosphate acyltransferase, phosphatidate phosphatase, and diacylglycerol kinase showed higher expression under 10% DO. Similarly, the gene encoding succinyl CoA synthetase, which catalyzes the reversible reaction for substrate-level phosphorylation in the TCA cycle, also exhibited higher expression. The upregulation of this gene possibly suggests a higher flux of succinate to succinyl CoA and ADP formation during growth under a lower oxygen supply. Overall, the above results suggest that a low DO supply induces transcriptional regulation of central carbon metabolism, which may lead to the reduced flux of precursors needed for biomass formation.

Furthermore, the transcriptomics results revealed upregulation of certain genes involved in amino acid metabolism at 10% DO level (Fig. 3). These genes were found to be particularly involved in the formation of amino acids such as threonine, glycine, serine, tryptophan, alanine, lysine, arginine, proline, leucine, and isoleucine from aspartate. In addition, the argI gene encoding ornithine carbamoyltransferase, an enzyme catalyzing the conversion of glutamate to arginine, showed an increased expression. Interestingly, the genes encoding branched-chain amino acid aminotransferase and hydroxymethylglutaryl-CoA lyase, involved in the acetyl-CoA biosynthesis from amino acids, exhibited higher expressions under 10% DO. Taken together, these results suggest that PKU#Mn16 cells mainly derive acetyl-CoA from amino acids when central carbon metabolism gets negatively impacted by low oxygen supply.

Fig. 3
figure 3

Differential gene expressions in amino acid metabolism between 10 and 30% DO levels at 24 h of fermentation. Genes in red or blue are supported by at least one transcript with significant regulation (|log2FC|≥ 1 and FDR < 0.05). Red: upregulated; Blue: downregulated. trpA tryptophan synthase alpha chain, SHMT glycine hydroxymethyltransferase, GCAT glycine C-acetyltransferase, lysC aspartate kinase, hom homoserine dehydrogenase, ALT alanine transaminase, argl ornithine carbamoyltransferase, OAT ornithine–oxo-acid transaminase, proC pyrroline-5-carboxylate reductase, LEU1 3-isopropylmalate dehydratase, ilvD dihydroxy-acid dehydratase, ilvE branched-chain amino acid aminotransferase, IMDH 3-isopropylmalate dehydrogenase, hmgL hydroxymethylglutaryl-CoA lyase, P4HA prolyl 4-hydroxylase

Transcriptional regulation of fatty acid metabolism

The fermentation results indicated that a 10% DO level can significantly decrease the TFA and DHA accumulation in PKU#Mn16 strain. To further understand the corresponding transcriptional regulation of fatty acid biosynthesis, the expression levels of FAS and PKS pathway genes at 10% DO were compared with those at 30% DO (Fig. 4). The comparative results revealed the downregulation of key genes involved in the FAS and PKS pathways at 10% DO level. These genes included acetyl-CoA carboxylase (ACC), ketoacyl synthase (KS), and ketoacyl reductase (KR). ACC catalyzes the carboxylation of acetyl-CoA to malonyl-CoA and therefore its downregulation suggests that the flux of acetyl-CoA to malonyl-CoA can significantly drop when oxygen supply is low. Similarly, the downregulation of genes encoding KS, KR, fatty acid synthase, and dehydratase indicated that low oxygen supply can negatively affect the formation of long-chain acyl-CoA. Furthermore, very-long-chain 3-ketoacyl-CoA synthase and very-long-chain 3-ketoacyl-CoA reductase were also downregulated, which suggests that a low oxygen supply can reduce the flux of saturated fatty acids into DHA.

Fig. 4
figure 4

Differential gene expressions in fatty acid metabolism between 10 and 30% DO levels at 24 h of fermentation. Genes in blue are supported by at least one transcript with significant downregulation (|log2FC|≥ 1 and FDR < 0.05). ACC acetyl-CoA carboxylase, KS ketoacyl synthase, KR ketoacyl reductase, FAS2 fatty acid synthase 2, ELO very-long-chain 3-ketoacyl-CoA synthase, KAR very-long-chain 3-ketoacyl-CoA reductase, DH dehydratase, MAT malonyl-CoA ACP acyltransferase, ER enoyl reductase, ACADSB short-chain 2-methyl acyl-CoA dehydrogenase, HADH enoyl-CoA hydratase, fadA acetyl-CoA acyltransferase, SA Stearic acid, ALA α-linolenic acid, SDA Stearidonic acid, ERA Eicosatrienoic acid, ETA Eicosatetraenoic acid, EPA Eicosapentaenoic acid, DPA Docosapentaenoic acid, DHA Docosahexaenoic acid

In addition to fatty acid biosynthetic genes, the comparative transcriptomics results showed that certain genes involved in fatty acid oxidation were downregulated at a 10% DO level. These DEGs encode short-chain 2-methyl acyl-CoA dehydrogenase, enoyl-CoA hydratase, and acetyl-CoA acyltransferase, which are involved in the β-oxidation of palmitoyl-CoA to acetyl-CoA (Fig. 4). This finding suggests that low oxygen supply can reduce the flux of fatty acids into acetyl-CoA and thereby affect the TCA cycle. Overall, the comparative transcriptomic results indicate that low oxygen supply to PKU#Mn16 cells can significantly affect the transcriptional regulation of fatty acid biosynthesis and oxidation.

Transcriptional regulation of antioxidative systems

To reveal the effect of low oxygen supply on the antioxidative systems of PKU#Mn16 strain, the expressions of genes involved in enzymatic and non-enzymatic antioxidant defense mechanisms were compared between 10 and 30% DO levels (Table 1). The comparative result showed that genes encoding superoxide dismutase and ascorbate peroxidase in the enzymatic antioxidant system were downregulated at a 10% DO level. Similarly, in the non-enzymatic antioxidant system, the genes encoding HMGCS and HMGCR of the terpenoid biosynthesis pathway, inositol-phosphate phosphatase and galactose dehydrogenase of ascorbate biosynthesis pathway, and carotenoid epsilon hydroxylase of carotenoid biosynthesis pathway were downregulated. These results suggest that low oxygen supply alleviates antioxidative systems through transcriptional regulation.

Table 1 Fold change of differentially-expressed genes (|log2FC|≥ 1 and FDR < 0.05) in antioxidative systems

Intracellular metabolite profiles

To confirm the effects of DO levels on the metabolism of PKU#Mn16, the intracellular metabolites accumulated under 10% and 30% DO levels were compared using the metabolomics approach. The comparative results revealed several long-chain fatty acids such as DHA, docosapentaenoic acid (DPA), dodecanoic acid (LA), tetradecanoic acid (MA), and hexadecenoic acid (PA) that exhibited lower accumulation at 10% DO level (Fig. 5a). The lower accumulation of these fatty acids supports the downregulation of FAS and PKS pathway genes evident in this study. In addition, certain amino acids like aspartate, glutamate, histidine, valine, isoleucine, and tryptophan showed higher accumulation at 10% DO level (Fig. 5b), which was also correlated with the upregulated genes of the amino acid biosynthesis pathway. Furthermore, certain other metabolites, namely citrate, 2-oxoglutarate, AMP, ADP, glycerol, and gluconate-6-phosphate showed more accumulation at 10% DO level (Fig. 5c). As these metabolites are intermediates of central carbon metabolism, their higher accumulation suggests decreased carbon flow through glycolysis, the TCA cycle, and other pathways of glucose and glycerol metabolism, which generates reducing equivalents and metabolite precursors for biomass growth.

Fig. 5
figure 5

Differential metabolites between 10 and 30% DO levels at 24 h of fermentation. The bar indicates the fold change (log2FC) of a fatty acids, b amino acids, and c other metabolites. Metabolites with |log2FC|≥ 1 and FDR < 0.05 are shown. DHA Docosahexaenoic acid (C22:6), DPA docosapentaenoic acid (C22:5), LA Dodecanoic acid (C12:0), MA Tetradecanoic acid (C14:0), PA Hexadecenoic acid (C16:1)


Impact of dissolved oxygen on DHA fermentation

DO has been reported to play a crucial role in regulating cell growth and lipid accumulation in thraustochytrids [19]. Previous studies suggested that oxygen limitation strategies can result in increased DHA percentage but decreased biomass productivity [13, 26, 27]. In this study, three different DO levels, namely 10%, 30%, and 50% were evaluated on PKU#Mn16 culture to reveal the oxygen-dependent dynamics of cell growth and DHA production. The growth dynamics of the PKU#Mn16 strain observed in our study showed that DO levels can significantly affect the growth rate. A low DO level (10%) can markedly reduce the growth rate in the early stages of fermentation, resulting in less biomass production than 30% and 50% DO levels (Additional file 1: Table S1). However, on prolonged fermentation, until the cells reached their stationary growth phase, a 10% DO level yielded more biomass than that of higher DO levels. A similar pattern of biomass accumulation was reported when Schizochytrium sp. S31 was cultivated in fed-batch mode at three different oxygen uptake rates [13]. The higher biomass accumulation at the 10% DO level observed in our study can be explained by the fact that PKU#Mn16 culture consumed more amount of the mixed carbon source (glucose and glycerol) compared with that consumed at higher DO levels (Additional file 1: Fig. S1). Our findings suggest that a 10% DO level is optimal for achieving high biomass accumulation, which could be a result of better consumption of the carbon source by PKU#Mn16 strain in the stationary growth phase.

Previous studies have suggested that stepwise aeration strategies for the cultivation of thraustochytrids can yield the best DHA concentrations compared with that in single-stage cultivation [11, 17]. In this study, the DHA fermentation results indicated that single-stage cultivation of PKU#Mn16 strain with a constant supply of 10% DO level for a period of 144 h can achieve a high amount of DHA. However, with higher DO levels (30% and 50%), the DHA concentrations were relatively lower (Additional file 1: Table S1), which indicates that these high levels of DO saturation provided an oxygen-rich environment that either led to the oxidation of DHA or activated the FAS pathway that is considered to synthesize short-chain fatty acids [28]. Overall, our findings suggest that a constant supply of 10% DO level can be an efficient strategy to improve the overall growth and metabolism of PKU#Mn16 strain and prevent oxidative stress, ultimately resulting in better DHA production.

Role of dissolved oxygen in carbon and amino acid metabolism

Central carbon metabolism uses a complex series of enzymatic steps to convert sugars into metabolic precursors that ultimately generate the entire cell biomass [29]. A previous study reported that high oxygen supply conditions can significantly upregulate certain genes related to central carbon metabolism, including glycolysis, the pentose phosphate pathway, and the TCA cycle [23]. In this study, the results of comparative transcriptomics revealed transcriptional regulation of central carbon metabolism, suggesting reduced fluxes of metabolic precursors needed for biomass formation at 10% DO level (Fig. 2). This was further supported by the reduced levels of certain precursor metabolites generated in the central carbon metabolism (Fig. 5c). The accumulation of citrate, 2-oxoglutarate, AMP, ADP, glycerol, and gluconate-6-phosphate possibly suggests decreased carbon flow through glycolysis, the TCA cycle, and other pathways of glucose and glycerol metabolism. As these pathways of the central carbon metabolism network generate reducing equivalents and metabolic precursors for biomass formation, their downregulation concurs with the reduced biomass formation in the early stages of fermentation at a 10% DO level.

Apart from the regulation of central carbon metabolism, our study revealed that amino acid metabolism can also get regulated when a 10% DO level is provided to PKU#Mn16 cells. Particularly, the upregulation of ALT, ilvE, and hmgL genes suggests increased conversion of alanine, valine, leucine, and isoleucine to acetyl-CoA. It can thus be speculated that these amino acids can be the main sources of acetyl-CoA at low DO levels. Furthermore, previous studies have reported that some of these amino acids can also play an important role in tolerance and cell viability when microbes combat environmental stress [30]. According to previous studies, tryptophan biosynthetic genes are related to ethanol tolerance of Saccharomyces cerevisiae [31]. The increased supplementation of histidine, leucine, and alanine was proposed to improve the stress tolerance of Saccharomyces cerevisiae against an excess glucose concentration [32]. When Schizochytrium sp. was exposed to an oxygen limitation condition, alanine acted as a protective molecule to cope with the variable conditions due to stress responses [24]. It has been suggested that alanine can be an essential regulatory substance for thraustochytrids to resist the environmental stress of oxygen deficiency. Similarly, glutamate has been reported to play essential roles in cellular metabolism in microorganisms [33]. Glutamine is also found to act as an oxygen stress protectant and can transform into glutamate [24]. In this study, the glutamate accumulation at 10% DO level was significantly increased compared with that at 30% DO level (Fig. 5b). Thus, this result suggests the role of glutamate in stress response against oxygen deficiency during early stages of fermentation. As abiotic stress often leads to the accumulation of amino acids [34, 35], our findings support the fact that certain amino acids can serve as potential stress mitigators of fatty acid biosynthesis in thraustochytrids under low oxygen supply conditions.

Effects of dissolved oxygen on fatty acid biosynthesis

The analysis of PKU#Mn16 transcriptome revealed several genes encoding fatty acid desaturases involved in the FAS pathway (Fig. 4). However, the critical Δ9, Δ12, Δ15, Δ6, and Δ5 desaturase genes essential for DHA synthesis were not found, suggesting incomplete FAS pathway in this strain. This was in agreement with previous reports [36, 37]. Interestingly, our study found genes related to the PKS pathway, which revealed that the PKS pathway played a more important role in PKU#Mn16 strain. Similarly, the transcriptome profiling of Aurantiochytrium sp. PKU#SW7 revealed that DHA biosynthesis occurs via the PKS pathway rather than the FAS pathway [36].

Molecular oxygen is required for fatty acid biosynthesis via the FAS pathway [22]. Whereas the PKS pathway only utilizes the PKS gene cluster to synthesize PUFAs from acyl-CoA de novo [38,39,40]. A high oxygen supply condition has been shown to upregulate the FAS gene expression in Schizochytrium sp. HX-308 [23], suggesting the role of DO in the FAS pathway. In our study, the comparative transcriptomics and metabolomics results indicated that low oxygen supply to PKU#Mn16 cells can significantly affect both fatty acid biosynthesis and oxidation (Fig. 4). ACC, a key enzyme in the fatty acid biosynthesis, produces malonyl-CoA that is used in both FAS and PKS pathways [41, 42]. Our results showed the downregulation of the ACC gene, which indicates that malonyl-CoA supply to these pathways may reduce at a 10% DO level, resulting in decreased flux of malonyl-ACP to DHA. The downregulation of KS, KR, and DH genes, and decreased accumulation of several fatty acids, including DHA (Fig. 5a), further support this inference.

Although an optimal oxygen supply condition improves lipid accumulation, a high level of DO can lead to lipid peroxidation [43]. This study found that the fatty acid oxidation pathway was downregulated at the level of 10% DO (Fig. 4), suggesting that a low oxygen supply condition is appropriate for limiting lipid peroxidation. Furthermore, for oleaginous microorganisms such as Schizochytrium, a high oxygen supply can result in its conversion into reactive oxygen species (ROS). Interestingly, it was found that some genes involved in antioxidative systems were downregulated at a 10% DO level, which suggests that ROS generation was less at low oxygen supply conditions. A low ROS generation indirectly supports the higher DHA accumulation at the 10% DO level observed in our study.


Central carbon, amino acid, and fatty acid metabolism are significantly downregulated in the early fermentation stage of Aurantiochytrium sp. PKU#Mn16 at low oxygen supply conditions (10% DO). This study revealed that a 10% DO supply to PKU#Mn16 culture can achieve the highest biomass and DHA production at the end of fermentation compared with that achieved at higher DO levels (30% and 50%). However, it was noted that in the early stages of fermentation, the biomass and DHA production could be significantly lower with a 10% DO level than that with higher DO levels. The findings of this study suggest that a constant supply of a low DO level in a variable volume fed-batch fermentation can be an efficient strategy to achieve optimal growth and DHA production.



A previously isolated thraustochytrid [44], Aurantiochytrium limacinum strain PKU#Mn16, was used in the present study. The PKU#Mn16 strain was subcultured every two weeks on M4 solid medium (glucose, 20 g/L; peptone, 1.5 g/L; yeast extract, 1 g/L; sea salt, 33 g/L; and agar, 20 g/L). The seed culture was prepared by cultivating a single colony of PKU#Mn16 from the solid M4 medium into flasks (100 mL) containing 30 mL of M4 medium at 28 °C on an orbital shaker at 170 rpm for 24 h. The 24 h grown seed culture (10%, v/v) was then transferred into flasks (1000 mL) containing 300 mL of M4 medium and cultivated at 28 °C on an orbital shaker at 170 rpm for 24 h. The resulting seed culture was used as the inoculum (10%, v/v) for further fermentation experiments.

Fermentation experiments

Experiments were conducted in a 5 L fermenter (model: SY-9000-V9, Shanghai DongMing Industrial Co. Ltd., Shanghai, China) equipped with DO and pH electrodes, a temperature sensor, an impeller, and an air pump. The working volume of the fermenter was 3 L. The aeration rate was maintained at 300 L/h and three different DO saturation levels, i.e., 10%, 30%, and 50%, were individually studied by altering the agitation speed. A 100% DO saturation level was achieved by setting the agitation at 800 rpm and aeration rate at 500 L/h. The seed culture (300 mL), as described in “Microorganism” section, was inoculated into the fermentation medium (glucose, 30 g/L; glycerol, 30 g/L; yeast extract, 15 g/L; KH2PO4, 0.25 g/L; sea salt, 33 g/L; and defoamer, 1 mL/L). The fermenter temperature was maintained at 28 ℃ throughout the fermentation. A stock solution (800 g/L) containing equal amounts of glucose and glycerol was prepared and used as a feed to replenish the substrate concentration in the fermenter. About 80–150 g/L of this feed was added to the fermentation broth when the substrate concentration dropped below 25 g/L (Additional file 1: Fig. S1). Samples were taken from the fermenter every 12 h for the analyses of residual carbon source, dry cell weight (DCW), and lipid content. The specific growth rate (µ) of the PKU#Mn16 culture under different DO levels was calculated using Eq. (1) [45].

$${B}_{t}={B}_{0}{e}^{\mu t}$$

where Bt is the DCW at a given time t, B0 is the DCW at t = 0, and µ is the specific growth rate.

Biochemical analyses

DCW, TFA, and fatty acid composition were determined following methods described in our previous study [46]. The concentrations of glucose and glycerol were determined by enzymatic method using commercial kits (Applygen Technologies Inc., Beijing, China).

Transcriptome profiling

To investigate the cell response to oxygen saturation, the transcriptional regulation in PKU#Mn16 strain between 10 and 30% saturation levels was analyzed. Three samples (15 mL each) were collected at 24 h of fermentation and analyzed in parallel for each saturation level. The cells were harvested from each sample, frozen with liquid nitrogen, and stored at −80 °C for RNA sequencing. Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's instructions. The RNA fragments were reverse-transcribed to create the cDNA library following the mRNA-Seq Sample Preparation Kit (Illumina, San Diego, USA) protocol. The average insert size for the paired-end libraries was 300 bp (± 50 bp). Total RNA extraction, quality, integrity, and quantity checks were conducted by Lianchuan Bio (Hangzhou, China).

The paired-end sequencing was performed on an Illumina NovaSeq™ 6000 platform (Illumina Inc., USA) at Lianchuan Bio (Hangzhou, China). Cutadapt [47] and in-house Perl scripts were used to remove the reads that contained adaptor contamination, low-quality bases, and undetermined bases. Then, the quality of the clean sequencing reads was verified using the FastQC tool (, including the Q20, Q30, and GC content. All downstream analyses were based on high-quality clean data. De novo assembly of the transcriptome was performed with Trinity 2.4.0 [48]. All assembled unigenes were aligned against the NCBI non-redundant (NR) protein database, Gene ontology (GO), SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), Pfam, and eggNOG databases using DIAMOND [49] with a threshold E value < 0.00001. Salmon [50] tool was used to quantify the expression of the transcripts by calculating TPM [51]. The differentially-expressed genes (DEGs) were selected with log2 (fold change) > 1 or log2 (fold change) < −1 and with statistical significance at alpha = 0.05 by R package edgeR [52]. The raw sequencing reads are deposited in NCBI under the BioProject PRJNA901145.

Metabolome analysis

The changes at the level of intracellular metabolites of PKU#Mn16 strain, in response to different (10% and 30%) DO levels, were evaluated by extracting the metabolites and analyzing them by mass spectrometry. Six samples (total OD660 of each sample = 24) were collected at 24 h of fermentation and analyzed in parallel for each DO level. The samples were centrifuged at 10,000 g for 5 min at −4 °C, the supernatant was discarded, and the cell pellets were washed thrice with 20 mL of PBS buffer (pH ~ 7.4). The resulting cell pellets were quickly frozen in liquid nitrogen and then stored at −80 °C until further processing at Lianchuan Bio (Hangzhou, China).

To extract the metabolites, the frozen cell pellets were thawed on ice, and 20 µL of each was quenched with 120 µL of precooled methanol (50% v/v), and then centrifugated at 4000g for 20 min. The supernatant was stored at −80 °C before the LC‐MS analysis. All samples were analyzed using a TripleTOF 5600 Plus high-resolution tandem mass spectrometer (SCIEX, Warrington, UK). Chromatographic separation was performed using an ultra-performance liquid chromatography (UPLC) system (SCIEX, UK). An ACQUITY UPLC T3 column (100 mm × 2.1 mm, 1.8 µm, Waters, UK) was used for the reversed-phase separation. The TripleTOF 5600 Plus system was used to detect the metabolites eluted from the column.

The open-access databases, KEGG and HMDB, were used to annotate the metabolites by matching the exact molecular mass data (m/z) to those from the database within a threshold of 10 ppm. Data normalization was performed on all samples using the probabilistic quotient normalization algorithm. The P value was analyzed by Student’s t‐test and adjusted for multiple tests using FDR (Benjamini–Hochberg).

Statistical analysis

The mean and standard deviation for each measured parameter and the test of significance (Student’s t-test) were computed in Microsoft Excel (Microsoft Corporation, 2016).

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.


  1. Patel A, Karageorgou D, Katapodis P, Sharma A, Rova U, Christakopoulos P, Matsakas L. Bioprospecting of thraustochytrids for omega-3 fatty acids: a sustainable approach to reduce dependency on animal sources. Trends Food Sci Technol. 2021;115:433–44.

    Article  CAS  Google Scholar 

  2. Orozco Colonia BS, de Melo Pereira GV, Soccol CR. Omega-3 microbial oils from marine thraustochytrids as a sustainable and technological solution: a review and patent landscape. Trends Food Sci Technol. 2020;99:244–56.

    Article  CAS  Google Scholar 

  3. Lu Q, Li H, Xiao Y, Liu H. A state-of-the-art review on the synthetic mechanisms, production technologies, and practical application of polyunsaturated fatty acids from microalgae. Algal Res. 2021;55: 102281.

    Article  Google Scholar 

  4. Adarme-Vega TC, Lim DKY, Timmins M, Vernen F, Li Y, Schenk PM. Microalgal biofactories: a promising approach towards sustainable omega-3 fatty acid production. Microb Cell Fact. 2012;11:96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Du F, Wang YZ, Xu YS, Shi TQ, Liu WZ, Sun XM, Huang H. Biotechnological production of lipid and terpenoid from thraustochytrids. Biotechnol Adv. 2021;48: 107725.

    Article  CAS  PubMed  Google Scholar 

  6. Gupta A, Barrow CJ, Puri M. Multiproduct biorefinery from marine thraustochytrids towards a circular bioeconomy. Trends Biotechnol. 2021;40:448–62.

    Article  PubMed  Google Scholar 

  7. Saini RK, Prasad P, Sreedhar RV, Akhilender Naidu K, Shang X, Keum YS. Omega-3 Polyunsaturated Fatty Acids (PUFAs): emerging plant and microbial sources, oxidative stability, bioavailability, and health benefits-a review. Antioxidants. 2021;10:1627.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Chi G, Xu Y, Cao X, Li Z, Cao M, Chisti Y, He N. Production of polyunsaturated fatty acids by Schizochytrium (Aurantiochytrium) spp. Biotechnol Adv. 2022;55: 107897.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang S, He Y, Sen B, Wang G. Reactive oxygen species and their applications toward enhanced lipid accumulation in oleaginous microorganisms. Bioresour Technol. 2020;307: 123234.

    Article  CAS  PubMed  Google Scholar 

  10. Sun XM, Ren LJ, Ji XJ, Chen SL, Guo DS, Huang H. Adaptive evolution of Schizochytrium sp. by continuous high oxygen stimulations to enhance docosahexaenoic acid synthesis. Bioresour Technol. 2016;211:374–81.

    Article  CAS  PubMed  Google Scholar 

  11. Chi Z, Liu Y, Frear C, Chen S. Study of a two-stage growth of DHA-producing marine algae Schizochytrium limacinum SR21 with shifting dissolved oxygen level. Appl Microbiol Biotechnol. 2009;81:1141–8.

    Article  CAS  PubMed  Google Scholar 

  12. Guo DS, Ji XJ, Ren LJ, Li GL, Yin FW, Huang H. Development of a real-time bioprocess monitoring method for docosahexaenoic acid production by Schizochytrium sp. Bioresour Technol. 2016;216:422–7.

    Article  CAS  PubMed  Google Scholar 

  13. Chang G, Wu J, Jiang C, Tian G, Wu Q, Chang M, Wang X. The relationship of oxygen uptake rate and kLa with rheological properties in high cell density cultivation of docosahexaenoic acid by Schizochytrium sp. S31. Bioresour Technol. 2014;152:234–40.

    Article  CAS  PubMed  Google Scholar 

  14. Chang G, Gao N, Tian G, Wu Q, Chang M, Wang X. Improvement of docosahexaenoic acid production on glycerol by Schizochytrium sp. S31 with constantly high oxygen transfer coefficient. Bioresour Technol. 2013;142:400–6.

    Article  CAS  PubMed  Google Scholar 

  15. Huang TY, Lu WC, Chu IM. A fermentation strategy for producing docosahexaenoic acid in Aurantiochytrium limacinum SR21 and increasing C22:6 proportions in total fatty acid. Bioresour Technol. 2012;123:8–14.

    Article  CAS  PubMed  Google Scholar 

  16. Guo DS, Ji XJ, Ren LJ, Li GL, Huang H. Improving docosahexaenoic acid production by Schizochytrium sp. using a newly designed high-oxygen-supply bioreactor. AIChE J. 2017;63:4278–86.

    Article  CAS  Google Scholar 

  17. Ren L, Ji X, Huang H, Qu L, Feng Y, Tong Q, Ouyang P. Development of a stepwise aeration control strategy for efficient docosahexaenoic acid production by Schizochytrium sp. Appl Microbiol Biotechnol. 2010;87:1649–56.

    Article  CAS  PubMed  Google Scholar 

  18. Qu L, Ji XJ, Ren LJ, Nie ZK, Feng Y, Wu WJ, Ouyang PK, Huang H. Enhancement of docosahexaenoic acid production by Schizochytrium sp. using a two-stage oxygen supply control strategy based on oxygen transfer coefficient. Lett Appl Microbiol. 2011;52:22–7.

    Article  CAS  PubMed  Google Scholar 

  19. Sun X, Geng L, Ren L, Ji X, Hao N, Chen K, Huang H. Influence of oxygen on the biosynthesis of polyunsaturated fatty acids in microalgae. Bioresour Technol. 2018;250:868–76.

    Article  CAS  PubMed  Google Scholar 

  20. Ren L, Sun X, Ji X, Chen S, Guo D, Huang H. Enhancement of docosahexaenoic acid synthesis by manipulation of antioxidant capacity and prevention of oxidative damage in Schizochytrium sp. Bioresour Technol. 2017;223:141–8.

    Article  CAS  PubMed  Google Scholar 

  21. Morabito C, Bournaud C, Maës C, Schuler M, Cigliano RA, Dellero Y, Maréchal E, Amato A, Rébeillé F. The lipid metabolism in thraustochytrids. Prog Lipid Res. 2019;76: 101007.

    Article  CAS  PubMed  Google Scholar 

  22. Xie Y, Wang G. Mechanisms of fatty acid synthesis in marine fungus-like protists. Appl Microbiol Biotechnol. 2015;99:8363–75.

    Article  CAS  PubMed  Google Scholar 

  23. Bi Z-Q, Ren L-J, Hu X-C, Sun X-M, Zhu S-Y, Ji X-J, Huang H. Transcriptome and gene expression analysis of docosahexaenoic acid producer Schizochytrium sp. under different oxygen supply conditions. Biotechnology Biofuels. 2018;11:249.

    Article  CAS  Google Scholar 

  24. Li J, Ren LJ, Sun GN, Qu L, Huang H. Comparative metabolomics analysis of docosahexaenoic acid fermentation processes by Schizochytrium sp. under different oxygen availability conditions. OMICS. 2013;17:269–81.

    Article  CAS  PubMed  Google Scholar 

  25. Geng L, Chen S, Sun X, Hu X, Ji X, Huang H, Ren L. Fermentation performance and metabolomic analysis of an engineered high-yield PUFA-producing strain of Schizochytrium sp. Bioprocess Biosyst Eng. 2019;42:71–81.

    Article  CAS  PubMed  Google Scholar 

  26. Jakobsen AN, Aasen IM, Josefsen KD, Strøm AR. Accumulation of docosahexaenoic acid-rich lipid in Thraustochytrid Aurantiochytrium sp. strain T66: effects of N and P starvation and O2 limitation. Appl Microbiol Biotechnol. 2008;80:297–306.

    Article  CAS  PubMed  Google Scholar 

  27. Bailey RB, Don D, Hansen JM, Mirrasoul PJ, Ruecker CM, Veeder GTI, Kaneko T, Barclay WR. Enhanced Production of Lipids Containing Polyenoic Fatty Acid by Very High Density Cultures of Eukaryotic Microbes in Fermentors. In US Patent No 6607900, 6607900. US2003

  28. Półbrat T, Konkol D, Korczyński M. Optimization of docosahexaenoic acid production by Schizochytrium sp.—a review. Biocatal Agric Biotechnol. 2021;35:102076.

    Article  Google Scholar 

  29. Noor E, Eden E, Milo R, Alon U. Central carbon metabolism as a minimal biochemical walk between precursors for biomass and energy. Mol Cell. 2010;39:809–20.

    Article  CAS  PubMed  Google Scholar 

  30. Zhao XQ, Bai FW. Mechanisms of yeast stress tolerance and its manipulation for efficient fuel ethanol production. J Biotechnol. 2009;144:23–30.

    Article  CAS  PubMed  Google Scholar 

  31. Hirasawa T, Yoshikawa K, Nakakura Y, Nagahisa K, Furusawa C, Katakura Y, Shimizu H, Shioya S. Identification of target genes conferring ethanol stress tolerance to Saccharomyces cerevisiae based on DNA microarray data analysis. J Biotechnol. 2007;131:34–44.

    Article  CAS  PubMed  Google Scholar 

  32. Pham TK, Wright PC. The proteomic response of Saccharomyces cerevisiae in very high glucose conditions with amino acid supplementation. J Proteome Res. 2008;7:4766.

    Article  CAS  PubMed  Google Scholar 

  33. Chen G, Fan KW, Lu FP, Li Q, Aki T, Chen F, Jiang Y. Optimization of nitrogen source for enhanced production of squalene from thraustochytrid Aurantiochytrium sp. New Biotechnol. 2010;27:382–9.

    Article  CAS  Google Scholar 

  34. Maroli AS, Nandula VK, Duke SO, Tharayil N. Stable isotope resolved metabolomics reveals the role of anabolic and catabolic processes in glyphosate-induced amino acid accumulation in Amaranthus palmeri biotypes. J Agric Food Chem. 2016;64:7040–8.

    Article  CAS  PubMed  Google Scholar 

  35. Chen H, Wang Q. Regulatory mechanisms of lipid biosynthesis in microalgae. Biolog Rev. 2021;96:2373.

    Article  CAS  Google Scholar 

  36. Liang Y, Liu Y, Tang J, Ma J, Cheng JJ, Daroch M. Transcriptomic profiling and gene disruption revealed that two genes related to PUFAs/DHA biosynthesis may be essential for cell growth of Aurantiochytrium sp. Mar Drugs. 2018;16:310.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Liu L, Hu Z, Li S, Yang H, Li S, Lv C, Zaynab M, Cheng CHK, Chen H, Yang X. Comparative transcriptomic analysis uncovers genes responsible for the DHA enhancement in the mutant Aurantiochytrium sp. Microorganisms. 2020;8:529.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Metz JG, Roessler P, Facciotti D, Levering C, Dittrich F, Lassner M, Valentine R, Lardizabal K, Domergue F, Yamada A, et al. Production of polyunsaturated fatty acids by polyketide synthases in both prokaryotes and eukaryotes. Science. 2001;293:290–3.

    Article  CAS  PubMed  Google Scholar 

  39. Hauvermale A, Kuner J, Rosenzweig B, Guerra D, Diltz S, Metz JG. Fatty acid production in Schizochytrium sp.: involvement of a polyunsaturated fatty acid synthase and a type I fatty acid synthase. Lipids. 2006;41:739–47.

    Article  CAS  PubMed  Google Scholar 

  40. Lippmeier JC, Crawford KS, Owen CB, Rivas AA, Metz JG, Apt KE. Characterization of both polyunsaturated fatty acid biosynthetic pathways in Schizochytrium sp. Lipids. 2009;44:621–30.

    Article  CAS  PubMed  Google Scholar 

  41. Zu XY, Zhong J, Luo DX, Tan JJ, Zhang QH, Wu Y, Liu JH, Cao RX, Wen GB, Cao DL. Chemical genetics of acetyl-CoA carboxylases. Molecules. 2013;18:1704–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Chen D, Yuan X, Liang LM, Liu K, Ye HY, Liu ZY, Liu YF, Huang LQ, He WJ, Chen YQ, et al. Overexpression of acetyl-CoA carboxylase increases fatty acid production in the green alga Chlamydomonas reinhardtii. Biotechnol Lett. 2019;41:1133–45.

    Article  CAS  PubMed  Google Scholar 

  43. Guichardant M, Chen P, Liu M, Calzada C, Colas R, Véricel E, Lagarde M. Functional lipidomics of oxidized products from polyunsaturated fatty acids. Chem Phys Lipid. 2011;164:544–8.

    Article  CAS  Google Scholar 

  44. Liu Y, Singh P, Sun Y, Luan S, Wang G. Culturable diversity and biochemical features of thraustochytrids from coastal waters of Southern China. Appl Microbiol Biotechnol. 2014;98:3241–55.

    Article  CAS  PubMed  Google Scholar 

  45. Kirchman DL. Calculating microbial growth rates from data on production and standing stocks. Mar Ecol Prog Ser. 2002;233:303–6.

    Article  Google Scholar 

  46. Wang Q, Sen B, Liu X, He Y, Xie Y, Wang G. Enhanced saturated fatty acids accumulation in cultures of newly-isolated strains of Schizochytrium sp. and Thraustochytriidae sp. for large-scale biodiesel production. Sci Total Environ. 2018;631–632:994–1004.

    Article  PubMed  Google Scholar 

  47. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.

    Article  Google Scholar 

  48. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.

    Article  CAS  PubMed  Google Scholar 

  50. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  52. Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139.

    Article  PubMed  Google Scholar 

Download references


This research was partially funded by the National Natural Science Foundation of China (32170063).

Author information

Authors and Affiliations



LL designed the study, carried out most of the experiments and wrote the manuscript. XZ assisted in the biochemical analysis of culture samples. HY assisted in the fermentation experiments. YW participated in some growth experiments. BS assisted in the conceptualization and design of the study, and contributed to editing and revising the manuscript. GW conceived the study, contributed to editing the manuscript, and received the funding. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Biswarup Sen or Guangyi Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Comparison of the fermentation parameters of Aurantiochytrium sp. PKU#Mn16 under different dissolved oxygen levels. Table S2. Annotation summary of de novo assembled unigenes. Table S3. Differentially-expressed genes in various metabolic pathways between 10% and 30% oxygen saturation levels. Figure S1. Feeding regime for the variable-volume fed-batch culture to evaluate the effects of different oxygen saturation levels (10%, 30%, and 50%) on growth and lipid production of PKU#Mn16 strain. The arrows indicate the point of feeding. Figure S2. The profiles of TFA production rates of PKU#Mn16 culture under different oxygen saturation levels. Rate was calculated by dividing the difference in TFA concentration by the change in time between two consecutive data. Figure S3. Bar plot of major fatty acids produced by PKU#Mn16 strain under different oxygen saturation levels (10% and 30%). C15, C16, DPA, DHA, and TFA stand for pentadecanoic acid, palmitic acid, docosapentaenoic acid, docosahexaenoic acid, and total fatty acids, respectively. Figure S4. Sequence length distribution of transcripts and unigenes

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, L., Zhu, X., Ye, H. et al. Low dissolved oxygen supply functions as a global regulator of the growth and metabolism of Aurantiochytrium sp. PKU#Mn16 in the early stages of docosahexaenoic acid fermentation. Microb Cell Fact 22, 52 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Thraustochytrids
  • Transcriptomics
  • Metabolomics
  • Central carbon metabolism
  • Amino acid metabolism
  • Fatty acids