Meta-analysis of transcriptomic profiles in Dunaliella tertiolecta reveals molecular pathway responses to different abiotic stresses
Bahman Panahi
A * , Mohammad Farhadian B , Nahid Hosseinzadeh Gharajeh C , Seyyed Abolghasem Mohammadi C and Mohammad Amin Hejazi D
A
B
C
D
Abstract
Microalgae are photosynthetic organisms and a potential source of sustainable metabolite production. However, different stress conditions might affect the production of various metabolites. In this study, a meta-analysis of RNA-seq experiments in Dunaliella tertiolecta was evaluated to compare metabolite biosynthesis pathways in response to abiotic stress conditions such as high light, nitrogen deficiency and high salinity. Results showed downregulation of light reaction, photorespiration, tetrapyrrole and lipid-related pathways occurred under salt stress. Nitrogen deficiency mostly induced the microalgal responses of light reaction and photorespiration metabolism. Phosphoenolpyruvate carboxylase, phosphoglucose isomerase, bisphosphoglycerate mutase and glucose-6-phosphate-1-dehydrogenase (involved in central carbon metabolism) were commonly upregulated under salt, light and nitrogen stresses. Interestingly, the results indicated that the meta-genes (modules of genes strongly correlated) were located in a hub of stress-specific protein–protein interaction (PPI) network. Module enrichment of meta-genes PPI networks highlighted the cross-talk between photosynthesis, fatty acids, starch and sucrose metabolism under multiple stress conditions. Moreover, it was observed that the coordinated expression of the tetrapyrrole intermediated with meta-genes was involved in starch biosynthesis. Our results also showed that the pathways of vitamin B6 metabolism, methane metabolism, ribosome biogenesis and folate biosynthesis responded specifically to different stress factors. Since the results of this study revealed the main pathways underlying the abiotic stress, they might be applied in optimised metabolite production by the microalga Dunaliella in future studies. PRISMA check list was also included in the study.
Keywords: central carbon metabolism, meta-analysis, network analysis, RNA-seq, transcriptome, translation-related genes, tetrapyrrole biosynthesis.
References
Adams C, Godfrey V, Wahlen B, Seefeldt L, Bugbee B (2013) Understanding precision nitrogen stress to optimize the growth and lipid content tradeoff in oleaginous green microalgae. Bioresource Technology 131, 188-194.
| Crossref | Google Scholar | PubMed |
Alkayal F, Albion RL, Tillett RL, Hathwaik LT, Lemos MS, Cushman JC (2010) Expressed sequence tag (EST) profiling in hyper saline shocked Dunaliella salina reveals high expression of protein synthetic apparatus components. Plant Science 179, 437-449.
| Crossref | Google Scholar | PubMed |
Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B (Methodological) 57, 289-300.
| Crossref | Google Scholar |
Bienz M (2006) The PHD finger, a nuclear protein-interaction domain. Trends in Biochemical Sciences 31, 35-40.
| Crossref | Google Scholar | PubMed |
Binte Safie SR, Ng YK, Yao L, Lee YK (2018) Growth bottlenecks of microalga Dunaliella tertiolecta in response to an up-shift in light intensity. European Journal of Phycology 53, 509-519.
| Crossref | Google Scholar |
Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120.
| Crossref | Google Scholar | PubMed |
Bray NL, Pimentel H, Melsted P, Pachter L (2016) Near-optimal probabilistic RNA-seq quantification. Nature Biotechnology 34, 525-527.
| Crossref | Google Scholar | PubMed |
Breuer G, Evers WAC, de Vree JH, Kleinegris DMM, Martens DE, Wijffels RH, Lamers PP (2013) Analysis of fatty acid content and composition in microalgae. Journal of Visualized Experiments 80, e50628.
| Crossref | Google Scholar |
Cihlář J, Füssy Z, Horák A, Oborník M (2016) Evolution of the tetrapyrrole biosynthetic pathway in secondary algae: conservation, redundancy and replacement. PLoS ONE 11, e0166338.
| Crossref | Google Scholar | PubMed |
Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Szcześniak MW, Gaffney DJ, Elo LL, Zhang X, Mortazavi A (2016) A survey of best practices for RNA-seq data analysis. Genome Biology 17(1), 1-19.
| Crossref | Google Scholar |
Courchesne NMD, Parisien A, Wang B, Lan CQ (2009) Enhancement of lipid production using biochemical, genetic and transcription factor engineering approaches. Journal of Biotechnology 141, 31-41.
| Crossref | Google Scholar | PubMed |
Daneshafrooz N, Bagherzadeh Cham M, Majidi M, Panahi B (2022a) Identification of potentially functional modules and diagnostic genes related to amyotrophic lateral sclerosis based on the WGCNA and LASSO algorithms. Scientific Reports 12(1), 20144.
| Crossref | Google Scholar | PubMed |
Daneshafrooz N, Joghataei MT, Mehdizadeh M, Alavi A, Barati M, Panahi B, Teimourian S, Zamani B (2022b) Identification of let-7f and miR-338 as plasma-based biomarkers for sporadic amyotrophic lateral sclerosis using meta-analysis and empirical validation. Scientific Reports 12(1), 1373.
| Crossref | Google Scholar | PubMed |
Del Campo JA, García-González M, Guerrero MG (2007) Outdoor cultivation of microalgae for carotenoid production: current state and perspectives. Applied Microbiology and Biotechnology 74, 1163-1174.
| Crossref | Google Scholar | PubMed |
Dellomonaco C, Fava F, Gonzalez R (2010) The path to next generation biofuels: successes and challenges in the era of synthetic biology. Microbial Cell Factories 9, 3.
| Crossref | Google Scholar | PubMed |
Fang L, Qi S, Xu Z, Wang W, He J, Chen X, Liu J (2017) De novo transcriptomic profiling of Dunaliella salina reveals concordant flows of glycerol metabolic pathways upon reciprocal salinity changes. Algal Research 23, 135-149.
| Crossref | Google Scholar |
Farhadian M, Rafat SA, Panahi B, Mayack C (2021) Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process. Scientific Reports 11(1), 2367.
| Crossref | Google Scholar | PubMed |
Feng D, Chen Z, Xue S, Zhang W (2011) Increased lipid production of the marine oleaginous microalgae Isochrysis zhangjiangensis (Chrysophyta) by nitrogen supplement. Bioresource Technology 102, 6710-6716.
| Crossref | Google Scholar | PubMed |
Ghahramani N, Shodja J, Rafat SA, Panahi B, Hasanpur K (2021) Integrative systems biology analysis elucidates mastitis disease underlying functional modules in dairy cattle. Frontiers in Genetics 12, 712306.
| Crossref | Google Scholar | PubMed |
Grünberger F, Knüppel R, Jüttner M, Fenk M, Borst A, Reichelt R, Hausner W, Soppa J, Ferreira-Cerca S, Grohmann D (2020) Exploring prokaryotic transcription, operon structures, rRNA maturation and modifications using Nanopore-based native RNA sequencing [Preprint]. BioRxiv. doi:10.1101/2019.12.18.880849
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li BO, Lieber M, MacManes MD (2013) De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nature protocols 8(8), 1494-1512.
| Crossref | Google Scholar |
Havaux M, Ksas B, Szewczyk A, Rumeau D, Franck F, Caffarri S, Triantaphylidès C (2009) Vitamin B6 deficient plants display increased sensitivity to high light and photo-oxidative stress. BMC Plant Biology 9, 130.
| Crossref | Google Scholar | PubMed |
Ho S-H, Nakanishi A, Kato Y, Yamasaki H, Chang J-S, Misawa N, Hirose Y, Minagawa J, Hasunuma T, Kondo A (2017) Dynamic metabolic profiling together with transcription analysis reveals salinity-induced starch-to-lipid biosynthesis in alga Chlamydomonas sp. JSC4. Scientific Reports 7, 45471.
| Crossref | Google Scholar | PubMed |
Ishizaki A, Fleming GR (2012) Quantum coherence in photosynthetic light harvesting. Annual Review of Condensed Matter Physics 3(1), 333-361.
| Crossref | Google Scholar |
Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Research 28, 27-30.
| Crossref | Google Scholar | PubMed |
Kim M, Park S, Polle JEW, Jin E (2010) Gene expression profiling of Dunaliella sp. acclimated to different salinities. Phycological Research 58, 17-28.
| Crossref | Google Scholar |
Klok AJ, Martens DE, Wijffels RH, Lamers PP (2013) Simultaneous growth and neutral lipid accumulation in microalgae. Bioresource Technology 134, 233-243.
| Crossref | Google Scholar | PubMed |
Kobayashi K, Masuda T (2016) Transcriptional regulation of tetrapyrrole biosynthesis in Arabidopsis thaliana. Frontiers in Plant Science 7, 1811.
| Crossref | Google Scholar | PubMed |
Kong W, Zhong H, Gong Z, Fang X, Sun T, Deng X, Li Y (2019) Meta-analysis of salt stress transcriptome responses in different rice genotypes at the seedling stage. Plants 8, 64.
| Crossref | Google Scholar |
Liang S, Qi Y, Zhao J, Li Y, Wang R, Shao J, Liu X, An L, Yu F (2017) Mutations in the Arabidopsis AtMRS2-11/AtMGT10/VAR5 gene cause leaf reticulation. Frontiers in Plant Science 8, 2007.
| Crossref | Google Scholar | PubMed |
Lopez D, Casero D, Cokus SJ, Merchant SS, Pellegrini M (2011) Algal functional annotation tool: a web-based analysis suite to functionally interpret large gene lists using integrated annotation and expression data. BMC Bioinformatics 12, 282.
| Crossref | Google Scholar | PubMed |
Moreno JI, Martín R, Castresana C (2005) Arabidopsis SHMT1, a serine hydroxymethyltransferase that functions in the photorespiratory pathway influences resistance to biotic and abiotic stress. The Plant Journal 41, 451-463.
| Crossref | Google Scholar | PubMed |
Nogaj LA, Srivastava A, van Lis R, Beale SI (2005) Cellular levels of glutamyl-tRNA reductase and glutamate-1-semialdehyde aminotransferase do not control chlorophyll synthesis in Chlamydomonas reinhardtii. Plant Physiology 139, 389-396.
| Crossref | Google Scholar | PubMed |
Obata T, Fernie AR, Nunes-Nesi A (2013) The central carbon and energy metabolism of marine diatoms. Metabolites 3(2), 325-346.
| Crossref | Google Scholar | PubMed |
Omidbakhshfard MA, Omranian N, Ahmadi FS, Nikoloski Z, Mueller-Roeber B (2012) Effect of salt stress on genes encoding translation-associated proteins in Arabidopsis thaliana. Plant Signaling & Behavior 7, 1095-1102.
| Crossref | Google Scholar | PubMed |
Panahi B, Hejazi MA (2020) Integrative analysis of gene expression and alternative splicing in microalgae grown under heterotrophic condition. PLoS ONE 15(6), e0234710.
| Crossref | Google Scholar | PubMed |
Panahi B, Hejazi MA (2021) Weighted gene co-expression network analysis of the salt-responsive transcriptomes reveals novel hub genes in green halophytic microalgae Dunaliella salina. Scientific Reports 11(1), 1607.
| Crossref | Google Scholar | PubMed |
Panahi B, Shahi A (2024) Trancriptome data mining in combination with co-expression network analysis identifies the functional modules and critical regulators in Hordeum vulgare L. in response to cold stress. Biochemistry and Biophysics Reports 37, 101620.
| Crossref | Google Scholar | PubMed |
Panahi B, Mohammadi SA, Ebrahimie E (2013) Identification of miRNAs and their potential targets in halophyte plant Thellungiella halophila. BioTechnologia 94(3), 285-290.
| Crossref | Google Scholar |
Panahi B, Frahadian M, Dums JT, Hejazi MA (2019a) Integration of cross species RNA-Seq meta-analysis and machine-learning models identifies the most important salt stress–responsive pathways in microalga Dunaliella. Frontiers in Genetics 10, 752.
| Crossref | Google Scholar | PubMed |
Panahi B, Mohammadi SA, Ruzicka K, Abbasi Holaso H, Zare Mehrjerdi M (2019b) Genome-wide identification and co-expression network analysis of nuclear factor-Y in barley revealed potential functions in salt stress. Physiology and Molecular Biology of Plants 25, 485-495.
| Crossref | Google Scholar | PubMed |
Panahi B, Farhadian M, Hejazi MA (2020) Systems biology approach identifies functional modules and regulatory hubs related to secondary metabolites accumulation after transition from autotrophic to heterotrophic growth condition in microalgae. PLoS ONE 15(2), e0225677.
| Crossref | Google Scholar | PubMed |
Perrine Z, Negi S, Sayre RT (2012) Optimization of photosynthetic light energy utilization by microalgae. Algal Research 1(2), 134-142.
| Crossref | Google Scholar |
Prochnik SE, Umen J, Nedelcu AM, et al. (2010) Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science 329, 223-226.
| Crossref | Google Scholar | PubMed |
Qu D, Show P-L, Miao X (2021) Transcription factor ChbZIP1 from alkaliphilic microalgae Chlorella sp. BLD enhancing alkaline tolerance in transgenic Arabidopsis thaliana. International Journal of Molecular Sciences 22(5), 2387.
| Crossref | Google Scholar |
Raschke M, Boycheva S, Crèvecoeur M, Nunes-Nesi A, Witt S, Fernie AR, Amrhein N, Fitzpatrick TB (2011) Enhanced levels of vitamin B6 increase aerial organ size and positively affect stress tolerance in Arabidopsis. The Plant Journal 66, 414-432.
| Crossref | Google Scholar | PubMed |
Rau A, Marot G, Jaffrézic F (2014) Differential meta-analysis of RNA-seq data from multiple studies. BMC Bioinformatics 15, 91.
| Crossref | Google Scholar | PubMed |
Razeghi J, Pishtab PA, Fathi P, Panahi B, Hejazi MA (2021) The feasibility of microalgae Dunaliella identification based on conserved regions of mitochondrial cytochrome b and cytochrome oxidase genes. Cytology and Genetics 55, 558-565.
| Crossref | Google Scholar |
Rest JS, Wilkins O, Yuan W, Purugganan MD, Gurevitch J (2016) Meta-analysis and meta-regression of transcriptomic responses to water stress in Arabidopsis. The Plant Journal 85, 548-560.
| Crossref | Google Scholar | PubMed |
Rismani-Yazdi H, Haznedaroglu BZ, Bibby K, Peccia J (2011) Transcriptome sequencing and annotation of the microalgae Dunaliella tertiolecta: pathway description and gene discovery for production of next-generation biofuels. BMC Genomics 12(1), 148.
| Crossref | Google Scholar |
Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140.
| Crossref | Google Scholar |
San Pedro A, González-López CV, Acién FG, Molina-Grima E (2013) Marine microalgae selection and culture conditions optimization for biodiesel production. Bioresource Technology 134, 353-361.
| Crossref | Google Scholar | PubMed |
Scarano D, Rao R, Corrado G (2017) In silico identification and annotation of non-coding RNAs by RNA-seq and de novo assembly of the transcriptome of tomato fruits. PLoS ONE 12(2), e0171504.
| Crossref | Google Scholar | PubMed |
Shang C, Pang B, Yu H, Gan S, Li Y (2022) Identification of targets of transcription factor WRINKLED1-Like related to lipid biosynthesis from marine microalga Dunaliella parva. Frontiers in Marine Science 8, 807493.
| Crossref | Google Scholar |
Sharifi S, Pakdel A, Ebrahimi M, Reecy JM, Fazeli Farsani S, Ebrahimie E (2018) Integration of machine learning and meta-analysis identifies the transcriptomic bio-signature of mastitis disease in cattle. PLoS ONE 13(2), e0191227.
| Crossref | Google Scholar | PubMed |
Shi T-Q, Wang L-R, Zhang Z-X, Sun X-M, Huang H (2020) Stresses as first-line tools for enhancing lipid and carotenoid production in microalgae. Frontiers in Bioengineering and Biotechnology 8, 610.
| Crossref | Google Scholar | PubMed |
Shin H, Hong S-J, Kim H, Yoo C, Lee H, Choi H-K, Lee C-G, Cho B-K (2015) Elucidation of the growth delimitation of Dunaliella tertiolecta under nitrogen stress by integrating transcriptome and peptidome analysis. Bioresource Technology 194, 57-66.
| Crossref | Google Scholar | PubMed |
Spolaore P, Joannis-Cassan C, Duran E, Isambert A (2006) Commercial applications of microalgae. Journal of Bioscience and Bioengineering 101, 87-96.
| Crossref | Google Scholar | PubMed |
Szklarczyk D, Franceschini A, Wyder S, et al. (2015) STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Research 43, D447-D452.
| Crossref | Google Scholar | PubMed |
Takahashi S, Okubo R, Kanesaki Y, Zhou B, Takaya K, Watanabe S, Tanaka K, Imamura S (2021) Identification of transcription factors and the regulatory genes involved in triacylglycerol accumulation in the unicellular red alga Cyanidioschyzon merolae. Plants 10(5), 971.
| Crossref | Google Scholar |
Tan KWM, Lin H, Shen H, Lee YK (2016) Nitrogen-induced metabolic changes and molecular determinants of carbon allocation in Dunaliella tertiolecta. Scientific Reports 6, 37235.
| Crossref | Google Scholar | PubMed |
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M (2004) MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. The Plant Journal 37, 914-939.
| Crossref | Google Scholar | PubMed |
Trapnell C, Pachter L, Salzberg SL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105-1111.
| Crossref | Google Scholar | PubMed |
Tseng GC, Ghosh D, Feingold E (2012) Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Research 40(9), 3785-3799.
| Crossref | Google Scholar | PubMed |
Wang P, Liang F-C, Wittmann D, Siegel A, Shan S-O, Grimm B (2018) Chloroplast SRP43 acts as a chaperone for glutamyl-tRNA reductase, the rate-limiting enzyme in tetrapyrrole biosynthesis. Proceedings of the National Academy of Sciences 115, E3588-E3596.
| Crossref | Google Scholar |
Wang C, Wang K, Ning J, Luo Q, Yang Y, Huang D, Li H (2021) Transcription factors from Haematococcus pluvialis involved in the regulation of astaxanthin biosynthesis under high light-sodium acetate stress. Frontiers in Bioengineering and Biotechnology 9, 650178.
| Crossref | Google Scholar |
Wattebled F, Ral J-P, Dauvillee D, Myers AM, James MG, Schlichting R, Giersch C, Ball SG, D’Hulst C (2003) STA11, a Chlamydomonas reinhardtii locus required for normal starch granule biogenesis, encodes disproportionating enzyme. Further evidence for a function of α-1,4 glucanotransferases during starch granule biosynthesis in green algae. Plant Physiology 132, 137-145.
| Crossref | Google Scholar | PubMed |
Yao L, Tan KWM, Tan TW, Lee YK (2017) Exploring the transcriptome of non-model oleaginous microalga Dunaliella tertiolecta through high-throughput sequencing and high performance computing. BMC Bioinformatics 18(1), 122.
| Crossref | Google Scholar | PubMed |
Zhang L, Zhang Z, Zhang X, Yao Y, Wang R, Duan B, Fan S (2019) Comprehensive meta-analysis and co-expression network analysis identify candidate genes for salt stress response in Arabidopsis. Plant Biosystems – An International Journal Dealing with all Aspects of Plant Biology 153, 367-377.
| Crossref | Google Scholar |
Zhang X, Yao Y, Li X, Zhang L, Fan S (2020) Transcriptomic analysis identifies novel genes and pathways for salt stress responses in Suaeda salsa leaves. Scientific Reports 10, 4236.
| Crossref | Google Scholar | PubMed |
Zheng Y, Jiao C, Sun H, et al. (2016) iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Molecular Plant 9, 1667-1670.
| Crossref | Google Scholar |
