Research article | Open | Published:
Closely related fungi employ diverse enzymatic strategies to degrade plant biomass
Biotechnology for Biofuelsvolume 8, Article number: 107 (2015)
Plant biomass is the major substrate for the production of biofuels and biochemicals, as well as food, textiles and other products. It is also the major carbon source for many fungi and enzymes of these fungi are essential for the depolymerization of plant polysaccharides in industrial processes. This is a highly complex process that involves a large number of extracellular enzymes as well as non-hydrolytic proteins, whose production in fungi is controlled by a set of transcriptional regulators. Aspergillus species form one of the best studied fungal genera in this field, and several species are used for the production of commercial enzyme cocktails.
It is often assumed that related fungi use similar enzymatic approaches to degrade plant polysaccharides. In this study we have compared the genomic content and the enzymes produced by eight Aspergilli for the degradation of plant biomass. All tested Aspergilli have a similar genomic potential to degrade plant biomass, with the exception of A. clavatus that has a strongly reduced pectinolytic ability. Despite this similar genomic potential their approaches to degrade plant biomass differ markedly in the overall activities as well as the specific enzymes they employ. While many of the genes have orthologs in (nearly) all tested species, only very few of the corresponding enzymes are produced by all species during growth on wheat bran or sugar beet pulp. In addition, significant differences were observed between the enzyme sets produced on these feedstocks, largely correlating with their polysaccharide composition.
These data demonstrate that Aspergillus species and possibly also other related fungi employ significantly different approaches to degrade plant biomass. This makes sense from an ecological perspective where mixed populations of fungi together degrade plant biomass. The results of this study indicate that combining the approaches from different species could result in improved enzyme mixtures for industrial applications, in particular saccharification of plant biomass for biofuel production. Such an approach may result in a much better improvement of saccharification efficiency than adding specific enzymes to the mixture of a single fungus, which is currently the most common approach used in biotechnology.
Plant biomass is a highly attractive substrate for the production of biofuels and biochemicals as it is a renewable resource with a lower carbon footprint as compared to fossil substrates. It is also the predominant carbon source for most fungi and consists largely of polymeric compounds, of which polysaccharides are the main components [1, 2]. In addition, lignin encrusts the polysaccharides and acts as a physical barrier that impedes fungal enzymes from gaining access to them. Fungi cannot take up intact polysaccharides, but need to degrade them extracellularly to monomeric and oligomeric compounds using diverse enzymatic mixtures . Plant polysaccharide degradation by fungi has been a topic of study for many decades due to its relevance in many industrial applications, such as paper and pulp, food and feed, beverages, textiles and detergents. The increasing interest in the production of alternative fuels and chemicals from plant biomass has provided an even greater push for research into fungal decomposition of plant biomass. In contrast to the earlier applications, production of biofuels and biochemicals would ideally involve complete depolymerization of the polysaccharides to monomers, putting a much higher demand on the efficiency of the enzymatic mixtures employed by industry.
Analysis of an increasing number of fungal genome sequences has demonstrated the fundamental differences in the plant polysaccharide degrading machinery of fungi [3–8]. In addition, the regulatory systems that control plant biomass degradation also differ strongly among fungi, although they are largely conserved among different Aspergillus species [9–16]. Results from a previous study on the utilization of polysaccharides by three Aspergilli  suggest that related fungal species may have developed different approaches to plant biomass degradation. In nature, biomass-degrading fungi live in mixed communities with other microorganisms. It can be expected that different species target distinct components of the substrate and degrade them using dissimilar enzyme combinations. An enhanced understanding of these strategies will not only increase our knowledge of fungal biodiversity, but will help in the design of efficient industrial enzyme mixtures for plant biomass degradation. In this study, we compared the plant biomass degradation potential and approaches of eight Aspergillus species: A. clavatus, A. fischeri, A. flavus, A. fumigatus, A. nidulans, A. niger, A. oryzae and A. terreus (Additional file 1: Table S1). The main aim was to evaluate if these related fungi have significant differences in their approach to degrade plant biomass and if this could provide leads to improve the saccharification efficiency of commercial enzyme cocktails. To do this we compared the genomic potential of these fungi as well as the enzymes sets they produce during growth on two common feedstocks, wheat bran and sugar beet pulp, that differ significantly in their composition (Table 1).
Genomic potential of the studied Aspergilli related to plant biomass utilization
Based on the Carbohydrate-Active enZymes (CAZy)  annotation pipeline, total numbers of glycoside hydrolases (GH), polysaccharide lyases (PL) and carbohydrate esterases (CE) vary among the species (Fig. 1; Table 2). The percentage of GH genes related to plant polysaccharide degradation (PPD) is 58–66% for all genomes, except that A. clavatus has 20–30% less GH genes than the others (Fig. 1), largely due to a reduction in pectinases (GH28, GH54, GH78, GH88) (Table 2). A. clavatus also contains the lowest percentage of PPD-related PL genes (71% as compared to >86%), which are also all related to pectin degradation. The variations in CAZy content are relatively small compared to previous studies with a more diverse set of fungal species [3–8]. This can be explained by their close phylogenetic relationships and their similar habitats, which would push genome evolution in a similar direction.
When the Aspergillus genomes were compared for individual CAZy families, significant differences in numbers of genes were observed (Additional file 2: Table S2A, B). Variations in gene numbers are particularly obvious in certain CAZy families involved in the degradation of mannan (GH26), pectin (GH28, GH53, GH78, GH88, GH93, PL1, PL3, CE8 and CE12), xyloglucan (GH29 and GH74), starch (GH31), sucrose/inulin (GH32), cellulose (GH45 and AA9), and xylan (GH115 and CE15). Genes encoding lignin-modifying peroxidases are not present in any of these genomes, but significant differences are found in the number of laccases and other oxidoreductase enzymes, which may play a role in lignin or polysaccharide degradation (Additional file 2: Table S2B). A. niger is richest in laccases (13 in its genome), while the other species have two to nine.
Orthologous clustering of the CAZymes showed that only 14.7% of the genes encoding hydrolytic enzymes are shared by all species (Additional file 2: Table S3A). In contrast, 27.5% of the genes are unique to a single species, with the largest number in A. nidulans, A. niger and A. terreus. For the oxidative enzymes, 10.8% of the genes are shared by all species, while 40.8% of the genes are unique to a single species, with again the largest number in A. nidulans, A. niger and A. terreus (Additional file 2: Table S3B). In general, the CAZyme distribution among the species follows their phylogenetic relationship. In total, this means that only 70 genes are shared by all species, while the number of unique genes differs strongly by species (Fig. 1).
Growth on plant biomass related substrates
Growth of the eight Aspergillus species was evaluated on 35 plant biomass related carbon sources (Fig. 2, full profiles are available at www.fung-growth.org). Two isolates per species were tested to check that the differences are species specific and not strain specific. The general growth speed differed between the strains of a species, but for most species no significant carbon source related differences were observed between the strains. An exception to this is A. niger CBS513.88 that grew poorly on all pure carbon sources and was shown to have an amino acid auxotrophy (unpublished data), which explains this phenotype. Apparently, both wheat bran (WB) and sugar beet pulp (SBP) contain sufficient protein/amino acids to supplement this deficiency. All other strains grew well on MM + glucose and glucose was therefore used as an internal reference to compare the strains, to avoid misleading differences caused by general differences in growth speed between the species. Growth on the other substrates relative to growth on glucose was then compared between the species.
Growth on pure cellulose was zero to very poor for all species. Most fungi had similar growth on glucose, maltose, starch and xylan, with the exception of A. nidulans, for which poor growth on maltose was observed for one strain. Growth on sucrose was similar to growth on glucose for nearly all strains, but strong differences were observed on inulin. Good growth on this substrate was observed for A. niger ATCC1015 and A. fischeri, while for all other species growth was reduced compared to sucrose. This does not correlate with the number of putative inulin-targeting genes in the genomes as A. fischeri has only two (Table 2), while A. niger has six just like A. terreus, which grows poorly on inulin. Good growth compared to the other species was observed for A. niger, A. nidulans, A. fumigatus and A. flavus on guar gum (galactomannan). While A. nidulans has the highest number of galactomannan-targeting genes (Table 2), A. niger has in fact the lowest number and A. terreus again grows poorly, even though it has the second highest number of galactomannan-targeting genes. Growth of A. clavatus was particularly poor on pectin which correlates well with its very low number of pectin-targeting genes (only 30 putative genes), which was less than half of the other species (Additional file 2: Table S2A, B).
All species grew well on wheat bran and also, with only A. clavatus having reduced growth, on sugar beet pulp. These substrates were therefore selected to analyze their enzymatic ability in more detail. Composition analysis (Table 1) showed that wheat bran contains mainly cellulose and (arabino)xylan, with xyloglucan and pectin as minor components. In contrast, sugar beet pulp contains mainly cellulose, xyloglucan and pectin, which explains the reduced growth of A. clavatus.
Enzyme profiles during growth on wheat bran and sugar beet pulp
A preliminary test in which eight plant biomass degrading enzymes were measured at day two, three and four of cultivation, demonstrated that at day three activities were maximal for all fungi (data not shown). This time point was therefore selected for the full enzymatic analysis. Nineteen extracellular, lignocellulose-active enzyme activities of the liquid cultures were measured (Additional file 1: Figures S1–S3). Comparison of these profiles demonstrated strong differences among the species, not only in the quantities of the activities, but also in the induction of specific enzymes. For instance, the highest activity levels for most enzymes of A. terreus were observed during growth on sugar beet pulp, while wheat bran resulted in higher levels of most enzymes for A. flavus (Additional file 1: Figure S1). When the individual activities were compared across the species, specific differences became noticeable. Wheat bran consists mainly of cellulose and arabinoxylan and the main regulator controlling degradation of these polysaccharides is XlnR, which is present in all Aspergilli . Endoxylanase and β-xylosidase were mainly produced on wheat bran, and levels were particularly high for A. niger (Additional file 1: Figure S1). Endoarabinanase, α-rhamnosidase, pectate lyase and endogalactanases, all related to pectin degradation, were mainly produced on sugar beet pulp, but rarely were all four activities produced by one species. Sugar beet pulp contains mainly cellulose and pectin and therefore pectinases and cellulases would be expected to be the main enzymes produced on this substrate, which is confirmed by our data.
Mass spectrometric analysis of the extracellular proteins confirmed the activity measurements with respect to the enzymes that were detected (Additional file 2: Table S4A–D). Figure 3 shows the presence of orthologous enzymes involved in the degradation of different polysaccharides in wheat bran and sugar beet pulp. This analysis demonstrates the high degree of diversity among the species in the production of orthologous enzymes. Only a few orthologous enzymes are produced by all or most species and in most cases they are produced on both wheat bran and sugar beet pulp (Fig. 3) although often with significantly different levels (Additional file 2: Table S4A–D). These data highlight the different enzymatic approaches used by the eight species to degrade plant biomass.
Correlation of CAZy profiles, taxonomy and enzyme activity of the eight Aspergilli
Figure 4 shows the correlation of the species for genome content, enzyme activity and production of individual enzymes. Correlating the number of genes per CAZy family demonstrated that with respect to genome content, closely related species (A. oryzae – A. flavus, A. fischeri – A. fumigatus) cluster together (Fig. 4a). This indicates that the evolution of their genome content related to plant biomass degradation follows the evolutionary history of the species. However, more distantly related species, A. nidulans and A. terreus, can display similar CAZy content. A possible explanation for this finding is that natural habitat exerts a stronger influence on genome evolution than phylogenetic relatedness.
No clear correlation was observed between the enzyme activities produced in response to complex substrates and evolutionary relatedness (Fig. 4b), possibly due to the range of non-plant substrates some species are known to consume (e.g. collagen for A. terreus, A. flavus and A. nidulans and insect larvae for A. clavatus), resulting in a varying biotope range and dependence on plant biomass. The composition of wheat bran and sugar beet pulp is different and they should elicit different activity profiles. For six of the tested species, the wheat bran and sugar beet pulp activity profiles diverge strongly. Unexpectedly, the sugar beet pulp and wheat bran activity profiles clustered together for A. flavus and A. oryzae. Two of the three tested A. niger strains (N402 and ATCC 1015) clustered together for both substrates, while the third (CBS 513.88) was strongly divergent in the enzyme activity profile. These results show that strains of the same species (CBS 513.88 and ATCC 1015) with near identical genomic content can use significantly different sets of enzymes to hydrolyze complex biomass. It should be noted that genome sequence analysis suggests that strains ATCC 1015 and N402 are likely descended from the same isolate (A. Tsang and co-workers, unpublished data), which explains the clustering of their activity profiles.
Correlation of the proteomics data did not follow the activity correlation (Fig. 4b, c), which can be explained by the production of non-orthologous enzymes for the same general activity by different species (Additional file 2: Table S4A–D). This adds an additional dimension to the highly divergent strategies of these Aspergilli. Considering the fairly similar genome content of these species, we conclude that the differences in their plant biomass degrading strategies are mainly at the regulatory level. More detailed studies into the regulation of orthologous CAZyme encoding genes in several species could reveal whether this is due to different sets of target genes of the main regulators or whether additional unknown regulators modulate the influence of the main regulators.
In this study, we compared eight Aspergilli with respect to plant polysaccharide degradation. The variations in CAZyme content between these species were relatively low as compared to previous studies in which a more diverse set of fungal species was compared [3–8]. This can be explained by the close phylogenetic relationships and/or by and the highly similar habitats of these Aspergilli, which would direct genome evolution in a similar direction. Human use of and/or interaction with the species differs markedly; with A. niger and A. oryzae being widely used industrial fungi, A. fumigatus one of the most significant opportunistic fungal human pathogens, and A. flavus a plant pathogen. However, all these species are common inhabitants of soil and stored agricultural products, and their spores are widespread in both indoor and outdoor environments. Although some of the sequenced strains are domesticated and not recent natural isolates, the comparison to a second strain that is a natural isolate showed that growth on 35 carbon sources is nearly identical for two strains of the same species. This demonstrates that the sequenced isolates have maintained their natural ability to use the tested carbon sources.
Hierarchical clustering of the plant polysaccharide degrading enzymes of these species demonstrated that in general the species with the most similar CAZome are also taxonomically close.
The number of unique genes (Additional file 2: Table S2A, B) per species also correlates well with the phylogenetic distance of the species. The lowest number was found for A. oryzae, A. flavus, A. fischeri and A. fumigatus. As the first two and the last two species, respectively, are closely related, their high similarity explains this low number. The more distant species (A. nidulans, A. terreus, A. niger, A. clavatus) have higher numbers of unique genes.
A high level of variation was detected in the enzyme activities of the tested species during growth on sugar beet pulp and wheat bran. While all species grew well on these substrates, with the exception of somewhat less growth of A. clavatus on sugar beet pulp, the enzyme profiles of the species showed strong differences. Wheat bran consists mainly of cellulose and arabinoxylan, while sugar beet pulp contains mainly cellulose and pectin. Enzymes able to degrade these different combinations of polysaccharides would therefore be expected to be prominent in the culture filtrate of all species grown on these substrates. Our study confirmed this as nearly all enzyme activities detected on wheat bran in all species are involved in xylan and cellulose degradation, while mainly pectinolytic and cellulolytic enzymes were detected on sugar beet pulp. The main regulator controlling the production of xylanolytic and cellulolytic enzymes in Aspergillus is XlnR, which has been studied in detail in A. niger, A. oryzae and A. nidulans [19–28]. XlnR activates the expression of xylanolytic and cellulolytic genes in response to the presence of xylan or xylose, the latter being the actual inducer. Indications for similar regulation have been reported for the other species [29–35], although the range of target genes can differ per species [36, 37]. Regulation studies in A. niger have previously demonstrated that pectinolytic genes are induced by galacturonic acid, rhamnose, polygalacturonic acid or pectin [38–40]. Differences in pectinolytic gene content between A. niger, A. nidulans and A. oryzae may be influenced by the pH of their natural habitat . An acidic pH favors pectin hydrolases, while a neutral to alkaline pH favors pectin lyases, supported by the finding that all fungal GH28 pectin hydrolases have activity optimum between pH 2 and pH 5, while pectin lyases have optimum between pH 7 and pH 10 (https://mycoclap.fungalgenomics.ca) . The pH of most samples was 7 except for A. nidulans on sugar beet pulp (pH 8), A. niger on wheat bran (pH 5.5) and sugar beet pulp (pH 4.5), and A. clavatus on sugar beet pulp (pH 6). The pH in the sugar beet pulp cultures correlates well with the pectin hydrolase and lyase activities and with the proteomics results (Additional file 2: Table S4A–D). Therefore, the differences in enzyme levels are likely caused by regulatory variation. Since the major regulators are shared by all tested species , their function or range of target genes in the tested Aspergilli is different and/or additional non-shared regulators are involved in the utilization of complex biomass. A difference in the function of the arabinanolytic regulator AraR in A. niger and A. nidulans was recently described , and the inducers for activation of AmyR also appear to differ between A. niger, A. nidulans and A. oryzae [44–47]. More detailed analysis of the set of target genes, and function and mechanism of the polysaccharide related regulators in the other species will be required to understand the mechanism responsible for these differences. Interestingly, the production of several cellulases appears to be conserved among the species, suggesting that this may be a core-activity for all species. In contrast, the production of hemicellulases is highly varied, suggesting specific adaptations of the species in their biomass degrading approach.
Laccase activity was detected for most species, with the highest activity on wheat bran for A. flavus and on sugar beet pulp for A. fumigatus. This does not correlate with the numbers of putative laccases detected in the genomes, suggesting significantly different regulation of the production of these enzymes among the species. Induction of laccase-encoding genes has been mainly studied in basidiomycetes in which transcription is modulated by metal ions (Cu2+, Ag+, Mn2+), aromatic compounds, nitrogen and carbon sources (nature and ratio) . In ascomycetes, regulatory elements such as heat shock elements, metal response elements and nitrogen metabolite regulation elements (NIT-2 like) were identified in the promoter region of laccase-like multicopper oxidase . In addition, laccases are also involved in other biological processes, such as spore pigment formation , and not only in lignin degradation, so the total number of laccases likely does not reflect the number of laccases which play a role in plant biomass degradation.
Although the fungi tested in this study produce diverse enzyme sets, they all grow well on the crude plant biomass substrates. This suggests that different strategies for the degradation of plant biomass may be equally efficient (as measured by fungal growth). In biotechnological applications, such as biofuel production, complete hydrolysis of the plant biomass is difficult to achieve with currently available enzyme cocktails. This may in part be explained by the absence of specific activities in these mixtures. The data obtained in this study show the existence of distinctly different enzymatic approaches to degrade biomass. A judicious mix of these approaches is likely to result in improved enzyme cocktails for biomass hydrolysis. Recently it was shown that addition of Podospora anserina hydrolases increases the efficiency of a Trichoderma reesei enzyme mixture . In this study we provide indications that similar results could be obtained with more closely related fungi. The advantage of using enzymes from other Aspergilli to improve enzyme cocktails of A. niger or A. oryzae is that heterologous production of these enzymes is not likely to cause problems due to the high similarity in gene structure of these species.
In this study we demonstrated that closely related fungi use highly diverse enzymatic strategies for the degradation of the same substrates, but with similar efficiency. It can be assumed that in nature not all fungi target the same parts of the substrate. Therefore the identification of the enzyme sets employed by different fungi can be used to design efficient commercial enzyme cocktails by combining these enzyme sets. This could significantly improve the saccharification efficiency of industrial enzyme cocktails. Combining the complete enzymes set of two fungi is likely to improve saccharification efficiency more than adding specific enzymes to the cocktail produced by a single species.
Media, growth conditions and chemicals
The fungal strains used in this study are listed in Additional file 1: Table S1. Aspergillus minimal medium was described previously . All monomeric and oligomeric carbon sources were added to a final concentration of 25 mM, while pure polymeric substrates and crude substrates were added to a final concentration of 1 and 3%, respectively. The pH of the medium was adjusted to 6.0. For plate growth, the center of the plates was inoculated with 2 μl of a suspension of 500 spores/μl and plates were incubated for 5 days at 30°C. All eight species were grown on minimal medium with 35 carbon sources including crude plant biomass, pure plant polysaccharides, oligosaccharides, monosaccharides and control substrates (casein, lignin) (Additional file 1: Figure S1, www.fung-growth.org). To confirm that the detected differences were species specific, a second isolate of each species was examined along with the sequenced strain. Growth on 25 mM glucose was used as a reference because the tested strains grow at different rates and glucose, among the monosaccharides, supported the fastest growth for all species. Growth on the other substrates relative to growth on glucose was then compared among the species. Growth on plates was analyzed by visual inspection by two authors independently after which these were compared and discussed.
Liquid cultures were inoculated with 106 spores/ml (final concentration) and incubated at 250 rpm for 3 days. All cultures were incubated at 30°C and performed in duplicate. Two to three strains of all species were grown in liquid cultures with 1% wheat bran or 1% sugar beet pulp. Culture filtrates after 3 days of cultivation were analyzed for the presence of free monomeric sugars, but no glucose, xylose, galacturonic acid, rhamnose or fructose was detected. SDS-PAGE analysis of the extracellular proteins revealed nearly identical profiles for strains of the same species (Additional file 1: Figure S6), indicating that enzyme production is highly conserved within a species. Detailed analysis of the produced enzymes was therefore only performed on a single strain.
Glucose, maltose, sucrose, inulin, beechwood xylan, Guar gum, apple pectin and all p-nitrophenyl-substrates were from Sigma–Aldrich. Soluble starch was from Difco. Red Debranched Arabinan (S-RDAR), Azo-CM-cellulose (S-ACMC), Azo-galactan (S-AGALP) and AZ-rhamnogalacturonan (S-AZRH), Azo-wheat arabinoxylan (S-AWAXP) and polygalacturonic acid (PGA) were from Megazyme International Ireland.
Composition analysis of plant biomass substrates
Sugar composition was determined by analyzing the sugars as their alditol acetate derivatives using GC-FLD as described previously .
The identification step of CAZymes followed the procedures previously described  where sequences are subject to BlastP analysis  against a library composed of modules derived from the CAZy database, the positive hits are then subjected to a modular annotation procedure that maps the individual modules onto the peptide using hits against libraries of catalytic and carbohydrate models derived from CAZy using BlastP or Hidden Markov models [54, 55]. The functional annotation step involves BlastP comparisons against a library of modules derived from biochemically characterized enzymes .
Orthology and synteny analysis
Genome scale protein ortholog clusters were constructed using OrthoMCL  by inflation factor 1, E value cutoff 1E−3, percentage match cutoff 60% as for identification of distant homologs . The orthologs clusters were further split according to the synteny detected by the Sybil algorithm  at www.aspgd.org. Sequences of genes were manually double checked by multiple sequence alignments with MAFFT  and potential errors of gene models were corrected.
All exo-acting CAZy enzyme activities were performed in microtiter plates. Reactions were carried out in 100 µl volumes containing 25 mM sodium acetate (pH 5), 0.01% substrate and suitably diluted culture filtrate. The mixture was incubated at 30°C for 2 h and the reaction was terminated by the addition of 100 µl 250 mM sodium carbonate. Enzyme activities (α-arabinofuranosidase, cellobiohydrolase, α-galactosidase, β-galactosidase, α-glucosidase, β-glucosidase, glucoamylase (α-maltosidase), β-mannosidase, α-rhamnosidase and β-xylosidase) were determined spectrophotometrically at 405 nm by measuring the release of p-nitrophenol (pNP) from their appropriate pNP-substrates and standardized against a known concentration of p-nitrophenol (pNP). Activities were expressed as nmol pNP/ml sample/min.
Endoarabinanase, endo-1,4-β-glucanase (cellulase), endo-1,4-β-galactanase and rhamnogalacturonanase activities were measured using 20 mg/ml of Red Debranched Arabinan (S-RDAR), Azo-CM-cellulose (S-ACMC), Azo-galactan (S-AGALP) and AZ-rhamnogalacturonan (S-AZRH), respectively. Endo-1,4-β-xylanase activity was measured using 10 mg/ml Azo-wheat arabinoxylan (S-AWAXP). 100 μl reactions were carried out containing equal volumes of buffered substrate (pH 4.5) and suitably diluted culture filtrate which were then incubated at 40°C for 1 h in the case of the endoarabinanase, endo-1,4-β-glucanase and endo-1,4-β-xylanase activities and 16 h for the endo-1,4-β-galactanase and rhamnogalacturonan activities. Endoarabinanase reactions were terminated with the addition of 400 μl 95% ethanol, endo-1,4-β-galactanase, rhamnogalacturonanase and endo-1,4-β-xylanase reactions with 250 μl 95% ethanol and endo-1,4-β-glucanase reactions with a 250 μl solution of sodium acetate trihydrate (40 mg/ml) and zinc acetate (4 mg/ml) in 76% ethanol. Precipitated reactions were then centrifuged at 1,000×g for 10 min and optical density of supernatants was measured at 590 nm. Endoarabinanase reactions were measured at 520 nm. Endo-acting enzyme activities are expressed as amount of dye released (absorbance change)/ml sample/min.
Pectate lyase activity was assayed using polygalacturonic acid (PGA). Reaction mixtures contained equal volumes of 50 mM N-cyclohexyl-3-aminopropanesulfonic acid (CAPS) (pH 10.0) and 2.5 mg/ml PGA, to which suitably diluted culture filtrate was added. Changes in absorbance at 235 nm were measured for approximately 30 min at 40°C.
Laccase activity was assayed using 2,2′-azino-di-(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS). Reaction mixtures contained 700 µl H2O, 100 µl 0.5 M glycine–HCl (pH 3.0), 100 µl culture filtrate and 100 µl 14 mM ABTS. The reaction was monitored by measuring the change in absorption at 436 nm at 30°C. The extinction coefficient of 29,300/M/cm was used for oxidized ABTS. Activity is expressed as is in nmol/min/ml.
Feruloyl esterase activities were determined spectrophotometrically (Shimadzu PharmaSpec UV-1700) at 37°C in 100 mM MOPS (3-(N-morpholino)propanesulfonic acid) buffer (pH 6.0). Methyl caffeate (MC), methyl ferulate (MF), methyl p-coumarate (MpC) and methyl sinapate (MS) (1.18 mM stock solutions in 100 mM MOPS, pH 6.0) were used as substrates. Reaction mixture contained 100 µl of culture liquid, 870 µl MOPS (3-(N-morpholino)propanesulfonic acid) buffer and the reaction was started by the addition of 30 µl substrate. Absorbance was monitored for 5 min at 308 nm for MpC (ε308 = 20,390/M/cm), 320 nm for MF (ε320 = 29,680/M/cm) and MS (ε320 = 15,890/M/cm), and 322 nm for MC (ε322 = 14,720/M/cm). FAE activities were expressed as nkat/l (10−9 mol/s/l).
Determination of monomeric sugars in the cultures
Presence of monomeric sugars in wheat bran and sugar beet pulp liquid cultures was measured by using Megazyme’s Assay kits for glucose and fructose (K-FRUGL), xylose (K-XYLOSE), glucuronic acid (K-URONIC) and rhamnose (K-RHAMNOSE) using the provider’s instruction. All measurements were done with two biological replicates.
Protein profiles were obtained by combining 25 μl of culture supernatant supplemented with 5 μl of 5 × Laemmli Loading Buffer (50 mM Tris–HCl pH 6.8, 2% SDS, 10% glycerol, 0.1 M dithiothreitol, 0.2 mg/ml Bromophenol Blue) and separating this on 12% SDS-PAGE gels. Proteins were visualized by silver staining and a PageRuler Unstained Protein Ladder (Thermo Scientific) was used as protein marker.
Proteins from 3 ml of culture filtrate were precipitated with cold TCA/acetone and the amount of protein recovered was determined using the RCDC kit assay (BioRad, Mississauga, ON, Canada). Five micrograms of protein were digested with trypsin and an aliquot analyzed by LC–MS/MS as previously described  on a Velos LTQ-Orbitrap mass spectrometer (Thermo-Fisher, San Jose, CA, USA). MS/MS data were processed using Proteome Discoverer Quant 1.3 (Thermo-Fisher) and spectral data were searched against Aspergillus protein databases downloaded from the Aspergillus Genome Database (AspGD). Search parameters used were 0.80 Da for fragment ion tolerance and 10.0 ppm for parent ion tolerance, fixed iodoacetamide cysteine modification and variable methionine oxidation. Protein and peptide identification confidence filters were applied to satisfy a 1% false discovery rate at the Peptide and Protein level. Protein grouping was applied so as to satisfy the principles of parsimony. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD000982.
Hierarchical clustering and correlation analysis
Matrix files of presence/absence and activity of CAZyme encoding genes and protein production measured by proteomics experiments were generated. The hierarchical clustering of CAZyme encoding genes were created using R  using the Euclidean distance with complete linkage and visualized by iTOL [62, 63]. The dendrogram and heatmaps of protein abundance were created and visualized using Genesis  with Pearson’s correlation and complete linkage.
Sequences of the RPB1, RPB2 (RNA polymerase II genes), Tsr1 (putative ribosome biogenesis protein), Cct8 (putative chaperonin complex component TCP-1) and AguA (α-glucuronidase) genes were downloaded from the full genome data sets and aligned using the Muscle software in the MEGA5 package . After aligning, the data sets were combined and maximum likelihood analysis was performed using RAxML version 7.2.8 . Each locus was treated as a separate partition. The number of bootstrap replicates was set on 1,000 replicates. Sequences of Penicillium chrysogenum Wisconsin 54–1,255 were used as outgroup.
plant polysaccharide degradation
Red Debranched Arabinan
sugar beet pulp
de Vries RP, Nayak V, van den Brink J, Vivas Duarte AD, Stalbrand H (2012) Fungal degradation of plant oligo- and polysaccharides. In: Grunwald P (ed) Carbohydrate modifying biocatalysts. Pan Stanford Publishing Pte Ltd., Singapore, pp 693–759
de Vries RP, Visser J (2001) Aspergillus enzymes involved in degradation of plant cell wall polysaccharides. Microb Mol Biol Rev 65:497–522
Amselem J, Cuomo CA, van Kan JA, Viaud M, Benito EP, Couloux A et al (2011) Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet 7:e1002230
Battaglia E, Benoit I, van den Brink J, Wiebenga A, Coutinho PM, Henrissat B et al (2011) Carbohydrate-active enzymes from the zygomycete fungus Rhizopus oryzae: a highly specialized approach to carbohydrate degradation depicted at genome level. BMC Genom 12:38
Berka RM, Grigoriev IV, Otillar R, Salamov A, Grimwood J, Reid I et al (2011) Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris. Nat Biotechnol 29:922–927
Coutinho PM, Andersen MR, Kolenova K, vanKuyk PA, Benoit I, Gruben BS et al (2009) Post-genomic insights into the plant polysaccharide degradation potential of Aspergillus nidulans and comparison to Aspergillus niger and Aspergillus oryzae. Fungal Genet Biol 46(Suppl 1):S161–S169
Eastwood DC, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A et al (2011) The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. Science 333:762–765
Espagne E, Lespinet O, Malagnac F, Da Silva C, Jaillon O, Porcel BM et al (2008) The genome sequence of the model ascomycete fungus Podospora anserina. Genome Biol 9:R77
Aro N, Pakula T, Penttila M (2005) Transcriptional regulation of plant cell wall degradation by filamentous fungi. FEMS Microbiol Rev 29:719–739
Baldrian P, Valaskova V (2008) Degradation of cellulose by basidiomycetous fungi. FEMS Microbiol Rev 32:501–521
Kowalczyk JE, Benoit I, de Vries RP (2014) Regulation of plant biomass utilization in Aspergillus. Adv Appl Microbiol 88:31–56
Kubicek CP, Starr TL, Glass NL (2014) Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. Ann Rev Phytopathol. 52:427–451
MacDonald J, Suzuki H, Master ER (2012) Expression and regulation of genes encoding lignocellulose-degrading activity in the genus Phanerochaete. Appl Microbiol Biotechnol 94:339–351
Stricker AR, Mach RL, de Graaff LH (2008) Regulation of transcription of cellulases- and hemicellulases-encoding genes in Aspergillus niger and Hypocrea jecorina (Trichoderma reesei). Appl Microbiol Biotechnol 78:211–220
Tani S, Kawaguchi T, Kobayashi T (2014) Complex regulation of hydrolytic enzyme genes for cellulosic biomass degradation in filamentous fungi. Appl Microbiol Biotechnol 98:4829–4837
Todd RB, Zhou M, Ohm RA, Leeggangers HA, Visser L, de Vries RP (2014) Prevalence of transcription factors in ascomycete and basidiomycete fungi. BMC Genom 15:214
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B (2014) The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42:D490–D495
Battaglia E, Visser L, Nijssen A, van Veluw J, Wösten HAB, de Vries RP (2011) Analysis of regulation of pentose utilisation in Aspergillus niger reveals evolutionary adaptations in the Eurotiales. Stud Mycol 69:31–38
de Vries RP, Visser J, de Graaff LH (1999) CreA modulates the XlnR-induced expression on xylose of Aspergillus niger genes involved in xylan degradation. Res Microbiol 150:281–285
Delmas S, Pullan ST, Gaddipati S, Kokolski M, Malla S, Blythe MJ et al (2012) Uncovering the genome-wide transcriptional responses of the filamentous fungus Aspergillus niger to lignocellulose using RNA sequencing. PLoS Genet 8:e1002875
Hasper AA, Trindade LM, van der Veen D, van Ooyen AJ, de Graaff LH (2004) Functional analysis of the transcriptional activator XlnR from Aspergillus niger. Microbiology 150:1367–1375
Marui J, Kitamoto N, Kato M, Kobayashi T, Tsukagoshi N (2002) Transcriptional activator, AoXlnR, mediates cellulose-inductive expression of the xylanolytic and cellulolytic genes in Aspergillus oryzae. FEBS Lett 528:279–282
Marui J, Tanaka A, Mimura S, de Graaff LH, Visser J, Kitamoto N et al (2002) A transcriptional activator, AoXlnR, controls the expression of genes encoding xylanolytic enzymes in Aspergillus oryzae. Fungal Genet Biol 35:157–169
Noguchi Y, Sano M, Kanamaru K, Ko T, Takeuchi M, Kato M et al (2009) Genes regulated by AoXlnR, the xylanolytic and cellulolytic transcriptional regulator, in Aspergillus oryzae. Appl Microbiol Biotechnol 85:141–154
Noguchi Y, Tanaka H, Kanamaru K, Kato M, Kobayashi T (2011) Xylose triggers reversible phosphorylation of XlnR, the fungal transcriptional activator of xylanolytic and cellulolytic genes in Aspergillus oryzae. Biosci Biotechnol Biochem 75:953–959
Tamayo EN, Villanueva A, Hasper AA, de Graaff LH, Ramon D, Orejas M (2008) CreA mediates repression of the regulatory gene xlnR which controls the production of xylanolytic enzymes in Aspergillus nidulans. Fungal Genet Biol 45:984–993
van Peij N, Gielkens MMC, de Vries RP, Visser J, de Graaff LH (1998) The transcriptional activator XlnR regulates both xylanolytic and endoglucanase gene expression in Aspergillus niger. Appl Environ Microbiol 64:3615–3619
van Peij NN, Visser J, de Graaff LH (1998) Isolation and analysis of xlnR, encoding a transcriptional activator co-ordinating xylanolytic expression in Aspergillus niger. Mol Microbiol 27:131–142
Bakri Y, Masson M, Thonart P (2010) Isolation and identification of two new fungal strains for xylanase production. Appl Biochem Biotechnol 162:1626–1634
de Souza CG, Girardo NS, Costa MA, Peralta RM (1999) Influence of growth conditions on the production of xylanolytic enzymes by Aspergillus flavus. J Basic Microbiol 39:155–160
Elshafei AM, Hassan MM, Haroun BM, Abdel-Fatah OM, Atta HM, Othman AM (2009) Purification and properties of an endoglucanase of Aspergillus terreus DSM 826. J Basic Microbiol 49:426–432
Gawande PV, Kamat MY (1999) Production of Aspergillus xylanase by lignocellulosic waste fermentation and its application. J Appl Microbiol 87:511–519
Mellon JE, Cotty PJ, Callicott KA, Abbas H (2011) Identification of a major xylanase from Aspergillus flavus as a 14-kD protein. Mycopathologia 172:299–305
Peixoto-Nogueira Sde C, Michelin M, Betini JH, Jorge JA, Terenzi HF, Polizeli Mde L (2009) Production of xylanase by Aspergilli using alternative carbon sources: application of the crude extract on cellulose pulp biobleaching. J Ind Microbiol Biotechnol 36:149–155
Tani S, Kanamasa S, Sumitani J-I, Arai M, Kawaguchi T (2012) XlnR-independent signaling pathway regulates both cellulase and xylanase genes in response to cellobiose in Aspergillus aculeatus. Curr Genet 58:93–104
Klaubauf S, Narang HM, Post H, Zhou M, Brunner K, Mach-Aigner AR et al (2014) Similar is not the same: Differences in the function of the (hemi-)cellulolytic regulator XlnR (Xlr1/Xyr1) in filamentous fungi. Fungal Genet Biol 72:73–81
Kunitake E, Tani S, Sumitani J, Kawaguchi T (2013) A novel transcriptional regulator, ClbR, controls the cellobiose- and cellulose-responsive induction of cellulase and xylanase genes regulated by two distinct signaling pathways in Aspergillus aculeatus. Appl Microbiol Biotechnol 97:2017–2028
de Vries RP, Jansen J, Aguilar G, Parenicová L, Benen JAE, Joosten V et al (2002) Expression profiling of pectinolytic genes from Aspergillus niger. FEBS Lett 530:41–47
Gruben BS, Zhou M, Wiebenga A, Ballering J, Overkamp KM, Punt PJ et al (2014) Aspergillus niger RhaR, a regulator involved in l-rhamnose release and catabolism. Appl Microbiol Biotechnol 98:5531–5540
Martens-Uzunova ES, Schaap PJ (2009) Assessment of the pectin degrading enzyme network of Aspergillus niger by functional genomics. Fungal Genet Biol 46(Suppl 1):S170–S179
Pel HJ, de Winde JH, Archer DB, Dyer PS, Hofmann G, Schaap PJ et al (2007) Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88. Nat Biotechnol 25:221–231
Murphy C, Powlowski J, Wu M, Butler G, Tsang A (2011) Curation of characterized glycoside hydrolases of fungal origin. Database J Biol Databases Curation 2011:bar020
Battaglia E, Hansen SF, Leendertse A, Madrid S, Mulder H, Nikolaev I et al (2011) Regulation of pentose utilisation by AraR, but not XlnR, differs in Aspergillus nidulans and Aspergillus niger. Appl Microbiol Biotechnol 91:387–397
Makita T, Katsuyama Y, Tani S, Suzuki H, Kato N, Todd RB et al (2009) Inducer-dependent nuclear localization of a Zn(II)(2)Cys(6) transcriptional activator, AmyR, in Aspergillus nidulans. Biosci Biotechnol Biochem 73:391–399
Murakoshi Y, Makita T, Kato M, Kobayashi T (2012) Comparison and characterization of alpha-amylase inducers in Aspergillus nidulans based on nuclear localization of AmyR. Appl Microbiol Biotechnol 94:1629–1635
Suzuki K, Tanaka M, Konno Y, Ichikawa T, Ichinose S, Hasegawa-Shiro S et al (2015) Distinct mechanism of activation of two transcription factors, AmyR and MalR, involved in amylolytic enzyme production in Aspergillus oryzae. Appl Microbiol Biotechnol 99:1805–1815
Vankuyk PA, Benen JA, Wosten HA, Visser J, de Vries RP (2012) A broader role for AmyR in Aspergillus niger: regulation of the utilisation of d-glucose or d-galactose containing oligo- and polysaccharides. Appl Microbiol Biotechnol 93:285–293
Piscitelli A, Giardina P, Lettera V, Pezzella C, Sannia G, Faraco V (2011) Induction and transcriptional regulation of laccases in fungi. Curr Genomics 12:104–112
Levasseur A, Saloheimo M, Navarro D, Andberg M, Pontarotti P, Kruus K et al (2010) Exploring laccase-like multicopper oxidase genes from the ascomycete Trichoderma reesei: a functional, phylogenetic and evolutionary study. BMC Biochem 11:32
Tsai HF, Wheeler MH, Chang YC, Kwon-Chung KJ (1999) A developmentally regulated gene cluster involved in conidial pigment biosynthesis in Aspergillus fumigatus. J Bacteriol 181:6469–6477
Couturier M, Haon M, Coutinho PM, Henrissat B, Lesage-Meessen L, Berrin JG (2010) Podospora anserina hemicellulases potentiate the Trichoderma reesei secretome for saccharification of lignocellulosic biomass. Appl Environ Microbiol 77:237–246
de Vries RP, Burgers K, van de Vondervoort PJI, Frisvad JC, Samson RA, Visser J (2004) A new black Aspergillus species, A. vadensis, is a promising host for homologous and heterologous protein production. Appl Environ Microbiol. 70:3954–3959
Kuhnel S, Schols HA, Gruppen H (2011) Aiming for the complete utilization of sugar-beet pulp: examination of the effects of mild acid and hydrothermal pretreatment followed by enzymatic digestion. Biotechnol Biofuels 4:14
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W et al (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389–3402
Eddy SR (1998) Profile hidden Markov models. Bioinformatics 14:755–763
Li L, Stoeckert CJ Jr, Roos DS (2003) OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 13:2178–2189
Boekhorst J, Snel B (2007) Identification of homologs in insignificant blast hits by exploiting extrinsic gene properties. BMC Bioinform 8:356
Crabtree J, Angiuoli SV, Wortman JR, White OR (2007) Sybil: methods and software for multiple genome comparison and visualization. Methods Mol Biol 408:93–108
Katoh K, Kuma K-I, Toh H, Miyata T (2005) MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res 33:511–518
Ozturkoglu Budak S, Zhou M, Brouwer C, Wiebenga A, Benoit I, Di Falco M et al (2014) A genomic survey of proteases in Aspergilli. BMC Genomics 15:523
Team RDC (2009) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna
Letunic I, Bork P (2007) Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23:127–128
Letunic I, Bork P (2011) Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res 39:W475–W478
Sturn A, Quackenbush J, Trajanoski Z (2002) Genesis: cluster analysis of microarray data. Bioinformatics 18:207–208
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28:2731–2739
Stamatakis A, Hoover P, Rougemont J (2008) A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol 57:758–771
IB, GAO, EB, BSG, EM, HMN, BTA and AW carried out the growth profiling studies. IB, HC, JvdB, KSH, LAJB, MRM, VM, BMC and HBOEB performed or supervised the enzyme assays. MZ, CPJMB, PAvK, PMC, AL, BH, AT and RPdV performed the bioinformatics studies. MDF and AT performed the proteomic studies. OB performed the sugar analysis. JH performed the phylogenetic analysis. IB, HC, MZ, AT and RPdV drafted the manuscript. RPdV conceived of, designed and coordinated the study. All authors read and approved the final manuscript.
IB, GAO, EB and BSG were supported by grants of the Dutch Technology Foundation STW, Applied Science division of NWO and the Technology Program of the Ministry of Economic Affairs UGC 07063 & 07938 to RPdV. IB was also supported by a grant of Kluyver Centre for Genomics of Industrial Fermentation KCII 2.2.27 to RPdV. GAO was also supported by a grants from DGAPA-UNAM number IN225710 and IN219813. JvdB was supported by a grant from the Netherlands Organisation for Scientific Research (NWO) JSTP.10.005 to RPdV. MZ was supported by a grant from the Netherlands Organisation for Scientific Research (NWO) and the Netherlands Genomics Initiative 93511035 to RPdV. The comparative genomics study was carried out on the Dutch national e-infrastructure with the support of SURF Foundation (e-infra130078). The proteomic work was supported by Genome Canada and Génome Québec.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Isabelle Benoit, Helena Culleton and Miaomiao Zhou contributed equally to this work
Combines Additional File Figures S1–S4 and Table S1. Figure S1: Hydrolytic enzyme activity profiles of the eight species. Figure S2: Laccase activity of the eight species. Figure S3: Differences in feruloyl esterase production. Figure S4: Conserved SDS-PAGE profiles for isolates of the same species. Table S1. Strains used in this study.
Combines Additional file Tables S2–S4 (excel format). Table S2A. Numbers of putative genes per CAZy family for the 10 genomes addressed in this study. Table S2B. Numbers of putative genes per Plant Polysaccharide Degradation-related CAZy family for the 10 genomes addressed in this study. Table S3A. Orthology clusters of feruloyl esterase (SF), glycoside hydrolase (GH), carbohydrate esterase and polysaccharide lyase (PL) families. Table S3B. Orthology clusters of auxiliary activities (AA). Table S4A. Detection of proteins in cultures grown on wheat bran sorted by CAZy family. Table S4B. Detection of proteins in cultures grown on sugar beet pulp sorted by CAZy family. Table S4C. Detected proteins in cultures grown on wheat bran sorted by number of species that contain an orthologue. Table S4D. Detected proteins in cultures grown on sugar beet pulp sorted by number of species that contain an orthologue.