Identification and characterization of an efficient acyl-CoA: diacylglycerol acyltransferase 1 (DGAT1) gene from the microalga Chlorella ellipsoidea

Background Oil in the form of triacylglycerols (TAGs) is quantitatively the most important storage form of energy for eukaryotic cells. Diacylglycerol acyltransferase (DGAT) is considered the rate-limiting enzyme for TAG accumulation. Chlorella, a unicellular eukaryotic green alga, has attracted much attention as a potential feedstock for renewable energy production. However, the function of DGAT1 in Chlorella has not been reported. Results A full-length cDNA encoding a putative diacylglycerol acyltransferase 1 (DGAT1, EC 2.3.1.20) was obtained from Chlorella ellipsoidea. The 2,142 bp open reading frame of this cDNA, designated CeDGAT1, encodes a protein of 713 amino acids showing no more than 40% identity with DGAT1s of higher plants. Transcript analysis showed that the expression level of CeDGAT1 markedly increased under nitrogen starvation, which led to significant triacylglycerol (TAG) accumulation. CeDGAT1 activity was confirmed in the yeast quadruple mutant strain H1246 by restoring its ability to produce TAG. Upon expression of CeDGAT1, the total fatty acid content in wild-type yeast (INVSc1) increased by 142%, significantly higher than that transformed with DGAT1s from higher plants, including even the oil crop soybean. The over-expression of CeDGAT1 under the NOS promoter in wild-type Arabidopsis thaliana and Brassica napus var. Westar significantly increased the oil content by 8–37% and 12–18% and the average 1,000-seed weight by 9–15% and 6–29%, respectively, but did not alter the fatty acid composition of the seed oil. The net increase in the 1,000-seed total lipid content was up to 25–50% in both transgenic Arabidopsis and B. napus. Conclusions We identified a gene encoding DGAT1 in C. ellipsoidea and confirmed that it plays an important role in TAG accumulation. This is the first functional analysis of DGAT1 in Chlorella. This information is important for understanding lipid synthesis and accumulation in Chlorella and for genetic engineering to enhance oil production in microalgae and oil plants. Electronic supplementary material The online version of this article (doi:10.1186/s12870-017-0995-5) contains supplementary material, which is available to authorized users.

In eukaryotes, three types of DGATs have been reported: the endoplasmic reticulum (ER)-localized DGAT1, DGAT2 and the soluble cytosolic DGAT3. Among them, DGAT1 and DGAT2 are responsible for the bulk of TAG synthesis in most organisms [9]. It has been proposed that these two enzymes have no redundant functions in TAG biosynthesis [10]. DGAT1 plays a dominating role in the determination of oil accumulation and fatty acid composition in seed oils [6], and DGAT2 may influence the content and composition of some plant seed oils containing unusual fatty acids (e.g., epoxy and hydroxyl) [11][12][13]. The role of the cytosolic DGAT3 has not yet been determined.
DGAT1s are ER membrane-bound proteins and possess six to nine transmembrane domains [14]. The most variable region of DGAT1 is the hydrophilic N terminus, which is quite unique for each DGAT1 and might serve distinct functions in different organisms [15]. Several conserved motifs, including acyl-CoA binding motif, DAG binding motif, the fatty acid-binding protein signature and a putative C-terminal ER retrieval motif, have been identified in DGAT1 [16]. Recently, site-directed mutagenesis was used to demonstrate the importance of some conserved residues in DGAT1s. For instance, mutagenesis at P216 and F439 in Tropaeolum majus DGAT1 resulted in a total loss of DGAT1 activity, while the substitution of S197 with alanine in a putative SnRK1 target site resulted in a strong increase in DGAT1 activity in the range of 38% to 80% [6].
The first eukaryotic DGAT1 gene was cloned from mouse [17], followed by isolation from other organisms [10,16,[18][19][20][21][22][23][24][25][26]. Many studies have investigated DGAT1s because of their important roles in TAG synthesis and have tried to use them to alter the quality and quantity of storage lipids in higher plants. For instance, the AS11 mutant of Arabidopsis, having reduced DGAT activity, showed a 75% reduction in seed lipids, but the expression of Arabidopsis DGAT1 in the AS11 mutant restored the wild-type levels of TAG and very-long-chain fatty acid content [27]. Moreover, the over-expression of AtD-GAT1 can greatly enhance the TAG content of transformed tobacco [19,28]. Subsequently, Tropaeolum majus DGAT1 significantly contributes to seed oil biosynthesis in wild-type Arabidopsis and Brassica napus by over-expression [6]. The co-expression of an epoxygenase from Stokesia laevis, SlEPX, and VgDGAT1 or VgDGAT2 from Vernonia galamensis greatly increased the accumulation of vernolic acid in both petunia leaves and soybean somatic embryos [13]. The over-expression of DGAT1 from Jatropha curcas showed an enhanced total oil content in seeds but did not show any phenotypic differences [25].
Unlike oil crops, microalgae have higher biomass production rates and many are exceedingly rich in oil. Therefore, microalgae have been regarded as potential resources for producing biodiesel, especially neutral lipids (e.g., triacyglycerols; TAGs) [29][30][31]. Many microalgal strains have the ability to accumulate substantial amounts of lipids in the form of TAGs under stress conditions, such as nitrogen starvation [30]. So far, several DGATs have been cloned and functionally characterized from microalgae. For instance, DGAT1-like [16] and DGAT2B [32] from the diatom Phaeodactylum tricornutum have been functionally characterized in a TAGdeficient mutant in the yeast Saccharomyces cerevisiae. Furthermore, the over-expression of PtDGAT2 in P. tricornutum resulted in a 35% increase in the neutral lipid content, and the fatty acid composition showed a significant increase in the proportion of polyunsaturated fatty acids [33]. Two DGAT2s (OtDGAT2A and OtDGAT2B) have been identified and characterized from Ostreococcus tauri, and OtDGAT2B possesses broad substrate specificity [34]. TpDGAT2 from the marine diatom Thalassiosira pseudonana significantly affects the fatty acid profile of TAG [35]. In Chlamydomonas reinhardtii, homology searches identified five DGAT2, encoded by DGTT1-DGTT5 [36]. Among them, DGTT1 and DGTT3 are active in TAG synthesis following nitrogen deprivation [37]. The expression of CrDGTT2 in Arabidopsis increased the leaf TAG content, with some molecular species containing very-long-chain fatty acids [38]. A gene encoding DGAT1 was also identified in C. reinhardtii after the transcript-based correction of gene models [39]. Other putative DGAT genes have been annotated in the genomes of some microalgae, such as Chlorella variabilis, Coccomyxa sp. C-169, Volvox carteri f. nagariensis, Ostreococcus lucimarinus, Fragilariopsis cylindrus and so on [40]. However, to date, there are few reports on the function of DGAT from the unicellular eukaryotic green alga, Chlorella, which is a desirable resource for producing biodiesel. Therefore, research on DGAT from Chlorella will advance our understanding of the molecular mechanisms underlying lipid metabolism during oil accumulation and will also provide a new means to improve the oil quality and content of microalgae and oil crops.
In the present study, we isolated a DGAT1 gene (CeD-GAT1) from C. ellipsoidea, a unicellular eukaryotic green alga that can be easily cultured under either autotrophic or heterotrophic conditions, and characterized its function in yeast and higher plants (Arabidopsis and B. napus). Compared with DGAT1s of higher plants, such as Glycine max, Arabidopsis and Brassica oleracea, CeDGAT1 could more effectively enhance fatty acid accumulation in the wild-type yeast (INVSc1). The over-expression of CeD-GAT1 can significantly enhance the seed oil content and seed weight in Arabidopsis and B. napus. Furthermore, the expression pattern of the isolated DGAT1 gene was investigated. This study would be helpful for understanding the function of DGAT from microalgae and for improving oil production in B. napus.
To examine the relationships among different sources of DGAT, a phylogenetic tree was generated from an alignment of the deduced amino acid sequences of CeD-GAT1 with 44 DGAT homologues from other species (including membrane-bound DGAT1 and DGAT2 and cytosolic DGAT3). The dendrogram grouped CeDGAT1 with members of the DGAT1 family. DGAT2 and DGAT3 formed a separate cluster different from DGAT1 branches (Fig. 1). In detail, CeDGAT1 clustered together with APDGAT1 (Auxenochlorella protothecoides) and CvDGAT1 (Chlorella vulgaris) and formed a distinct lineage. Furthermore, CeDGAT1 is more closely related to DGAT1s from animals and fungi than to DGAT1s from plants.
Protein analysis with the TMpred program [41] predicted nine strongly hydrophobic transmembrane regions, which is consistent with the nine transmembrane domains predicted for B. napus DGAT1, T. majus DGAT1 and Arabidopsis DGAT1 [6,9,20]. Using the PROSITE database [42], a number of putative functional motifs, including N-glycosylation, cAMP-and cGMPdependent protein kinase phosphorylation, protein kinase C phosphorylation, casein kinase II phosphorylation, tyrosine kinase phosphorylation, N-myristoylation and amidation sites were identified in CeDGAT1 (Additional file 1: Table S1). Compared with AtDGAT1, N-glycosylation and amidation sites were found only in CeDGAT1, while a leucine zipper pattern was detected only in AtDGAT1. It remains to be determined whether these sites are important for the functional regulation of the enzyme in vivo.
As previously reported [43,44], a consensus amino acid sequence for an acyl-CoA binding motif (Fig. 2, I) and a conserved sequence for the DAG-binding motif (Fig. 2, III) of DGAT were also found in CeDGAT1. There was also a fatty acid-binding protein signature spanning residues Ala571 to Asn587 (Fig. 2, II) containing a putative tyrosine phosphorylation site: Tyr582 [6]. The CeD-GAT1 protein contained an invariant proline (Pro381), which is thought to participate in presenting the fatty acyl group to the active site for esterification to (diacyl) glycerol and is critical for DGAT1 activity [45]. A highly conserved serine residue (Ser410) was essential for the activity of acyl-CoA: cholesterol acyltransferase, an enzyme closely related to DGAT1 [24,46]. CeDGAT1 showed a leucine zipper motif with only one conserved leucine (Leu386) in the sequence. A visual examination of CeDGAT1 also revealed the sequence of a putative C-terminal ER retrieval motif (YYHDW, Fig. 2, IV), which is similar to other DGAT1 proteins in plants [11].
Lipid analysis and expression pattern of CeDGAT1 in C. ellipsoidea To understand the relationship between CeDGAT1 expression and lipid synthesis, we investigated the expression pattern of CeDGAT1, the variation in biomass, and the total fatty acid (TFA) and TAG contents in heterotrophic cultures of C. ellipsoidea on nitrogen-replete and nitrogen-depleted (1/4 N) media (Fig. 3). As shown in Fig. 3a, the algal cell biomass increased more slowly in nitrogen-depleted cultures than in nitrogen-replete cultures, but both the TFA content and TAG content increased rapidly from 19% and 9% to 44% and 23%, respectively, from 36 to 108 h ( Fig. 3b and c). On the contrary, TFA and TAG contents were not significantly changed in the growth process of nitrogen sufficient condition. So there was a 108% and 212% increase in TFA content and TAG content, respectively at the 108 th hour under nitrogen-depleted culture condition than under nitrogen sufficient culture condition.
Quantitative real-time PCR was performed to examine the expression profiles of CeDGAT1 in C. ellipsoidea cells under nitrogen-replete and nitrogen-depleted (1/ 4 N) conditions (Fig. 3d). 18S rRNA was used as an internal reference control. We noted that CeDGAT1, which catalyses the last committed step in TAG biosynthesis, was downregulated under nitrogen-replete conditions. Nevertheless, CeDGAT1 showed transient upregulation, with its transcript level peaking at 84 h following the onset of nitrogen depletion and declining thereafter. The upregulation of CeDGAT1 was concomitant with the increase in the TFA and TAG contents under nitrogen deprivation, suggesting that CeDGAT1 was highly induced by nitrogen deprivation and that its increased expression coupled with lipid content change may play an important role in TAG accumulation.

CeDGAT1 can recover the TAG synthesis of the quadruple mutant yeast strain H1246
To verify the diacylglycerol acyltransferase activity of CeDGAT1, the CeDGAT1 gene was heterologously expressed in the TAG-deficient S. cerevisiae quadruple mutant strain H1246 [47], which lacks the four genes DGA1, LRO1, ARE1 and ARE2 encoding DGAT, PDAT (phosphatidylcholine: diacylglycerol acyltransferase), ASAT1 (acyl-CoA: sterol acyltransferase 1) and ASAT2 (acyl-CoA: sterol acyltransferase 2), respectively. These four genes are essential for the formation of neutral lipids. Lipid bodies can be formed by the expression of at least one of four genes. INVSc1 and H1246 cells harbouring an empty pYES2.0 vector were used as positive and negative controls, respectively.
Nile Red was used to stain lipid bodies in yeast cells. Lipid bodies were present in the wild-type yeast strain INVSc1 or the mutant strain transformed with the CeD-GAT1 gene but were undetectable in the quadruple mutant strain carrying the empty expression vector (pYES2.0) (Fig. 4a). The total lipids were extracted from yeast cells and then subjected to TLC (Thin-Layer Chromatography) analysis. Upon expression of CeDGAT1, a prominent band corresponding to TAG appeared on the TLC plate as expected, whereas no TAGs were identified in mutant yeast cells lacking the endogenous yeast DGAT and PDAT activities (Fig. 4b), which is consistent with the results from Nile Red staining (Fig. 4a). These results suggest that CeDGAT1 can successfully restore the ability of the quadruple mutant strain H1246 to form neutral lipids and confirm that CeDGAT1 encodes a functional protein capable of catalysing the last step of TAG biosynthesis.
Heterologous expression of CeDGAT1 can more significantly increase the total fatty acid content in the wild-type yeast than DGAT1s from some higher plants We separately transferred CeDGAT1 and another three DGAT1 genes from higher plants, including oil crop G. max, A. thaliana and B. oleracea, into the wild-type yeast (INVSc1). RT-PCR results showed that the DGAT1 genes had nearly identical expression patterns in transgenic yeast (Additional file 2: Figure S1). The fatty acid contents of yeast carrying DGAT1 genes were measured by GC. As shown in Fig. 5, the total fatty acid contents of yeast significantly increased due to the expression of different DGAT1s when compared to the yeast transformed with pYES2.0. In detail, the total fatty acid content in yeast carrying AtDGAT1, GmDGAT1, BoDGAT1 and CeDGAT1 was 234.7 μg/mg, 243.3 μg/mg, 258.5 μg/ mg and 290.0 μg/mg, respectively. Among the yeast expressing different DGAT1 genes, the total fatty acid content in the yeast transformed with CeDGAT1 increased most remarkably, by 142%. In contrast (compared with the yeast expressing CeDGAT1), the total fatty acid contents of the yeasts transformed with the three DGAT1 genes from higher plants significantly decreased by 19% (AtDGAT1), 16% (GmDGAT1) and 11% (BoDGAT1), respectively. These results suggest that CeDGAT1 may Over-expression of CeDGAT1 enhances the seed oil content and seed weight in higher plants To explore CeDGAT1 as a tool to manipulate acyl-CoA pools and to engineer TAGs in plants, CeDGAT1 was over-expressed in Arabidopsis and B. napus var. Westar under the control of the constitutive NOS promoter. In Arabidopsis, three independent homozygous lines were selected for advancement to the T4 generation and used for detailed analysis. RT-PCR results showed that the CeDGAT1 transcript was expressed in transgenic lines over-expressing CeDGAT1 (Additional file 3: Figure S2). Identical amino acid residues are highlighted in black. Conserved motifs or putative signatures are boxed, such as the acyl-CoA binding signature (I); the fatty acid protein signature (II), which contains a tyrosine phosphorylation site (★); the DAG-binding site (III); and the putative endoplasmic reticulum retrieval motif in the C-terminus (IV). The highly conserved proline and serine residues are marked by asterisks. The region containing a conserved leucine repeat (L) is also marked The transgenic lines did not show any visible morphological difference from untransformed control plants (data not shown). GC analysis revealed that the transformation of wild-type Arabidopsis with CeDGAT1 leads to a higher seed oil content (Fig. 6a). The average total fatty acid content in wild-type Arabidopsis seeds was 31.5 mg/100 mg of seed, but increased to 33.9-43.1 mg/ 100 mg in the transgenic lines. Thus, the transgenic seeds displayed an approximately 8-37% higher oil content than that of wild-type plants. In addition, the 1,000seed weight of transgenic plants was 9-15% greater than that of the control (Fig. 6b).
In B. napus, twelve transgenic plants were identified by PCR (Additional file 4: Figure S3), and four of the twelve independent T3 transgenic B. napus lines were chosen for further analysis. The RT-PCR results showed that the CeDGAT1 transcript was expressed in transgenic lines (Additional file 3: Figure S2). Again, transgenic B. napus lines did not show any visible morphological difference from wild-type plants (data not shown). The fatty acid content in the seeds of wild-type B. napus and transgenic lines was measured using GC. As shown in Fig. 7a, the average total fatty acid content in wild-type B. napus seeds was 39.6 mg/100 mg of seed. In lines expressing CeDGAT1 under the NOS promoter, the total fatty acid content increased to 44.5-46.8 mg/100 mg of seed, representing an increase of 12-18% over the control. Moreover, the average 1,000-seed weight in the CeDGAT1 transgenic lines increased by 6-29% compared to that of wild-type plants (Fig. 7b). These results indicated that the CeDGAT1 gene can stimulate fatty acid biosynthesis and enhance seed weight.

Subcellular localization of CeDGAT1
EGFP-tagged CeDGAT1 was expressed in tobacco BY-2 suspension cells and then examined under confocal laserscanning microscopy. As shown in Fig. 8, the distribution pattern of EGFP-CeDGAT1 was similar to that of endogenous ER stained with ER-tracker™ Red, and the signal was visualized around the nucleus and thread-like ER networks, indicating typical ER localization.

Discussion
TAGs are quantitatively the most important storage form of energy for eukaryotic cells. The synthesis of TAG from DAG by DGAT is believed to be the major flux control step in oil biosynthesis. Much research has focused on DGAT because it is an enzyme unique to TAG synthesis in plants. However, the function of DGAT1 from Chlorella has not been reported.
In this study we cloned and characterized a novel DGAT1 gene (CeDGAT1) from C. ellipsoidea. Proteinprotein BLAST showed that CeDGAT1 shared no more than 40% identity with DGAT1s of higher plants, which resulted in a difference in the predicted three-dimensional structures (Additional file 5: Figure S4). Functional characterization in yeast showed that CeDGAT1 can increase the TAG content more than can AtDGAT1, GmDGAT1 and BoDGAT1, resulting in a significant increase in the total lipid content of yeast of 142%. Further investigations of the relationships between the CeDGAT1 activity and structure are needed. Its higher activity provides a scientific and economic basis for the use of C. ellipsoidea as an oil-producing alga to produce more oil in a short time.
In higher plants, the expression of DGAT generally correlates with oil deposition in developing seeds [21]. For soybeans, a stronger expression of DGAT1 was found in developing seeds than in other tissues [13]. However, DGAT1 transcripts were also detected in other plant tissues, e.g., AtDGAT was expressed in a wide range of tissues but most strongly in developing embryos and flower petals [18]. DGAT1 is also highly expressed during pollen development, presumably contributing to TAG accumulation in the pollen grain [48]. These findings suggested that these DGAT enzymes may be related to other physiological processes in addition to seed oil synthesis [15]. For unicellular eukaryotic green algae, all physiological processes take place within a cell; thus, the expression of DGAT1 can directly reflect the dynamics of TAG accumulation.
The effect of nutrition pattern alteration on algal cell growth, lipid accumulation, and cellular component   Fig. 6 The detection of the fatty acid content and the 1,000-seed weight of CeDGAT1 transgenic Arabidopsis. a Average total fatty acid content in dry seeds of different homozygous (T4) transgenic Arabidopsis lines. b Average 1,000-seed weight (expressed as milligrams of weight/1,000 seeds) of transgenic Arabidopsis T4 seeds. Col-0, wild type Arabidopsis; CeDGAT1-1, CeDGAT1-9, and CeDGAT1-12, transgenic Arabidopsis lines expressing CeDGAT1; C16:0, palmitic acid; C18:0, stearic acid; C18:1, oleic acid; C18:2, linoleic acid; C18:3, linolenic acid; C20:1, eicosenoic acid; and total, the sum of the C16:0, C18:0, C18:1, C18:2, C18:3 and C20:1 contents. The bars are the standard deviations (SDs) of three technical repeats. For the same fatty acid component, numbers with different letters are statistically significant (P < 0.05) changes has been analysed in several studies [49][50][51][52][53][54]. Upon nitrogen starvation, both starches and lipids increased greatly within C. zofingiensis [51] and Nannochloropsis oceanica cells [54], and N-deficiency plus P-repletion was a promising lipid trigger to motivate lipid accumulation in C. protothecoides cells [50]. In C. reinhardtii, three genes encoding acyltransferases, DGAT1, DGTT1, and PDAT1, were induced by nitrogen starvation and are likely to play a role in TAG accumulation based on their patterns of expression [39]. At the transcript level in N. oceanica, enhanced TAG synthesis under N-depleted conditions primarily involved the upregulation of seven putative DGAT genes and the downregulation of six other DGAT genes [55]. However, the expression patterns of lipid biosynthesisrelated genes, including DGAT1 in Chlorella, have not been extensively studied in these processes. Our results revealed that the upregulation of CeDGAT1 was closely related to the significant increase in the TFA and TAG contents under nitrogen deprivation, suggesting that CeDGAT1 plays an important role in TAG accumulation. Our findings contribute to the understanding of the microalgal response to element deprivation and the mechanism of lipid synthesis and accumulation in Chlorella, but much remains to be elucidated regarding the precise contribution of N starvation to microalgal metabolism.
Several previous studies have reported that the genetic manipulation of DGAT can lead to increased oil biosynthesis in Arabidopsis and B. napus. For instance, the seed-specific over-expression of A. thaliana DGAT1 in wild-type Arabidopsis can increase the seed oil content by 11-28% and the seed weight by 2.5-32.3% [43]. Similarly, the seed-specific expression of TmDGAT1 in wildtype Arabidopsis resulted in a 10-33% net increase in the seed oil content and a 15% increase in the 1,000-seed weight in transgenic Arabidopsis. Furthermore, the seedspecific expression of TmDGAT1 in high-erucic acid B. napus led to a net increase of 11-15% in the seed oil content of transgenic plants [6]. In addition, the overexpression of AtDGAT1 and BnDGAT1 in canola under the control of the napin promoter led to an increase of 2.5-7% in the oil content [56]. The over-expression of JcDGAT1 in Arabidopsis under both CaMV35S promoter and a seed specific promoter resulted in a 30-41% increase in the seed oil content [25]. Our studies showed that the expression of CeDGAT1 in Arabidopsis and B. napus under the NOS promoter does indeed increase oil biosynthesis in transgenic seeds by approximately 8-37% and 12-18% over that of the control. In addition, neither the Arabidopsis nor the B. napus CeDGAT1 transformants showed significant changes in fatty acid composition. In some studies, however, there were alterations in the oil composition through DGAT expression. The over-expression of JcDGAT1 in Arabidopsis resulted in a significant decrease in oleic acid (C18:1) and an increase in linolenic acid (C18:3) [25], and the transgenic expression of Sesamum indicum DGAT1 in Arabidopsis led to an increase in eicosenoic acid (C20:1) and a reduction in oleic acid (C18:1) in seed oil [26]. More importantly, the expression of CeDGAT1 in Arabidopsis and B. napus under the NOS promoter also led to a significant increase in the average 1,000-seed weight in CeDGAT1 transgenic lines, by 9-15% and 6-29% in Arabidopsis and B. napus, respectively, and thus there was no decrease in the 1,000-seed weight caused by the oil content increase. Considering the constitutive expression of CeDGAT1 under the NOS promoter, larger increments in seed oil biosynthesis and seed weight can probably be expected when using a seed specific promoter. Furthermore, there was a difference between CeDGAT1 transgenic Arabidopsis and B. napus, with respect to effects on seed oil biosynthesis and seed weight. The seed oil content increased more in transgenic Arabidopsis than in B. napus, but the average seed weight increase was greater in B. napus. Interestingly, the increase in the oil content on a per-1,000-seed basis was similar between transgenic Arabidopsis and B. napus, at approximately 25-50%. To date, there has been no report that the over-expression of DGAT1 can significantly increase the seed weight in the oil plant B. napus, although this effect has been reported in Arabidopsis. In addition, transgenic plants showed no other phenotypic differences. Therefore, CeDGAT1 should have great potential for increasing the net oil production of the oil plant oilseed rape.

Conclusion
We cloned a novel DGAT1 gene (CeDGAT1) from C. ellipsoidea. CeDGAT1 is novel protein, sharing a low identity (≤40%) with DGAT1s from higher plants. The expression of CeDGAT1 is highly related to rapid lipid accumulation in C. ellipsoidea under nitrogen deprivation. In yeast, the expression of CeDGAT1 can significantly increase the lipid content and shows greater ability for improving the lipid synthesis than DGAT1s from some higher plants, including that from soybean. Moreover, the expression of CeDGAT1 in Arabidopsis and oilseed rape can lead to a net increase in the 1,000-seed lipid content of transgenic plants of 25-50%. These findings should be helpful for understanding the function of DGAT from microalgae and the mechanism of lipid synthesis and accumulation and may also provide technology for enhancing lipid production in microalgae and oil plants.

Strains and growth conditions
Chlorella ellipsoidea was initially cultured mixotrophically in 1 L flasks containing 500 mL of sterilized Endo medium [57] and incubated at 25°C under illumination (100 μmol photons/m 2 /s) for one week with shaking at 160 r/min. These pre-cultured cells were transferred to Erlenmeyer flasks (3 L), each containing 1 L of fresh medium to a final volume of 1.5 L, and incubated for another 4 days. These pre-cultured cells, after centrifugation and washing with sterilized water, were sampled as the starting point (0 h). Then, the collected cells were resuspended at a density of approximately 9.5 g/L in a 20 L BioFlo 415 fermentor (New Brunswick Scientific, USA) containing 14 L of modified Endo medium, in which urea was used to replace KNO 3 in the original Endo medium. The media containing urea at 0.4 g/L and 1.6 g/L were named N-depleted (1/4 N) and Nrepletion media, respectively. The culture conditions were maintained at 25°C, and a thermocirculator was used to maintain a constant temperature in the fermentor by circulating water through the jacket. The fermentor was aerated with filtered ambient air at a flow rate of 0.5 vvm, and the pH was maintained at 6.8 using 1 M KOH. The cultures were sampled 24, 36, 60, 84 and 108 h after they were initiated using N-depleted (1/4 N) or N-repletion medium. The cells were harvested by centrifugation (5,000 g at 4°C for 5 min). Aliquots for RNA analysis and gene cloning were frozen in liquid nitrogen and stored at −80°C if not immediately used, while those for lipid analysis were washed with water and freeze-dried.
Cloning of a cDNA encoding DGAT1 from C. ellipsoidea Total RNA was isolated from algae cells of the exponential growth phase of C. ellipsoidea using the EasySpin RNA Extraction Kit (Aidlab Biotech, Beijing, China), and cDNA was prepared from 5 μg of total RNA-template with the ReverTra Ace qPCR RT Kit (Toyobo, Osaka, Japan). The coding sequence of CeDGAT1 was amplified using the gene-specific primers P1 and P2 (Additional file 6: Table S2)  conditions for PCR were as follows: denatured at 95°C for 10 min, followed by 30 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 2 min; and a final extension step of 72°C for 10 min. The amplified cDNA was cloned into the pEASY-Blunt vector (TransGen Biotech, Beijing, China), and the corresponding clones were verified by PCR and DNA sequencing.

Yeast expression vector construction and transformation
The full-length CeDGAT1 ORF was amplified using primers P3 and P4 (Additional file 6: Table S2) and subcloned between the Hind III and EcoR I sites of the pYES2.0 yeast expression vector. The Saccharomyces cerevisiae strains (the wild-type strain INVSc1 and the mutant strain H1246) were transformed using the LiAc method [58]. We also separately transferred another three DGAT1 genes from higher plants, including the oil crop Glycine max (accession no. AY496439.1), Arabidopsis thaliana (accession no. NM_127503.2) and Brassica oleracea (unpublished data from our laboratory) into the wild-type yeast (INVSc1). The primers that were used for DGAT1 gene cloning are shown in Additional file 6: Table S2. Transformants were selected on synthetic complete medium lacking uracil (SC-ura). For heterologous expression studies, the yeast strains were transferred into liquid SC-ura medium containing 2% (w/v) glucose at 30°C overnight and then induced by adding 2% (w/v) galactose and 1% (w/v) tergitol NP-40 (Sigma, Taufkirchen, Germany) for an additional 72 h at 20°C. The expression of DGAT1s in transgenic yeast was verified at the transcript level by RT-PCR (for the RT-PCR primers see Additional file 6: Table S2).

Plant vector construction and transformation
The complete CeDGAT1 was cloned into the plant expression vector pCAMBIA2301 under the control of the nopaline synthase (NOS) promoter and nos terminator, yielding pCAMBIA2301-CeDGAT1 (Additional file 7: Figure S5). The final binary vector was verified and then transferred into Agrobacterium tumefaciens strain GV3101 by the freeze-thaw method [59]. Arabidopsis plants were transformed by vacuum infiltration [60]. Brassica napus var. Westar was transformed using hypocotyl explants and the modified method of DeBlock et al. [61]. T1 generation seeds were selected on kanamycin (50 mg/L), and then the selected transformed plants were transferred to soil. T3 transgenic B. napus lines and homozygous T4 transgenic Arabidopsis lines were used for seed and oil analyses. Genomic DNA was isolated from B. napus var. Westar leaf material. The stable integration of the NOS: CeDGAT1: nos cassette into the genome of transgenic B. napus was checked by PCR amplification using the specific primers P21 and P22 (Additional file 6: Table S2). GUS histochemical staining of the leaves from the transgenic lines was also conducted as described by Jefferson et al. [62]. In the meantime, the expression of CeDGAT1 in Arabidopsis and B. napus was detected by RT-PCR using the CeDGAT1specific primers P11 and P12 (Additional file 6: Table S2). The Arabidopsis housekeeping gene actin (primers P23 and P24) and the B. napus housekeeping gene GAPDH (primers P25 and P26) were used as internal controls.

Nile Red staining and microscopy
The Nile Red staining was used to visualize the intracellular lipid bodies as an indicators of TAG formation [63]. For yeast cell staining, a 500 μL suspension of yeast cells in the culture medium was stained with 5 μL of Nile Red (1 mg/mL in acetone stock), incubated in the dark for 5 min, and immediately used for microscopic analysis.

Lipid analysis by TLC and GC-MS
For the analysis of lipids from yeast and C. ellipsoidea, the cells were harvested by centrifugation, and the resulting cell pellets were ground to a fine powder under liquid nitrogen and subsequently treated with isopropanol at 80°C for 10-15 min to stop the lipolytic activity. Isopropanol was evaporated under nitrogen gas before lipid extraction. The total lipids were extracted according to a modified version of the Bligh and Dyer method [64], and TAG was separated from the total lipids by thin-layer chromatography (TLC) on Silica Gel 60 plates (Merck, Darmstadt, Germany). The solvents that were used were hexane/diethyl ether/glacial acetic acid (70:30:1, v/v). The lipids were visualized by spraying Primuline (Sigma, 10 mg/100 mL acetone: water (60:40 v/v)) and exposing the plate to UV. Triolein (Sigma) was used as the standard. TAGs were recovered from the TLC plates and then trans-esterified with 5% H 2 SO 4 in methanol at 85°C for 1 h. The fatty acid methyl esters (FAMEs) were extracted with hexane and analysed by GC-MS following the methods described in the following section.

Fatty acid analysis
Cellular fatty acids were extracted by incubating 10 mg of dried seeds of control and transformed plants or 50 mg of yeast powder and freeze-dried algae powder in 3 mL of 7.5% (w/v) KOH in methanol for saponification at 70°C for 4 h. After the pH was adjusted to 2.0 with HCl, the fatty acid were subjected to methylesterification with 2 mL of 14% (w/v) boron trifluoride in methanol at 70°C for 1.5 h. A phase separation was produced by adding 1 mL of 0.9% (w/v) NaCl and 4 mL of hexane. The upper phase was dried under a nitrogen gas flow and resuspended in 0.3 mL of acetic ether prior to GC analysis. An analysis of fatty acid methyl esters (FAME) was performed by GC-MS (gas chromatography-mass spectrometry, TurboMass, PerkinElmer, USA) equipped with a capillary column (BPX-70, 30 m × 0.25 mm × 0.25 μm). Hydrogen was used as the carrier gas at a flow rate of 1.0 mL/min. The injector and detector temperatures were held at 250°C. The column oven was temperatureprogrammed from 100 to 190°C at 15°C/min, where the temperature was held for 1 min increased to 220°C at 10°C/min, and then held for 4 min. The total FA content was quantified using heptadecanoic acid (C17:0) (Sigma) as an internal standard added to samples prior to extraction.

Dry weight determination
For dry weight determination, the algal cells were collected by filtering the culture through pre-weighed Whatman GF/C filter paper (1.2 μm pore size). Then, the filter paper was dried at 80°C in an oven until the weight was constant.
The C. ellipsoidea biomass concentration (w/v) was equivalent to a specific value of the cell dry weight (DW) that was determined by OD 540 Table S2).
To normalize the transcript levels in each sample, 18S rRNA was used as the internal standard (primers P29 and P30). The relative expression was computed following the formula (Cta-Ctb), where Cta and Ctb are the average Ct values of the reference and target genes, respectively.

Subcellular localization of CeDGAT1 in tobacco BY-2 Cells
To determine its subcellular location, CeDGAT1 without the stop codon was amplified from cDNA, cloned into the GATEWAY donor vector pGWC according to the method described by Chen et al. [67] and sequenced. The CeDGAT1 was then introduced in the destination vector pGWB5 with infusion with EGFP under the CaMV 35S promoter by LR reaction following the manufacturer's instructions (Invitrogen). The final construct CeDGAT1::EGFP was transferred into Agrobacterium tumefaciens strain GV3101. The Agrobacterium-mediated transformation of tobacco BY2 cells was performed according to An [68] and Genschik et al. [69]. Briefly, BY2 cells were incubated for 3 days at 27°C in the dark (without shaking) with Agrobacterium GV3101 containing CeDGAT1::EGFP. Subsequently, cells were plated on medium containing two antibiotics: timentin (500 mg/L) to kill off Agrobacterium and kanamycin (100 mg/L) to select for transformed cells. Transformed cells appeared after one month as a callus on plates and were then transferred to fresh plates once a month. A suspension culture was obtained by the addition of small transformed callus clumps to liquid culture medium containing kanamycin. The transformed BY-2 suspension cells were stained with ER-tracker™ Red following the manufacturer's instructions (Invitrogen) and then observed and photographed using a Leica TCS SP5 confocal laser scanning microscope (Leica Microsystems, Germany).

Statistical analysis
All of the experimental data were statistically compared using a one-way analysis of variance (ANOVA) with the software Statistical Product and Service Solutions (SPSS) v19.0, followed by a post-hoc test to determine the significant difference among the treatment means.

Additional files
Additional file 1: Table S1. Putative functional motifs in CeDGAT1 and AtDGAT1. (DOCX 17 kb) Additional file 2: Figure S1. RT-PCR detection of DGAT1 genes in transgenic yeast (INVSc1). The yeast actin was used as an internal control. 1, The yeast transformed with pYES2.0; 2-5, the yeast expressing AtDGAT1, GmDGAT1, BnDGAT1 and CeDGAT1, respectively. (DOCX 55 kb) Additional file 3: Figure S2. A schematic map of the pCAMBIA2301-NOS-CeDGAT1-nos plasmid. The pCAMBIA2301-NOS-CeDGAT1-nos vector contained an expression box of the CeDGAT1 gene from C. ellipsoidea under the control of the NOS promoter and nos terminator; an expression box of the GUS gene controlled by the CaMV35S promoter and nos