Transcriptomic analysis reveals hub genes and pathways in response to acetic acid stress in Kluyveromyces marxianus during high-temperature ethanol fermentation

Abstract

The thermotolerant yeast Kluyveromyces marxianus is known for its potential in high-temperature ethanol fermentation, yet it suffers from excess acetic acid production at elevated temperatures, which hinders ethanol production. To better understand how the yeast responds to acetic acid stress during high-temperature ethanol fermentation, this study investigated its transcriptomic changes under this condition. RNA sequencing (RNA-seq) was used to identify differentially expressed genes (DEGs) and enriched gene ontology (GO) terms and pathways under acetic acid stress. The results showed that 611 genes were differentially expressed, and GO and pathway enrichment analysis revealed that acetic acid stress promoted protein catabolism but repressed protein synthesis during high-temperature fermentation. Protein–protein interaction (PPI) networks were also constructed based on the interactions between proteins coded by the DEGs. Hub genes and key modules in the PPI networks were identified, providing insight into the mechanisms of this yeast's response to acetic acid stress. The findings suggest that the decrease in ethanol production is caused by the imbalance between protein catabolism and protein synthesis. Overall, this study provides valuable insights into the mechanisms of K. marxianus's response to acetic acid stress and highlights the importance of maintaining a proper balance between protein catabolism and protein synthesis for high-temperature ethanol fermentation.

Keywords

Kluyveromyces marxianus / Acetic acid / Transcriptomics / Protein–protein interaction network

References

[1] AnJ, KwonH, KimE, LeeYM, KoHJ, ParkH, ChoiI-G, KimS, KimKH, KimW, ChoiW. Tolerance to acetic acid is improved by mutations of the TATA-binding protein gene. Environ Microbiol, 2015, 17(3):656-669

[2] ArneborgN, JespersenL, JakobsenM. Individual cells of Saccharomyces cerevisiae and Zygosaccharomyces bailii exhibit different short-term intracellular pH responses to acetic acid. Arch Microbiol, 2000, 174(1–2):125-128

[3] BaderGD, HogueCW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinf, 2003, 4: 2

[4] BardJAM, GoodallEA, GreeneER, JonssonE, DongKC, MartinA. Structure and Function of the 26S Proteasome. Annu Rev Biochem, 2018, 87: 697-724

[5] CarrozzaMJ, UtleyRT, WorkmanJL, CôtéJ. The diverse functions of histone acetyltransferase complexes. Trends Genet, 2003, 19(6):321-329

[6] CasalM, CardosoH, LeaoC. Mechanisms regulating the transport of acetic acid in Saccharomyces cerevisiae. Microbiology (reading), 1996, 142(Pt 6):1385-1390

[7] ChinCH, ChenSH, WuHH, HoCW, KoMT, LinCY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol, 2014, 8(Suppl 4):S11

[8] DongY, HuJ, FanL, ChenQ. RNA-Seq-based transcriptomic and metabolomic analysis reveal stress responses and programmed cell death induced by acetic acid in Saccharomyces cerevisiae. Sci Rep, 2017, 7(1):42659

[9] FleckCB, BrockM. Re-characterisation of Saccharomyces cerevisiae Ach1p: Fungal CoA-transferases are involved in acetic acid detoxification. Fungal Genet Biol, 2009, 46(6):473-485

[10] FonsecaGG, HeinzleE, WittmannC, GombertAK. The yeast Kluyveromyces marxianus and its biotechnological potential. Appl Microbiol Biotechnol, 2008, 79(3):339-354

[11] FuX, LiP, ZhangL, LiS. Understanding the stress responses of Kluyveromyces marxianus after an arrest during high-temperature ethanol fermentation based on integration of RNA-Seq and metabolite data. Appl Microbiol Biotechnol, 2019, 103(6):2715-2729

[12] GengP, ZhangL, ShiGY. Omics analysis of acetic acid tolerance in Saccharomyces cerevisiae. World J Microbiol Biotechnol, 2017, 33(5):94

[13] GoldembergJ. Ethanol for a sustainable energy future. Science, 2007, 315(5813):808-810

[14] KimD, LangmeadB, SalzbergSL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods, 2015, 12(4):357-360

[15] Kosaka T, Tsuzuno T, Nishida S, Pattanakittivorakul S, Murata M, Miyakawa I, Lertwattanasakul N, Limtong S, Yamada M (2022) Distinct Metabolic Flow in Response to Temperature in Thermotolerant Kluyveromyces marxianus. Appl Environ Microbiol 88(6):e0200621 doi:https://doi.org/10.1128/aem.02006-21

[16] KurdistaniSK, GrunsteinM. Histone acetylation and deacetylation in yeast. Nat Rev Mol Cell Biol, 2003, 4(4):276-284

[17] LertwattanasakulN, KosakaT, HosoyamaA, SuzukiY, RodrussameeN, MatsutaniM, MurataM, FujimotoN, SuprayogiTK, LimtongS, FujitaN, YamadaM. Genetic basis of the highly efficient yeast Kluyveromyces marxianus: complete genome sequence and transcriptome analyses. Biotechnol Biofuels, 2015, 8: 47

[18] LiB-Z, YuanY-J. Transcriptome shifts in response to furfural and acetic acid in Saccharomyces cerevisiae. Appl Microbiol Biotechnol, 2010, 86(6):1915-1924

[19] LiP, TanX, FuX, DangY, LiS. Metabolomic analysis reveals Kluyveromyces marxianus's stress responses during high-temperature ethanol fermentation. Process Biochem, 2021, 102: 386-392

[20] LimtongS, SringiewC, YongmanitchaiW. Production of fuel ethanol at high temperature from sugar cane juice by a newly isolated Kluyveromyces marxianus. Bioresour Technol, 2007, 98(17):3367-3374

[21] LoveMI, HuberW, AndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014, 15(12):550

[22] MaoX, CaiT, OlyarchukJG, WeiL. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics, 2005, 21(19):3787-3793

[23] MartynovaJ, KokinaA, KibildsJ, LiepinsJ, ScerbakaR, VigantsA. Effects of acetate on Kluyveromyces marxianus DSM 5422 growth and metabolism. Appl Microbiol Biotechnol, 2016, 100(10):4585-4594

[24] McBrideHM, NeuspielM, WasiakS. Mitochondria: More Than Just a Powerhouse. Curr Biol, 2006, 16(14):R551-R560

[25] NonklangS, Abdel-BanatBM, Cha-aimK, MoonjaiN, HoshidaH, LimtongS, YamadaM, AkadaR. High-temperature ethanol fermentation and transformation with linear DNA in the thermotolerant yeast Kluyveromyces marxianus DMKU3-1042. Appl Environ Microbiol, 2008, 74(24):7514-7521

[26] Ohmen JD, Kloeckener-Gruissem B, McEwen JE (1988) Molecular cloning and nucleotide sequence of the nuclear PET122 gene required for expression of the mitochondrial COX3 gene in S.cerevisiae. Nucleic Acids Research 16(22):10783–10802 https://doi.org/10.1093/nar/16.22.10783

[27] Orlandi I, Casatta N, Vai M (2012) Lack of Ach1 CoA-Transferase Triggers Apoptosis and Decreases Chronological Lifespan in Yeast. Frontiers in Oncology 2 https://doi.org/10.3389/fonc.2012.00067

[28] PapaFR, HochstrasserM. The yeast DOA4 gene encodes a deubiquitinating enzyme related to a product of the human tre-2 oncogene. Nature, 1993, 366(6453):313-319

[29] PapaFR, AmerikAY, HochstrasserM. Interaction of the Doa4 deubiquitinating enzyme with the yeast 26S proteasome. Mol Biol Cell, 1999, 10(3):741-756

[30] Petibon C, Malik Ghulam M, Catala M, Abou Elela S (2021) Regulation of ribosomal protein genes: An ordered anarchy. WIREs RNA 12(3):e1632 https://doi.org/10.1002/wrna.1632

[31] PrasaiK, RobinsonLC, TatchellK, HarrisonL. Saccharomyces cerevisiae Mhr1 can bind Xho I-induced mitochondrial DNA double-strand breaks in vivo. Mitochondrion, 2018, 42: 23-32

[32] RosalenyLE, AntúnezO, Ruiz-GarcíaAB, Pérez-OrtínJE, TorderaV. Yeast HAT1 and HAT2 deletions have different life-span and transcriptome phenotypes. FEBS Lett, 2005, 579(19):4063-4068

[33] RugthawornP, MurataY, MachidaM, ApiwatanapiwatW, HirookaA, ThanapaseW, DangjareanH, UshiwakaS, MorimitsuK, KosugiA, AraiT, VaithanomsatP. Growth inhibition of thermotolerant yeast, Kluyveromyces marxianus, in hydrolysates from cassava pulp. Appl Biochem Biotechnol, 2014, 173(5):1197-1208

[34] SalvoA, BritoJ, ArtaxoP, GeigerFM. Reduced ultrafine particle levels in Sao Paulo's atmosphere during shifts from gasoline to ethanol use. Nat Commun, 2017, 8(1):77

[35] ScardoniG, PetterliniM, LaudannaC. Analyzing biological network parameters with CentiScaPe. Bioinformatics, 2009, 25(21):2857-2859

[36] Scully MJ, Norris GA, Falconi TMA, MacIntosh DL (2021) Carbon intensity of corn ethanol in the United States: state of the science. Environmental Research Letters 16(4) doi:https://doi.org/10.1088/1748-9326/abde08

[37] ShannonP, MarkielA, OzierO, BaligaNS, WangJT, RamageD, AminN, SchwikowskiB, IdekerT. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 2003, 13(11):2498-2504

[38] SzklarczykD, GableAL, NastouKC, LyonD, KirschR, PyysaloS, DonchevaNT, LegeayM, FangT, BorkP, JensenLJ, von MeringC. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res, 2021, 49(D1):D605-D612

[39] WangD, WuD, YangX, HongJ. Transcriptomic analysis of thermotolerant yeast Kluyveromyces marxianus in multiple inhibitors tolerance. RSC Adv, 2018, 8: 14177-14192

[40] YoungMD, WakefieldMJ, SmythGK, OshlackA. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol, 2010, 11(2):R14

[41] ZengL, HuangJ, FengP, ZhaoX, SiZ, LongX, ChengQ, YiY. Transcriptomic analysis of formic acid stress response in Saccharomyces cerevisiae. World J Microbiol Biotechnol, 2022, 38(2):34

Funding

National Undergraduate Training Program for Innovation and Entrepreneurship(202198039); Beijing Municipal Education Commission