Research article | Open | Published:
Carbohydrate-active enzymes from the zygomycete fungus Rhizopus oryzae: a highly specialized approach to carbohydrate degradation depicted at genome level
BMC Genomicsvolume 12, Article number: 38 (2011)
Rhizopus oryzae is a zygomycete filamentous fungus, well-known as a saprobe ubiquitous in soil and as a pathogenic/spoilage fungus, causing Rhizopus rot and mucomycoses.
Carbohydrate Active enzyme (CAZy) annotation of the R. oryzae identified, in contrast to other filamentous fungi, a low number of glycoside hydrolases (GHs) and a high number of glycosyl transferases (GTs) and carbohydrate esterases (CEs). A detailed analysis of CAZy families, supported by growth data, demonstrates highly specialized plant and fungal cell wall degrading abilities distinct from ascomycetes and basidiomycetes. The specific genomic and growth features for degradation of easily digestible plant cell wall mono- and polysaccharides (starch, galactomannan, unbranched pectin, hexose sugars), chitin, chitosan, β-1,3-glucan and fungal cell wall fractions suggest specific adaptations of R. oryzae to its environment.
CAZy analyses of the genome of the zygomycete fungus R. oryzae and comparison to ascomycetes and basidiomycete species revealed how evolution has shaped its genetic content with respect to carbohydrate degradation, after divergence from the Ascomycota and Basidiomycota.
The phylum Zygomycota is a primitive and early diverging group of fungi. Fungal species belonging to this basal fungal lineage are characterized by sexual reproduction via zygospores, asexual reproduction by uni-to multispored sporangia and, in most species, nonseptate (i.e coenocytic) hyphae. The Zygomycota is divided into two classes, the Zygomycetes and the Trichomycetes. The Zygomycetes are an ecologically diverse class of fungi, including both saprobes and pathogens of plants, animals (including humans) and other fungi. Rhizopus oryzae is a member of the order Mucorales, genus of Rhizopus. The order of Mucorales consists of genera such as Mucor, Phycomyces, Rhizopus, Rhizomucor. Although Rhizopus and Mucor are more closely related to each other than either of them to Phycomyces, they are not closely related genera . R. orzyae has been reported as the predominant human-pathogen causing zygomycosis, a highly destructive and lethal infection on immune-compromised hosts . In industry, it is one of the main Rhizopus spp. used in several traditional Asian fermented foods such as tempeh . Furthermore, it is the well-known producer of metabolites such as organic acids (e.g. lactic acid, fumaric acid), ethanol and hydrolytic enzymes (e.g. glucoamylases, polygalacturonases) from pentose sugars and agricultural wastes such as barley, cassava, corn, potato pulp, oats and rice [4–9].
R. oryzae has a world-wide distribution with a high prevalence in tropical and subtropical regions. It has been isolated from many substrates, including a wide variety of soils, decaying vegetation, fruit, vegetables, seeds and dung . As a fast growing filamentous fungus, it is known as a primary or secondary colonizer , invading quickly on easily accessible and digestible (e.g. rich in simple sugars) substrates. Recently the genome sequence of R. oryzae strain 99-880 has been published . The genome paper focussed mainly on the evolution of the genome, which was shown to include a whole genome duplication event. Expansions of gene families have been observed for specific cellular processes related to the human pathogenic lifestyle such as cell growth and virulence. Little attention was given, however, to the ability of R. oryzae to degrade plant and fungal polysaccharides.
Several recent studies have demonstrated a strong relationship between the repertoire of carbohydrate active enzymes (CAZymes, http://www.cazy.org)  in fungal genomes and their saprophytic lifestyle [14–17]. These studies focussed in particular on those CAZymes involved in polysaccharide degradation. However, for R. oryzae only two CAZy enzyme families, GH18 (chitinases)  and GH28 (polygalacturonases) , have been studied in detail. Here we report the CAZy analysis of the R. oryzae and link this annotation to the ability of this fungus to use plant and fungal polysaccharides as carbon source. This will not only provide a better understanding of the ecological role of the R. oryzae, but also of sugar consumption and degradation by zygomycete fungi in comparison to ascomycete and basidiomycete fungi.
Results and Discussion
CAZy annotation reveals different numbers of CAZy families for R. oryzae in comparison to ascomycete and basidiomycete fungi
A total of 317 CAZymes have been identified in the R. oryzae genome. This value is similar to the number of CAZymes found in known plant cell wall (PCW) degrading basidiomycetes, higher than for yeasts, but significantly lower than for filamentous ascomycetes (Table 1). A total of 116 glycoside hydrolases were found in the R. oryzae genome, which is lower to what has been found in other filamentous fungi. In contrast, the total number of 130 glycosyl transferases is a very high value by comparison to other fungi (Table 1). Average numbers of polysaccharide lyases (PL), carbohydrate esterases (CE) and carbohydrate binding modules (CBM) were identified in R. oryzae.
In the next sections a more detailed comparison between the putative CAZymes in R. oryzae and ascomycete and basidiomycete genomes related to the degradation of PCW polysaccharides (cellulose, xyloglucan, β-1,3-1,4-glucan, xylan, galactomannan, starch, inulin and pectin) is described (Table 2). This analysis will be complemented by that of enzyme activities that may be involved in the metabolism of endogenous or exogenous fungal cell wall (FCW) polysaccharides. Although displayed for comparative purposes, the numbers of CAZy modules of the yeast S. cerevisiae are not included in the discussion, as this organism lacks most of the enzymes involved in the degradation of PCW polysaccharides. As a foreword, the tentative assignment of putative enzymes to specific metabolic pathways found below is meant to reveal the overall ability to deal individually with different polysaccharides. A more rigorous accounting, if ever possible, would only be possible following experimental characterization of many gene products. It should be noted that the CAZy annotations are based upon the protein models derived from the annotation of the genomes. As the pipelines from the JGI and the Broad Institute are not identical, this may cause small differences in the numbers per family that affected the comparison.
CAZy annotation of the R. oryzae genome indicates highly specialized plant polysaccharide utilisation
Degradation of cellulose is achieved by the synergistic action of endoglucanases (endo-β-D-1,4-glucanases) and cellobiohydrolases, assisted by β-glucosidases. It should be noted that β-glucosidases also have other functions than to assist cellulose degradation and therefore not all β-glucosidases may be involved in cellulose breakdown. A total of 21 and 19 candidate cellulolytic enzymes are found in R. oryzae and in U. maydis genomes, respectively, lower than what is observed for the other species (Table 2). P. anserina had the largest set of candidate enzymes (72) related to cellulose degradation, as shown previously .
Fungal cellobiohydrolases are classified in families GH6 and GH7 and β-glucosidases are found in families GH1 and GH3. The R. oryzae genome contains no protein models belonging to families GH1, GH6, GH7. Moreover, only six GH3 genes were identified in R. oryzae compared while filamentous ascomycetes typically contain more than ten (Table 2). None of the GH3 members were predicted to be β-glucosidases (Additional File 1).
R. oryzae contains the largest number of putative enzymes assigned to family GH45 (5) whose characterised members are typically endoglucanases (Additional Files 1 and 2). Three out of these five GH45 ORFs contain a N-terminal CBM1 module (Additional File 1) and have previously been shown to encode functional endoglucanases (rce1, rce2 and rce3) [20, 21]. GH45 endoglucanases with a similar N-terminal CBM1 module have been identified in other members of the Mucorales , including MCE1 and MCE2 from Mucor circinelloides and PCE1 from Phycomyces nitens. Interestingly, the R. oryzae genome contains four putative protein models of GH9 (Supp. Table 2), a family of endoglucanases mostly found in bacteria, plants, and occasionally in animals . GH9 members are absent in all ascomycete filamentous fungi, while basidiomycetes generally encode a single GH9 protein harbouring a C-terminal membrane spanning region (Additional File 2). The role of fungal GH9s remains unclear, but is probably unrelated to cellulose degradation . Finally, R. oryzae does not appear to encode any family GH61 protein, which have been shown to boost cellulose breakdown by cellulases [27, 28].
Altogether the analysis suggests that the main cellulose-degrading enzymes produced by R. oryzae are GH45 endoglucanases and that the fungus uses them solely to access the plant material. Cellulose is therefore probably not a nutritional source for R. oryzae, as enzymes able to fully degrade cellulose into glucose have not been found.
The R. oryzae genome is poor in protein models related to xyloglucan degradation (3) (Table 2). Only U. maydis had a similarly low number, while the other basidiomycete and ascomycete fungi contain 8 to 17 protein models (Additional File 2). The genome appears to contain no xyloglucan-active β-1,4-D-endoglucanase (GH12, 74) nor any side-chain cleaving enzymes such as α-L-fucosidases of families GH29 and GH95 (Additional File 2). The 3 members of family GH31 are candidate α-glucosidases and not likely α-D-xylosidases (Additional File 1).
The R. oryzae genome also encodes only a relatively small number of enzymes for the breakdown of β-1,3-1,4-glucans (Table 2). 1,3-1,4-β-Glucans, which are abundant components of the cell walls of grasses , are typically hydrolysed by lichenases and endo-1,3-1,4-β-glucanases of family GH16. Eight GH16 protein models were identified in R. oryzae, a low number compared to ascomycete and basidiomycete species (Table 2).
Endoxylanases, commonly found in families GH10 and GH11, cleave the xylan backbone into smaller oligosaccharides, which are degraded further by β-xylosidases to xylose . The R. oryzae genome contains no putative ORFs belonging to either family GH10 or GH11. However, it has been demonstrated that R. oryzae CBS 112.07 produces a low molecular weight endoxylanase . Interestingly, family GH3 or GH43 candidate β-xylosidases were detected in the genome (Additional File 2). The two GH43 members were functionally annotated as candidates for arabinan endo-1, 5-α-L-arabinosidases (Additional File 1). R. oryzae appears to be unable to breakdown xylan significantly as it lacks essential xylanolytic enzymes.. With the exception of R. oryzae and the mycorrhizal basidiomycete L. bicolor (Figure 1, Additional File 2 and ), all other analyzed fungal species contain putative endoxylanases from families GH10 and GH11. The latter two fungi appear to share the inability to degrade xylan with the oomycete P. ultimum[17, 32].
The R. oryzae genome contains 13 candidate CAZymes related to the degradation of the main chain and side-chains of galactomannan. This value is similar to that of U. maydis, but significantly lower than that found in the other analyzed species (Table 2). No clear candidate β-endomannanase (EC 126.96.36.199) or β-mannosidase (EC 188.8.131.52) from families GH2, GH5 and GH26 were found in R. oryzae. Among the few GH5s found in R. oryzae, only one was found to be distantly related to bacterial β-mannosidases (Additional File 1). This constitutes presently the unique candidate enzyme for a possible role in the breakdown of the (gluco)mannan backbone. In A. niger, the presence of a single characterized β-mannosidase (MndA) and β-mannanase (Man5A) has been shown to be sufficient for an efficient degradation of mannan from different sources [33, 34].
The genome R. oryzae contains a total of six protein models in families GH27 and GH36, mostly described as α-galactosidases (Additional File 1). The enrichment of α-galactosidases in the genome suggests galactose release by R. oryzae, probably for nutritional use, from structural or storage galactoglucomannan, oligosaccharides such as melibiose, raffinose, stachyose, and/or xylan or xyloglucan.
Starch and inulin degradation
The R. oryzae genome contains an average number of candidate CAZymes involved in starch degradation compared to the ascomycetes and basidiomycetes (Table 2). R. oryzae contains the highest number of GH15 enzymes (6), whereas few members from families GH13 (4) and GH31 (3) were detected (Table 2). Four of the six GH15 members were annotated as candidate glucoamylases or related enzymes (Additional File 1). Glucoamylase (AmyA; RO3G_00082.3) has been characterized previously , in part due its common usage in the starch industry because of its high activity and stability . One of them has been previously characterized as a glucoamylase (AmyA; RO3G_00082.3)  and this enzyme has been commonly used in industry because of high activity and stability . AmyA consists of a C-terminal catalytic domain, connected to an N-terminal starch-binding module of CBM21 by an O-glycosylated linker . Interestingly, a very high number of CBM21-containing protein models (8) containing a CBM21 module were observed in R. oryzae compared to ascomycetes and basidiomycetes (1 - 4) (Additional File 1). Except for AmyA the actual function of these putative proteins is unknown. No catalytic modules were detected for the other seven ORFs, but all eight ORFs contain a putative phosphatase regulatory subunit.
Candidate starch-binding CBM20-containing protein models are found in all ascomycetes and basidiomycetes, excepting U. maydis. No CBM20 module is present in the R. oryzae genome, a feature probably functionally compensated by the eight CBM21-containing proteins. Interestingly, most CBM21-containing proteins are found among Eukaryota, contrasts with CBM20s that are equally dispersed among Bacteria and Eukaryota. As expected, none of the other CBM families containing starch-binding modules (CBM25, CBM26, CBM34 and CBM41) were found in the species analysed in this study. They are only found in prokaryotes (mostly Eubacteria) and the only one in Eukaryotes is in picophytotoplanktonic green algae .
No clear candidate invertases or inulinases could be detected in the genome of R. oryzae. In the R. orzyae NBRC 4785 strain, a glucoamylase (similar to R. oryzae NRRL 395 amyB) with sucrose-hydrolyzing activity has been however reported . AmyB has no starch-binding module and no apparent glucoamylase activity, but it was highly similar to the catalytic module of AmyA (glucoamylase). Only the equivalent AmyA mentioned earlier is found in the genome of R. oryzae strain 99-88 suggesting strain variations within family GH15. Interestingly, Aspergillus starch-degrading glucoamylases have been known for their promiscuity in degrading substrates containing α-linked D-glucosides  but no other reports of an ability to degrade sucrose are known.
In the R. oryzae genome, pectin degradation appears to be the focus of the highest number of putative CAZymes (27 predicted enzymes in Table 2). In comparison to the ascomycete and basidiomycete species, this is an average value. Less pectin-degrading CAZymes were predicted for U. maydis and L. bicolor, while the highest potential for pectin degradation is found in A. nidulans and A. oryzae (with an estimation of 94 and 109 CAZymes each, respectively). The R. oryzae genome is especially rich in family GH28 protein models, with twelve and three genes encoding candidate endopolygalacturonases (EC 184.108.40.206) and exopolygalacturonases (EC 220.127.116.11), respectively . Endopolygalacturonases particularly hydrolyse the α-1,4 linked D-galacturonic acid within the homogalacturonan chain and exo-polygalacturonases hydrolyse this chain from the nonreducing end . Our analysis shows that the enrichment of main-chain degrading enzymes of pectin was accompanied by low numbers of accessory enzymes (Table 2). A possible explanation is that R. oryzae is limited in the degradation of unbranched pectins, unlike most ascomycete and basidiomycete fungi that can degrade a variety of lowly and highly branched pectins.
This view fits well with the pre- and/or post harvest pathogenic lifestyle of R. oryzae. As a causal agent of Rhizopus rot, R. oryzae produces endopolygalacturonases for maceration of plant seedlings, fruits and vegetables and plant structures such as flowers, bulbs and tubers [42–44]. The pectin structure in berries, black currants and apples has been shown to contain a high amount of D-galacturonic acid and a small amount of rhamnose residues, indicating the presence of unbranched pectin [45, 46]. Rhizopus rot is characterized by the dissolution of the PCW middle lamella in fruits and vegetables, that are typically rich in unbranched pectin (homogalacturonans) with a low degree of methyl esterification . Such dissolution would facilitate accessibility into plant tissue. The natural ability of R. oryzae to degrade unbranched pectin has been exploited in industry, where GH28 endopolygalacturonases has been shown to be very useful in flax-retting (i.e. a process that degrades the smooth regions of pectin for separation of plant fibres) [48–50].
In summary, the genome contains mainly enzymes involved in degradation of storage polysaccharides (galactomannan and starch) and the backbone PCW structural polysaccharides (cellulose and pectin). The apparent absence of (i) exo-acting enzymes for cellulose, (ii) accessory enzymes for pectin degradation, and (iii) any enzymes for degradation of xyloglucan, β-1,3-1,4-glucan, xylan and inulin, suggest that these polysaccharides do not serve R. oryzae as a major carbon source. The degrading capacity of R. oryzae allows rapid invasion on and into the substrate, outcompeting secondary and late colonizers that will have to cope with a more extensive degradation of more resistant complex PCW polysaccharides.
The hydrolytic potential of the R. oryzae genome correlates with growth on polysaccharides and related monosaccharides
The CAZy annotation results were compared to the growth profile of two R. oryzae strains (the sequenced isolate 99-880 and the type strain of the species CBS 112.07) and a selection of species from the Fung-Growth database http://www.fung-growth.org on a number of carbon sources (Figure 1). Overall growth of the type isolate was better than the sequenced isolate. A strong preference for citrus pectin and guar gum (a galactomannan) was observed for R. oryzae 99-880, followed by inulin and starch. This pattern was largely the same for the type isolate, but this strain grew better on starch than the sequenced strain. Good growth on pectin, galactomannan and starch correlates with the presence of related CAZymes in the genome and the occurrence of R. oryzae on fruits, vegetables and seeds . Fruits and vegetables are rich in pectin, whereas seeds are enriched in storage polysaccharides such as galactomannan, starch and inulin. Both R. orzyae strains were able to grow on inulin, despite the absence of known genes required for its degradation. Genome annotation of the oomycete Pythium ultimum initially also identified no candidate invertase . An invertase gene was detected after Blast analysis of the raw nucleotide sequence, indicating that this was missed by the model-building pipeline. In light of this we used the same approach, but failed to identify any putative invertase in the sequenced strain of R. oryzae. Alternative explanations for the growth on inulin could be the presence of traces of other carbon sources in the commercial inulin preparation, or the presence of yet to discover fructofuranosidases.
No to poor growth (similar to growth on agar without carbon source) was observed on xylan (Figure 1) and cellulose (data not shown). Xylan was utilized by all other tested species with the exception of L. bicolor (Figure 1). This is in agreement with the absence of genes required to degradation of xylan in the genomes of R. oryzae and L. bicolor, and presence in all other tested genomes.
Growth on citrus pectin, guar gum, starch and inulin correlated well with growth on their associated monosaccharide components. R. oryzae showed average to good growth on D-galactose and D-mannose, the main components of galactomannan (Figure 1). The pectin main chain is mainly composed of D-galacturonic acid residues with small amounts of L-rhamnose residues, while the side chains consist of L-arabinose and D-galactose . This correlates well with the ability of R. oryzae to grow on D-galacturonic acid (Figure 1) and L-rhamnose (data not shown). No growth of P. chrysogenum, L. bicolor, P. anserina and M. grisea was observed on D-galacturonic acid (Figure 1). In general, fungi do not prefer D-galacturonic acid as a carbon source (de Vries et al, unpublished results). No to poor growth occurred on L-arabinose or D-ribose (data not shown), which strongly suggests that pentose sugars are not preferred by R. oryzae. Growth on D-xylose was observed for R. oryzae CBS 112.07, which confirms a previous study . However, no growth was observed for 99-880 on this substrate, demonstrating a significant difference between the isolates for xylose utilisation. Part of this difference may be due to a poor ability of some R. oryzae strains to germinate on D-xylose .
CAZy annotation reveals specific differences in the chitinolytic and glucanolytic enzyme system for R. oryzae
The following sections describe the comparison of the number of CAZy family members between R. oryzae and the ascomycete and basidiomycete species that are implicated in the turnover/modification of the fungal cell wall (FCW) polysaccharides chitin, chitosan, α-1,3-glucan, β-1,3-glucan and β-1,6-glucan (Table 3). R. oryzae appears to contain over 40 CAZymes for chitosan degradation/biosynthesis, more than any of the species analyzed so far. The chitin recycling machinery found in R. oryzae, is similar in size to that observed in filamentous ascomycete and basidiomycete species (Table 3). Significantly, low numbers of CAZymes for the degradation of α-1,3-glucan, β-1,3-glucan and β-1,6-glucan were found in R. oryzae when compared to ascomycetes and basidiomycetes.
Chitosan and chitin degradation
Chitin deacetylases catalyze the deacetylation of chitin to chitosan. In R. oryzae, a very large set of chitin deacetylases of family CE4 is found (Additional File 2). With a total of 34 CE4 models, this is more than three times the previous highest number of CE4s observed in the basidiomycete L. bicolor or in the ascomycete M. grisea. Most of the CE4 proteins in R. oryzae appear to be GPI-anchored (Additional File 1). The biosynthesis of cell wall chitosan, the main cell wall component of zygomycete species [52, 53], proceeds by coordinated action of both chitin synthase and chitin deacetylases . The many chitin deacetylases in the R. oryzae genome complement a large set of chitin synthases, previously reported as an expansion compared to ascomycete genomes . The physiological role of chitin deacetylases in R. oryzae is presently unknown but they have been suggested to play roles in FCW biosynthesis in Mucor rouxii and Absidia coerulea and in plant-pathogen interactions in Colletotrichum lindemuthianum and Aspergillus nidulans. Alternatively chitin deacetylases may also be involved in the digestion of chitin, whereby the resulting chitosan would be cleaved by chitosanases. It has been proposed that marine fungi use chitin deacetylases for the decomposition of chitin  according to the chitosan pathway .
Chitosanases (EC 18.104.22.168) have been described in families GH5, GH7, GH8, GH46, GH75, and GH80, the last three containing exclusively chitosanases. All presently characterized GH75 proteins are fungal chitosanases . Members from this family are present in all filamentous ascomycete species, but absent in R. oryzae, in the yeast S. cerevisiae and in all described basidiomycete fungi (Additional File 2). In R. oryzae, one and two proteins from families GH8 and GH46 were identified, respectively (Additional File 2). Our analysis of the corresponding sequences from R. oryzae indicates that they are only distantly related to bacterial chitosanases (Additional File 1). Bacterial chitosanases are known to play a role in the degradation and utilization of exogenous chitosan, whereas fungal chitosanases are assumed to recycle chitosan from the FCW. Production of chitosanases has been demonstrated during autolysis in Mucor rouxi.
Fungal chitinases from family GH18 are essential for FCW remodelling during growth and development . In R. oryzae, a total of 14 GH18 members was found, a value similar to what is found in known ascomycete and basidiomycete species (Additional File 2). All candidate proteins belonging to family GH18 are likely chitinases in R. orzyae, as previously suggested by phylogenetic analysis . The R. oryzae genome contains four protein models belonging to family GH20 and comparison to functionally characterized enzymes revealed that these members are related to β-N-acetylglucosaminidases (EC 22.214.171.124) or β-N-acetylhexosaminidases (EC 126.96.36.199). Family CE9 which includes N-acetylglucosamine-6-phosphate deacetylases (EC 188.8.131.52) is important for the metabolism of chitin. R. oryzae has only one family CE9 member, a number similar to ascomycete and basidiomycete fungi. Fungal chitinases have been shown to have diverse roles, such as remodelling of their own cell wall and release of nutrients . It remains unknown, whether these chitinolytic genes of R. oryzae also have a role in nutritional processes. Rhizopus and Mucor species have been shown to use chitinases for hyphal growth and autolysis [61, 62]. For other zygomycete fungi such as Mortierella spp, chitinolytic activity has been linked to utilisation of nutrients .
Both α-1,3-glucosidase and glucan endo-α-1,3-glucosidase are involved in hydrolysis of α-1, 3-glucan. α-1,3-Glucosidases are classified in families GH31 and GH63 and glucan endo-α-1,3-glucosidase in families GH71 and 87. R. oryzae encodes five proteins likely to breakdown α-1,3-glucan degradation, a lower value compared to the filamentous ascomycetes and basidiomycetes, excepting U. maydis. Functional analysis suggests that three members of family GH31 and two of family GH63 are candidate α-glucosidases, but it is unknown if these α-glucosidases could deal specifically with α-1,3-glycoside linkages. Both families are known to contain different enzyme activities, such as α-glucosidases that hydrolyse α-1, 4-linkages. No proteins were found in families GH71 and GH87, known for bearing glucan endo-α-1,3-glucosidases, unlike several ascomycete and basidiomycete genomes that contain GH71 members.
β-1,3-Glucan and β-1,6-glucan degradation
β-1,3-Glucanases can be divided into exo-β-1,3-glucanases and endo-β-1,3-glucanases. Four of the seven GH5 members from R. oryzae were functionally annotated as a candidate β-glucosidase related to exo-1,3-β-glucanases (Additional File 1), while five of the six protein models of family GH3 were functionally annotated as a candidate β-glucosidase or exo-1,3-β-glucosidase. R. oryzae also encodes two proteins belonging to family GH72, a family typically containing GPI-anchored β-1,3-glucanosyltransferases. These enzymes are known to play a central role in the cross linking of cell wall β-1,3-glucans to other cell wall β-glucans . The two R. oryzae GH72s are indeed related to characterized β-1,3-glucanosyltransglycosylases, and have a C-terminal GPI-anchor. The cell wall of zygomycetes consists of primarily of chitin and chitosan . However, there is evidence that the cell wall of zygomycetes could contain small amounts of β-1,3-glucan. β-1,3-Glucans have been suggested to participate in the regulation of cell wall morphology in the zygomycete M. rouxii. We have identified three putative 1,3-β-D-glucan synthases (family GT48) encoded by the R. oryzae genome. The likely presence of β-1,3-glucan might explain the recently detected susceptibility of R. oryzae to caspofungin, an antifungal agent which inhibits the 1,3-β-D-glucan synthase . The presence of reduced amounts of β-1,3-glucans in the cell wall does not correlate with the significant numbers of 1,3-β-glucanases in the genome. This could indicate a possible role for the candidate β-1,3-glucanases in the degradation of cell walls of ascomycetes and basidiomycetes, whose cell walls are enriched in chitin and β-1,3-glucans [52, 66, 67].
β-1,6-Glucanases may be found in families GH5 and GH30. The R. orzyae genome contains seven protein models belonging to these two GH families (Table 3), but comparison to biochemically characterized enzymes allowed no clear identification of candidate β-1,6-glucanases in R. orzyae (Additional File 1). This agrees well with the lack of reports on the presence of β-1,6-glucan in the cell wall of R. oryzae.
The hydrolytic potential of R. oryzae genome correlates with growth on fungal cell wall fractions and polysaccharides
The CAZy annotation results are consistent with growth data of R. oryzae on fungal cell wall polysaccharides and fractions thereof (Figure 2). Growth of R. oryzae on β-1,3-glucan, chitin and chitosan was comparable to the colony size and density on D-glucose. Significantly improved growth was observed on fungal cell wall fractions of S. cerevisiae, A. niger and A. bisporus (aggregated mycelium) (Figure 2), but no growth was observed on cell wall fractions of R. oryzae and P. blakesleeanes (data not shown). The zygomycete M. circinelloides showed a similar growth phenotype on cell wall polysaccharides and fractions. In contrast, poor growth of A. niger was observed on β-1,3-glucan, chitin, chitosan and cell wall fractions. All three species showed good growth on N-acetylglucosamine (Figure 2).
The chitinolytic and β-1,3-glucanolytic enzyme system of R. oryzae could have a role in either nutrition or fungal interactions. Fungal interactions may include inhibition of growth of coexisting fungi, a potential defence system for mycoparasitic competitors such as ascomycete Trichoderma spp. that are known to parasitize a range of fungi . Degradation of chitin, chitosan and/or β-1,3 glucan in the environment can serves as a carbon and nitrogen source for fungi. Recycling of cell wall material of former fungal colonizers has been observed for primary and secondary basidiomycete species of wood during succession . Growth of R. oryzae on chitin and chitosan suggest the ability to digest these polymers in the environment. Together with previous reported locations of isolation, e.g. the soil of salt-marshes , our observations are compatible with an adaptation of R. orzyae to inhabit coastal environments.
The R. oryzae genome encodes a repertoire of CAZymes distinct of that of ascomycete and basidiomycete fungi, suggesting strong evolutionally differences and adaptation to the environment at the genomic level. In contrast to most ascomycetes and basidiomycetes, the set of CAZymes in the R. oryzae genome - supported by the growth profile - reveals the ability to use easily digestible sugars accompanied by an inability to degrade complex plant cell wall polysaccharides. The chitinolytic and glucanolytic enzyme system and the ability to degrade chitin, chitosan, β-1,3-glucan and ascomycete and basidiomycete cell walls suggest a role in both morphological processes and non-plant based nutrition. The differences in hydrolytic potential identified for R. oryzae suggest the development of unique strategies for polysaccharide degradation in the Zygomycota when compared to Basidiomycota and Ascomycota.
Strains and growth conditions
The fungal strains and culture conditions are listed in Suppl. Table 3. Carbon sources were added to the minimal medium (MM) containing 1,5% agar (Merck,101614) at the following concentrations: 1% (w/v) for soluble starch, guar gum, citrus pectin, inulin and beechwood xylan and 25 mM for D-glucose, D-mannose, D-galactose, D-galacturonic acid, D-fructose and D-xylose (Additional File 3). The pH of the medium was adjusted to 6.0. The plates were inoculated with 2 μl of 500 spores/μl or a mycelium plug (1 mm diameter). For growth comparison on fungal cell wall fractions and polysaccharides, of Rhizopus oryzae, Mucor circinelloides and Aspergillus niger were grown on plates containing Aspergillus minimal medium  and 1,5% agar with 25 mM glucose, or 0.5% (w/v) N-acetylglucosamine, chitin, chitosan, β-1,3-glucan (Additional File 3), or 0.5% (w/v) dried cell wall fractions of Saccharomyces cerevisiae, Aspergillus niger, Agaricus bisporus. Plates were inoculated with 5 μl of 500 A. niger spores/μl, or with a plug of R. oryzae or M. circinelloides mycelium. Cultures were grown for 3 days at 25°C.
All 17,459 protein-encoding ORFs from the R. oryzae RA 99-880 genome were submitted to analysis using the CAZy annotation pipeline in a two-step procedure of identification and annotation. The identification step of CAZymes follows the procedures previously described  where sequences are subject to Blastp analysis  against a library composed of modules derived from CAZy, the positive hits are then subject to a modular annotation procedure that maps the individual modules against on the peptide using hits against libraries of catalytic and carbohydrate models derived from CAZy using BlastP or profile Hidden Markov models [71, 72]. The results are augmented with signal peptide, transmembrane, and GPI predictions by human curators [73–75]. The fragmentary models and all models suspected of splicing prediction errors are identified. The functional annotation step involves BlastP comparisons against a library modules derived from biochemically characterized enzymes . This manual comparison to identified functions yields three levels of detailed annotation: i) candidate activity, reflects a non-ambiguous assignment thanks to high levels of similarity and/or a large functional homogeneity of the hits; ii) related to activity, reflects a more distant similarity and/or a functional variability of the results; iii) distantly related to activity, reflects a weak degree of similarity to known activities. In most cases a generic activity (ex: α-glycosidase) was equally assigned.
Fungal cell wall extractions
Cell wall fractions were isolated as described by Damveld et al..
Budziszewska J, Wilk M, Wrzosek M: Taxonomic revision of the genus Rhizomucor. IMC9: The Biology of Fungi: 2010; Edinburgh, UK. 2010
Ribes JA, Vanover-Sams CL, Baker DJ: Zygomycetes in human disease. Clin Microbiol Rev. 2000, 13 (2): 236-301. 10.1128/CMR.13.2.236-301.2000.
Nout MJR, Rombouts FM: Recent developments in tempe research. J Appl Bacteriol. 1990, 69: 609-633.
Takahashi T, Tsuchida Y, Irie M: Purification and some properties of three forms of glucoamylase from a Rhizopus species. J Biochem. 1978, 84 (5): 1183-1194.
Hang R-cYaYD: Kinetics of direct fermentation of agricultural commodities to L(+) lactic acid by Rhizopus oryzae. Biotechnology Letters. 1989, 11 (8): 597-600. 10.1007/BF01040043.
Carta FS, Soccol CR, Ramos LP, Fontana JD: Production of fumaric acid by fermentation of enzymatic hydrolysates derived from cassava bagasse. Bioresource Technology. 1999, 68 (1): 23-28. 10.1016/S0960-8524(98)00074-1.
Saito K, Kawamura Y, Oda Y: Role of the pectinolytic enzyme in the lactic acid fermentation of potato pulp by Rhizopus oryzae. J Ind Microbiol Biotechnol. 2003, 30 (7): 440-444. 10.1007/s10295-003-0071-z.
Maas RH, Bakker RR, Eggink G, Weusthuis RA: Lactic acid production from xylose by the fungus Rhizopus oryzae. Appl Microbiol Biotechnol. 2006, 72 (5): 861-868. 10.1007/s00253-006-0379-5.
Thongchul N, Navankasattusas S, Yang ST: Production of lactic acid and ethanol by Rhizopus oryzae integrated with cassava pulp hydrolysis. Bioprocess Biosyst Eng. 2010, 33 (3): 407-416. 10.1007/s00449-009-0341-x.
Domsch KH, Gams W, Anderson TH: Compendium of Soil Fungi. Academic Press, London. 1980, 1: 703-709.
Richardson M: The ecology of the Zygomycetes and its impact on environmental exposure. Clin Microbiol Infect. 2009, 15 (Suppl 5): 2-9. 10.1111/j.1469-0691.2009.02972.x.
Ma LJ, Ibrahim AS, Skory C, Grabherr MG, Burger G, Butler M, Elias M, Idnurm A, Lang BF, Sone T, et al: Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication. PLoS Genet. 2009, 5 (7): e1000549-10.1371/journal.pgen.1000549.
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 2009, D233-238. 10.1093/nar/gkn663. 37 Database
Coutinho PM, Andersen MR, Kolenova K, vanKuyk PA, Benoit I, Gruben BS, Trejo-Aguilar B, Visser H, van Solingen P, Pakula T, et al: Post-genomic insights into the plant polysaccharide degradation potential of Aspergillus nidulans and comparison to Aspergillus niger and Aspergillus oryzae. Fungal Genet Biol. 2009, 46 (Suppl 1): S161-S169. 10.1016/j.fgb.2008.07.020.
Espagne E, Lespinet O, Malagnac F, Da Silva C, Jaillon O, Porcel BM, Couloux A, Aury JM, Segurens B, Poulain J, et al: The genome sequence of the model ascomycete fungus Podospora anserina. Genome Biol. 2008, 9 (5): R77-10.1186/gb-2008-9-5-r77.
Martinez D, Challacombe J, Morgenstern I, Hibbett D, Schmoll M, Kubicek CP, Ferreira P, Ruiz-Duenas FJ, Martinez AT, Kersten P, 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 (6): 1954-1959. 10.1073/pnas.0809575106.
Martin F, Aerts A, Ahren D, Brun A, Danchin EG, Duchaussoy F, Gibon J, Kohler A, Lindquist E, Pereda V, et al: The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature. 2008, 452 (7183): 88-92. 10.1038/nature06556.
Karlsson M, Stenlid J: Comparative Evolutionary Histories of the Fungal Chitinase Gene Family Reveal Non-Random Size Expansions and Contractions due to Adaptive Natural Selection. Evol Bioinform Online. 2008, 4: 47-60.
Mertens JA, Burdick RC, Rooney AP: Identification, biochemical characterization, and evolution of the Rhizopus oryzae 99-880 polygalacturonase gene family. Fungal Genet Biol. 2008, 45 (12): 1616-1624. 10.1016/j.fgb.2008.09.009.
Moriya T, Murashima K, Nakane A, Yanai K, Sumida N, Koga J, Murakami T, Kono T: Molecular cloning of endo-beta-D-1,4-glucanase genes, rce1, rce2, and rce3, from Rhizopus oryzae. J Bacteriol. 2003, 185 (5): 1749-1756. 10.1128/JB.185.5.1749-1756.2003.
Murashima K, Nishimura T, Nakamura Y, Koga J, Moriya T, Sumida N, Yaguchi T, Kono T: Purification and characterization of new endo-1,4-β-D-glucanases from Rhizopus oryzae. 2002, 30 (3): 319-326.
Shimonaka A, Koga J, Baba Y, Nishimura T, Murashima K, Kubota H, Kono T: Specific characteristics of family 45 endoglucanases from Mucorales in the use of textiles and laundry. Biosci Biotechnol Biochem. 2006, 70 (4): 1013-1016. 10.1271/bbb.70.1013.
Baba Y, Shimonaka A, Koga J, Kubota H, Kono T: Alternative splicing produces two endoglucanases with one or two carbohydrate-binding modules in Mucor circinelloides. J Bacteriol. 2005, 187 (9): 3045-3051. 10.1128/JB.187.9.3045-3051.2005.
Shimonaka A, Baba Y, Koga J, Nakane A, Kubota H, Kono T: Molecular cloning of a gene encoding endo-beta-D-1,4-glucanase PCE1 from Phycomyces nitens. Biosci Biotechnol Biochem. 2004, 68 (11): 2299-2305. 10.1271/bbb.68.2299.
Ni J, Takehara M, Watanabe H: Identification of activity related amino acid mutations of a GH9 termite cellulase. Bioresour Technol. 2010, 101 (16): 6438-6443. 10.1016/j.biortech.2010.03.045.
Wymelenberg AV, Denman S, Dietrich D, Bassett J, Yu X, Atalla R, Predki P, Rudsander U, Teeri TT, Cullen D: Transcript analysis of genes encoding a family 61 endoglucanase and a putative membrane-anchored family 9 glycosyl hydrolase from Phanerochaete chrysosporium. Appl Environ Microbiol. 2002, 68 (11): 5765-5768. 10.1128/AEM.68.11.5765-5768.2002.
Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, et al: Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry. 2010, 49 (15): 3305-3316. 10.1021/bi100009p.
Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sorlie M, Eijsink VG: An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010, 330 (6001): 219-222. 10.1126/science.1192231.
Burton RA, Fincher GB: (1,3;1,4)-beta-D-glucans in cell walls of the poaceae, lower plants, and fungi: a tale of two linkages. Mol Plant. 2009, 2 (5): 873-882. 10.1093/mp/ssp063.
de Vries RP, Visser J: Aspergillus enzymes involved in degradation of plant cell wall polysaccharides. Microbiol Mol Biol Rev. 2001, 65 (4): 497-522. 10.1128/MMBR.65.4.497-522.2001. table of contents.
Bakir U, Yavascaoglub S, Guvencb F, Ersayinb A: An endo-β-1,4-xylanase from Rhizopus oryzae: production, partial purification and biochemical characterization. Enzyme and Microbial Technology. 2001, 29 (6-7): 328-334. 10.1016/S0141-0229(01)00379-9.
Levesque CA, Brouwer H, Cano L, Hamilton JP, Holt C, Huitema E, Raffaele S, Robideau GP, Thines M, Win J, et al: Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire. Genome Biol. 2010, 11 (7): R73-10.1186/gb-2010-11-7-r73.
Ademark P, de Vries RP, Hagglund P, Stalbrand H, Visser J: Cloning and characterization of Aspergillus niger genes encoding an alpha-galactosidase and a beta-mannosidase involved in galactomannan degradation. Eur J Biochem. 2001, 268 (10): 2982-2990. 10.1046/j.1432-1327.2001.02188.x.
Do BC, Dang TT, Berrin JG, Haltrich D, To KA, Sigoillot JC, Yamabhai M: Cloning, expression in Pichia pastoris, and characterization of a thermostable GH5 mannan endo-1,4-beta-mannosidase from Aspergillus niger BK01. Microb Cell Fact. 2009, 8: 59-10.1186/1475-2859-8-59.
Ashikari T, Nakamura N, Tanaka Y, Kiuchi N, ShiBano Y, Tanaka T, Amachi T, Yoshizumi H: Rhizopus raw-starch-degrading glucoamylase: Its cloning and expression in yeast. Agric Biol Chem. 1986, 50: 957-964.
Norouzian D, Akbarzadeh A, Scharer JM, Moo Young M: Fungal glucoamylases. Biotechnol Adv. 2006, 24 (1): 80-85. 10.1016/j.biotechadv.2005.06.003.
Chou WI, Pai TW, Liu SH, Hsiung BK, Chang MD: The family 21 carbohydrate-binding module of glucoamylase from Rhizopus oryzae consists of two sites playing distinct roles in ligand binding. Biochem J. 2006, 396 (3): 469-477. 10.1042/BJ20051982.
Ral JP, Derelle E, Ferraz C, Wattebled F, Farinas B, Corellou F, Buleon A, Slomianny MC, Delvalle D, d'Hulst C, et al: Starch division and partitioning. A mechanism for granule propagation and maintenance in the picophytoplanktonic green alga Ostreococcus tauri. Plant Physiol. 2004, 136 (2): 3333-3340. 10.1104/pp.104.044131.
Mertens J, Skory C: Isolation and characterization of a second glucoamylase gene without a starch binding domain from Rhizopus oryzae. Enzyme Microbiol Tech. 2007, 40: 874-880. 10.1016/j.enzmictec.2006.07.003.
Watanabe T, Oda Y: Comparison of sucrose-hydrolyzing enzymes produced by Rhizopus oryzae and Amylomyces rouxii. Biosci Biotechnol Biochem. 2008, 72 (12): 3167-3173. 10.1271/bbb.80346.
Meagher MM, Reilly PJ: Kinetics of the hydrolysis of di- and trisaccharides with Aspergillus niger glucoamylases I and II. Biotechnol Bioeng. 1989, 34 (5): 689-693. 10.1002/bit.260340513.
Yoshida S, Tsuyumu S, Tsukiboshi T: Macerating Enzymes Produced by Rhizopus oryzae in Infected Mulberry Roots. J Phytopathology. 2003, 151: 436-441. 10.1046/j.1439-0434.2003.00747.x.
Yoshida S, Murakam R, Watanabe T, Koyama A: Rhizopus Rot of Mulberry-grafted Saplings Caused by Rhizopus oryzae. J Gen Plant Pathol. 2001, 67: 291-293. 10.1007/PL00013033.
Shtienberg D: Rhizopus head rot of confectionery sunflower: effects on yield quantity and quality and implications for disease management. Phytopathology. 1997, 87 (12): 1226-1232. 10.1094/PHYTO.19184.108.40.2066.
de Vries JA, Rombouts FM, Voragen AGJ, Pilnik W: Enzymic degradation of apple pectins. Carbohydrate Polymers. 1982, 2: 25-33. 10.1016/0144-8617(82)90043-1.
Hilz H, Bakx EJ, Schols HA, Voragen GJ: Cell wall polysaccharides in black currants and bilberries - characterisation in berries, juice, and press cake. Carbohydrate Polymers. 2005, 59 (4): 477-488. 10.1016/j.carbpol.2004.11.002.
Steele NM, McCann MC, Roberts K: Pectin Modification in Cell Walls of Ripening Tomatoes Occurs in Distinct Domains. Plant Physiol. 1997, 114 (1): 373-381.
Xiao Z, Wang S, Bergeron H, Zhang J, Lau PC: A flax-retting endopolygalacturonase-encoding gene from Rhizopus oryzae. Antonie Van Leeuwenhoek. 2008, 94 (4): 563-571. 10.1007/s10482-008-9274-7.
Zhang J, Henriksson H, Szabo IJ, Henriksson G, Johansson G: The active component in the flax-retting system of the zygomycete Rhizopus oryzae sb is a family 28 polygalacturonase. J Ind Microbiol Biotechnol. 2005, 32 (10): 431-438. 10.1007/s10295-005-0014-y.
Zhang J, Henriksson G, Johansson G: Polygalacturonase is the key component in enzymatic retting of flax. J Biotechnol. 2000, 81 (1): 85-89. 10.1016/S0168-1656(00)00286-8.
Maas RH, Springer J, Eggink G, Weusthuis RA: Xylose metabolism in the fungus Rhizopus oryzae: effect of growth and respiration on L+-lactic acid production. J Ind Microbiol Biotechnol. 2008, 35 (6): 569-578. 10.1007/s10295-008-0318-9.
Bartnicki-Garcia S: Cell wall chemistry, morphogenesis, and taxonomy of fungi. Annu Rev Microbiol. 1968, 22: 87-108. 10.1146/annurev.mi.22.100168.000511.
Zhao Y, Park RD, Muzzarelli RA: Chitin deacetylases: properties and applications. Mar Drugs. 2010, 8 (1): 24-46. 10.3390/md8010024.
Davis L, Bartinicki-Garcia S: Chitosan synthesis by the tandem action of chitin synthetase and chitin deacetylase from M. rouxii. Biochemistry. 1984, 23: 1065-1073. 10.1021/bi00301a005.
Tsigos I, Martinou A, Kafetzopoulos D, Bouriotis V: Chitin deacetylases: new, versatile tools in biotechnology. Trends Biotechnol. 2000, 18 (7): 305-312. 10.1016/S0167-7799(00)01462-1.
Gooday GW: Physiology of microbial degradation of chitin and chitosan. Biodegradation. 1990, 1: 177-190. 10.1007/BF00058835.
Cheng CY, Chang CH, Wu YJ, Li YK: Exploration of glycosyl hydrolase family 75, a chitosanase from Aspergillus fumigatus. J Biol Chem. 2006, 281 (6): 3137-3144. 10.1074/jbc.M512506200.
Somashekar D, Joseph R: Chitosanases -properties and applications: a review. Bioresource Technology. 1996, 55: 35-45. 10.1016/0960-8524(95)00144-1.
Seidl V: Chitinases of filamentous fungi: a large group of diverse proteins with multiple physiological functions. Fungal Biology Reviews. 2008, 22 (1): 36-42. 10.1016/j.fbr.2008.03.002.
Gortari MC, Hours RA: Fungal chitinases and their biological role in the antagonism onto nematode eggs. A review. Mycol Progress. 2008, 7: 221-238. 10.1007/s11557-008-0571-3.
Rast DM, Horsch M, Furter R, Gooday GW: A complex chitinolytic system in exponentially growing mycelium of Mucor rouxii: properties and function. J Gen Microbiol. 1991, 137 (12): 2797-2810.
Yanai K, Takaya N, Kojima N, Horiuchi H, Ohta A, Takagi M: Purification of two chitinases from Rhizopus oligosporus and isolation and sequencing of the encoding genes. J Bacteriol. 1992, 174 (22): 7398-7406.
Hartland RP, Fontaine T, Debeaupuis JP, Simenel C, Delepierre M, Latge JP: A novel beta-(1-3)-glucanosyltransferase from the cell wall of Aspergillus fumigatus. J Biol Chem. 1996, 271 (43): 26843-26849. 10.1074/jbc.271.43.26843.
Argimon S, Galello F, Pereyra E, Rossi S, Moreno S: Mucor rouxii Rho1 protein; characterization and possible role in polarized growth. Antonie Van Leeuwenhoek. 2007, 91 (3): 237-251. 10.1007/s10482-006-9113-7.
Ibrahim AS, Bowman JC, Avanessian V, Brown K, Spellberg B, Edwards JE, Douglas CM: Caspofungin inhibits Rhizopus oryzae 1,3-beta-D-glucan synthase, lowers burden in brain measured by quantitative PCR, and improves survival at a low but not a high dose during murine disseminated zygomycosis. Antimicrob Agents Chemother. 2005, 49 (2): 721-727. 10.1128/AAC.49.2.721-727.2005.
Klis FM, Boorsma A, De Groot PW: Cell wall construction in Saccharomyces cerevisiae. Yeast. 2006, 23 (3): 185-202. 10.1002/yea.1349.
Latge JP: The cell wall: a carbohydrate armour for the fungal cell. Mol Microbiol. 2007, 66 (2): 279-290. 10.1111/j.1365-2958.2007.05872.x.
Harman GE, Howell CR, Viterbo A, Chet I, Lorito M: Trichoderma species--opportunistic, avirulent plant symbionts. Nat Rev Microbiol. 2004, 2 (1): 43-56. 10.1038/nrmicro797.
Lindahl BD, Finlay RD: Activities of chitinolytic enzymes during primary and secondary colonization of wood by basidiomycetous fungi. New Phytol. 2006, 169 (2): 389-397. 10.1111/j.1469-8137.2005.01581.x.
de Vries RP, Frisvad JC, van de Vondervoort PJ, Burgers K, Kuijpers AF, Samson RA, Visser J: Aspergillus vadensis, a new species of the group of black Aspergilli. Antonie Van Leeuwenhoek. 2005, 87 (3): 195-203. 10.1007/s10482-004-3194-y.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14 (9): 755-763. 10.1093/bioinformatics/14.9.755.
Eisenhaber B, Schneider G, Wildpaner M, Eisenhaber F: A sensitive predictor for potential GPI lipid modification sites in fungal protein sequences and its application to genome-wide studies for Aspergillus nidulans, Candida albicans, Neurospora crassa, Saccharomyces cerevisiae and Schizosaccharomyces pombe. J Mol Biol. 2004, 337 (2): 243-253. 10.1016/j.jmb.2004.01.025.
Fankhauser N, Maser P: Identification of GPI anchor attachment signals by a Kohonen self-organizing map. Bioinformatics. 2005, 21 (9): 1846-1852. 10.1093/bioinformatics/bti299.
Kall L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method. J Mol Biol. 2004, 338 (5): 1027-1036. 10.1016/j.jmb.2004.03.016.
Damveld RA, Arentshorst M, VanKuyk PA, Klis FM, van den Hondel CA, Ram AF: Characterisation of CwpA, a putative glycosylphosphatidylinositol-anchored cell wall mannoprotein in the filamentous fungus Aspergillus niger. Fungal Genet Biol. 2005, 42 (10): 873-885. 10.1016/j.fgb.2005.06.006.
Kamper J, Kahmann R, Bolker M, Ma LJ, Brefort T, Saville BJ, Banuett F, Kronstad JW, Gold SE, Muller O, et al: Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis. Nature. 2006, 444 (7115): 97-101. 10.1038/nature05248.
Martinez D, Larrondo LF, Putnam N, Gelpke MD, Huang K, Chapman J, Helfenbein KG, Ramaiya P, Detter JC, Larimer F, et al: Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78. Nat Biotechnol. 2004, 22 (6): 695-700. 10.1038/nbt967.
Ohm RA, de Jong JF, Lugones LG, Aerts A, Kothe E, Stajich JE, de Vries RP, Record E, Levasseur A, Baker SE, et al: Genome sequence of the model mushroom Schizophyllum commune. Nat Biotechnol. 2010, 28 (9): 957-963. 10.1038/nbt.1643.
Coleman JJ, Rounsley SD, Rodriguez-Carres M, Kuo A, Wasmann CC, Grimwood J, Schmutz J, Taga M, White GJ, Zhou S, et al: The genome of Nectria haematococca: contribution of supernumerary chromosomes to gene expansion. PLoS Genet. 2009, 5 (8): e1000618-10.1371/journal.pgen.1000618.
Cuomo CA, Guldener U, Xu JR, Trail F, Turgeon BG, Di Pietro A, Walton JD, Ma LJ, Baker SE, Rep M, et al: The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science. 2007, 317 (5843): 1400-1402. 10.1126/science.1143708.
Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H, et al: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434 (7036): 980-986. 10.1038/nature03449.
Nierman WC, Pain A, Anderson MJ, Wortman JR, Kim HS, Arroyo J, Berriman M, Abe K, Archer DB, Bermejo C, et al: Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature. 2005, 438 (7071): 1151-1156. 10.1038/nature04332.
Machida M, Asai K, Sano M, Tanaka T, Kumagai T, Terai G, Kusumoto K, Arima T, Akita O, Kashiwagi Y, et al: Genome sequencing and analysis of Aspergillus oryzae. Nature. 2005, 438 (7071): 1157-1161. 10.1038/nature04300.
Pel HJ, de Winde JH, Archer DB, Dyer PS, Hofmann G, Schaap PJ, Turner G, de Vries RP, Albang R, Albermann K, et al: Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88. Nat Biotechnol. 2007, 25 (2): 221-231. 10.1038/nbt1282.
van den Berg MA, Albang R, Albermann K, Badger JH, Daran JM, Driessen AJ, Garcia-Estrada C, Fedorova ND, Harris DM, Heijne WH, et al: Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum. Nat Biotechnol. 2008, 26 (10): 1161-1168. 10.1038/nbt.1498.
Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, Basturkmen M, Spevak CC, Clutterbuck J, et al: Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature. 2005, 438 (7071): 1105-1115. 10.1038/nature04341.
This research was supported by a NWO Van Gogh grant (VGP 89-185) to R.P de Vries. BH and PMC acknowledge financial support from the ANR to project E-Tricel (grant ANR-07-BIOE-006). The CAZy database is funded in part by GIS-IBiSA.
All authors read and approved the final version of the paper. EB participated in the growth profiling and CAZy annotation and was the main author of the paper. IB contributed to interpretation of the results. JvdB and AW contributed to the experimental work. PMC and BH participated and supervised the CAZy annotation and contributed to the manuscript. RPdV supervised the experimental work and contributed to the manuscript.