- Open Access
Proteomic analysis reflects an environmental alkalinization-coupled pH-dependent mechanism of regulating lignocellulases in Trichoderma guizhouense NJAU4742
Biotechnology for Biofuels volume 13, Article number: 6 (2020)
Filamentous fungi have the ability to efficiently decompose plant biomass, and thus are widely used in the biofuel and bioprocess industries. In process, ambient pH has been reported to strongly affect the performance of the applied functional filamentous fungi. In this study, Trichoderma guizhouense NJAU4742 was investigated under the fermentation of rice straw at different initial pH values for a detailed study.
The results showed that NJAU4742 strain could tolerate ambient pH values ranging from 3.0 to 9.0, but had significantly higher growth speed and extracellular enzyme activities under acidic conditions. At low ambient pH (< 4), NJAU4742 strain achieved rapid degradation of rice straw by elevating the ambient pH to an optimal range through environmental alkalinization. Further proteomic analysis identified a total of 1139 intracellular and extracellular proteins during the solid-state fermentation processes, including the quantified 190 carbohydrate-active enzymes (CAZymes) responsible for rice straw degradation, such as 19 cellulases, 47 hemicellulases and 11 chitinases. Meanwhile, the analysis results clearly showed that the secreted lignocellulases had a synergistic trend in distribution according to the ambient pH, and thus led to a pH-dependent classification of lignocellulases in T. guizhouense NJAU4742.
Most functional lignocellulases were found to be differently regulated by the ambient pH in T. guizhouense NJAU4742, which had the ability of speeding up biomass degradation by elevating the ambient pH through environmental alkalinization. These findings contribute to the theoretical basis for the biodegradation of plant biomass by filamentous fungi in the biofuel and bioprocess industries.
Plant biomass, as the most abundant natural material on earth, mainly consists of cellulose (40–50%), hemicellulose (25–30%) and lignin (15–20%) . For a long time, scientists have tried their best to overcome the structural recalcitrance and time-consuming degradation process of plant biomass [2,3,4]. Detailed in structure, cellulose and hemicellulose are the most abundant natural polysaccharides, where the former exists mostly as a crystalline structure [5, 6], and the latter exists as a mixture of xylan, glucomannan, arabinogalactan, xyloglucan and galactoglucomannan . Hemicellulose is usually decorated with various side-chain substituents, such as acetyl, α-arabinosyl, α-glucuronyl and 4-O-methylglucuronycyl residues . As a complex phenolic polymer, lignin is composed of three monolignol monomers that can be methoxylated to various degrees: p-coumaryl alcohol, coniferyl alcohol, and sinapyl alcohol. These lignols are incorporated into lignin in the form of the phenylpropanoid p-hydroxyphenyl (H), guaiacyl (G), and syringyl (S) units, respectively . Finally, by physical and chemical forces, all three different polymers form a 3D microstructure in the plant cell wall, which is very stable and shows strong resistance to hydration .
So far, plant biomass is widely used in two areas: biofuel and bioorganic fertilizer , in which filamentous fungi are employed as the plant biomass-degraders for their highly efficient secretion of many synergistic enzymes, such as cellulases, hemicellulases, pectinases and ligninases. However, it is still a challenge to completely degrade plant biomass in an economical amount of time just with the secreted enzymes from filamentous fungi. In the biofuel industry, acid or alkaline pretreatment is used and will partly destroy and significantly loosen the entire stable structure of plant biomass, thus providing a less recalcitrant and much more accessible material for the hydrolytic enzymes . Similarly, in the field of bioorganic fertilizer, amino acids with animal origins are added into plant biomass for highly efficient fermentation by functional filamentous fungi, aiming to supply nutrients, slightly pretreat and prevent the growth of other microbes . Certainly, all of these treatments would significantly change the subsequent environmental pH of the next step regardless of biofuel or bioorganic fertilizer production. Additionally, biological organic fertilizer containing the functional filamentous fungi would be applied to different farmlands, which have fluctuant natural pH values ranging from 4.0 to 9.0 . These facts drew our attention to investigate the effect of pH on the growth of filamentous fungi, or more specifically, focus on the effect of pH on the regulation of lignocellulases in functional filamentous fungi.
pH regulation of lignocellulases in filamentous fungi has already been investigated previously. For example, in Aspergillus nidulans, XlnB (xylanase) and AbfB (α-l-arabinofuranosidase) were acid-expressed, while XlnA (xylanase) was alkaline-expressed . In Botrytis cinerea, protein gel identification showed that three hydrolases (α-1,2-mannosidase, α-l-arabinofuranosidase and β-1,3-glucosidase) were up or specifically accumulated at pH 4.0, while four other hydrolases (exocellulase, arabinogalactan endo-1,4-β-galactosidase, rhamnogalacturonan acetyl esterase and exoarabinanase) were up or specifically accumulated at pH 6.0 . Transcriptome analysis of Trichoderma reesei also indicated that a series of genes encoding cellulases and hemicellulases were selectively induced under various pH values . Furthermore, deeper studies revealed that a pH regulatory system could be mediated by the transcription factor of PacC, which was responsible for differential protein expression at different pH values [14, 17]. It was reported that knocking out the pacC gene could significantly affect both cellulase and xylanase activities in two model filamentous fungi, Trichoderma reesei and Neurospora crassa [18, 19]. These previous results showed that a PacC-mediated mechanism likely exists as a general way of regulating lignocellulases in filamentous fungi. However, only portions of enzymes or sporadic indirect pieces of evidence (transcriptome sequencing) were provided in those previous studies, and very few studies have tried to entirely investigate the change in all proteins in response to pH in a filamentous fungus.
Trichoderma guizhouense NJAU4742 strain, which was isolated from compost, was widely used as a functional active component of biological organic fertilizer. In this study, NJAU4742 strain was fermented at different initial pH values to investigate the effects on degradation of the sole carbon source (rice straw). In the whole process, the growth conditions and extracellular lignocellulose activities were analyzed in detail. Simultaneously, the proteomic analysis using the Sequential Window Acquisition of all Theoretical fragment ions (SWATH) technique provided a specific description of how pH affected the intracellular and extracellular proteins synthesized by NJAU4742 strain.
Pre-evaluating the ability of T. guizhouense NJAU4742 to degrade rice straw at different pH values in the submerged fermentation
In this study, spores of NJAU4742 strain (1 × 107/L) were inoculated into triangular flasks containing liquid minimal medium (MM) with a pH ranging from 2.0 to 10.0 and 2% (w/v) rice straw as the sole carbon source. During cultivation at 28 °C and 150 rpm, pH changes were determined, and results are shown in Fig. 1a. The growth of NJAU4742 strain could significantly regulate the ambient pH; for example, initial pH values of 3.0, 4.0, 5.0, 6.0 and 7.0 were all increased by an average increment of approximately one pH value after the cultivation of NJAU4742 strain for 7 days, but the initial pH values of 2.0, 9.0 and 10.0 remained stable throughout the entire cultivation period. Meanwhile, endoglucanase and xylanase activities were also evaluated for cellulose and hemicellulose degradation abilities, respectively, during the biodegradation process (Fig. 1b, c). Consistent with the pH changes, no enzyme activities were detected in the cultivations with the initial pH values of 2.0, 9.0 and 10.0. Surely, both could be explained by the fact that the initial high concentration of H+ or OH− could actually inhibit spore germination and further growth of NJAU4742 strain in the submerged fermentation (data not shown). In addition, both the highest xylanase (8.23 ± 0.20 U/mL) and exoglucanase (2.19 ± 0.06 U/mL) activities were obtained in the treatment with initial pH values of 5.0 in the third and 5th day, respectively. In general, enzyme activities, in fermentations with the initial pH values ranging from 3.0 to 7.0, had no statistical differences between each other (p > 0.05), and were all kept in a relatively high level from the 3rd day (the average xylanase activities: 3.84–5.02 U/mL, and the average exoglucanase activities: 0.87–1.41 U/mL). All these results indicated that T. guizhouense NJAU4742 strain preferred acidic environment, and the ambient pH was one of the critical parameters for the lignocellulose utilization by this strain.
The performance of NJAU4742 strain under solid-state fermentations of rice straw at different pH values
Solid-state fermentation is an important method for cultivating filamentous fungi when used in agriculture or industry. Based on the data obtained from the submerged fermentations, pH values of 2.0, 3.0, 6.0, 8.0 and 9.0 were chosen as the initial ambient pH values for the solid-state fermentations of NJAU4742 strain using rice straw as the sole carbon source. During the fermentation course of 7 days, the spores germinated and grew until large amounts of green mature spores were produced (a complete growth cycle, Additional file 1: Figure S1). In solid-state fermentation, an initial pH of 6.0 could significantly accelerate the growth of NJAU4742 strain compared to other pH values, and initial pH values of 3.0 and 8.0 could also normally incubate this strain, but with some delay in the growth when compared to the former treatment. Similar to the submerged fermentation, an initial pH of 2.0 completely inhibited the spore germination of NJAU4742 strain, and a relatively small amount of mycelia were observed in the treatment with the initial pH value of 9.0.
The diverse growth conditions showed that different ambient pH values had strong effects on the NJAU4742 strain during the solid-state fermentation process of rice straw, and thus, it was very interesting to investigate the extracellular conditions of lignocellulose degradation in this process. Therefore, the treatments with initial pHs of 3.0, 6.0 and 9.0 (SSF3.0, SSF6.0 and SSF9.0) were chosen for more systematic analysis in detail. First, we detected the releasing amount of carbon dioxide (RCD) under these conditions, and the results (Fig. 2a) showed that the RCDs of SSF3.0, SSF6.0 and SSF9.0 all increased rapidly in the first 5 days. Then, the RCDs of SSF6.0 and SSF9.0 dropped rapidly to a similar final amount (3096 ppm and 2235 ppm, respectively) on the 7th day, but the decrease in RCD of SSF3.0 happened after the 6th day. For the values of RCDs, 3.91- or 1.41-fold changes were observed between the highest values of SSF6.0 and SSF9.0 or SSF3.0. These results revealed the following facts: (1) NJAU4742 strain grew the best and had the largest amount of biomass in SSF6.0 (Additional file 1: Figure S1); (2) sporulation was the predominant reason for the decrease in RCDs; (3) the growth cycle of SSF3.0 was clearly delayed when compared to that of SSF6.0 even though the biomass was equivalent. In addition, NJAU4742 strain in SSF9.0 also had a normal growth cycle but a limited increase in biomass, and this was assumed to be led by the low utilization efficiency of rice straw. The electron microscope analysis (Fig. 2b) indicated that rice straw in SSF9.0 still largely maintained a smooth surface similar to the rice straw without bio-fermentation (control check, CK), in comparison to the total destruction of the rice straw surface structure both in SSF3.0 and SSF6.0. Together with the fact that SSF9.0 also owned the lowest enzyme activities in all treatments (Fig. 3a–c, p < 0.01), including xylanase, endoglucanase and exoglucanase activities, and these results indicated that low utilization efficiency of rice straw in SSF9.0 seemed to be resulted by the low enzymatic hydrolysis of rice straw. To verify this viewpoint as mentioned above, we detected the concentrations of extracellular proteins (CEP) in all three treatments (Fig. 3d), and SSF9.0 indeed had the lowest values, but the differences (79% and 67% for SSF3.0 and SSF6.0, respectively) are not enough to completely support the hypothesis that a low amount of secreted enzymes led to an inefficient destruction of rice straw in SSF9.0. Therefore, it was hypothesized that the ambient pH of SSF9.0 was another main reason that largely inhibited all enzyme activities, as it is well known that most secreted cellulases and hemicellulases can only maintain approximately 20% of enzyme activity at pH 9.0 [20,21,22]. The crude enzymes of NJAU4742 strain induced by rice straw were detected to only have 19.13 ± 0.31% of xylanase activity and 1.12 ± 0.07% endoglucanase activity at pH 9.0 when compared to pH 5.0. Therefore, it is necessary to detect the ambient pH in SF9.0, and the results showed that there was no obvious change in pH during the 7-day fermentation (Fig. 3e). This high pH of 9.0 extremely affected extracellular enzyme activities and stabilities, and thereby strongly contributed to the growth limitation of high pH. Interestingly, 3.0 was also a value far from the optimal pH of most lignocellulases, but the rice straw in SSF3.0 was effectively utilized. The change profile of pH value in SSF3.0 illustrated that the pH value was rapidly and largely increased as the fermentation continued until reaching a final value of 5.8 (Fig. 3e), which was very close to the optimal pH for most lignocellulases secreted by acidic filamentous fungi. Finally, the pH increases in SSF6.0 was also observed, which indicated that the growth of NJAU4742 strain could significantly eliminate the high concentration of H+ in two ways: (1) neutralized by the intermediate products released from the rice straw; (2) absorbed and utilized by NJAU4742 strain, or offset by the secretions from NJAU4742 strain. In fact, the results in Fig. 3f excluded the first mechanism, and thus indicated that the second mechanism is responsible for the pH increase during the solid-state fermentation process.
Proteomic analysis for solid-state fermentations of NJAU4742 strain at pH 3.0, 6.0 and 9.0
The different growth conditions and enzymatic properties above showed that diverse regulations must exist in the solid-state fermentations of NJAU4742 strain at different initial pH values, thus proteomic analysis was used to investigate the details of responding to ambient pH by NJAU4742 strain. After protein extraction, the mycelial proteins and extracellular proteins were mixed based on the protein concentration (mg g−1 dw, protein/substrates), meaning that both were extracted from the same amount of the fermented materials (see “Methods”). SDS-PAGE analysis showed that the mixed total proteins were repeatable in each treatment and that significant differences existed among SSF3.0, SSF6.0 and SSF9.0 (Additional file 1: Figure S2). Consistent with the highest enzyme activities, total proteins of SSF6.0 had the most abundant protein bonds, as shown in Additional file 1: Figure S2. After quality control, total proteins from SSF3.0, SSF6.0 and SSF 9.0 were separately detected and analyzed using the SWATH technique, which resulted in a total number of 1139 proteins when using the whole protein sequences of NJAU4742 strain as the database.
For all the determined proteins, their coverage rate of peptides was calculated and seemed to fit a noncentral F-distribution (Additional file 1: Figure S3A), with approximately 90% of identified proteins having a coverage rate of peptides lower than 37%. The protein abundances in all samples (SSF3.0, SSF6.0 and SSF9.0) conformed to a normal distribution with a range of 103–108 (Additional file 1: Figure S3B). Together with the considerable repeatability of each sample (Additional file 1: Figure S3C) and objective statistical information produced with SWATH analysis (data not shown), these data of the identified proteins allowed the deep analysis of the responses to different ambient pH values by NJAU4742 strain.
Protein location analysis
All 1139 identified proteins were classified based on their locations (Fig. 4 and Additional file 2), in which 921 proteins were annotated as intracellular proteins and 346 proteins therein were located in different detailed intracellular regions, such as the ribosome (58 proteins), endoplasmic reticulum (31 proteins), cytoplasm/cytosol (82 proteins), golgi apparatus (19 proteins), mitochondria (59 proteins), and nucleus (63 proteins). Additionally, 265 proteins were annotated as extracellular proteins, including 94 glycoside hydrolases, 11 carbohydrate esterases, 11 auxiliary activities enzymes, 27 proteases, 21 small secreted cysteine-rich proteins (SSCRPs), 7 nucleases, 44 functional proteins and 50 unknown proteins. These identified proteins had a very uneven distribution of expression abundances. For example, though less than a quarter of the total number of proteins, the 265 extracellular proteins amazingly shared > 58.2% of total protein abundances, and did not conform to the normal distribution. In contrast, most intracellular organelles fitted a much greater distribution of spindle morphology in protein abundance when compared to extracellular proteins (Fig. 4). These findings indicated that the secreted proteins from NJAU4742 strain did not follow intracellular rules, but continually secreted enzymes to achieve much more nutrition by degrading rice straw in a resource-poor environment. From this aspect, glycoside hydrolases and proteases evidently played important roles in their > 72% of extracellular protein abundances, which showed the critical need for carbon and nitrogen sources by strain NJAU4742 strain.
Rice straw, a natural lignocellulose, mainly consisted of cellulose, hemicellulose and lignin. In solid-state fermentation, rice straw was efficiently deconstructed by the secreted enzymes, which were usually classified into different families in the CAZyme database, including glycoside hydrolases (GHs), polysaccharide lyases (PLs), carbohydrate esterases (CEs), auxiliary activities (AAs) and carbohydrate-binding modules (CBMs). Unfortunately, no PLs were identified in this experiment, even though they have positive functions in degrading rice straw. In Fig. 5, Additional file 1: Figure S4 and Additional file 3, 190 identified CAZymes (GHs, CEs, AAs and CBMs), and 510 other KEGG pathway-annotated proteins were used to construct a protein network, which significantly showed the functions of pH in regulating these identified proteins. For enzymes in the CEs and AAs, the treatment of SSF6.0 had a relatively higher inductive effects when compared to those of the treatments of SSF3.0 and SSF9.0. For enzymes in the GHs, however, the result became complex. When comparing SSF3.0 and SSF6.0, differential proteins were divided into two numerically comparable groups (the up-regulated and the down-regulated, 1.5-fold change and p < 0.05), which meant that SSF3.0 or SSF6.0 could selectively enhance or inhibit the inductions of some enzymes in GHs. This phenomenon did not happen in the comparison between SSF9.0 and SSF6.0, where nearly all differential proteins had high abundances in SSF6.0.
The identified CAZymes included cellulases, hemicellulases, chitinases and other functional enzymes (Fig. 6). Among cellulases, 5 endoglucanases (families: GH5, GH7, GH12 and GH16), 2 exoglucanases (GH6 and GH7), 8 β-glucosidases (GH1, GH3 and GH5) and 3 polysaccharide monooxygenases (AA9 and AA11) were identified. For hemicellulases, 8 endoxylanases (GH10, GH11 and GH30), 6 β-xylosidases (GH3, GH39 and GH43), 8 hydrolytic enzymes of xylan side chains (GH54, GH67, CE5 and CE15), 10 mannanases (GH2, GH5, GH38 and GH92) and 14 other functional hemicellulases (GH2, GH27, GH28, GH30, GH35, GH36, GH37, GH62, GH65, GH74 and GH115) were identified. Additionally, 11 chitinases (GH18, GH20 and GH75) and 113 other enzymes (including glucanases, glucosidases, esterases, and unknown proteins) were also identified. As shown in Fig. 6, the up different inductions by pH values were also reflected in these subgroups (cellulases and hemicellulases) of GHs when comparing SSF3.0 and SSF6.0. That is, some cellulases/hemicellulases were highly expressed in SSF3.0, and some other cellulases/hemicellulases were highly expressed in SSF6.0. For those cellulases/hemicellulases down-regulated in SSF3.0, they were also mostly down-regulated in SSF9.0 when compared to SSF6.0 (1.5-fold change and p < 0.05). However, it is different for chitinases, where SSF3.0 seemed to be the best induction condition in all treatments. For these rice straw-degrading-related CAZymes, different induction performances evidently existed, including (1) highest expression abundances in SSF3.0; (2) highest expression abundances in SSF6.0; (3) lowest expression abundances in SSF6.0 and (4) no obvious differences at all pH values. All these conditions clearly indicated that different ambient pH values functioned to regulate the secretion of various lignocellulases.
When focusing on the protein abundances of those secreted CAZymes, very significant differences existed (Fig. 7). Among cellulases, the most abundant enzymes, including 2 exoglucanases, 1 β-glucosidase and 1 polysaccharide monooxygenase, composed 15.6% of all extracellular proteins in protein abundance. Among hemicellulases, four abundant enzymes (3 endoxylanases and 1 α-N-arabinofuranosidase) composed 16.2% of all extracellular proteins in protein abundance. Meanwhile, a xyloglucanase and a α-1,3-glucanase had percentages of 4.1% and 2.4%, respectively. These results reflected that some enzymes in CAZymes were the universal choices among filamentous fungi, while their largely unequal secretions possibly satisfied the need for rapid degradation of rice straw. Certainly, protein abundances could not decide the enzyme functions, but their differential distribution here actually showed that a natural proportion of these secreted functional enzymes was reasonable for degradation of the lignocellulose.
In solid-state fermentation, NJAU4742 strain had adapted to a wide environmental pH condition, ranging from 3.0 to 9.0, and showed a distinct preference for the specific pH of 6.0. Beyond adaptation, NJAU4742 strain could also significantly reconstruct the environmental pH condition by alkalinization, in which the pH value was elevated from 3.0 to 6.0 in 5 days. Therefore, it is very interesting to discuss NJAU4742’s alkalinization, which is hence hypothesized to be the result of three possible mechanisms: (I) the released components from the degradation of rice straw neutralized the acidic environmental pH; (II) NJAU4742 strain absorbed extracellular H+ and (III) the secretions of NJAU4742 strain elevated the environmental pH. For hypothesis I, the main components of rice straw contained polysaccharides (cellulose, xylan and mannan) and lignin, however, the released basic units, such as oligosaccharides, glucose, xylose, mannose and aromatic alcohols, actually had little influence on the ambient pH. The components that could alter the ambient pH were pectin and those substituent groups of hemicellulose, of which the complete hydrolysis would release various organic acids, including galacturonic acid, glucuronic acid, ferulic acid, p-coumalic acid, and acetic acid . Meanwhile, protein degradation by extracellular proteases could also release various amino acids from the substrate. Apparently, these produced organic acids will function as a reverse to the environmental alkalinization by NJAU4742 strain, and were probably absorbed and utilized rapidly by NJAU4742 strain, much like the produced sugars, and this behavior functioned in regulating ambient pH, which belonged to part of hypothesis II (absorbing extracellular acids). Additionally, the original high concentration of H+ in SSF3.0 strongly drew attention to an enzyme called the H+-ATPase, which was identified and significantly regulated by pH in this study (OPB40853, p < 0.05, Additional file 2). This enzyme was generally located on the plasma membrane of plants and fungi, and belonged to the P-type cation-translocating ATPase family. Many previous studies showed that H+-ATPases were always regulated by glucose, an acidic pH and weak organic acids , which existed in solid-state fermentations of NJAU4742 strain and resulted in the highest expression level of the H+-ATPase in SSF6.0. The complex regulation of the H+-ATPase in the SSF of NJAU4742 strain suggested its relationship with NJAU4742 strain alkalinization, but may also work as a counter for its reported function of pumping H+ out of the cell to create a transmembrane electrochemical proton gradient for nutrient uptake and stabilizing intracellular pH [24, 25]. Actually, under such a low pH of 3.0 condition (high natural potential across plasma membrane), it was hard to know whether the H+-ATPase was still needed to work or if there were any possible patterns of the H+-ATPase helping NJAU4742 strain absorb H+ in some way, which still required a more in-depth study. For hypothesis III, many studies have proposed that fungi could modify the pH of their environment by secreting acids or alkalis. For example, Sclerotinia sclerotiorum and Botrytis sp. could produce oxalic acid , while Penicillium sp. and Aspergillus sp. can secrete mainly gluconic and citric acids [27, 28]. Although pH was slightly dropped in SSF9.0, it was still difficult to judge whether NJAU4742 strain secreted acids in this study, but at least this was not the main determinant for the acidic environment. Some other fungi, such as Colletotrichum gloeosporioides , Dictyostelium discoideum  and Debaryomyces hansenii , and especially for Candida albicans , could alkalinize their environment by secreting ammonia, a byproduct during the amino acid metabolic process. By a BLAST search, the homologs of an ammonium transporter of MEPB and three ammonia transporters (ATOs) in this process were found in strain NJAU4742, but no homologs were detected for the critical regulator (Stp2p) and amino acid catabolic enzyme (Ach1p). Among the four detected homologs, only one ATO protein (OPB44235) was identified by SWATH, but its highest expression level in SSF6.0 could not be explained logically. These previous reports of other fungi’s alkalinization supported that NJAU4742 strain in this study might elevate the ambient pH by secreting another kind of alkali (not ammonia), which was entirely possible because C. albicans also shared multiple mechanisms of environmental alkalinization not only ammonia secretion [32,33,34]. Certainly, H2O2 and secondary metabolites were also considered as the possible functional parts of the complex environmental alkalinization of NJAU4742 strain.
In SSF3.0, the growth inhibition of NJAU4742 strain by the low pH was relieved as the ambient pH increased, which indicated that environmental alkalinization determined the growth ability of NJAU4742 strain at pH 3.0. This fact reflects that a complex regulatory system exists for detecting different pH values and further, responding to the pH. The results here clearly showed that extracellular proteins, such as cellulases, hemicellulases and chitinases, were under the regulation of this possible system. Similarly, some cellulase and xylanase genes were also reported to be controlled by the ambient pH in T. reesei , Humicola grisea  and Neurospora crassa , and thus together showed the universality of pH regulation of lignocellulases in filamentous fungi. Actually, there are already some reports that the integrity of seven proteins, designated PacC, PalA, PalB, PalC, PalF, PalH and PalI, were responsible for gene regulation at different pH values . Under acidic growth conditions, PacC72 (72 kDa) is inactive without any known role in transcriptional regulation, but it is converted to PacC53 (53 kDa) and further to PacC27 (27 kDa), which is subsequently active and functions in the activation of alkaline-expressed genes and repression of acid-expressed genes. PacC orthologs were found to modulate the expression of cellulolytic genes in response to ambient pH in filamentous fungi, but the effects varied depending on the organisms. For example, deletion of the pacC gene in T. reesei significantly increased the cellulolytic enzyme activities at neutral pH , but the transcriptome analysis showed only a few cellulase- and hemicellulase-encoding genes were clearly under PacC regulation . In Aspergillus nidulans , PacC significantly upregulated the expression of cellulolytic genes at alkaline pH, but had a reverse observation in N. crassa . The homologs of the seven proteins in the whole PacC-mediated system could be detected and encoded in the genome of NJAU4742 strain, and thus suggested that those pH-regulated lignocellulases were possibly under the control of this system. Although none of the homologs appeared in the 1139 identified proteins, further study of PacC-mediated and other possible systems in regulating extracellular enzymes by the ambient pH, will be very interesting.
Obviously, filamentous fungi regulated the secretion of extracellular enzymes by the ambient pH, however, we need to understand the biological significance behind this behavior. The organisms are highly evolved, and thus the secretion adjustment of extracellular enzymes at different pH values is believed to be an adaptive mechanism to ensure better survival in that environment. That means, NJAU4742 strain possibly achieved more efficient degradation of rice straw by adjusting the differential secretion of cellulases and hemicellulases when exposed under the different pH values. Logically, the adjustment is to achieve the optimal synergy of extracellular enzymes under different pH conditions, which possibly could be explained from the different enzymatic properties of cellulases and hemicellulases. And, it is really the fact that cellulases and xylanases from the same filamentous fungus have different enzymatic properties, such as thermostability, optimal temperature, specific activity and functional domains . Just like extracellular proteases, they were distinguished as acidic, neutral and alkaline proteases and regulated by the PacC-mediated system . The secreted proteases of NJAU4742 strain also showed significant pH regulation in this study (Additional file 2). Besides, Xiong et al.  reported that xylanase I was most secreted and active at low pH value (4.0), while xylanase III was most secreted and active at high pH value (6.0) in T. reesei. Xylanase II was clearly produced at both pH values. These results and analysis clearly support such a judgment that the pH-dependent regulation of extracellular enzymes is probably for achieving higher enzymatic synergy by secreting the optimal enzymes at different pH values, so as to obtain a most efficient degradation of lignocellulose. Certainly, more detailed and in-depth researches need to be carried out to explain this pH-dependent synergy of lignocellulases.
Additionally, we also noticed the stable differential distribution of extracellular enzyme abundances, which was believed to be the result of evolution, as well as a kind of synergistic performance of multiple enzymes. For cellulose degradation, exoglucanase attacked the natural and LPMO-produced ends of crystalline cellulose, and then the products could be hydrolyzed into glucose by β-glucosidase. The high abundance of these three functional enzymes suggested the specific enhancement of the main steps in the degradation of crystalline cellulose, but not including the endoglucanases, which seemed to have difficulties when attacking the middle of the cellulose chain in the crystalline structure. In contrast, endoxylanase was the largest abundant secreted enzyme among all xylanases, corresponding to the fact that endoxylanase could easily attack the whole chain of xylan for its non-crystalline structure. Facing the complex 3D structure of lignocellulose, the differential distribution of extracellular enzymes is indeed likely to be for greater synergy efficiency, and further analyzing this synergy will be very interesting and valuable in biofuel industry.
This study described the performances of NJAU4742 strain during the submerged and solid-state fermentation processes at different initial pH values in detail. The results showed that ambient pH dominated the growth, sporulation, extracellular enzyme activities and decomposition rate of rice straw by NJAU4742 strain. Proteomic analysis deeply revealed the intracellular and extracellular proteins during the solid-state fermentation processes and clearly quantified 190 functional enzymes for rice straw degradation especially for their extracellular synergetic distribution. More importantly, the pH-regulated lignocellulose-degrading enzymes in filamentous fungi were firstly distinguished and classified in detail. All these results contribute to the theoretical basis for the degradation of plant biomass by filamentous fungi in the biofuel and biological organic fertilizer industries.
Strains and culture conditions
Trichoderma guizhouense NJAU4742 was isolated from a mature compost sample and stored in 15% glycerol at − 80 °C in our lab, and its genome sequence was already published in the NCBI database (Accession No. LVVK00000000.1). Murashige and Skoog (MS) medium (1.36 g/L KH2PO4, 2.13 g/L Na2HPO4, 0.2 g/L MgSO4·7H2O, and trace elements) was used as the liquid medium for cultivating NJAU4742 in both submerged and solid-state fermentation, but varied in percentages of rice straw, which were 2% (w/v) or 25% (w/v), respectively. The rice straw used in this study was obtained from the local farmland (Nanjing, China) and then washed with running water for 48 h, after which it was dried entirely at 40 °C for 5 h and ball-milled into powder using a ball-miller (Retsch MM400, Germany) with ZrO2 balls (15 mm, diameter) and vessel (35 mL, volume), and the parameters were set as follows: 15 Hz, 1 min milling followed by a 30-s pause, with a total time of 5 min, after which the rice straw powder was screened by a 30 mesh sieve with the average size less than 0.6 mm. The initial pH values of the submerged or solid-state fermentations were adjusted by adding 0.5 M H2SO4 or NaOH and monitored with a pH meter (PB-10, Sartorius, Germany).
Spores of NJAU4742 strain were collected from the incubated PDA plates by washing followed by the filtration through 4 layers of gauze. After two repeated washing by the sterilized ddH2O, the spore concentration was subsequently adjusted to be 1 × 107/mL through the method of hemocytometer counting. For the submerged fermentation, a total of 5 × 106 spores were inoculated into 50 mL MS medium with 2% (w/v) rice straw in a 250-mL triangular flask and grown under the condition of 28 °C and 150 rpm. For the solid-state fermentation, a total of 1 × 107 spores were inoculated into a 500-mL triangular flask containing 75 mL MS medium and 25 g rice straw, and grown under the condition of 28 °C. All treatments were performed for total 7 days, and sampled at the incubation times of 1, 2, 3, 4, 5, 6 and 7 days. Three biological repeats existed for each sampling point.
Enzyme activity assays
Endoglucanase, exoglucanase and xylanase activities were measured using carboxymethyl cellulose sodium (CMC-Na) (Sinopharm Chemical Reagent Co., China), p-nitrophenyl-β-d-cellobiose (pNPC) (Sigma, USA), and oat spelts xylan (Sigma, USA) as the substrate, respectively, according to the methods described previously [20, 41]. One unit of enzyme activity was defined as the amount of enzyme required to release 1 μmol reducing sugars or pNP from the substrate in 1 min. Protein concentrations were determined using the Bradford protein assay kit (ThermoFisher, USA) based on the product description . All the enzyme activities were reported as the means of at least three replicates.
Electron microscope analysis and respiratory rate determination
After fixation with 2.5% glutaraldehyde for 2 days, the samples from different treatments were immediately dehydrated in 10 mL of ice-cold 200 proof ethanol (Sigma-Aldrich), point-dried in a critical point dryer (HCP-2, Hitachi High-Technologies Corporation, Japan) and coated with 60% Au/Pd in a sputter coater (Sputter Coater Baltec SCD500, Bal-Tel) . The surface morphologies of different samples were observed by using an S-4800 II field emission scanning electron microscope (Hitachi, Japan). The respiratory rate of NJAU4742 strain during the solid-state fermentation process was evaluated by tracking the release of CO2. Briefly, the triangular flasks of solid-state fermentations were sealed with a parafilm 2 h before sampling. When sampling, 20 mL of gas in the triangular flask was extracted with a syringe and then injected into the gas sampling bag (E-SWITCH, China), and the contents of CO2 was determined by a gas chromatography (Agilent 7890A) equipped with Porapak Q column and a flame ionization detector (FID) according to Zhang et al.  with some modification: 500 µL of gas were injected through a septum with the temperature of 80 °C, and the carrier gas (N2) flow-rate was 30 mL min−1; meanwhile, the temperature of the methanizer and detector were set as 450 °C and 250 °C, respectively.
Protein extraction and SWATH analysis
Samples of each treatment in solid-state fermentations at the 4th day were taken and used to extract proteins for proteome analysis. Extracellular protein extraction was carried out as follows: 150 mL deionized ddH2O was added into triangular flasks, which were then shaken at 180 rpm for 1 h, and the supernatants were collected by centrifugation at 10,000 rpm for 10 min and filtration through a 0.22 µm sterile membrane. The precipitates of the former step were used to prepare the mycelial intracellular proteins, which were extracted by using the NoviPure® Soil Protein Extraction Kit (MOBIO, USA) according to the manufacturer’s protocol. After that, extracellular and intracellular proteins, both extracted from the same amount of fermentation materials samples (mg g−1 dw, protein/substrates), were directly mixed and concentrated by lyophilization and then stored at − 80 °C. Two biological repeats existed for each condition in proteomic analysis.
The whole SWATH analysis was commissioned to a specialized company (GeneCreate, Wuhan, China), and the analysis is described briefly as follows: 100 μg of proteins in each sample were digested by trypsin (Promega, USA) at 37 °C overnight, separated by the Strata-X C18 pillar and dried using a vacuum concentration meter. LC–ESI-MS/MS analysis was performed on an AB SCIEX nanoLC–MS/MS (Triple TOF 5600 plus) system in two phases: data-dependent acquisition (DDA) was followed by SWATH acquisition, where all samples were mixed and detected in DDA, and the resulting data were used as a library for analysis of each sample by SWATH. Spectral library generation and SWATH data processing were performed using skyline version 3.5 software  complying with the following rules: (i) peptides containing modifications and/or shared by different protein entries/isoforms were excluded from selection; (ii) peptides having a ProteinPilot-identified confidence less than 95% were excluded; (iii) up to five fragment ions that were ranked by intensity were chosen; (iv) fragment ions within the SWATH isolation window were excluded from selection; and (v) to control false discovery rate, a random mass shift of Q1 and Q3 m/z strategy  was used to create a decoy spectra library for targeted peak extraction, a mass-to-charge tolerance of at most 10 ppm (0.05 m/z) was allowed for both peptide precursor and fragment ion. Fragment ion areas that belonged to one peptide were added to obtain a peptide’s abundance, and the total abundance of peptides for a given protein was determined to obtain the protein’s abundance. To eliminate the random errors and sample bias, all the data among samples were normalized using the median normalization method , and the mProphet algorithm was used both for sample normalization and for assessing the data confidence .
CAZymes annotation was carried out using the hmmscan tool from HMMER software  to search the dbCAN database , and the output was processed by the script of hmmscan-parser.sh (https://github.com/carden24/Bioinformatics_scripts/blob/master/hmmscan-parser.sh). KEGG pathway annotation was performed using KEGG Orthology-Based Annotation System 2.0 (KOBAS) . Gene location was manually annotated by combining information from GO Ontology, KEGG annotation, NCBI gene descriptions and some previous articles. Signal peptide was predicted using SignalP 4.0 software . Circos  and Cytoscape  were used for data visualization. Violin plot was performed using the R package vioplot .
Availability of data and materials
the releasing amount of carbon dioxide
the concentrations of extracellular proteins
Kubicek CP, Kubicek EM. Enzymatic deconstruction of plant biomass by fungal enzymes. Curr Opin Chem Biol. 2016;35:51–7.
Himmel ME, Ding S-Y, Johnson DK, Adney WS, Nimlos MR, Brady JW, Foust TD. Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science. 2007;315:804–7.
Mosier N, Wyman C, Dale B, Elander R, Lee YY, Holtzapple M, Ladisch M. Features of promising technologies for pretreatment of lignocellulosic biomass. Bioresour Technol. 2005;96:673–86.
Limayem A, Ricke SC. Lignocellulosic biomass for bioethanol production: current perspectives, potential issues and future prospects. Prog Energy Combust Sci. 2012;38:449–67.
Nishiyama Y, Sugiyama J, Chanzy H, Langan P. Crystal structure and hydrogen bonding system in cellulose Iα from synchrotron X-ray and neutron fiber diffraction. J Am Chem Soc. 2003;125:14300–6.
Nishiyama Y, Langan P, Chanzy H. Crystal structure and hydrogen-bonding system in cellulose Iβ from synchrotron X-ray and neutron fiber diffraction. J Am Chem Soc. 2002;124:9074–82.
Gírio FM, Fonseca C, Carvalheiro F, Duarte LC, Marques S, Bogel-Łukasik R. Hemicelluloses for fuel ethanol: a review. Bioresour Technol. 2010;101:4775–800.
Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010;153:895–905.
Gupta VK, Kubicek CP, Berrin J-G, Wilson DW, Couturier M, Berlin A, Filho EXF, Ezeji T. Fungal enzymes for bio-products from sustainable and waste biomass. Trends Biochem Sci. 2016;41:633–45.
Chen L, Yang X, Raza W, Luo J, Zhang F, Shen Q. Solid-state fermentation of agro-industrial wastes to produce bioorganic fertilizer for the biocontrol of Fusarium wilt of cucumber in continuously cropped soil. Bioresour Technol. 2011;102:3900–10.
Zheng Y, Zhao J, Xu F, Li Y. Pretreatment of lignocellulosic biomass for enhanced biogas production. Prog Energy Combust Sci. 2014;42:35–53.
Liu H, Chen D, Zhang R, Hang X, Li R, Shen Q. Amino acids hydrolyzed from animal carcasses are a good additive for the production of bio-organic fertilizer. Front Microbiol. 2016;7:1290.
Guo JH, Liu XJ, Zhang Y, Shen JL, Han WX, Zhang WF, Christie P, Goulding KWT, Vitousek PM, Zhang FS. Significant acidification in major Chinese croplands. Science. 2010. https://doi.org/10.1126/science.1182570.
Peñalva MA, Tilburn J, Bignell E, Arst HN Jr. Ambient pH gene regulation in fungi: making connections. Trends Microbiol. 2008;16:291–300.
Li B, Wang W, Zong Y, Qin G, Tian S. Exploring pathogenic mechanisms of Botrytis cinerea secretome under different ambient pH based on comparative proteomic analysis. J Proteome Res. 2012;11:4249–60.
Häkkinen M, Sivasiddarthan D, Aro N, Saloheimo M, Pakula TM. The effects of extracellular pH and of the transcriptional regulator PACI on the transcriptome of Trichoderma reesei. Microb Cell Factories. 2015;14:63.
Arst HN Jr, Peñalva MA. pH regulation in Aspergillus and parallels with higher eukaryotic regulatory systems. Trends Genet. 2003;19:224–31.
He R, Ma L, Li C, Jia W, Li D, Zhang D, Chen S. Trpac1, a pH response transcription regulator, is involved in cellulase gene expression in Trichoderma reesei. Enzyme Microb Technol. 2014;67:17–26.
Antoniêto ACC, Pedersoli WR, dos Santos Castro L, da Silva Santos R, da Silva Cruz AH, Nogueira KMV, Silva-Rocha R, Rossi A, Silva RN. Deletion of pH regulator pac-3 affects cellulase and xylanase activity during sugarcane bagasse degradation by Neurospora crassa. PLoS ONE. 2017;12:e0169796.
Miao Y, Li J, Xiao Z, Shen Q, Zhang R. Characterization and identification of the xylanolytic enzymes from Aspergillus fumigatus Z5. BMC Microbiol. 2015;15:126.
Miao Y, Kong Y, Li P, Li G, Liu D, Shen Q, Zhang R. Effect of CBM1 and linker region on enzymatic properties of a novel thermostable dimeric GH10 xylanase (Xyn10A) from filamentous fungus Aspergillus fumigatus Z5. AMB Express. 2018;8:44.
Liu D, Zhang R, Yang X, Wu H, Xu D, Tang Z, Shen Q. Thermostable cellulase production of Aspergillus fumigatus Z5 under solid-state fermentation and its application in degradation of agricultural wastes. Int Biodeterior Biodegrad. 2011;65:717–25.
Portillo F. Regulation of plasma membrane H+-ATPase in fungi and plants. Biochim Biophys Acta BBA Rev Biomembr. 2000;1469:31–42.
Palmgren M, Morsomme P. The plasma membrane H+-ATPase, a simple polypeptide with a long history. Yeast. 2018;36(4):201–10.
Hesse SJ, Ruijter GJG, Dijkema C, Visser J. Intracellular pH homeostasis in the filamentous fungus Aspergillus niger. Eur J Biochem. 2002;269:3485–94.
Manteau S, Abouna S, Lambert B, Legendre L. Differential regulation by ambient pH of putative virulence factor secretion by the phytopathogenic fungus Botrytis cinerea. FEMS Microbiol Ecol. 2003;43:359–66.
Prusky D, Yakoby N. Pathogenic fungi: leading or led by ambient pH? Mol Plant Pathol. 2003;4:509–16.
Ruijter GJ, Visser J. Characterization of Aspergillus niger phosphoglucose isomerase. Use for quantitative determination of erythrose 4-phosphate. Biochimie. 1999;81:267–72.
Shnaiderman C, Miyara I, Kobiler I, Sherman A, Prusky D. Differential activation of ammonium transporters during the accumulation of ammonia by Colletotrichum gloeosporioides and its effect on appressoria formation and pathogenicity. Mol Plant Microbe Interact. 2013;26:345–55.
Schindler J, Sussman M. Ammonia determines the choice of morphogenetic pathways in Dictyostelium discoideum. J Mol Biol. 1977;116:161–9.
Gori K, Mortensen HD, Arneborg N, Jespersen L. Ammonia production and its possible role as a mediator of communication for Debaryomyces hansenii and other cheese-relevant yeast species. J Dairy Sci. 2007;90:5032–41.
Vylkova S, Carman AJ, Danhof HA, Collette JR, Zhou H, Lorenz MC. The fungal pathogen Candida albicans autoinduces hyphal morphogenesis by raising extracellular pH. MBio. 2011;2:e00055-11.
Danhof HA, Vylkova S, Vesely EM, Ford AE, Gonzalez-Garay M, Lorenz MC. Robust extracellular pH modulation by Candida albicans during growth in carboxylic acids. MBio. 2016;7:e01646-16.
Naseem S, Araya E, Konopka JB. Hyphal growth in Candida albicans does not require induction of hyphal-specific gene expression. Mol Biol Cell. 2015;26:1174–87.
Bailey MJ, Buchert J, Viikari L. Effect of pH on production of xylanase by Trichoderma reesei on xylan- and cellulose-based media. Appl Microbiol Biotechnol. 1993;40:224–9.
Mello-de-Sousa TM, Silva-Pereira I, Poças-Fonseca MJ. Carbon source and pH-dependent transcriptional regulation of cellulase genes of Humicola grisea var. thermoidea grown on sugarcane bagasse. Enzyme Microb Technol. 2011;48:19–26.
Kunitake E, Hagiwara D, Miyamoto K, Kanamaru K, Kimura M, Kobayashi T. Regulation of genes encoding cellulolytic enzymes by Pal-PacC signaling in Aspergillus nidulans. Appl Microbiol Biotechnol. 2016;100:3621–35.
Miao Y, Li P, Li G, Liu D, Druzhinina IS, Kubicek CP, Shen Q, Zhang R. Two degradation strategies for overcoming the recalcitrance of natural lignocellulosic xylan by polysaccharides-binding GH10 and GH11 xylanases of filamentous fungi. Environ Microbiol. 2017;19:1054–64.
Peñalva MA, Arst HN. Recent advances in the characterization of ambient pH regulation of gene expression in filamentous fungi and yeasts. Annu Rev Microbiol. 2004;58:425–51.
Xiong H, von Weymarn N, Leisola M, Turunen O. Influence of pH on the production of xylanases by Trichoderma reesei Rut C-30. Process Biochem. 2004;39:731–6.
Liu D, Li J, Zhao S, Zhang R, Wang M, Miao Y, Shen Y, Shen Q. Secretome diversity and quantitative analysis of cellulolytic Aspergillus fumigatus Z5 in the presence of different carbon sources. Biotechnol Biofuels. 2013;6:149.
Stoscheck CM.  Quantitation of protein. In: Deutscher MP, editor. Methods in enzymology. Boca Raton: Academic Press; 1990. p. 50–68.
Wang M, Ma L, Kong Z, Wang Q, Fang L, Liu D, Shen Q. Insights on the aerobic biodegradation of agricultural wastes under simulated rapid composting conditions. J Clean Prod. 2019;220:688–97.
Zhang L, Wang H, Zou J, Rogers WE, Siemann E. Non-native plant litter enhances soil carbon dioxide emissions in an invaded annual grassland. PLoS ONE. 2014;9:e92301.
Sterling C, Crouch R, Russell DJ, Calderón AI. 1H-NMR quantification of major saccharides in açaí raw materials: a comparison of the internal standard methodology with the absolute intensity qNMR method. Phytochem Anal. 2013;24:631–7.
Röst HL, Rosenberger G, Navarro P, Gillet L, Miladinović SM, Schubert OT, Wolski W, Collins BC, Malmström J, Malmström L, Aebersold R. OpenSWATH enables automated, targeted analysis of data-independent acquisition MS data. Nat Biotechnol. 2014;32:219–23.
Oberg AL, Mahoney DW. Statistical methods for quantitative mass spectrometry proteomic experiments with labeling. BMC Bioinform. 2012;13:S7.
Reiter L, Rinner O, Picotti P, Hüttenhain R, Beck M, Brusniak M-Y, Hengartner MO, Aebersold R. mProphet: automated data processing and statistical validation for large-scale SRM experiments. Nat Methods. 2011;8:430.
Finn RD, Clements J, Arndt W, Miller BL, Wheeler TJ, Schreiber F, Bateman A, Eddy SR. HMMER web server: 2015 update. Nucleic Acids Res. 2015;43:W30–8.
Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40:W445–51.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li C, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785.
Krzywinski MI, Schein JE, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.
Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2010;27:431–2.
Hintze JL, Nelson RD. Violin plots: a box plot-density trace synergism. Am Stat. 1998;52:181–4.
This research was financially supported by the National Natural Science Foundation of China (31801935), the Jiangsu Provincial Natural Science Foundation (BK20180538), and Special Fund for Agro-scientific Research in the Public Interest (201503110). QS was supported by the Innovative Research Team Development Plan of the Ministry of Education of China (IRT_17R56).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
SWATH results of the identified proteins in SSF3.0, SSF6.0 and SSF9.0. Figure S1. Growth conditions of NJAU4742 strain in SSF2.0, SSF3.0, SSF6.0, SSF8.0 and SSF9.0. Figure S2. 2D SDS-PAGE of the total proteins in SSF3.0, SSF6.0 and SSF9.0 for SWATH analysis. Figure S3. Statistical analysis of SWATH results. After SWATH detection, Frequency distribution of protein coverages and protein abundances were shown in (A) and (B), respectively. Protein abundance from each repeat was plotted and compared in (C). Figure S4. Protein expression differences between SSF9.0 and SSF6.0 visualized as a Cytoscape interaction network.
The identification and quantification results of all the proteins from different treatments.
The raw data for the Cytoscape interaction network analysis.
About this article
Cite this article
Miao, Y., Chen, X., Li, T. et al. Proteomic analysis reflects an environmental alkalinization-coupled pH-dependent mechanism of regulating lignocellulases in Trichoderma guizhouense NJAU4742. Biotechnol Biofuels 13, 6 (2020). https://doi.org/10.1186/s13068-020-1651-0
- Plant biomass
- Environmental alkalinization