- Open Access
Differential β-glucosidase expression as a function of carbon source availability in Talaromyces amestolkiae: a genomic and proteomic approach
Biotechnology for Biofuels volume 10, Article number: 161 (2017)
Genomic and proteomic analysis are potent tools for metabolic characterization of microorganisms. Although cellulose usually triggers cellulase production in cellulolytic fungi, the secretion of the different enzymes involved in polymer conversion is subjected to different factors, depending on growth conditions. These enzymes are key factors in biomass exploitation for second generation bioethanol production. Although highly effective commercial cocktails are available, they are usually deficient for β-glucosidase activity, and genera like Penicillium and Talaromyces are being explored for its production.
This article presents the description of Talaromyces amestolkiae as a cellulase-producer fungus that secretes high levels of β-glucosidase. β-1,4-endoglucanase, exoglucanase, and β-glucosidase activities were quantified in the presence of different carbon sources. Although the two first activities were only induced with cellulosic substrates, β-glucosidase levels were similar in all carbon sources tested. Sequencing and analysis of the genome of this fungus revealed multiple genes encoding β-glucosidases. Extracellular proteome analysis showed different induction patterns. In all conditions assayed, glycosyl hydrolases were the most abundant proteins in the supernatants, albeit the ratio of the diverse enzymes from this family depended on the carbon source. At least two different β-glucosidases have been identified in this work: one is induced by cellulose and the other one is carbon source-independent. The crudes induced by Avicel and glucose were independently used as supplements for saccharification of slurry from acid-catalyzed steam-exploded wheat straw, obtaining the highest yields of fermentable glucose using crudes induced by cellulose.
The genome of T. amestolkiae contains several genes encoding β-glucosidases and the fungus secretes high levels of this activity, regardless of the carbon source availability, although its production is repressed by glucose. Two main different β-glucosidases have been identified from proteomic shotgun analysis. One of them is produced under different carbon sources, while the other is induced in cellulosic substrates and is a good supplement to Celluclast in saccharification of pretreated wheat straw.
Second generation bioethanol represents an efficient alternative to conventional energy supply, involving the exploitation of renewable sources, usually disposable wastes from other industrial or agricultural activities. In this process, polysaccharides from lignocellulosic biomass are hydrolyzed to glucose and xylose units that are further converted to ethanol . When raw materials with hardly available polysaccharides are used, a pretreatment is needed in order to facilitate its accessibility. Steam-explosion is maybe the most commonly applied method to perform this pretreatment, generating a biomass slurry that usually contains by-products which adversely affect downstream steps like enzymatic hydrolysis or ethanol fermentation [2, 3]. Although all stages are relevant for the process, the improvement of the composition and dosage of the hydrolytic enzymatic cocktails remains a subject of great interest to industry. An efficient hydrolysis of carbohydrates, especially cellulose, the major polysaccharide from plant biomass, is indispensable to recover fermentable sugars from feedstock. Cellulose is a homopolysaccharide of β-1,4-linked d-glucose residues whose enzymatic conversion into free glucose requires the coordinated work of the cellulase system, consisting of three different kinds of hydrolases: (1) cellobiohydrolases (CBH) named also as exoglucanases, which cleave cellobiose units from the chain ends, (2) β-1,4-endoglucanases (EG), which hydrolyze the polymers internally, reducing their degree of polymerization and (3) β-glucosidases (BGLU), which convert cellobiose and soluble cellodextrins into glucose .
The main industrial producers of cellulases are Trichoderma reesei and Aspergillus niger [5, 6], but the enzymatic crudes released by these microorganisms, in particular those from Trichoderma sp., are deficient in β-glucosidase , and the saccharification cocktails must be supplemented with this activity to increase the efficiency of cellulose degradation.
The current approach for production of lignocellulosic bioethanol relies on combining different cellulase cocktails, which is very effective but involves a significant increase of the production costs. Then, many efforts are being devoted to discover microorganisms that secrete high amounts of glycosidases and to search for novel robust enzymes, which efficiently transform cell-wall polysaccharides, even in the presence of the undesirable by-products found in the process streams. In this sense, Penicillium sp. and its perfect states (Talaromyces or Eupenicillium) have attracted much attention in the last years for their high cellulase activity [8, 9]. In addition, these fungi could have special interest for cellulose and hemicelluloses transformation in the context of bioethanol production or other industrial applications, since they also produce high levels of xylanases [10,11,12].
This work reports the characterization of T. amestolkiae as cellulase producer, and the sequencing and assembly of its genome. The cellulolytic activities secreted upon fungal growth in several carbon sources were quantified, and the pool of extracellular proteins in each one of the secretomes was analyzed by massive peptide analysis. The data obtained were comprehensively related to the number of genes codifying for members of different glycosyl hydrolases families (GH) in T. amestolkiae genome, showing that β-glucosidases are highly represented from DNA to proteins secreted in this fungus. The role of these enzymes in saccharification of wheat straw acid slurry is discussed.
Fungal strain and culture media
The fungus was isolated from cereal wastes and deposited in the IJFM culture collection at “Centro de Investigaciones Biológicas” (Madrid, Spain), with the reference A795. This isolate was identified on the basis of molecular and phenotypic analyses (see Additional file 1), and hereinafter, it will be denominated as T. amestolkiae CIB. Fungal cultures were grown on 2% agar-malt Petri dishes at 28 °C and incubated for 7 days to obtain spore suspensions. One agar plug of about 1 cm2 was cut from actively growing mycelium and deposited in 15 mL Falcon tubes with 5 mL of 1% NaCl with 0.1% Tween 80, shaken, and 200 μL withdrawn to inoculate 250 mL Erlenmeyer flasks with 50 mL of CSS medium (28 °C, 250 rpm for 7 days). CSS medium (pH 5.6) contained (L−1): 40 g glucose, 0.4 g FeSO4 × 7H2O, 9 g (NH4)2SO4, 4 g K2HPO4, 26.3 g corn steep solids, 7 g CaCO3, and 2.8 mL soybean oil. Two mL of these cultures were used as inocula for the production of cellulase activities in Mandels medium . Its components are (L−1): 2.0 g KH2PO4, 1.3 g (NH4)2SO4, 0.3 g urea, 0.3 g MgSO4·7H2O, 0.3 g CaCl2, 5 mg FeSO4·7H2O, 1.6 mg MnSO4·H2O, 1.4 mg ZnSO4·7H2O, and 1 g Bacto Peptone. The medium was supplemented with 1% Avicel, glucose, acid slurry from wheat straw (acid-catalyzed steam-exploded wheat straw provided by Abengoa), maltose, fructose, glycerol, xylose or cellobiose, or 2% beechwood xylan (reported as the best concentration for xylanases’ production ) as carbon sources. In all cases, the cultures were carried out in 250 mL Erlenmeyer flasks with 50 mL of culture medium and incubated at 28 °C and 250 rpm. Samples were taken periodically from three replicate flasks and the mycelium separated from the culture liquid by centrifugation at 13,000 × g and 4 °C for 5 min.
Enzyme and protein assays and other determinations
Proteins were quantified according to the Bradford method using the Bio-Rad reagent and bovine serum albumin as the standard. Avicelase (total microcrystalline-cellulose hydrolyzing activity, as indicative of exocellulase activity) and β-1,4-endoglucanase activities were measured by determining the release of reducing sugars by the Somogyi–Nelson method . The standard enzymatic assays were performed in 50 mM sodium citrate buffer, pH 5.0, containing appropriately diluted crudes and 1% Avicel (Merck) or 2% low viscosity carboxymethylcellulose (CMC, Sigma), as substrates for Avicelase and β-1,4-endoglucanase activities, respectively. Linearity of the enzymatic assays was checked by performing the reaction at two different incubation times (5 and 10 min). β-glucosidase activity was assayed spectrophotometrically following p-nitrophenol (pNP) release (ε 410 = 15,200 M−1 cm−1) from p-nitrophenyl-β-d-glucopyranoside (Sigma-Aldrich), in 100 mM sodium citrate buffer pH 5, using 1.4% sodium carbonate to stop the reaction. One unit of activity was defined as the amount of enzyme releasing 1 μmol of reducing sugars or pNP per minute under the above conditions. Direct quantification of glucose and cellobiose was accomplished by high-performance liquid chromatography (HPLC) on an Agilent 1200 series system equipped with a refractive index detector. Aliquots of 100 μL were loaded onto a SUPELCOGEL C-G610H column (Sigma) equilibrated with 5 mM H2SO4. The column was previously calibrated by injecting 100 μL of glucose or cellobiose in a concentration range from 0.5 to 60 mM. From the area under the peaks, a calibration curve was calculated for each compound. Peaks were identified from their retention times, by comparison with those of the commercial standards, and their concentrations were calculated from the calibration curves.
Genome sequencing and assembly
DNA was extracted with the DNeasy Plant minikit (Qiagen). Sequencing was performed on an Illumina HiSeqTM 2000 system. 90 bp long reads were obtained by paired-end sequencing of a DNA library composed of 537 ± 396 bp inserts. The Illumina GA Pipeline version 1.5 was used for a first removal of adapter sequences, contamination and low-quality reads from raw reads. Then A5-miseq pipeline  was used in order to carry out a subsequent filtering of the reads, the assembly, and the scaffolding. Assembly was performed de novo and evaluated in terms of N50 and L50. This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession MIKG00000000. The version described in this paper is version MIKG01000000.
Gene prediction, general function annotation, and CAZyME prediction
Gene models were predicted by AUGUSTUS trained with Aspergillus fumigatus . The obtained open reading frames (ORFs) were used as query to interrogate the KEGG and KOG databases using KAAS  and WebMGA  servers, respectively, in order to assign general protein functions profiles. The specific annotation of CAZy families [20, 21] was accomplished by submitting the predicted genes to dbCAN and filtering the results by applying an E value <1·10−20 as cutoff. Then, the assignment of β-glucosidase function was carried out by doing a BLASTP of the predicted glycosyl hydrolases from T. amestolkiae against the characterized GHs from the CAZy database and annotating the function of the best hit with an E value <1·10−20. The presence or absence of signal peptide in all the predicted proteins was analyzed using Phobius, a combined transmembrane topology and signal peptide predictor .
Shotgun analysis of the secretomes
For differential proteomic analysis, samples of the extracellular pool of proteins from culture supernatants of T. amestolkiae grown in the different substrates were independently analyzed. Two mL of the liquid supernatants from 7-day-old cultures were centrifuged at 14,000×g for 5 min to remove insoluble material and freeze-dried. Then, samples containing 5 μg of proteins were first dissolved in sample buffer, denatured, loaded in a 12% SDS-gel and allowed to run for 10 min in order to remove non-protein compounds before the proteomic analysis. The protein band was horizontally divided into two similar portions and each one was excised in small pieces and destained with 50 mM ammonium bicarbonate/50% ACN (acetonitrile), dehydrated with ACN and dried. Samples were reduced by adding DTT to a final concentration of 10 mM and alkylated with iodoacetamide to a final concentration of 50 mM. Then, gel pieces were dried, rehydrated with 12.5 ng/µL trypsin in 50 mM ammonium bicarbonate, and incubated overnight at 30 °C. Peptides were extracted at 37 °C using 100% ACN and 0.5% TFA, respectively. The peptide pool was dried and then cleaned using ZipTip with 0.6 µL C18 resin (Millipore) and reconstituted in 5 µL 0.1% formic acid/2% ACN (v/v) prior to MS analysis.
All peptide separations were carried out on a NanoEasy HPLC (Proxeon Biosystems) coupled to a nanoelectrospray ion source (Proxeon Biosystems). For each analysis, the sample peptides were loaded onto a C18-A1 ASY-Column 2 cm precolumn (Thermo Scientific) and then eluted onto a Biosphere C18 column (C18, inner diameter 75 µm, 15 cm long, 3 µm particle size, Nano Separations). The mobile phase flow rate was 250 nL/min using 0.1% formic acid, 2% ACN in water (solvent A) and 0.1% formic acid and 100% ACN (solvent B). The gradient profile was set as follows: 0–35% solvent B for 120 min, 35–45% solvent B for 20 min, 45–95% solvent B for 9 and 5 min isocratically at 95%. Six microliters of each sample were injected. Full-scan MS spectra (m/z 300–1800) were acquired in the LTQ-Orbitrap Velos in the positive ion mode with a target value of 1,000,000 at a resolution of 60,000 and the 15 most intense ions were selected for collision-induced dissociation (CID) fragmentation in the LTQ with a target value of 10,000 and normalized collision energy of 35%.
Mass spectra *.raw files were searched against an in—house specific database of the T. amestolkiae genome (10408 sequences, 5662098 residues) using the SEQUEST search engine through Proteome Discoverer (version 184.108.40.206, Thermo). Precursor and fragments mass tolerance were set to 10 ppm and 0.5 Da, respectively. Search parameters included a maximum of two missed cleavages allowed, carbamidomethylation of cysteines as a fixed modification and oxidation of methionine as a variable modification. Peptides were validated through the algorithm Percolator  and only those with high confidence were admitted (FDR 0.01). Unless otherwise specified, protein identifications were accepted if they contained at least two identified peptides. Results were inferred from data obtained from two technical replicates from two different biological samples. Relative low precision quantification of the proteins in the samples analyzed was estimated from the sum of the number of peptide spectrum matches (PSMs) corresponding to each protein [24,25,26]. These values were used to calculate the percentage of individual proteins or categories (% of GHs, % of BGLUs) in the whole protein pool. The data provided in tables and figures correspond to the mean value of the %PSMs from the biological replicates.
Wheat straw saccharification
Celluclast 1.5 L FG (Novozymes) was used as commercial cocktail for saccharification. An amount of Celluclast containing 0.5 U of BGLU activity was added to 66.6 mg/mL of acid slurry from wheat straw (40% content in cellulose) in 1.5 mL sodium acetate 100 mM, pH 4, at 50 °C. Samples were taken periodically and the glucose released was quantified using the Glucose TR kit from Spinreact.
Celluclast was supplemented with 0.5 U of BGLU activity from N50010 (Novozyme), TAM377 (T. amestolkiae enzymatic crude produced using glucose as carbon source) or TAM3821 (T. amestolkiae enzymatic crude produced using Avicel as carbon source).
Results and discussion
Cellulase production by T. amestolkiae CIB
Different carbon sources, namely Avicel (microcrystalline cellulose), beechwood xylan, and wheat straw slurry, were tested as inducers of the cellulase activities. In addition, a culture with glucose, which has been described as cellulases’ inhibitor [27, 28], was grown under the same conditions. According to the results shown in Fig. 1, Avicel seemed to be the best substrate for cellulase production. Cellobiohydrolases and β-1,4-endoglucanase activities (Fig. 1a, b) were mainly detected in media containing Avicel or acid wheat straw slurry. Nevertheless, as this last substrate is a mixture of lignocellulose components, the production of cellulolytic enzymes was slightly lower. β-1,4-endoglucanase activity was strongly induced by Avicel (over 10 U/mL after 8 days), being also measurable in cultures containing slurry. BGLU activity reached similar values regardless of the carbon source used (1.4–1.8 U/mL), even with glucose (Fig. 1c). It is well known that, in general, Avicelase and β-1,4-endoglucanase production are repressed in the presence of glucose , but little is known about BGLU production in the presence of easily available carbon sources. Non-inducible BGLU activity has been described in some fungal species [30, 31] but, in most cases, the secretion of this enzyme is induced by cellulose . In addition, the regulation and expression of cellulolytic enzymes can be affected by the presence of inducers, derived from cellulose metabolism, as oligosaccharides or transglycosylation products, or experience catabolic repression by glucose, as reported in T. reesei . Then, the effect of small molecules as cellobiose, glycerol, fructose, and xylose on BGLU production was also assayed (Fig. 2a). In virtually all cases, the released activity was very similar (around 1.4 U/mL) and comparable to that obtained before for polymeric cellulosic substrates and glucose, getting lower values for cellobiose, the only disaccharide evaluated (1 U/mL), which needs to be transformed in glucose before being consumed by the fungus. These data demonstrate that this T. amestolkiae strain does not require specific inducers for BGLU production.
Monitoring glucose consumption in cultures containing this monosaccharide as carbon source revealed the complete depletion of the sugar after 24 h (Fig. 2b), which may indicate that, in this case, carbon starvation is triggering BGLU production. Daily addition of 0.5 or 1% glucose to cultures in Mandels medium containing either 1% glucose or Avicel (Fig. 2c) caused only basal production of BGLU, confirming that the production of this enzyme is repressed by glucose. Thus, when easily assimilable carbon sources are exhausted, carbon starvation might induce BGLU production in T. amestolkiae, while the production of CBH and EG is induced only when cellulosic substrates are present.
However, these data raised the question of whether the BGLU secreted under carbon starvation and in the presence of cellulose are the same enzyme or not. Preliminary results from zimograms after isoelectric focusing suggested that they could be different proteins (not shown). The production of basal levels of cellulases, usually Avicelase and β-1,4-endoglucanase, has been described and interpreted as a fungal strategy intended to take advantage of any cellulosic material present in the surrounding medium . On the contrary, the production of BGLU in a carbon source-independent way could be due to its role in the regulation of cellulase production. In this sense, BGLU from T. reesei can convert cellobiose into sophorose, a strong cellulase inducer .
In this work, the analysis of cellulase diversity in T. amestolkiae and its differential production in specific conditions was tackled according to several strategies, as reported above.
General features of T. amestolkiae CIB genome
In a first approach, the cellulolytic potential of the fungus was analyzed through the sequencing and annotation of its genome. The draft genome sequence of T. amestolkiae CIB was based on high-throughput sequencing system (~65-fold coverage) and de novo assembling. The 24,330,860 pair-ended reads were incorporated into 215 scaffolds, among which 132 comprised more than 1 kb. Genome size was determined to be 33.7 Mb and N50 and L50 statistics displayed values of 1,486,010 and 9, respectively. The assembled genome resulted in the prediction of 10,408 ORFs (Additional file 2). As seen below, 342 of these ORFs were supported by the secretome data. The prediction of GH enzymes was performed by submitting the putative ORFs to the dbCAN server and filtering the results with the selected cutoff. By this way, 325 glycosyl hydrolases were identified, a number close to those described for Fusarium verticillioides (332), Aspergillus flavus (334), and Aspergillus oryzae (317) .
β-glucosidases are highly represented in the genome of T. amestolkiae
Glycosyl hydrolases involved in hydrolysis of lignocellulose are widely extended across the CAZy families . Figure 3 depicts the profile of the main families of CAZymes implicated in plant cell-wall degradation encoded by the genomes of T. amestolkiae and six microorganisms of remarkable lignocellulolytic interest. The selection of these species was made attending to the criteria of relevance and availability of accurately annotated genomic data. T. reesei is the main industrial source of the cellulases and hemicellulases added to commercial cocktails for lignocellulose saccharification [6, 35]. Other fungi, as Aspergillus niger and Penicillium oxalicum (previously classified as Penicillium decumbens) are also industrial producers of cellulolytic enzymes [36, 37]. Thielavia terrestris secretes thermostable hydrolases , and Clostridium thermocellum is a well-known Firmicutes capable of directly converting cellulose into ethanol and other value-added products [39, 40].
Data analysis revealed specific CAZyme profiles for each species although, as expected, the bacterium displayed the biggest differences. Distinctive CAZyme profiles among plant cell-wall-degrading fungi have been linked to nutritional mode adaptations and may imply different strategies for attacking lignocellulose at the enzymatic level [34, 41, 42]. Regarding T. amestolkiae, the number of encoded CAZymes was significantly higher than in the other organisms. In particular, the high amount of genes encoding enzymes from the GH3 family, which comprises many of the reported fungal β-glucosidases, was intriguing and it could be related to the findings of BGLU production reported above. The presence of this activity in the culture supernatants of this organism, regardless of the carbon source, may suggest a central role of these enzymes in the metabolic strategy adopted by T. amestolkiae along evolution.
To assess this possibility, we carried out a deeper comparison focusing on GH1, GH3, GH5, and GH30, the main GH families encoding fungal β-glucosidases (Fig. 4a). The total number of genes for these enzymes in the genome of T. amestolkiae was annotated by running BLASTP against the characterized GHs from the CAZy database, and compared with those reported for each of the cellulolytic species used as reference (Fig. 4b). According to our BLASTP analysis, 5 BGLUs from GH1 family, 18 from GH3, and 1 from GH5 are codified in the genome of this isolate (Table 1).
The CAZymes types and abundance represented in a given fungal genome may condition a certain strategy for hydrolyzing lignocellulose . The high number of β-glucosidases annotated in T. amestolkiae pointed them as key factors during lignocellulose degradation by this fungus. This is especially relevant considering that commercial mixtures based on T. reesei crudes (the main industrial source of this type of enzymes) are deficient in β-glucosidase activity . In this sense, this unexpectedly high number of cellulases and hemicellulases in T. amestolkiae opens the possibility of improving lignocellulose saccharification through the synergistic addition of several CAZymes.
The expression of β-glucosidases by T. amestolkiae in the presence of different lignocellulosic substrates was validated, beyond the genomic level, by proteomic analysis of the extracellular proteins secreted by this fungus.
Differential shotgun analysis of proteins secreted by T. amestolkiae
The enzymatic pools released under the four assayed conditions were compared by sequential tryptic digestion and LC–MS/MS of the whole peptide mixtures produced (Additional file 3). The number of proteins identified in each secretome differed according to the complexity of the carbon source used for fungal growth. The maximum number of identifications came from samples produced in the two polymeric substrates, xylan (184) and slurry (144), while only 119 and 104 proteins were recognized in the secretomes from glucose and Avicel cultures, respectively. Xylan and slurry are structurally complex materials that contain heterogeneous and branched polymers. Hence, it is expectable that the amount and variety of extracellular enzymes required to accomplish their degradation is higher than those needed to metabolize simpler substrates.
Table 2 lists the ten extracellular proteins identified with maximal confidence (with the highest scores) in each condition. Mainly CBH, followed by xylanases and BGLU, was produced when the fungus grew in the medium with Avicel, and the same distribution was observed in cultures with wheat straw slurry. The presence of xylanases in the secretome of P. chrysosporium grown in cellulose has already been reported . β-xylosidases were detected only in xylan cultures, as previously described . Surprisingly, a glutaminase was expressed in all conditions. This protein might be released upon induction of protein catabolism by early carbon starvation , and its production in Thermomyces lanuginosus SSBP cultures in corncob has been previously described . Other extracellular proteins with different functionalities were also detected in glucose cultures. Some of them are probably involved in the metabolism of fungal carbohydrates, as an hexosaminidase . For example, cellulase regulation in H. jecorina proceeds through deglycosylation catalyzed by glucosaminidases, which can modulate enzyme activity by this mechanism . Under starvation stress, the increased production of other enzymes involved in fungal cell-wall degradation, as β-1,3-glucanases, has also been observed. The heterogeneity of the proteins identified in the extracellular pool of proteins from glucose cultures is in agreement with the high variability detected by KOG functional analysis (see below). The PSM values of two independent experiments were taken as low precision, semi-quantitative records for further evaluation of the shotgun analysis. These values account for the number of scans identified for each protein, which are roughly related to the amount of a given protein in the sample , being analyzed from different perspectives: functional analysis (KOG attribution), GHs diversity, and relative abundance of cellulases and hemicellulases.
Functional analysis of the proteins secreted in the presence of different carbon sources
Functional analysis based on KOG categories (Table 3) showed that most proteins detected in all conditions tested were implicated in carbohydrate metabolism and transport (70.5–91.8%). The second most abundant group clustered enzymes involved in amino acid metabolism and transport. The rest of proteins from each one of the secretomes studied belonged to a different number of KOG categories. For each condition, the categories that group more than 2% of the total proteins identified were considered to be the most representative. These outlines describe perfectly the results deduced for Avicel cultures, in which more than 95% of the extracellular proteins are categorized in the two main groups. The enzymes in the other three secretomes studied share significant amounts of proteins involved in signal transduction and transport, energy production and conversion. The extracellular proteins from glucose cultures had the lowest representation of enzymes involved in carbohydrate metabolism and a significantly high amount of proteins related to cell-wall biogenesis, as compared to the other samples. Similarly, the notable representation of enzymes engaged in post-translational modification, protein turnover, and with chaperone functions in the proteins from glucose, xylan, and slurry supernatants merits especial attention. As commented before, the number of proteins detected in the shotgun analysis of the last two secretomes is considerably high due to the need of using a battery of enzymes to metabolize these complex and heterogeneous lignocellulosic materials, and this fact would also justify the over-representation of enzymes involved in the synthesis and secretion of the proteins required to accomplish this task. Although these are mainly intracellular proteins, they can be detected in the culture supernatant due to fungal autolysis. In fact, when the sequences of the proteins identified were examined for the presence of signal peptide, it turned out that more than 30% of proteins in glucose, xylan, and slurry supernatants were intracellular, which contrasts with the 14% of intracellular proteins in Avicel supernatants. These data coincide with the results from semi-quantitative analysis (% PSM) that indicated that more than 90% of the total protein amount from Avicel supernatants corresponds to extracellular enzymes, against 74–79% of extracellular proteins in xylan, glucose, and slurry. Similarly, the classification of the carbohydrolases identified into GH families reflected 88% of extracellular enzymes in the medium with Avicel (95% PSM), and 74, 72 and 82% (89, 78 and 89% PSM) in the supernatants from glucose, xylan, and slurry, respectively. More than 5% of the proteins detected in slurry supernatants were categorized as energy production and conversion enzymes, and may have been secreted due to the presence of lignin in the culture medium, an heterogeneous substrate that comes from plant cell wall and is catabolized by oxidases . The number of different proteins identified in the extracellular medium recovered upon fungal growth in the four carbon sources was analyzed by a Venn representation (Fig. 5). According to these data, 37 proteins were produced in all conditions evaluated, and 21, 5, 44, and 15 enzymes were exclusive from Avicel, glucose, xylan, and slurry supernatants, respectively. These data agree with those accounted for total proteins, since the maximum number of total proteins was also found in xylan cultures.
In the present work, the proteomic analysis was carried out in samples from 7-day-old cultures. As already mentioned, glucose was depleted after 24 h growth, which corresponds to severe carbon starvation in these cultures. On the other hand, the scarce solubility of xylan and slurry components may lead to an overall pseudo-starvation, as oligo- and monosaccharides are presumably slowly released and quickly metabolized by the fungus. Similar cellulase induction has been reported after depletion of the carbon source and during early exposure to cellulose in Neurospora crassa and A. niger cultures . Our data are consistent with those obtained for the extracellular proteins detected in carbon-starved cultures of the ascomycete Aspergillus niger that also cluster into multiple KOG categories .
GHs are the main extracellular enzymes secreted by T. amestolkiae
The major enzymes secreted in all conditions tested are GHs, which constitute a very complex group of carbohydrate-active enzymes. Their current classification into families is entirely based on the similarity of proteins’ sequence, and two enzymes with the same activity can be included in different families.
As explained before, low precision quantitative data for GHs were also based on the number of PSMs (mean value from the two biological replicates). Total GH abundance was first calculated as the sum of PSMs from all GHs identified in each individual sample, and then the percentage of each enzyme was obtained. The GHs identified in the proteomic analysis were categorized into 31 families (Table 4). Proteins from family GH3, associated to degradation of plant cell-wall polysaccharides (β-glucosidases, β-xylosidases, arabinofuranosidases and exo-1,3-1,4-glucanases), are produced in all carbon sources tested. In terms of relative abundance (% PSM), these enzymes are the main extracellular GHs found in the xylan-containing medium (>25%), representing 16.3% and 13.7% in glucose and slurry media and less than 10% in the Avicel crude. GH7 proteins, which include endo-β-1,4-glucanases and reducing end-acting cellobiohydrolases, were profusely produced in the two cellulose-containing media, Avicel (47.1% of the total GHs) and slurry (29.5%), while GH6 endoglucanases and non-reducing end cellobiohydrolases were especially abundant in the secretome from the Avicel medium (15.1%). Proteins from these families were induced by xylan to a lesser extent (1.5% for GH6 and 5.6% for GH7). Family GH10 groups enzymes with endo-xylanase activity and is well represented in Avicel, slurry, and xylan crudes. In addition, α-carbohydrolases belonging to GH15 and GH31 families were identified in the samples analyzed. The huge amount of GH15 glucoamylases/glucodextranases detected in glucose cultures (28%) can be related to fungal autolysis, since α-(1-3) and α-(1-4) glucans have been reported as cell-wall components of several Talaromyces species . On the other hand, GH31 proteins act on α-linked xylo-, manno-, galacto-, or gluco-oligosaccharides contributing to the final breakdown of oligomers released from plant hemicellulose or fungal polysaccharides.
The percentages of the different types of cellulases and hemicellulases secreted in each medium are summarized in Fig. 6.
β-Glucosidases produced in different culture conditions
Further analysis of the proteomic data allowed the evaluation of the β-glucosidases induced in the carbon sources tested (Fig. 7). An approximation to the relative abundance of each one of them, calculated as the percentage of total β-glucosidases, was analyzed. The number of hypothetical BGLUs detected varied from 6 different enzymes in Avicel and glucose media to 8 and 10 for slurry and xylan cultures, respectively. Protein g377 (see Additional file 1 for sequence) was the most abundant in media containing slurry, glucose, or xylan, representing 72, 55, and 45% of the total BGLUs, respectively, while in Avicel, the main BGLU was protein g3821 (70%). However, low amounts of this enzyme were also produced in media with glucose (<2%), xylan (5%), or slurry (13%). These findings suggest that BGLU g3821 is strongly induced by cellulose but constitutively expressed in basal levels in the other conditions studied. On the contrary, the production of BGLU g377 seems to be carbon source-independent, although its synthesis and secretion may be triggered upon starvation. Preliminary RT-PCR experiments indicated that g3821 is overexpressed in Avicel cultures and repressed when glucose was used as the carbon source, while g377 is overexpressed in both media, which is compatible with a constitutive expression, independent of the presence of cellulose in the medium (data not shown).
The enzyme g377 could be responsible for the BGLU activity detected when T. amestolkiae was cultivated using non-cellulosic materials as carbon sources (mono and disaccharides). Similarly to this, one of the two different β-glucosidases produced by Stachybotrys atra is non-inducible . Its role could also be related to cell-wall remodeling, as seen in some plants .
T. amestolkiae crude supernatants as BGLU sources for saccharification
Preliminary saccharification assays were performed using Celluclast 1.5 L FG or T. amestolkiae crudes from Avicel as the main sources of cellulases. Released glucose was 33% higher when Celluclast 1.5 L FG was used, probably for its high content in CBH and EG compared to T. amestolkiae crudes, rich in BGLU.
For that reason, enzymatic hydrolysis of wheat straw slurry was assayed using Celluclast 1.5 L FG as the main source of cellulases. As BGLU supplements for saccharification, 0.5 units of this activity were incorporated to the saccharification mixture using T. amestolkiae crudes from glucose (TAM377) or Avicel (TAM3821) culture media, or the commercial preparation N50010 from Novozymes (Fig. 8).
The treatment of wheat straw slurry with only Celluclast released 50% glucose after 96 h. In the same time period, glucose release enhanced significantly when Celluclast was supplemented with TAM3821 (80%), N50010 (73%) or TAM377 (61%).
These data show that these enzymatic crudes from T. amestolkiae are very efficient as supplements for lignocellulosic biomass saccharification, especially TAM3821 that provided superior glucose solubilization than N50010. This could be due to a high affinity of BGLU g3821 for cellulosic substrates, since this protein was only induced by cellulose. Nevertheless, the existence of synergistic catalysis between this BGLU and the EGs and CBHs induced by cellulosic substrates cannot be ruled out. Further purification and characterization of both BGLU will be carried out in order to ascertain their different catalytic constants, substrate specificities, and efficiencies.
A strong extracellular β-glucosidase activity was detected in T. amestolkiae cultures grown in four carbon sources of different composition and complexity. Genome sequencing and annotation disclosed a high number of genes encoding BGLUs, which suggests the relevance of these enzymes for this organism. Differential semi-quantitative proteomic analysis of the secretomes recovered from each culture revealed the production of two main proteins with BGLU activity. One of them was constitutively released regardless the carbon source used, although it was repressed by glucose, and showed to be less efficient than the commercial BGLU preparation used for comparison in saccharification experiments. The other was strongly induced by cellulose and revealed to be an excellent supplement to Celluclast for saccharification of pretreated wheat straw.
cellobiohydrolase, Avicelase, exocellulase
liquid chromatography coupled to tandem mass spectrometry
high performance liquid chromatography
Kyoto Encyclopedia of Genes and Genomes
EuKaryotic Orthologous Groups
open reading frames
peptide spectrum match
culture supernatant of T. amestolkiae growing on Mandels plus glucose, expressing more g377 BGLU
culture supernatant of T. amestolkiae growing on Mandels plus Avicel, expressing more g3821 BGLU
Naik SN, Goud VV, Rout PK, Dalai AK. Production of first and second generation biofuels: a comprehensive review. Renew Sustain Energy Rev. 2010;14:578–97.
Jurado M, Prieto A, Martinez-Alcala A, Martinez AT, Jesus Martinez M. Laccase detoxification of steam-exploded wheat straw for second generation bioethanol. Bioresour Technol. 2009;100:6378–84.
Alvira P, Tomas-Pejo E, Ballesteros M, Negro M. Pretreatment technologies for an efficient bioethanol production process based on enzymatic hydrolysis: a review. Bioresour Technol. 2010;101:4851–61.
Zhang YHP, Lynd LR. Toward an aggregated understanding of enzymatic hydrolysis of cellulose: noncomplexed cellulase systems. Biotechnol Bioeng. 2004;88:797–824.
Duff SJB, Murray WD. Bioconversion of forest products industry waste cellulosics to fuel ethanol: a review. Bioresour Technol. 1996;55:1–33.
Martinez D, Berka RM, Henrissat B, Saloheimo M, Arvas M, Baker SE, Chapman J, Chertkov O, Coutinho PM, Cullen D, Danchin EGJ, Grigoriev IV, Harris P, Jackson M, Kubicek CP, Han CS, Ho I, Larrondo LF, de Leon AL, Magnuson JK, Merino S, Misra M, Nelson B, Putnam N, Robbertse B, Salamov AA, Schmoll M, Terry A, Thayer N, Westerholm-Parvinen A, et al. Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina). Nat Biotechnol. 2008;26:553–60.
Maeda RN, Serpa VI, Rocha VAL, Mesquita RAA, Santa Anna LMM, de Castro AM, Driemeier CE, Pereira N, Polikarpov I. Enzymatic hydrolysis of pretreated sugar cane bagasse using Penicillium funiculosum and Trichoderma harzianum cellulases. Process Biochem. 2011;46:1196–201.
Bhiri F, Chaabouni SE, Limam F, Ghrir R, Marzouki N. Purification and biochemical characterization of extracellular beta-glucosidases from the hypercellulolytic pol6 mutant of Penicillium occitanis. Appl Biochem Biotechnol. 2008;149:169–82.
Szijarto N, Szengyel Z, Liden G, Reczey K. Dynamics of cellulase production by glucose grown cultures of Trichoderma reesei rut-C30 as a response to addition of cellulose. Appl Biochem Biotechnol. 2004;113:115–24.
Adsul M, Bastawde K, Varma A, Gokhale VD. Strain improvement of Penicillium janthinellum NCIM 1171 for increased cellulase production. Bioresour Technol. 2007;98:1467–73.
Chávez R, Bull P, Eyzaguirre J. The xylanolytic enzyme system from the genus Penicillium. J Biotechnol. 2006;123:413–33.
Ustinov BB, Gusakov AV, Antonov AI, Sinitsyn AP. Comparison of properties and mode of action of six secreted xylanases from Chrysosporium lucknowense. Enzyme Microb Technol. 2008;43:56–65.
Mandels M, Weber J. Production of cellulases. Adv Chem Ser. 1969;95:391–414.
Nieto-Domínguez M, de Eugenio LI, Barriuso J, Prieto A, Fernández de Toro B, Canales-Mayordomo Á, Martínez MJ. Novel pH-stable glycoside hydrolase family 3 β-Xylosidase from Talaromyces amestolkiae: an enzyme displaying regioselective transxylosylation. Appl Environ Microbiol. 2015;81:6380–92.
Somogyi M. A new reagent for the determination of sugars. J Biol Chem. 1945;160:61–8.
Tritt A, Eisen JA, Facciotti MT, Darling AE. An integrated pipeline for de novo assembly of microbial genomes. PLoS ONE. 2012;7:e42304.
Hoff KJ, Stanke M. WebAUGUSTUS—a web service for training AUGUSTUS and predicting genes in eukaryotes. Nucleic Acids Res. 2013;41(Web Server issue):W123–8.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35(Web Server issue):W182–5.
Wu S, Zhu Z, Fu L, Niu B, Li W. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genom. 2011;12:444.
Carbohydrate active enzymes database. http://www.cazy.org/.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(Database issue):D490–5.
Käll L, Krogh A, Sonnhammer EL. A combined transmembrane topology and signal peptide prediction method. J Mol Biol. 2004;338:1027–36.
Käll L, Canterbury JD, Weston J, Noble WS, MacCoss MJ. Semi-supervised learning for peptide identification from shotgun proteomics datasets. Nat Methods. 2007;4:923–5.
Liu H, Sadygov RG, Yates JR. A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem. 2004;76:4193–201.
Schulze WX, Usadel B. Quantitation in mass-spectrometry-based proteomics. Annu Rev Plant Biol. 2010;61:491–516.
Salvachúa D, Katahira R, Cleveland NS, Khanna P, Resch MG, Black BA, Purvine SO, Zink EM, Prieto A, Martínez MJ, Martínez AT, Simmons BA, Gladden JM, Beckham GT, Martínez AT, Speranza M, Ruiz-Dueñas F, Ferreira P, Camarero S, Guillén F, Martínez MJ, Gutiérrez A, del Río J, Salvachúa D, Prieto A, Martínez AT, Martínez MJ, Picart P, de María PD, Schallmey A, et al. Lignin depolymerization by fungal secretomes and a microbial sink. Green Chem. 2016;18:6046–62.
Chen H. Biological fundamentals for the biotechnology of lignocellulose. In: Chen H, editor. Biotechnology of lignocellulose—theory and practice. Dordrecht: Springer; 2014. p. 73–141.
Hanif A, Yasmeen A, Rajoka MI. Induction, production, repression, and de-repression of exoglucanase synthesis in Aspergillus niger. Bioresour Technol. 2004;94:311–9.
Amore A, Giacobbe S, Faraco V. Regulation of cellulase and hemicellulase gene expression in fungi. Curr Genomics. 2013;14:230–49.
Shewale JG. Beta-Glucosidase: its role in cellulase synthesis and hydrolysis of cellulose. Int J Biochem. 1982;14:435–43.
Lee J, Kwon K-S, Hah YC. Regulation of β-glucosidase biosynthesis in Aspergillus nidulans. FEMS Microbiol Lett. 1996;135:79–84.
Ilmén M. Saloheimo a, Onnela ML, Penttilä ME: Regulation of cellulase gene expression in the filamentous fungus Trichoderma reesei. Appl Environ Microbiol. 1997;63:1298–306.
Fowler T, Brown RD. The bgl1 gene encoding extracellular beta-glucosidase from Trichoderma reesei is required for rapid induction of the cellulase complex. Mol Microbiol. 1992;6:3225–35.
Zhao Z, Liu H, Wang C, Xu J-R, Cantarel B, Coutinho P, Rancurel C, Bernard T, Lombard V, Henrissat B, Ospina-Giraldo M, Griffith J, Laird E, Mingora C, Murphy C, Powlowski J, Wu M, Butler G, Tsang A, Couturier M, Navarro D, Olive C, Chevret D, Haon M, Favel A, Lesage-Meessen L, Henrissat B, Coutinho P, Berrin J, King B, et al. Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 2013;14:274.
Häkkinen M, Arvas M, Oja M, Aro N, Penttilä M, Saloheimo M, Pakula TM, Sánchez C, McKendry P, Jordan D, Bowman M, Braker J, Dien B, Hector R, Lee C, Mertens J, Wagschal K, Sukumaran R, Singhania R, Mathew G, Pandey A, Cantarel B, Coutinho P, Rancurel C, Bernard T, Lombard V, Henrissat B, Schuster A, Schmoll M, Teeri T, et al. Re-annotation of the CAZy genes of Trichoderma reesei and transcription in the presence of lignocellulosic substrates. Microb Cell Fact. 2012;11:134.
Liu G, Zhang L, Wei X, Zou G, Qin Y, Ma L, Li J, Zheng H, Wang S, Wang C, Xun L, Zhao G-P, Zhou Z, Qu Y. Genomic and secretomic analyses reveal unique features of the lignocellulolytic enzyme system of Penicillium decumbens. PLoS ONE. 2013;8:e55185.
Pel HJ, deWinde JH, Archer DB, Dyer PS, Hofmann G, Schaap PJ, Turner G, deVries RP, Albang R, Albermann K, Andersen MR, Bendtsen JD, Benen JAE, van den Berg M, Breestraat S, Caddick MX, Contreras R, Cornell M, Coutinho PM, Danchin EGJ, Debets AJM, Dekker P, van Dijck PWM, van Dijk A, Dijkhuizen L, Driessen AJM, d’Enfert C, Geysens S, Goosen C, Groot GSP, et al. Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88. Nat Biotechnol. 2007;25:221–31.
Berka RM, Grigoriev IV, Otillar R, Salamov A, Grimwood J, Reid I, Ishmael N, John T, Darmond C, Moisan M-C, Henrissat B, Coutinho PM, Lombard V, Natvig DO, Lindquist E, Schmutz J, Lucas S, Harris P, Powlowski J, Bellemare A, Taylor D, Butler G, de Vries RP, Allijn IE, van den Brink J, Ushinsky S, Storms R, Powell AJ, Paulsen IT, Elbourne LDH, et al. Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris. Nat Biotechnol. 2011;29:922–7.
Stevenson DM, Weimer PJ. Expression of 17 genes in Clostridium thermocellum ATCC 27405 during fermentation of cellulose or cellobiose in continuous culture. Appl Environ Microbiol. 2005;71:4672–8.
Roberts SB, Gowen CM, Bruks J, Fong SS. Genome-scale metabolic analysis of Clostridium thermocellum for bioethanol production. BMC Systems Biol. 2010;4:31.
Znameroski EA, Glass NL. Using a model filamentous fungus to unravel mechanisms of lignocellulose deconstruction. Biotechnol Biofuels. 2013;6:6.
Martinez D, Challacombe J, Morgenstern I, Hibbett D, Schmoll M, Kubicek CP, Ferreira P, Ruiz-Duenas FJ, Martinez AT, Kersten P, Hammel KE, Vanden Wymelenberg A, Gaskell J, Lindquist E, Sabat G, Bondurant SS, Larrondo LF, Canessa P, Vicuna R, Yadav J, Doddapaneni H, Subramanian V, Pisabarro AG, Lavin JL, Oguiza JA, Master E, Henrissat B, Coutinho PM, Harris P, Magnuson JK, et al. Genome, transcriptome, and secretome analysis of wood decay fungus Postia placenta supports unique mechanisms of lignocellulose conversion. Proc Natl Acad Sci USA. 2009;106:1954–9.
Kumar R, Singh S, Singh OV. Bioconversion of lignocellulosic biomass: biochemical and molecular perspectives. J Ind Microbiol Biotechnol. 2008;35:377–91.
Hori C, Igarashi K, Katayama A, Samejima M. Effects of xylan and starch on secretome of the basidiomycete Phanerochaete chrysosporium grown on cellulose. FEMS Microbiol Lett. 2011;321:14–23.
van Munster JM, Daly P, Delmas S, Pullan ST, Blythe MJ, Malla S, Kokolski M, Noltorp ECM, Wennberg K, Fetherston R, Beniston R, Yu X, Dupree P, Archer DB. The role of carbon starvation in the induction of enzymes that degrade plant-derived carbohydrates in Aspergillus niger. Fungal Genet Biol. 2014;72:34–47.
Winger AM, Heazlewood JL, Chan LJG, Petzold CJ, Permaul K, Singh S. Secretome analysis of the thermophilic xylanase hyper-producer Thermomyces lanuginosus SSBP cultivated on corn cobs. J Ind Microbiol Biotechnol. 2014;41:1687–96.
Gooday G. What are the roles of chitinases in the growing fungus? FEMS Microbiol Lett. 1992;100:387–91.
Stals I, Karkehabadi S, Kim S, Ward M, Van Landschoot A, Devreese B, Sandgren M. High resolution crystal structure of the endo-N-acetyl-β-d-glucosaminidase responsible for the deglycosylation of Hypocrea jecorina cellulases. PLoS ONE. 2012;7:e40854.
Roberts TC, Johansson HJ, McClorey G, Godfrey C, Blomberg KEM, Coursindel T, Gait MJ, Smith CIE, Lehtiö J, El Andaloussi S, Wood MJA. Multi-level omics analysis in a murine model of dystrophin loss and therapeutic restoration. Hum Mol Genet. 2015;24:6756–68.
Nitsche BM, Jørgensen TR, Akeroyd M, Meyer V, Ram AFJ. The carbon starvation response of Aspergillus niger during submerged cultivation: insights from the transcriptome and secretome. BMC Genomics. 2012;13:380.
Prieto A, Bernabé M, Leal JA. Isolation, purification and chemical characterization of alkali-extractable polysaccharides from the cell walls of Talaromyces species. Mycol Res. 1995;99:69–75.
De Gussem RL, Aerts GM, Claeyssens M, De Bruyne CK. Purification and properties of an induced β-d-glucosidase from Stachybotrys atra. Biochim Biophys Acta Enzymol. 1978;525:142–53.
Ketudat Cairns JR, Mahong B, Baiya S, Jeon J-S. β-Glucosidases: multitasking, moonlighting or simply misunderstood? Plant Sci. 2015;241:246–59.
LE, JM, JG, JB, and MN contributed to design the research and participate in the experiments. LA and MN analyzed and assembled fungal genome. LE and MN drafted the manuscript. MJM and AP co-coordinated the research, and reviewed and edited the manuscript. All authors read and approved the final manuscript.
The authors thank the Proteomics and Genomics facility at CIB and also Ph.D. Alexander Holm Viborg for kindly providing the bioinformatics tool (http://research.ahv.dk/cazy) used for characterized GHs sequences extraction from CAZy database.
The authors declare that they have no competing interests.
Availability of supporting data
This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession MIKG00000000. The version described in this paper is MIKG01000000.
This work has been funded by projects RTC-2014-1777-3 and BIO2015-68387-R from MINECO and S2013/MAE-2907 from Comunidad de Madrid. M. Nieto and J.A. Méndez thank their fellowships, FPU of the MINECO and Tatiana Pérez de Guzmán el Bueno Foundation, respectively.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Identification of the fungal isolate. This file contains three figures, the materials and methods associated to the information presented, and a brief discussion of the data presented. Figure S1. Maximum likelihood phylogenetic analysis of RPB1 (A), ITS (B), and BT2 (C) regions from different Talaromyces strains. Figure S2. Agar colonies of T. amestolkiae, Penicillium purpurogenum var. rubrisclerotium and Penicillium rubrum. Figure S3. SEM micrography of conidiophores and hyphae from the three fungal species.
Additional file 2. Gene prediction from T. amestolkiae genome. The file contains the amino acid sequences of the proteins encoded by the genome of T. amestolkiae.
Additional file 3. Proteins in secretomes of T. amestolkiae growing with different carbon sources: 3A) Avicel; 3B) Glucose; 3C) Slurry; 3D) Xylan. The file contains four tables with the list of proteins identified in each condition.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
de Eugenio, L.I., Méndez-Líter, J.A., Nieto-Domínguez, M. et al. Differential β-glucosidase expression as a function of carbon source availability in Talaromyces amestolkiae: a genomic and proteomic approach. Biotechnol Biofuels 10, 161 (2017). https://doi.org/10.1186/s13068-017-0844-7
- Lignocellulosic biomass