Weighted Gene Co-expression Network Analysis Identifies Critical Genes for the Production of Cellulase and Xylanase in Penicillium oxalicum

Genes involved in cellular processes undergo environment-dependent co-regulation, but the co-expression patterns of fungal cellulase and xylanase-encoding genes remain unclear. Here, we identified two novel carbon sources, methylcellulose and 2-hydroxyethyl cellulose, which efficiently induced the secretion of cellulases and xylanases in Penicillium oxalicum. Comparative transcriptomic analyses identified carbon source-specific transcriptional patterns, mainly including major cellulase and xylanase-encoding genes, genes involved in glycolysis/gluconeogenesis and the tricarboxylic acid cycle, and genes encoding transcription factors, transporters and G protein-coupled receptors. Moreover, the weighted correlation network analysis of time-course transcriptomes, generated 17 highly connected modules. Module MEivory, comprising 120 members, included major cellulase and xylanase-encoding genes, genes encoding the key regulators PoxClrB and PoxXlnR, and a cellodextrin transporter POX06051/CdtC, which were tightly correlated with the filter-paper cellulase, carboxymethylcellulase and xylanase activities in P. oxalicum. An expression kinetic analysis indicated that members in MEivory were activated integrally by carbon sources, but their expressional levels were carbon source- and/or induction duration-dependent. Three uncharacterized regulatory genes in MEivory were identified, which regulate the production of cellulases and xylanases in P. oxalicum. These findings provide insights into the mechanisms associated with the synthesis and secretion of fungal cellulases and xylanases, and a guide for P. oxalicum application in biotechnology.