Low-temperature effects on docosahexaenoic acid biosynthesis in Schizochytrium sp. TIO01 and its proposed underlying mechanism
Biotechnology for Biofuels volume 13, Article number: 172 (2020)
Schizochytrium species are known for their abundant production of docosahexaenoic acid (DHA). Low temperatures can promote the biosynthesis of polyunsaturated fatty acids (PUFAs) in many species. This study investigates low-temperature effects on DHA biosynthesis in Schizochytrium sp. TIO01 and its underlying mechanism.
The Schizochytrium fatty acid biosynthesis pathway was evaluated based on de novo genome assembly (contig N50 = 2.86 Mb) and iTRAQ-based protein identification. Our findings revealed that desaturases, involved in DHA synthesis via the fatty acid synthase (FAS) pathway, were completely absent. The polyketide synthase (PKS) pathway and the FAS pathway are, respectively, responsible for DHA and saturated fatty acid synthesis in Schizochytrium. Analysis of fatty acid composition profiles indicates that low temperature has a significant impact on the production of DHA in Schizochytrium, increasing the DHA content from 43 to 65% of total fatty acids. However, the expression levels of PKS pathway genes were not significantly regulated as the DHA content increased. Further, gene expression analysis showed that pathways related to the production of substrates (acetyl-CoA and NADPH) for fatty acid synthesis (the branched-chain amino acid degradation pathway and the pentose phosphate pathway) and genes related to saturated fatty acid biosynthesis (the FAS pathway genes and malic enzyme) were, respectively, upregulated and downregulated. These results indicate that low temperatures increase the DHA content by likely promoting the entry of relatively large amounts of substrates into the PKS pathway.
In this study, we provide genomic, proteomic, and transcriptomic evidence for the fatty acid synthesis pathway in Schizochytrium and propose a mechanism by which low temperatures promote the accumulation of DHA in Schizochytrium. The high-quality and nearly complete genome sequence of Schizochytrium provides a valuable reference for investigating the regulation of polyunsaturated fatty acid biosynthesis and the evolutionary characteristics of Thraustochytriidae species.
Docosahexaenoic acid (DHA, C22:6), a major ω-3 polyunsaturated fatty acid (PUFA), is widely distributed among phospholipids in the human brain and retina, playing a vital role in human health . A variety of marine microorganisms are rich in DHA [2, 3]. As one of the most promising DHA-producing species, Schizochytrium has the ability to accumulate DHA at more than 40% of the total fatty acid (TFA) concentration and ~ 40% of the dry cell weight . On account of this characteristic, the PUFA biosynthesis pathway in Schizochytrium has attracted increasing attention and is currently being exploited by several companies .
The entire PUFA synthesis pathway conventionally occurs via two routes. The first route includes the standard fatty acid synthase (FAS) pathway, which involves serial desaturation and elongation of short saturated fatty acids (SSFAs, C16:0 or C18:0). The second route includes the polyketide synthase (PKS) pathway . In 2002, Metz et al. reported that the PKS pathway is responsible for PUFA synthesis and does not rely on elongation and desaturation of the FAS pathway in Schizochytrium through the use of labeling experiments . This result was subsequently validated by Lippmeier et al. who showed that knock-outs of the PKS gene led to PUFA auxotrophic behavior . Although these biochemical and genetic studies on Schizochytrium have suggested that the PKS pathway is responsible for DHA biosynthesis [7,8,9], recent studies revealed that enhancement of the FAS pathway contributes significantly to increased DHA productivity in Schizochytrium [10, 11]. Based on these results, questions arise as to which route is employed by Schizochytrium for DHA synthesis. Genome sequencing is an ideal approach for investigating this issue, however, information on high-quality whole genomes of Schizochytrium currently available in public databases is limited , preventing identification of the Schizochytrium PUFA biosynthesis pathway.
Low temperatures have been shown to impact PUFA productivity in many PUFA-producing organisms [13,14,15,16,17]. However, few studies have investigated the mechanism underlying these effects. Ma et al. pioneered in-depth studies in Aurantiochytrium, showing that the upregulation of the PKS pathway and downregulation of the FAS pathway play important roles in promoting DHA accumulation [13, 18]. Min et al. reported that significant upregulation of the fatty acid desaturases associated with the FAS pathway leads to PUFA accumulation in Bangia fuscopurpurea . In Shewanella piezotlerans and Photobacterium profundum, upregulation of PKS pathway genes was not detected despite the increase in PUFA accumulation at a low temperature [20, 21]. In Schizochytrium, Zeng et al. reported that low temperatures have a significant impact on DHA production, increasing DHA levels up to 50% of the TFA level . The mechanism underlying this low-temperature-induced DHA accumulation remains unknown in Schizochytrium.
Previously, we isolated a Schizochytrium strain from sea water (designated Schizochytrium sp. TIO01) . In this study, the effects of low temperature on fatty acid synthesis in Schizochytrium sp TIO01 were analyzed and the underlying mechanism investigated. Here, we propose a reconstructed fatty acid biosynthesis pathway of Schizochytrium sp. TIO01 based on de novo genome assembly and iTRAQ-based LC-MS-MS. The effects of low temperature on fatty acid synthesis were investigated by differential gene expression analysis.
Low temperatures promote DHA accumulation in Schizochytrium
The fatty acid composition profiles of Schizochytrium cultured under both cold and normal temperatures were analyzed. Samples used to analyze fatty acids were cultured and collected at an optical density at 650 nm of approximately 0.7 at 28 or 16 °C, as indicated in Fig. 1a. Generally, Schizochytrium sp. TIO01 had similar total lipid levels (16 °C: 53.53% ± 2.18% (weight/cell dry weight); 28 °C: 51.05% ± 3.05% (weight/cell dry weight)) and lipid content profiles at the two temperatures. Pentadecanoic acid (C15:0), palmitic acid (C16:0), heptadecanoic acid (C17:0), docosapentaenoic acid (DPA, C22:5, ω-6), and DHA (C22:6, ω-3) were all detected under both conditions. Tetracosanoic acid (C24:0) was detected only when Schizochytrium was cultured at 28 °C. At low temperatures, Schizochytrium exhibited reduced saturated fatty acid synthesis and a preference for increased DHA synthesis. DHA constituted ~ 65% of the total fatty acid content at 16 °C compared to ~ 43% at 28 °C (Fig. 1b).
Schizochytrium genome assembly, assessment and annotation
To reveal the fatty acid biosynthesis pathways and further improve our understanding of the underlying mechanism by which low temperatures promote DHA accumulation in Schizochytrium, we performed de novo whole-genome assembly and annotation. Using 13.9 Gb of PacBio RS II subreads (217 × genomic data) and 31 Gb of Illumina PE250 clean data (480 × genomic data), a genome size of 64 Mb with a 45% GC base ratio containing 34 scaffolds and 3 circular contigs (one of the circular contigs is the mtDNA genome, which was previously reported ) was obtained. The N50 scaffold and N50 contig were 5.83 Mb and 2.86 Mb, respectively. Table 1 shows detailed information regarding the assembled Schizochytrium genome and the genome for Thraustochytriaceae species. To inspect the accuracy of de novo assembly, multiple independent sources of reads were used to assess the assembled genome. Among the 62,457,386 paired Illumina PE250 reads (approximately 480 × coverage of the genome), over 99% could be aligned to the genome. The overall alignment rate of the transcriptomic reads from 14 independent culture conditions (approximately 900 × coverage of the genome) ranged from 95 to 97%, as shown in Additional file 1: Table S1). These results suggest that the assembled Schizochytrium genome was of high-quality and nearly complete. We also evaluated the completeness of the genome assembly using BUSCO-v3.0 (database: eukaryota_odb9) . The results showed that 90.4% of the core genes were detected (C: 90.4% [S: 89.1%, D: 1.3%], F: 2.3%, M: 7.3%, 274/303). Based on the integration of RNA-seq data, protein alignment, and de novo predictions, 12,392 protein-coding genes were predicted with an average transcript size of 1876.3 bp and a mean of 1.7 exons per gene. In all, 6796 out of the 12,392 predicted proteins (54.84%) could be classified into families according to their putative functions. The proportions for each database were as follows: KEGG (4362, 35.20%), SwissProt (5078, 40.98%), TrEMBL (6336, 51.13%), Nr (6063, 48.93%), GO (1235, 9.97%), and InterPro (6242, 50.37%).
iTRAQ-based protein identification
To identify Schizochytrium proteins, six samples of Schizochytrium cells were grown at both normal temperatures and low temperatures and subjected to iTRAQ-based proteomic analysis. Tandem mass spectra were searched against the Schizochytrium protein database containing the genomics-predicted proteins and transcriptomics-predicted novel proteins. A total of 4008 proteins were quantified (detailed information regarding the iTRAQ-identified proteins is provided in Additional file 2: Table S2), of which 3196 were annotated by the reference genome and 812 were predicted by transcriptomic analysis. Using KEGG pathway analysis, 2939 proteins among the 4008 proteins were annotated in 45 pathways. These identified proteins are involved in global and overview maps, signal transduction, carbohydrate metabolism, lipid metabolism, amino acid metabolism, and energy metabolism, among others. The top 20 KEGG pathways of iTRAQ-identified proteins are shown in Fig. 2a.
Fatty acid biosynthesis pathway
Based on genome annotations and protein identification, the fatty acid biosynthesis pathways of Schizochytrium, including the saturated fatty acid synthesis pathway and the polyunsaturated fatty acid synthesis pathway were reconstructed (Fig. 3). Annotations of these proteins and their detailed LC-MS-MS information are provided in Additional file 3: Table S3 and Additional file 2: Table S2, respectively.
As shown in Fig. 3, the FAS pathway and the PKS pathways are, respectively, responsible for saturated fatty acid and polyunsaturated fatty acid synthesis. There are two principal pathways involved in short saturated fatty acid (SSFA) synthesis in Schizochytrium. Similar to fungi that use type I fatty acid synthase (FAS) to synthesize short saturated fatty acids in the cytosol, Schizochytrium utilizes a single large multifunctional polypeptide to synthesize short saturated fatty acids, which contains 4230 amino acids and shares approximately 99% protein identity with the fatty acid synthase from Aurantiochytrium mangrovei (AKV56231.1).
In addition to type I fatty acid synthase, the type II short saturated fatty acid biosynthesis pathway is also believed to be able to synthesize short saturated fatty acids in the cytosol. Type II short saturated fatty acid biosynthesis pathways are usually found in archaea and bacteria, which contain a set of discrete monofunctional enzymes . Most of the genes related to the type II short saturated fatty acid biosynthesis pathway were annotated and identified in Schizochytrium (Fig. 3, Additional file 2: Table S2 and Additional file 3: Table S3), except for 3-hydroxyacyl-ACP dehydratase (FabZ). 3-Hydroxyacyl-ACP dehydratase catalyzes the third step of short saturated fatty acid biosynthesis, dehydrating 3-hydroxyacyl-ACP to form trans-2-enoyl-ACP. These results suggest that the type II short saturated fatty acid pathway in Schizochytrium sp. TIO01 is incomplete and that type I FAS is responsible for the synthesis of C16:0 and odd-length short saturated fatty acids (C15:0 and C17:0) . In Schizochytrium, the long-chain fatty acid elongation process uses C16:0 as a substrate to synthesize long-chain saturated fatty acids (C24:0). This procedure consists of four sequential reactions: condensation, reduction, dehydration, and reduction. Genes participating in the elongation process were all annotated (Additional file 3: Table S3).
Similar to other Thraustochytriidae species that use the PKS pathway to synthesize polyunsaturated fatty acids (PUFAs) , Schizochytrium contains three large multifunctional PKS pathway genes, namely, pfaA, pfaB and pfaC. These three genes are highly homologous (degrees of protein identify > 99%) to proteins from another PUFA-producing species, Aurantiochytrium sp. L-BL10. In addition to the PKS pathway, the FAS pathway, which involves serial desaturation and elongation of short saturated fatty acids, is also believed to be able to synthesize PUFAs . Several elongases associated with the FAS pathway for PUFA synthesis were annotated and identified in Schizochytrium (Fig. 3, Additional file 2: Table S2 and Additional file 3: Table S3), however, the desaturase associated with the FAS pathway for PUFA synthesis was completely absent. The intermediates of PUFA compounds involved in the FAS pathway for DHA synthesis, such as C18:1, C18:2, C18:3, and C20:3, were also not detected by GC. These results indicate that the FAS pathway for DHA synthesis is incomplete in Schizochytrium and that the PKS pathway is responsible for DHA synthesis.
Transcriptomic profiling under cold temperatures
To improve the understanding of the mechanism by which low temperatures promote DHA accumulation, transcriptomic analysis was performed on six samples at both normal and lower temperatures as described previously. The obtained reads represented an average of 207.11 times the Schizochytrium genome length. Among these reads, 11,215 expressed genes were detected, including 9660 genomics predicted genes and 1555 novel predicted genes. A total of 1546 genes among 11,215 expressed genes were significantly associated with the low-temperature response (|fold change|> = 2 and adjusted p value < = 0.001), including 1237 downregulated genes and 309 upregulated genes. Using KEGG pathway analysis, 846 of the 1546 differentially expressed genes were annotated in 45 pathways. These differentially expressed genes were involved in amino acid metabolism, carbohydrate metabolism, lipid metabolism, global and overview maps, signal transduction, the endocrine system, and so on. The top 20 KEGG pathways of differentially expressed genes are shown in Fig. 2b.
Significant differentially expressed genes involved in fatty acid synthesis
As shown in Fig. 1b, Schizochytrium exhibited reduced saturated fatty acid content and a preference for increased DHA content when cultured under cold temperatures. Transcriptomic and q-PCR analyses indicated that genes involved in the FAS pathway responsible for saturated fatty acid synthesis (such as those encoding fatty acid synthase (FAS), long-chain fatty acid elongase (ELO), and very-long-chain 3-oxoacyl-CoA reductase (VLCR)) were significantly downregulated (Fig. 4 and Table 2). Gene expression was also significantly downregulated for malic enzyme, which catalyzes the conversion of malate to pyruvate to produce NADPH, which enters exclusively into the FAS pathway for saturated fatty acid synthesis [33, 34]. These results were consistent with the reduced saturated fatty acid content (C15:0, C17:0 and C24:0) at low temperatures. Generally, downregulation of the FAS pathway and malic enzyme is responsible for reduced saturated fatty acid content at low temperatures, however, the production of C16:0 was not affected in this study (Fig. 1b). Further analysis showed that two genes involved in fatty acid metabolism in organelles were significantly regulated (Fig. 4). One of these genes was significantly downregulated and encodes acyl-CoA dehydrogenase (ACADM), which catalyzes the first step of fatty acid beta-oxidation degradation in peroxisomes. The other gene was significantly upregulated and encodes enoyl-CoA hydratase (ECHS), which catalyzes the third step of fatty acid synthesis in mitochondria. These results suggest that fatty acid metabolism in organelles plays an important role in C16:0 synthesis and complements the biosynthesis of C16:0 under low-temperature conditions. Although low temperature has a significant impact on the production of DHA, genes involved in the PKS pathway responsible for the biosynthesis of DHA were not significantly regulated according to the transcriptomic analysis (Fig. 4) and q-PCR results (Table 2). This suggests that the increased DHA content was not a result of the overexpression of PKS pathway genes.
Significant differentially expressed genes associated with glycolysis and the pentose phosphate pathway
Glycolysis and the pentose phosphate pathway directly produce the substrates (acetyl-CoA and NADPH) for fatty acid synthesis. When Schizochytrium was cultured under cold conditions, genes involved in glycolysis such as those encoding glycolysis hexokinase (HK), 6-phosphofructokinase (PFKA), and phosphoglycerate mutase (PGAM) were significantly downregulated (|fold change|> = 2 and adjusted p value < 0.001) (Fig. 4 and Table 2), indicating that the glycolysis pathway was inhibited to some extent under cold stress. This result is consistent with the glucose consumption analysis (Fig. 1a).
Glucose-6-phosphate 1-dehydrogenase (PGD), which plays a critical role in microbial and algal lipid accumulation [31, 35], is involved in the pentose phosphate pathway and catalyzes the conversion of gluconate-6P to D-ribulose-5P to produce NADPH. Our findings show that this enzyme was significantly upregulated (Fig. 4 and Table 2). Song et al. and Liu et al. have reported that the NADPH produced by PGD exclusively enters the PKS pathway for PUFA biosynthesis [33, 34]. Taken together, these results indicate that glucose is degraded via the pentose phosphate pathway to produce a relatively large amount of NADPH entering the PKS pathway for DHA biosynthesis at low temperatures.
Significant differentially expressed genes involved in branched-chain amino acid metabolism and the TCA cycle
Amino acid degradation leads to the production of acetyl-CoA which is an important substrate for fatty acid synthesis. According to the results of the differentially expressed gene analysis, low temperatures induced significant regulation changes for 20 genes involved in branched-chain amino acid degradation (upregulation of 16 and downregulation of 4) (Fig. 5 and Table 2). Four significantly downregulated genes included genes encoding hydroxymethylglutaryl-CoA synthase (HMGCS), acyl-CoA dehydrogenase (ACADM) and ketol-acid reductoisomerase (BCATC). HMGCS catalyzes the reverse reaction of Leu degradation, leading to the consumption of acetyl-CoA. ACADM participates in Ile degradation and the fatty acid beta-oxidation pathway. There are six copies of ACADM genes in Schizochytrium, two of which were significantly downregulated. It is uncertain whether these two downregulated ACADM genes play roles in Ile degradation or in the fatty acid beta-oxidation pathway (Fig. 4). Generally, low temperatures enhanced the degradation of branched-chain amino acids leading to increased production of acetyl-CoA for fatty acid synthesis and TCA intermediate compounds (fumarate and succinyl-CoA).
Most of the genes that participate in the TCA cycle were not significantly regulated at lower temperatures except succinyl-CoA synthetase (LSC). Downregulated succinyl-CoA synthetase reduces the production of the TCA intermediate compound succinate. However, the levels of succinate and the downstream compound fumarate could be increased by enhanced branched-chain amino acid degradation (Tyr and Ile degradation pathways) (Figs. 4 and 5). This indicates that the TCA intermediate compounds (fumarate and succinyl-CoA) produced by branched-chain amino acid degradation could be further degraded to produce the acetyl-CoA precursor citrate through the TCA cycle. Taken together, more acetyl-CoA for fatty acid synthesis could be produced by branched-chain amino acid degradation and the TCA cycle at low temperatures.
Schizochytrium species are well known in the field of single-cell oil research for their ability to produce the health-benefitting docosahexaenoic acid. However, the PUFA biosynthesis pathway in Schizochytrium has long been controversial [10, 11, 36]. Based on the high-quality genome sequence of Schizochytrium sp TIO01, we found that the desaturase involved in the FAS pathway for PUFA synthesis was absent in Schizochytrium sp. TIO01, and the FAS pathway and the PKS pathway were, respectively, responsible for SFA and PUFA synthesis. This result was consistent with previous biochemical and genetic studies on Schizochytrium [7,8,9]. The PKS pathway responsible for PUFA synthesis also exists in other Thraustochytriidae species including one of the closest relatives of Schizochytrium, Thraustochytrium. In this species, desaturases involved in the FAS pathway for PUFA synthesis are completely absent. Evidence that the PKS pathway is responsible for PUFA synthesis includes a genomic survey , heterologous candidate gene expression  and an in vivo feeding experiment . In another close relative, Aurantiochytrium, the PKS pathway and the FAS pathway have been proposed to, respectively, produce PUFAs and SFAs based on fatty acid feeding experiments . Taken together, these results suggest that the PKS pathway, not the FAS pathway, is responsible for PUFA synthesis in Thraustochytriidae.
Low temperature has a significant impact on PUFA accumulation in many organisms that use the PKS pathway for PUFA synthesis. In this study, we found that low temperatures reduced saturated fatty acid biosynthesis and significantly increased DHA synthesis in Schizochytrium. Furthermore, gene expression analysis revealed that FAS pathway genes were significantly downregulated. In contrast, the expression of PKS pathway genes that are responsible for DHA synthesis remained unchanged relative to expression levels at normal temperature. These results suggest that significantly increased DHA content is not attributable to upregulation of PKS pathway genes. This phenomenon also occurs in Shewanella piezotlerans and Photobacterium profundum, where PUFA content is increased under low-temperature conditions and is not correlated with the overexpression of PKS pathway genes [20, 21]. In one of the closest relatives of Schizochytrium, Aurantiochytrium, Ma et al. reported that the upregulated PKS pathway genes and downregulated FAS pathway play an important role in promoting DHA accumulation [13, 18]. These results suggest that different mechanisms for regulating PUFA synthesis exist in different organisms.
NADPH and acetyl-CoA are generally accepted as the main substrates for fatty acid biosynthesis. Therefore, fatty acid production is promoted by the overexpression of genes involved in their production (malic enzyme, glucose-6-phosphate dehydrogenase, acetyl-CoA synthetase, citrate lyase, the pyruvate dehydrogenase complex) or by exogenous supplementation with an acetyl-CoA precursor (citrate, acetate) . The two major contributors to NADPH production include malic enzyme and glucose-6-phosphate 1-dehydrogenase (PGD) in the oxidative pentose phosphate pathway. When Schizochytrium was cultured at a low temperature, downregulated malic enzyme and upregulated glucose-6-phosphate 1-dehydrogenase (PGD) were detected. Although downregulation of malic enzyme reduces the production of NADPH, recent studies have shown that the NADPH produced by malic enzyme and glucose-6-phosphate 1-dehydrogenase enters the FAS pathway for SFA biosynthesis and the PKS pathway for PUFA biosynthesis, respectively [33, 34]. Downregulation of malic enzyme coupled with decreased saturated fatty acid content, and upregulation of glucose-6-phosphate 1-dehydrogenase coupled with increased polyunsaturated fatty acid content were also observed in this study. Generally, the upregulated pentose phosphate pathway is expected to increase the amount of NADPH entering the PKS pathway for DHA synthesis at low temperatures.
Glycolysis and branched-chain amino acid degradation pathways directly produce acetyl-CoA. When Schizochytrium was cultured under low-temperature conditions, glycolysis expression levels were downregulated while branched-chain amino acid degradation pathway levels were upregulated. Although downregulated glycolysis leads to reduced acetyl-CoA production, the upregulated branched-chain amino acid degradation pathway could increase the production of acetyl-CoA. In addition, the downregulated FAS pathway would decrease the consumption of acetyl-CoA, leading to decreased saturated fatty acid content. Generally, a relatively large amount of acetyl-CoA would be able to enter the PKS pathway for DHA synthesis at low-temperature conditions.
Taken together, upregulation of the pentose phosphate pathway and branched-chain amino acid degradation pathway (increasing the production of the substrates for DHA synthesis) and downregulation of the FAS pathway and malic enzyme (reducing the consumption of the substrates) likely play a key role in DHA accumulation in Schizochytrium under low temperatures.
In the present study, multi-omics approaches were used to investigate the effects of low temperatures on DHA biosynthesis and the underlying mechanism in Schizochytrium sp. TIO01. Our findings provide evidence that the PKS and FAS pathways are, respectively, responsible for DHA and saturated fatty acid synthesis in Schizochytrium. In addition, our work supports the findings that under low-temperature conditions, increased DHA content is not correlated with PKS pathway genes, but rather, increased substrate levels induced by low-temperature conditions entering the PKS pathway.
Growth, glucose consumption, and fatty acid analysis
The Schizochytrium sp. TIO01 strain used in this study was isolated previously and deposited in the China General Culture Collection Center (CGMCC no. 4603). Modified YPD medium (glucose 40 g/L, yeast extract 15 g/L, peptone 5 g/L, pH 7.0) with a salinity equivalent to 50% that of seawater was used to culture Schizochytrium.
Growth analysis: Schizochytrium was cultured overnight at 28 °C in YPD medium and then inoculated into fresh medium at 180 rpm and 28 or 16 °C. Samples were collected at certain intervals, and the optical density was analyzed at 650 nm. The Glucose (GO) Assay Kit (Sigma, GAGO20-1KT) was used to analyze the glucose concentration according to the manufacturer’s protocol. Samples used to analyze the fatty acids were cultured at 28 °C or 16 °C and collected when the OD650 reached approximately 0.7. The total lipid content was extracted from freeze-dried samples according to the Bligh–Dyer method  and analyzed by an Agilent 7890B GC instrument with H2 as the carrier gas using an Agilent J&W CP-Sil 88 column (100 m × 0.25 mm i.d., 0.2 μm film thickness; Agilent Technologies). The operating conditions were as follows: the injector and detector temperatures were kept at 250 °C and 260 °C, respectively; the temperature program was as follows: initial temperature 100 °C for 5 min, increased at 8 °C/min to 180 °C and holding at this temperature for 9 min, followed by increasing at 1 °C/min to 230 °C and holding at this temperature for 15 min. The carrier gas was H2 with a flow rate of 40 mL/min and the detector gas flow rate was 40 mL/min for H2, and 400 mL/min for air. The FAME mix (C4–C24) consisting of fatty acid standards was purchased from Sigma (cat. no. 18919-1AMP).
Genome assembly, prediction, and annotation
Genome sequencing sample preparation: Schizochytrium was cultured at 28 °C and 180 rpm and collected when the OD650 reached approximately 0.7. Total genomic DNA was extracted, and whole-genome shotgun sequencing was performed using the PacBio RS II or Illumina PE250 sequencing platform at the Beijing Genomic Institute (BGI, Shenzhen, China).
Transcriptomics sample preparation for gene prediction: 14 samples were cultured in 7 different media at 28 or 16 °C (the medium components are shown in Additional file 4: Table S4). Harvesting was performed, respectively, when the OD650 reached 0.7. Total RNA was extracted and sequencing on each sample was performed on the BGISEQ-500 (PE100) or Illumina platform (PE150) at the BGI, Shenzhen, China.
Genome assembly: Falcon  (fc_env_180425) was used to assemble the PacBio long reads multiple times to produce the longest contig N50. A two-step polishing strategy was used to improve accuracy. Initial polishing was performed by Arrow using PacBio long reads, and further corrective polishing was performed by Pilon-v1.22  with highly accurate Illumina PE250 reads. SSPACE-LongRead-v1.1  was used to construct a scaffold using PacBio reads longer than 10,000 bp. Contigs assembled by Platanus-v1.2.4  using Illumina PE250 reads were used to fill the gaps by GMcloser-v1.6.2 . Illumina PE250 reads and transcriptomic reads were analyzed by Bowtie2-v3.4.1  and Histat2-v2.10 , respectively, to assess genome completeness. BUSCO-v3.0  (database: eukaryota_odb9) was used to evaluate the completeness of the genome.
Gene prediction and annotation: MAKER-v3.01.02  was used to predict genes by combining results obtained from ab initio prediction (AUGUSTUS-v3.3.1  and SNAP-v2006-07–28 ), protein homology search (against protein sequences from the NCBI database: Aurantiochytrium sp. FCC1311 (GCA_002897355.1), Blastocystis hominis (GCA_000151665.1), Fistulifera solaris (GCA_002217885.1), Achlya hypogyna (GCA_002081595.1), Phytophthora parasitic (GCA_000247585.2), and Thalassiosira pseudonana (GCA_000149405.2); FGENESH-predicted proteins using the Phaeodactylum tricornutum parameters were also used as homologous proteins (https://www.softberry.com/berry.phtml?topic=fgenesh&group=programs&subgroup=gfind)), and transcriptomics (transcripts from 14 diverse conditions were generated using Trinity-v2.6.6 de novo assembly  and used as EST evidence). Predicted genes were annotated by alignment to the nonredundant protein sequence (Nr), Kyoto Encyclopedia of Genes and Genomes (KEGG), SwissProt, TrEMBL and InterPro databases using the Basic Local Alignment Search Tool (BLAST) -v2.7.1 + , with an E-value cutoff of 1E-5. Gene Ontology (GO) terms were assigned to the genes using the BLAST2GO pipeline .
iTRAQ labeling, LC–MS/MS and data analysis
iTRAQ labeling and MS/MS: six samples were cultured in YPD medium at two different temperatures (28 and 16 °C). Two samples were collected when the OD650 reached 0.7 at 28 °C (28log); two samples were collected when the OD650 reached 0.7 at 16 °C (16log); and two samples were collected when the OD650 reached 1.4 at 16 °C (16sta). After cell sonication, 100 μg of the proteins was reduced with 10 mM DTT and alkylated with 15 mM iodoacetamide (IAM). Samples were digested at 37 °C with sequencing-grade modified trypsin solution. Desalted peptides were labeled with iTRAQ reagents (6-plex kit, Applied Biosystems, USA) according to the manufacturer's protocol for the 8-plex iTRAQ Kit. Labels 113 and 116 were used to label 16log; 114 and 117 were used to label 16sta; and 115 and 118 were used to label 28log. The iTRAQ-labeled peptide mixture was dissolved and separated into 20 fractions, and each fraction was loaded by an autosampler onto a 2-cm C18 trap column with buffer A (5% ACN, 0.1%, FA, with a linear gradient from 5 to 35%) and buffer B (98% ACN, 0.1% FA) at a flow rate of 300 nL/min, followed by a wash with up to 80% buffer B. Intact peptides were detected in a Thermo Q-Exactive Orbitrap instrument at a resolution of 70,000. The electrospray voltage applied was 1.6 kV. The AGC target was 3e-6 for full MS and 1e-5 for MS2. For MS scans, the m/z scan range was 350 to 2000 Da.
iTRAQ data analysis: MS/MS data were processed using MaxQuant (v. 22.214.171.124). Tandem mass spectra were searched against the Schizochytrium protein database, which contains genomics-predicted proteins and transcriptomics-predicted novel proteins. Enzyme specificity was set as full cleavage by trypsin, with two maximum missed cleavage sites permitted. The precursor and fragment ion mass tolerances were 5 ppm and 0.02 Da, respectively. The precursor charged states allowed ranged from 1 to 5. Carbamidomethylation (Cys), iTRAQ 8plex (K) and iTRAQ 8plex (N-term) were set as fixed modifications, whereas oxidation (Met) and acetylation (protein N-terminal) were set as variable modifications. For all experiments, only unique peptides were considered for protein quantification. The estimated false discovery rate (FDR) thresholds for peptides and proteins were specified at a maximum of 1%. The minimum peptide length was set at 6.
Transcriptomics sample preparation and q-PCR validation
Transcriptomics sample preparation for analyzing gene expression under two different culture conditions: six samples were cultured at 28 and 16 °C in YPD medium. Harvesting was performed when the OD650 reached approximately 0.7. Total RNA was extracted and sequencing on each sample was performed on the Illumina platform (PE150) at BGI, Shenzhen, China. Data were analyzed using the Hisat2-StringTie pipeline .
q-PCR validation: RNA was extracted using the TRIzol method followed by DNase treatment. The Prime Script II 1st Strand cDNA Synthesis Kit (Takara Biotech Co., Ltd., Dalian, China) was used for cDNA synthesis. q-PCR was performed on a CFX Connect real-time system (Bio-Rad) using Bio-Rad SYBR Green, and 18S rDNA was used as an internal standard; the 2−ΔΔCT method was used to analyze the results . All primers were designed by Primer3Web (https://bioinfo.ut.ee/primer3/), and sequence information is provided in Additional file 5: Table S5.
Availability of data and materials
The raw sequencing reads used for genome assembly and annotation have been deposited in the SRA (Sequence Read Archive) database under Bioproject ID PRJNA483186. The Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession SMSO00000000. Transcriptomic data for analysis of differential gene expression under the two conditions have been deposited in the SRA database under Bioproject ID PRJNA517294. Protein LC-MS-MS data have been deposited in the IPRoX database (https://www.iprox.org/) under the identifier IPX0001511000.
Total fatty acid
Saturated fatty acid
Short saturated fatty acid
Polyunsaturated fatty acid
Fatty acid synthase
3-Oxoacyl-ACP synthase II
Long-chain fatty acid elongase
Very long-chain 3-oxoacyl-CoA reductase
Very long-chain (3R)-3-hydroxyacyl-CoA dehydratase
Polyketide synthase subunit A
Polyketide synthase subunit B
Polyketide synthase subunit C
Beta-ketoacyl-acyl-carrier-protein synthase domain
Ketoacyl reductase domain
Branched-chain amino acid aminotransferase
2-Oxoisovalerate dehydrogenase E1 component
2-Oxoisovalerate dehydrogenase E2 component
3-Methylcrotonyl-CoA carboxylase beta subunit
Next generation sequencing
Isobaric tags for relative and absolute quantitation
Innis SM. Dietary omega 3 fatty acids and the developing brain. Brain Res. 2008;1237:35–43.
Ward OP, Singh A. Omega-3/6 fatty acids: alternative sources of production. Process Biochem. 2005;40(12):3627–52.
Gao M, Song X, Feng Y, Li W, Cui Q. Isolation and characterization of Aurantiochytrium species: high docosahexaenoic acid (DHA) production by the newly isolated microalga, Aurantiochytrium sp. SD116. J Oleo Sci. 2013;62(3):143–51.
Barclay WR. Fermentation process for producing long chain omega-3 fatty acids with euryhaline microorganisms. U.S. Patent 451567B1, Sep 17, 2002
Marchan LF, Chang KJL, Nichols PD, Mitchell WJ, Polglase JL, Gutierrez T. Taxonomy, ecology and biotechnological applications of thraustochytrids: a review. Biotechnol Adv. 2018;36(1):26–46.
Zhao X, Dauenpen M, Qu C, Qiu X. Genomic analysis of genes involved in the biosynthesis of very long chain polyunsaturated fatty acids in Thraustochytrium sp. 26185. Lipids. 2016;51(9):1065–75.
Metz JG, Roessler P, Facciotti D, Levering C, Dittrich F, Lassner M, Valentine R, Lardizabal K, Domergue F, Yamada A. Production of polyunsaturated fatty acids by polyketide synthases in both prokaryotes and eukaryotes. Science. 2001;293(5528):290–3.
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(7):621–30.
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(8):739–47.
Chen W, Zhou PP, Zhang M, Zhu YM, Wang XP, Luo XA, Bao ZD, Yu LJ. Transcriptome analysis reveals that up-regulation of the fatty acid synthase gene promotes the accumulation of docosahexaenoic acid in Schizochytrium sp. S056 when glycerol is used. Algal Res. 2016;15:83–92.
Hoang MH, Nguyen C, Pham HQ, Van Nguyen L, Hai TN, Ha CH, Dai Nhan L, Anh HTL, Quynh HTH, Ha NC. Transcriptome sequencing and comparative analysis of Schizochytrium mangrovei PQ6 at different cultivation times. Biotech Lett. 2016;38(10):1781–9.
Ji XJ, Mo KQ, Ren LJ, Li GL, Huang JZ, Huang H. Genome sequence of Schizochytrium sp. CCTCC M209059, an effective producer of docosahexaenoic acid-rich lipids. Genome Announc. 2015. https://doi.org/10.1128/genomeA.00819-15.
Ma Z, Tan Y, Cui G, Feng Y, Cui Q, Song X. Transcriptome and gene expression analysis of DHA producer Aurantiochytrium under low temperature conditions. Sci Rep. 2015;5:14446.
Zhou PP, Lu MB, Li W, Yu LJ. Microbial production of docosahexaenoic acid by a low temperature-adaptive strain Thraustochytriidae sp. Z105: screening and optimization. J Basic Microbiol. 2010;50(4):380–7.
Aussant J, Guihéneuf F, Stengel DB. Impact of temperature on fatty acid composition and nutritional value in eight species of microalgae. Appl Microbiol Biotechnol. 2018;102:5279–97.
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(5):269–81.
Ren LJ, Feng Y, Li J, Qu L, Huang H. Impact of phosphate concentration on docosahexaenoic acid production and related enzyme activities in fermentation of Schizochytrium sp. Bioprocess Biosyst Eng. 2013;36(9):1177–83.
Ma Z, Tian M, Tan Y, Cui G, Feng Y, Cui Q, Song X. Response mechanism of the docosahexaenoic acid producer Aurantiochytrium under cold stress. Algal Res. 2017;25:191–9.
Cao M, Wang D, Mao Y, Kong F, Bi G, Xing Q, Weng Z. Integrating transcriptomics and metabolomics to characterize the regulation of EPA biosynthesis in response to cold stress in seaweed Bangia fuscopurpurea. PLoS ONE. 2017;12(12):e0186986.
Allen EE, Bartlett DH. Structure and regulation of the omega-3 polyunsaturated fatty acid synthase genes from the deep-sea bacterium Photobacterium profundum strain SS9. Microbiology. 2002;148(6):1903–13.
Wang F, Xiao X, Ou HY, Gai Y, Wang F. Role and regulation of fatty acid biosynthesis in the response of Shewanella piezotolerans WP3 to different temperatures and pressures. J Bacteriol. 2009;191(8):2574–84.
Zeng Y, Ji XJ, Lian M, Ren LJ, Jin LJ, Ouyang PK, Huang H. Development of a temperature shift strategy for efficient docosahexaenoic acid production by a marine fungoid protist, Schizochytrium sp. HX-308. Appl Biochem Biotechnol. 2011;164(3):249–55.
Cheng RB, Lin XZ, Wang ZK, Yang SJ, Rong H, Ma Y. Establishment of a transgene expression system for the marine microalga Schizochytrium by 18S rDNA-targeted homologous recombination. World J Microbiol Biotechnol. 2011;27(3):737–41.
Wang Z, Lou S, Hu F, Wu P, Yang L, Li H, He L, Lin X. Complete mitochondrial genome of a DHA-rich protist Schizochytrium sp. TIO1101. Mitochondrial DNA Part B. 2016;1(1):126–7.
Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.
Ganuza E, Yang S, Amezquita M, Giraldo-Silva A, Andersen RA. Genomics, biology and phylogeny Aurantiochytrium acetophilum sp. nov. (Thraustrochytriaceae), including first evidence of sexual reproduction. Protist. 2019;170(2):209–32.
Iwasaka H, Koyanagi R, Satoh R, Nagano A, Watanabe K, Hisata K, Satoh N, Aki T. A possible trifunctional β-Carotene synthase gene identified in the draft genome of Aurantiochytrium sp. Strain KH105. Genes. 2018;9(4):200.
Liu B, Ertesvåg H, Aasen IM, Vadstein O, Brautaset T, Heggeset TMB. Draft genome sequence of the docosahexaenoic acid producing thraustochytrid Aurantiochytrium sp. T66. Genomics Data. 2016;8:115–6.
Seddiki K, Godart F, Cigliano RA, Sanseverino W, Barakat M, Ortet P, Rébeillé F, Maréchal E, Cagnac O, Amato A. Sequencing, de novo assembly, and annotation of the complete genome of a new Thraustochytrid species, strain CCAP_4062/3. Genome Announc. 2018;6(11):e01335-e11317.
Jenke-Kodama H, Sandmann A, Müller R, Dittmann E. Evolutionary implications of bacterial polyketide synthases. Mol Biol Evol. 2005;22(10):2027–39.
Zhao X, Qiu X. Analysis of the biosynthetic process of fatty acids in Thraustochytrium. Biochimie. 2018;144:108–14.
Meesapyodsuk D, Qiu X. Biosynthetic mechanism of very long chain polyunsaturated fatty acids in Thraustochytrium sp. 26185. J Lipid Res. 2016;57:1854–64.
Song X, Tan Y, Liu Y, Zhang J, Liu G, Feng Y, Cui Q. Different impacts of short-chain fatty acids on saturated and polyunsaturated fatty acid biosynthesis in Aurantiochytrium sp. SD116. J Agric Food Chem. 2013;61(41):9876–81.
Liu B, Liu J, Sun P, Ma X, Jiang Y, Chen F. Sesamol enhances cell growth and the biosynthesis and accumulation of docosahexaenoic acid in the microalga Crypthecodinium cohnii. J Agric Food Chem. 2015;63(23):5640–5.
Xue J, Balamurugan S, Li DW, Liu YH, Zeng H, Wang L, Yang WD, Liu JS, Li HY. Glucose-6-phosphate dehydrogenase as a target for highly efficient fatty acid biosynthesis in microalgae by enhancing NADPH supply. Metab Eng. 2017;41:212–21.
Zhao B, Li Y, Li C, Yang H, Wang W. Enhancement of Schizochytrium DHA synthesis by plasma mutagenesis aided with malonic acid and zeocin screening. Appl Microbiol Biotechnol. 2018;102(5):2351–61.
Sun XM, Ren LJ, Zhao QY, Ji XJ, Huang H. Enhancement of lipid accumulation in microalgae by metabolic engineering. Biochim Biophys Acta Mol Cell Biol Lipids. 2019;1864(4):552–66.
Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. CanJ Biochem Physiol. 1959;37(8):911–7.
Chin CS, Peluso P, Sedlazeck FJ, Nattestad M, Concepcion GT, Clum A, Dunn C, O’Malley R, Figueroa-Balderas R, Morales-Cruz A, et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016;13(12):1050–4.
Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, Cuomo CA, Zeng Q, Wortman J, Young SK, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE. 2014;9(11):e112963.
Boetzer M, Pirovano W. SSPACE-LongRead: scaffolding bacterial draft genomes using long read sequence information. BMC Bioinformatics. 2014;15:211.
Kajitani R, Toshimoto K, Noguchi H, Toyoda A, Ogura Y, Okuno M, Yabana M, Harada M, Nagayasu E, Maruyama H, et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 2014;24(8):1384–95.
Kosugi S, Hirakawa H, Tabata S. GMcloser: closing gaps in assemblies accurately with a likelihood-based selection of contig or long-read alignments. Bioinformatics. 2015;31(23):3733–41.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT StringTie and Ballgown. Nat Protoc. 2016;11(9):1650–67.
Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12(1):491.
Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008;24(5):637–44.
Korf I. Gene finding in novel genomes. BMC Bioinformatics. 2004;5(1):59.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Xiong Q, Feng J, Li ST, Zhang GY, Qiao ZX, Chen Z, Wu Y, Lin Y, Li T, Ge F, et al. Integrated transcriptomic and proteomic analysis of the global response of Synechococcus to high light stress. Mol Cell Proteomics. 2015;14(4):1038–53.
The authors would like to thank Dr. Xuanmin Guang (BGI, Shenzhen, China) for his helpful advice on NGS data analysis.
This work was financially supported by the National Natural Science Foundation of China (no. 31700028); Scientific Research Foundation of Third Institute of Oceanography, SOA (no. 2016008); Natural Science Foundation of Fujian Province (no. 2016J05100).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of data and read alignments used in genome assembly and assessment.
Detailed information regarding the iTRAQ-identified proteins.
List of predicted genes involved in fatty acid synthesis.
Constitution of media used in RNA_seq sample preparation for protein-coding gene prediction.
List of q-PCR primer sequences used in this study.
About this article
Cite this article
Hu, F., Clevenger, A.L., Zheng, P. et al. Low-temperature effects on docosahexaenoic acid biosynthesis in Schizochytrium sp. TIO01 and its proposed underlying mechanism. Biotechnol Biofuels 13, 172 (2020). https://doi.org/10.1186/s13068-020-01811-y
- Low temperature
- Polyunsaturated fatty acid
- De novo genome assembly
- Differential expressed genes