Frontiers reaches 6.4 on Journal Impact Factors

Data Report ARTICLE Provisionally accepted The full-text will be published soon. Notify me

Front. Bioeng. Biotechnol. | doi: 10.3389/fbioe.2018.00019

Phosphoproteome of the oleaginous green alga, Chlorella vulgaris UTEX 395, under nitrogen replete and deplete conditions.

  • 1National Renewable Energy Laboratory (DOE), United States

The unicellular green alga, Chlorella vulgaris UTEX 395, represents a promising biocatalyst for renewable biofuel production due to its relatively rapid growth rate and high lipid accumulation capacity (Guarnieri et al., 2011, Guarnieri et al., 2012, Gerken et al., 2013, Griffiths et al., 2014, Zuniga et al., 2016). Prior analyses have unveiled the global proteome dynamics of C. vulgaris following nitrogen depletion, which induces a high lipid accumulation phenotype (Guarnieri et al. 2011, Guarnieri et al., 2013). More recently, we have reported a draft genome, genome scale model, and nitrosoproteome for this alga (Zuniga et al. 2016, Guarnieri et al., 2017, Henard et al., 2017) providing further insight into lipid biosynthetic-, nutrient response-, and post-transcriptional- regulatory mechanisms. To further our understanding of these regulatory mechanisms and expand the knowledge base surrounding this organism, comparative phosphoproteomic analyses were conducted under nitrogen replete and deplete conditions to identify differentially phosphorylated proteins that will aid in the evaluation of the potential role of phospho-regulation in lipogenesis.

Algal cultivation: Chlorella vulgaris UTEX 395 was cultivated in biological triplicate, as described previously (Guarnieri et al. 2011, Guarnieri et al. 2013). Briefly, algae were inoculated at OD750 = 0.05 in 1L Roux bottles using modified Bold's Basal Media (mBBM). Cultures were maintained at 25°C±1°C, with continuous (24 hr) white fluorescent light illumination (200 µE m−2 s−1). Cultures were supplemented with 2% CO2/air and mixed with a magnetic stir bar at 500 rpm. 50mL of nitrogen replete cell culture was harvested at OD750 = 2, centrifuged for 5min at 5,000× g, and immediately quenched on liquid nitrogen. To induce nitrogen deprivation, the remaining cell culture was centrifuged for 5 min at 5,000× g, washed once in nitrogen-free mBBM, and resuspended in nitrogen-free mBBM for continued growth to OD750 = 4, followed by harvest, as described above. Cell pellets were immediately quenched in liquid nitrogen.

Protein Isolation: Cell pellets were thawed and solubilized on ice in 2 mL of lysis buffer (50 mM Tris, pH 8.0, 150 mM NaCl, 1 mM DTT, 10% glycerol, supplemented with 1× PhosSTOP Phosphatase Inhibitor Cocktail
(Roche Diagnostics Corporation). The cells were then sonicated on ice at 4°C, at 90% power setting for 30 seconds×6 cycles, with a one minute cool-down period between sonication cycles using a Braun-Sonic-L ultrasonicator. Lysates were cleared via two cycles of centrifugation at 16,000× g at 4°C for 30 minutes, and the supernatants were isolated for use in subsequent phosphoproteomic analysis.

Peptide Digestion and Phosphopeptide Enrichment: 500µg of lysate material was digested in solution with the following concentration of buffer and enzymes: 100mM DTT, 20mM HEPES, 100mM iodoacetamide, 25ng/µL trypsin (final E:S ratio of 1:20), and 20.5% TFA. 200µg of peptide were desalted using solid phase extraction on a Waters C18 Sep Pak 1cc cartridge (part no. WAT054960, as follows: column was activated with 1mL 100% acetonitrile followed by 1mL 50% acetonitrile, and equilibrated with 3 x 1mL 0.1% TFA. The sample was loaded on column and washed 2x with 1mL 0.1% TFA, followed by elution with 0.5mL 60% acetonitrile, 0.1% TFA. The eluate was dried to completion and phosphopeptides were enriched using the GL Sciences TiO2 kit (5010-21312), following manufacturer's instructions. Samples were dried to completion and resuspended in 120µL 0.1% TFA.

Phosphoproteomic analysis: Each enriched sample was analyzed by nano LC/MS/MS with a Waters NanoAcquity HPLC system interfaced to a ThermoFisher Q Exactive. The acquisition order was randomized. In each case peptides were loaded on a trapping column and eluted over a 75μm analytical column at 350nL/min; both columns were packed with Jupiter Proteo resin (Phenomenex), employing a 2h gradient, as follows: Buffer A - 0.1% formic acid in water; Buffer B - 0.1% formic acid in acetonitrile; t0 = 2%B, t1min = 5% B, t95 = 25% B, t110 = 35% B, t112 = 90% B, t113 = 2%B, t120 = 2%B. The mass spectrometer was operated in data-dependent mode, with MS performed in the Orbitrap at 70,000 FWHM resolution and MS/MS performed using HCD and product ions detected in the Orbitrap at 17,500 FWHM resolution. The fifteen most abundant ions were selected for MS/MS.

Data Processing: Data were processed through the MaxQuant software v1.5.1.0 ( for MS data recalibration, database filtering to 1% false discovery rate (FDR), calculation of peak areas (label-free), and assignment of phosphosite localization probability. Data were searched using a local copy of Andromeda with the following parameters: i) Enzyme: Trypsin/P ii) Database: Cv395_Maker_7100 (concatenated forward and reverse plus common contaminants), iii) Fixed modification: Carbamidomethyl (C), iv) Variable modifications: Oxidation (M), Acetyl (Protein N-term), Phospho (STY), v) Fragment Mass Tolerance: 20ppm. The Phospho(STY)Sites output was further processed using Excel. A total of 691 phosphosites (considering positional variances within the same peptide) were detected at 1% FDR.

Data Analysis: All analysis of intensities utilized the R language ( Proteins with no blastx identification were removed from the analysis. Intensities were filtered for presence in 2 of 3 replicates to remove erroneous measures and then compared for strict presence/absence in the two conditions. The filtered dataset of raw intensities was then analyzed for differential expression using the edgeR package (Robinson et al., 2010) and the top differentially expressed genes (p < 0.05) were analyzed further for gene ontology (GO). GO identifications for each protein were obtained in tab-separated format from the UniProt web interface ( and parsed in R to compare the top GO terms in the control replete cells versus the depleted cells. GO terms were split into molecular function, cellular component, and biological process categories when comparing between the treatments.

Phosphoproteome data deposition: A master file with the full list of phosphosites and associated localization and statistical data is deposited at the FigShare repository, along with differential expression and GO analyses. All files are accessible at:
Author Contributions
Project planning and wet lab execution were conducted by EK and MG. AG and CH directed phosphoproteome data processing and handling. MG and AG wrote the manuscript.

Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

This work was supported by National Renewable Energy Laboratory, Laboratory-Directed Research and Development (LDRD) grants 06511103 and 06511301 and the U.S. Department of Energy, Office of Energy Efficiency and Renewable Energy (EERE) under Agreement No. 22000. We kindly thank Drs. Richard Jones and David Allen for assistance in phosphoproteomic data acquisition.

Guarnieri, M. T., A. Nag, S. L. Smolinski, A. Darzins, M. Seibert and P. T. Pienkos (2011). Examination of triacylglycerol biosynthetic pathways via de novo transcriptomic and proteomic analyses in an unsequenced microalga. PLoS ONE 6: e25851. doi:10.1371/journal.pone.0025851

Guarnieri, M. T., L. M. L. Laruens, E. P. Knoshaug, Y. C. Chou, B. S. Donohoe and P. T. Pienkos (2012). Complex systems engineering: A case study for an unsequenced microalga, John Wiley & Sons, Inc. doi: 10.1002/9781118433034.ch8

Gerken, H. G., B. Donohoe and E. P. Knoshaug (2013). Enzymatic cell wall degradation of Chlorella vulgaris and other microalgae for biofuels production. Planta 237: 239-253. doi:10.1007/s00425-012-1765-0

Griffiths, M. J., R. P. van Hille and S. T. L. Harrison (2014). The effect of nitrogen limitation on lipid productivity and cell composition in Chlorella vulgaris. Appl Microbiol Biotechnol 98: 2345-2356. doi:10.1007/s00253-013-5442-4

Zuniga, C., C. T. Li, T. Huelsman, J. Levering, D. C. Zielinski, B. O. McConnell, C. P. Long, E. P. Knoshaug, M. T. Guarnieri, M. R. Antoniewicz, M. J. Betenbaugh and K. Zengler (2016). Genome-scale metabolic model for the green alga Chlorella vulgaris UTEX395 accurately predicts phenotypes under autotrophic, heterotrophic, and mixotrophic growth conditions. Plant Phys. doi:10.1104/pp.16.00593

Guarnieri, M. T., A. Nag, S. Yang and P. T. Pienkos (2013). Proteomic analysis of Chlorella vulgaris: Potential targets for enhanced lipid accumulation. J Proteomics 93: 245-253. doi: 10.1016/j.jprot.2013.05.025

Guarnieri, M. T., J. Levering, K. Zengler, C. A. Henard, J. L. Boore and E. P. Knoshaug (2017). Draft genome sequence of the oleaginous green alga, Chlorella vulgaris UTEX 395. Frontiers Plant Sci submitted.

Henard, C. A., M. T. Guarnieri and E. P. Knoshaug (2017). The Chlorella vulgaris S-Nitrosoproteome under nitrogen-replete and -deplete conditions. Front Bioeng Biotechnol Jan 17; 4:100. doi: 10.3389/fbioe.2016.00100

Robinson, M. D., D. J. McCarthy and G. K. Smyth (2010). edgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139-140.

Keywords: Microalgae, Biofuels, phosphoproteome, Lipids, Chlorella vulgaris

Received: 04 Oct 2017; Accepted: 12 Feb 2018.

Edited by:

Qiang Wang, Institute of Hydrobiology (CAS), China

Reviewed by:

Xianhai Zeng, Xiamen University, China
Yingchun Wang, Chinese Academy of Sciences, China  

Copyright: © 2018 Guarnieri, Gerritsen, Henard and Knoshaug. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

* Correspondence: Dr. Michael T. Guarnieri, National Renewable Energy Laboratory (DOE), 15013 Denver West Parkway, Golden, Golden, 80401, CO, United States,