Transcriptional regulatory networks of methanol-independent protein expression in Pichia pastoris under the AOX1 promoter with trans-acting elements engineering

Abstract

To explore the differences in the intracellular transcriptional mechanism in carbon-derepressed and wild-type Pichia pastoris strains fed with three different carbon sources. RNA in carbon-derepressed (Δmig1Δmig2Δnrg1-Mit1; Mut) and wild-type (WT) P. pastoris fed with three different carbon sources (dextrose, glycerol, and methanol) were sequenced. Differentially expressed genes (DEGs) associated with these carbon sources were obtained and clustered into modules using weighted gene co-expression network analysis (WGCNA). Signaling pathway enrichment analysis was performed using KEGG, and protein to protein interaction (PPI) network was also constructed. A total of 2536 DEGs were obtained from three intersections, and some of them were enriched in carbon sources and involved in carbon metabolism, secondary metabolisms, and amino acid biosynthesis. Two modules, MEgreenyellow (involved in protease, oxidative phosphorylation, endoplasmic reticulum protein processing, folate carbon pool, and glycerol phospholipid metabolism pathways) and MEmidnightblue (involved in protease, endocytosis, steroid biosynthesis, and hippo signaling pathways) were significantly correlated with the strain type. Eight hub genes and two sub-networks were obtained from PPI network. Sub-network A enriched in proteasomes pathway while sub-network B enriched in ribosome pathway. The genes involved in carbon metabolism, secondary metabolic, and amino acid biosynthesis pathways changed significantly under different carbon sources. The changes in proteasome and ribosome activities play roles in carbohydrate metabolism in the methanol-free P_AOX1 start-up Mut strain.

Keywords

Pichia pastoris / RNA-seq / Carbon sources / AOX1 promoter

References

[1] Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol

[2] Anders S, Huber W (2012) Differential expression of RNA-Seq data at the gene level-the DESeq package. Genome Biol

[3] Bandettini WP, . MultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study. J Cardiovasc Magn Reson, 2012, 14: 83.

[4] Barabási AL. Scale-free networks: a decade and beyond. Science, 2009, 325(5939): 412-413.

[5] Beilharz TH, Preiss T. Translational profiling: the genome-wide measure of the nascent proteome. Brief Funct Genomic Proteomic, 2004, 3(2): 103-111.

[6] Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser, 1995, 57(1): 289-300.

[7] Brierley RA, Bussineau C, Kosson R, Melton A, Siegel RS. Fermentation development of recombinant Pichia pastoris expressing the heterologous gene: bovine lysozyme. Ann N Y Acad Sci, 1990, 589: 350-362.

[8] Chen HX, Chu J, Zhang SL, . Intracellular expression of Vitreoscilla hemoglobin improves S-adenosylmethionine production in a recombinant Pichia pastoris. Appl Microbiol Biotechnol, 2007, 74(6): 1205-1212.

[9] Florea L, Song L, Salzberg SL. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. F1000 Res, 2013, 2: 188.

[10] Inan M, Meagher MM. Non-repressing carbon sources for alcohol oxidase (AOX1) promoter of Pichia pastoris. J Biosci Bioeng, 2001, 92(6): 585-589.

[11] Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res, 2000, 28(1): 27-30.

[12] Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol, 2013, 14: R36.

[13] Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics, 2008, 9(1): 559.

[14] Liang S, Wang B, Pan L, Ye Y, He M, Han S, . Comprehensive structural annotation of Pichia pastoris transcriptome and the response to various carbon sources using deep paired-end RNA sequencing. BMC Genomics, 2014, 9(4): 511-525.

[15] Lin J, Panigraphy D, Trinh L, Folkman J, Shiloach J. Production process for recombinant human angiostatin in Pichia pastoris. J Ind Microbiol Biotechnol, 2000, 24(1): 31-35.

[16] Lin-Cereghino GP, Godfrey L, de la Cruz BJ, Johnson S, Khuongsathiene S, Tolstorukov I, . Mxr1p, a key regulator of the methanol utilization pathway and peroxisomal genes in Pichia pastoris. Mol Cell Biol, 2006, 26(3): 883-897.

[17] Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol, 2007, 25(1): 117-124.

[18] Polupanov AS, Nazarko VY, Sibirny AA. Gss1 protein of the methylotrophic yeast Pichia pastoris is involved in glucose sensing, pexophagy and catabolite repression. Int J Biochem Cell Biol, 2012, 44(11): 1906-1918.

[19] Portnoy T, Margeot A, Linke R, Atanasova L, Fekete E, Sándor E, . The CRE1 carbon catabolite repressor of the fungus Trichoderma reesei: a master regulator of carbon assimilation. BMC Genomics, 2011, 12: 269.

[20] Prielhofer R, Maurer M, Klein J, Wenger J, Kiziak C, Gasser B, Mattanovich D. Induction without methanol: novel regulated promoters enable high-level expression in Pichia pastoris. Microb Cell Fact, 2013, 12(1): 5.

[21] Prielhofer R, Cartwright SP, Graf AB, Valli M, Bill RM, Mattanovich D, Gasser B. Pichia pastoris regulates its gene-specific response to different carbon sources at the transcriptional, rather than the translational, level. BMC Genomics, 2015, 16: 167.

[22] Russmayer H, Buchetics M, Gruber C, Valli M, Grillitsch K, Modarres G, Guerrasio R, Klavins K, Neubauer S, Drexler H, Steiger M, Troyer C, Al Chalabi A, Krebiehl G, Sonntag D, Zellnig G, Daum G, Graf AB, Altmann F, Koellensperger G, Hann S, Sauer M, Mattanovich D, Gasser B. Systems-level organization of yeast methylotrophic lifestyle. BMC Syst Biol, 2015, 13: 80.

[23] Santt O, Pfirrmann T, Braun B, Juretschke J, Kimmig P, Scheel H, Hofmann K, Thumm M, Wolf DH. The yeast GID complex, a novel ubiquitin ligase (E3) involved in the regulation of carbohydrate metabolism. Mol Biol Cell, 2008, 19(8): 3323-3333.

[24] Shannon P, . Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 2003, 13(11): 2498-2504.

[25] Szklarczyk D, et al (2014) STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res:1003

[26] Tang Y, . CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks. BioSystems, 2015, 127: 67-72.

[27] Tomko RJ Jr, Hochstrasser M. Order of the proteasomal ATPases and eukaryotic proteasome assembly. Cell Biochem Biophys, 2011, 60(1–2): 13-20.

[28] Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, . Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc, 2012, 7: 562-578.

[29] Trapnell CR, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc, 2014, 9(10): 2513.

[30] van den Hazel HB, Kielland-Brandt MC, Winther JR. Biosynthesis and function of yeast vacuolar proteases: review. Yeast, 1996, 12: 1-16.

[31] Vogl T, Glieder A. Regulation of Pichia pastoris promoters and its consequences for protein production. New Biotechnol, 2013, 30(4): 385-404.

[32] Wang XY, Yan X (2017) The metabolic difference and mechanism of Komagataella phaffii in glycerol and methanol culture based on transcriptomic analysis. Acta Microbiol Sin:1–17

[33] Wang Z, Wang Y, Zhang D, Li J, Hua Z, Du G, Chen J. Enhancement of cell viability and alkaline polygalacturonate lyase production by sorbitol co-feeding with methanol in Pichia pastoris fermentation. Biores Technol, 2010, 101(4): 1318-1323.

[34] Wang J, Wang X, Shi L, Qi F, Zhang P, Zhang Y, Zhou X, Song Z, Cai M. Methanol-independent protein expression by AOX1 promoter with trans-acting elements engineering and glucose–glycerol-shift induction in Pichia pastoris. Sci Rep, 2017, 7: 41850.

[35] Wu J, Mao X, Cai T, Luo J, Wei L. KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res, 2006, 34: 720-724.

[36] Xie CEA. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res, 2011, 39: 316-322.

[37] Zhan C, Wang S, Sun Y, Dai X, Liu X, Harvey L, McNeil B, Yang Y, Bai Z. The Pichia pastoris transmembrane protein GT1 is a glycerol transporter and relieves the repression of glycerol on AOX1 expression. FEMS Yeast Res, 2016

[38] Zhang P, Zhang W, Zhou X, Bai P, Cregg JM, Zhang Y. Catabolite repression of Aox in Pichia pastoris is dependent on hexose transporter PpHxt1 and pexophagy. Appl Environ Microbiol, 2010, 76(18): 6108.

[39] Zhang CL, Qin ZJ, Wang GZ, Ji ZB, Wang JM. Transcriptome and RNA-Seq technology. Biotechnol Bull, 2012, 12: 51-56.