Combining manipulation of transcription factors and overexpression of the target genes to enhance lignocellulolytic enzyme production in Penicillium oxalicum
Biotechnology for Biofuels volume 10, Article number: 100 (2017)
Lignocellulolytic enzymes are the main enzymes to saccharify lignocellulose from renewable plant biomass in the bio-based economy. The production of these enzymes is transcriptionally regulated by multiple transcription factors. We previously engineered Penicillium oxalicum for improved cellulase production via manipulation of three genes in the cellulase expression regulatory network. However, the potential of combinational engineering of multiple regulators and their targets at protein abundance and activity levels has not been fully explored.
Here, we verified that a point mutation XlnRA871V in transcription factor XlnR enhanced the expression of lignocellulolytic enzymes, particularly hemicellulases, in P. oxalicum. Then, overexpression of XlnRA871V with a constitutive PDE_02864 promoter was combined with the overexpression of cellulase transcriptional activator ClrB and deletion of carbon catabolite repressor CreA. The resulted strain RE-7 showed 8.9- and 51.5-fold increased production of cellulase and xylanase relative to the starting strain M12, respectively. Further overexpression of two major cellulase genes cbh1-2 and eg1 enabled an additional 13.0% improvement of cellulase production. In addition, XlnRA871V led to decreased production of β-glucosidase and amylase, which could be attributed to the reduced transcription of corresponding enzyme-encoding genes.
The results illustrated that combinational manipulation of the involved transcription factors and their target genes was a viable strategy for efficient production of lignocellulolytic enzymes in filamentous fungi. The striking negative effect of XlnRA871V mutation on amylase production was also highlighted.
Lignocellulosic biomass is the most abundant renewable and sustainable resource on earth and has been regarded as a resource of fermentable sugars for the production of liquid biofuels and biochemicals for a long time . As the main agents of plant decay in the environment, filamentous fungi such as Trichoderma reesei [2, 3], Aspergillus spp. [4, 5], and Neurospora crassa [6, 7] can express a complex mixture of enzymes to synergistically deconstruct the polysaccharides in plant cell walls. Commercial biofuel production remains hindered by the high cost of enzyme production for biomass conversion . Thus, construction of lignocellulolytic enzyme high-producing strains is important to improve the economy of bioconversion of lignocellulosic materials.
Penicillium oxalicum wild-type strain 114-2 has been studied for cellulase production for more than 30 years in China . Three main cellobiohydrolases (CBHs), fifteen endoglucanases (EGs), eleven β-glucosidases (BGLs), and fifty-one hemicellulases were predicted to be encoded in its genome [9, 10]. Owing to these enzymes, the lignocellulolytic enzyme system of P. oxalicum is more diverse than that of the main industrial strain T. reesei . A P. oxalicum mutant JU-A10-T, which has enhanced cellulase expression after multiple rounds of mutagenesis and screening, has been applied for commercial cellulase production for 20 years . Comparative genomic analysis of JU-A10-T and 114-2 and subsequent functional verification indicated that a frameshift mutation in the gene encoding carbon catabolite repressor CreA in JU-A10-T strain was responsible for its cellulase hyper-production . The similar case was reported in T. reesei, where the truncation of CRE1 (ortholog of CreA) also dramatically increased cellulase expression in high-producing strain RUT-C30 . The function of CreA/CRE1 in regulating various biological processes has been described in many fungi [13,14,15]. Besides, ClrB/CLR-2, as a pathway-specific transcription factor, is essential for cellulase expression in both N. crassa and Aspergillus nidulans . In our previous work, a transcription factor gene deletion library was constructed in P. oxalicum and several regulators including ClrB, CreA, XlnR, and AmyR were proved to regulate cellulase expression . A strategy of genetically modifying these regulators through gene overexpression and deletion was used to efficiently improve the production of lignocellulolytic enzymes . However, the potential of this strategy has not been fully explored regarding the manipulation of transcription factors at the activity level.
Hemicellulose, as the second most abundant component of lignocellulosic biomass, is a group of heterogeneous polysaccharides including xylan and mannan. . The structural heterogeneity and complex constituents made them require a complex set of enzymes for efficient degradation. The removal of hemicellulose could make cellulose more accessible to cellulolytic enzymes. Thus, improving the production of hemicellulase is a potential strategy for more efficient deconstruction of lignocellulosic biomass . As the major transcriptional activator of xylan degradation and xylose utilization in filamentous fungi, XlnR (orthologs named XYR1 or XLR-1 in different species) has different roles in the regulation of cellulose degradation . In Aspergillus niger and T. reesei, XlnR/XYR1 was proved to be essential for the expression of both cellulase and hemicellulase genes [20, 21]. However, XLR-1 did not play a significant role in cellulase gene expression in N. crassa . In P. oxalicum, we also found that XlnR was indispensible for hemicellulase induction and to a less extent for cellulase expression . A point mutation in XYR1 (A824V) in T. reesei was responsible for the strong deregulation (i.e., inducer-independent high expression) of both cellulase and xylanase gene expressions . The homologous mutation on XLR-1 in N. crassa also resulted in constitutive xylanase expression [22, 23]. These findings offer a potential target for us to further engineer the regulatory pathway for high production of lignocellulolytic enzymes in P. oxalicum.
In this study, a point mutation A871V in XlnR was verified to increase the expression of lignocellulolytic enzyme genes in P. oxalicum. Furthermore, this mutation was combined with the manipulation of genes encoding CreA and ClrB as well as two main cellulases to evaluate the strategy for improving lignocellulolytic enzyme productions through engineering both regulatory proteins and their targets.
Enhanced expression of lignocellulolytic enzyme genes due to the point mutation XlnRA871V
In our previous study, we identified that overexpression of xlnR resulted in significantly improved xylanase production and a slight increase in cellulase production in P. oxalicum . Considering the above-mentioned mutation of XYR1/XLR-1 that markedly enhanced the production of cellulase or xylanase in T. reesei and N. crassa [22, 23], we examined whether the homologous point mutation XlnRA871V in P. oxalicum had the similar effect. Consequently, an xlnR A871V cassette under the control of its native promoter and terminator was constructed (Fig. 1a) and transformed into P. oxalicum wild-type strain 114-2. Several stable transformants of xlnR A871V (Table 1) were obtained after two rounds of conidial separation and purification. One of these transformants was identified by PCR amplification (Fig. 1a), Southern blot (Additional file 1A), and DNA sequencing. These results indicated that the xlnR A871V cassette was solely integrated at the native xlnR locus in this strain. Neither the promoter and terminator regions nor the coding region of xlnR bores any other mutations.
Given the previously reported enhanced production of xylanase and cellulase in gpdA(p)::xlnR mutant where xlnR was expressed under the control of the A. nidulans gpdA promoter , we compared the enzyme production abilities of xlnR A871V mutant with those in the wild-type strain and gpdA(p)::xlnR mutant. The xlnR A871V mutant displayed nearly identical phenotypes on glucose and wheat bran plates to wild-type strain, while gpdA(p)::xlnR mutant displayed an obvious hydrolysis halo on cellulose (Fig. 2). To examine whether the transcriptional levels of lignocellulolytic genes in P. oxalicum respond to XlnRA871V, the strains were transferred to wheat bran medium containing xylan and cellulose as inducers. The transcription of major xylanase gene xyn10A in the xlnR A871V mutant showed a 91.6-fold increase at 4 h and 24.8-fold increase at 24 h compared to those in the wild-type strain (Fig. 3a, b). Transcriptional levels of major cellulase genes cbh1-2 and eg1 also increased by 1.1 to 4.2-folds (Fig. 3a, b). In contrast, the transcription of major β-glucosidase gene bgl1 [24, 25] showed decrease in both xlnR A871V and gpdA(p)::xlnR than those in wild-type strain (Fig. 3a, b). Neither the overexpression nor the mutation of xlnR affected the transcription of creA and clrB (Fig. 3c, d). The above results demonstrated that the mutation XlnRA871V facilitated the expression of xylanolytic enzymes, and moderately the expression of cellulolytic enzymes, in P. oxalicum.
To investigate whether the enzyme production levels correspond to the change of transcriptional levels in xlnR A871V mutant, various lignocellulolytic enzyme activities were analyzed when the strains were grown in wheat bran medium. The xlnR A871V mutant showed significant increases in FPase (filter paper enzyme), pNPCase (CBH), CMCase (EG), xylanase, pNPXase (β-xylosidase) activities, and extracellular protein concentration compared to wild-type strain after 120-h cultivation (Fig. 4a–c, e–g). Particularly, the production of xylanase increased by 8.9-fold in xlnR A871V (Fig. 4e). On the other hand, the level of pNPGase (BGL) was significantly reduced in the mutant (Fig. 4d). The changes in the enzyme production levels were confirmed by zymography analyses of culture supernatants (Additional file 2A–C), which suggested at least four xylanase components were up-regulated in xlnR A871V. Notably, the xlnR A871V mutant showed greater increases of cellulase and hemicellulase production than the gpdA(p)::xlnR strain, particularly in the later phase of cultivation (Fig. 4a–c, e, f; Additional file 3A-C).
When the culture supernatants of xlnR A871V mutant were analyzed by SDS-PAGE, differential protein patterns were clearly observed compared with those in wild-type strain (Fig. 4h). Five of the obviously differential protein bands between the strains were selected and subsequently analyzed by MS–MS (Table 2). The four proteins with higher abundances in xlnR A871V mutant were identified as xylanases (PDE_04478 and PDE_03573) and β-xylosidases (PDE_00049 and PDE_08037), respectively, while the protein with dramatically decreased abundance was found to be glucoamylase Amy15A (PDE_09417). Further qRT-PCR analysis showed that xlnR A871V mutant had 6.9-, 0.9-, 10.6-, and 20.7-fold increases in the transcription of PDE_00049, PDE_03573, PDE_04478, and PDE_08037, respectively, compared with those in wild-type strain in cellulose medium (Additional file 4). These data implied that genes encoding xylanases and xylosidases are among the efficiently regulated targets of XlnRA871V in P. oxalicum.
Identification of PDE_02864(p) as a strong promoter for manipulation of the cellulase transcriptional regulatory pathway
Previous result showed the dose-controlled regulation of cellulase production by ClrB in P. oxalicum . Considering the significant effect of XlnRA871V on the expression of lignocellulolytic enzyme genes, we supposed overexpression of this mutant using promoters stronger than its native promoter would more efficiently improve the enzyme production levels. Four candidate genes including the expansin gene PDE_02102, the xylanase gene PDE_04478, the 40S ribosomal protein S8 gene PDE_02864, and an unannotated gene PDE_07106 were picked out in consideration of their high transcriptional levels in cellulose medium based on the RNA-Seq results  for further exploration of their promoters.
Expressions of the above four genes in wild-type strain were analyzed in the media with glucose or cellulose as sole carbon source. Transcriptional levels of PDE_02102 and PDE_04478 showed typical induced expression on cellulose compared with those on glucose (Fig. 5a, b). Transcription of PDE_02864 showed high and sustained levels under both glucose and cellulose conditions (Fig. 5c). For PDE_07106, the transcription decreased at 24 and 48 h on cellulose (Fig. 5d). When the promoters of these four genes were used to drive the overexpression of clrB, PDE_02864(p)::clrB strain showed the greatest increase in the ratio of hydrolysis halo to colony diameter in either wild-type or gpdA(p)::clrB background  (Fig. 5e, f). As expected, constitutive expression of clrB was achieved in strains PDE_02864(p)::clrB and PDE_07106(p)::clrB (Additional file 5A). Since the FPase activity in PDE_02864(p)::clrB strain was higher than that in PDE_07106(p)::clrB strain in wheat bran medium (Additional file 5B), the PDE_02864 promoter was used for the overexpression of xlnR A871V in the following strain engineering studies.
Enhancement of lignocellulolytic enzyme production by combining XlnRA871V overexpression with the manipulation of creA and clrB
We have previously identified the synergistic effect on improving cellulolytic and xylanolytic gene expressions by overexpression of clrB and deletion of creA in P. oxalicum wild-type strain 114-2 through two-step genetic manipulations [10, 17]. Here, we first validated the efficiency of this combinatorial manipulation in M12 strain , a uracil auxotrophic mutant derived from wild-type 114-2. M12 showed similar production of cellulase and xylanase to wild-type strain (Additional file 6A–C). Consequently, a ΔcreA-gpdA(p)::clrB cassette with gpdA(p)::clrB cassette flanked by creA upstream and downstream sequences was constructed (Fig. 1b) and transformed into M12. Southern blot (Additional file 1B) and PCR (Fig. 1b) showed one of the transformants contained a single-copy integration of ΔcreA-gpdA(p)::clrB cassette at the creA gene locus, which was named RE-6 (Table 1) and further characterized. As expected, RE-6 mutant showed increased hydrolysis halo on cellulose plate (Fig. 6a) and significantly enhanced transcription of lignocellulolytic enzyme-encoding genes (Fig. 7a, b) as well as the production of corresponding enzymes (Fig. 8a–f) than those in M12. Notably, increased transcription of bgl1 (Fig. 7a, b) and production of β-glucosidase (Fig. 8d) were observed in RE-6 compared with those in M12.
Given the synergistic effect of the overexpression of ClrB and XlnR on cellulase expression , we assumed the overexpression of xlnR A871V in RE-6 might further enhance cellulase expression. Correspondingly, PDE_02864(p)::xlnR A871V cassette (Fig. 1c) was constructed and transformed into RE-6, producing a mutant named RE-7 (Table 1). PCR (Fig. 1c) and Southern blot analyses (Additional file 1C) showed the PDE_02864(p)::xlnR A871V expression cassette was integrated into the genome in RE-7. RE-7 showed larger halo on cellulose plate compared with RE-6 (Fig. 6a). In wheat bran medium, RE-7 showed 4.2-, 4.2-, and 23.8-fold increases in transcriptional levels of cbh1-2, eg1, and xyn10A compared with those in RE-6 at 4 h (Fig. 7a), respectively. The FPase, pNPCase, CMCase, xylanase and pNPXase activities, and extracellular protein concentration in RE-7 mutant increased by 0.6-, 1.0-, 0.7-, 3.0-, 1.6-, and 0.3-folds, respectively, when compared with those in RE-6 mutant at 120 h (Fig. 8a–c, e–g). The increases in cellulase and xylanase were also observed in another independent transformant RE-7-2 (Additional file 7A-B). Again, decreased bgl1 transcription at 24 h (Fig. 7b) and pNPGase activity (Fig. 8d) were observed in RE-7 relative to RE-6, which was also revealed by zymography analyses of the culture supernatants (Additional file 2C). Additionally, SDS-PAGE results showed that much more protein bands were detected in RE-7 compared with those in RE-6 mutant (Fig. 8h). Comparison of individual protein bands between strains suggested higher concentrations of xylanases PDE_03573 and PDE_04478 and decreased concentration of amylase Amy 15A in RE-7 mutant than those in RE-6 (Fig. 8h).
Combinational engineering of cellulolytic transcription factors and main cellulase genes conferred further improved cellulase production
To explore whether engineering of cellulase transcription factors and direct tuning the expression of their target genes could be combined to improve cellulase production, cellulase genes cbh1-2 and eg1 were chosen for overexpression together with the manipulation of creA, clrB, and xlnR. These two genes cbh1-2 and eg1 were tightly regulated by ClrB, CreA, and XlnR (Fig. 7) and their products were identified as the major components in the secretome of P. oxalicum .
In order to construct large DNA fragment containing multi-gene expression cassettes, we used RecET direct cloning technology . The final 14.1-kb fragment obtained contained four independent gene expression cassettes including PDE_02864(p)::xlnR A871V, cbh1-2(p)::cbh1-2, eg1(p)::eg1, and pyrG (as the selection marker), as shown in Fig. 1d. The multi-gene fragment was transformed into RE-6 and generated the RE-8 mutant (Table 1) as verified by PCR (Fig. 1d) and Southern blot analyses (Additional file 1D). RE-8 showed a similar colony on wheat bran but a larger halo on cellulose plate compared with those of parental strain RE-6 (Fig. 6a). As expected, RE-8 displayed 1.8- and 2.2-fold greater cbh1-2 transcriptional levels and 3.5- and 5.9-fold higher eg1 transcriptional levels than those in RE-7 at 4 h and 24 h, respectively (Fig. 7a, b). In accordance with the mRNA levels, RE-8 mutant showed 0.8- and 0.8-fold increased pNPCase and CMCase activities relative to RE-7 (Fig. 8b, c). While the transcriptional level of xyn10A was similar between RE-7 and RE-8 (Fig. 7a), xylanase production after 72 h of cultivation was lower in the latter strain (Fig. 8e). The higher level of CMCase production in RE-8 relative to RE-7 was confirmed by zymography analysis (Additional file 2B). Compared with the original strain M12, the quintuple mutant RE-8 showed 10.2-, 18.1-, 30.6-, 29.7-, and 3.9-fold higher FPase, pNPCase, CMCase, xylanase, and pNPXase activities, respectively (Fig. 8a–c, e–f). Moreover, specific activities of FPase and xylanase (U/mg total protein) in RE-8 at 120 h increased by 1.8- and 6.5-folds, respectively, relative to M12. That is to say, the performance of the enzyme mixture produced by RE-8 might have been enhanced compared with M12 with the same protein loading in lignocelluloses hydrolysis.
XlnRA871V negatively regulated amylase production
Cellulase and amylase are the major components in the secretome of P. oxalicum , and the key activator of amylase expression, AmyR, was identified as a negative regulator of cellulase expression . Given the result of SDS-PAGE that the xlnR A871V mutant had a remarkable decrease in the production of the main amylase Amy15A (Fig. 4h), we compared amylase activity and the transcription level of amy15A between wild-type and xlnR A871V strains (Fig. 9a–c). Up to 98% decrease in amylase activity was observed at 120 h in xlnR A871V relative to wild-type (Fig. 9a). Consistently, zymography analysis showed that two protein bands with amylase activity were hardly detected in the mutant (Fig. 9b). Besides, the transcriptional levels of amy15A and amyR in xlnR A871V mutant were about 37 and 18% lower than those in wild-type strain at 4 h and 40 and 41% lower at 24 h, respectively (Fig. 9c). The similar changes were observed by comparing RE-7 carrying xlnR A871V overexpression with its parent strain RE-6 (Fig. 9d–f). While amylase production was enhanced in RE-6 relative to M12 (Fig. 9d, e), further introduction of the xlnR A871V mutation in RE-7 dramatically reduced amylase production (Fig. 9d, e) as well as the transcription of amy15A and amyR (Fig. 9f) to similar or even lower levels than those in M12.
Suppression of the phenotype of creA deletion by overexpression of XlnRA871V
The cre1 knock-out strain showed smaller colony and reduced hyphal growth in T. reesei . Similar phenomenon was also observed in creA deletion mutant in P. oxalicum . In this study, the colony of RE-6 mutant showed a 50% decrease in diameter than that of M12 strain after 84-h cultivation on glucose plate (Fig. 6b). Interestingly, the growth of RE-7 mutant with the xlnR A871V overexpression showed a recovery on glucose but not on wheat bran (Fig. 6a, b). The recovered growth was also observed in RE-8 (Fig. 6a, b) and other two transformants carrying the XlnRA871V mutation constructed from RE-6 (Additional file 7C). To further investigate this effect of xlnR A871V on growth, conidia of the strains were inoculated in liquid glucose medium and cultivated for 48 h. Similar to the results on glucose plates, the decreased hyphal growth of RE-6 mutant relative to M12 was retrieved with XlnRA871V overexpression (Fig. 6c, d).
In this study, we found that XlnRA871V significantly improved the expression of lignocellulolytic enzyme genes in P. oxalicum. The xlnR A871V mutant with a lower xlnR level (Fig. 3c, d) had a greater improvement in enzyme production than gpdA(p)::xlnR, suggesting the high efficiency of engineering the activity but not the abundance of transcription factors. Considering the enhancement of cellulase production by XlnRA871V is less than that of hemicellulase production, combination of XlnRA871V with other regulators of cellulase expression including CreA and ClrB was shown to be efficient for the coordinative production of lignocellulolytic enzymes.
Introduction of XlnRA871V resulted in dramatically reduced amylase production, which could be attributed to the decreased expression of amy15A and amyR. However, other mechanisms (e.g., changes at protein synthesis and secretion levels) cannot be excluded, as the decrease of amy15A expression was less severe than that of amylase activity (Fig. 9a, c). Since AmyR was proved to negatively regulate cellulase expression in cellulose medium in P. oxalicum , the down-regulation of amyR in xlnR A871V mutant could also contribute to the increased cellulase expression. That is to say, XlnRA871V might enhance the expression of lignocellulolytic enzyme genes through both direct and indirect (mediated by AmyR) manners. In addition, the decreased amylase activity in xlnR A871V should release less glucose from starch than wild-type strain in the wheat bran medium, and thus reduced the carbon catabolite repression on lignocellulolytic enzyme expression. Combination of XlnRA871V mutation and deletion of amyR is expected to be conducted in the future to elucidate the genetic interaction between the two regulators and may further improve the production of lignocellulolytic enzymes.
One trade-off effect of manipulating the regulatory pathway for improving cellulase production is that the performance of the secreted enzymes could be affected due to changes in the composition of enzyme mixture. For example, introduction of XlnRA871V reduced the production of β-glucosidase (Figs. 4d, 8d), which is a rate-limiting enzyme in the degradation of lignocellulosic materials . Here, the introduction of extra copies of cbh1-2 and eg1 genes enabled increased production of corresponding enzymes, suggesting the possibility to strengthen the production of specific enzyme components following the engineering of regulatory pathway. Given the fact that efficient degradation of different lignocellulosic materials requires enzyme mixtures with different compositions , direct adjustment of the expression of genes encoding individual enzymes should be critically considered in the genetic engineering of strains for the production of optimized enzyme cocktails. It should be pointed out that the overexpression of these enzyme-encoding genes could produce an amplified effect when being performed in strains with engineered regulatory pathways (e.g., RE-6).
Efficient genetic engineering tools, e.g., recycling of selection marker and assembly of large DNA fragments, are required for manipulating multiple genes in the same strain. Here, we used the uracil auxotrophic strain M12  for genetic engineering, where both the presence and absence of selection marker pyrG could be easily selected . To ensure pyrG being efficiently excised, the six sequences recognized by β-recombinase [30, 31] were added outside the pyrG cassette in this study (Fig. 1c, d), which allowed further manipulation of other targets using the same marker. In addition, we utilized a RecET protocol in the construction of strain RE-8, which is particularly suitable for fast acquisition of gene cluster containing multiple gene expression cassettes .
In this article, we verified the crucial regulatory function of XlnRA871V on lignocellulolytic enzyme expression in P. oxalicum, and combined the manipulation of three regulators (deletion of creA and overexpression of xlnR A871V and clrB) and two main cellulase genes (overexpressions of cbh1-2 and eg1) in one engineered strain. The pentagenic recombinant strain showed significantly improved cellulolytic and hemicellulolytic enzyme production. Our results signify that redesigning the regulation network by combinational manipulation of the activity and abundance of regulators as well as their target genes is a promising strategy for developing lignocellulolytic enzyme hyper-producing strains. Considering the functions of these regulators are conserved in many ascomycete fungi, the strategy should also be suitable to other fungal species.
Strains and culture media
The wild-type strain (CGMCC5302) and all mutants in this article were cultured on wheat bran extract slants to obtain fresh conidia. Vogel’s salts offered the necessary trace elements for hyphal growth in most media unless otherwise noted. Then, 1 L 50× Vogel’s salt solution was prepared. The solution contained 125 g Na3Citrate·2H2O, 250 g KH2PO4, 100 g NH4NO3, 10 g MgSO4·7H2O, 5 g CaCl2·2H2O, 0.25 mg biotin, 5 mL trace element solution (5% citric acid·H2O, 5% ZnSO4·7H2O, 1% Fe(NH4)2(SO4)2·6H2O, 0.25% CuSO4·5H2O, 0.05% MnSO4·H2O, 0.05% H3BO3, and 0.05% Na2MoO4·2H2O, wt/vol), and 755 mL distilled water. The medium for transformation contained 1 M d-sorbitol, 2% glucose, and 1× Vogel’s salt solution. The solid medium for hyphal morphology observation included 1.5% agar and 1× Vogel’s salt solution supplemented with 2% glucose or 1% cellulose, or complex medium referred to simply as wheat bran medium (see below) in this study. The liquid medium containing 2% glucose was used for hyphal growth to produce pre-cultures. The medium containing 2% glucose or 2% cellulose, or the wheat bran medium (2% corn con residue, 0.6% Avicel, 4.66% wheat bran, 1.0% soybean cake powder, 0.2% (NH4)2SO4, 0.28% NaNO3, 0.1% urea, 0.3% KH2PO4, and 0.05% MgSO4, w/v), was employed for the mRNA extraction and cellulase enzyme activity analysis. All plates were incubated in a 30 °C incubator, and all liquid cultures were grown in 300 mL flasks at 30 °C and 200 rpm.
Fresh spore solution was added into wheat bran plates and incubated at 30 °C for 14–15 h. The mycelia were collected, and solution S1 (1.2 M d-sorbitol, 0.1 M KH2PO4, pH 5.6) with 3‰ lysing enzyme from Trichoderma harzianum (Sigma-Aldrich, USA) was added to degrade the mycelial cell wall at 30 °C for 3 h. The solution was filtered into a 50 mL centrifuge tube, and the suspension was centrifuged at 2500 rpm and 4 °C for 10 min. The supernatant was discarded, and the sediment was resuspended with 5 mL solution S2 (1 M d-sorbitol, 50 mM CaCl2, and 10 mM Tris–HCl at pH 7.5). The suspension was centrifuged at 2500 rpm and 4 °C for 10 min. The supernatant was discarded, and the precipitate was resuspended with 400 μL solution S2 and stored on ice. The transformation system, including 5 μL DNA fragment (concentration ≥ 100 ng/μL), 25 μL solution S3 (25% PEG 6000, 50 mM CaCl2, and 10 mM Tris–HCl at pH 7.5), and 100 μL protoplast suspension, was kept on ice for 20 min. Then, 1 mL solution S3 was added to the system, and the system was allowed to stand for another 5 min at room temperature. Afterward, 2 mL solution S2 was added to terminate the reaction. The solution mixture was added into growth medium (1 M d-sorbitol, 2% glucose, 1× Vogel’s salt solution, 0.7% agarose, hygromycin B or pyrithiamine for transformant selection) at 55 °C for plate pouring. After 4-days cultivation at 30 °C, transformants were obtained.
Construction of xlnR A871V, RE-6, RE-7, and RE-8 mutants
xlnR A871V: The xlnR A871V expression cassette was divided into two parts because of the 871 point mutation. The fragments were, respectively, amplified from the genomic DNA of wild-type strain with primer pairs xmU-F + xmU-R and xmCDS-F + xmCDS-R (Additional file 8). The sequence GCC encoding alanine was replaced by GTC encoding valine in primer xmCDS-F, and GGC was replaced by GAC in primer xmU-R. The hph selection maker cassette was amplified from the pSilent-1 plasmid  with primer pair hphS-F + hphS-R (Additional file 8). The xlnR downstream fragment was amplified with xmD-F + xmD-R (Additional file 8). In accordance with the double-joint PCR method , 25 bp homology sequences were, respectively, designed on xmU-R, xmCDS-R, and xmD-F for a four-fragment fusion. Primer pair xmU-F + xmDR was used to amplify the 8.4 kb xlnR A871V cassette. The obtained overexpression cassette was transformed into the wild-type strain.
RE-6: The gpdA(p)::clrB cassette was constructed as previously described . The primer pair gpdA-F2 + ptrA-R1 was used to amplify the gpdA(p)::clrB cassette in this article. The upstream and downstream fragments of creA were amplified with primer pairs creAU-F +creAU-R and creAD-F +creAD-R, respectively. The three fragments were fused in accordance with the double-joint PCR method, and the primer pair creAU-F + creAD-R was used to amplify the final 9.1 kb cassette. The obtained overexpression cassette was transformed into the M12 strain.
RE-7: The xlnR A871V native promoter was replaced with the PDE_02864 promoter. The PDE_02864 promoter was amplified with primer pair PDE_02864-F + PDE_02864-R. The xlnR A871V open reading frame and 3′ untranslated region were amplified with the primer pair xlnR-F + xlnR-R. The pyrG selection marker cassette from A. nidulans with two six sequences on both sides was amplified with primer pair six-pyrG-F + six-pyrG-R. Three fragments were fused in accordance with the double-joint PCR method. Primer pair PDE_02864-F + six-pyrG-R was used to amplify the 6.7 kb PDE_02864(p)::xlnR A871V cassette. The obtained overexpression cassette was transformed into the RE-6 mutant.
RE-8: PDE_02864(p)::xlnR A871, cbh1-2(p)::cbh1-2, and eg1(p)::eg1 were each amplified with PDE_02864-F2 + xlnR-R, cbh1-F + cbh1-R, and eg1-F + eg1-R primer pairs. The pyrG selection marker cassette from A. nidulans with six sequences on both sides was amplified with six-pyrG-F2 + six-pyrG-R2 primer pair. The RecET direct cloning method was employed to join the four fragments. The 20 μL ligation system was prepared using the following components: 0.13 μL T4 DNA polymerase, 0.2 μL purified bovine serum albumin 100×, 2 μL NE buffer 2, 5 μL linearized p15A vector , and 12.67 μL of all amplified fragments. The system was subjected to the following program: incubation at 25 °C for 30 min, 75 °C for 20 min, and 50 °C for 30 min. The ligation system was spotted on a nitrocellulose membrane (Merck Millipore, Ireland) floating on water for 40 min to remove the metal ions and then stored at −20 °C. E. coli GBdir-gyrA462 competent cells were prepared and transformed in accordance with modified methods as described by Wang et al. . The obtained overexpression cassette was transformed into the RE-6 mutant.
Equal conidial suspension (107/mL, 1 μL) was spotted on each plate for phenotype observations. The plates were kept at 30 °C for 4 to 8 days.
Equal spores (107) of strains were inoculated into 100 mL glucose media, and mycelia were gathered at 48 and 72 h. Mycelia were dried at 80 °C for 3 h to constant weight for biomass measurements.
Strains were pre-cultured in 2% glucose medium at 30 °C and 300 rpm for 20 h. Equal amounts of mycelia were transferred to 1× Vogel’s salt solution and cultured for 2 h. Then, the mycelia were transferred to wheat bran/cellulose media. The mycelia were ground in liquid nitrogen and added to 1 mL TRIzol solution. mRNA extraction and cDNA synthesis were strictly performed in accordance with the manufacturer’s instructions of RNAiso™ reagent (TaKaRa, Japan) and PrimeScript RT Reagent Kit (TaKaRa, Japan). The qRT-PCR reaction mixture contained 10 μL SYBR II, 1 μL cDNA template, 1.6 μL primer pair, and 7.4 μL DEPC-treated water. The reaction program was as follows: 95 °C for 2 min, 40 cycles at 95 °C for 10 s, and 61 °C for 30 s. The fluorescence signal was gathered at the end of each extension step at 80 °C, and a melting curve program with the temperature gradient of 0.1 °C/s (from 65 to 95 °C) was performed. The transcriptional level of actin (PDE_01092) was selected as an internal reference for data normalization. The primers are shown in Additional file 8.
Southern blot analysis
Three primer pairs (listed in Additional file 8) were designed to amplify probes xlnR, clrB, and eg1, respectively. The genomic DNA of wild-type, xlnR A871V, and gpdA(p)::xlnR strains were, respectively, digested by restriction enzyme PstΙ and hybridized with probe xlnR for xlnR A871Vcassette identification. HindΙΙΙ was used to digest the genomic DNA of the RE-6 and RE-7 strains for PDE_02864(p)::xlnR A871V cassette identification. ScaΙ was used to digest the genomic DNA of M12, RE-6, RE-7, and RE-8 for ΔcreA-gpdA(p)::clrB cassette identification. BamHΙ was used to digest the genomic DNA of RE-6 and RE-8 strains for PDE_02864(p)::xlnR A871V-cbh1-2(p)::cbh1-2-eg1(p)::eg1 cassette identification. DIG-High Prime DNA Labeling and Detection starter Kit I (Roche, Switzerland) was used for probe preparation, fragment hybridization, and immunological detection according to the manufacturer’s instructions.
Enzyme activity assays
Substrates, including Whatman No. 1 filter paper (Shanghai, China), 10 mg/mL beechwood xylan (Sigma, USA) solution, 1 mg/mL pNPC (Sigma, USA) solution, 10 mg/mL CMC-Na solution (Sigma, USA), 1 mg/mL pNPG (Sigma, USA) solution, and 10 mg/mL starch (Sigma, USA) solution, were prepared for FPase, xylanase, pNPCase, CMCase, pNPGase, and amylase activity assays, respectively. For the FPase activity assay, the reaction system, including a 0.05 g Whatman No. 1 filter paper, 1.5 mL acetate buffer (pH 4.8), and 0.5 mL diluted culture supernatant, was maintained at 50 °C for 60 min. Then, 3 mL DNS  was added to the reaction system. The reaction system was placed in boiling water for 10 min and 20 mL distilled water was added. Absorbance of the reaction system at 540 nm was determined. Up to 1.5 mL beechwood xylan, CMC-Na, or starch solution was added for xylanase, CMCase, or amylase activity assays, respectively. Reaction at 50 °C for 30 min was offered for xylanase and CMCase activity assays, and reaction at 40 °C for 10 min was offered for amylase activity assay. The concentration of reducing sugars was measured with DNS method. pNPCase, pNPXase, and pNPGase activity assays were as follows: the reaction system including 50 μL pNPC/pNPG/pNPX solution and 100 μL diluted culture supernatant was kept at 50 °C for 30 min, and 150 μL 10% Na2CO3 was added to stop the reaction. Absorbance of the reaction system at 420 nm was determined. One unit of enzyme activity was defined as the amount of enzyme required to release 1 μmol product (glucose/xylose/p-nitrophenyl) from the substrate per minute under the assayed conditions. The concentration of extracellular protein was measured by the Bradford Kit (Sangon, Shanghai, China).
Up to 12.5% gel was prepared for protein separation. The mixture of 15 μL culture supernatant and 3 μL 5× loading buffer was placed in boiling water for 10 min and was loaded into the gel. Coomassie brilliant blue R250 (Sangon, Shanghai, China) was used to color the gel for 30 min. Destainer (glacial acetic acid:ethanol:distilled water, 1:1:8,v/v/v) was used to destain the gel for another 2–4 h.
Protein analysis using MS/MS
Protein bands were cut from SDS-PAGE gels and soaked in 30% acetonitrile/100 mM NH4HCO3 solution for color removal. Then, the supernatant was discarded and protein sample was incubated in 100 mM NH4HCO3 at room temperature for 15 min. The supernatant was discarded and the samples were freeze-dried with liquid nitrogen. The samples were digested overnight in 12.5 ng/μL trypsin in 25 mM NH4HCO3 and then extracted thrice with 60% acetonitrile/0.1% TFA. The extracts were pooled and dried completely, and then 2 μL 20% acetonitrile was added to dissolve the extracts and 1 μL sample solution was spotted to the MALDI target and was naturally dried. Up to 0.5 μL 60% acetonitrile/0.1% TFA solution was added and the sample solution was naturally dried. The 5800 MALDI-TOF/TOF (Applied Biosystems, USA) was used to analyze protein peptides. The Nd:YAG laser with 355 nm wave length of exciting light and the 2.2 kV electrospray voltage was applied. MS/MS queries were performed using Mascot search engine 2.2 (Matrix Science, Ltd.). Proteins were identified by searching the P. oxalicum protein database (http://genome.jgi.doe.gov/Penox1/Penox1.home.html).
Xylan (0.25%), CMC-Na (0.25%), and starch (0.25%) were, respectively, added to the gels for SDS-PAGE to analyze the xylanase, CMCase, and amylase activities of samples. The gels were soaked in 5% Triton X-100 for 1 h to remove SDS. The gels were incubated in sodium acetate buffer (0.2 M, pH 4.8) at 50 °C (40 °C for amylase activity assay) for 30 min for enzymatic reaction. The gel containing xylan or CMC-Na was soaked in the solution containing 0.25% Congo Red with 10% ethanol for 30 min for substrate staining and then the gels were transferred to destaining solution (1 M NaCl) for 20 min. Iodine solution (6% KI and 0.6% I2) was used to detect the starch degradation for amylase assay. For β-glucosidase detection, native-PAGE without SDS addition was performed , and the gel was soaked in solution including 0.1% aesculin solution and 0.5% FeCl3 at 50 °C for 10 min.
Statistical tests were done with one-tailed homoscedastic (equal variance) t test in Microsoft Excel 20 **.
The GenBank accession numbers for the five proteins manipulated in this study are as follows: XlnR, EPS32714; ClrB, EPS31045; CreA, EPS28222; CBH1, EPS32984; EG1, EPS32968.
Filter paper enzyme
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Glass NL, Schmoll M, Cate JH, Coradetti S. Plant cell wall deconstruction by ascomycete fungi. Annu Rev Microbiol. 2013;67:477–98.
Kubicek CP, Mikus M, Schuster A, Schmoll M, Seiboth B. Metabolic engineering strategies for the improvement of cellulase production by Hypocrea jecorina. Biotechnol Biofuels. 2009;2:19.
Peterson R, Nevalainen H. Trichoderma reesei RUT-C30–thirty years of strain improvement. Microbiology. 2012;158:58–68.
van Peij NNME, Gielkens MMC, de Vries RP, Visser J, de Graaff LH. The transcriptional activator XlnR regulates both xylanolytic and endoglucanase gene expression in Aspergillus niger. Appl Environ Microbiol. 1998;64(10):3615–9.
Watanabe J, Tanaka H, Mogi Y, Yamazaki T, Suzuki K, Watanabe T, et al. Loss of Aspergillus oryzae amyR function indirectly affects hemicellulolytic and cellulolytic enzyme production. J Biosci Bioeng. 2011;111(4):408–13.
Tian C, Beeson WT, Iavarone AT, Sun J, Marletta MA, Cate JH, et al. Systems analysis of plant cell wall degradation by the model filamentous fungus Neurospora crassa. Proc Natl Acad Sci USA. 2009;106(52):22157–62.
Coradetti ST, Xiong Y, Glass NL. Analysis of a conserved cellulase transcriptional regulator reveals inducer-independent production of cellulolytic enzymes in Neurospora crassa. Microbiol Open. 2013;2(4):595–609.
Himmel ME, Ding SY, Johnson DK, Adney WS, Nimlos MR, Brady JW, et al. Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science. 2007;315(5813):804–7.
Liu G, Zhang L, Wei X, Zou G, Qin Y, Ma L, et al. Genomic and secretomic analyses reveal unique features of the lignocellulolytic enzyme system of Penicillium decumbens. PLoS One. 2013;8(2):e55185.
Li Z, Yao G, Wu R, Gao L, Kan Q, Liu M, et al. Synergistic and dose-controlled regulation of cellulase gene expression in Penicillium oxalicum. PLoS Genet. 2015;11(9):e1005509.
Liu G, Zhang L, Qin Y, Zou G, Li Z, Yan X, et al. Long-term strain improvements accumulate mutations in regulatory elements responsible for hyper-production of cellulolytic enzymes. Sci Rep. 2013;3:1569.
Ilmen M, Thrane C, Penttila M. The glucose repressor gene cre1 of Trichoderma: isolation and expression of a full-length and a truncated mutant form. Mol Genet Genom. 1996;251(4):451–60.
Strauss J, Mach RL, Zeilinger S, Hartler G, Stoffler G, Wolschek M, et al. Cre1, the carbon catabolite repressor protein from Trichoderma reesei. FEBS Lett. 1995;376(1–2):103–7.
Mach-Aigner AR, Omony J, Jovanovic B, van Boxtel AJ, de Graaff LH. d-Xylose concentration-dependent hydrolase expression profiles and the function of CreA and XlnR in Aspergillus niger. Appl Environ Microbiol. 2012;78(9):3145–55.
Sun J, Glass NL. Identification of the CRE-1 cellulolytic regulon in Neurospora crassa. PLoS One. 2011;6(9):e25654.
Coradetti ST, Craig JP, Xiong Y, Shock T, Tian C, Glass NL. Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi. Proc Natl Acad Sci USA. 2012;109(19):7397–402.
Yao G, Li Z, Gao L, Wu R, Kan Q, Liu G, et al. Redesigning the regulatory pathway to enhance cellulase production in Penicillium oxalicum. Biotechnol Biofuels. 2015;8:71.
Sun J, Tian C, Diamond S, Glass NL. Deciphering transcriptional regulatory mechanisms associated with hemicellulose degradation in Neurospora crassa. Eukaryot Cell. 2012;11(4):482–93.
Klaubauf S, Narang HM, Post H, Zhou M, Brunner K, Mach-Aigner AR, et al. Similar is not the same: differences in the function of the (hemi-)cellulolytic regulator XlnR (Xlr1/Xyr1) in filamentous fungi. Fungal Genet Biol. 2014;72:73–81.
Raulo R, Kokolski M, Archer DB. The roles of the zinc finger transcription factors XlnR, ClrA and ClrB in the breakdown of lignocellulose by Aspergillus niger. AMB Express. 2016;6(1):5.
Mach-Aigner AR, Pucher ME, Steiger MG, Bauer GE, Preis SJ, Mach RL. Transcriptional regulation of xyr1, encoding the main regulator of the xylanolytic and cellulolytic enzyme system in Hypocrea jecorina. Appl Environ Microbiol. 2008;74(21):6554–62.
Derntl C, Gudynaite-Savitch L, Calixte S, White T, Mach RL, Mach-Aigner AR. Mutation of the Xylanase regulator 1 causes a glucose blind hydrolase expressing phenotype in industrially used Trichoderma strains. Biotechnol Biofuels. 2013;6(1):62.
Craig JP, Coradetti ST, Starr TL, Glass NL. Direct target network of the Neurospora crassa plant cell wall deconstruction regulators CLR-1, CLR-2, and XLR-1. MBio. 2015;6(5):e01452-15.
Chen M, Qin Y, Cao Q, Liu G, Li J, Li Z, et al. Promotion of extracellular lignocellulolytic enzymes production by restraining the intracellular beta-glucosidase in Penicillium decumbens. Bioresour Technol. 2013;137:33–40.
Yao GS, Wu RM, Kan QB, Gao LW, Liu M, Yang P, et al. Production of a high-efficiency cellulase complex via beta-glucosidase engineering in Penicillium oxalicum. Biotechnol Biofuels. 2016;9:11.
Qin YQ, Zheng K, Liu GD, Chen M, Qu YB. Improved cellulolytic efficacy in penicillium decumbens via heterologous expression of hypocrea jecorina endoglucanase II. Arch Biol Sci. 2013;65(1):305–14.
Portnoy T, Margeot A, Linke R, Atanasova L, Fekete E, Sandor E, et al. The CRE1 carbon catabolite repressor of the fungus Trichoderma reesei: a master regulator of carbon assimilation. BMC Genom. 2011;12:269.
Liu GD, Qin YQ, Li ZH, Qu YB. Development of highly efficient, low-cost lignocellulolytic enzyme systems in the post-genomic era. Biotechnol Adv. 2013;31(6):962–75.
Ouedraogo JP, Arentshorst M, Nikolaev I, Barends S, Ram AFJ. I-SceI-mediated double-strand DNA breaks stimulate efficient gene targeting in the industrial fungus Trichoderma reesei. Appl Microbiol Biotechnol. 2015;99(23):10083–95.
Szewczyk E, Kasuga T, Fan ZL. A new variant of self-excising beta-recombinase/six cassette for repetitive gene deletion and homokaryon purification in Neurospora crassa. J Microbiol Methods. 2014;100:17–23.
Hartmann T, Dumig M, Jaber BM, Szewczyk E, Olbermann P, Morschhauser J, et al. Validation of a self-excising marker in the human pathogen Aspergillus fumigatus by employing the beta-rec/six site-specific recombination system. Appl Environ Microbiol. 2010;76(18):6313–7.
Wang H, Li Z, Jia R, Hou Y, Yin J, Bian X, et al. RecET direct cloning and Redαβ recombineering of biosynthetic gene clusters, large operons or single genes for heterologous expression. Nat Protoc. 2016;11(7):1175–90.
Yu JH, Hamari Z, Han KH, Seo JA, Reyes-Dominguez Y, Scazzocchio C. Double-joint PCR: a PCR-based molecular tool for gene manipulations in filamentous fungi. Fungal Genet Biol. 2004;41(11):973–81.
Zhang X, Liu N, Yang F, Li J, Wang L, Chen G, et al. In situ demonstration and quantitative analysis of the intrinsic properties of glycoside hydrolases. Electrophoresis. 2012;33(2):280–7.
YBQ, XS, ZHL, and LWG designed the study. LWG, CQX, ML, PY, and LLY performed the experiments. LWG and ZHL analyzed the data and wrote the manuscript. YBQ, XS, and ZHL revised the manuscript. All authors read and approved the final manuscript.
We thank Xuezhi Li, Guodong Liu, Yuqi Qin, Jian Zhao, and Yao hua Zhong from Shandong University for helpful advice on this work. We thank Youming Zhang from Shandong University for advice on RecET direct cloning experiment.
The authors declare that they have no competing interests.
Availability of supporting data
All data supporting the conclusions of this article are included within the manuscript and in the additional information.
Consent for publication
All authors agreed to publish this article.
Ethics approval and consent to participate
This study was funded by the National Natural Science Foundation of China (Nos. 30970052 and 31270089 and 31670079), and supported by State Key Laboratory of Microbial Technology Open Projects Fund (No. M2016-07).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Fig. S1 Southern bolt analysis of the constructed strains. (A) Probe xlnR was used to detect the copy numbers of xlnR gene in the wild-type, xlnR A871V, and gpdA(p)::xlnR strains. (B) Probe clrB was used to detect the copy numbers of clrB in M12, RE-6, RE-7, and RE-8 strains. (C) Probe xlnR was used to detect the copy numbers of xlnR in RE-6 and RE-7 strains. (D) Probe eg1 was used to detect the copy numbers of eg1 in RE-6 and RE-8 strains.
Additional file 2. Fig. S2 Activity analysis of lignocellulolytic enzymes by zymography. (A-C) Xylanase, CMCase and β-glucosidase activity analyses of the constructed strains, respectively.
Additional file 3. Fig. S3 Enzyme activity analysis of wild-type, xlnR A871V and gpdA(p)::xlnR strains at 24 h and 48 h. FPase (A), CMCase (B) and xylanase activities (C) of supernatants from wild-type, xlnR A871V and gpdA(p)::xlnR strains on wheat bran were determined. Error bars represent the standard deviations. Statistical significance of the differences between wild-type and each mutant were calculated. *P < 0.05, **P < 0.01, ***P < 0.001.
Additional file 4. Fig. S4 Transcription analysis of xylanase and xylosidase genes in wild-type and xlnR A871V strains. Error bars represent the standard deviations. Statistical significance of the difference between wild-type and each mutant were calculated. *P < 0.05, **P < 0.01, ***P < 0.001.
Additional file 5. Fig. S5 Transcription of clrB and cellulase activity analyses of clrB overexpression mutants. (A) Transcriptional levels of clrB in wild-type and four clrB overexpression strains on glucose and cellulose. (B) FPase activity analysis of gpdA(p)::clrB, PDE_02864(p)::clrB and PDE_07106(p)::clrB strains on wheat bran. Error bars represent the standard deviations. Statistical significance of the difference between gpdA(p)::clrB and each of the other two mutants were calculated. *P < 0.05, **P < 0.01, ***P < 0.001.
Additional file 6. Fig. S6 Lignocellulolytic enzyme activity analysis of wild-type and M12 strains. FPase (A), pNPCase (B) and xylanase activities (C) on wheat bran were determined. Error bars represent the standard deviations.
Additional file 7. Fig. S7 Enzyme activity analysis and phenotype observation of the transformant RE-7-2. FPase (A) and xylanase activities (B) of RE-6 and RE-7-2 strains on wheat bran at 72 h, 96 h and 120 h were determined. Error bars represent the standard deviations. (C) Phenotypic analysis of RE-6 and three parallel transformants (RE-7, RE-7-1 and RE-7-2) on glucose plate after 4 days’ cultivation.
About this article
Cite this article
Gao, L., Li, Z., Xia, C. et al. Combining manipulation of transcription factors and overexpression of the target genes to enhance lignocellulolytic enzyme production in Penicillium oxalicum . Biotechnol Biofuels 10, 100 (2017). https://doi.org/10.1186/s13068-017-0783-3