Jing Chen1,2,3
Yuanhao Chen1,2,3
Weiling He1,2,3Honghao Liang1,2,3Ting Hong1,2,3
Tangcheng Li1,2,3
Hong Du1,2,3*- 1Guangdong Provincial Key Laboratory of Marine Biotechnology, STU-UNIVPM Joint Algal Research Center, Institute of Marine Sciences, Shantou University, Shantou, Guangdong, China
- 2Southern Marine Science and Engineering Guangdong Laboratory, Guangzhou, Guangdong, China
- 3Guangdong Provincial Key Laboratory of Marine Disaster Prediction and Prevention, Shantou University, Shantou, Guangdong, China
Background: The green alga Chlamydomonas reinhardtii can grow photoautotrophically utilizing light and CO2, and heterotrophically utilizing acetate. The physiological and biochemical responses of autotrophy and heterotrophy are different in C. reinhardtii. However, there is no complete understanding of the molecular physiology between autotrophy and heterotrophy. Therefore, we performed biochemical, molecular and transcriptome analysis of C. reinhardtii between autotrophy and heterotrophy.
Results: The cell growth characterization demonstrated that heterotrophic cell had enhanced growth rates, and autotrophic cell accumulated more chlorophyll. The transcriptome data showed that a total of 2,970 differentially expressed genes (DEGs) were identified from photoautotrophy 12h (P12h) to heterotrophy 12h (H12h). The DEGs were involved in photosynthesis, the tricarboxylic acid cycle (TCA), pyruvate and oxidative phosphorylation metabolisms. Moreover, the results of qRT-PCR revealed that the relative expression levels of malate dehydrogenase (MDH), succinate dehydrogenase (SDH), ATP synthase (ATPase), and starch synthase (SSS) were increased significantly from P12h and H12h. The protein activity of NAD-malate dehydrogenase (NAD-MDH) and succinate dehydrogenase (SDH) were significantly higher in the H12h group.
Conclusion: The above results indicated that the high growth rate observed in heterotrophic cell may be the effects of environmental or genetic regulation of photosynthesis. Therefore, the identification of novel candidate genes in heterotrophy will contribute to the development of microalga strains with higher growth capacity and better performance for biomass production.
Introduction
As we all know, microalgae are found in oceans and lakes around the world. They are abundant in vitamins, fatty acids, proteins, and several nutrients, which are important food sources for fish, shellfish, and other animals (Spolaore et al., 2006; Wulff et al., 2009). Microalgae are typically autotrophic organisms that depend on organic and inorganic carbon sources for cell growth (Li et al., 2020). Phototrophic microalgal species can produce energy and carbon for growth using inorganic compounds and light. Heterotrophic microalgal species are wholly dependent on organic material for nutrition (Caron, 2016). In addition, a small proportion of microalgae can be both photoautotrophic and heterotrophic (Markou, 2020). Some microalgal species cannot utilize organic carbon to grow. A deficiency of glucose-concentrating capacity was demonstrated in three cases: Nostoc sp. (PCC7118) and two Synechocystis strains (PCC6308 and PCC7008) (Zhang et al., 1998). Spirulina sp. was able to grow photoautotrophically, heterotrophically, and mixotrophically (Chojnacka and Noworyta, 2004). Chlorella vulgaris was able to grow in both autotrophic and heterotrophic conditions. While autotrophic growth provided higher cellular lipid content, the lipid productivity was much lower compared with heterotrophic growth (Liang et al., 2009). The green alga C. reinhardtii can grow photoautotrophically utilizing CO2, heterotrophically utilizing acetate, and mixotrophically utilizing both carbon sources (Heifetz et al., 2000). C. reinhardtii is an excellent model organism in study photosynthesis, cell physiology, and photoprotective mechanisms (Wakao et al., 2021). C. reinhardtii was photoautotrophic through light-related metabolic pathways, such as photosynthesis, porphyrin, chlorophyll metabolism, and carotenoid biosynthesis. When C. reinhardtii was cultured under heterotrophic conditions, they mainly showed dark hypoxic metabolism and acetic acid metabolism (Boyle and Morgan, 2009a).
With the rapid development of omics, physiological data and omics-based technologies, such as transcriptomics, metabolomics, and proteomics have been increasingly used cooperatively to provide mechanistic insights into growth differences of microalga in different cultivation conditions (Vidotti et al., 2020). Transcriptome sequencing (RNA-seq) is a cost-effective and time-saving approach to survey the putative differentially expressed genes (DEGs) and contribute to understanding biological processes and metabolic pathways of microalgae from autotrophic to heterotrophic conditions (Vidotti et al., 2020; Ma et al., 2023). Transcriptome analysis contributed to finding target genes and pathways that could be the focus of further physiological and molecular studies (Dong et al., 2023). Recent studies have summarized advances in transcriptomic research on heterotrophic and autotrophic protists, including physiology and metabolism, development and life cycles, environmental and ecological studies (Anderson, 2022a).
Previous studies demonstrated that C. reinhardtii mainly fed on cells during photoautotrophy through light-related metabolic pathways, such as photosynthesis, porphyrin, chlorophyll metabolisms, and carotenoid biosynthesis. When C. reinhardtii was cultured under heterotrophic conditions, they mainly showed dark hypoxic metabolism and acetic acid metabolism (Heifetz et al., 2000). However, little information is available on dark oxygen metabolism in C. reinhardtii. Dark culture conditions can accumulate more triacylglycerol (a type of fat) in C. reinhardtii, which have potential applications for biofuel production (Liang, 2013). Therefore, it is critical to investigate mechanisms by which C. reinhardtii can adapt to these transitions, with focus on how the cell alter pathways, such as glycolysis (gluconeogenesis), the TCA cycle, lipid metabolism, energy storage, and cell growth (Liang, 2013; Hu et al., 2018). Additionally, during changes in cellular physiology of C. reinhardtii, several enzymes are notable for playing key roles in the complex network of dark metabolism between photoautotrophic and heterotrophic cultures (Eriksen, 2008; Morales-Sánchez et al., 2015). However, there is limited research that focus on the transcriptome to compare growth differences of C. reinhardtii under autotrophic and heterotrophic conditions. This study attempted to identify molecular alterations by combining physiological, transcriptomic, and molecular methods from autotrophy to heterotrophy. The ultimate goal was to elucidate the molecular basis for growth differences between photoautotrophy and heterotrophy in C. reinhardtii and to provide a theoretical basis for the in-depth research on photoautotrophy and heterotrophy in alga. Additionally, candidate genes for transcriptome screening in heterotrophy will contribute to development of microalga strains with higher growth capacity and better performance for biomass production.
Materials and methods
Cell strain and culture conditions
C. reinhardtii strain CC-125 was obtained from the Chlamydomonas Resource Center (CRC, University of Minnesota). For photoautotrophy, CC-125 cell was continuously cultured in a high-salt (HS) medium (Sueoka, 1960), and heterotrophy was cultured in tris–acetate–phosphate (TAP) medium (Gorman and Levine, 1965) in 500 mL Erlenmeyer flasks (with 150 mL medium). The air was bubbled into the medium at 0.4 L/min using an LZB-3WB float-type flowmeter (Nanjing Shunlaida Measurement and Control Equipment Co., Ltd), which was required for photoautotrophy. Photoautotrophy was cultured under a temperature of 25 °C with a light intensity of 80 μmol m-2 s-1, a 24:0 h (light:dark) photoperiod, and constant shaking at 120 rpm. By contrast, the culture condition of heterotrophy was the same as above except for the light time, which was a 0:24 h (light:dark) photoperiod.
Growth measurement and morphological observation
The cell density of C. reinhardtii was measured by cell counting or optical density methods between autotrophic and heterotrophic cultures. Chlorophyll was extracted using an 80% acetone solution and the concentrations of chlorophyll a (Chl a), chlorophyll b (Chl b), and total chlorophyll were measured using spectrophotometry as described by Lichtenthaler et al (Harmut, 1987). The following equations (Lichtenthaler and Wellburn, 1983) (Unit: μg/mL) were used: Equations 1–3 were used to calculate chlorophyll a chlorophyll b and total chlorophyll content.
RNA quantification and qualification
In the transcriptome analysis, 2×10^7 algal cells cultured in photoautotrophy for 12 hours were collected (for the P12h group). Meanwhile, 2×10^7 algal cells cultured in heterotrophy for 12 hours were collected (for the H12h group). Above cells were transferred into liquid nitrogen for quick freezing and transferred to a -80°C refrigerator for storage. Total RNA was extracted and sequenced by previous methods (Fu et al., 2022). The integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioana-lyzer 2100 system (Agilent Technologies, CA, USA). The transcriptome experiments were performed by Novogene (https://www.novogene.cn/).
Transcriptome analysis
The clustering of index-coded samples was performed on the cBot Cluster Generation System using the Tru Seq PE Cluster Kit v3-cBot-HS (Illumina). After cluster generation, libraries were prepared for sequencing on the Illumina Nova seq platform, generating 150 bp paired-end reads. Clean data (clean reads) were obtained by removing reads containing read adapters, reads containing ploy-N, and low-quality reads (the number of bases in Ophred<= 20 is more than 50% of the entire read length). FPKM was the expected per kilobase of transcript per million fragments mapped. The statistical methods were provided by the DESeq2 R software (1.20.0). The Benjamini and Hochberg methods for reducing the false discovery rate at 5% were used to modify the obtained P-values. Genes identified by DESeq2 as having differential expression were those with an adjusted P-value of ≤ 0.05.
Cluster Profiler R package was used to implement a gene ontology (GO) enrichment analysis of the differentially expressed genes (DEGs) that were adjusted for genes length bias. We used a corrected P-value of ≤ 0.05 and a |log2 (fold change) | ≥ 1 to evaluate significant genes expression differences. GO keywords were deemed to be significantly enriched by DEGs (Young et al., 2010). The KEGG (Kanehisa et al., 2008) database was used to understand high-level biological system functions and utilities based on molecular data, particularly sizable molecular datasets produced by genome sequencing and other high-throughput experimental technologies (http://www.genome.jp/kegg/).
Determination of key enzyme activities
The 2×10^7 cells were collected by centrifugation at 3000 g for 10 min and stored at -80 °C. For enzyme activities, we used 5 million cells in 1 mL PBS and were disrupted by sonication. After centrifugation with 10000 g at 4°C for 10 min, the supernatant was removed and kept on ice for testing. Then, we transferred the supernatant into a test tube containing 10 mL of distilled water and shook before use. Key enzyme activities were then detected using the following kits: Rubisco-Bisphosphate Carboxylase/Oxygenase (Rubisco) Activity Assay Kit; NAD-Malate Dehydrogenase (NAD-MDH) Activity Assay Kit; Soluble Starch Synthase (SSS) Activity Assay Kit; Pyruvate Kinase (PK) Activity Assay Kit; Succinate Dehydrogenase (SDH) Activity Assay Kit; ATP Activity Assay Kit.
RNA extraction and RT-qPCR
The 2×10^7 algal cells were collected and the total RNA was extracted with TRIzol (Takara, Shiga, Japan). The concentration and purity of RNA were assessed using a Nanodrop microspectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). After determining the RNA quality using agarose gel electrophoresis, complementary cDNA synthesis was carried out using Rever Tra Ace qPCR RT Master Mix with gDNA Remover (TOYOBO). Quantitative real-time PCR (qRT-PCR) was performed with an Applied Biosystems 7300 Real Time PCR System and a Roche Light Cycler 96 system, using TB Green® Premix Ex Taq™ II to measure the relative transcript levels of genes associated with metabolic changes. The ubiquitin ligase UBC8 (Phytozomev10.2, cre03.g159200.t1.1) was used as the endogenous control (Díaz-Santos et al., 2016). Cycling conditions were 10 min at 95°C with 40 cycles for melting (30 s at 95°C), annealing (30 s at 60°C), and extension (30 s at 72°C). The genes of interest and the primers were listed in Supplementary Table S6. The data was quantitatively analyzed by 2−△△CT.
Statistical analysis
The experiments were performed with three or four biological replicates and presented as means ± standard deviation (SD). Statistical analysis was performed using the statistical package “SPSS 20.0 for Windows”. The results of qRT-PCR were calculated as described above. PowerPoint, GraphPad Prism 8, and Paint 3D (Windows 10) were used to draw and modify figures. All data were firstly tested for homoscedasticity and normality. The significance of different treatment groups was analyzed with one-way analysis of variance (ANOVA) and two-way ANOVA, both of which were followed by Tukey’s test, with significance determined at P< 0.05 (*) and P< 0.01(**). Bars are means ± SD (n = 3 or 4).
Results
Effects of photoautotrophy and heterotrophy on the growth of C. reinhardtii
As shown in Figure 1, the growth characteristics of C. reinhardtii in autotrophic and heterotrophic conditions were different, including cell numbers, OD750, Fv/Fm variation curves and total chlorophyll content. The cell numbers increased significantly during ten days in photoautotrophy or heterotrophy (Figure 1A). The changes in OD750 of C. reinhardtii were investigated (Figure 1B). It increased significantly to a maximum of 0.66 in photoautotrophy from day 0 to day 4, and then declined slightly over the next few days. By contrast, the heterotrophy increased to its highest point of 0.671 on the eighth day, and then it decreased slowly. The Fv/Fm of photoautotrophy and heterotrophy showed the same downward trend, but heterotrophy was falling faster (Figure 1C). The total chlorophyll content of photoautotrophy rose considerably from day 0 to day 4 and then declined gradually in the following days. By contrast, the total chlorophyll content of photoautotrophy did not change abruptly, with variation from 2.4 μg/mL to 2.6 μg/mL (Figure 1D). To observe the differences in above indicators between P12h to H12h, we showed bar graphs as in Figures 1E–H. The P12h group was significantly higher than the H12h group in the cell numbers, OD750, Fv/Fm, and total chlorophyll content. Photoautotrophy and heterotrophy significantly affected the growth of C. reinhardtii.
Figure 1 The growth characteristics of C. reinhardtii in photoautotrophy and heterotrophy. (A) Number of cells cultured under photoautotrophy and heterotrophy at 0, 2, 4, 6, 8, and 10 days; (B) OD750 variation curves; (C) Fv/Fm variation curves; (D) Total chlorophyll content; (E) Number of cells between P12h and H12h; (F) OD750 from P12h to H12h; (G) Fv/Fm between P12h and H12h; (H) Total chlorophyll content between P12h and H12h. Significance was determined by t-test analysis: ** indicated P < 0.01.
Analysis of transcriptome in P12h and H12h
For a better insight into differences of growth, biochemical composition, and photosynthetic physiology of photoautotrophy and heterotrophy in C. reinhardtii, transcriptome analysis was used in this study. We constructed total RNA libraries of P12h and H12h by RNA-sequencing (Figure 2). According to the Illumina sequencing, the transcriptome sequencing of C. reinhardtii produced 39.84–47.81 million raw reads (Illumina NovaSeq 6000). Fastp (version 0.19.7) was used to filter the raw data, obtaining over 37.57 million clean reads (Supplementary Table S1). The dependability and credibility of data were shown by a strong correlation coefficient (R2 > 0.909) of gene expression between biological replicates (Supplementary Figure S1). To assess differences in the sequencing data between P12h and H12h, holistic principal component analysis (PCA) was performed. The first and second principal components explained 44.61% and 19.38% of variation, respectively (Figure 2A). The Venn diagram showed the number of uniquely expressed genes in P12h or H12h, and overlapping areas showed the number of co-expressed genes in both groups (Figure 2B). The differential genes were screened with p-value ≤ 0.05 and the absolute value of log2 (FC) ≥ 1 as the screening condition. A total of 2970 differential genes (1317 upregulated, 1653 downregulated) were identified between the P12h and the H12h groups (Figure 2C). The heat map indicated expression levels of DEGs in the P12h and the H12h groups (Figure 2D). DEGs clustered with the R package Mclust demonstrated that the group had similar expression profiles.
Figure 2 RNA-seq data expression of C. reinhardtii in photoautotrophy and heterotrophy. (A) PCA analysis, PC1 indicated the first principal component and PC2 indicated the second principal; (B) Venn diagram showing expressed genes in P12h vs. H12h; (C) The volcanic map of DEGs between P12h and H12h; (D) Heatmap analysis of the DEGs in H12h and P12h.
Functional annotations and classification of DEGs
To explore the regulatory role of DEGs between P12h and H12h, GO annotation and KEGG pathway analyses were performed. In the comparison of DEGs from H12h vs. P12h, the GO analysis showed enrichment of three categories, including Biological Process (BP), Cellular Component (CC), and Molecular Function (MF). They were mainly enriched in ATP hydrolysis and energy transport, including sulfuric ester hydrolase activity, cofactor binding, energy-coupled proton transmembrane transport, and coenzyme binding (Supplementary Figure S2). In P12h vs. H12h, the top 20 GO enrichments with DEGs were mainly enriched in membrane, ATP hydrolysis, energy transport cofactor binding, transporter activity and cation binding (Figure 3A). Additionally, KEGG pathway enrichment analysis revealed that the majority of DEGs were involved in secondary metabolism (Supplementary Figure S3). The top 20 KEGG pathway enrichment analysis demonstrated that DEGs in P12h vs. H12h were mainly enriched in several metabolic pathways, including biosynthesis of secondary metabolites, carbon metabolism, biosynthesis of amino acids, oxidative phosphorylation, citrate cycle, chlorophyll metabolism, and others (Figure 3B).
Figure 3 Functional analysis of DEGs. (A) Top 20 GO enrichments with DEGs; (B) Top 20 KEGG pathways enriched with DEGs.
Analyses of the DEGs related to the photosynthesis and carbon fixation metabolic pathways in P12h and H12h
According to the KEGG enrichment analysis, the photosynthesis and carbon fixation metabolic pathways were one of the representative pathways in P12h and H12h. Based on previous studies and transcriptomic data, the simplified photosynthetic carbon fixation metabolic pathway was constructed in Figure 4A. The ten DEGs were identified related to the carbon fixation pathway (Supplementary Table S2), four of these genes were involved in C4-Dicarboxylic acid cycle, including AST, MDN2, MME, and AAT. The expression levels of these genes increased from P12h to H12h according to the FPKM values of transcriptome (Figure 4B). MDN1 and MME were involved in the crassulacean pathway (CAM), and their expression levels were increased in P12h vs. H12h according to the FPKM values. Ribulose-1,5-bisphosphate carboxylase (Rubisco) as a key enzyme in the process of the Calvin cycle, its expression pattern was decreased (Figure 4B).
Figure 4 Differentially expressed genes (DEGs) between P12h and H12h were involved in the photosynthesis and carbon fixation metabolic pathways. (A) Schematic overview of photosynthesis and carbon fixation in chloroplast of C. reinhardtii, PEP: carboxylate phosphoenolpyruvate, PEPC: phosphoenolpyruvate carboxylase, OAA: oxalo acetic acid, Rubisco: ribulose bisphosphate carboxylase, PGA: 3-phosphoglyceric acid, PGAL: 3-phosphoglyceraldehyde, CAM: crassulacean pathway, red words indicated up-regulated DEGs and blue words indicated down-regulated DEGs; (B) The heatmap of relative expression of DEGs.
Analyses of the DEGs related to glycolysis and the TCA cycle metabolic pathways in P12h and H12h
A total of 20 DEGs related to the glycolysis and the TCA cycle metabolic pathways were screened (Supplementary Table S3). A proposed glycolysis and the TCA cycle metabolic pathways were built based on these identified DEGs (Figure 5A). The expression levels of these DEGs were analyzed according to the FPKM values of transcriptome (Figure 5B). The expression levels of two key enzymes in glycolysis pathway phosphoglucomutase (GPM1) and glyceraldehyde 3-phosphate dehydrogenase (GAPN1) were up-regulated. However, the expression levels of three genes (ALD5, PYK6, and GAP2) were down-regulated from P12h to H12h (Figure 5B). Pyruvate, an intermediate of glycolysis, enters the tricarboxylic acid cycle (TCA). The expression levels of 11 key enzymes (PDC3, OGD2, ACH1, SDH2, FUM1, MDN1, SDH3, MDN3, ACLA1, CIS1 and PDC1) were significantly increased between P12h and H12h (Figure 5). The expression of most DEGs significantly increased from P12h to H12h in glycolysis and the TCA cycle metabolic pathways.
Figure 5 Differentially expressed genes (DEGs) between P12h and H12h were involved in glycolysis and the TCA cycle metabolic pathways. (A) Scheme of glycolysis and the TCA cycle in mitochondria of C. reinhardtii. 1,3-DPG, 1,3-diphosphoglyceric acid; Ac-CoA, acetyl-CoA; (B) The heatmap of relative expression of DEGs.
Analyses of the DEGs related to pyruvate metabolic pathway in P12h and H12h
To understand the roles of DEGs in regulating pyruvate metabolic pathway, correlation analysis between pyruvate metabolic network and the expression profiles of the DEGs were performed (Figure 6). In this study, a total of 8 DEGs were identified to be implicated in pyruvate metabolic pathway, including ALD5, PYK6, MME1, PDC1, ACS, MDN1, MDN3 and FUM1 (Figure 6A and Supplementary Table S4). According to the FPKM values of transcriptome, the expression of ALD5 and PYK6 gradually decreased from P12h to H12h. By contrast, the expression levels of MME1, PDC1, ACS, MDN1, MDN3 and FUM1 were dramatically up-regulated (Figure 6B).
Figure 6 Differentially expressed genes (DEGs) between P12h and H12h were involved in pyruvate metabolic pathway. (A) Scheme of pyruvate metabolic network in mitochondria of C. reinhardtii. (B) The heatmap of relative expression of DEGs.
Analyses of the DEGs related to oxidative phosphorylation in P12h and H12h
According to transcriptome results, the 17 DEGs involved in oxidative phosphorylation were significantly up-regulated from P12h to H12h (Figure 7; Supplementary Table S5). Based on previous research we constructed a diagram of complex V (ATP synthase) which was a key enzyme in oxidative phosphorylation (Figure 7A). Moreover, the 12 genes were summarized together with ATPase for catalyze ATP synthesis (Figure 7B). Transcriptome data revealed that expression levels of these genes were significantly increased in P12h vs. H12h by the FPKM values of transcriptome (Figure 7B).
Figure 7 Differentially expressed genes (DEGs) between P12h and H12h were involved in oxidative phosphorylation. (A) Schematic diagram showing ATP synthase. ATP synthase structural model was constructed (Ge et al., 2021). (B) The heatmap of relative expression of DEGs related to oxidative phosphorylation.
The relative expressions of marker genes
To further explore the relationship between DEGs and the major metabolic pathways mentioned above, the relative expression levels of RBCL, MDH, SDH, Atpase, and SSS were detected by qRT-PCR (Figure 8). Detailed analysis indicated that ribulose-1,5-bisphosphate carboxylase (RBCL) transcripts were increased from P12h to H12h (Figure 8A). The relative expression patterns of MDH, SDH, ATPase and SSS were increased significantly between P12h and H12h. As shown in Figure 8B, relative expression levels of NAD-malate dehydrogenase (MDH) in the H12h group were significantly higher than in the P12h group. Figure 8C demonstrated that the expression of succinate dehydrogenase (SDH) in the P12h group was 1.0120 and dramatically increased to 12.2199 in the H12h group. In P12h vs. H12h, the expression of ATP synthase (ATPase) increased significantly (P< 0.05) (Figure 8D). The expression of starch synthase (SSS) increased from only 1.0319 in the P12h group to 13.0615 in the H12h group (Figure 8E). These results were consistent with RNA-seq.
Figure 8 qPCR expression profile of marker genes RBCL, MDH, SDH, Atpase, and SSS between the P12h and the H12h groups. (A) Expression changes of RBCL; (B) Expression changes of MDH; (C) Expression changes of SDH; (D) Expression changes of ATPase; (E) Expression changes of SSS. Data are representative of three independent experiments. Significance was determined by t-test analysis: ** indicated P < 0.01.
The changes of key enzyme activities in P12h and H12h
To comprehend the roles of key DEGs in the regulation of main metabolic pathways, the changes of key enzyme activities were examined. The activity of Rubisco showed a significant decline (p< 0.01) from the P12h group to the H12h group (Figure 9A). In contrast, the NAD-MDH activity of the H12h group was significantly higher than the P12h group (Figure 9B). The changes of SSS activity and pyruvate kinase activity (PK) were no significant among two groups, presented in Figures 9C, 9D. Activity of the SDH was markedly increased (p< 0.01) from the P12h group to the H12h group (Figure 9E).
Figure 9 The key enzyme activities Rubisco, NAD-MDH, SSS, PK, and SDH in the P12h and the H12h groups. (A) Rubisco activity in P12h and H12h. (B) NAD-MDH activity; (C) SSS activity; (D) PK activity; (E) SDH activity. Data are representative of 3 independent experiments. Significance was determined by t-test analysis: * indicated P < 0.05.
Discussion
Photoautotrophy and heterotrophy affect the growth of C. reinhardtii
The photoautotrophic growth is an attractive cultivation because it can be used for natural resources, for example, sunlight and carbon dioxide. Even though photosynthetic autotrophy has these advantages, the maximum biomass concentrations can be negatively correlated with light penetration (Richmond, 2013). The heterotrophy of microalga is cultivated in the dark where the cell growth and reproduction are supported by organic carbons (Chen and Jiang, 2017). It has been reported that heterotrophic growth of the diatom Nitschia laevis, the dinoflagellate Crypthecodinium cohnii, and the rhodophyte Galdieria sulphuraria accumulated much higher cell densities and biomass productivities than autotrophic growth cell (De Swaaf et al., 2003; Wen and Chen, 2003; Schmidt et al., 2005; Graverholt and Eriksen, 2007). The previous results indicated that cultivation conditions had significant effects on C. reinhardtii biomass productivity and bioproduct accumulation (Pessoa et al., 2023). It demonstrated that the biomass, carbohydrates, and chlorophyll b production were higher in photoautotrophic growth. In contrast, lipids and chlorophyll concentrations were higher in heterotrophic growth (Boyle and Morgan, 2009b). Our results represented that Fv/Fm and total chlorophyll content of photoautotrophy were significantly higher than heterotrophy (Figure 1). Our results have the same findings as previous studies. The previous studies reported that the chlorophyll content of Nannochloropsis and Chlorella were reduced in heterotrophic conditions, compared with autotrophic conditions (Cheirsilp and Torpee, 2012).
Transcriptomics reveals the expression of DEGs from autotrophs to heterotrophs
In previous studies, the transcriptome was widely used in heterotrophic and autotrophic microalga, including physiology and metabolism, development and life cycles, environmental and ecological studies (Anderson, 2022b). In this study, the transcriptomics analysis was estimated for samples from photoautotrophy (P12h) and heterotrophy (H12h) (Figure 2). Heatmap showed the FPKM expression patterns that were normalized, mean-centered, and scaled by Z-score method. In the bubble chart of KEGG enrichment, amino acid metabolic pathways and some glycol-metabolism pathways were enriched in H12h vs. P12h (Figure 3). Above results suggested that these metabolic pathways presumably were involved in the growth of C. reinhardtii under heterotrophy.
Functional analysis of DEGs in main metabolic pathways
Further functional evaluation of DEGs was performed by analyzing the important metabolic pathways, including photosynthesis carbon fixation, glycolysis, the TCA cycle, pyruvate metabolism, and oxidative phosphorylation. The Calvin cycle is an important process of carbon fixation in photosynthesis which chloroplasts convert air-borne CO2 into carbohydrates (Yang et al., 2017). In this study, there were 9 DEGs involved in the photosynthesis and carbon fixation metabolic pathways, which were increased in heterotrophy, as shown in Figure 5. Heterotrophic cultures accelerated the photosynthetic carbon fixation process, providing cell with more 3-phosphoglycerate for cellular metabolic processes that may consume excess ATP and NADPH than autotrophic conditions.
Pyruvate is the main product of cytosolic glycolysis and the primary respiratory substrate for energy production to support cell growth in most organisms (Le et al., 2022). It is first oxidized by the mitochondrial pyruvate dehydrogenase complex (PDC) then produces acetyl-CoA which enters the tricarboxylic acid cycle (TCA), finally, generates reductant equivalent for ATP production (Day and Hanson, 1977; Le et al., 2021). The transcriptome data analysis revealed alterations of gene expression of the enzyme coding genes from major metabolic pathways, including Glycolysis, Pyruvate metabolism, and Tricarboxylic acid cycle (TCA) (Figures 5, 6). Several alterations in the expression of enzymes point that TCA may exert an important role in accelerating or increasing the production of citrate and malate in heterotrophy.
Oxidative phosphorylation of cell is the main form of energy generation to maintain vital life activities. There are four main membrane-bound complex I (Hatefi et al., 1979), complex II (Lemire and Oyedotun, 2002), complex III (Tsukihara et al., 1996; Zhang et al., 1998), and complex IV (Poyton et al., 1995) participate in mitochondrial oxidative phosphorylation. In bacteria, ATP synthase consists of eight subunits, α3β3γδϵab2cn, and F1 corresponds to α3β3γδϵ and F0 to ab2cn (Weber and Senior, 2003). Transcriptomic results demonstrated that 12 genes were involved in ATP synthesis and the expression levels were dramatically raised from P12h to H12h in Figure 7. These differences between photoautotrophic and heterotrophic conditions indicated that the cell may transport excess electrons through cyclic electron transport chain under heterotrophic conditions. The genes involved in ATPase were activated in heterotrophic conditions, indicating that ATP was required for cell phosphorylation.
Molecular mechanism through which DEGs regulate between P12h and H12h
In this study, the activity and expression levels of key enzymes were analyzed in the above major metabolic pathways. RBCL is the major soluble protein in plants, which catalyzes the initial step in Calvin’s reductive pentose phosphate cycle (Miziorko and Lorimer, 1983). MDH plays a crucial role in several metabolic pathways, including the TCA cycle, glyoxylate bypass, amino acid synthesis, and so on (Musrati et al., 1998). SDH catalyzes the oxidation of succinate to fumarate and thereby link the TCA cycle (Jardim-Messeder et al., 2015). The transcripts of MDH, SDH, Atpase, and SSS were significantly increased during the H12h treatment, as shown in Figure 8. These differences between autotrophic and heterotrophic conditions indicated that glycolysis and the TCA cycle, oxidative phosphorylation, and pyruvate metabolisms were activated in heterotrophic conditions.
According to previous studies, the autotrophic, mixotrophic, and heterotrophic phenotypes in the microalga Chlorella vulgaris was compared by using proteomics and transcriptomics, indicating that myo-inositol may exert a signaling and regulatory role over the cell growth performance of C. vulgaris in mixotrophy (Vidotti et al., 2020). However, there were limited publications that discussed the potential growth differences of C. reinhardtii under autotrophic and heterotrophic conditions. In present study, we examined pivotal indicators related to the growth of C. reinhardtii under autotrophic and heterotrophic cultures, in which the P12h group had significantly higher cell numbers than the H12h group. Moreover, the transcriptome and molecular analysis indicated that the expression of genes in major metabolisms (glycolysis, tricarboxylic acid cycle, pyruvate metabolism, and oxidative phosphorylation) were activated in heterotrophy. To gain more insight into the regulation of metabolic pathways, we will identify novel candidate genes in heterotrophy and photoautotrophy to explain regulatory mechanisms of cell growth.
Conclusion
Heterotrophic cultures of C. reinhardtii had a significantly higher cell count than autotrophic conditions. Our findings showed enhanced photosynthetic carbon fixation, the TCA cycle, pyruvate and oxidative phosphorylation metabolisms, which ensured an efficient supply of energy and carbon skeleton for rapid growth of C. reinhardtii under heterotrophic conditions. Therefore, the results of this study provided a theoretical basis and an experimental basis for the growth of C. reinhardtii in photoautotrophy and heterotrophy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
Chlamydomonas reinhardtii strain CC-125 was obtained from the Chlamydomonas Resource Center (CRC, University of Minnesota). The study complies with all relevant institutional, national, and international guidelines and legislation in the field.
Author contributions
JC: Writing – original draft, Writing – review & editing. YC: Investigation, Resources, Software, Validation, Visualization, Writing – review & editing. WH: Data curation, Formal analysis, Writing – review & editing. HL: Investigation, Methodology, Writing – original draft. TH: Data curation, Methodology, Writing – review & editing. TL: Supervision, Writing – original draft. HD: Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. We are grateful for the funding of Project supported by Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (SML2023SP204); Science and Technology Plan Projects of Guangdong Province (No. 2021B1212050025); the Program for University Innovation Team of Guangdong Province (2022KCXTD008); the Start- Up funding of Shantou University (NTF23030T).
Conflict of interest
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.
The reviewer HC declared a past co-authorship with the author HD to the handling editor.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1407915/full#supplementary-material
Abbreviations
P12h, photoautotrophy 12h; H12h, heterotrophy 12h; DEGs, Differentially expressed genes; TCA, cycle tricarboxylic acid cycle; RNA-seq, Ribonucleic acid sequencing; PBS, phosphate buffered saline; FPKM, Fragments per kilobase of transcript per million mapped reads; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; NAD-MDH, NAD-Malate Dehydrogenase; SSS, Soluble Starch Synthase; PK, Pyruvate Kinase; PCR, Polymerase chain reaction; qRT-PCR, Quantitative real-time PCR; PCA, Principal component analysis.
References
Anderson, O. R. (2022a). Recent advances in application of transcriptomics: research on heterotrophic and autotrophic protists. Acta Protozoologica 61, 47–75. doi: 10.4467/16890027AP
Anderson, O. R. (2022b). Recent advances in application of transcriptomics: research on heterotrophic and autotrophic protists. Acta Protozoologica 61, 47–75. doi: 10.4467/16890027AP
Boyle, N. R., Morgan, J. A. (2009a). Flux balance analysis of primary metabolism in Chlamydomonas reinhardtii. BMC Syst. Biol. 3, 1–14. doi: 10.1186/1752-0509-3-4
Boyle, N. R., Morgan, J. A. (2009b). Flux balance analysis of primary metabolism in Chlamydomonas reinhardtii. BMC Syst. Biol. 3, 1–14. doi: 10.1186/1752-0509-3-4
Caron, D. A. (2016). Mixotrophy stirs up our understanding of marine food webs. Proc. Natl. Acad. Sci. 113, 2806–2808. doi: 10.1073/pnas.1600718113
Cheirsilp, B., Torpee, S. (2012). Enhanced growth and lipid production of microalgae under mixotrophic culture condition: effect of light intensity, glucose concentration and fed-batch cultivation. Biores. Technol. 110, 510–516. doi: 10.1016/j.biortech.2012.01.125
Chen, H.-H., Jiang, J.-G. (2017). Lipid accumulation mechanisms in auto- and heterotrophic microalgae. J. Agric. Food Chem. 65, 8099–8110. doi: 10.1021/acs.jafc.7b03495
Chojnacka, K., Noworyta, A. (2004). Evaluation of Spirulina sp. growth in photoautotrophic, heterotrophic and mixotrophic cultures. Enzyme Microbial. Technol. 34, 461–465. doi: 10.1016/j.enzmictec.2003.12.002
Day, D. A., Hanson, J. B. (1977). Pyruvate and malate transport and oxidation in corn mitochondria. Plant Physiol. 59, 630–635. doi: 10.1104/pp.59.4.630
De Swaaf, M. E., Sijtsma, L., Pronk, J. T. (2003). High-cell-density fed-batch cultivation of the docosahexaenoic acid producing marine alga Crypthecodinium cohnii. Biotechnol. Bioeng. 81, 666–672. doi: 10.1002/bit.10513
Díaz-Santos, E., Vila, M., Vigara, J., León, R. (2016). A new approach to express transgenes in microalgae and its use to increase the flocculation ability of Chlamydomonas reinhardtii. J. Appl. Phycol. 28, 1611–1621. doi: 10.1007/s10811-015-0706-2
Dong, Q., Wu, Y., Li, B., Chen, X., Peng, L., Sahito, Z. A., et al. (2023). Multiple insights into lignin-mediated cadmium detoxification in rice (Oryza sativa). J. Hazardous Mater. 458, 131931. doi: 10.1016/j.jhazmat.2023.131931
Eriksen, N. T. (2008). The technology of microalgal culturing. Biotechnol. letters 30, 1525–1536. doi: 10.1007/s10529-008-9740-3
Fu, A., Zheng, Y., Lv, Y., Watkins, C. B., Bai, C., Ma, L., et al. (2022). Multi-omics analysis reveals specific modifications associated with reduced chilling injury in bell pepper fruit by methyl jamonate. Postharvest Biol. Technol. 185, 111799. doi: 10.1016/j.postharvbio.2021.111799
Ge, Q., Jia, D., Cen, D., Qi, Y., Shi, C., Li, J., et al. (2021). Micropeptide ASAP encoded by LINC00467 promotes colorectal cancer progression by directly modulating ATP synthase activity. J. Clin. Invest. 131 (22). doi: 10.1172/JCI152911
Gorman, D. S., Levine, R. P. (1965). Cytochrome f and plastocyanin: their sequence in the photosynthetic electron transport chain of Chlamydomonas reinhardi. Proc. Natl. Acad. Sci. 54, 1665–1669. doi: 10.1073/pnas.54.6.1665
Graverholt, O. S., Eriksen, N. T. (2007). Heterotrophic high-cell-density fed-batch and continuous-flow cultures of Galdieria sulphuraria and production of phycocyanin. Appl. Microbiol. Biotechnol. 77, 69–75. doi: 10.1007/s00253-007-1150-2
Harmut, A. (1987). Chlorophylls and carotenoids: pigments of photosynthetic biomembranes. Methods Enzymol. 148, 350–382. doi: 10.1016/0076-6879(87)48036-1
Hatefi, Y., Galante, Y. M., Stiggall, D. L., Ian Ragan, C. (1979). Proteins, polypeptides, prosthetic groups, and enzymic properties of complexes I, II, III, IV, and V of the mitochondrial oxidative phosphorylation system. Methods Enzymol. 56, 577–602. doi: 10.1016/0076-6879(79)56056-X
Heifetz, P. B., Förster, B., Osmond, C. B., Giles, L. J., Boynton, J. E. (2000). Effects of acetate on facultative autotrophy in chlamydomonas reinhardtii assessed by photosynthetic measurements and stable isotope analyses. Plant Physiol. 122, 1439–1446. doi: 10.1104/pp.122.4.1439
Hu, J., Nagarajan, D., Zhang, Q., Chang, J. S., Lee, D. J. (2018). Heterotrophic cultivation of microalgae for pigment production: a review. Biotechnol. Advances 36, 54–67. doi: 10.1016/j.biotechadv.2017.09.009
Jardim-Messeder, D., Caverzan, A., Rauber, R., de Souza Ferreira, E., Margis-Pinheiro, M., Galina, A. (2015). Succinate dehydrogenase (mitochondrial complex II) is a source of reactive oxygen species in plants and regulates development and stress responses. New Phytologist 208, 776–789. doi: 10.1111/nph.13515
Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, D480–D484. doi: 10.1093/nar/gkm882
Le, X. H., Lee, C. P., Millar, A. H. (2021). The mitochondrial pyruvate carrier (MPC) complex mediates one of three pyruvate-supplying pathways that sustain Arabidopsis respiratory metabolism. Plant Cell. 33, 2776–2793. doi: 10.1093/plcell/koab148
Le, X. H., Lee, C. P., Monachello, D., Millar, A. H. (2022). Metabolic evidence for distinct pyruvate pools inside plant mitochondria. Nat. Plants 8, 694–705. doi: 10.1038/s41477-022-01165-3
Lemire, B. D., Oyedotun, K. S. (2002). The Saccharomyces cerevisiae mitochondrial succinate:ubiquinone oxidoreductase. Biochim. Biophys. Acta (BBA)-Bioenergetics 1553, 102–116. doi: 10.1016/S0005-2728(01)00229-8
Li, T., Yang, F., Xu, J., Wu, H., Mo, J., Dai, L., et al. (2020). Evaluating differences in growth, photosynthetic efficiency, and transcriptome of Asterarcys sp. SCS-1881 under autotrophic, mixotrophic, and heterotrophic culturing conditions. Algal. Res. 45, 101753. doi: 10.1016/j.algal.2019.101753
Liang, Y. (2013). Producing liquid transportation fuels from heterotrophic microalgae. Appl. Energy 104, 860–868. doi: 10.1016/j.apenergy.2012.10.067
Liang, Y., Sarkany, N., Cui, Y. (2009). Biomass and lipid productivities of Chlorella vulgaris under autotrophic, heterotrophic and mixotrophic growth conditions. Biotechnol. Letters 31, 1043–1049. doi: 10.1007/s10529-009-9975-7
Lichtenthaler, H. K., Wellburn, A. R. (1983). Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents. Biochem. Soc. Trans. 11, 591–592. doi: 10.1042/bst0110591
Ma, R., Zhang, Z., Fang, H., Liu, X., Ho, S.-H., Xie, Y., et al. (2023). Unveiling the underlying molecular mechanisms of high lutein production efficiency in Chlorella sorokiniana FZU60 under a mixotrophy/photoautotrophy two-stage strategy by transcriptomic, physiological, and biochemical analyses. Biotechnol. Biofuels Bioproducts 16, 1–15. doi: 10.1186/s13068-023-02300-8
Markou, G. (2020). “Overview of microalgal cultivation, biomass processing and application,” in Handbook of Algal Science, technology and medicine (Academic Press), 343–352. doi: 10.1016/B978-0-12-818305-2.00021-8
Miziorko, H. M., Lorimer, G. H. (1983). Ribulose-1,5-bisphosphate carboxylase-oxygenase. Annu. Rev. Biochem. 52, 507–535. doi: 10.1146/annurev.bi.52.070183.002451
Morales-Sánchez, D., Martinez-Rodriguez, O. A., Kyndt, J., Martinez, A. (2015). Heterotrophic growth of microalgae: metabolic aspects. World J. Microbiol. Biotechnol. 31, 1–9. doi: 10.1007/s11274-014-1773-2
Musrati, R. A., Kollarova, M., Mernik, N., Mikulasova, D. (1998). Malate dehydrogenase: distribution, function and properties. Gen. Physiol. biophys. 17, 193–210.
Pessoa, J. S., de Oliveira, C. F. M., Mena-Chalco, J. P., de Carvalho, J. C. M., Ferreira‐Camargo, L. S.. (2023). Trends on Chlamydomonas reinhardtii growth regimes and bioproducts. Biotechnol. Appl. Biochem. 70, 1830–1842. doi: 10.1002/bab.2486
Poyton, R. O., Goehring, B., Droste, M., Sevarino, K. A., Allen, L. A., Zhao, X.-J. (1995). Cytochrome-c oxidase from Saccharomyces cerevisiae. Methods Enzymol. 260, 97–116. doi: 10.1016/0076-6879(95)60133-3
Richmond, A. (2013). “Biological principles of mass cultivation of photoautotrophic microalgae,” in Handbook of microalgal culture: applied phycology and biotechnology, 169–204. doi: 10.1002/9781118567166
Schmidt, R. A., Wiebe, M. G., Eriksen, N. T. (2005). Heterotrophic high cell-density fed-batch cultures of the phycocyanin-producing red alga Galdieria sulphuraria. Biotechnol. Bioeng. 90, 77–84. doi: 10.1002/bit.20417
Spolaore, P., Joannis-Cassan, C., Duran, E., Isambert, A. (2006). Commercial applications of microalgae. J. Biosci. Bioeng. 101, 87–96. doi: 10.1263/jbb.101.87
Sueoka, N. (1960). Mitotic replication of deoxyribonucleic acid in chlamydomonas reinhardi. Proc. Natl. Acad. Sci. 46, 83–91. doi: 10.1073/pnas.46.1.83
Tsukihara, T., Aoyama, H., Yamashita, E., Tomizaki, T., Yamaguchi, H., Shinzawa-Itoh, K., et al. (1996). The Whole Structure of the 13-Subunit Oxidized Cytochrome c Oxidase at 2.8 Å. Science 272, 1136–1144. doi: 10.1126/science.272.5265.1136
Vidotti, A. D. S., Riaño-Pachón, D. M., Mattiello, L., Giraldi, L. A., Winck, F. V., Franco, T. T. (2020). Analysis of autotrophic, mixotrophic and heterotrophic phenotypes in the microalgae Chlorella vulgaris using time-resolved proteomics and transcriptomics approaches. Algal. Res. 51, 102060. doi: 10.1016/j.algal.2020.102060
Wakao, S., Shih, P. M., Guan, K., Schackwitz, W., Ye, J., Patel, D., et al. (2021). Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii. PloS Genet. 17, 1009725. doi: 10.1371/journal.pgen.1009725
Weber, J., Senior, A. E. (2003). ATP synthesis driven by proton transport in F1F0-ATP synthase. FEBS Letters 545, 61–70. doi: 10.1016/S0014-5793(03)00394-6
Wen, Z.-Y., Chen, F. (2003). Heterotrophic production of eicosapentaenoic acid by microalgae. Biotechnol. Advances 21, 273–294. doi: 10.1016/S0734-9750(03)00051-X
Wulff, A., Iken, K., Quartino, M. L., Al-Handal, A., Wiencke, C., Clayton, M. N. (2009). Biodiversity, biogeography and zonation of marine benthic micro and macroalgae in the Arctic and Antarctic. Botanica Marina 52, 491–507. doi: 10.1515/BOT.2009.072
Yang, B., Liu, J., Ma, X., Guo, B., Liu, B., Wu, T., et al. (2017). Genetic engineering of the Calvin cycle toward enhanced photosynthetic CO2 fixation in microalgae. Biotechnol. Biofuels 10, 1–13. doi: 10.1186/s13068-017-0916-8
Young, M. D., Wakefield, M. J., Smyth, G. K., Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, 1–12. doi: 10.1186/gb-2010-11-2-r14
Zhang, Z., Huang, L., Shulmeister, V. M., Chi, Y.-I., Kim, K. K., Hung, L.-W., et al. (1998). Electron transfer by domain movement in cytochrome bc1. Nature 392, 677–684. doi: 10.1038/33612
Keywords: C. reinhardtii, photoautotrophy, heterotrophy, transcriptome, metabolic pathways
Citation: Chen J, Chen Y, He W, Liang H, Hong T, Li T and Du H (2024) Transcriptome analysis reveals the molecular mechanism of differences in growth between photoautotrophy and heterotrophy in Chlamydomonas reinhardtii. Front. Plant Sci. 15:1407915. doi: 10.3389/fpls.2024.1407915
Received: 29 March 2024; Accepted: 20 May 2024;
Published: 19 June 2024.
Edited by:
Yan Lu, Western Michigan University, United StatesCopyright © 2024 Chen, Chen, He, Liang, Hong, Li and Du. 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(s) 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: Hong Du, aGR1QHN0dS5lZHUuY24=