Skip to main content

Metabolomic, proteomic and lactylated proteomic analyses indicate lactate plays important roles in maintaining energy and C:N homeostasis in Phaeodactylum tricornutum



Phaeodactylum tricornutum accumulates lipids while the growth also increases under high CO2, shedding light on its potential application in the reduction of CO2 emissions and at the same time acquiring biodiesel raw materials. However, the sensing and transducing of high C:N signals and the related response mechanism(s) remained unknown.


In this study, a multiple omics analysis was performed with P. tricornutum under low nitrogen (LN) and high CO2 (HC) conditions. The results indicated that 2-oxoglutarate was significantly increased under both LN and HC. Meanwhile, proteins involved in carbon concentration mechanism decreased, indicated that 2-oxoglutarate might regulate C:N balance through suppressing carbon fixation. Lactate, which acts in energy metabolism, signal transduction and ‘LactoylLys’ modification on proteins, was the most upregulated metabolite under both LN and HC conditions. Meanwhile, proteins involved in carbon, nitrogen and energy metabolisms were significantly regulated. Western blotting analysis suggested that non-histone L-lactylation modification was enhanced under LN and HC. Moreover, lactylated proteins were enriched in photosynthesis, central carbon metabolism, nitrogen metabolism, fatty acid synthesis and oxidative phosphorylation.


It is suggested that lactate might play important roles in energy homeostatic maintenance and C:N balance regulation in P. tricornutum through protein lactylation modification.


Diatoms are important phytoplankton which are responsible for approximately 20% of global primary productivity [1, 2] and thus play an important role in carbon fixation and the material cycle on earth [3, 4]. The great success of diatoms is mainly due to their rapid response to fluctuations in nitrate in the ocean. Diatoms possess a complete urea cycle, which might facilitate the reassimilation of nitrogen in catabolic N-containing components [5], making them response quickly to nitrogen-limited conditions. In contrast to carbon regulation in green algae and higher plants, which increase carbon stores under nitrogen depletion conditions, diatoms accelerate glycolysis and the TCA cycle to provide carbon skeletons for nitrogen redistribution [6]. This response of diatoms under nitrogen-limited conditions is a result of the superior efficiency of intracellular carbon store breakdown as a source of carbon for the reassimilation of nitrogen than for photosynthesis [6]. Because of the importance of diatoms and their unique characteristics, the response of diatoms to nitrogen depletion has received much attention in the past decade.

The pennate diatom Phaeodactylum tricornutum (P. tricornutum) is considered a model organism for the study of physiology, evolution and biochemistry in single-cell organism due to its clear genomic background [7], universal molecular toolbox [8] and stable transgene expression system [9, 10]. In the past decade, the physiological response [11], gene expression and metabolic pathway regulation mechanisms of P. tricornutum under nitrogen deficiency have been studied [12,13,14,15]. Under nitrogen deficiency, P. tricornutum accumulates lipids at the expense of depressed growth, photosynthesis, chlorophyll biosynthesis, the Calvin cycle and chrysolaminarin biosynthesis, and increased recycling of N compounds, such as amino acids, proteins, chlorophyll and nucleic acids; high turnover of cell components (e.g. soluble sugar, membrane lipids and phosphatidylcholines) that reallocates cellular nutrients to glycolysis and the TCA cycle, which are accelerated to provide carbon skeletons for redistribution of nitrogen, was all observed [12, 16, 17]. However, the mechanisms of sensing and transducing nitrogen deficiency signals and the related responses have not been elucidated.

In cyanobacteria, 2-oxoglutarate (2-OG), the metabolic intermediate of TCA cycle, played an important role in the C:N balance regulation by binding to the key regulators PII and NtcA, two receptors in the nitrogen starvation signalling pathways. Our previous study indicated that 2-OG content was upregulated 22-fold in P. tricornutum under nitrogen-deficient conditions [18]. As an important metabolite of central carbon metabolism, 2-OG provides a carbon skeleton for the assimilation of nitrogen through glutamate synthase, connecting intracellular carbon metabolism and nitrogen metabolism [19, 20]. Thus, the accumulation of 2-OG indicates a high carbon: nitrogen (C:N) ratio, and accumulated 2-OG emits a nitrogen starvation signal, playing an important role in the regulation of carbon and nitrogen metabolism in cyanobacteria [21]. When the C:N ratio is high, 2-OG suppresses carbon fixation by downregulating genes involved in the Calvin cycle and the carbon concentration mechanism (CCM), which leads to elevated CO2 concentrations near Rubisco to enhance CO2 fixation [20, 21]. The accumulation of 2-OG in P. tricornutum under nitrogen-deficient conditions, together with changes in carbon metabolism in response to nitrogen starvation in diatoms bears a greater resemblance to the response of cyanobacteria than to that of higher plants or green algae, triggering our interest in investigating the role of 2-OG in P. tricornutum.

If 2-OG indeed plays a role in C:N balance regulation in P. tricornutum, it might also act under other conditions which can result in a high intracellular C:N ratio. Besides nitrogen depletion, high inorganic carbon concentration, i.e. ocean acidification induced by increasing of atmospheric CO2, can also result in a high intracellular C:N ratio. In fact, our previous study suggests some similar response of P. tricornutum to high CO2 concentration as to low nitrogen, including accumulation of lipid, increase in glycolysis and TCA cycle activities and remodelling of cellular components [22, 23]. We propose that nitrogen deficiency or high CO2 concentration essentially disrupts the intracellular C:N balance and energy balance, thus triggering a series of signal transduction and regulation processes to maintain cell homeostasis. In contrast to cyanobacteria, P. tricornutum is an oleaginous microalga in which lipids accumulate under high C:N ratio conditions. This lipid accumulation is accompanied with remodelling of intracellular components and redistributing of metabolites. Hence, carbon and nitrogen level regulation may not be limited to the carbon and nitrogen assimilation; it may also be related to the degradation and redirection of intracellular catabolic products.

Although researchers have investigated the response to nitrogen deficiency in P. tricornutum at various levels [13,14,15,16,17], no studies have been conducted on signal sensing and transduction or their regulatory pathways. In this study, we characterized the proteome, and metabolome of P. tricornutum under normal (NC), low nitrogen (LN) and high CO2 (HC) conditions. Western blotting was performed to assess lactylation in proteins. Lactylated proteomic analysis for NC was conducted, revealing the probable roles of lactylation in P. tricornutum. The results reveal that lactate might play an important role in C:N balance regulation and energy homeostasis maintaining in P. tricornutum through signal transduction as well as protein lactylation modification.


Change in intracellular metabolite pools in response to HC and LN conditions

Table 1 shows the ratios of key cellular metabolites detected by NMR under different culture conditions (NC, LN and HC). In total, 41 metabolites, including alcohols (1), amines (2), amino acids and their derivatives (23), ammonium compounds (3), food and drug compounds (1), nucleic acid components (3), organic acids (7) and sugar (1), were identified. Most metabolites in LN and HC cultures showed similar regulation patterns, among which lactate was the most significantly upregulated metabolite, under both LN and HC conditions, compared to the NC condition, while 2-OG ranked second and third among significantly upregulated metabolites under LN and HC conditions, respectively. Of the 23 amino acids and their derivatives, the levels of two (Phe and Tyr) did not change; the levels of two (Gly and sarcosine) were increased, and the levels of the remaining 19 were decreased under LN and HC conditions, indicating N depletion. The levels of Gly and sarcosine were increased, indicating their possible importance in N redistribution. The levels of all the nucleic acid components detected were decreased, confirming N depletion. The levels of most organic acids were decreased, while increases in 2-OG and lactate reflected their important roles in C:N balance regulation. Glucose levels increased, likely because of increased carbohydrate catabolism. Some metabolites showed opposite expression patterns under HC and LN conditions. The levels of two amines (dimethylamine and methylamine) were increased under LN conditions and decreased under HC conditions, while the level of one ammonium compound (choline) decreased under LN conditions and increased under HC conditions.

Table 1 Ratios of abundances of key cellular metabolites under different culture conditions determined by NMR

Differential expression of proteins under HC and LN conditions

Overview of the proteomic data

Based on label-free proteome quantification technology, 5958 proteins were identified, among which 5325 proteins were quantified in P. tricornutum under NC, LN and HC conditions. Most peptides consisted of 7–20 amino acids (Additional file 6: Fig. S1a), which conforms to general rules based on enzymatic hydrolysis and mass spectrometry fragmentation. Most proteins corresponded to more than two peptides (Additional file 6: Fig. S1b). The coverage of most proteins was less than 30% (Additional file 6: Fig. S1c). The molecular weights of the identified proteins were evenly distributed (Additional file 6: Fig. S1d). A full range of functional annotations of the identified proteins was obtained on the basis of Gene Ontology (GO), protein domain, KEGG pathway, COG functional classification and subcellular structure positioning. The details of the annotations results are presented in Additional file 2: Table S1.

Identification and enrichment analysis of differentially expressed proteins (DEPs)

To understand the proteome responses to C:N ratio fluctuation in P. tricornutum, proteomic analysis based on label-free proteome quantification technology was performed to identify the DEPs in P. tricornutum under NC, LN and HC conditions. According to the criteria of a 1.5-fold increase (or 0.67-fold decrease) and a P value < 0.05 [24], a total of 1471 DEPs, including 694 differentially upregulated proteins (DUPs) and 777 differentially downregulated proteins (DDPs), were identified under LN conditions (Additional file 3: Table S2), and 327 DEPs, including 215 DUPs and 112 DDPs, were identified under HC conditions (Additional file 4: Table S3).

The KEGG pathway-based enrichment analysis revealed significant enrichment of DEPs in pathways involved in carbon metabolism, nitrogen metabolism, lipid metabolism, protein and amino acids metabolism, energy metabolism and signalling, including the enrichment of DUPs under LN conditions in TCA cycle, nitrogen metabolism, phenylpropanoid biosynthesis, phosphatidylinositol signalling system, glycerophospholipid metabolism and oxidative phosphorylation (Fig. 1a); the enrichment of DDPs under LN conditions in photosynthesis, biosynthesis of amino acids, carbon fixation in photosynthetic organisms and nitrogen metabolism (Fig. 1a), the enrichment of DUPs under HC conditions in photosynthesis—antenna proteins, phenylpropanoid biosynthesis and glycolysis / gluconeogenesis (Fig. 1b), and the enrichment of DDPs under HC conditions in biosynthesis of amino acids, pyrimidine metabolism and nitrogen metabolism (Fig. 1b).

Fig. 1
figure 1

Enriched KEGG pathways analysis of differentially expressed proteins (all P < 0.05). a LN/NC. b HC/NC. Red, pathways terms of upregulated proteins. Green, pathways terms of downregulated proteins. X axis, − log10 (Fisher' exact test P value)

Similar expression patterns of DEPs under LN and HC conditions indicated their functions in C:N balance regulation

As showed in Table 2, the third most upregulated DEP was predicted to be a transcription factor (UniProt B7G9R3), with fold changes of 44.9 and 17.1 under LN and HC conditions, respectively. The two most downregulated DEPs under LN conditions, UniProt B7G9B2 (0.0024-fold change) and UniProt B7G9B0 (0.0062-fold change), are two cell surface proteins and were found to be the second and fourth most downregulated proteins under HC conditions with a 0.0192-fold change and 0.0296-fold change, respectively. The third and fourth most downregulated DEPs under LN conditions, UniProt B7G9B1 (0.0115-fold change) and UniProt B7FYL2 (0.0145-fold change), are both activated upon iron starvation [7]. They were also downregulated under HC conditions, ranking first and seventh among DDPs with a 0.0114-fold change and 0.0859-fold change, respectively. Many DEPs (20 under LN and 2 under HC conditions) are members of the solute carrier family (Additional file 3: Table S2 and Additional file 4: Table S3), indicating a frequent exchange of metabolites. In addition, many DEPs are transcription factors or involved in signal transduction (Additional file 3: Table S2 and Additional file 4: Table S3), indicating that in P. tricornutum, the C:N balance might be regulated through signalling pathways. This similar expression pattern of DEPs under LN and HC conditions indicated they might play a role in C:N balance regulation.

Table 2 Different expressed proteins under LN/NC and HC/NC

DEPs involved in phenylpropanoid biosynthesis and nitrogen metabolism

The most upregulated DEP under both LN (790.8-fold change) and HC (287.8-fold change) was predicted to be peroxiredoxin (Tables S2 & S3), which is involved in phenylpropanoid biosynthesis. The second and fifth (LN) or fourth (HC) most upregulated DEPs were predicted to be transcription factors (UniProt B7GDM6 and UniProt B7FUS4) involved in the regulation of phenylpropanoid biosynthesis, with fold changes of 167.6 and 26.4 under LN conditions and 50.0 and 14.4 under HC conditions (Tables S2 & S3), respectively. A total of 11 DEPs were involved in nitrogen metabolism under LN, while those of HC were three (Table 2, Fig. 2), indicating that though P. tricornutum experienced nitrogen deficiency stress under both LN and HC conditions, it is more severe under LN.

Fig. 2
figure 2

Different expressed proteins involved in photosynthesis, carbon concentration mechanism (CCM), glycolysis, Calvin cycle, TCA cycle, pentose phosphate pathway, and oxidation phosphorylation. The number in the circle indicates the no. of proteins listed in Table 2 where the full information for the proteins can be obtained. Square, LN/NC; circle, HC/NC. Red, upregulation of protein content. Green, downregulation of protein contents. Blue, the regulation mood of protein content. Grey, proteins do not detect in the experiment group (LN or HC)

DEPs involved in photosynthesis, carbon concentration mechanism (CCM), central carbon metabolism and chrysolaminarin degradation

Two DUPs and four DDPs in LN were light-harvesting complex-related proteins. A total of 18 DDPs and two DUPs under LN were involved in photosynthesis, indicating a severe downregulated of photosynthesis under LN (Table 2, Fig. 2). Five light-harvesting complex-related proteins were upregulated under HC condition, indicating enhance of photosynthesis under HC (Table 2, Fig. 2). In addition, proteins involved in the CCM were significantly downregulated under LN conditions, including two Na+-independent Cl/HCO3 exchangers (UniProt B7S435 and B7S437); four carbonic anhydrases (CA, UniProt B7FNU0, B7GA80, B7FXP8, and B5Y401); pyruvate, orthophosphate dikinase (PPDK, UniProt B7G585); and phosphoenolpyruvate carboxykinase (PEPCK, UniProt B7GA05) (Table 2, Fig. 2). A CA (B7FNT2) and a PPDK (UniProt B7G585) were significantly downregulated under HC condition (Table 2, Fig. 2). Under LN, proteins involved in Calvin cycle, glycolysis, and pentose phosphate pathway were mainly downregulated, while that of TCA cycle increased (Table 2, Fig. 2). Under HC there were not significantly change in most proteins involved in the central carbon metabolism, except that some proteins involved in glycolysis increased. Diatoms store carbon in the form of 1,3-β-d-glucan (chrysolaminarin) or lipids. We observed an increase in β-glucosidase under both LN (UniProt B7GBX3, UniProt B7FWJ8 and UniProt B7GB76) and HC (UniProt B7FWJ8) conditions (Table 2), indicating the degradation of chrysolaminarin.

DEPs involved in phosphatidylinositol signalling system

Five proteins involved in phosphatidylinositol signalling system, including diacylglycerol kinase (DGK), 1-phosphatidylinositol-4-phosphate 5-kinase (PIP5K), phosphatidylinositol phospholipase C (PLCD), calmodulin, phosphatidylinositol 4-kinase B (PI4KB) and inositol polyphosphate 5-phosphatase (INPP5B/F) were all significantly up-regulated under LN condition (Table 2), indicating phosphatidylinositol signalling system might play an important role in nitrogen-limited response.

DEPs involved in lactate metabolism and lysine lactylation

D-lactate dehydrogenase (ldhA), which catalyses the dehydration of pyruvate to generate D-lactate, was upregulated under both LN (5.6-fold change) and HC (4.5-fold change) conditions, a finding consistent with lactate accumulation under both LN and HC conditions. Besides, a l-lactate permease (LCTP), which was reported to expressed under nitrogen-limited conditions [7], and enzymes related to D-lactylation, namely, lactolglutathione lyase (GLXI) and glyoxalase (GLO1 and GLO2) [25], were also detected in the P. tricornutum proteome analysis, although their levels did not vary significantly upon different treatments (Table 2).

Protein validation by PRM and biochemical analysis

We performed PRM to eliminate most of the possible errors though the label-free data in this report had undergone rigorous statistical and bioinformatics analysis. Analysis on FbaC5, ZEP3, bglX, nirA, PPdK, NR, PCCA, LDHA, GLO_2, GLO_1, Formidase, PDH1, PEPCK1, GltX and acnBto was performed to validate the label-free results. The results of target proteins detected by PRM were consistent with the label-free analysis (Fig. 3), which indicated that our proteomics results based on label free were highly reliable and reproducible.

Fig. 3
figure 3

The comparison of protein expression by label free and PRM. a LN/NC. b HC/NC. FbaC5 fructose-bisphosphate aldolase, ZEP3 zeaxanthin epoxidase, bglX beta-glucosidase, nirA ferredoxin-nitrite reductase, PPdK pyruvate, orthophosphate dikinase, NR nitrate reductase, PCCA carboxylase propionyl-CoA carboxylase, LDHA D-lactate dehydrogenase, GLO_2 glyoxalase OS, GLO_1 glyoxalase, PDH1 precursor of dehydrogenase pyruvate dehydrogenase E1, PEPCK1 phosphoenolpyruvate carboxykinase, GltX ferredoxin-dependent glutamate synthase, acnB aconitate hydratase 2 / 2-methylisocitrate dehydratase

Pan-acetylation and pan-lactyllysine Western blotting analyses

The overall regulation of acetylation and lactylation patterns of all proteins was determined by Western blotting using a pan anti-acetyllysine antibody and a pan-lactyllysine antibody, respectively. As shown in Fig. 4, many proteins were either acetylated or lactylated under NC, LN and HC conditions, but the overall acetylation and lactylation levels differed. Increased lactylation was observed on larger proteins under LN and HC conditions, consistent with the increase in lactate content implied by the metabolomic analysis results. Under LN and HC conditions, increased acetylation was observed on proteins with a molecular weight of 12–15 kDa, which were most likely members of the histone family. These findings confirmed that protein lactylation proceeds through temporal dynamics that differ from those of protein acetylation.

Fig. 4
figure 4

Pan-acetylation and pan-lactyllysine Western blotting analyses. a Coomassie brilliant blue staining. b Western blotting with pan anti-acetyllysine antibody. c Western blotting with pan anti-lactyllysine antibody

Lactylated proteomic analysis under NC in P. tricornutum

Global profiling of lactylated proteome

To study the probable function of lysine lactylation in P. tricornutum, we performed global profiling of the lysine lactylome for NC. A total of 463 lactyllysine sites across 212 proteins were identified (Additional file 5: Table S4). To further understand the biological regulations and functions of lactylated proteins in C/N balance in P. tricornutum, we carried out Clusters of Orthologous Groups (COG) annotation and KEGG pathway enrichment analysis. The results of COG functional annotation analysis showed that various lactylated proteins were engaged with crucial biological processes, such as posttranslational modification, protein turnover and chaperones; carbohydrate transport and metabolism; translation, ribosomal structure and biogenesis; energy production and conversion; amino acid transport and metabolism; and lipid transport and metabolism (Fig. 5).

Fig. 5
figure 5

Clusters of Orthologous Groups (COG) annotation for lactylated proteins

Lactylated proteins are enriched on photosynthesis, central carbon metabolism and fatty acid biosynthesis pathways

Protein domain analysis revealed that lactylated proteins were concentrated in ATP synthase, GADPH, core histone H2A/H2B/H3/H4, phosphoglycerate kinase, fructose-bisphosphate aldolase, carbonic anhydrase, enolase and so on (Fig. 6a). KEGG pathway enrichment analysis illustrated that lactylated proteins were concentrated in photosynthesis, carbon fixation in photosynthetic organisms, photosynthesis-antenna proteins, glycolysis/gluconeogenesis, oxidative phosphorylation, fatty acid biosynthesis, fructose and mannose metabolism and pentose phosphate pathway (Fig. 6b), indicating that lactylation played an important role in energy balance and carbon metabolism.

Fig. 6
figure 6

Protein domain enrichment analysis (a) and KEGG pathway enrichment analysis (b) for lactylated proteins

Lactylated proteins involved in glycolysis/gluconeogenesis, Calvin cycle, pentose phosphate pathway and fatty acid synthesis are shown in Fig. 7. Besides, nitrogen metabolism-related proteins, including ammonium transporter, glutamine synthetase and carbonic anhydrase, were also detected (Fig. 7), indicating that lactylation also played an important role in nitrogen metabolism. Simultaneously, components of oxidative phosphorylation and photosynthetic phosphorylation were also detected, indicating that lactylation participates in energy metabolism.

Fig. 7
figure 7

Lactylated proteins involved in photosynthesis, glycolysis, Calvin cycle, TCA cycle, pentose phosphate pathway, fatty acid synthesis and oxidation phosphorylation. FBA fructose-bisphosphate aldolase, TPI triosephosphate isomerase, GAPDH glyceraldehyde 3-phosphate dehydrogenase (phosphorylating), GAP d-glyceraldehyde 3-phosphate, PGK phosphoglycerate kinase, BPG 3-phospho-d-glyceroyl phosphate, 3PG 3-phospho-d-glycerate, 2PG phospho-d-glycerate, PGAM 2,3-bisphosphoglycerate-dependent phosphoglycerate mutase, ENO enolase, X5P d-xylulose 5-phosphate, TKL transketolase, E4P d-erythrose 4-phosphate, S7P sedoheptulose 7-phosphate, SBP sedoheptulose 1,7-bisphosphate, PRK phosphoribulokinase, MCOA malonyl-CoA, ACACA acetyl-CoA carboxylase / biotin carboxylase 1, Malonyl-[acp] malonyl-[acyl-carrier protein], FabD [acyl-carrier-protein] S-malonyltransferase, ACSL long-chain acyl-CoA synthetase, FabF 3-oxoacyl-[acyl-carrier-protein] synthase II, FabI enoyl-[acyl-carrier protein] reductase I, CA carbonic anhydrase, AMT ammonium transporter, LHCA light-harvesting complex I chlorophyll a/b binding protein 1, ATPF1A F-type H+/Na+-transporting ATPase subunit alpha, ATPF1B F-type H+/Na+-transporting ATPase subunit beta, ATPF1G F-type H+-transporting ATPase subunit gamma, ATPF0B F-type H+-transporting ATPase subunit b, NUOA NADH-quinone oxidoreductase subunit A, NDUFAB1 NADH dehydrogenase (ubiquinone) 1 alpha/beta subcomplex 1, acyl-carrier protein, ATPF0A F-type H+-transporting ATPase subunit a, PPA inorganic pyrophosphatase, ATPeF0O F-type H + -transporting ATPase subunit O, ATPeV1A V-type H+-transporting ATPase subunit A

Interestingly, lactylation mainly occurred in proteins that located in cytoplasm and chloroplast, but rarely in proteins located in mitochondria, implying that protein lactylation modifications are affected by cell compartmentalization. It is probable that there is lactate permease on chloroplast membrane, but not on mitochondria membrane.

Influence of lactate on growth and total lipid content of P. tricornutum

As shown in Fig. 8, lactate at 1 mM, 5 mM and 10 mM could increase (P < 0.05) growth of P. tricornutum, with 5 mM being the most significant. Total lipid content was determined to verify if lactate influenced total lipid content. The results showed that lactate can significantly (P < 0.05) promote the accumulation of lipids in P. tricornutum, and the total lipid content increases with the increase of lactate concentration, indicating that lactate might play an important role in lipid metabolism in P. tricornutum. Yet as lactate can also function as a carbon source, further studies, such as 13C labelling, are needed to verify the function of lactate in P. tricornutum.

Fig. 8
figure 8

Influence of lactate on growth (a) and total lipid content (b) of P. tricornutum


2-OG might regulate P. tricornutum C:N balance by suppressing CCM

In recent decades, the response of P. tricornutum to nitrogen deficiency has received considerable attention. Research mainly concerns the physiological response, gene expression and metabolic pathway shift or lipid accumulation process [12,13,14, 26]. However, the mechanisms behind P. tricornutum responses, e.g. sensing nitrogen deficiency and transmitting relevant signals to induce the response mechanism, have not been elucidated. To some extent, nitrogen deficiency (LN) and a high CO2 concentration (HC) indicate a high C:N ratio [23, 27], which might trigger the activation of a series of signalling pathways. For example, 2-OG indicates a high C:N ratio and plays an important role in the regulation of carbon and nitrogen metabolism in cyanobacteria [21]. By depressing the expression of genes involved in CCM, 2-OG represses carbon fixation and thus maintains the intracellular C:N balance. In this study, we observed 14-fold and 12-fold upregulation of 2-OG expression under LN and HC conditions (Table 1). In addition, proteins involved in the CCM were significantly downregulated under LN conditions (Table 2). The accumulation of 2-OG and the downregulation of CCM-related genes under both LN and HC conditions indicated that 2-OG might play an important role in sensing and transducing high C:N ratio signals and act as a ligand that interacts with carbon fixation regulators. This supposition is consistent with the idea that carbon metabolism induced in diatoms in response to nitrogen starvation bears closer resemblance to the response of cyanobacteria than that of higher plants or green algae [6].

Probable role of lactate in P. tricornutum C:N balance regulation

In addition to an increase in 2-OG, we also observed a dramatic increase in lactate content under both LN and HC conditions. Lactate is used to be considered as a metabolism waste by-product of glycolysis. Currently, lactate is considered a crucial regulatory nexus for energy metabolism [28], a signalling molecule participating in a series of metabolic processes [29, 30] and a substrate for posttranslational modifications (PTMs) through the ‘LactoylLys’ modification of proteins [25, 31]. The increase in lactate under both LN and HC conditions indicated its likely roles in the following processes:

  1. 1.

    Regulation of energy balance. In organisms, the main source of ATP is oxidative phosphorylation, in which electrons in NADH (produced by cell metabolism) are delivered to oxygen via the respiratory chain with the concomitant production of ATP and NAD+. Under hypoxic conditions, oxidative phosphorylation is repressed, NADH accumulates and NAD+ regeneration is blocked. To maintain the energy supply under hypoxic conditions, plants switch from respiration to fermentative metabolism by activating lactate dehydrogenase, which catalyses the reduction of pyruvate to generate lactate, accompanied by the consumption of NADH. Thus, lactate dehydrogenase can regulate cellular energy balance. In this study, we observed a significant decrease in AMP under both LN and HC conditions (Table 1), indicating a high ATP:AMP ratio; e.g., cells are in high energy mode. This outcome might be a result of a decrease in the NADH consumption process, i.e. nitrogen reduction. We observed a decrease expression in nitrate reductase (UniProt B7G997) and ferredoxin-nitrite reductase (UnitProt B7G0L4) under both LN and HC conditions (Table 2, Fig. 2). These results might indicate decreased NADH consumption. Furthermore, processes in which NADH is produced, i.e., the TCA cycle under LN conditions and pyruvate dehydrogenase and photosynthesis under HC conditions, were accelerated. This acceleration might disrupt the energy balance and activate lactate dehydrogenase, resulting in the accumulation of lactate. In this study, we observed a 5.57-fold and 4.46-fold increase in lactate dehydrogenase (ldhA, B7S4E4) under LN and HC conditions (Table2), indicating the upregulation of lactate reductant from pyruvate. Moreover, the shuttling of L-lactate between different organs and cells creates a major circulatory mechanism of carbohydrate sources and NADH/NAD+ [32]. In this study, we also detected L-lactate permease (LCTP, Gene ID: 7,195,398) under LN condition (Table 2 and Additional file 2: Table S1), suggesting the important role of lactate in the C:N balance and energy balance regulation.

  2. 2.

    Regulation of lipid metabolism. In contrast to cyanobacteria, P. tricornutum is an oleaginous microalga that accumulates lipids under nitrogen-limited conditions. In this process, P. tricornutum remodels intracellular components and redistributes metabolites [14]. Therefore, its C:N balance regulation may not be limited to the assimilation process of carbon and nitrogen and may also be related to the degradation and redirection of intracellular catabolic products. Lactate itself can act as a signal molecule in the cAMP pathway and regulate a series of cellular processes. In mammalian adipose tissue, for example, lactate can combine with its sensor, an orphan G protein-coupled receptor, activating it and suppressing lipolysis in fat cells [30]. Despite its categorization as a photosynthetic organism, P. tricornutum has an unusual evolutionary history, with secondary endosymbiotic origin and mosaic genome that contains ‘animal-like’, ‘plant-like’ and ‘bacteria-like’ genes [33,34,35,36,37]. Some of its genes are more similar to those of animals than to their photosynthetic counterparts [38]. Recently, it was reported that in P. tricornutum, phosphatidylinositol 3-kinase (PI3K), which is involved in lipid accumulation, is more closely associated with mammalian homologues than those in higher plants [39]. The accumulation of both lactate and lipids in P. tricornutum under LN and HC conditions indicated that lactate might also suppress lipolysis in P. tricornutum, as it does in mammalian cells, and thus facilitate lipid accumulation in P. tricornutum under LN and HC conditions.

  3. 3.

    ‘LactoylLys’ modification on proteins (lactylation). Researchers recently found that lactate can regulate protein functions through posttranslational modifications (PTMs) by conjugating to protein Lys residues and generating a ‘LactoylLys’ modification on proteins, revealing non-metabolic functions of lactate [31]. Protein lactylation is a new topic currently studied mainly in mammals, as protein lactylation plays an important role in diseases, such as neoplasia, sepsis and autoimmune diseases. Studies related to protein lactylation in algae have not been reported. Protein lactylation is positively correlated with intracellular lactate concentration [25]. In this study, lactate was the most significantly upregulated metabolite under both LN and HC conditions. The WB results (Fig. 4) indicated proteins in P. tricornutum had undergone lactylation, which was particularly evident under LN and HC conditions. Lactylated proteomics under NC condition revealed that lactylated proteins were significantly involved in photosynthesis, oxidative phosphorylation, central carbon metabolism including Calvin cycle, glycolysis and pentose phosphate pathway, fatty acid synthesis and nitrogen metabolism, confirming protein lactylation in P. tricornutum played an important role in C:N balance regulation, energy homeostasis maintenance and fatty acid synthesis.

Probable role of the IP3K signalling pathway under LN conditions

In this study, we found that genes involved in the phosphatidylinositol signalling system were upregulated at the protein level (Table 2) under LN conditions. It has been reported that in P. tricornutum, the expression of genes involved in MAPK signalling was upregulated, while that in lipid metabolism was decreased under treatment with a phosphatidylinositol 3-kinase inhibitor, demonstrating a role for phosphatidylinositol 3-kinase in the allocation of carbon, including as a factor of lipid reduction in P. tricornutum [39]. In a study of the green alga Chlamydomonas reinhardtii, it was reported that the phosphatidylinositol 3-kinase signalling system is a master regulator of energy and carbon metabolism [40], regulating the homeostasis of membrane lipids, TAGs, starch, free fatty acids (FFAs) and ATP by influencing other components of the signalling network. Here, we also found that the expression of a carboxypeptidase (B5Y4V6), which negatively regulates membrane lipid metabolism, was upregulated under both LN and HC conditions. The enrichment of upregulated proteins in the phosphatidylinositol signalling system and related genes indicated that this pathway might play important regulatory roles (energy and carbon balance regulation and homeostasis maintenance) in the P. tricornutum response to LN conditions. In addition, genes involved in autophagy, plant–pathogen interactions, endocytosis and phagosomes were also enriched, indicating a complex regulatory pattern through signal transduction in P. tricornutum under LN conditions.

Probable role of phenylpropanoids in N recycling under LN and HC conditions

The most upregulated DUP under both LN (790.8-fold change) and HC (287.8-fold change) conditions was PRDX6, which is involved in phenylpropanoid biosynthesis, followed by two predicted transcription factors involved in the regulation of phenylpropanoid biosynthesis ranked second and fifth (LN) and second and fourth (HC) (Tables S2 & S3), indicating an important role for the phenylpropanoid biosynthesis pathway under LN and HC conditions. The phenylpropanoid biosynthesis pathway is involved in the degradation of phenylalanine and tyrosine and thus produces abundant NH4+, facilitating the recycling of nitrogen under N-limited conditions, while nitrate inhibits large sectors of phenylpropanoid metabolism [41]. The upregulation of the phenylpropanoid biosynthesis pathway might facilitate nitrogen reallocation under both LN and HC conditions.

Overall response to a high C:N ratio in P. tricornutum under HC and LN conditions

Considering the combined results of the proteome, metabolome and WB analysis and lactylated proteome with respect to P. tricornutum under HC and LN conditions, we propose the following responses of P. tricornutum to C:N imbalance (Fig. 9): a high C:N ratio stimulates 2-OG accumulation. The resulting nitrogen reduction and nitrite reduction were reduced as a result of N limitation, thereby reducing the consumption of NADH. In addition, accelerated glycolysis, photosynthesis (HC conditions) and the TCA cycle (LN conditions) result in increases in NADH and ATP (Fig. 2). A low AMP:ATP ratio (Table 1) depresses the activity of NADH dehydrogenase, and the accumulation of NADH triggers the activity of ldhA (Table 2), resulting in the accumulation of lactate (Table 1). 2-OG and lactate then serve as signalling molecules and regulate a series of pathways. 2-OG senses and transduces signals when the C:N ratio is high, regulates genes involved in the CCM and reduces the extent carbon fixation (Table 2, Fig. 2). Moreover 2-OG can act on the PHD and E3 ubiquitin ligation complex and induce protein degradation through ubiquitin-mediated proteolysis, providing carbon skeletons and NH4+ for the reallocation of nitrogen. Lactate can regulate the function of proteins through PTM and regulate cellular components through a series of signalling pathways, i.e. IP3K signalling pathways. Lactylated proteins were involved in central carbon metabolism, photosynthesis, oxidative phosphorylation and nitrogen metabolism, thus can maintain C:N balance and energy homeostasis.

Fig. 9
figure 9

Overall response to high C/N ratio in P. tricornutum under HC and LN conditions. Red, upregulation of metabolites or protein content of genes. Green, downregulation of metabolites or protein contents of genes. Blue, the regulation mood of metabolites or protein content of genes remain unknown. Full line, the regulation of metabolites or proteins are verified by omics data. Dashes lines, the regulation of metabolites or proteins are speculated by related metabolites or proteins. NDUFA12 NADH dehydrogenase (ubiquinone) 1 alpha subcomplex subunit 12, NR nitrate reductase (NAD(P)H), nirA ferredoxin-nitrite reductase, ldhA d-lactate dehydrogenase. PYR pyruvate, CCM carbon concentration mechanism


Multiple omics analysis was conducted in P. tricornutum under HC and LN conditions. The results indicated 2-OG and lactate was significantly up-regulated under both LN and HC conditions. Genes involved in CCM were down-regulated, indicating 2-OG regulated C:N balance by suppressing CCM-related genes. Genes involved in lactate redundant and energy metabolism are significant regulated. By WB analysis it was proved that lactylation modification occurred in P. tricornutum. What is more, the non-histone lactylation modification was enhanced under LN and HC. Lactylated proteomic analysis revealing that lactylation might play an important role in energy and carbon balance regulation and homeostasis maintenance in P. tricornutum. Our result reveal C:N regulation mechanism under LN and HC conditions, shedding light on the potential application of P. tricornutum to the reduction of CO2 emissions from industrial flue gas streams and at the same time acquiring biodiesel raw materials.

Materials and methods

Cell culture and treatments

P. tricornutum Bohlin was obtained from the Microalgae Culture Center (MACC) at Ocean University of China. Algal cells were cultured using NaHCO3-free sterilized artificial seawater supplemented with f/2 [42] under cool white fluorescent light (~ 100 μmol m–2 s–1) at 20 °C and a 12:12 dark:light cycle. A UV/visible spectrophotometer (UV-1800, Shimadzu, Japan) was used to measure the absorbance at 730 nm (A730 nm) to monitor cell growth. For high CO2 (HC, ~ 2000 μatm) or normal CO2 (NC, ~ 400 μatm) treatment, the culture was bubbled with ~ 2,000 μatm of CO2 or ambient air with a flow rate of 500 ml min–1. For the low nitrogen (LN) treatment, a nitrogen-free f/2 medium was produced with artificial seawater and the cultures were bubbled with ambient air at a flow rate of 500 ml min–1. For experiments testing the influence of lactate on growth and lipid content, ( ±)-lactate (Xilong Chemical Co., Ltd., Shantou, Guangdong province, China) was added to a final concentration of 1, 5 and 10 mM. Each treatment was administered in triplicate sets. After 8 days, cell pellets were collected and then washed using distilled water at 5,000 g for 4 min. The pellets were frozen in liquid nitrogen and stored at − 80 °C.

Lipid analysis and metabolomic analysis

Total lipids were extracted and weighed as described by previously [43] with minor modifications. Metabolomic analysis was performed as described [18]. Full details are given in Additional file 1: Methods S1.

Proteomic analysis

Algal cells were ground into powder in liquid nitrogen, and proteins were extracted using phenol buffer (containing 10 mM dithiothreitol and a 1% protease inhibitor) and lysed by ultrasound. After adding Tris-balanced phenol and centrifuging (4℃, 5500g, 10 min), the supernatant was maintained overnight to allow precipitation with ammonium acetate/methanol, and then, it was washed with methanol and acetone. The protein precipitate was dissolved with urea, and the concentration was determined with a BCA protein assay kit (Beyotime). Equal amounts of protein were hydrolysed by trypsin and analysed by LC MS/MS as described in Additional file 1: Methods S1.

The MS/MS raw data were administered by operating MaxQuant ( with an integrated Andromeda search engine (v1.6.15.0). Full detail parameters are given in Additional file 1: Methods S1. Significantly differentially expressed proteins were defined by the criteria of a fold change (FC) ≥ 1.5 or ≤ 0.67 (P < 0.05). Functional annotation of the proteins was performed using non-redundant protein (NR), Swiss Protein (SwissProt), Gene Ontology (GO) and KEGG databases.

Protein validation by parallel reaction monitoring (PRM)

To verify the protein expression levels obtained by label-free analysis, the expression levels of 15 selected proteins, including fructose-bisphosphate aldolase (FbaC5), zeaxanthin epoxidase (ZEP3), beta-glucosidase (bgl), ferredoxin-nitrite reductase (nirA), pyruvate, orthophosphate dikinase (PPdK), nitrate reductase (NR), propionyl-CoA carboxylase (PCCA), d-lactate dehydrogenase (LDHA), glyoxalase (GLO_2 and GLO_1), formidase, precursor of dehydrogenase pyruvate dehydrogenase E1 (PDH1), phosphoenolpyruvate carboxykinase (PEPCK1), ferredoxin-dependent glutamate synthase (GltX) and aconitate hydratase 2 / 2-methylisocitrate dehydratase (acnB), were further quantified by PRM analysis using the original protein samples (NC, LN, and HC) as described in Sect. 2.3. PRM validation was performed as described previously [44]with minor revisions.

Pan anti-acetylation and pan anti-lactyllysine as determined by Western blotting

Proteins used in Western blotting were extracted according to the same operation manual as used in proteomic analysis except that 3 μM TSA (trichostatin A) and 50 mM nicotinamide were added. Protein concentration was determined by BCA reagent. Western blotting samples were prepared as follows: equal amounts of protein (20 µg/sample) were mixed with SDS-PAGE 4 × loading buffer with a ratio of 3:1 (v/v), heating for 10 min at 95 °C. The same volume of 20% protein marker, boiled samples and 4 × loading buffer were loaded into wells. Gel electrophoresis was then performed as follows: at 80 V for 30 min, then at 120 V until loading dye reach the bottom edge of SDS-PAGE gel. Following transfer the protein to nitrocellulose membrane at 100 V for 1 h at 4 °C, the membranes were blocked in 1 × TBST with 5% non-fat milk with gentle shaking for 1 h at RT. By rinsing with 1 × TBST for three times (10 min each), the membrane was incubated in primary antibody (anti-acetyllysine antibody, PTM-101, Lot: 12838533L303, 1:1000 dilution in 1 × TBST with 2.5% BSA; and pan anti-lactyllysine, PTM-1401RM, Lot: L011121, 1:500 dilution in 1 × TBST with 2.5% BSA) at 4 °C with gentle shaking overnight, rinsing with1 × TBST for 10 min three times, incubating for 1 h at RT with secondary antibody (Thermo, Pierce, goat anti-mouse IgG, H + L, peroxidase-conjugated, 31,430 and 31,460, respectively) diluted 1:10,000 in in 1 × TBST with 5% non-fat milk), rinsing with 1 × TBST for 10 min three times, and staining with ECL western blot detection reagent (Beyotime, Beijing, China) for 2 min, followed by quantification using multiple exposure time to acquire the optimal picture.

Lactylated proteomic analysis

Protein was extracted and hydrolysed by trypsin as described in 2.3, following by pan-antibody-based PTM enrichment. Briefly, tryptic peptides dissolved in NETN buffer (100 mM NaCl, 1 mM EDTA, 50 mM Tris–HCl, 0.5% NP-40, pH 8.0) were incubated with pre-washed antibody beads (Lot: L011121, PTM Bio) at 4 °C overnight with gentle shaking. After that the beads were washed for four times with NETN buffer and twice with H2O. Then the bound peptides were eluted from the beads with 0.1% trifluoroacetic acid. Finally, the eluted fractions were combined and vacuum-dried and desalted with C18 ZipTips (Millipore) according to the manufacturer’s instructions.

LC MS/MS and database search were conducted as described in 2.3. Full detail parameters are given in Methods S1. For motif analysis, an adjacent sequence model of lactyllysine was analysed by motif-x (, integrating 10 amino acids upstream and downstream of the modification sites (lactyl-21-mers). When the number of a lactyl-21-mer is greater than 20 and P < 0.000001, the lactyl-21-mer is considered to be a motif of modified peptides. Protein sequences for all P. tricornutum proteins were utilized as background databases.

Statistical analysis

The significance of differences between samples was assessed by two-tailed student t-tests and one-way analysis of variance. The statistical significance level was set at P < 0.05.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


  1. Field C, Behrenfeld M, Randerson J, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Sci. 1998;281(5374):237–40.

    Article  CAS  Google Scholar 

  2. Falkowski P, Barber R, Smetacek V. Biogeochemical controls and feedbacks on ocean primary production. Sci. 1998;281(5374):200–6.

    Article  CAS  Google Scholar 

  3. Treguer P, Nelson D, Van Bennekom A, DeMaster D, Leynaert A, Queguiner B. The silica balance in the world ocean: a reestimate. Sci. 1995;268(5209):375.

    Article  CAS  Google Scholar 

  4. Werner D. Silicate metabolism. The biology of diatoms 1977, 149.

  5. Allen AE, Dupont CL, Obornik M, Horak A, Nunes-Nesi A, McCrow JP, Zheng H, Johnson DA, Hu H, Fernie AR, et al. Evolution and metabolic significance of the urea cycle in photosynthetic diatoms. Nature. 2011;473(7346):203.

    Article  CAS  Google Scholar 

  6. Hockin NL, Mock T, Mulholland F, Kopriva S, Malin G. The response of diatom central carbon metabolism to nitrogen starvation is different from that of green algae and higher plants. Plant Physiol. 2012;158(1):299–312.

    Article  CAS  Google Scholar 

  7. Bowler C, Allen AE, Badger JH, Grimwood J, Jabbari K, Kuo A, Maheswari U, Martens C, Maumus F, Otillar RP. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature. 2008;456(7219):239–44.

    Article  CAS  Google Scholar 

  8. Siaut M, Heijde M, Mangogna M, Montsant A, Coesel S, Allen A, Manfredonia A, Falciatore A, Bowler C. Molecular toolbox for studying diatom biology in Phaeodactylum tricornutum. Gene. 2007;406(1–2):23–35.

    Article  CAS  Google Scholar 

  9. De Risco V, Raniello R, Maumus F, Rogato A, Bowler C, Falciatore A. Gene silencing in the marine diatom Phaeodactylum tricornutum. NAR. 2009;37(14): e96.

    Article  Google Scholar 

  10. Stukenberg D, Zauner S, Dell’Aquila G, Maier UG. Optimizing CRISPR/cas9 for the diatom Phaeodactylum tricornutum. Front Plant Sci. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Li F, Beardall J, Collins S, Gao K. Decreased photosynthesis and growth with reduced respiration in the model diatom Phaeodactylum tricornutum grown under elevated CO2 over 1800 generations. Glob Change Biol. 2016.

    Article  Google Scholar 

  12. Burrows EH, Bennette NB, Carrieri D, Dixon JL, Brinker A, Frada M, Baldassano SN, Falkowski PG, Dismukes GC. Dynamics of lipid biosynthesis and redistribution in the marine diatom Phaeodactylum tricornutum under nitrate deprivation. Bioenergy Res. 2012;5(4):876–85.

    Article  CAS  Google Scholar 

  13. Levitan O, Dinamarca J, Zelzion E, Lun DS, Guerra LT, Kim MK, Kim J, Van Mooy BAS, Bhattacharya D, Falkowski PG. Remodeling of intermediate metabolism in the diatom Phaeodactylum tricornutum under nitrogen stress. Proc Natl Acad Sci U S A. 2015;112(2):412–7.

    Article  CAS  Google Scholar 

  14. Remmers IM, D’Adamo S, Martens DE, de Vos RCH, Mumm R, America AHP, Cordewener JHG, Bakker LV, Peters SA, Wijffels RH, et al. Orchestration of transcriptome, proteome and metabolome in the diatom Phaeodactylum tricornutum during nitrogen limitation. Algal Res Biomass Biofuels Bioprod. 2018;35:33–49.

    Google Scholar 

  15. Yang Z-K, Ma Y-H, Zheng J-W, Yang W-D, Liu J-S, Li H-Y. Proteomics to reveal metabolic network shifts towards lipid accumulation following nitrogen deprivation in the diatom Phaeodactylum tricornutum. J Appl Phycol. 2014;26(1):73–82.

    Article  CAS  Google Scholar 

  16. Alipanah L, Rohloff J, Winge P, Bones AM, Brembu T. Whole-cell response to nitrogen deprivation in the diatom Phaeodactylum tricornutum. JExB. 2015;66(20):6281–96.

    CAS  Google Scholar 

  17. Yang Z-K, Niu Y-F, Ma Y-H, Xue J, Zhang M-H, Yang W-D, Liu J-S, Lu S-H, Guan Y, Li H-Y. Molecular and cellular mechanisms of neutral lipid accumulation in diatom following nitrogen deprivation. Biotechnol Biofuels. 2013;6(1):67.

    Article  CAS  Google Scholar 

  18. Huang A, Liu L, Yang C, Wang G. Phaeodactylum tricornutum photorespiration takes part in glycerol metabolism and is important for nitrogen-limited response. Biotechnol Biofuels. 2015;8(1):73.

    Article  Google Scholar 

  19. Huergo LF, Dixon R. The emergence of 2-oxoglutarate as a master regulator metabolite. Microbiol Mol Biol Rev. 2015;79(4):419–35.

    Article  CAS  Google Scholar 

  20. Selim KA, Ermilova E, Forchhammer K. From cyanobacteria to archaeplastida: new evolutionary insights into PII signalling in the plant kingdom. N Phytol. 2020;227(3):722–31.

    Article  Google Scholar 

  21. Jiang Y, Wang X, Sun H, Han S, Li W, Cui N, Lin G, Zhang J, Cheng W, Cao D-D, et al. Coordinating carbon and nitrogen metabolic signaling through the cyanobacterial global repressor NdhR. Proc Natl Acad Sci U S A. 2018;115(2):403–8.

    Article  CAS  Google Scholar 

  22. Huang A, Wu S, Gu W, Li Y, Xie X, Wang G. Provision of carbon skeleton for lipid synthesis from the breakdown of intracellular protein and soluble sugar in Phaeodactylum tricornutum under high CO2. BMC Biotechnol. 2019;19(1):53–53.

    Article  Google Scholar 

  23. Wu S, Huang A, Zhang B, Huan L, Zhao P, Lin A, Wang G. Enzyme activity highlights the importance of the oxidative pentose phosphate pathway in lipid accumulation and growth of Phaeodactylum tricornutum under CO2 concentration. Biotechnol Biofuels. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Sun L, Lin C, Li X, Xing L, Huo D, Sun J, Zhang L, Yang H. Comparative phospho- and acetyl proteomics analysis of posttranslational modifications regulating intestine regeneration in sea cucumbers. Front Physiol. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Gaffney DO, Jennings EQ, Anderson CC, Marentette JO, Shi T, Oxvig A-MS, Streeter MD, Johannsen M, Spiegel DA, Chapman E, et al. Non-enzymatic lysine lactoylation of glycolytic enzymes. Cell Chem Biol. 2020;27(2):206.

    Article  CAS  Google Scholar 

  26. Yang Z, Zheng J, Niu Y, Yang W, Liu J, Li H. Systems-level analysis of the metabolic responses of the diatom Phaeodactylum tricornutum to phosphorus stress. Environ Microbiol. 2014;16(6):1793–807.

    Article  CAS  Google Scholar 

  27. Li W, Gao K, Beardall J. Interactive effects of ocean acidification and nitrogen-limitation on the diatom Phaeodactylum tricornutum. PLoS ONE. 2012;7(12): e51590.

    Article  CAS  Google Scholar 

  28. Sola-Penna M. Metabolic regulation by lactate. IUBMB Life. 2008;60(9):605–8.

    Article  CAS  Google Scholar 

  29. Ballou ER, Avelar GM, Childers DS, Mackie J, Bain JM, Wagener J, Kastora SL, Panea MD, Hardison SE, Walker LA, et al. Lactate signalling regulates fungal beta-glucan masking and immune evasion. Nat Microbiol. 2017.

    Article  Google Scholar 

  30. Liu C, Wu J, Zhu J, Kuei C, Yu J, Shelton J, Sutton SW, Li X, Yun SJ, Mirzadegan T, et al. Lactate inhibits lipolysis in fat cells through activation of an orphan G-protein-coupled receptor, GPR81. J Biol Chem. 2009;284(5):2811–22.

    Article  CAS  Google Scholar 

  31. Zhang D, Tang Z, Huang H, Zhou G, Cui C, Weng Y, Liu W, Kim S, Lee S, Perez-Neut M, et al. Metabolic regulation of gene expression by histone lactylation. Nature. 2019;574(7779):575.

    Article  CAS  Google Scholar 

  32. Adeva-Andany M, Lopez-Ojen M, Funcasta-Calderon R, Ameneiros-Rodriguez E, Donapetry-Garcia C, Vila-Altesor M, Rodriguez-Seijas J. Comprehensive review on lactate metabolism in human health. Mitochondrion. 2014;17:76–100.

    Article  CAS  Google Scholar 

  33. Delwiche CF, Palmer JD. The origin of plastids and their spread via secondary symbiosis. Plant Syst Evol 1997, pp 53–86.

  34. Gibbs S. The chloroplasts of some algal groups may have evolved from endosymbiotic eukaryotic algae. NYASA. 1981;361:193–208.

    CAS  Google Scholar 

  35. McFadden G, van Dooren G. Evolution: red algal genome affirms a common origin of all plastids. Curr Biol. 2004;14(13):514–6.

    Article  Google Scholar 

  36. Medlin LK, Kooistra W, Schmid AMM. A review of the evolution of the diatoms—a total approach using molecules, morphology and geology. The origin and early evolution of the diatoms: fossil, molecular and biogeographical approaches 2000:13–35.

  37. Nisbet R, Kilian O, McFadden G. Diatom genomics: genetic acquisitions and mergers. Curr Biol. 2004;14(24):1048–50.

    Article  Google Scholar 

  38. Scala S, Carels N, Falciatore A, Chiusano ML, Bowler C. Genome properties of the diatom Phaeodactylum tricornutum. Plant Physiol. 2002;129(3):993–1002.

    Article  CAS  Google Scholar 

  39. Zhang H, Yin W, Ma D, Bao Q, Liu X, Xu K, Zhao C, Miao Y, Liu J. Transcriptional reprogramming of intermediate metabolism gene induced by phosphatidylinositol 3-kinase in Phaeodactylum tricornutum. Algal Res Biomass Biofuels Bioprod. 2020.

    Article  Google Scholar 

  40. Ramanan R, Tran Q-G, Cho D-H, Jung J-E, Kim B-H, Shin S-Y, Choi S-H, Liu K-H, Kim D-S, Lee S-J, et al. The ancient phosphatidylinositol 3-kinase signaling system is a master regulator of energy and carbon metabolism in algae. Plant Physiol. 2018;177(3):1050–65.

    Article  CAS  Google Scholar 

  41. Fritz C, Palacios-Rojas N, Feil R, Stitt M. Regulation of secondary metabolism by the carbon-nitrogen status in tobacco: nitrate inhibits large sectors of phenylpropanoid metabolism. Plant J Cell Mol Biol. 2006;46(4):533–48.

    Article  CAS  Google Scholar 

  42. Guillard R. Culture of phytoplankton for feeding marine invertebrates. Culture Mar Invertebr Anim. 1975;26–60.

    Article  Google Scholar 

  43. Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37(8):911–7.

    Article  CAS  Google Scholar 

  44. Guo H, Chen T, Liang Z, Fan L, Shen Y, Zhou D. iTRAQ and PMR-based comparative proteomic profiling in gills of white shrimp Litopenaeus vannamei under copper stress. Chemosphere. 2021.

    Article  PubMed  Google Scholar 

Download references


The authors would like to thank Jingjie PTM BioLab (Hangzhou) Co. Inc. for quantitative proteomic analysis assistance.


The authors are grateful for the financial support from National Natural Science Foundation of China (41876158), Natural Science Foundation of Hainan Province (420QN219), and the Scientific Research Foundation of Hainan University (KYQD(ZR)20060, KYQD(ZR)1967 and KYQD(ZR)21005).

Author information

Authors and Affiliations



AYH carried out the experiments, performed the data analysis and drafted the manuscript. YXL drafted the manuscript. JWD and SYG cultured the P. tricornutum. XNC, XZ, HL and WR analysed data. ZYX conceived of the study, and drafted the manuscript. All the authors read and approved the manuscript.

Corresponding author

Correspondence to Zhenyu Xie.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Methods S1.

Detail parameters for materials and methods.

Additional file 2: Table S1.

Annotation combines of proteomics data.

Additional file 3: Table S2.

Differentially expressed proteins under LN.

Additional file 4: Table S3.

Differentially expressed proteins under HC.

Additional file 5: Table S4.

Annotation combines of lactylated proteomics data.

Additional file 6: Fig. S1.

Overview of proteomics data. (a) Peptide length distribution. (b) Peptide number distribution. (c) Protein coverage distribution. (d) Protein molecular weight distribution.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Huang, A., Li, Y., Duan, J. et al. Metabolomic, proteomic and lactylated proteomic analyses indicate lactate plays important roles in maintaining energy and C:N homeostasis in Phaeodactylum tricornutum. Biotechnol Biofuels 15, 61 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: