- Open Access
Predominance of Trichoderma and Penicillium in cellulolytic aerobic filamentous fungi from subtropical and tropical forests in China, and their use in finding highly efficient β-glucosidase
Biotechnology for Biofuels volume 7, Article number: 107 (2014)
Cellulose is the most abundant biomass on earth. The major players in cellulose degradation in nature are cellulases produced by microorganisms. Aerobic filamentous fungi are the main sources of commercial cellulase. Trichoderma reesei has been explored extensively for cellulase production; however, its major limitations are its low β-glucosidase activity and inefficiency in biomass degradation. The aim of this work was to isolate new fungal strains from subtropical and tropical forests in China, which produce high levels of cellulase in order to facilitate development of improved commercial cellulases.
We isolated 305 fungal strains from 330 samples collected from subtropical and tropical virgin forests in China. Of these, 31 strains were found to have Avicelase activity of more than 0.2 U/ml in liquid batch cultivation. Molecular analyses of the 31 strains based on internal transcribed spacer sequences revealed that 18 were Trichoderma and 13 were Penicillium species. The best-performing isolate was Trichoderma koningiopsis FCD3-1, which had similar Avicelase activity to T. reesei Rut-C30. Most interestingly, strain FCD3-1 exhibited extracellular β-glucosidase activity of 1.18 U/ml, which was approximately 17 times higher than that of Rut-C30. One β-glucosidase secreted by FCD3-1 was purified, and its gene was cloned and identified. The β-glucosidase belonged to glycosyl hydrolase (GH) family 3, sharing the highest identity of 94% with a GH family 3 protein from Trichoderma atroviride IMI 206040, and was designated TkBgl3A. The optimal pH and temperature of TkBgl3A were 4.5 and 65°C, respectively. The enzyme retained over 90% activity for 360 hours at pH 4.0 and 30°C, which are the usual conditions used for simultaneous saccharification and fermentation (SSF) of cellulose to ethanol. The enzyme showed significantly higher specific activity toward natural substrate cellobiose (141.4 U/mg) than toward artificial substrate p-nitrophenyl-beta-D-glucopyranoside (108.0 U/mg).
Strains of Trichoderma and Penicillium were the predominant cellulolytic fungi in subtropical and tropical forests in China. T. koningiopsis FCD3-1 was the most efficient producer of cellulase, and also produced a high level of β-glucosidase. The high specific activity toward cellobiose and stability under SSF conditions of the purified β-glucosidase from FCD3-1 indicates its potential application in SSF of cellulose to bioethanol.
Cellulose, a type of homogeneous polysaccharide that exists as units of cellobiose connected by β-1,4-glycosidic bonds, is the most abundant renewable biomass in nature, with 7.2 × 1010 tons synthesized by photosynthesis every year . Cellulose can be hydrolyzed by cellulase into D-glucose, which can be further transformed into renewable fuel and various valuable chemicals . These characteristics have stimulated studies on the degradation and utilization of cellulose. However, natural cellulose is difficult to be degraded to glucose owing to its structure and the fact that it is often found in combination with lignin and hemicellulose . Cellulose degradation and utilization are part of the global carbon cycle, with the major degradation power of cellulose in nature derived from cellulases produced by microorganisms [4–7].
Hydrolysis of cellulose to glucose requires three groups of enzymes: endo-1,4-β-glucanase (EC 188.8.131.52), cellobiohydrolase (EC 184.108.40.206), and β-glucosidase (EC 220.127.116.11). Endo-1,4-β-glucanase randomly cuts cellulose chains, yielding cellulose ends, while cellobiohydrolase attacks the cellulose chains from the reducing or non-reducing ends, yielding cellobiose, and β-glucosidase hydrolyzes cellobiose and cello-oligosaccharides to form glucose .
The mechanisms of microorganisms producing cellulases are very diverse. Among these microorganisms are the aerobic filamentous fungi, which produce an extracellular non-composite cellulase system [2, 4]. They are the main sources of commercial cellulase, with Trichoderma reesei being the most widely used .
In the 1970s, Mandels and Sternberge assessed the cellulase production of over 14,000 collected strains, and found that T. reesei produced the highest cellulase activity. Various T. reesei strains with higher yields of cellulase were then developed by mutagenesis. Of these, T. reesei Rut-C30 was one of the most characteristic strains with effective cellulase yield, and has been selected as the reference strain for many studies .
However, the hyperproducing mutant strains of T. reesei produce cellulase with all essential components, but lack the optimum amount of β-glucosidase for efficient cellulose hydrolysis, and the β-glucosidase it produces is also glucose-sensitive . As β-glucosidase catalyzes the final step of cellulose hydrolysis thus releasing glucose from the inhibitory cellobiose, it plays an important role in the biological conversion of cellulose to glucose .
Based on sequence similarities, β-glucosidases are mostly placed in either family 1 or family 3 of GH . Product inhibition and thermal inactivation of β-glucosidases constitute two major barriers to commercial development of enzymatic hydrolysis of cellulose . Therefore, the production costs of cellulosic ethanol would be significantly reduced by improving the glucose tolerance and thermostablity of β-glucosidases . In addition, β-glucosidase activities are generally measured using artificial substrates such as p-nitrophenyl-β-D-glycopyranoside (pNPG) or 4-methylumbelliferyl-β-D-glucuronide. When using these artificial substrates, most enzymes show higher catalytic activity toward artificial substrates than toward cellobiose [16–20]. However, UeBgl3A from Ustilago esculenta and a recombinant β-glucosidase from Caldicellulosiruptor saccharolyticus are exceptions, with higher catalytic activity toward cellobiose than toward artificial substrate [21, 22].
Guangxi and Yunnan Provinces are in southwest China, with the Tropic of Cancer passing through them. There is abundant sunlight, bountiful rainfall, and large areas of subtropical and tropical forests. These forests are well preserved because of low levels of human activity, and soils are mainly weakly acidic. These conditions provide abundant biomass and microbial resources. Different climatic types and natural vegetation characteristics have led to the evolution of diverse cellulolytic organisms capable of degrading cellulose .
In this work, novel fungal strains producing high cellulase activity were isolated, screened, and identified from soil samples systematically collected from virgin forests in subtropical and tropical climates in Guangxi and Yunnan Provinces in China, with a predominance of Trichoderma and Penicillium. Trichoderma koningiopsis FCD3-1 was the most efficient producer of cellulases, and also had extracellular β-glucosidase activity at levels approximately 17 times higher than that of T. reesei Rut-C30. One β-glucosidase purified from FCD3-1 was stable under the usual conditions of simultaneous saccharification and fermentation (SSF), and showed high activity toward the natural substrate cellobiose, indicating its potential application for saccharification of cellulosic biomass for bioethanol production.
Isolation, screening, and identification of fungal strains with higher cellulase activity from forests in subtropical and tropical regions of China
In total, 330 soil and rotten wood samples were collected from sampling sites in virgin forests in 8 natural reserves in Guangxi and Yunnan Provinces, China (Table 1). From these 330 samples, 305 fungal strains were isolated, using agar medium plates containing Whatman No. 1 filter paper as the sole carbon source. Of the 305 fungal strains, 231 could be cultured by shake-flask cultivation, using liquid medium containing microcrystalline cellulose (Avicel PH101) as the sole carbon source. Of these 231 fungal strains (Table 2), 51 produced Avicelase activity of greater than 0.05 U/ml 3 days after cultivation under these conditions (Table 2).
The optimum pH and temperature (Table 3) of the cellulases produced by the 51 fungal strains were measured against Avicel. Following this, the Avicelase activity of cellulases produced by the 51 strains at peak enzyme production was measured under conditions of the respective optimum pH and temperature of the crude enzyme. Of the 51 strains, 31 produced Avicelase activity of greater than 0.2 U/ml (Tables 2 and 3), compared with 0.39 U/ml for T. reesei Rut-C30 under the same cultivation conditions. Two fungal strains, FCD3-1 and BM48-3, had Avicelase activity of 0.37 and 0.36 U/ml, respectively, which did not differ statistically from that of Rut-C30. Strains HP35-3 and HP29-3 also produced relatively high Avicelase activity at 0.34 and 0.33 U/ml, respectively.
The 31 strains were further identified by sequence analysis of internal transcribed spacer (ITS) and observation of morphological characteristics. Thus, 18 strains were identified as belonging to genus Trichoderma and 13 to Penicillium (Table 3).
Notably, 9 of the 31 strains were isolated from Huaping National Nature Reserve in Guangxi, which was the sample location that gave the highest number of strains with Avicelase activity of greater than 0.2 U/ml, and 5 of these 9 strains were identified as Trichoderma harzianum (Table 3). Of the 31 strains, 8 were isolated from Xishuangbanna National Nature Reserve in Yunnan, which was the sample location with the second highest number of strains producing Avicelase activity greater than 0.2 U/ml, and 3 of these 8 strains were identified as T. harzianum (Table 3).
Enzymatic characterization of cellulases from fungal strains producing higher Avicelase activity
The cellulase produced by the four highest-producing strains, FCD3-1, BM48-3, HP35-3, and HP29-3, were further measured against substrates of filter paper, carboxylmethylcellulose sodium salt (CMC-Na), and pNPG, in addition to Avicel. Of the four strains tested, T. koningiopsis strain FCD3-1 had the most efficient enzyme activity toward Avicel, filter paper, and CMC-Na (Figure 1). Interestingly, strain FCD3-1 produced a high level of β-glucosidase activity at 1.18 U/ml, which was approximately 17 times higher than that of Rut-C30 (Figure 1).
The ITS of strain FCD3-1 was amplified by PCR and sequenced (GenBank accession number KJ619600). Using the BLAST algorithm, the ITS sequence of FCD3-1 was compared with sequences in the GenBank database, and the best BLAST match was found to be the ITS of T. koningiopsis strain T-1 (GenBank accession number KC884758) at 100% identity. The morphological identification of FCD3-1 (see Additional file 1) and the molecular identification of the strain were consistent. Therefore, fungal strain FCD3-1 was identified as T. koningiopsis.
Glucose was the sole hydrolytic product of Avicel hydrolyzed by the crude cellulase produced by strain FCD3-1 (Figure 2). In comparison, glucose and cellobiose were the major hydrolytic products of Avicel hydrolyzed by the crude cellulases produced by the fungal strains BM48-3, HP35-3, and HP29-3, and the reference strain Rut-C30 (Figure 2; see Additional file 2).
Purification and characterization of one β-glucosidase secreted by T. koningiopsis FCD3-1
To determine the β-glucosidase responsible for the high level of β-glucosidase activity by T. koningiopsis FCD3-1, the extracellular proteins secreted by this strain were subjected to protein purification. One β-glucosidase secreted by FCD3-1 was purified by using ammonium sulfate precipitation, followed by hydrophobic interaction chromatography and ion exchange chromatography. The homogeneity of the purified enzyme was confirmed by non-denaturing polyacrylamide gel electrophoresis (PAGE) analysis, which showed a single protein band after staining with Coomassie blue (Figure 3A), coincident with a β-glucosidase activity band detected using esculin as the substrate (Figure 3A). A single protein band was also detected after electrophoresis under denaturing conditions (SDS-PAGE), with molecular mass of about 100 kDa (Figure 3A).
The 100-kDa band on SDS-PAGE was excised, subjected to trypsin digestion, and analyzed by liquid chromatography–tandem mass spectrometry (LC/MS/MS). Mascot (A software from Matrix Science, London, UK) search showed that eight sequences in a GH family 3 protein (GenBank accession number EHK21862) from Trichoderma virens Gv29-8 matched (Table 4). The gene encoding the β-glucosidase was amplified from the genomic DNA of strain FCD3-1 and sequenced. This gene, with four introns, encoded a polypeptide with 866 amino acids, which contained the eight matched peptide sequences in the GH family 3 protein from T. virens Gv29-8, confirming that the cloned gene encoded the purified β-glucosidase.
The encoded polypeptide by the cloned gene shared its highest identity (94%) with a GH3 protein (GenBank accession number EHK46786) from Trichoderma atroviride IMI 206040. Simple Modular Architecture Research Tool (SMART) analysis showed that the polypeptide composed of 4 domains (aa residues of 1 to 21 belonged to the signal peptide, aa 88 to 313 belonged to the GH3 N-terminal domain, aa 413 to 652 belonged to the GH3 C-terminal domain, and aa 788 to 856 was described as a fibronectin type III-like domain), indicating that the polypeptide was a GH3 β-glucosidase. Thus, the enzyme was named TkBgl3A.
Characteristics of β-glucosidase TkBgl3A from T. koningiopsis FCD3-1
The purified β-glucosidase TkBgl3A was most active toward pNPG at pH 4.5 (Figure 3B), and this enzyme was stable in the pH range of 3.0 to 9.5 (Figure 3B). The β-glucosidase activity was optimal at 65°C (Figure 3C); and after 1 hour of incubation, TkBgl3A retained more than 50% activity at 60°C (Figure 3C). The enzyme retained more than 90% activity for 360 hours at pH 4.0 and 30°C, which are the usual conditions for SSF of cellulose to produce ethanol (Figure 3D). The Michaelis–Menten constant (Km) and maximal velocity (Vmax) for hydrolysis of pNPG were 1.21 mM and 314 μmol/min/mg, respectively. The catalytic constant (Kcat) was 523/s. The enzyme was competitively inhibited by glucose, with an inhibition constant (Ki) value of 3.29 mM.
Specific activity of TkBgl3A toward cellobiose was 141.4 U/mg, which was significantly higher than that toward artificial substrate pNPG (108 U/mg) (Table 5). It is noteworthy that high-performance liquid chromatography (HPLC) analysis showed that the purified enzyme TkBgl3A could also hydrolyze cello-oligosaccharides (cellotriose, cellotetraose, cellopentaose, cellohexaose) and mainly produced glucose with high efficiency (see Additional file 3).
Metal ions had no significant effect on TkBgl3A when using pNPG as substrate, whereas the surfactant SDS significantly inhibited the activity of TkBgl3A. The reaction mix containing 0.25% SDS retained 13.5% of enzyme activity relative to the untreated sample (see Additional file 4).
We carried out large-scale isolation and screening of secretomes with Avicelase activity from fungal strains isolated from tropical and subtropical forests in China. Kaisa et al. isolated approximately 600 mesophilic lignocellulolytic fungi from soil samples and decaying plant materials, and found that the majority of these strains were Trichoderma spp., with two Penicillium strains also found to be good cellulase producers . The predominance of Trichoderma and Penicillium strains in hydrolyzing cellulose in Guangxi and Yunnan Provinces indicated that fungal strains from these two genera may play important roles in cellulose degradation and therefore in carbon recycling in tropical and subtropical forest soils in China.
Several experimental studies have concluded that effect of the physical properties of cellulose on cellulase binding and substrate accessibility by the enzyme is responsible for the rate of hydrolysis of lignocellulose by cellulase. Physical properties of cellulose include crystallinity, degree of polymerization, and accessible surface area [26–28]. Crystalline cellulose is the main form of cellulose in most plant cell walls, and crystallinity is a key factor affecting hydrolysis of cellulose because the glycosidic bonds in crystalline regions are difficult to be hydrolyzed compared with those in amorphous regions [27, 28]. Avicel is one form of crystalline cellulose with a high degree of crystallinity, and it is very difficult to be hydrolyzed .
Desrchers et al. reported a fungal strain that can produce an enzyme complex with high carboxymethycellulase and cellobiase activity but low Avicelase activity (0.02 U/ml) . Whereas the Avicelase activity produced by an actinomycete strain isolated by Li et al. was only 0.008 U/ml , Douglas et al. isolated a fungal strain with higher Avicelase activity (0.12 to 0.17 U/ml) . In the present study, 17% of the strains isolated using Whatman No. 1 filter paper as the sole carbon source had significant Avicelase activity (greater than 0.05 U/ml), and 10% of the strains had Avicelase activity greater than 0.2 U/ml, which is a high proportion of strains producing high levels of Avicelase activity.
T. koningiopsis FCD3-1 was the best-performing isolate, producing a similar level of Avicelase activity to T. reesei Rut-C30. Interestingly, strain FCD3-1 produced about 17 times as much β-glucosidase as T. reesei Rut-C30. Although T. reesei has been extensively studied and reported to produce a highly active extracellular cellulase system, most of its strains produce β-glucosidase at low levels compared with other cellulolytic fungi such as Aspergillus species, because most of its β-glucosidase is bound to its cell wall. Exogenous supplementation of β-glucosidase in T. reesei cellulase preparations has been used as an alternative strategy to increase the rate and extent of cellulose hydrolysis . Nevertheless, wild-type Trichoderma strains producing significant quantities of β-glucosidase have been reported. Kov́acs et al. screened more than 150 wild-type Trichoderma strains isolated from 30 countries based on the overall cellulase production (filter paper activity). Trichoderma atroviride TUB F-1505 was the best wild-type extracellular cellulase producer, producing significant quantities of β-glucosidase, but it had 10 to 25% less extracellular cellulase activity than T. reesei Rut-C30. TUB F-1505 was isolated from a tropical rain forest soil sample collected near Manaus, close to the Rio Negro, Brazil . Cellulase production by six species of Trichoderma cultivated on medicinal plant processings was comparatively evaluated, and Trichoderma citrinoviride MTCC No. 2418 was found to be the most efficient producer of cellulases and also produced a high level of β-glucosidase. T. citrinoviride produced five times more β-glucosidase than produced by T. reesei MTCC No. 164 . However, purification of β-glucosidase and cloning of the gene encoding the enzyme from T. atroviride and T. citrinoviride were not reported.
In the present study, the purified β-glucosidase TkBgl3A was most active at pH 4.5. Most fungal β-glucosidases have an optimal pH range of 3.5 to 6.5 . The optimal pH for the β-glucosidase was similar to the commercial β-glucosidase Novozyme S188 . TkBgl3A was stable at pH 3.0 to 9.5, and its activity did not vary substantially within this range; its pH stability was notably better than that of previously reported β-glucosidases . TkBgl3A had an optimum temperature of 65°C; in general, the optimum temperatures for fungal β-glucosidases are in the range of 50 to 70°C . In addition, TkBgl3A displayed good thermostability, similar to Bgl3A from P. funiculosum NCL1 .
Most microbial β-glucosidases have glucose Ki ranging from as low as 0.5 mM to 100 mM or more . Although glucose acted as a competitive inhibitor of pNPG hydrolysis by TkBgl3A, with a Ki of 3.29 mM, the glucose inhibition of β-glucosidase could be overcome by an SSF process in which glucose is rapidly converted to ethanol by yeast . The purified TkBgl3A was stable under SSF conditions, which facilitates the use of the enzyme in SSF.
TkBgl3A exhibited a similar affinity toward the synthetic substrate pNPG (1.21 mM) as that of the commercially available β-glucosidase Novozyme S188 (1.03 mM). However, the specific activity of TkBgl3A for cellobiose was approximately four times higher than that of Novozyme S188 . The higher activity of TkBgl3A toward cellobiose and cello-oligosaccharides again supports the potential application of the enzyme in the final step of cellulose saccharification [40, 41].
Large-scale isolation and screening of secretomes with Avicelase activity for fungal strains from tropical and subtropical forests in China led to the isolation and identification of 31 fungal strains producing Avicelase activity of greater than 0.2 U/ml in liquid cultivation. All 31 fungal strains were identified as Trichoderma or Penicillium species, indicating the predominance of these genera in cellulose degradation in subtropical and tropical forests in China. T. koningiopsis FCD3-1 was the most efficient producer of cellulases and also had a high level of extracellular β-glucosidase. The purified β-glucosidase TkBgl3A from FCD3-1 showed high specific activity toward natural substrate cellobiose and good stability under SSF conditions, indicating its potential application in SSF of cellulose to produce bioethanol.
Collection of soil and rotten wood samples
Soil and rotten wood in virgin forests were collected in 2007 and 2008 (July to October) from eight natural reserves in both tropical and subtropical areas in different biotopes in China. The sampling route was set at intervals of 100 meters, with one sample collected per interval. Surface soil at sampling sites was removed, and the soil at a depth of 50 to 200 mm was collected, with about 50 g per sample.
Cultivation and isolation of fungal strains
The suspensions of soil samples were serially diluted. Three dilutions of each sample were chosen, and 0.1 ml of each sample was spread onto the surface of Medium P agar plates containing Whatman No. 1 filter paper as the sole carbon source. The composition of Medium P was as follows: 2 g/l NaNO3, 0.5 g/l MgSO4 · 7H2O, 0.5 g/l NaCl, 0.01 g/l FeSO4 · 7H2O, 1.0 g/l KH2PO4, 0.4 g/l yeast, extract, 10 g/l Whatman No. 1 filter paper, and 15.0 g/l agar, at pH 5.0. After incubation at 28°C for 3 to 5 days, the single colonies of fungi were picked using an inoculating hook and inoculated onto a new Medium P agar plate.
Shake-flask cultivation of fungal strains
The medium used for the shake-flask cultivation was named Medium A, and contained microcrystalline cellulose powder (Avicel PH101, Sigma, St. Louis, MO, USA) as the sole carbon source. The composition of Medium A was as follows: 10 g/l Avicel (Sigma, USA), 2 g/l KH2PO4, 1.4 g/l (NH4)2SO4, 0.3 g/l MgSO4 · 7H2O, 0.3 g/l CaCl2, 1 g/l CaCO3, 1 g/l peptone (Difco, USA), 1 g/l Tween-80, 0.3 g/l urea, 10 mg/l FeSO4 · 7H2O, 3.2 mg/l MnSO4, 6.9 mg/l ZnSO4 · 7H2O, and 0.4 mg/l CoCl2 · 6H2O. All materials were supplied by Sinopharm Chemical Reagent Co., Ltd., Shanghai, China unless otherwise stated. The pH was adjusted to 5.0 before sterilization.
Shake-flask cultivation was carried out in 500 ml cotton-plugged Erlenmeyer flasks containing 100 ml of Medium A. After autoclaving at 121°C for 20 minutes, the flasks were inoculated by removing the mycelia from a single colony of fungi on a Medium P agar plate. Flask cultivation was performed at 28°C on a rotary shaker at 200 rpm. After 3 days of cultivation, samples were removed and centrifuged at 4°C and 9,000 × g (Centrifuge Model:himac CF 16RX, Rotor type: T15A29, Hitachi Koki Co., Ltd., Tokyo, Japan) for 10 minutes, and the clear supernatants were removed and used for enzyme assays.
Avicelase activity produced by isolated fungal strains was determined using the 3,5-dinitrosalicylic acid (DNS) method . The reaction system was prepared by mixing 50 μl of crude enzyme (appropriately diluted) with 450 μl of 1.1% (w/v) Avicel. The buffer used for resuspending the substrate was 0.1 M sodium citrate buffer at appropriate pH. The mixtures were incubated at 50°C for 1 hour, and then the reactions were stopped by adding 1 ml of DNS reagent. All mixtures were heated in boiling water for 5 minutes for color development. Subsequently, 200 μl of each sample was transferred to a 96-well microplate, and absorbance was measured at 540 nm. 1 U of Avicelase activity was defined as the amount of enzyme that released 1 μmol of reducing sugars equivalent to glucose per minute during the reaction.
To determine the optimal pH of the crude enzyme produced by isolated fungal strains, 50 μl of supernatant from flask cultivation were incubated with 450 μl of Avicel (1.1% w/v) at 50°C with a range of different pH levels (3.0 to 6.0 at intervals of 0.5). To determine the optimal temperature of crude enzyme produced by isolated fungal strains, 50 μl of supernatant from flask cultivation were incubated with 450 μl of Avicel (1.1%, w/v) at the optimal pH and a temperature range of 30 to 65°C with an interval of 5°C. The maximum Avicelase activity obtained at different pH or temperature levels was considered to be 100%.
Hydrolysis products of Avicel by the crude enzymes produced by fungal strains FCD3-1, BM48-3, HP35-3, HP29-3, and Rut-C30 were analyzed by HPLC. After 5 days of shaking cultivation at 28°C in 500 ml cotton-plugged Erlenmeyer flasks containing 100 ml of Medium A, the cultures of fungal strains were centrifuged at 3,000 × g (Centrifuge Model:himac CF 16RX, Rotor type: T5SS31, Hitachi Koki Co., Ltd., Tokyo, Japan) for 10 minutes, and the filter paper activity of each collected supernatant was measured. Each enzyme solution was respectively added to 100 mg of Avicel in a test tube at a dose of 10 U/g Avicel and the total volume of the mixture was adjusted to 5 ml with 0.1 M citrate–phosphate buffer (CPB) at pH 5.0. The hydrolysis was performed in a water bath shaker with shaking at 200 rpm at 45°C. Hydrolysis was stopped after reactions were conducted for 0, 2, 4, 8, 16 and 32 hours by boiling for 10 minutes. After the cooled sample was filtered through a 0.22 μm filter (Nylon 66 Syringe Filter, Jinteng, China), the collected filtrate was analyzed with HPLC (LC-10A VP Plus, Shimadzu Kyoto, Japan) equipped with a refractive index detector (Shimadzu). Sugars were separated using a Hypersil NH2 column (Dalian Elite Analytical Instruments Co., Ltd., China) at 30°C, with 75% acetonitrile as eluent, at a flow rate of 1.0 ml/min.
Identification of isolated fungal strains
Cell mycelia of fungal isolates were obtained from pure culture in Yeast Mold (YM) broth (10 g/l glucose, 3 g/l yeast extract, 3 g/l malt extract, and 5 g/l peptone). Genomic DNA of the isolates was extracted and used as PCR template. PCR was performed to amplify the ITS (approximately 600 bp, including ITS1, the 5.8 S rRNA gene, and ITS2 region) of fungal ribosomal DNA with primers ITS1 5′-TCCGTAGGTGAACCTGCGG-3′ and ITS4 5′-TCCTCCGCTTATTGATATGC-3′. The 25 μl mixtures were composed of 1 μl of template DNA, 12.5 μl of 2 × Taq PCR Master Mix (purchased from Tiangen Biotech, Beijing, China; containing 0.5 U Taq DNA polymerase, 500 μM of each dNTP, 20 mM Tris–HCl pH 8.3, 100 mM KCl, 3 mM MgCl2, and bromophenol blue), 1 μl of each primer (10 μM) and 9.5 μl of double-distilled H2O. The PCR procedure was as follows: primary denaturation for 5 minutes at 94°C, followed by 30 cycles of denaturation at 94°C for 30 seconds, annealing at 55°C for 30 seconds, and extension at 72°C for 60 seconds, then an additional reaction for 10 minutes at 72°C. The PCR products were separated in 0.8% agarose gels to confirm purity, quantity, and size, and sent to Sangon Biotech (Shanghai) Co., Ltd. (China) for sequencing. The ITS sequences were compared with other ITS sequences in GenBank using the BlastN program (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Media and culture conditions for cultivation of T. koningiopsis FCD3-1
The media used for cultivation of T. koningiopsis FCD3-1 were potato dextrose agar (PDA), Avicel agar medium, or fungal cultivation medium. The Avicel agar medium contained: 3 g Avicel PH-101 (Sigma), 2.0 g KH2PO4, 1.4 g (NH4)2SO4, 0.3 g MgSO4 · 7H2O, 0.3 g CaCl2, 5 mg FeSO4 · 7H2O, 1.6 mg MnSO4, 1.7 mg ZnCl2, 2 mg CoCl2, and 20 g agar in 1 litre ultra-pure water. The fungal cultivation medium contained: 4.0 g KH2PO4, 2.8 g (NH4)2SO4, 0.6 g MgSO4 · 7H2O, 0.6 g CaCl2, 0.6 g urea, 10 mg FeSO4 · 7H2O, 3.2 mg MnSO4, 2.8 mg ZnSO4, 2.4 mg CoCl2, 20 g wheat bran, and 30 g Avicel PH-101 in 1 litre of ultra-pure water. Each medium was adjusted to pH 5.0.
To cultivate the T. koningiopsis FCD3-1 for high-level production of β-glucosidase, FCD3-1 was firstly inoculated onto PDA plates at 28°C for 72 hours. The hyphae were then incubated in Avicel agar medium at 28°C for 96 hours, and incubated in fungal cultivation medium at 28°C on a rotary shaker at 200 rpm for 1 week to assess extracellular β-glucosidase production.
Purification and characterization of one β-glucosidase secreted by T. koningiopsis FCD3-1
The culture of FCD3-1 was centrifuged for 20 minutes at 9,000 × g (Centrifuge Model:himac CF 16RX, Rotor type: T15A29, Hitachi Koki Co., Ltd., Tokyo, Japan) and 4°C after cultivation. The supernatant was carefully transferred to a sanitized container. Purification of β-glucosidase was mainly based on the methods of Parry et al. . The supernatant (crude enzyme) was subjected to 40% (NH4)2SO4 precipitation. After centrifugation, the supernatant was collected and then subjected to 60% (NH4)2SO4 precipitation overnight. The precipitated protein was dissolved in 1.5 M (NH4)2SO4 solution and loaded onto a HiPrep 16/10 Phenyl FF (high sub) column (GE Healthcare, Amersham Place, Buckinghamshire, UK). The eluted β-glucosidase was then loaded on an anion column (5 ml HiTrap™ Q FF, GE Healthcare, UK).
The purified protein sample was separated by SDS-PAGE and the protein band was excised. The band pieces were washed with 100 μl dH2O for 15 minutes and were further successively washed for 10 minutes with 100 μl CH3CN and 200 μl 50 mM NH4HCO3/CH3CN (50:50 v/v). After the band pieces were treated with 100 μl of CH3CN to dehydrate the proteins for 10 minutes, the band pieces were treated with 20 μl of 10 ng/μl modified trypsin in 20 mM NH4HCO3 for 30 minutes, then more 20 mM NH4HCO3 was added to cover the band pieces. After overnight incubation at 37°C, an equal volume of CH3CN was added. After incubation at room temperature for 30 minutes and centrifugation, the supernatant was transferred to a new clean Eppendorf tube, and then lyophilized by vacuum centrifugation and stored at -80°C.
Liquid chromatography was performed on a nano Acquity UPLC system (Waters Corporation, Milford, USA) connected to a LTQ Orbitrap XL mass spectrometer (Thermo Scientific, Bremen, Germany) equipped with an online nano-electrospray ion source (Michrom Bioresources, Auburn, USA). Peptides were resuspended with 25 μl of solvent A (5% acetonitrile and 0.1% formic acid in water), then 20 μl of peptide solution, was loaded onto a Captrap Peptide column (0.5 mm internal diameter × 2.0 mm; Michrom Bioresources) using a 20 μl/min flow rate of solvent A for 5 minutes, and then separated on a Magic C18AQ reverse phase column (100 μm internal diameter × 15 cm; Michrom Bioresources) with a linear gradient, starting from 5% solvent B (90% acetonitrile and 0.1% formic acid in water) to 45% B (that is, from 95% A to 55% A) in 70 minutes. The column flow rate was maintained at 500 nl/min and column temperature at 35°C. An electrospray voltage of 1.4 kV versus the inlet of the mass spectrometer was used.
A LTQ Orbitrap XL mass spectrometer was operated in the data-dependent mode to switch automatically between MS and MS/MS acquisition. Survey full-scan MS spectra with one microscan (350 to 1600 m/z) was acquired in the Obitrap with a mass resolution of 100,000 at 400 m/z, followed by MS/MS of the eight most intense peptide ions in the LTQ analyzer. The automatic gain control (AGC) was set to 1,000,000 ions, with maximum accumulation times of 500 ms. For MS/MS, we used an isolation window of 2 m/z, and the AGC of LTQ was set to 20 000 ions, with maximum accumulation time of 120 ms. Single charge state was rejected, and dynamic exclusion was used with two microscans using exclusion durations of 10 and 90 seconds. For MS/MS, precursor ions were activated using 35% normalized collision energy at the default activation q of 0.25 and an activation time of 30 ms. The spectra were recorded with Xcalibur (v2.0.7) software.
The mass spectra were searched using the Mascot Daemon software (v2.3.0; Matrix Science, London, UK) based on the Mascot algorithm. We searched the National Center for Biotechnology Information database (downloaded on 27 March 2013, taxonomy is Trichoderma with 39.001 entries). To reduce false-positive identification results, a decoy database containing the reverse sequences was appended to the database. The searching parameters were set up as follows: full trypsin cleavage with two missed cleavages was considered; oxidation on methionine and aceytlation of the protein N-terminus were set as variable modifications.; the peptide mass tolerance was 10 ppm; and the fragment ion tolerance was 1.0 Da. Peptides for which Mascot scores exceeded the 99.9% confidence level score (ions score ≥ 40) were accepted as correct matches.
Protein and enzyme activity assay of β-glucosidase secreted by T. koningiopsis FCD3-1
Determination of the total protein was performed according to the method of Branford . Standard curves were prepared with bovine serum albumin. For enzyme assay, each reaction was performed at 50°C and pH 5.0 (in triplicate if not otherwise mentioned). The reaction mixture was 116 μl of 0.1 M CPB (pH 5.0), 14 μl of substrate solution (25 mM pNPG) and 10 μl of enzyme solution. After addition of enzyme, the mixture was incubated for 10 minutes at 50°C and stopped with 70 μl of 0.4 M Na2CO3. The developed color was read at 410 nm on a microplate reader (BioTek, USA). One unit of β-glucosidase activity was expressed as the amount of enzyme required to release 1 μmol of p-nitrophenol (pNP) per minute.
The optimal pH of β-glucosidase activity was determined by incubating the purified enzyme at 50°C for 10 minutes in CPB with pH values in the range of 3.0 to 7.0 at intervals of 0.5. To determine the optimal temperature, the enzyme was incubated in CPB (0.1 M and pH 4.5) for 10 minutes at different temperatures in the range of 30 to 80°C at intervals of 5°C .
The pH stability of purified enzyme was assessed by incubating 500 μl of purified enzyme at 4°C for 24 hours in various buffers adjusted to different pH values, followed by checking the residual activity . To determine the thermostability of β-glucosidase activity, the purified enzyme was incubated at different temperatures (30, 35, 40, 45, 50, 55, 60, and 65°C) in the absence of substrate. After incubation for 1 hour, the residual β-glucosidase activity was determined as described above. Stability of the purified enzyme under SSF conditions (30°C and pH 4.0) was assessed by mixing 100 μl of purified enzyme with CPB pH 4.0 to a volume of 5 ml, followed by incubating the enzyme mixture at 30°C for varying times before checking the remaining activity.
The influence of metal ions on enzyme activity was determined by mixing the enzyme with CPB (pH 4.5) that contained various metal ions and reagents.
The kinetic parameters (Vmax and Km) in hydrolysis of pNPG by the purified β-glucosidase were determined under optimal conditions, and the values for Km and Vmax were estimated by applying a nonlinear curve fit using GraphPad Prism4 from GraphPad Software (San Diego, CA, USA). Catalytic constants (Kcat) and catalytic efficiency ratios (Kcat/Km) were determined from the obtained kinetic parameter values .
The activity of the purified β-glucosidase against pNPG and cello-oligosaccharides (degree of polymerization from 2 to 6) was measured in 0.1 M CPB (pH 4.5) at 65°C for 10 minutes. For the hydrolysis of the pNPG, 250 μl reaction mixtures were set up by adding 0.1 μg purified enzyme and pNPG stock solution to final concentration of 1 mM. For hydrolysis of cello-oligosaccharides, 250 μl reaction mixtures were set up by adding 0.28 μg purified enzyme and cello-oligosaccharide to final concentration of 1 mM. The amount of released glucose was determined under the same conditions as described above using a glucose oxidase kit (Maker Co., Ltd., China). One unit of β-glucosidase activity was expressed as the amount of enzyme required to release 1 μmol of pNP or glucose per minute . In addition, the hydrolysates of cello-oligosaccharides produced by the purified β-glucosidase TkBgl3A were analyzed by HPLC using a refractive index detector (RI detector k-2301; Knauer, Germany). Sugars were separated using a Hypersil NH2 column (Dalian Elite Analytical Instrument Co., Ltd., China) at 25°C, with 70% acetonitrile as eluent, at a flow rate of 1.0 mL/min.
Cloning of Tkbgl3A gene
For isolation of genomic DNA from T. konigiopsis strain FCD3-1, the fungal strain was cultivated on PDA plates for 3 days. Genomic DNA was isolated from the lyophilized mycelia using the hexadecyl trimethyl ammonium bromide (CTAB) extraction method . To amplify the TkBgl3A gene, the forward primer 5′-ATGTTGTCAAATTCAATCATCGCCGTTG-3′ and reverse primer 5′-CTAAGGAAGAACCTGGCTCAAGGGTAAA-3′ were used, which were designed based on the coding sequence (GenBank accession number ABDF02000059) of GH family 3 protein from T. virens GV29-8 using vector NTI 10 software. The extracted genomic DNA from T. konigiopsis FCD3-1 was used as template. The PCR amplification was carried out using the following condition: one cycle at 95°C for 4 minutes, followed by 30 cycles at 95°C for 30 seconds, 56°C for 30 seconds, and 72°C for 1.5 minutes, and then one cycle at 72°C for 10 minutes. The PCR product was cloned into the T-vector, and then transformed into Escherichia coli DH5α. Three of the transformants were sequenced by BGI TechSolutions Co., Ltd. (Shenzhen, China). The sequence analyses were carried out using the Lasergene DNASTAR software package, Vector NTI 10 software, NCBI-BLAST, and SMART (http://smart.embl-heidelberg.de).
Nucleotide sequence accession numbers
All the ITS sequences of the 31 fungal strains reported in this paper were deposited in the GenBank database under accession numbers KJ619589 to KJ619593 and KJ619599 to KJ619642. The β-glucosidase gene TkBgl3A is under GenBank accession number KJ623247.
Malhi Y: Carbon in the atmosphere and terrestrial biosphere in the 21st century. Philos Trans A Math Phys Eng Sci 2002, 360: 2925-2945.
Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev 2002, 66: 506-577.
Boer W, Folman LB, Summerbell RC, Boddy L: Living in a fungal world: impact of fungi on soil bacterial niche development. FEMS Microbiol Rev 2005, 29: 795-811.
Wilson DB: Microbial diversity of cellulose hydrolysis. Curr Opin Microbiol 2011, 14: 259-263.
Weimer PJ, Russell JB, Muck RE: Lessons from the cow: what the ruminant animal can teach us about consolidated bioprocessing of cellulosic biomass. Bioresour Technol 2009, 100: 5323-5331.
Distel DL, Beaudoin DJ, Morrill W: Coexistence of multiple proteobacterial endosymbionts in the gills of the wood-boring Bivalve Lyrodus pedicellatus (Bivalvia: Teredinidae). Appl Environ Microbiol 2002, 68: 6292-6299.
Ohkuma M: Termite symbiotic systems: efficient bio-recycling of lignocellulose. Appl Microbiol Biotechnol 2003, 61: 1-9.
Béguin P, Aubert JP: The biological degradation of cellulose. FEMS Microbiol Rev 1994, 13: 25-58.
Gusakov AV: Alternatives to Trichoderma reesei in biofuel production. Trends Biotechnol 2011, 29: 419-425.
Le Crom S, Schackwitz W, Pennacchio L, Magnuson JK, Culley DE, Collett JR, Martin J, Druzhinina IS, Mathis H, Monot F, Seiboth B, Cherry B, Rey M, Berka R, Kubicek CP, Baker SE, Margeot A: Tracking the roots of cellulase hyperproduction by the fungus Trichoderma reesei using massively parallel DNA sequencing. Proc Natl Acad Sci U S A 2009, 106: 16151-16156.
Singhania RR, Patel AK, Sukumaran RK, Larroche C, Pandey A: Role and significance of beta-glucosidases in the hydrolysis of cellulose for bioethanol production. Bioresour Technol 2013, 127: 500-507.
Hong J, Tamaki H, Kumagai H: Cloning and functional expression of thermostable β-glucosidase gene from Thermoascus aurantiacus . Appl Microbiol Biotechnol 2007, 73: 1331-1339.
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res 2009, 37: D233-D238.
Yoon JJ, Kim KY, Cha CJ: Purification and characterization of thermostable β-glucosidases from the brown-rot basidiomycete Fomitopsis palustris grown on microcrystalline cellulose. J Microbiol 2008, 46: 51-55.
Liu W, Hong J, Bevan DR, Zhang YH: Fast identification of thermostable beta-glucosidase mutants on cellobiose by a novel combinatorial selection/screening approach. Biotechnol Bioeng 2009, 103: 1087-1094.
Del Pozo MV, Fernández-Arrojo L, Gil-Martínez J, Montesinos A, Chernikova TN, Nechitaylo TY, Waliszek A, Tortajada M, Rojas A, Huws SA, Golyshina OV, Newbold CJ, Polaina J, Ferrer M, Golyshin PN: Microbial β-glucosidases from cow rumen metagenome enhance the saccharification of lignocellulose in combination with commercial cellulase cocktail. Biotechnol Biofuels 2012, 5: 73.
Liu D, Zhang R, Yang X, Zhang Z, Song S, Miao Y, Shen Q: Characterization of a thermostable β-glucosidase from Aspergillus fumigatus Z5, and its functional expression in Pichia pastoris X33. Microb Cell Fact 2012, 11: 25.
Uchima CA, Tokuda G, Watanabe H, Kitamoto K, Arioka M: Heterologous expression and characterization of a glucose-stimulated β-glucosidase from the termite Neotermes koshunensis in Aspergillus oryzae . Appl Microbiol Biotechnol 2011, 89: 1761-1771.
Feng Y, Duan CJ, Liu L, Tang JL, Feng JX: Properties of a metagenome-derived beta-glucosidase from the contents of rabbit cecum. Biosci Biotechnol Biochem 2009, 73: 1470-1473.
Beloqui A, Nechitaylo TY, López-Cortés N, Ghazi A, Guazzaroni ME, Polaina J, Strittmatter AW, Reva O, Waliczek A, Yakimov MM, Golyshina OV, Ferrer M, Golyshin PN: Diversity of glycosyl hydrolases from cellulose-depleting communities enriched from casts of two earthworm species. Appl Environ Microbiol 2010, 76: 5934-5946.
Nakajima M, Yamashita T, Takahashi M, Nakano Y, Takeda T: Identification, cloning, and characterization of β-glucosidase from Ustilago esculenta . Appl Microbiol Biotechnol 2012, 93: 1989-1998.
Hong MR, Kim YS, Park CS, Lee JK, Kim YS, Oh DK: Characterization of a recombinant beta-glucosidase from the thermophilic bacterium Caldicellulosiruptor saccharolyticus . J Biosci Bioeng 2009, 108: 36-40.
Merino ST, Cherry J: Progress and challenges in enzyme development for biomass utilization. Adv Biochem Eng Biotechnol 2007, 108: 95-120.
Montenecourt BS, Eveleigh DE: Selective screening methods for the isolation of high yielding cellulase mutants of Trichoderma reesei . In Hydrolysis of Cellulose: Mechanisms of Enzymatic and Acid Catalysis, Advances in Chemistry Series, Volume 181. Edited by: Brown R, Jurasek L. Washington, DC: American Chemical Society; 1979:289-301.
Marjamaa K, Toth K, Bromann PA, Szakacs G, Kruus K: Novel Penicillium cellulases for total hydrolysis of lignocellulosics. Enzyme Microb Technol 2013, 52: 358-369.
Turon X, Rojas OJ, Deinhammer RS: Enzymatic kinetics of cellulose hydrolysis: a QCM-D study. Langmuir 2008, 24: 3880-3887.
Ahola S, Turon X, Osterberg M, Laine J, Rojas OJ: Enzymatic hydrolysis of native cellulose nanofibrils and other cellulose model films: effect of surface structure. Langmuir 2008, 24: 11592-11599.
Chang VS, Holtzapple MT: Fundamental factors affecting biomass enzymatic reactivity. Appl Biochem Biotechnol 2000, 84–86: 5-37.
Zhang YH, Lynd LR: Toward an aggregated understanding of enzymatic hydrolysis of cellulose: noncomplexed cellulase systems. Biotechnol Bioeng 2004, 88(7):797-824.
Desrochers M, Jurasek L, Paice MG: Production of cellulase, beta-glucosidase, and xylanase by schizophyllum commune grown on a cellulose-peptone medium. Dev Ind Microbiol 1981, 22: 675-684.
Li X, Gao P: Isolation and partial characterization of cellulose-degrading strain of Streptomyces sp. LX from soil. Lett Appl Microbiol 1996, 22: 209-213.
Mountfort DO, Asher RA: Production and regulation of cellulase by two strains of the rumen anaerobic fungus Neocallimastix frontalis . Appl Environ Microbiol 1985, 49: 1314-1322.
Tiwari P, Misra BN, Sangwan NS: β-Glucosidases from the fungus Trichoderma : an efficient cellulase machinery in biotechnological applications. Biomed Res In 2013, 2013: 203735.
Kov́acs K, Megyeri L, Szakacsa G, Kubicekc CP, Galbeb M, Zacchi G: Trichoderma atroviride mutants with enhanced production of cellulase and β-glucosidase on pretreated willow. Enzyme Microb Technol 2008, 43: 48-55.
Chandra M, Kalra A, Sharma PK, Sangwan RS: Cellulase production by six Trichoderma spp. fermented on medicinal plant processings. J Ind Microbiol Biotechnol 2009, 36: 605-609.
Park AR, Hong JH, Kim JJ, Yoon JJ: Biochemical characterization of an extracellular β-glucosidase from the fungus, Penicillium italicum , isolated from rotten citrus peel. Mycobiology 2012, 40: 173-180.
Chauve M, Mathis H, Huc D, Casanave D, Monot F, Lopes Ferreira N: Comparative kinetic analysis of two fungal beta-glucosidases. Biotechnol Biofuels 2010, 3: 3.
Ramani G, Meera B, Vanitha C, Rao M, Gunasekaran P: Production, purification, and characterization of a β-glucosidase of Penicillium funiculosum NCL1. Appl Biochem Biotechnol 2012, 167: 959-972.
Ohgren K, Bura R, Lesnicki , Saddler J, Zacchi G: A comparison between simultaneous saccharification and fermentation and separate hydrolysis and fermentation using steam-pretreated corn stover. Process Biochem 2007, 42: 834-839.
Bhatia Y, Mishra S, Bisaria VS: Purification and characterization of recombinant Escherichia coli -expressed Pichia etchellsii beta-glucosidase II with high hydrolytic activity on sophorose. Appl Microbiol Biotechnol 2005, 66: 527-535.
Chen H, Hayn M, Esterbauer H: Purification and characterization of two extracellular beta-glucosidases from Trichoderma reesei . Biochim Biophys Acta 1992, 1121: 54-60.
Miller GL: Use of Dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem 1959, 31: 426-428.
Ghose TK: Measurement of cellulase activities. Pure Appl Chem 1987, 59: 257-268.
Berghem LE, Pettersson LG: The mechanism of enzymatic cellulose degradation. Isolation and some properties of a beta-glucosidase from Trichoderma viride . Eur J Biochem 1974, 46: 295-305.
Parry NJ, Beever DE, Owen E, Vandenberghe I, Van Beeumen J, Bhat MK: Biochemical characterization and mechanism of action of a thermostable beta-glucosidase purified from Thermoascus aurantiacus . Biochem J 2001, 353: 117-127.
Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 1976, 72: 248-254.
Michlmayr H, Schümann C, da Silva NM, Kulbe KD, del Hierro AM: Isolation and basic characterization of a β-glucosidase from a strain of Lactobacillus brevis isolated from a malolactic starter culture. J Appl Microbiol 2010, 108: 550-559.
Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard RW: Ribosomal DNA spacer-length polymorphisms in barley: mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci 1984, 81: 8014-8018.
This work was supported by grants from the Guangxi Natural Science Foundation (grant number 2012GXNSFGA060005), the Bagui Scholar Program of Guangxi (grant number 2011A001), Guangxi Program for Science and Technology Development (grant number 11107008-5) and the National Natural Science Foundation of China (grant number 30960013).
The authors declare that they have no competing interests.
ZZ collected soil samples, isolated the fungal strains, assayed the fungal enzyme production, identified the fungal isolates, and was involved in drafting the manuscript. JLL and JYL purified and characterized the β-glucosidase TkBgl3A, cloned the encoding gene, and participated in drafting the manuscript. CJD participated in data analysis. QSM participated in conceiving and designing the study. JXF conceived, designed the study, and critically revised and finally approved the manuscript. All authors read and approved the final manuscript.
Zheng Zhang, Jun-Liang Liu, Jian-Yi Lan contributed equally to this work.
Electronic supplementary material
Additional file 1: T. koningiopsis FCD3-1.(PDF 236 KB)
Additional file 2:(PDF 108 KB)
The hydrolysates of cello-oligosaccharides by the purified β-glucosidase TkBgl3A at different time intervals.
Additional file 3:(PDF 84 KB)
Additional file 4: T. koningiopsis FCD3-1.(PDF 126 KB)
About this article
Cite this article
Zhang, Z., Liu, JL., Lan, JY. et al. Predominance of Trichoderma and Penicillium in cellulolytic aerobic filamentous fungi from subtropical and tropical forests in China, and their use in finding highly efficient β-glucosidase. Biotechnol Biofuels 7, 107 (2014). https://doi.org/10.1186/1754-6834-7-107