Skip to main content

Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes

Abstract

Background

Since its inception, the carbohydrate-active enzymes database (CAZy; http://www.cazy.org) has described the families of enzymes that cleave or build complex carbohydrates, namely the glycoside hydrolases (GH), the polysaccharide lyases (PL), the carbohydrate esterases (CE), the glycosyltransferases (GT) and their appended non-catalytic carbohydrate-binding modules (CBM). The recent discovery that members of families CBM33 and family GH61 are in fact lytic polysaccharide monooxygenases (LPMO), demands a reclassification of these families into a suitable category.

Results

Because lignin is invariably found together with polysaccharides in the plant cell wall and because lignin fragments are likely to act in concert with (LPMO), we have decided to join the families of lignin degradation enzymes to the LPMO families and launch a new CAZy class that we name “Auxiliary Activities” in order to accommodate a range of enzyme mechanisms and substrates related to lignocellulose conversion. Comparative analyses of these auxiliary activities in 41 fungal genomes reveal a pertinent division of several fungal groups and subgroups combining their phylogenetic origin and their nutritional mode (white vs. brown rot).

Conclusions

The new class introduced in the CAZy database extends the traditional CAZy families, and provides a better coverage of the full extent of the lignocellulose breakdown machinery.

Background

Decades of research into plant cell wall deconstruction have progressively revealed the enzyme activities involved in lignocellulolysis, for fundamental as well as for applied purposes. Lignin, the main non-carbohydrate structural component of the plant cell wall, is formed of an intricate network of phenolic compounds that constitutes a hard, hydrophobic and insoluble barrier. The plant cell wall carbohydrates form a complex matrix composed of crystalline and insoluble cellulose fibres intimately covered by an intricate network of more easily solubilizable carbohydrate polysaccharides, generically designated as hemicellulose, and, more distantly, by a gel-forming layer of galacturonic acid-containing polysaccharides, commonly designated as pectin. Upon maturation of the cell wall, the phenolic network of lignin is gradually deposited into this matrix and further stabilized through ester bonds to hemicellulose components. Even though the catabolism of each of these individual components can be independently studied, synergy in the breakdown of the various components is commonplace, because the accessibility of each constituent is limited by the presence and nature of the others. Therefore, a global view of the breakdown of all components of the plant cell wall, i.e. cellulose, hemicellulose, pectin and lignin is required to unravel the general mechanisms of lignocellulolysis and to explore the lignocellulolytic capabilities of organisms with sequenced genomes.

The assembly and breakdown of carbohydrate polymers and glycoconjugates are carried out by a diverse panel of Carbohydrate-Active enzymes, or ‘CAZymes’ [13]. The CAZymes classification system, based on amino acid sequence similarities, protein folds and enzymatic mechanisms, has been integrated and meticulously updated in the CAZy database (http://www.cazy.org) since 1998. CAZy currently incorporates more than 300 sequence families subdivided into the following classes: glycoside hydrolases (GH), glycosyltransferases (GT), polysaccharide lyases (PL), carbohydrate esterases (CE), and carbohydrate-binding modules (CBM). CAZy provides manually-curated information for each of these categories, making it now possible to cover all known families involved in cellulolysis, hemicellulolysis and pectinolysis.

Lignin forms a three-dimensional network with an aromatic backbone composed of dimethoxylated (syringyl, S), monomethoxylated (guaiacyl, G) and non-methoxylated (p-hydroxyphenyl, H) phenylpropanoid units. These phenylpropanoid units are linked by C-C and aryl-ether linkages. The structural complexity, insolubility and high molecular weight of lignin render its enzymatic degradation challenging but rewarding for microorganisms to gain access and convert the polysaccharides. Although lignin degradation is far from being fully deciphered, different microbial enzymatic strategies have already been evidenced to modify and/or degrade the recalcitrant lignin matrix [46]. The literature originally opposed two types of plant cell wall-degrading systems: the hydrolytic system responsible for carbohydrate degradation and the oxidative ligninolytic system responsible for lignin depolymerization. Contrary to the subtle specificity of carbohydrate-degrading enzymes, ligninolytic enzymes generate highly-reactive and non-specific free radicals that cleave carbon–carbon and ether inter-unit bonds. Thus, the enzymatic mechanism of ligninolysis has been described as an “enzymatic combustion” involving several non-specific oxidative species. However, such non-specific oxidative action is unlikely to be restricted to lignin digestion and all other plant cell wall components are also susceptible to undergoing oxidation. A number of carbohydrate oxidases have thus been identified in all kingdoms of life, although predominantly in fungi [7]. During microbial attack, these oxidases promote extended cleavage of plant cell wall components and yield reactive molecules (e.g. H2O2).

Moreover, recent evidence highlights the critical role of alternative enzymatic partners involved in the oxidation of other cell wall components. Originally described as cellulases, GH61 enzymes have recently been shown to be copper-dependent lytic polysaccharide monooxygenases (LPMO) that enhance cellulose degradation in concert with classical cellulases [811]. LPMO (GH61) catalyze the oxidative cleavage of cellulose using low molecular weight reducing agents such as ascorbate, gallate, reduced glutathione, and even lignin [12]. It has been also noted that the action of LPMO is strongly potentiated when combined with cellobiose dehydrogenase (CDH) [11, 13, 14]. Unlike traditional cellulases which require the presence of a loose, isolated substrate chain in the active site cleft or tunnel in order to perform cleavage, the strong oxidative mechanism of LPMO could allow oxidation at the surface of a cellulose crystalline microfibril, where chains are packed against each other. In fact, members of carbohydrate-binding module family CBM33 were first demonstrated to also be oxidative enzymes that cleave polysaccharide chains in crystalline chitin and cellulose [1517]. The NMR and isothermal titration calorimetry studies of CBM33 demonstrate that these enzymes are also copper-dependent LPMO [18, 19]. Therefore, LPMO are essential oxidative enzymes boosting the enzymatic conversion of recalcitrant polysaccharides. As a consequence, today’s modern consensus model for efficient enzymatic decomposition of cellulose should be extended to include not only the synergistic combination of processive and non-processive glycoside hydrolases but also its oxidative cleavage [20].

Cell wall carbohydrates and lignin co-occur in all land plants and are intimately interconnected and covalently cross-linked. We therefore believe that a strict boundary between (hemi)cellulolysis and ligninolysis is perhaps artificial, since functional inference often hinges on the choice of the substrates and analytical methods used in biochemical assays. Moreover, several publications have evidenced a concerted action of specific enzymes able to act on both carbohydrates and lignin [13, 21]. In the post-genomic era, where the quest for efficient enzymes for plant cell wall deconstruction has become a major research subject, the simultaneous identification of CAZymes and lignin-degrading enzymes is required to describe the full enzymatic repertoire necessary for plant cell wall deconstruction. Here, we describe the integration of a novel class in the CAZy database. This novel category, broadly termed “Auxiliary Activities” groups together the families of LPMO and the families of redox enzymes involved in lignin breakdown. Like the traditional CAZy families, the new families are based on sequence similarity with one or several biochemically-characterized founding member(s), ensuring that members of a given family share the same three-dimensional structure. Combined with the traditional CAZy families, this addition allows a complete description of the main actors involved in the degradation of all plant cell wall components. Like all enzyme classes in CAZy, the novel class features a hierarchical division in families and, where relevant, clans and subfamilies, all freely accessible at http://www.cazy.org.

Results and discussion

Classification of the auxiliary activities

We name “Auxiliary Activities” (AA) a widespread group of catalytic modules involved in plant cell wall degradation. The criteria for integration in the AA classification consists in the potential ability to help the original GH, PL and CE enzymes gain access to the carbohydrates encrusted in the plant cell wall. The AA category therefore encompasses a large class of modules not strictly restricted to a single catalytic reaction mechanism or to specific substrates. All AA families were created according to relevant biochemical characterizations retrieved from the literature for at least one member. The families and lists of AA proteins are continuously updated based on their similarity with the characterized models considered as functional anchors. Within families, subfamilies were manually defined according to their homology relationships between members of the focal family. The newly created AA class now encompasses, improves, replaces and extends a former classification dedicated to fungal ligninolytic enzymes [22]. Moreover, it aims to provide a global list of carbohydrate oxidases. This first release of the current AA class is divided into 10 families, and includes subfamilies for families AA1, AA3 and AA5 (Table 1). Subfamilies are designated using the family name plus a suffix indicating the subfamily. Each AA family and subfamily is concisely described below.

Table 1 Subfamily division and known activities in the AA families

Family AA1

The characterized AA1 enzymes are multicopper oxidases that use diphenols and related substances as donors with oxygen as the acceptor. The newly classified AA1 enzymes are widely distributed in a range of organisms from plants, fungi and insects to bacteria. The AA1 family is currently divided into 3 subfamilies as shown below:

Subfamily AA1_1: Laccases (EC 1.10.3.2) are blue copper oxidases that catalyze the one-electron oxidation of phenolics, aromatic amines and other electron-rich substrates via the reduction of oxygen to water [23]. These phenol oxidases require the presence of mediator molecules to perform the direct oxidation of high-redox-potential lignin components [24].

Subfamily AA1_2: Ferroxidases, also termed Fe(II):oxygen oxidoreductase, are enzymes present in this subfamily that catalyze the oxidization of iron II to iron III: 4 Fe2+ + 4 H+ + O2 = 4 Fe3+ + 2H2O. Ferroxidases are key enzymes involved in iron homeostasis and have been extensively studied in the yeast Saccharomyces cerevisiae (Fet3p and Fet5p) [25, 26]. Through Fe2+ oxidation, ferroxidases are potentially able to modulate Fenton reactions or reduce toxic levels of hydroxyl radicals, as suggested in Phanerochaete chrysosporium[27].

AA1_3: The laccase-like multicopper oxidases (LMCOs) present here are multicopper oxidases of unknown function. Fungal LMCOs are usually found in Ascomycota, and AA1_3 enzymes are suggested to be involved in diverse functions including lignin degradation, in addition to other physiological roles such as iron metabolism, pathogenic interactions, synthesis of conidial pigments, induction of fruiting bodies, and offensive/defensive fungal strategies [28, 29].

Family AA2

To date, this family contains class II lignin-modifying peroxidases that are part of the plant peroxidase superfamily [30, 31]. Initially described as true ligninases due to their high redox potential, class II peroxidases are secreted heme-containing enzymes that use hydrogen peroxide or organic peroxides as electron acceptors to catalyze a number of oxidative reactions in which two electrons are derived from substrate molecules to reduce the enzyme (via Compound I and Compound II intermediates) followed by a concomitant release of two water molecules [32]. Known AA2 enzymes include manganese peroxidase (MnP, EC 1.11.1.13), lignin peroxidase (LiP, EC 1.11.1.14) and versatile peroxidase (VP, EC 1.11.1.16). LiPs oxidize phenolic aromatic substrates and catalyze oxidative cleavages of C–C bonds and ether bonds in high-redox-potential non-phenolic aromatic substrates. MnPs oxidize Mn(II) to Mn(III) which in turn oxidize a variety of phenolic model compounds able to degrade and⁄ or modify lignin polymers. VPs are the third type of class II peroxidases included in the family AA2. VPs are the hybrid LIP-like MnP enzymes that combine the catalytic properties of LiP and MnP and are thus able to oxidize Mn(II), phenolic and non-phenolic substrates.

Family AA3

This family belongs to the glucose-methanol-choline (GMC) oxidoreductases family first outlined by Cavener [33]. AA3 enzymes are flavoproteins containing a flavin-adenine dinucleotide (FAD)-binding domain. Family AA3 can be divided into 4 subfamilies: AA3_1 (mostly covering cellobiose dehydrogenases), AA3_2 (including both aryl alcohol oxidase and glucose 1-oxidase), AA3_3 (alcohol oxidase) and AA3_4 (pyranose 2-oxidase).

Subfamily AA3_1: This subfamily includes cellobiose dehydrogenases (CDH, EC 1.1.99.18) that are extracellular hemoflavoenzymes produced under cellulolytic culture conditions by a number of lignocellulose-degrading fungi [34]. Although direct evidence is still missing, it has been proposed that CDH are involved in cellulose, hemicellulose and lignin biodegradation [35]. CDH employ a wide spectrum of electron acceptors to oxidize soluble cellodextrins to the corresponding lactones [21, 36]. They can generate hydroxyl radicals by Fenton-type reactions, thereby potentially oxidizing lignin. In the literature, the core CDH was originally described as a monomeric protein with a bipartite domain organization containing an N-terminal cytochrome domain of spectral class b connected to a C-terminal dehydrogenase domain containing the FAD redox cofactor. These two distinct modules in CDH evolved in parallel as fused genes [37]. Therefore, the modularity of CDH has been independently integrated in our classification as two separate modules, i.e. AA3_1 and AA8, respectively. AA3_1 represents the FAD-binding dehydrogenase module also designated cellobiose quinone oxidoreductase CBQ. The modularity of AA enzymes is discussed below.

Subfamily AA3_2: This family includes the two closely-related FAD-dependent enzymes, aryl-alcohol oxidase (AO, EC 1.1.3.7) and glucose 1-oxidase (GOX, EC 1.1.3.4). In our sequence-based classification, AO and GOX are strongly related and found in the same subfamily AA3_2. AO enzymes catalyze the oxidative dehydrogenation of several aromatic and aliphatic polyunsaturated alcohols with an α-carbon primary hydroxyl group, with the concomitant reduction of O2 to H2O2. GOX oxidizes the hydroxyl group at the C1 position of sugars. Both AO and GOX are H2O2-generating enzymes required for class II peroxidases to oxidize lignin. Lignin oxidation can alternatively be driven by reduction of H2O2 to hydroxyl free radical (OH) via ferrous iron recycling [38, 39]. GOX-generated H2O2 (intracellular enzymes) requires an export system that has not been reported to date [40]. In addition, related and more distant GMC oxidoreductases have been added in this subfamily based on sequence similarity to AO and GOX.

Subfamily AA3_3: Thus far, this subfamily contains members characterized as alcohol (methanol) oxidases (AlcOx; EC 1.1.3.13). AlcOx catalyze the oxidation of primary aliphatic alcohols to their respective aldehydes. They are able to use the methanol produced during lignin attack (demethoxylation) to supply the H2O2 necessary for hydroxyl radical generation in Fenton reactions [41].

Subfamily AA3_4: To date the only known activity in this subfamily is pyranose 2-oxidases (P2O, EC 1.1.3.10), which catalyze the C2-oxidation of D-glucopyranose and structurally analogous mono- and di-saccharides to the corresponding sugars while reducing O2 to H2O2[42].

Family AA4

The characterized members of this family are exclusively vanillyl-alcohol oxidases (VAO) that catalyze the conversion of a wide range of phenolic compounds bearing side chains at the para-position of the aromatic ring [43]. A broad spectrum of aromatic compounds is produced during the degradation of lignin, and AA4 enzymes are known to be active with a wide variety of these key intermediates.

Family AA5

The AA5 family belongs to the copper radical oxidases family and includes two subfamilies, namely AA5_1 containing characterized glyoxal oxidases and AA5_2 containing galactose oxidases. AA5_1 and AA5_2 share similar tertiary structures and virtually identical active sites despite different catalytic specificities and low sequence similarity.

Subfamily AA5_1: Glyoxal oxidases are copper-radical oxidases with a broad specificity able to oxidize simple aldehydes to the corresponding carboxylic acids [44]. Glyoxal oxidases are considered one of the central H2O2-generating enzymes. AA5_1 also integrates glyoxal-like sequences called cro1 to cro6 genes that share low overall sequence similarity to glyoxal oxidase but have conserved residues surrounding the catalytic site [45].

Subfamily AA5_2: Galactose oxidases (GAO) present here are copper-containing enzymes that catalyze a reaction comprising two separable half-reactions, oxidation of a primary alcohol, and reduction of O2 to H2O2[46]. GAO catalyze the oxidation of a wide range of carbohydrates (including galactose) but also primary alcohols into the corresponding aldehydes, with the reduction of O2 into H2O2.

Family AA6

This family contains members among which all experimentally characterized proteins are 1,4-benzoquinone reductases. These are intracellular enzymes involved in the biodegradation of aromatic compounds and in the protection of fungal cells from reactive quinone compounds. During the fungal degradation of aromatic compounds, quinones are key intermediate metabolites. Benzoquinone reductases are involved in a quinone redox cycle that generates extracellular Fenton reagents [4749].

Family AA7

The gluco-oligosaccharide oxidases (GOO) found in this family oxidize the C1 hydroxyl groups of β-1,4-linked sugars through two half-reactions featuring the oxidation of the reducing sugar to the lactone and spontaneous hydrolysis to the corresponding acid [7]. GOO are able to oxidize a variety of carbohydrates (D-glucose, maltose, lactose, cellobiose, malto- and cello-oligosaccharides) with the concomitant reduction of O2 to H2O2[50]. Known AA7 enzymes are potentially involved in the biotransformation or detoxification of lignocellulosic compounds.

Family AA8

Iron reductases consist in a cytochrome domain (protoheme IX) of spectral class b. Iron reductase was first described as the N-terminal hemic module found in the bipartite domain organization of the flavocytochrome CDH. In P. chrysosporium, AA8 stimulates the reduction of Fe(III) by acting as an electron sink or by directly reducing Fe(III) [51]. The physiological role of this family is unclear, but it could participate in the generation of highly reactive hydroxyl radicals (OH) via Fenton’s reaction. Consequently, AA8 proteins could be involved in the non-enzymatic cellulose chain breakage that involves hydroxyl radical attack [52]. Moreover, the C-terminal CBM1 might bring the generation of hydroxyl radicals by the Fenton’s reaction in the vicinity of cellulose.

Family AA9

This family includes the lytic polysaccharide monooxygenases previously classified as family GH61 glycoside hydrolases [8, 53]. AA9 enzymes have been recently characterized as copper-dependent polysaccharide monooxygenases containing highly conserved histidine residues implicated in a type-2 copper center. These enzymes are believed to act directly on cellulose, rendering it more accessible to more traditional CAZyme action [16]. The oxidative cleavages of cellulose catalysed by AA9 enzymes require reducing equivalents generated by the AA3_1 enzyme or alternatively low-molecular-weight reducing agents [8, 12]. AA9 enzymes are shown to oxidize the C1 carbon of the glucose ring structure but may also be able to oxidize C4 and C6 [9, 11, 14]. A number of AA9 members have been identified mainly in fungal genomes, indicating the widespread nature of this family, especially in the fungal wood decayers.

Family AA10

AA10 renames carbohydrate-binding module family CBM33. As for family AA9, the newly renamed AA10 members are copper-dependent monooxygenases that utilize molecular oxygen and external electron donors [18, 19]. Known AA10 enzymes are capable of oxidatively cleaving polysaccharide chains in crystalline chitin and cellulose to produce aldonic acids [15, 17]. AA10 action should increase the substrate accessibility for hydrolytic enzymes, and thus act synergistically with chitinases and cellulases. AA10 are structurally similar to AA9. Both families are devoid of classic substrate binding grooves or pockets yet they possess distinctive flat substrate-binding surfaces. Family AA10 also includes a conserved N-terminal histidine residue which is an integral component of the metal binding site. AA10 proteins are not limited to a specific kingdom, but are so far predominantly represented in bacteria [54].

Evolutionary relationships and modularity in AA families

Over the course of evolution, fungal genomes have been shaped by a series of complex processes leading to changes in biochemical function that are reflected in the amino acid sequence [5559]. Sequences evolve as a continuum in which positive, neutral and purifying selection drive the pace of the fixed mutations. Furthermore, new protein architectures (module gains, losses and shuffling) are frequently observed to occur and are thus recurrent evolutionary events shaping a given protein. As a rule, ligninolytic enzymes were found to be less modular than cellulolytic enzymes suggesting that close contact with polysaccharides is less compulsory. Notable exceptions were made for certain AA families in which various modular combinations were clearly evidenced (Figure 1). As recurrently observed in the established CAZy families, AA families have multimodular structures composed of catalytic modules associated to different ancillary modules like CBMs. This modular distribution is not uniformly shared across all AA families but remains strongly family-dependent. Families AA1, AA3, AA4, AA5, AA7, AA8, AA9 and AA10 are especially prone to modular organization.

Figure 1
figure 1

Examples of modular AA proteins. The additional modules represented include: carbohydrate-binding modules (CBMs; green), conserved modules of unknown function (X; red) and other domains (labelled grey boxes). Signal peptides (purple), proline/serine-rich linkers (cyan) and unassigned connecting regions (unlabelled grey boxes) are also shown.

The AA1 family was surprisingly found to be associated to various CBMs (CBM1, CBM18 and CBM20), which is singular for multicopper oxidases. In agreement with the diverse functions attributed to the AA1 family, the modular multicopper oxidases may be functionally linked to the modification of unexpected substrates such as chitin or even starch. Among the AA3 subfamilies, AA3_1 possesses a modular organization where the respective modules can be found either isolated or associated to the N-terminal cytochrome module AA8 and/or to a C-terminal CBM1. Despite different possible associations between AA8, AA3_1 and CBM1, their modular arrangement is strictly conserved. The glyoxal and glyoxal-like modules included in subfamily AA5_1 can contain one or more N-terminal WSC modules in tandem with the catalytic copper radical oxidase module. WSC domains are approx. 100 amino acid long with eight conserved cysteines, which are predicted to form disulfide bridges. The function of WSC domains is not known. Concerning subfamily AA5_2, N-terminal CBM13 and CBM32 modules can be found associated to the core catalytic module. A similar association was identified in AA7 in which associations to CBM12 and CBM18 modules are recurrent, which could point to a potential role of AA7 in the chitin modification. Interestingly, a shared modular organization was also evidenced between the AA4 and AA7 families that feature an N-terminal FAD-binding-4 module (pfam01565). The common modular organization between AA4 and AA7 reflects their shared ability to use FAD as a co-factor. As described above, AA8 is frequently associated to AA3_1, but AA8 can also be found fused to distinct modules. In fact, AA8 can be produced either as a single module or associated to AA3_1, CBM or other unknown modules. It is therefore possible that AA8 could be used as a key functional partner for the discovery of new candidate modules potentially involved in plant cell wall degradation.

Concerning the LPMO families AA9 and AA10 (formerly GH61 and CBM33, respectively), both are frequently multimodular and can bear a C-terminal CBM. To date, AA9 have only been found in association with CBM1 and CBM18, whereas AA10 are frequently fused to one or several modules of family CBM2, CBM5, CBM10, CBM12, and fibronectin type-III modules. Interestingly, the AA9 and AA10 enzymes always bear their ancilliary modules at the C-terminus, in agreement with the critical role of the N-terminal histidine for enzyme activity. Indeed, the first N-terminal histidine is involved in the coordination of the copper ion and the loss of this position would render the protein inactive.

In conclusion, diverse modular organizations are frequently detected in the AA families. The evolution of ligninolytic oxidoreductases makes use of classical evolutionary events in which shuffling events make it possible to associate different modules together. Gain and loss of modularity in AA families are both subsequently driven by common selective mechanisms (conservative, neutral or Darwinian selection), leading to different fates for “neo”-modular emergences. The modularity of ligninolytic modules has historically been underestimated. Additional modular architectures will doubtless be discovered in the near future by exploring novel and original taxonomic groups via systematic genome sequencing.

Occurrence of AA genes in fungal genomes

We selected 41 fungal genomes in order to compare their AA gene contents. The selected genomes mainly included several well-known white rot fungi (WR fungi), brown rot fungi (BR fungi) and plant pathogens from the Ascomycota (13 species) and Basidiomycota (28 species) divisions. A total of 1576 AA were identified in the 41 genomes. The AA distribution was not uniform across species (Figure 2). Coprinopsis cinerea, Stereum hirsutum and Trametes versicolor possess the highest number of AA (> 75 putative AA) while Malassezia globosa, Sporobolomyces roseus and Tremella mesenterica possess the lowest number of AA (≤ 5 putative AA). The well-known “industrial fungus” Trichoderma reesei is also particularly poor in AA enzymes as also demonstrated for the traditional CAZy families [60]. The distribution of the AA families was compared for each fungal genome, demonstrating that different enzymatic tactics are employed across the fungal kingdom (Figure 3, Additional files 1 and 2).

Figure 2
figure 2

Number of AA proteins in the 41 selected fungal genomes. Basidiomycetes are depicted in blue and ascomycetes are depicted in grey.

Figure 3
figure 3

Distribution of AA proteins in the 41 selected fungal genomes.

In order to organize the dataset into meaningful groups and identify coherent patterns, a hierarchical clustering of the number of ligninolytic AA family members was generated. Only ligninolytic families AA1 to AA8 were considered to focus exclusively on the lignin-degrading abilities of the 41 fungi including some important WR and BR representatives. Clustering of the AA dataset identified several major groups and subgroups of species that are discussed below according to phylogenetic origin and nutritional mode (Figure 4, Additional file 3).

Figure 4
figure 4

Comparison of the ligninolytic AA repertoires identified in the selected fungal genomes using double hierarchical clustering. Top tree: family number according to the AA classification. Left tree: Fungal genomes analyzed. A and B correspond to the Ascomycotina and Basidiomycotina divisions, respectively. The abundance of the different AA proteins within a family is represented by a colour scale from 0 (dark blue) to ≥ 10 occurrences (red) per species. Note that only ligninolytic AA families were selected (AA1 to AA8).

Phylogenetic origin

Two main groups appeared which corresponded to the Ascomycota and Basidiomycota phyla. In evolution, Ascomycota and Basidiomycota emerged around ~520 Myr after the split of the common ancestor (Dikarya) [61], with each group evolving its own set of ligninolytic enzymes. Ascomycota genomes contained no identifiable class II peroxidases (AA2) or ligninolytic laccase (AA1_1) whereas all Ascomycota possessed laccase-like multicopper oxidases (AA1_3). Surprisingly, AA1_3 (but not AA1_1) was also observed in the basidiomycete Auricularia delicata. This inconsistency between the current phylogeny including AA1_3 and the traditionally-accepted phylogeny of A. delicata poses a number of questions. A. delicata is one of the earliest branching Agaricomycete lineages, meaning that AA1_3 could either have been inherited and conserved from a common ancestor or acquired by horizontal transfer [62]. Ascomycota possess a higher number of AA3_1 (CDH) than Basidiomycota (averaging ~2.3 vs ~0.6, respectively). Only AA3_1 from Ascomycota could be fused to a C-terminal CBM1 module (~28% of total AA3_1) [34, 63]. The acquisition of a CBM1 may have occurred early before the diversification of the AA3_1 family in the Ascomycota division, or alternatively this modular organization may have been prematurely lost in Basidiomycota. As expected, Ascomycota are particularly prone to a high proportion of AA8 family members compared to Basidiomycota (~5 copies/genome vs ~1.5 copy/genome, respectively). Interestingly, AA8 can be individually fused to CBM1 at the C-terminal end (only), and this modular organization occurred in both fungal divisions. In terms of AA content, the non-ligninolytic Pichia stipitis stands out as containing only a few AA and clustering with low-AA-content genomes such as M. globosa. P. stipitis belongs to the subphylum Saccharomycotina, contrary to the other fungal Ascomycota genomes which belong to the subphylum Pezizomycotina. Basidiomycota presented significant differences in AA family distribution compared to Ascomycota. Subfamily AA1_1 ligninolytic laccases, like AA2 class II peroxidases, were only identified in Basidiomycota. AA5_1 was prevalent in Basidiomycota with a mean of ~5 copies per genome vs ~1 copy in Ascomycota.

Nutritional mode

Among the Basidiomycota sequenced to date, all BR fungi lack ligninolytic AA2 enzymes whereas WR fungi contain up to 26 AA2 proteins and average ~11 AA2 copies per genome. This high prevalence of AA2 enzymes in WR fungi re-enforced the central role of ligninolytic class II peroxidases in wood degradation, and argues for the use of AA2 as a family “marker” for predicting the wood-decaying activities and nutritional mode of different fungal groups (WR vs BR). Also worthy of mention is that known BR fungi such as F. pinicola, W. cocos and P. placenta possess “low redox potential peroxidases” lacking residues involved in Mn2+ binding. These “low redox potential peroxidases” may oxidize low redox-potential dyes and phenols, but their putative involvement in plant cell wall degradation needs to be experimentally studied before drawing more conclusions (Additional file 2) [64]. No AA2 enzyme was found in S. commune despite the fact that it was initially considered a WR fungus although its lignin-degrading capacity is limited [65]. Based on this similarity with BR fungi, S. commune clustered into the BR group. Concerning the AA3_2 family (AO and GOO), fewer enzymes were comparatively found among BR fungi, especially in F. pinicola, W. cocos, C. puteana, S. lacrymans and Dacryopinax sp. In WR fungi, AO enzymes are known to participate in the fungal degradation of lignin by providing the H2O2 required by ligninolytic peroxidases. This result is in agreement with the high AA2 content in WR fungi. Although subfamily AA3_1 was considered an important actor of lignocellulolysis [66], we found no significant correlation with the nutritional modes. However, the AA3_1 subfamily is prevalent in some “soft rot fungi” from the Ascomycota group, along with the AA1_3 subfamily.

Dacryopinax sp. (“jelly fungi” belonging to the Dacrymycetales order) has a limited number of AA families, including no member of subfamily AA1_1 among BR, and was clustered apart from other BR fungi. Based on the distribution of fungal species, the emergence of BR fungi seems to have followed a different evolutionary scenario. BR fungi may have appeared several times in evolution, accompanied by reductions in specific families. The chronology of divergences in fungal nutritional modes and the origin of lignin degradation were recently deciphered using comparative and functional genomics [52, 61]. BR fungi are able to perform an initial nonenzymatic attack on the wood cell wall, generating hydroxyl radicals (OH) extracellularly via the Fenton reaction. The generation of H2O2 could be achieved by oxidases, especially the high number of AA5_1 copper radical oxidases present in BR fungi (averaging ~3.4 copies per genome). In WR fungi, the species A. delicata, P. carnosa and P. chrysosporium belong to a separate cluster in which no subfamily AA1_1 was found. Although distant multicopper oxidases could be identified in these three WR species genomes, the absence of classical ligninolytic laccases revealed different enzymatic panels and, subsequently, different evolutionary groups within WR-mode lineages. Although the ability of laccases to remove lignin has been demonstrated during delignification [67], AA1_1 ligninolytic laccases are not compulsory sensu stricto in the WR mode of wood decay.

Co-occurrence of AA families in fungal genomes

In addition to identifying clusters of species, the AA families clustered into groups of biological significance (Figure 4). For instance, AA1_1, AA3_2, AA3_3, AA2 and AA5_1 clustered together. The AA2-family class II peroxidases require H2O2, and several extracellular enzymes have been postulated to be involved in extracellular H2O2 production. As expected, the main H2O2-producers included members of the families AA5_1 (known for glyoxal oxidase activity), AA3_2 (including AO and GOX) and AA3_3 (known as AlcOx). Interestingly, AA1_1 enzymes may also be involved synergistically with AA3_2 members to produce hydroxyl radicals, as suggested in this cluster. Indeed, laccases are able to produce hydroxyl radical (OH) during the oxidation of hydroquinones, and the addition of AO from subfamily AA3_2 strongly increases OH generation [68]. AA1_3, AA8 and AA1_3 also clustered together. AA8 and AA3_1 share a modular organization that means they are preferentially identified together. The laccase-like multicopper oxidases AA1_3 were grouped with the AA8 and AA3_1 families essentially based on the fact that all three families are found in high contents in Ascomycota. Interestingly, there is already evidence of synergistic action between these families, as it was demonstrated that the combined action of CDH and laccases resulted in substantial increases in the decolorization of synthetic dyes [6971].

Multigenicity in the AA families

The annotation of AA in the fungal genomes revealed a striking discrepancy between the numbers of members within families. The multiplicity of members mainly concerns families AA1_1, AA2, AA3_2 and AA9 in which up to 17, 26, 22 and 33 candidate proteins were assigned in individual genomes, respectively (Additional file 1). In contrast, other families seem particularly stable or less populated (i.e. AA3_1 ranged from 0 to 4 AA per genome). Interestingly, family expansions occurred in both ligninolytic AA (sub)families (AA1_1, AA2) and in AA9 which are likely involved in the oxidation of polysaccharides. From an evolutionary standpoint, the expansion of specific families is particularly intriguing. Family expansions are not randomly distributed, as only specific families are repeatedly targeted in several genomes, which is synonymous with convergence in evolution. A driving force for the emergence of multigenic families and gene novelty is the duplication event [72]. After duplication, the probability of fixation and maintenance of duplicates depends on numerous variables, and a combination of stochastic but also adaptive mechanisms take place. The remarkable expansion of the AA1_1, AA2 and AA9 families raises the question of the selective pressures constraining these families, and consequently the functional relevance at the organismal level (i.e. fitness). Different scenarios can be hypothesized to explain the multigenic AA families, ranging from i) functional redundancy to ii) functional diversification or iii) fine-tuned regulation of alternative genes. As fungi have to contend with the high complexity of the plant cell wall, the expansion of the AA families may reflect adaptations to the degradation of the substrates or to the use of different cofactors present in situ in the plant cell wall.

Conclusion

Plant cell wall degradation requires an arsenal of catalytic and non-catalytic proteins produced by lignocellulolytic microorganisms. Lignocellulolysis is a complex process, and fungi have evolved different enzymatic strategies to utilize these substrates as carbon sources. The growing number of published genome sequences and the allied “omics” data offers a unique opportunity for gaining deeper insight into the evolution of lignocellulolysis and the families involved in this biological machinery of high biotechnological potential. The novel AA classification provides a complementary insight into enzymatic lignocellulolysis by focusing on oxidative enzyme families. All known lignocellulolytic categories are now present in the continuously updated CAZy database. This expansion of the scope of the database is not restricted to the breakdown of plant cell wall as family AA10 (formerly CBM33) has been shown to act not only on cellulose, but also on another recalcitrant and crystalline polysaccharide, chitin [1517]. The presently released AA families are strongly biased towards the fungal kingdom. New kingdoms will be explored in the near future, especially in Bacteria where lignin breakdown may be more significant than previously thought [73]. Indeed, there are a number of literature reports of soil bacteria or aromatic-degrading bacteria from the digestive microbiota of animals that are able to depolymerise lignin (e.g. Nocardia, Rhodococcus, Sphingobium; Streptomyces) [7476]. The precise enzymology of bacterial lignin conversion is still in infancy and enzyme families involved in lignin breakdown will be added to the AA category in due time. It is hoped that the common annotation procedures and the dissemination provided by the CAZy database framework will allow cross fertilization of the various components of lignocellulose research. Rational integration of experimentally determined functional data as described earlier for traditional CAZymes [3] allied to biochemical expertise, will ultimately make it possible to provide a reliable and evolving functional annotation of all enzymes involved in the breakdown of lignocellulose and other recalcitrant biopolymers.

Methods

Family creation and propagation

The new families were created based on sequence similarity criteria analogous to those used in the past to define CAZy families [3]. Family creation starts with the identification of biochemically-characterized enzyme sequences that are decomposed into functional modules using available structural data whenever possible, followed by comparisons and manual curation. The modularity of these seeding examples is propagated to clearly identify homologous sequences using BLASTP [77] against GenBank sequences [78] in the initial stages of family definition. Families GH61 and CBM33 were just renamed (respectively AA9 and AA10) and proper information has been added to the headers of former families GH61 and CBM33 to indicate this fact. It is important to stress that inclusion of new sequences in the families is driven solely by sequence similarity and not by taxonomy or biochemistry considerations. In consequence, some families are obviously multi-functional and these were decomposed into subfamilies using criteria similar to those used in prior subfamily analyses [7981].

Hierarchical clustering

The number of each representative for each ligninolytic family (AA1 to AA8) was listed in order to compare the ligninolytic abilities of the 41 selected fungi. Data were ordered and grouped using the distance matrix computation and hierarchical clustering analysis calculated by MultiExperimentViewer [82]. A double clustering was carried out based on Pearson correlation as a distance matrix and complete linkage clustering as a linkage method.

Abbreviations

AA:

Auxiliary Activities

CAZy:

Carbohydrate-Active enZymes

AlcOx:

Alcohol oxidase

AO:

Aryl-alcohol oxidase

BR:

Brown rot

CBM:

Carbohydrate-binding modules

CDH:

Cellobiose dehydrogenase

CE:

Carbohydrate esterase

GAO:

Galactose oxidase

GH:

Glycoside hydrolase

GMC oxidoreductase:

glucose-methanol-choline oxidoreductase

GOO:

Gluco-oligosaccharide oxidase

GOX:

Glucose-1-oxidase

GT:

Glycosyltransferases

LMCO:

Laccase-like multicopper oxidase

LPMO:

Lytic polysaccharide monooxygenase

Myr:

Million years

PL:

Polysaccharide lyase

MnP:

Manganese peroxidase

LiP:

Lignin peroxidase

VAO:

Vanillyl-alcohol oxidase

VP:

Versatile peroxidase

WR:

White rot

References

  1. Henrissat B: A classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem J 1991, 280: 309-316.

    Article  CAS  Google Scholar 

  2. Henrissat B, Bairoch A: Updating the sequence-based classification of glycosyl hydrolases. Biochem J 1996, 316: 695-696.

    Article  Google Scholar 

  3. Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res 2009, 37: D233-D238. 10.1093/nar/gkn663

    Article  CAS  Google Scholar 

  4. Ruiz-Dueñas FJ, Martínez AT: Microbial degradation of lignin: how a bulky recalcitrant polymer is efficiently recycled in nature and how we can take advantage of this. Microb Biotechnol 2009, 2: 164-177. 10.1111/j.1751-7915.2008.00078.x

    Article  Google Scholar 

  5. Martínez AT, Ruiz-Dueñas FJ, Martínez MJ, Del Río JC, Gutiérrez A: Enzymatic delignification of plant cell wall: from nature to mill. Curr Opin Biotechnol 2009, 20: 348-357. 10.1016/j.copbio.2009.05.002

    Article  Google Scholar 

  6. Fernandez-Fueyo E, Ruiz-Dueñas FJ, Ferreira P, Floudas D, Hibbett DS, Canessa P, Larrondo LF, James TY, Seelenfreund D, Lobos S, Polanco R, Tello M, Honda Y, Watanabe T, Watanabe T, Ryu JS, Kubicek CP, Schmoll M, Gaskell J, Hammel KE, St John FJ, Vanden Wymelenberg A, Sabat G, Splinter BonDurant S, Syed K, Yadav JS, Doddapaneni H, Subramanian V, Lavín JL, Oguiza JA: Comparative genomics of Ceriporiopsis subvermispora and Phanerochaete chrysosporium provide insight into selective ligninolysis. Proc Natl Acad Sci U S A 2012, 109: 5458-5463. 10.1073/pnas.1119912109

    Article  CAS  Google Scholar 

  7. van Hellemond EW, Leferink NG, Heuts DP, Fraaije MW, van Berkel WJ: Occurrence and biocatalytic potential of carbohydrate oxidases. Adv Appl Microbiol 2006, 60: 17-54.

    Article  CAS  Google Scholar 

  8. Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JC, Johansen KS, Krogh KB, Jørgensen CI, Tovborg M, Anthonsen A, Tryfona T, Walter CP, Dupree P, Xu F, Davies GJ, Walton PH: Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci U S A 2011, 108: 15079-15084. 10.1073/pnas.1105776108

    Article  CAS  Google Scholar 

  9. Beeson WT, Phillips CM, Cate JH, Marletta MA: Oxidative cleavage of cellulose by fungal copper-dependent polysaccharide monooxygenases. J Am Chem Soc 2012, 134: 890-892. 10.1021/ja210657t

    Article  CAS  Google Scholar 

  10. Li X, Beeson WT, Phillips CM, Marletta MA, Cate JHD: Structural basis for substrate targeting and catalysis by fungal polysaccharide monooxygenases. Structure 2012, 20: 1051-1061. 10.1016/j.str.2012.04.002

    Article  Google Scholar 

  11. Bey M, Zhou S, Poidevin L, Henrissat B, Coutinho PM, Berrin JG, Sigoillot JC: Cello-oligosaccharide oxidation reveals differences between two lytic polysaccharide monooxygenases (family GH61) from Podospora anserina . Appl Environ Microbiol 2013, 79: 488-496. 10.1128/AEM.02942-12

    Article  CAS  Google Scholar 

  12. Westereng B, Ishida T, Vaaje-Kolstad G, Wu M, Eijsink VGH, Igarashi K, Samejima M, Ståhlberg J, Horn SJ, Sandgren M: The putative endoglucanase PcGH61D from Phanerochaete chrysosporium is a metal-dependent oxidative enzyme that cleaves cellulose. PLoS One 2011, 16: e27807.

    Article  Google Scholar 

  13. Langston JA, Brown K, Xu F, Borch K, Garner A, Sweeney MD: Cloning, expression, and characterization of a cellobiose dehydrogenase from Thielavia terrestris induced under cellulose growth conditions. Biochim Biophys Acta 1824, 2012: 802-812.

    Google Scholar 

  14. Phillips CM, Beeson WT, Cate JH, Marletta MA: Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa . ACS Chem Biol 2011, 6: 1399-1406. 10.1021/cb200351y

    Article  CAS  Google Scholar 

  15. Vaaje-Kolstad G, Horn SJ, van Aalten DMF, Synstad B, Eijsink VGH: The non-catalytic chitin-binding protein CBP21 from Serratia marcescens is essential for chitin degradation. J Biol Chem 2005, 280: 28492-28497. 10.1074/jbc.M504468200

    Article  CAS  Google Scholar 

  16. Vaaje-Kolstad G, Westereng B, Horn SJ, Liu ZL, Zhai H, Sørlie M, Eijsink VGH: An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science 2010, 330: 219-222. 10.1126/science.1192231

    Article  CAS  Google Scholar 

  17. Forsberg Z, Vaaje-Kolstad G, Westereng B, Bunæs AC, Stenstrøm Y, MacKenzie A, Sørlie M, Horn SJ, Eijsink VG: Cleavage of cellulose by a CBM33 protein. Protein Sci 2011, 20: 1479-1483. 10.1002/pro.689

    Article  CAS  Google Scholar 

  18. Vaaje-Kolstad G, Bøhle LA, Gåseidnes S, Dalhus B, Bjørås M, Mathiesen G, Eijsink VG: Characterization of the chitinolytic machinery of Enterococcus faecalis V583 and high-resolution structure of its oxidative CBM33 enzyme. J Mol Biol 2012, 416: 239-254. 10.1016/j.jmb.2011.12.033

    Article  CAS  Google Scholar 

  19. Aachmann FL, Sørlie M, Skjåk-Bræk G, Eijsink VG, Vaaje-Kolstad G: NMR structure of a lytic polysaccharide monooxygenase provides insight into copper binding, protein dynamics, and substrate interactions. Proc Natl Acad Sci U S A 2012, 109: 18779-18784. 10.1073/pnas.1208822109

    Article  CAS  Google Scholar 

  20. Mba Medie F, Davies GJ, Drancourt M, Henrissat B: Genome analyses highlight the different biological roles of cellulases. Nat Rev Microbiol 2012, 10: 227-234. 10.1038/nrmicro2729

    Article  Google Scholar 

  21. Henriksson G, Johansson G, Pettersson G: A critical review of cellobiose dehydrogenases. J Biotechnol 2000, 78: 93-113. 10.1016/S0168-1656(00)00206-6

    Article  CAS  Google Scholar 

  22. Levasseur A, Piumi F, Coutinho PM, Rancurel C, Asther M, Delattre M, Henrissat B, Pontarotti P, Asther M, Record E: FOLy: an integrated database for the classification and functional annotation of fungal oxidoreductases potentially involved in the degradation of lignin and related aromatic compounds. Fungal Genet Biol 2008, 45: 638-645. 10.1016/j.fgb.2008.01.004

    Article  CAS  Google Scholar 

  23. Baldrian P: Fungal laccases - occurrence and properties. FEMS Microbiol Rev 2006, 30: 215-242. 10.1111/j.1574-4976.2005.00010.x

    Article  CAS  Google Scholar 

  24. Cañas AI, Camarero S: Laccases and their natural mediators: biotechnological tools for sustainable eco-friendly processes. Biotechnol Adv 2010, 28: 694-705. 10.1016/j.biotechadv.2010.05.002

    Article  Google Scholar 

  25. Askwith C, Eide D, Van Ho A, Bernard PS, Li L, Davis-Kaplan S, Sipe DM, Kaplan J: The FET3 gene ofS. cerevisiaeencodes a multicopper oxidase required for ferrous iron uptake. Cell 1994, 76: 403-410. 10.1016/0092-8674(94)90346-8

    Article  CAS  Google Scholar 

  26. Kosman DJ: Molecular mechanisms of iron uptake in fungi. Mol Microbiol 2003, 47: 1185-1197. 10.1046/j.1365-2958.2003.03368.x

    Article  CAS  Google Scholar 

  27. Kersten P, Cullen D: Extracellular oxidative systems of the lignin-degrading Basidiomycete Phanerochaete chrysosporium . Fungal Genet Biol 2007, 44: 77-87. 10.1016/j.fgb.2006.07.007

    Article  CAS  Google Scholar 

  28. Levasseur A, Saloheimo M, Navarro D, Andberg M, Pontarotti P, Kruus K, Record E: Exploring laccase-like multicopper oxidase genes from the ascomycete Trichoderma reesei : a functional, phylogenetic and evolutionary study. BMC Biochem 2010, 11: 32. 10.1186/1471-2091-11-32

    Article  Google Scholar 

  29. Hölker U, Dohse J, Höfer M: Extracellular laccases in ascomycetes Trichoderma atroviride and Trichoderma harzianum . Folia Microbiol 2002, 47: 423-427. 10.1007/BF02818702

    Article  Google Scholar 

  30. Koua D, Cerutti L, Falquet L, Sigrist CJ, Theiler G, Hulo N, Dunand C: PeroxiBase: a database with new tools for peroxidase family classification. Nucleic Acids Res 2009, 37: D261-D266. 10.1093/nar/gkn680

    Article  CAS  Google Scholar 

  31. Lundell TK, Mäkelä MR, Hildén K: Lignin-modifying enzymes in filamentous basidiomycetes–ecological, functional and phylogenetic review. J Basic Microbiol 2010, 50: 5-20. 10.1002/jobm.200900338

    Article  CAS  Google Scholar 

  32. Hofrichter M, Ullrich R, Pecyna MJ, Liers C, Lundell T: New and classic families of secreted fungal heme peroxidases. Appl Microbiol Biotechnol 2010, 87: 871-897. 10.1007/s00253-010-2633-0

    Article  CAS  Google Scholar 

  33. Cavener DR: GMC oxidoreductases—a newly defined family of homologous proteins with diverse catalytic activities. J Mol Biol 1992, 223: 811-814. 10.1016/0022-2836(92)90992-S

    Article  CAS  Google Scholar 

  34. Zámocký M, Ludwig R, Peterbauer C, Hallberg BM, Divne C, Nicholls P, Haltrich D: Cellobiose dehydrogenase- a flavocytochrome from wood-degrading, phytopathogenic and saprotropic fungi. Curr Protein Pept Sci 2006, 7: 255-280. 10.2174/138920306777452367

    Article  Google Scholar 

  35. Kremer SM, Wood PM: Cellobiose oxidase from Phanerochaete chrysosporium as a source of Fenton’s reagent. Biochem Soc Trans 1992, 20: 110S.

    Article  CAS  Google Scholar 

  36. Turbe-Doan A, Arfi Y, Record E, Estrada-Alvarado I, Levasseur A: Heterologous production of cellobiose dehydrogenases from the basidiomycete Coprinopsis cinerea and the ascomycete Podospora anserina and their effect on saccharification of wheat straw . Appl Microbiol Biotechnol 2012. In press

    Google Scholar 

  37. Zámocký M, Hallberg M, Ludwig R, Divne C, Haltrich D: Ancestral gene fusion in cellobiose dehydrogenases reflects a specific evolution of GMC oxidoreductases in fungi. Gene 2004, 338: 1-14. 10.1016/j.gene.2004.04.025

    Article  Google Scholar 

  38. Hernández-Ortega A, Ferreira P, Martínez AT: Fungal aryl-alcohol oxidase: a peroxide-producing flavoenzyme involved in lignin degradation. Appl Microbiol Biotechnol 2012, 93: 1395-1410. 10.1007/s00253-011-3836-8

    Article  Google Scholar 

  39. Martínez AT, Speranza M, Ruiz-Dueñas FJ, Ferreira P, Camarero S, Guillén F, Martínez MJ, Gutiérrez A, del Río JC: Biodegradation of lignocellulosics: microbial, chemical, and enzymatic aspects of the fungal attack of lignin. Int Microbiol 2005, 8: 195-204.

    Google Scholar 

  40. Eriksson KE, Pettersson B, Volc J, Musílek V: Formation and partial characterization of glucose-2-oxidase, a H 2 O 2 producing enzyme in Phanerochaete chrysosporium . Appl Microbiol Biotechnol 1986, 23: 257-262.

    Article  CAS  Google Scholar 

  41. Daniel G, Volc J, Filonova L, Plihal O, Kubátová E, Halada P: Characteristics of Gloeophyllum trabeum alcohol oxidase, an extracellular source of H 2 O 2 in brown rot decay of wood. Appl Environ Microbiol 2007, 73: 6241-6253. 10.1128/AEM.00977-07

    Article  CAS  Google Scholar 

  42. Giffhorn F: Fungal pyranose oxidases: occurrence, properties and biotechnical applications in carbohydrate chemistry. Appl Microbiol Biotechnol 2000, 54: 727-740. 10.1007/s002530000446

    Article  CAS  Google Scholar 

  43. de Jong E, van Berkel WJ, van der Zwan RP, de Bont JA: Purification and characterization of vanillyl-alcohol oxidase from Penicillium simplicissimum . A novel aromatic alcohol oxidase containing covalently bound FAD. Eur J Biochem 1992, 208: 651-657. 10.1111/j.1432-1033.1992.tb17231.x

    Article  CAS  Google Scholar 

  44. Whittaker MM, Kersten PJ, Nakamura N, Sanders-Loehr J, Schweizer ES, Whittaker JW: Glyoxal oxidase from Phanerochaete chrysosporium is a new radical-copper oxidase. J Biol Chem 1996, 271: 681-687. 10.1074/jbc.271.2.681

    Article  CAS  Google Scholar 

  45. Vanden Wymelenberg A, Sabat G, Mozuch M, Kersten PJ, Cullen D, Blanchette RA: Structure, organization, and transcriptional regulation of a family of copper radical oxidase genes in the lignin-degrading basidiomycete Phanerochaete chrysosporium . Appl Environ Microbiol 2006, 72: 4871-4877. 10.1128/AEM.00375-06

    Article  CAS  Google Scholar 

  46. Whittaker JW: Free radical catalysis by galactose oxidase. Chem Rev 2003, 103: 2347-2363. 10.1021/cr020425z

    Article  CAS  Google Scholar 

  47. Brock BJ, Rieble S, Gold MH: Purification and characterization of a 1,4-Benzoquinone reductase from the basidiomycete Phanerochaete chrysosporium . Appl Environ Microbiol 1995, 61: 3076-3081.

    CAS  Google Scholar 

  48. Lee SS, Moon DS, Choi HT, Song HG: Purification and characterization of an intracellular NADH:quinone reductase from Trametes versicolor . J Microbiol 2007, 45: 333-338.

    CAS  Google Scholar 

  49. Jensen K Jr, Ryan Z, Wymelenberg A, Cullen D, Hammel K: An NADH:quinone oxidoreductase active during biodegradation by the brown-rot basidiomycete Gloeophyllum trabeum . Appl Environ Microbiol 2002, 68: 2699-2703. 10.1128/AEM.68.6.2699-2703.2002

    Article  CAS  Google Scholar 

  50. Fan Z, Oguntimein GB, Reilly PJ: Characterization of kinetics and thermostability of Acremonium strictum glucooligosaccharide oxidase. Biotechnol Bioeng 2000, 68: 231-237. 10.1002/(SICI)1097-0290(20000420)68:2<231::AID-BIT12>3.0.CO;2-D

    Article  CAS  Google Scholar 

  51. Henriksson G, Johansson G, Pettersson G: Is cellobiose oxidase from Phanerochaete chrysosporium a one-electron reductase? Biochim Biophys Acta 1993, 1144: 184-190. 10.1016/0005-2728(93)90171-B

    Article  CAS  Google Scholar 

  52. Eastwood DC, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A, Asiegbu FO, Baker SE, Barry K, Bendiksby M, Blumentritt M, Coutinho PM, Cullen D, de Vries RP, Gathman A, Goodell B, Henrissat B, Ihrmark K, Kauserud H, Kohler A, LaButti K, Lapidus A, Lavin JL, Lee YH, Lindquist E, Lilly W, Lucas S, Morin E, Murat C, Oguiza JA: The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. Science 2011, 333: 762-765. 10.1126/science.1205411

    Article  CAS  Google Scholar 

  53. Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, Xu F, Cherry J, Larsen S, Lo LL: Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: Structure and function of a large, enigmatic family. Biochemistry 2010, 49: 3305-3316. 10.1021/bi100009p

    Article  CAS  Google Scholar 

  54. Horn SJ, Vaaje-Kolstad G, Westereng B, Eijsink VG: Novel enzymes for the degradation of cellulose. Biotechnol Biofuels 2012, 5: 45. 10.1186/1754-6834-5-45

    Article  CAS  Google Scholar 

  55. Kimura M: The neutral theory of molecular evolution. Cambridge UK: Cambridge University Press; 1983.

    Book  Google Scholar 

  56. Ganfornina MD, Sanchez D: Generation of evolutionary novelty by functional shift. Bioessays 1999, 21: 432-439. 10.1002/(SICI)1521-1878(199905)21:5<432::AID-BIES10>3.0.CO;2-T

    Article  CAS  Google Scholar 

  57. Lynch M, Conery JS: The origins of genome complexity. Science 2003, 302: 1401-1404. 10.1126/science.1089370

    Article  CAS  Google Scholar 

  58. Levasseur A, Pontarotti P: The role of duplications in the evolution of genomes highlights the need for evolutionary-based approaches in comparative genomics. Biol Direct 2011, 6: 11. 10.1186/1745-6150-6-11

    Article  CAS  Google Scholar 

  59. Levasseur A, Orlando L, Bailly X, Milinkovitch MC, Danchin EG, Pontarotti P: Conceptual bases for quantifying the role of the environment on gene evolution: the participation of positive selection and neutral evolution. Biol Rev Camb Philos Soc 2007, 82: 551-572. 10.1111/j.1469-185X.2007.00024.x

    Article  Google Scholar 

  60. Martinez D, Berka RM, Henrissat B, Saloheimo M, Arvas M, Baker SE, Chapman J, Chertkov O, Coutinho PM, Cullen D, Danchin EG, Grigoriev IV, Harris P, Jackson M, Kubicek CP, Han CS, Ho I, Larrondo LF, de Leon AL, Magnuson JK, Merino S, Misra M, Nelson B, Putnam N, Robbertse B, Salamov AA, Schmoll M, Terry A, Thayer N, Westerholm-Parvinen A: Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina ). Nat Biotechnol 2008, 26: 553-560. 10.1038/nbt1403

    Article  CAS  Google Scholar 

  61. Floudas D, Binder M, Riley R, Barry K, Blanchette RA, Henrissat B, Martínez AT, Otillar R, Spatafora JW, Yadav JS, Aerts A, Benoit I, Boyd A, Carlson A, Copeland A, Coutinho PM, de Vries RP, Ferreira P, Findley K, Foster B, Gaskell J, Glotzer D, Górecki P, Heitman J, Hesse C, Hori C, Igarashi K, Jurgens JA, Kallen N, Kersten P: The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. Science 2012, 336: 1715-1719. 10.1126/science.1221748

    Article  CAS  Google Scholar 

  62. Rosewich UL, Kistler HC: Role of horizontal gene transfer in the evolution of fungi. Annu Rev Phytopathol 2000, 38: 325-363. 10.1146/annurev.phyto.38.1.325

    Article  CAS  Google Scholar 

  63. Harreither W, Sygmund C, Augustin M, Narciso M, Rabinovich ML, Gorton L, Haltrich D, Ludwig R: Catalytic properties and classification of cellobiose dehydrogenases from ascomycetes. Appl Environ Microbiol 2011, 77: 1804-1815. 10.1128/AEM.02052-10

    Article  CAS  Google Scholar 

  64. Martinez AT: Molecular biology and structure-function of lignin-degrading heme peroxidases. Enzym Microb Technol 2002, 30: 425-444. 10.1016/S0141-0229(01)00521-X

    Article  CAS  Google Scholar 

  65. Ohm RA, de Jong JF, Lugones LG, Aerts A, Kothe E, Stajich JE, de Vries RP, Record E, Levasseur A, Baker SE, Bartholomew KA, Coutinho PM, Erdmann S, Fowler TJ, Gathman AC, Lombard V, Henrissat B, Knabe N, Kües U, Lilly WW, Lindquist E, Lucas S, Magnuson JK, Piumi F, Raudaskoski M, Salamov A, Schmutz J, Schwarze FW, van Kuyk PA, Horton JS: Genome sequence of the model mushroomSchizophyllum commune. Nat Biotechnol 2010, 28: 957-963. 10.1038/nbt.1643

    Article  CAS  Google Scholar 

  66. Cameron MD, Aust SD: Cellobiose dehydrogenase-an extracellular fungal flavocytochrome. Enzyme Microb Technol 2001, 28: 129-138. 10.1016/S0141-0229(00)00307-0

    Article  CAS  Google Scholar 

  67. Gutiérrez A, Rencoret J, Cadena EM, Rico A, Barth D, del Río JC, Martínez AT: Demonstration of laccase-based removal of lignin from wood and non-wood plant feedstocks. Bioresour Technol 2012, 119: 114-122.

    Article  Google Scholar 

  68. Guillén F, Gómez-Toribio V, Martínez MJ, Martínez AT: Production of hydroxyl radical by the synergistic action of fungal laccase and aryl alcohol oxidase. Arch Biochem Biophys 2000, 383: 142-147. 10.1006/abbi.2000.2053

    Article  Google Scholar 

  69. Ciullini I, Tilli S, Scozzafava A, Briganti F: Fungal laccase, cellobiose dehydrogenase, and chemical mediators: combined actions for the decolorization of different classes of textile dyes. Bioresour Technol 2008, 99: 7003-7010. 10.1016/j.biortech.2008.01.019

    Article  CAS  Google Scholar 

  70. Enayatzamir K, Tabandeh F, Yakhchali B, Alikhani HA, Rodríguez Couto S: Assessment of the joint effect of laccase and cellobiose dehydrogenase on the decolouration of different synthetic dyes. J Hazard Mater 2009, 169: 176-181. 10.1016/j.jhazmat.2009.03.088

    Article  CAS  Google Scholar 

  71. Tilli S, Ciullini I, Scozzafava A, Briganti F: Differential decolorization of textile dyes in mixtures and the joint effect of laccase and cellobiose dehydrogenase activities present in extracellular extracts from Funalia trogii . Enzyme Microb Technol 2011, 49: 465-471. 10.1016/j.enzmictec.2011.08.002

    Article  CAS  Google Scholar 

  72. Ohno S: Evolution by gene duplication. New York: Springer; 1970.

    Book  Google Scholar 

  73. Zimmermann W: Degradation of lignin by bacteria. J Biotechnol 1990, 13: 119-130. 10.1016/0168-1656(90)90098-V

    Article  CAS  Google Scholar 

  74. Bugg TD, Ahmad M, Hardiman EM, Singh R: The emerging role for bacteria in lignin degradation and bio-product formation. Curr Opin Biotechnol 2011, 22: 394-400. 10.1016/j.copbio.2010.10.009

    Article  CAS  Google Scholar 

  75. Masai E, Katayama Y, Fukuda M: Genetic and biochemical investigations on bacterial catabolic pathways for lignin-derived aromatic compounds. Biosci Biotechnol Biochem 2007, 71: 1-15. 10.1271/bbb.60437

    Article  CAS  Google Scholar 

  76. Bugg TD, Ahmad M, Hardiman EM, Rahmanpour R: Pathways for degradation of lignin in bacteria and fungi. Nat Prod Rep 2011, 28: 1883-1896. 10.1039/c1np00042j

    Article  CAS  Google Scholar 

  77. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215: 403-410.

    Article  CAS  Google Scholar 

  78. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res 2005, 33: D34-D38. 10.1093/nar/gni032

    Article  CAS  Google Scholar 

  79. Stam MR, Danchin EG, Rancurel C, Coutinho PM, Henrissat B: Dividing the large glycoside hydrolase family 13 into subfamilies: towards improved functional annotations of alpha-amylase-related proteins. Protein Eng Des Sel 2006, 19: 555-562. 10.1093/protein/gzl044

    Article  CAS  Google Scholar 

  80. Lombard V, Bernard T, Rancurel C, Brumer H, Coutinho PM, Henrissat B: A hierarchical classification of polysaccharide lyases for glycogenomics. Biochem J 2010, 432: 437-444. 10.1042/BJ20101185

    Article  CAS  Google Scholar 

  81. Aspeborg H, Coutinho PM, Wang Y, Brumer H 3rd, Henrissat B: Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5). BMC Evol Biol 2012, 12: 186. 10.1186/1471-2148-12-186

    Article  CAS  Google Scholar 

  82. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques 2003, 34: 374-378.

    CAS  Google Scholar 

Download references

Acknowledgements

We wish to thank Corinne Rancurel for her technical help. The support of the French National Institute for Agricultural Research, INRA (Departments MICA and CEPIA) is gratefully acknowledged. This work was performed in NEMO, a project under EU’s 7th framework program, grant number 222699.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Anthony Levasseur or Bernard Henrissat.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AL and BH conceived the study. AL designed and carried out the bioinformatics studies. ED and VL contributed to the data integration into the database. PMC helped to conceive the study. AL drafted the manuscript. AL and BH wrote the final manuscript. All authors have read and approved the final manuscript for publication.

Electronic supplementary material

Additional file 1: Number of AA proteins in the 41 selected fungal genomes. (PDF 17 KB)

Additional file 2: List of AA proteins identified in the 41 selected fungal genomes. (XLSX 104 KB)

13068_2013_263_MOESM3_ESM.pdf

Additional file 3: Schematic distribution of ligninolytic AA in fungal genomes. Each genome is represented as a circle. The total number of AA is given in the middle of each circle. Groups 1 to 4 were defined according to the hierarchical clustering. (PDF 7277 kb) (PDF 7 MB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Levasseur, A., Drula, E., Lombard, V. et al. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels 6, 41 (2013). https://doi.org/10.1186/1754-6834-6-41

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1754-6834-6-41

Keywords