Open Access

Artificial miRNA inhibition of phosphoenolpyruvate carboxylase increases fatty acid production in a green microalga Chlamydomonas reinhardtii

Biotechnology for Biofuels201710:91

DOI: 10.1186/s13068-017-0779-z

Received: 26 January 2017

Accepted: 6 April 2017

Published: 13 April 2017



Nutrient limitation, such as nitrogen depletion, is the most widely used method for improving microalgae fatty acid production; however, these harsh conditions also inhibit algal growth significantly and even kill cells at all. To avoid these problems, we used artificial microRNA (amiRNA) technology as a useful tool to manipulate metabolic pathways to increase fatty acid contents effectively in the green microalga Chlamydomonas reinhardtii. We down-regulated the expression of phosphoenolpyruvate carboxylase (PEPC), which catalyzes the formation of oxaloacetate from phosphoenolpyruvate and regulates carbon flux.


amiRNAs against two CrPEPC genes were designed and transformed into Chlamydomonas cells and amiRNAs were induced by heat shock treatment. The transcription levels of amiRNAs increased 16–28 times, resulting in the remarkable decreases of the expression of CrPEPCs. In the end, inhibiting the expression of the CrPEPC genes dramatically increased the total fatty acid content in the transgenic algae by 28.7–48.6%, which mostly increased the content of C16–C22 fatty acids. Furthermore, the highest content was that of C18:3n3 with an average increase of 35.75%, while C20–C22 fatty acid content significantly increased by 85–160%.


Overall our results suggest that heat shock treatment induced the expression of amiRNAs, which can effectively down-regulate the expression of CrPEPCs in C. reinhardtii, resulting in an increase of fatty acid synthesis with the most significant increase occurring for C16 to C22 fatty acids.


Chlamydomonas reinhardtii Lipids Biodiesel amiRNA CrPEPC1 CrPEPC2


Biodiesel has become the focus of the world’s new energy research to replace petroleum diesel as it comes from a wide variety of sources, is biodegradable [1], has a high calorific value, and above all is environmentally friendly [2, 3]. Microalgae are the biodiesel source of choice with high growth rate, high oil production, and ease of cultivation [4, 5], which make them the most promising biodiesel material [6]. However, the lipid content in microalgae must be improved if biodiesel from microalgae is to be used for commercial scale production in the future [7, 8].

In recent years, increasing lipid content in microalgae through genetic engineering techniques has been tried. However, the expression of acetyl coenzyme A carboxylase (ACCase) and fatty acid synthase (FAS) from higher plants in microalgae did not achieve the aim of dramatically increasing the oil content [9, 10]. Meanwhile, previous studies showed that there was a negative correlation between the activities of phosphoenolpyruvate carboxylase (PEPC), a key enzyme of the amino acid metabolic pathway, and lipid accumulation because of sharing its ACCase with a common substrate, pyruvate [11, 12]. Therefore, it is safe to propose that the inhibition of the activity of PEPC may promote the flow of carbon to fatty acid synthesis, thereby enhancing the oil content of cells (Fig. 1).
Fig. 1

Inhibition of PEPC by miRNA promotes the synthesis of fatty acids

Among numerous gene-silencing systems, microRNAs (miRNAs) are receiving much attention because of their high efficiency and accuracy. The miRNAs guide post-transcriptional regulation by means of targeted RNA degradation and translational arrest [13, 14]. Hence, amiRNA (artificial microRNA)-mediated knockdown technology as a novel technology has opened up the possibility of specifically down-regulating a gene of interest, avoiding problems associated with silencing of off-targets by siRNAs [15, 16]. Chlamydomonas reinhardtii, a unicellular green alga from the phylum Chlorophyta, is an ideal model system for studying the role of miRNAs in regulating gene expression because it has well-established genetic engineering systems in nucleus, chloroplast, and mitochondria [17, 18]. For example, pyruvate formate lyase (PFL1) expression was reduced by 80–90% at the transcriptional level by a designed amiRNA in Chlamydomonas [19], indicating that the technique amiRNA broadens the application prospects in the field of gene regulation in microalgae.

It is well known that there were two subtypes of CrPEPC enzymes, named as CrPEPC1 and CrPEPC2, in C. reinhardtii [20, 21]. Here we designed amiRNAs, based on CrPEPC1/2 sequences, to reduce their expression and explored their potential for increasing fatty acid accumulation in C. reinhardtii. This is the first report on the use of novel amiRNAs to silence CrPEPCs and improve the accumulation of lipid in C. reinhardtii. The transcription levels of amiRNA-PEPC1 and amiRNA-PEPC2 were drastically increased when induced by heat shock (HS) and the transcription levels of CrPEPC1 and CrPEPC2 were reduced accordingly, leading to an increase in fatty acid content.


Screening and identification of transgenic mutants

The recombinant vectors pH-aRP1 and pH-aRP2 were introduced into the cell wall-deficient C. reinhardtii strain CC-849 by a “glass-bead” method, using zeomycin (10 μg/mL) as the selective antibiotic. The positive transformants were then verified by genomic DNA-PCR amplification and Sanger sequencing (see Additional file 1: Figure S1). Individual transformants were verified and selected for further study. Mutant strains with the names such as aRP1.1 and aRP1.2 were derived from individual amicroRNA-PEPC1 transformants, and aRP2.1, aRP2.2, and aRP2.3 were from amicroRNA-PEPC2 transformants. The growths of WT and transgenic algae were similar and the highest biomass was around 8 × 106 cells/mL for all strains, indicating that the transformed genes had no effect on the growth of transgenic algal cells (see Additional file 2: Figure S2).

Effect of heat shock (HS) on amiRNA transcription in transgenic algae

Since the amiRNA vector contains an Hsp70A heat shock-responsive promoter, amiRNA expression level after HS treatment was evaluated using quantitative RT-PCR, with U4 snoRNA as the internal reference. The 30-min HS increased the amicroRNA-PEPC1 (aRP1) transcription level to 17–19 times that of the control (Fig. 2a). When the algae were cooled down to normal conditions (25 °C), the aRP1 expression level declined, but still remained three times higher than the control. Similarly, heat shock also induced amicroRNA-PEPC2 (aRP2) expression, at 16–28 times higher than the control (Fig. 2b).
Fig. 2

Effects of heat shock on amicroRNA-PEPC1/2 expression levels. NT indicates a normal condition; HS indicates heat shock; aRP1.1 and aRP1.2 represent two different strains of amicroRNA-PEPC1 (a); aRP2.1, aRP2.2, and aRP2.3 represent three different strains of amicroRNA-PEPC2 (b) transformants; asterisk and double asterisk indicate significantly different amicroRNA-PEPC1/2 gene expression in transgenic algae when compared with NT

Effects of HS on CrPEPC genes’ expression in transgenic algae

After HS, the transcription level of CrPEPC1 was increased quickly in the WT, and after 2 h it reached the highest level when cooled down to the normal conditions (Fig. 3a). In contrast, in Chlamydomonas mutants with amicroRNA-PEPC1 (aRP1.1–aRP1.2), CrPEPC1 transcript levels were dramatically reduced to only 7–9.8% of the control after HS treatment and returned to normal conditions (Fig. 3a). These results indicated that aRP1 transformation in Chlamydomonas inhibits the expression of the target gene CrPEPC1.
Fig. 3

Heat shock-reduced transcription level of CrPEPC1/2. a aRP1.1 and aRP1.2 represent two different strains of amicroRNA-PEPC1 transformant; b aRP2.1, aRP2.2, and aRP2.3 represent three different strains of amicroRNA-PEPC2 transformant; NT indicates normal condition; HS indicates heat shock; asterisk and double asterisk indicate CrPEPC genes’ expression in transgenic algae compared to that in non-transgenic C. reinhardtii CC-849 on the gene expression level significantly

Interestingly, in the WT algae, CrPEPC2 expression under HS was strikingly different from CrPEPC1 expression, suggesting different HS-triggered regulation of CrPEPC1 and CrPEPC2 in Chlamydomonas. Similar to CrPEPC1 in transgenic aRP2 cells, CrPEPC2 gene expression levels under the same conditions were reduced to less than that of the control group. For instance, under normal conditions the CrPEPC2 expression level decreased by 80% in aRP2 strain compared to the WT. In addition, 30 min after the HS the CrPEPC2 expression level in transformant aRP2.3 was reduced by 90% compared to the control, indicating that the amiRNA aRP2 significantly reduced CrPEPC2 transcripts in Chlamydomonas. The HS treatment reduced CrPEPC2 in both WT and transgenic cells. After return to normal conditions, CrPEPC2 gene expression in transgenic algae aRP2.1, aRP2.2, and aRP2.3 was reduced to 30–70% that of the control group. Based on the data analysis, it was confirmed that CrPEPC2 gene expression was suppressed by overexpressed amiRNAs in the aRP2 strains (Fig. 3b).

Fatty acid composition increases in genetically modified algae after HS treatment

Fatty acids (FAs) were extracted from WT and all transgenic algal cells under normal conditions and after HS. In WT cells under routine culture condition, fatty acid composition analysis indicated that the major long-chain fatty acids C16:0, C16:4, and C18:3n3 are the most abundant in proportion to the total fatty acid (TFA), accounting for up to 80.97%. Under normal conditions, the TFA contents of genetically modified algae strains including aRP1 and aRP2 strains normally increase by about 3–23% (see Table 1 for more details). TFA were further increased after HS treatments. More HS treatments were carried out and more TFA were produced in algal cells. For example, the TFA content in WT was 85.5 ± 1.5, 93.8 ± 3.01, and 106.8 ± 12.2 mg/g (dry weight, DW), after zero, one, and three heat shocks, respectively (see Tables 1, 2, 3 for more details).
Table 1

Fatty acids and TFA contents (% cell dry weight) in transgenic strains with amicroRNA-PEPCs as analyzed by GC–MS

Fatty acids










0.046 ± 0.001

0.053 ± 0.001

0.059 ± 0.001

0.165 ± 0.029

0.046 ± 0.003



0.035 ± 0.031


0.234 ± 0.004

0.29 ± 0.012

0.293 ± 0.022

0.369 ± 0.021

0.233 ± 0.032

0.254 ± 0.015

0.246 ± 0.022

0.314 ± 0.055


15.46 ± 0.477

16.411 ± 0.685

18.063 ± 0.277a

18.602 ± 1.161a

15.054 ± 0.19

21.709 ± 1.005a

21.268 ± 1.344a

19.14 ± 1.667a


2.247 ± 0.033

3.162 ± 0.117b

3.231 ± 0.465

3.849 ± 0.377a

2.36 ± 0.212

2.571 ± 0.164

2.838 ± 0.117a

2.727 ± 0.32


2.12 ± 0.066

2.561 ± 0.085a

2.22 ± 0.194

2.539 ± 0.18a

2.08 ± 0.205

2.544 ± 0.165a

2.115 ± 0.247

2.232 ± 0.286


22.987 ± 0.097

25.535 ± 0.862a

25.67 ± 1.452a

25.852 ± 1.363a

24.25 ± 1.418

25.057 ± 1.338

27.57 ± 0.696b

28.82 ± 1.088b

C18:1 t

0.818 ± 0.041

1.233 ± 0.108b

1.769 ± 0.414

2.649 ± 0.633a

0.91 ± 0.089

1.283 ± 0.109b

1.235 ± 0.295

0.969 ± 0.165

C18:2 t

3.353 ± 0.343

4.253 ± 0.365a

4.219 ± 0.081a

4.242 ± 0.288a

3.448 ± 0.211

4.189 ± 0.667

3.575 ± 0.596

5.034 ± 0.783a


5.806 ± 0.178

7.345 ± 0.49b

5.696 ± 0.309

6.768 ± 0.358a

5.688 ± 0.614

8.376 ± 0.586b

6.059 ± 0.297

6.682 ± 1.092


30.777 ± 0.465

31.353 ± 1.037

34.298 ± 2.817

32.517 ± 1.715

31.633 ± 2.849

32.311 ± 1.103

34.343 ± 1.283a

36.74 ± 0.781


1.504 ± 0.033

2.054 ± 0.046b

1.555 ± 0.153

1.795 ± 0.095b

1.599 ± 0.204

2.133 ± 0.032b

1.946 ± 0.067a

1.977 ± 0.246b


0.012 ± 0.003

0.089 ± 0.029a

0.055 ± 0.074

0.086 ± 0.049

0.026 ± 0.026



0.065 ± 0.067


0.031 ± 0.001

0.06 ± 0.002b

0.042 ± 0.016

0.062 ± 0.004b

0.042 ± 0.017

0.005 ± 0.002b

0.007 ± 0.007a

0.047 ± 0.029



0.081 ± 0.057


0.052 ± 0.003b

0.026 ± 0.037

0.011 ± 0.003a

0.008 ± 0.002a

0.08 ± 0.094


0.085 ± 0.001

0.216 ± 0.014b

0.146 ± 0.067

0.242 ± 0.015b

0.141 ± 0.065

0.102 ± 0.004b

0.105 ± 0.012a

0.18 ± 0.067


0.016 ± 0.004

0.093 ± 0.018a

0.045 ± 0.061

0.089 ± 0.006b

0.033 ± 0.045

0.014 ± 0.001

0.014 ± 0.002

0.074 ± 0.051


85.497 ± 1.581

94.788 ± 3.312a

97.377 ± 3.06a

99.878 ± 5.496a

87.56 ± 3.208

100.561 ± 4.856b

101.33 ± 4.842b

105.193 ± 3.618b

The table shows the average of all the observations from 4 repeated tests, in the form of mean ± standard error of representation, wild-type C. reinhardtii CC-849 represents a control, TFA represents the total fatty acid content. All data are expressed in mg/g (dry weight)

The data represent the means ± SD of three replicate experiments and were analyzed by Student’s t-test (n = 3)

aP <0.05; bP <0.01. aRP1.1, aRP1.2, aRP1.3, aRP1.4 is transformed amicroRNA-PEPC1 into C. reinhardtii; aRP2.1, aRP2.2, aRP2.3 is transformed amicroRNA-PEPC2 into C. reinhardtii

Table 2

Various fatty acid compositions and contents of the transformants under conditions of heat shock for 30 min

Fatty acids










0.049 ± 0.001

0.083 ± 0.002

0.083 ± 0.003

0.079 ± 0.001

0.085 ± 0.002

0.083 ± 0.003

0.085 ± 0.004

0.086 ± 0.001


0.373 ± 0.028

0.467 ± 0.009

0.452 ± 0.004

0.44 ± 0.015

0.502 ± 0.007

0.416 ± 0.055

0.431 ± 0.009

0.435 ± 0.007


17.646 ± 0.179

21.968 ± 0.931b

23.709 ± 1.406b

21.003 ± 0.758b

23.445 ± 0.293b

20.94 ± 0.719b

22.76 ± 0.406b

21.256 ± 1.432b


2.666 ± 0.234

3.032 ± 0.145a

3.62 ± 0.072b

3.364 ± 0.038b

3.353 ± 0.044b

2.386 ± 0.293

3.092 ± 0.045a

2.789 ± 0.148


3.37 ± 0.272

3.748 ± 0.114a

3.563 ± 0.023

3.126 ± 0.067

3.487 ± 0.02

3.51 ± 0.177

3.338 ± 0.036

3.059 ± 0.141


20.756 ± 1.093

24.823 ± 0.314b

28.023 ± 1.839b

28.09 ± 0.848b

28.529 ± 0.146b

23.025 ± 0.781a

24.66 ± 1.453b

24.532 ± 0.254b

C18:1 t

2.393 ± 0.387

1.889 ± 0.246

2.625 ± 0.112

1.591 ± 0.06b

2.133 ± 0.183

1.958 ± 0.589

2.354 ± 0.05

2.011 ± 0.447

C18:2 t

9.163 ± 0.516

7.913 ± 0.555a

6.73 ± 0.123b

4.806 ± 0.261b

9.775 ± 0.354

7.605 ± 1.309

11.37 ± 0.673b

10.042 ± 1.465


6.651 ± 0.565

7.272 ± 0.3

7.352 ± 0.123

7.06 ± 0.116

8.598 ± 0.273b

8.503 ± 0.591b

7.663 ± 0.365a

7.244 ± 0.289


28.757 ± 1.697

33.663 ± 0.646b

36.882 ± 1.682b

35.235 ± 1.615b

39.886 ± 0.278b

31.335 ± 1.129a

35.924 ± 1.12b

35.765 ± 0.892b


1.475 ± 0.023

1.483 ± 0.023

1.493 ± 0.038

1.656 ± 0.063b

1.663 ± 0.034b

1.516 ± 0.117

1.376 ± 0.096

1.354 ± 0.055b


0.053 ± 0.009

0.093 ± 0.003b

0.096 ± 0.008b

0.095 ± 0.002b

0.093 ± 0.003b

0.098 ± 0.009b

0.097 ± 0.005b

0.095 ± 0.005b


0.071 ± 0.002

0.105 ± 0.001b

0.107 ± 0.004b

0.106 ± 0.003b

0.119 ± 0.003b

0.111 ± 0.003b

0.118 ± 0.003b

0.119 ± 0.006b


0.042 ± 0.003

0.081 ± 0.002b

0.084 ± 0.005b

0.079 ± 0.001b

0.083 ± 0.004b

0.098 ± 0.017b

0.085 ± 0.005b

0.082 ± 0.003b


0.215 ± 0.003

0.381 ± 0.013b

0.382 ± 0.009b

0.339 ± 0.064b

0.375 ± 0.016b

0.36 ± 0.012b

0.363 ± 0.006b

0.299 ± 0.063a


0.074 ± 0.005

0.145 ± 0.002b

0.158 ± 0.015b

0.142 ± 0.001b

0.142 ± 0.001b

0.143 ± 0.004b

0.142 ± 0.002b

0.144 ± 0.005b


93.754 ± 3.01

107.14 ± 3.306b

119.21 ± 6.55b

107.21 ± 3.7b

122.27 ± 0.86b

102.1 ± 1.15b

113.86 ± 3.41b

109.312 ± 4.99b

The data represent the means ± SD of three replicate experiments and were analyzed by Student’s t-test (n = 3)

aP <0.05; bP <0.01. aRP1.1, aRP1.2, aRP1.3, aRP1.4 is transformed amicroRNA-PEPC1 into C. reinhardtii; aRP2.1, aRP2.2, aRP2.3 is transformed amicroRNA-PEPC2 into C. reinhardtii

Table 3

Various fatty acid compositions and contents of the transformants under heat shock conditions for 30 min ×3

Fatty acids










0.046 ± 0.001

0.087 ± 0.001

0.074 ± 0.007

0.079 ± 0.002

0.08 ± 0.003

0.11 ± 0.031

0.084 ± 0.001

0.081 ± 0.001


0.444 ± 0.016

0.625 ± 0.04

0.531 ± 0.047

0.556 ± 0.015

0.526 ± 0.011

0.69 ± 0.006

0.538 ± 0.007

0.523 ± 0.004


20.208 ± 0.377

29.869 ± 0.812b

27.318 ± 2.346b

28.64 ± 0.535b

28.231 ± 0.61b

35.076 ± 1.781b

27.86 ± 0.999b

28.1 ± 0.558b


3.865 ± 0.068

4.823 ± 0.826

4.98 ± 0.501b

4.68 ± 0.178b

4.114 ± 0.104b

4.335 ± 0.476

4.258 ± 0.128b

4.679 ± 0.108b


3.257 ± 0.079

4 ± 0.177b

3.533 ± 0.156

4.014 ± 0.069b

3.84 ± 0.077b

5.094 ± 0.061b

3.27 ± 0.111

3.218 ± 0.054


26.971 ± 0.282

36.832 ± 0.68b

34.999 ± 1.856b

38.697 ± 1.17b

35.366 ± 0.708b

33.142 ± 3.959a

36.463 ± 0.69b

38.57 ± 0.913b

C18:1 t

1.696 ± 0.177

4.201 ± 1.199b

2.764 ± 0.497b

2.22 ± 0.088b

1.872 ± 0.187

2.965 ± 0.32b

1.551 ± 0.119

1.857 ± 0.264

C18:2 t

5.687 ± 0.048

8.655 ± 1.497b

7.507 ± 0.558b

7.548 ± 0.165b

6.248 ± 0.065b

13.539 ± 0.682b

6.988 ± 0.43b

7.789 ± 0.208b


8.042 ± 0.054

9.797 ± 0.304b

8.918 ± 0.353b

9.98 ± 0.229b

9.17 ± 0.472b

13.672 ± 1.253b

9.549 ± 0.318b

10.004 ± 0.31b


34.281 ± 0.647

48.215 ± 2.613b

44.218 ± 2.514b

48.18 ± 1.527b

45.387 ± 1.151b

46.848 ± 5.121a

45.256 ± 0.94b

47.64 ± 1.178b


1.793 ± 0.051

2.179 ± 0.348

1.751 ± 0.079

2.154 ± 0.05b

2.165 ± 0.089b

2.077 ± 0.276

2.099 ± 0.041b

2.161 ± 0.086b


0.048 ± 0.001

0.097 ± 0.003b

0.089 ± 0.01b

0.093 ± 0.001b

0.097 ± 0.005b

0.125 ± 0.05a

0.111 ± 0.034a

0.095 ± 0.005b


0.063 ± 0.001

0.124 ± 0.019b

0.104 ± 0.006b

0.11 ± 0.004b

0.107 ± 0.001b

0.145 ± 0.046b

0.105 ± 0.002b

0.105 ± 0.003b


0.041 ± 0.002

0.091 ± 0.009b

0.078 ± 0.009b

0.081 ± 0.002b

0.083 ± 0.002b

0.111 ± 0.046b

0.084 ± 0.002b

0.081 ± 0.003b


0.266 ± 0.012

0.455 ± 0.006b

0.4 ± 0.027b

0.415 ± 0.012b

0.415 ± 0.005b

0.532 ± 0.11b

0.389 ± 0.009b

0.386 ± 0.002b


0.073 ± 0.004

0.14 ± 0.001b

0.142 ± 0.023b

0.144 ± 0.003b

0.144 ± 0.003b

0.188 ± 0.085a

0.142 ± 0.002b

0.147 ± 0.004b


106.78 ± 12.24

150.19 ± 2.882b

137.406 ± 8.66b

140.212 ± 9.72b

137.845 ± 4.65b

158.65 ± 18.42b

138.743 ± 8.4b

145.44 ± 11.61b

For applying heat shock (HS), 400 mL cells (WT and transgenic algae) in mid-logarithmic phase were incubated at 40 °C and light intensity 20 μmol/m2/s for 30 min. As for triple heat shocks, a 4-h recovery period was utilized between each heat shock

The data represent the means ± SD of three replicate experiments and were analyzed by Student’s t-test (n = 3)

aP <0.05; bP <0.01. aRP1.1, aRP1.2, aRP1.3, aRP1.4 is transformed amicroRNA-PEPC1 into C. reinhardtii; aRP2.1, aRP2.2, aRP2.3 is transformed amicroRNA-PEPC2 into C. reinhardtii

After HS treatments for three times, the TFA of WT was increased by 24.9% in the control. The mutants had significantly higher TFA production than their respective control before heat shock with 58.3% (aRP1.1), 41% (aRP1.2), 47.6% (aRP1.3), 57.3% (aRP1.14), 57.8% (aRP2.1), 36.9% (aRP2.2), and 38.3% (aRP2.3) (Table 3). What is more, compared with WT cells that underwent HS treatments three times, transgenic algal cells had significantly higher TFA than the respective WT controls, with increases of 40.7% (aRP1.1), 28.7% (aRP1.2), 38.2% (aRP1.3), 29.1% (aRP1.14), 48.6% (aRP2.1), 29.9% (aRP2.2), and 36.2% (aRP2.3), respectively (Fig. 4). The results showed that the heat shock-induced overexpression of CrPEPCs amiRNAs significantly increases the TFA content of transgenic Chlamydomonas strains.
Fig. 4

Total fatty acid content in transformants: a aRP1 before heat shock; b aRP2 before heat shock; c aRP1 after ×1 heat shock; d aRP2 after ×1 heat shock; e aRP1 after ×2 heat shock; f aRP2 after ×2 heat shock; g aRP1 after ×3 heat shock; and h aRP2 after ×3 heat shock

The increase in TFA content of transgenic algae was contributed by increases of nearly all the different fatty acids (Table 3). Interestingly, C18:3n3, C16:0, and C16:4 were among the most significant constituents after heat shocks. The fatty acid C18:3n3 had the highest content which increased from 9.93 to 13.9 mg/g (DW) with an average increase of 35.75% in transgenic strains when comparing with 3 × HS WT; the C16:0 contents increased from 7.11 to the maximum 16.87 mg/g (DW); C16:4 content increased from 6.17 to 11.6 mg/g (DW), with a 43% increase in the transformant aRP2.3. Moreover, although the low contents of long-chain fatty acids such as C20:3, C22:2, and C22:6 were observed, they significantly increased from 85 to 160% according to GC–MS after 3 × HS (Table 3). These results suggest that HS treatment to induce the expression of amiRNAs can effectively down-regulate CrPEPC genes’ expression in C. reinhardtii, resulting in an increase of fatty acid synthesis with the most significant increase occurring for C16 to C22 fatty acids.


In this study, the WMD3 platform was employed to design amiRNAs against the Chlamydomonas reinhardtii CrPEPC1 and CrPEPC2 genes, and the amiRNAs were successfully induced by the heat shock protein (Hsp70A) promoter-driven expression system. Using 42 °C as the hHS, we observed significantly increased amiRNAs, aRP1 and aRP2 against CrPEPC1 and CrPEPC2, respectively, which were overexpressed at least 16 times after HS. Hsp70A promoter known as heat-inducible promoter was described a lot, showing that the transcription levels of foreign genes controlled by Hsp70A promoter were significantly increased after incubation at 42 °C for 30 min [22, 23]. The similar results were observed here that amiRNA expression reached the highest level after incubation at 42 °C for 30 min. We used multi-heat shocks to maximize amiRNAs’ expression. In this research, the expression level of amiRNAs was higher than the others [24, 25]. Till now, there was no report on the increase of TFA by amiRNA technology in microalgae. Our results showed that the expression of amiRNAs was increased about 28 times after HS, while that of CrPEPCs were decreased 90%, leading to a 48.6% increase in TFA. Here, we showed a highly effective miRNA expression system that was inducible and significantly increased TFA in C. reinhardtii.

It is well known that there are two PEPC isoenzymes in C. reinhardtii, CrPEPC1 and CrPEPC2 [21]. CrPEPC1 is proposed as a structurally homologous tetramer with maximum catalytic activity, while CrPEPC2 works as a multi-subunit complex whose expression pattern is more complex than CrPEPC1 [26, 27]. Activities of CrPEPC1 were decreased by 39–50% using RNAi technology [12, 28]. Here we observe the differential transcription regulation between CrPEPC1 and CrPEPC2 under HS. For example, after HS, a sharp rise of CrPEPC1 in the WT algal cells was observed. On the other hand, the gene expression of CrPEPC2 declined constantly after HS. Results from several labs proved that it was reliable to silence target gene using amiRNA in C. reinhardtii [19, 24, 25]. It was shown that the amiRNAs designed in our experiment showed high specificity and efficiency compared to CrPEPCs . After HS, the expression levels of aRP1 and aRP2 were increased 19–28 times, respectively, resulting in specific silencing of CrPEPCs. Although CrPEPC1 and CrPEPC2 had different response patterns to HS treatment, their effective down-regulation by amiRNAs under HS treatment compared with the WT was similar to significant reductions. For instance, aRP1 reduced the expression level of CrPEPC1 by ~93% compared to that in WT, and aRP2 down-regulated CrPEPC2 expression levels by up to 90% compared to WT. Furthermore, we did not detect any off-target phenomenon in all transgenic algae even after 1-year cultivation in the lab, indicating the stability of transgenic algae. Additionally, the inducible expression system for miRNA can reduce the cell damage caused by the expression of amiRNA. Therefore, our results showed that we succeeded in obtaining amiRNAs that were highly specific to silencing of CrPEPCs.

The observations demonstrate that we have constructed amiRNAs to inhibit the expression of CrPEPC genes in Chlamydomonas successfully. The designed amiRNAs own a high degree of complementarity with CrPEPC mRNAs to ensure that the RNA-induced silencing complex could bind to CrPEPC mRNAs resulting in the cleavage of them. It appears that decreasing CrPEPC mRNA abundance led to efficiently silencing their expression, redirecting carbon flux to fatty acid synthesis. As a result, the down-regulated CrPEPCs’ expression levels decreased the protein biosynthesis, and thus more carbon source flowed to the fatty acid synthesis pathway, resulting in the TFA increase. Although it was previously claimed that the activity of CrPEPC1 was inhibited by RNA interference, TAG levels only increased 20% as observed in a previous study [28], which never met the demand of industrial production. Moreover, using RNA interference technology requires a multi-step vector strategy with off-target effects so that it is likely not the best strategy for the manipulation of fatty acid synthesis [17]. Hence, the amiRNA with high specificity and efficiency was needed to improve TFA production in microalgae. In the present study, positive results were observed in the case of CrPEPC1 and CrPEPC2, the down-regulation of which resulted in significantly improved TFA content of microalgae, with a highest increase of 48.6% [158.65 mg/g (DW)] compared to that of WT.

We conducted GC–MS analyses to define the specifically changed FA components and found a significantly increase in long-chain fatty acid contents using our amiRNA method. Especially, the contents of C16 to C18 fatty acids generally increased by more than 30%. This is consistent with a previous report in oilseed rape seeds, with a substantial increase of C18 fatty acids following a sharp reduction of PEPC activity [29]. Moreover, a notable increase in the contents of long-chain fatty acids such as C20:3, C22:2, and C22:6 was also observed, which is useful for producing biodiesel.


The finding of this study indicated that overexpressing amiRNAs with an Hsp70A promoter remarkably inhibited the transcription of CrPEPCs, resulting in a significant increase of TFA in Chlamydomonas. This is the first report in which amiRNA was employed to increase TFA in Chlamydomonas. Although the expression levels of CrPEPC1 and CrPEPC2 are regulated in different ways, the amiRNAs remarkably inhibited them both at the transcription level, which made it possible to silence CrPEPC1 and CrPEPC2 at the same time. Moreover, TFA in transgenic algae was further increased with the help of multi-heat shocks and the contents of C16 and C22 fatty acids were significantly increased. In summary, this study demonstrates a high potential for the use of efficient CrPEPC amicroRNAs to improve the TFA in Chlamydomonas and provides a new possible way to increase the TFA production in other microalgae such as diatom for biodiesel.


Algal strains and culture conditions

A cell wall-deficient C. reinhardtii strain, CC-849, was obtained from the Chlamydomonas Genetic Center of Duke University (Duke University, Durham, NC, USA). Cells were cultured in TAP (Tris-acetate-phosphate) medium at 25 °C and under continuous cool-white fluorescent lamps (~20 μmol photons/m2 s) (normal condition, NT). For HS induction, 400 mL cells (WT and transgenic) in mid-logarithmic phase were incubated in a water bath at 42 °C for 30 min. After the treatment, algal cells were subjected to DNA/RNA extraction and GC–MS detection.

amiRNA constructs

The Web MicroRNA Designer (WMD3) website ( was used to design mature amicroRNA-PEPC1/2 to target the CrPEPC1/2 genes [17]. ami-PEPC1 (5´-TATTGGATTGAAAGGTCGCTA-3´) and ami-PEPC2 (5´-TTAACCAAACATTTTCGGCAC-3´) are complementary to nucleotides 117–137 bp of the coding region of CrPEPC1 and CrPEPC2, respectively. Then ami-PEPC1 and ami-PEPC2 replaced the mature sequence of C. reinhardtii endogenous miRNA-cre-MIR1162 (Accession Number: MI0006223) [24], to obtain amicroRNA-PEPC1 (aRP1) and amicroRNA-PEPC2 (aRP2), respectively (see Additional file 3). The constructed amicroRNA-PEPC1/2 (154 bp) was commercially synthesized in vitro and was cloned in plasmid pUC57 (Sangon Biotech Co., Ltd, Shanghai, China).

The expression vector pH124 containing the Ampicillin- and Zeomycin-resistant genes and a strong heat-inducible promoter (Hsp70A-RBCS2) was constructed in our lab [30]. Enzymatic digestion sites of NheI (GCTAGC) and PmaC I (CACGTG) were added at 5′- and 3′-end of aRP1/2, respectively. The DNA fragments of amicroRNA-PEPC1 (aRP1) and amicroRNA-PEPC2 (aRP2) were then ligated to the restriction enzyme sites, NheI and PmaCI, in the Chlamydomonas transformation vector pH124, resulting in pH-amicroRNA-PEPC1 (pH-aRP1) and pH-amicroRNA-PEPC2 (pH-aRP2), respectively (Fig. 5).
Fig. 5

Schematic diagram of the construct pH-amicroRNA-PEPC1/2. pH124 contained the ampicillin- and zeomycin-resistant genes and a strong heat-inducible promoter (Hsp70A-RBCS2). The DNA fragments of amicroRNA-PEPC1 (aRP1) and amicroRNA-PEPC2 (aRP2) were ligated to the restriction enzyme sites, NheI and PmaCI, within the Chlamydomonas transformation vector pH124, resulting in pH-amicroRNA-PEPC1 (pH-aRP1) and pH-amicroRNA-PEPC2 (pH-aRP2), respectively. The amiRNA precursors included amicroRNA-PEPC1/2, which were overexpressed under Hsp70A promoter

DNA-PCR analysis of transgenic alga

PCR analysis was utilized to confirm the successful transformation of our target amicroRNA-PEPC1/2. Genomic DNA was extracted from both WT and transgenic algae using the DNeasy kit (Takara, Japan). PCR was carried out according to the standard protocols to verify the presence of amicroRNA-PEPC1/2 on transgenic alga. PCR was performed using the primer pair pH124-F (5′TGACCTCCACTTTCAGCGACA3′) and pH124-R (5′ACTTGAGAGCAGTATCTTCCATCCA3′), which resulted in an amplicon of about 675 bp. PCR conditions were as follows: incubation at 95 °C for 2 min, followed by 25 cycles of 95 °C for 30 s, 60 °C for 15 s, and 72 °C for 15 s, plus a final extension for 7 min. All the amplified products were purified and verified by sequencing analyses (Sangon Biotech., Shanghai, China).

RNA extraction and RT-PCR

Total RNA was extracted using RNAiso Plus for Total RNA kit and DNA synthesis performed using PrimeScript™ Double-Strand cDNA Synthesis Kit (Takara Biotechnology Co., LTD, Dalian, China) according to the manufacturer’s instructions. The miRNA was extracted using RNAiso for Small RNA (Takara Biotechnology Co., LTD, Dalian, China), and reversed and polyadenylated using SYBR® PrimeScript™ miRNA RT-PCR Kit (Takara Biotechnology Co., LTD, Dalian, China) according to the manufacturer’s instructions. The amplified positive products were verified by sequencing analyses (Sangon Biotech., Shanghai, China).

Quantitative real-time PCR (qRT-PCR)

To quantitatively detect changes of amicroRNA-PEPC1/2 and CrPEPC1/2 expression in both wild-type and transgenic algae, qRT-PCR was performed with Applied Biosystems 7300 real-time PCR System (Framingham, MA, USA), and the primers used are listed in Additional file 4: Table S1. The standard protocol was applied to CrPEPC1/2 expression detection using SYBR Premix Ex Taq™II (Takara, Japan) according to the manufacturer’s instruction, while the expression of amicroRNA-PEPC1/2 was detected as described previously [31]. PCR conditions were as follows: one step of 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. Analysis of the melting curve of amplicons was used to test the specificity of the primers. The β-actin gene and U4 snoRNA were used as the reference genes in qRT-PCR detection of CrPEPCs and amicroRNA-PEPC1/2, respectively. Data with an R 2 above 0.998 were analyzed using the 2−ΔΔCt program [32]. Three technical replicates and two biological replicates were accustomed.

Nuclear transformation of C. reinhardtii

Genetic transformation of C. reinhardtii CC-849 was carried out using the “glass beads” method, according to Kindle [33] and Wang et al. [34].

Fatty acid methyl ester (FAME) transformation and FAME analysis

Total lipid extraction was performed as described by Lu et al. [35] with slight modifications. Briefly, 20 mg of lyophilized cells were suspended in 1 mL of 2 M NaOHCH3OH solution, shaken (100 rpm) for 1 h at 25 °C, and incubated at 75 °C for 15 min. After cooling down, the mixture was spiked with 1 mL of 4 M HCl–CH3OH and pH was set to below 2.0 with HCl, followed by incubation at 75 °C for 15 min. After that, FAMEs were extracted with 1 mL hexane, shaken by hand for 30 s, and then centrifuged at 4000g for 2 min. The hexane phase was collected and stored at −20 °C for further GC–MS analysis. 100 μg of C19 FA was added before extraction to estimate the recovery rate. Qualification and quantification of FAMEs were performed on a Thermo Trace GC Ultra Gas Chromatograph coupled to Thermo PolarisQ Mass Spectrometer which was equipped with an HP-5MS column (30 m × 0.25-mm id, film thickness 0.25 μm). The temperature of the injector was maintained at 250 °C. Helium was used as the carrier gas and ions were generated by a 70 eV electron beam and the mass range scanned was 50–650 m/z at a rate of 2 scan/s. The oven temperature for FAME analysis was initially maintained at 70 °C for 5 min followed by a temperature rate of 5 °C/min to 200 °C and then held for 5 min, 5 °C/min to 204 °C and then held for 2 min, 5 °C/min to 220 °C and then held for 3 min, and 5 °C/min to 255 °C and then held for 5 min. Peak identification was performed by matching the mass spectra of each compound with the National Institute of Standards and Technology mass spectral library. Automatic peak deconvolution was processed with MassLynx software (V4.1, Waters Corp., USA) [32]. The datasets of FAME profiling for further analysis were obtained by normalizing with the internal standards in the same chromatograms, respectively.

Statistical analysis

All exposure experiments were repeated three times independently, and data were reported as the mean with standard deviation (SD). For gene expression experiments, quantitative real-time PCR analysis was performed using the BioRAD iQ5 software. For each gene, the fold change expressed as mean ± SD (% control) was calculated using the (standard curve) approximation corrected for primer efficiency and normalized to housekeeping gene β-actin expression values. Statistical analyses were performed using the Student’s t test and Pearson Correlation analysis (SPSS13.0). For all of the data analysis, a p value <0.05 was considered statistically significant.





phosphoenolpyruvate carboxylase


pyruvate kinase


oxaloacetic acid


tricarboxylic acid cycle


fatty acid synthase


acetyl coenzyme A carboxylase


pyruvate decarboxylase


malonyl-CoA:ACP transacylase


artificial microRNA


acetyl coenzyme A carboxylase


total fatty acid


heat shock


normal condition


dry weight


Authors’ contributions

CGW, XC, and ZLH participated in the design of experiments, collected the data, and drafted the manuscript. CGW, XC, and JXW participated in data collection. CGW, HL, and JXW participated in the design of experiments and helped in writing the manuscript. All authors read and approved the final manuscript.


We would like to thank the anonymous reviewers for their constructive suggestions.

Competing interests

The authors declare that they have no competing interests.

Availability of supporting data

Data will be made available from the corresponding author on reasonable request.


This work was supported by the National Natural Science Foundation of China (Grant Nos. 31470389, 31470431), Guangdong Natural Science Foundation (2014A030308017, 2016A030313052), and Shenzhen Grant Plan for Science and Technology (JCYJ20160422171614147,JSGG20130411160539208, NYSW20140327010012).

Publisher's Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Guangdong Technology Research Center for Marine Algal Bioengineering, College of Life Sciences, Shenzhen University
Guangdong Key Laboratory of Plant Epigenetics, College of Life Sciences, Shenzhen University
Nanshan District Key Lab for Biopolymers and Safety Evaluation, Shenzhen University


  1. Demirbaş A. Diesel fuel from vegetable oil via transesterification and soap pyrolysis. Energy Sour. 2002;24(9):835–41.View ArticleGoogle Scholar
  2. Khan SA, Rashmi Hussain MZ, Prasad S, Banerjee UC. Prospects of biodiesel production from microalgae in India. Renew Sustain Energy Rev. 2009;13(9):2361–72.View ArticleGoogle Scholar
  3. Sharma T, Gour RS, Kant A, Chauhan RS. Lipid content in Scenedesmus species correlates with multiple genes of fatty acid and triacylglycerol biosynthetic pathways. Algal Res. 2015;12:341–9.View ArticleGoogle Scholar
  4. Minowa T, Yokoyama SY, Kishimoto M, Okakura T. Oil production from algal cells of Dunaliella tertiolecta by direct thermochemical liquefaction. Fuel. 1995;74(12):1735–8.View ArticleGoogle Scholar
  5. Janaun J, Ellis N. Perspectives on biodiesel as a sustainable fuel. Renew Sustain Energy Rev. 2010;14(4):1312–20.View ArticleGoogle Scholar
  6. Mata TM, Martins AA, Caetano NS. Microalgae for biodiesel production and other applications: a review. Renew Sustain Energy Rev. 2010;14(1):217–32.View ArticleGoogle Scholar
  7. Banerjee C, Singh PK, Shukla P. Microalgal bioengineering for sustainable energy development: recent transgenesis and metabolic engineering strategies. Biotechnol J. 2016;41(3):355–63.Google Scholar
  8. Radakovits R, Jinkerson RE, Darzins A, Posewitz MC. Genetic engineering of algae for enhanced biofuel production. Eukaryot Cell. 2010;9(4):486–501.View ArticleGoogle Scholar
  9. Dunahay TG, Jarvis EE, Roessler PG. Genetic transformation of the diatoms Cyclotella cryptica and Navicula saprophila. J Phycol. 1995;31(6):1004–12.View ArticleGoogle Scholar
  10. Dehesh K, Tai H, Edwards P, Byrne J, Jaworski JG. Overexpression of 3-ketoacyl-acyl-carrier protein synthase IIIs in plants reduces the rate of lipid synthesis. Plant Physiol. 2011;25(2):1103–14.Google Scholar
  11. Pathirana MS, Samac DA, Roeven R, Yoshioka H, Vance CP, Gantt JS. Analyses of phosphoenolpyruvate carboxylase gene structure and expression in alfalfa nodules. Plant Journal. 1997;12(2):293–304.View ArticleGoogle Scholar
  12. Deng X, Li Y, Fei X. The mRNA abundance of pepc2, gene is negatively correlated with oil content in Chlamydomonas reinhardtii. Biomass Bioenerg. 2011;35(5):1811–7.View ArticleGoogle Scholar
  13. Schwab R, Ossowski S, Riester M, Warthmann N, Weigel D. Highly specific gene silencing by artificial microRNAs in Arabidopsis. Plant Cell. 2006;18(51):1121–33.View ArticleGoogle Scholar
  14. Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000;403(6772):901–6.View ArticleGoogle Scholar
  15. Yu S, Pilot G. Testing the efficiency of plant artificial microRNAs by transient expression in Nicotiana benthamiana reveals additional action at the translational level. Front Plant Sci. 2014;5:622.Google Scholar
  16. Carbonell A, Fahlgren N, Mitchell S, Cox KLJ, Reilly KC, Mockler TC, Carrington JC. Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA, precursors. Plant J. 2015;82(6):1061–75.View ArticleGoogle Scholar
  17. Molnár A, Bassett A, Thuenemann E, Schwach F, Karkare S, Ossowski S, Weigel D, Baulcombe D. Highly specific gene silencing by artificial microRNAs in the unicellular alga Chlamydomonas reinhardtii. Plant J. 2009;58(1):165–74.View ArticleGoogle Scholar
  18. Molnár A, Schwach F, Studholme DJ, Thuenemann EC, Baulcombe DC. miRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii. Nature. 2007;447(7148):1126–9.View ArticleGoogle Scholar
  19. Burgess SJ, Tredwell G, Molnàr A. Artificial microRNA-mediated knockdown of pyruvate formate lyase (PFL1) provides evidence for an active 3-hydroxybutyrate production pathway in the green alga Chlamydomonas reinhardtii. J Biotechnol. 2012;162(1):57–66.View ArticleGoogle Scholar
  20. Rivoal J, Trzos S, Gage DA, Plaxton WC, Turpin DH. Two unrelated phosphoenolpyruvate carboxylase polypeptides physically interact in the high molecular mass isoforms of this enzyme in the unicellular green alga Selenastrum minutum. J Biol Chem. 2001;276(16):12588–97.View ArticleGoogle Scholar
  21. O’Leary B, Park J, Plaxton WC. The remarkable diversity of plant PEPC (phosphoenolpyruvate carboxylase): recent insights into the physiological functions and post-translational controls of non-photosynthetic PEPCs. Biochem J. 2011;436(1):15–34.View ArticleGoogle Scholar
  22. Schroda M, Blöcker D, Beck CF. The HSP70A promoter as a tool for the improved expression of transgenes in Chlamydomonas. Plant J. 2000;21(2):121–31.View ArticleGoogle Scholar
  23. Schroda M, Beck CF, Vallon O. Sequence elements within an HSP70 promoter counteract transcriptional transgene silencing in Chlamydomonas. Plant J. 2002;31(4):445–55.View ArticleGoogle Scholar
  24. Zhao T, Li G, Mi S, Li S, Hannon GJ, Wang XJ, Qi Y. A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev. 2007;21(10):1190–203.View ArticleGoogle Scholar
  25. Godman JE, Molnár A, Baulcombe DC, Balk J. RNA silencing of hydrogenase (-like) genes and investigation of their physiological roles in the green alga Chlamydomonas reinhardtii. Biochem J. 2010;431(3):345–51.View ArticleGoogle Scholar
  26. Mamedov TG, Moellering ER, Chollet R. Identification and expression analysis of two inorganic C- and N-responsive genes encoding novel and distinct molecular forms of eukaryotic phosphoenolpyruvate carboxylase in the green microalga Chlamydomonas reinhardtii. Plant J. 2005;42(6):832–43.View ArticleGoogle Scholar
  27. Moellering ER, Ouyang Y, Mamedov TG, Chollet R. The two divergent PEP-carboxylase catalytic subunits in the green microalga Chlamydomonas reinhardtii, respond reversibly to inorganic-N supply and co-exist in the high-molecular-mass, hetero-oligomeric Class-2 PEPC complex. FEBS Lett. 2007;581(25):4871–6.View ArticleGoogle Scholar
  28. Deng X, Cai J, Li Y, Fei X. Expression and knockdown of the pepc1, gene affect carbon flux in the biosynthesis of triacylglycerols by the green alga Chlamydomonas reinhardtii. Biotech Lett. 2014;36(11):449–568.View ArticleGoogle Scholar
  29. Sebei K, Ouerghi Z, Kallel H, Boukhchina S. Evolution of phosphoenolpyruvate carboxylase activity and lipid content during seed maturation of two spring rapeseed cultivars (Brassica napus L.). CR Biol. 2006;329(9):719–25.View ArticleGoogle Scholar
  30. Wu JX, Hu ZL, Wang CG, Li SF, Lei AP. Efficient expression of green fluorescent protein (gfp) mediated by a chimeric promoter in Chlamydomonas reinhardtii. Chin J Oceanol Limnol. 2008;26(3):242–7.View ArticleGoogle Scholar
  31. Shu L, Hu Z. Characterization and differential expression of microRNAs elicited by sulfur deprivation in Chlamydomonas reinhardtii. BMC Genom. 2012;13(1):1–12.View ArticleGoogle Scholar
  32. Lei AP, Chen H, Shen GM, Hu ZL, Chen L, Wang JX. Expression of fatty acid synthesis genes and fatty acid accumulation in Haematococcus pluvialis under different stressors. Biotechnol Biofuels. 2012;5(18):1–11.Google Scholar
  33. Kindle KL. High-frequency nuclear transformation of Chlamydomonas reinhardtii. Proc Natl Acad Sci USA. 1990;87(3):1228–32.View ArticleGoogle Scholar
  34. Wang CG, Peng XC, Wang JX, Lei AP, Li H, Hu ZL. A β-carotene ketolase gene (bkt1) promoter regulated by sodium acetate and light in a model green microalga Chlamydomonas reinhardtii. Algal Res. 2016;20:61–9.View ArticleGoogle Scholar
  35. Lu S, Wang J, Niu Y, Jie Y, Jian Z, Yuan Y. Metabolic profiling reveals growth related FEAM productivity and quality of Chlorella sorokiniana with different inoculum sizes. Biotechnol Bioeng. 2012;109(7):1651–62.View ArticleGoogle Scholar


© The Author(s) 2017