Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5)
BMC Evolutionary Biology volume 12, Article number: 186 (2012)
The large Glycoside Hydrolase family 5 (GH5) groups together a wide range of enzymes acting on β-linked oligo- and polysaccharides, and glycoconjugates from a large spectrum of organisms. The long and complex evolution of this family of enzymes and its broad sequence diversity limits functional prediction. With the objective of improving the differentiation of enzyme specificities in a knowledge-based context, and to obtain new evolutionary insights, we present here a new, robust subfamily classification of family GH5.
About 80% of the current sequences were assigned into 51 subfamilies in a global analysis of all publicly available GH5 sequences and associated biochemical data. Examination of subfamilies with catalytically-active members revealed that one third are monospecific (containing a single enzyme activity), although new functions may be discovered with biochemical characterization in the future. Furthermore, twenty subfamilies presently have no characterization whatsoever and many others have only limited structural and biochemical data. Mapping of functional knowledge onto the GH5 phylogenetic tree revealed that the sequence space of this historical and industrially important family is far from well dispersed, highlighting targets in need of further study. The analysis also uncovered a number of GH5 proteins which have lost their catalytic machinery, indicating evolution towards novel functions.
Overall, the subfamily division of GH5 provides an actively curated resource for large-scale protein sequence annotation for glycogenomics; the subfamily assignments are openly accessible via the Carbohydrate-Active Enzyme database at http://www.cazy.org/GH5.html.
Carbohydrates, in the form of mono-, di-, oligo-, and polysaccharides, as well as glycoconjugates, play fundamental roles in all forms of life . Beyond their role in energy storage, carbohydrates are central to diverse biological processes such as host-pathogen interactions, signal transduction, inflammation, intracellular trafficking, diseases, and differentiation/development. Not least, as structural components of terrestrial biomass, carbohydrates comprise approximately 75% of the carbon fixed annually by primary production . Sugar-rich plant cell walls, seeds, and tubers thus represent a renewable material with significant potential to address energy and material needs.
A striking feature of carbohydrates is their remarkable structural complexity, due to a rich diversity of monosaccharide building blocks, and the possibility of numerous stereo- and regiospecific linkages , which give rise to both simple linear and complex, highly branched molecules . A decade of investments in genomics and proteomics has greatly improved our interpretation of the molecular language of the cell, but deciphering the complex carbohydrate-based information in the biomolecular landscape is still in its infancy. Indeed, glycomics has been identified both as “the last frontier of molecular and cellular biology”  as well as an “emerging technology that will change the world” .
Functional analysis of glycans and glycoconjugates is complicated by the fact that they are not direct genetic products, but are instead synthesized, recognized, modified, and degraded by a plethora of carbohydrate-active enzymes (CAZymes) and binding proteins. In the synthetic direction, phosphosugar-dependent glycosyltransferases (GTs) catalyze the formation of glycosidic linkages, whereas their breakdown is mediated by glycoside hydrolases (GHs) and polysaccharide lyases (PLs), with the assistance of carbohydrate esterases (CEs). The structural diversity of carbohydrates is reflected in an abundance of CAZyme-encoding genes, which comprise 1-3% of the genome of most organisms . Expanding and harnessing knowledge of the complexity of the “CAZome” is thus essential to understanding the complexity of the glycome.
The protein sequence-based classification of CAZymes was initiated in 1991 as a complement to the long-standing Enzyme Commission (EC) number system , which is based solely on enzyme activities . Given the prevalence of convergent evolution of enzymes that cleave glycosidic bonds, as well as the demonstrable catalytic promiscuity of individual enzymes, sequence-based classification has proven to be a robust way to unify information on enzyme structure, specificity, and mechanism, which provides enormous predictive power . Initially motivated by a need to delineate cellulases (EC 126.96.36.199) into distinct structural families , the first incarnation of the GH family classification, as such, comprised 35 GH families . The number of families increased steadily with the growing interest in Glycobiology so that, as of August 2012, 130 sequence-based families of GHs have been defined in the continuously updated CAZy database .
Presently, one of the largest GH families is GH5, historically known as “cellulase family A” as it was the first cellulase family described . GH5 exemplifies a family with a large variety of specificities: it currently contains close to 20 experimentally determined enzyme activities denoted with an EC number. The abundance of GH5 enzymes in different ecological niches has been highlighted by their frequent identification in metagenomes of diverse microbial communities [12–14], as well as the genomes of individual organisms . As with other CAZyme families , GH5 members are commonly found to be encoded as parts of multi-modular polypeptide chains containing other catalytic, substrate-binding, and functionally unidentified or yet to be described modules.
Within the large GH5 family, a discernible diversity of sequences was observed soon after its creation. The first five subfamilies of GH5 (A1-A5) were identified as early as 1990 . Subfamily A6 was introduced in 1997  and the following year eukaryotic and prokaryotic β-mannanases were assigned to A7 and A8, respectively . Subsequently, subfamily A9 was introduced in a study, which notably also suggested the merger of A5 and A6 . Finally, A10 was the most recently defined GH5 subfamily , while new subfamilies that presently lack a unique identifier have also been suggested [21, 22]. Family GH5 belongs to clan GH-A, which presently groups 19 GH families to form the largest set of evolutionarily related GH families described in CAZy thus far (a clan is a group of families that arise from a common but very distant ancestor; despite weak sequence similarity, clan members share conserved protein fold and catalytic machinery).
Families such as GH5 were originally defined with a very small number of sequences. With the accumulation of an increasing body of sequence data, the relationship between the original families has sometimes changed enough to merit reexamination of family membership. Very recently, detailed three-dimensional structural analysis led to the reclassification of several GH5 sequences into family GH30 based on the organization of secondary structural elements around the conserved (β/α)8 fold of the catalytic module .
Given the continuing expansion in sequence numbers and the partial GH5/GH30 reclassification, it is clear that a global re-analysis of the subfamily division of GH5 is now needed. The rapid accumulation of genomic data in the past decade revealed a complex and varied sequence space, with the consequence that a substantial portion of GH5 family members are currently not assigned to any subfamily. This situation will only become worse as the rate of (meta)genomic sequencing continues to increase with phenomenal rapidity. Further, this flood of data will cause an increasing reliance on computer-based annotation, which necessarily requires a robust framework to produce meaningful functional predictions. The division of CAZyme families into subfamilies based on phylogenetic analysis has been applied as a successful approach to meet this challenge: Subfamily classification of GH13, GH30 and all of the PL families has demonstrated that the majority of the defined subfamilies were monospecific, thus indicating a significantly better correlation of substrate specificity between sequences at the subfamily level than the family level [23–25]. Significantly, the division into subfamilies allows the identification of currently uncharacterized subfamilies that can subsequently be analyzed biochemically and structurally to potentially unveil new activities.
Hence, we present here an improved, robust subfamily classification for GH5 by employing a large-scale analysis of all publicly available sequences. Our intention is that the introduction of this additional hierarchical level across this important GH family will serve to guide enzyme discovery, structure-function analysis, and biocatalyst improvement in post-genomic efforts. Not least, many enzyme activities relevant to biomass analysis and conversion are found in GH5 (e.g., cellulases, mannanases, xylanases, galactanases, and xyloglucanases), as are enzymes with biomedical applications . Significantly, the present analysis unveiled a large number of sparsely or incompletely characterized subfamilies that may still hide a number of unsuspected activities and singular structural features.
Results and discussion
Our bioinformatics approach allowed the division of close to 2300 GH5 catalytic modules into 51 distinct subfamilies, as shown in the global phylogenetic tree (Figure 1 and Additional file 1: Figure S1); subfamily information is summarized in Table 1. Subfamily naming follows the procedure devised for GH13, where the family number is followed by an Arabic numeral that reflects the order of creation : GH5_1 to GH5_53. This series is essentially continuous, with a few exceptions due to historical reasons: All of the previously described subfamilies (A1-A10) have been re-identified in the current investigation except for A3 and A4, which are merged into a single subfamily GH5_4 and A5 and A6 which are unified in subfamily GH5_5 (Figure 1). To maintain consistency with earlier literature, the re-identified historical subfamilies have retained the original Arabic numeral. For example, the subfamily formerly known as A2 is hereby designated GH5_2. The absence of subfamilies GH5_3 and GH5_6 reflects the two fusion events involving the historical subfamilies described above .
In addition to the new and historical designations, the taxonomical range of the included sequences, experimentally determined enzyme activities and representative 3-D structures are presented in Table 1 for each subfamily. Notably, all of the 33 enzymes with a solved 3-D structure have been assigned to a subfamily, resulting in thirteen individual subfamilies out of 51 with at least one structural representative. Genes that encode GH5 enzymes are present in most organisms ranging from Archaea and Eubacteria to Eukaryotes, e.g. fungi and plants. From an anthropocentric perspective, a GH5 member is notably lacking in the human genome. Examples of metazoan GH5 genes are also notably scarce and are limited to nematodes, mollusks, and arthropods, likely resulting from horizontal transfer. For example, several independent horizontal events of transfer of cellulase and xylanase genes from Bacteria to nematodes have been described  and the transfer of a bacterial β-mannanase to an insect was recently documented . The taxonomical range at the subfamily level is, naturally, more restricted. A few smaller subfamilies are currently specific to certain types of organisms (Table 1). For example, eight subfamilies (GH5_15, GH5_16, GH5_23, GH5_24, GH5_30, GH5_31, GH5_49 and GH5_50) contain only fungal sequences. Subfamily GH5_14 contains only plant members, whereas members of GH5_20 come exclusively from Stramenopiles. These limited taxonomic distributions may represent biological (catalytic) specialization, if not biased by a still incomplete genome sequencing of organisms.
The definition of subfamilies was restricted to phylogenetic clades with five or more members from different organisms available in the public protein databases (GenBank and UniProt) in order to capture sufficient diversity for a robust subfamily definition. Using these criteria, the overall success rate of subfamily grouping was of approximately 80%, i.e., about 20% of the analyzed GH5 sequences could not be classified into subfamilies having at least five public members. The sequences that have not yet been assigned to subfamilies will likely define new subfamilies as the pool of available sequences continues to increase. In the future, these subfamilies will be gradually released with new identifiers when they have been sufficiently populated. Compared to the GH13 subfamily classification , the GH5 subdivision has resulted in both a higher number of subfamilies and a higher number of uncharacterized subfamilies, suggesting that family GH5 is comparatively less well explored.
To further refine the global subfamily analysis, maximum likelihood phylogenetic analysis has been performed on each subfamily (Additional file 1: Figures S1: GH5_1 – GH5_53). In addition to database accession numbers, information about substrate specificity (as provided by EC numbers) has been included, and 3-D structures highlighted in each subfamily tree. The 51 subfamilies (numbered GH5_1 to GH5_53 as explained elsewhere) have been categorized based on available enzyme activity data, to aid their individual descriptions, below. Thus, the first group, comprised of “Enzymatically-Active Subfamilies”, contains subfamilies with at least one member whose enzyme activity has been shown. The extent of the documented characterization varies substantially within the group, from simple information obtained from enzymatic assays insufficient to assign a particular EC number (these enzymes are denoted with a star in Additional file 1: Figure S1), to detailed enzyme specificity and kinetics studies. Among the better characterized subfamilies, monospecific subfamilies are distinguished by the presence of only one EC number for one or more members, whereas in polyspecific subfamilies, two or more enzymatic activities have been observed in different members. “Uncharacterized Subfamilies” comprise the second major group; these subfamilies currently lack documented enzymatic activity altogether.
In this context, it is worth noting that the majority of the large, well characterized subfamilies were polyspecific. These highly populated subfamilies were also the first ones to be identified and described, and often, despite observed polyspecificity, one particular activity predominates. For example, 22 of 24 characterized enzymes in GH5_1 are reported to be endo-glucanases, whereas one protein is a documented as a cellobiohydrolase, and another was described as displaying licheninase activity. It is, however, difficult to draw far-reaching conclusions based on these observations. On one hand, it may be that the acquisition of a new specificity within a subfamily is a rare event; alternatively, the observation of one or few activities in specific subfamilies may simply be a consequence of differences in the range of substrates tested experimentally.
Another significant aspect of many GH5 family proteins is that their protein sequence may include additional modules with different functions, and in particular CBMs . A great variety of modular structures may be found throughout the family and in a number of individual subfamilies. An analysis of all the complexity of modular structures found in the family goes beyond the objectives of this study, and some aspects of this diversity are illustrated in Figure 2. For instance, many members of subfamily GH5_8 are modular and reveal two major trends: (i) the addition of one or of multiple CBMs (see Figure 2a) is more common and may be associated not only to the nature of the main substrate of the corresponding catalytic domain, particularly in complex substrates; and (ii) the combinations with other catalytic modules to form bifunctional enzymes (see Figure 2b but also 2c), are more rare but particular useful to reveal interacting or synergistic enzyme activities of some catalytic modules. Numerous modular arrangements can also be found in other large subfamilies like GH5_1, GH5_2 and GH5_7. CBMs can be located on the N- or C-terminal side of the GH5 module (as illustrated in GH5_8 in Figure 2a). The combination of catalytic domains may target different tissue components. Some may, for instance, target cellulose and cellulose associated substrates (as in Figure 2b) but bifunctional enzymes likely targeting hemicellulose may also be found (several examples in Figure 2c).
Subfamilies with identified active enzymes
A number of subfamilies exhibit a single activity among their characterized enzyme members. Multiple individual examples within a subfamily improve the degree of confidence regarding subfamily monospecificity, while subfamilies with only a single characterized representative may be subject to reinterpretation in the future as the breadth of biochemical data increases.
GH5_5 endo-β-1,4-glucanases (EC 188.8.131.52)
The largest subfamily that contains only a single EC number is GH5_5, which is primarily composed of secreted bacterial and fungal enzymes (Additional file 1: Figure S1: GH5_5). All investigated enzymes in GH5_5 display endo-β-1,4-glucanase activity (EC 184.108.40.206). One crystal structure has been determined for a Thermoascus aurantiacus endo-glucanase . In fungi, about half of the GH5 proteins harbor a CBM1 module at the N- or C-terminus, which is compatible with an active role on cellulose. No modular proteins are found among the bacterial members of the subfamily.
GH5_8, GH5_10 and GH5_17: endo-β-1,4-mannanase (EC 220.127.116.11)
Members of the large GH5_8 subfamily are all extracellular mannan endo-β-1,4-mannosidases (EC 18.104.22.168) according to available biochemical characterization, and this subfamily was historically described as the bacterial mannanase subfamily A8 . Structural analysis has highlighted distinctive features of alkaline β-mannanases . The subfamily now contains a single eukaryotic enzyme from the beetle Hypothenemus hampai, resulting from a horizontal gene transfer from bacteria .
In the closely related but distinct subfamilies GH5_10 and GH5_17, only extracellular enzymes with endo-β-1,4-mannanase activity (EC 22.214.171.124) have equally been reported. Metazoan sequences, including sequences originating from mollusks and arthropods, as well as bacterial sequences compose the GH5_10 subfamily, whereas subfamily GH5_17 only harbors bacterial enzymes. Currently, none of the bacterial GH5_10 members have a documented enzyme activity, although a Fibrobacter succinogenes enzyme is active on AZCL-galactomannan . The bacterial enzymes in GH5_17 are all cellulosomal components from the genus Clostridium.
GH5_14 : (plant) exo-β-1,3-glucosidase (EC 126.96.36.199)
Subfamily GH5_14 comprises exclusively plant enzymes and relies on a single functional characterization at present. Recombinant expression of the rice OsGH5BG showed that the enzyme had glucan β-1,3-glucosidase activity (EC 188.8.131.52) . A unique feature of GH5_14 is the fascin-like module inserted after β-strand 1. Fascin is a human actin-binding protein, but the function of the plant fascin-like domain is unknown. Members of GH5_14 are well represented throughout the plant kingdom, but most interestingly, a representative is absent in the leading plant model organism Arabidopsis thaliana.
GH5_15 : (fungal) (endo-)β-1,6-glucanases (EC 184.108.40.206)
Subfamily GH5_15 is a small but well characterized subfamily composed of secreted fungal enzymes. The identified β-1,6-glucanase activity (EC 220.127.116.11) is important for the mycoparasitic activity and probably cell wall recycling by some fungi. The GH5_15 phylogenetic tree displays two major clades (Additional file 1: Figure S1: GH5_15). The largest clade is formed by enzymes issued from fungi from the class of Sordariomycetes. Eurotiomycetes are found the second subgroup.
GH5_16 : (fungal) endo-β-1,6-galactanase (EC 18.104.22.168)
GH5_16 is another example of a monospecific subfamily of secreted enzymes where only a single fully sequenced biochemically characterized β-1,6-galactanase (EC 22.214.171.124) is currently known , although partial N-terminal sequence of an Aspergillus enzyme, closely related to sequences of other members of the subfamily, yields the same activity . Several other enzymes with β-1,6-galactanase activity have been moved from GH5 to GH30_5 , but subfamily GH5_16 clearly remains within family GH5. The known β-1,6-galactanase (EC 126.96.36.199) is involved in larch wood arabinogalactan degradation .
GH5_21 : (bacteroidetes) endo-β-1,4-xylanase (EC 188.8.131.52)
Xylanase (EC 184.108.40.206) activity has been recently established  for a number of xylanolytic bacteroidetes enzymes belonging to this subfamily. These GH5_21 endo-xylanases integrate xylan utilization gene clusters found in Prevotella and Bacteroides species and are all apparently secreted. Significant differences in their mode of action have been observed, despite the inclusion in the same subfamily. Different GH5_21 enzymes were shown to release different products from wheat arabinoxylan .
GH5_22 GH5_31, GH5_34, GH5_39, and GH5_53 : single β-glycanase characterizations
In addition to subfamilies GH5_14 and GH5_16 described above, five other subfamilies are distinguished by harboring only a single experimentally characterized enzyme. Endo-β-1,4-glucanase activity (EC 220.127.116.11) has been determined in subfamilies GH5_22 and GH5_39. GH5_31 is a small subfamily currently restricted to secreted fungal proteins. A β-1,3-(gluco)mannanase activity (EC 3.2.1.-)a has been reported for an enzyme from Paecilomyces lilacinus. In the small subfamily GH5_34, composed of extracellular modular proteins from bacterial origin, there is a single enzymatically and structurally characterized enzyme from Clostridium thermocellum. Notably, although this is the first reported enzyme with arabinoxylanase activity (EC 3.2.1.-)b it was designated CtXyl5A in spite of its inability to attack unsubstituted xylans . Finally, the small subfamily GH5_53 is a modular extracellular subfamily that contains a single characterized cellodextrinase (EC 18.104.22.168) .
GH5_27 GH5_28 and GH5_29 : endo-glycosylceramidases (EC 22.214.171.124)
GH5_27, GH5_28 and GH5_29 are subfamilies exclusively containing extracellular endo-glycosylceramidases. Subfamily GH5_27 is formed of sequences of eukaryotic origin while the small subfamilies GH5_28 and GH5_29 are bacterial. Interestingly, the first subfamily is found among the four subfamilies that contain metazoan GH5 enzymes. All but one of the enzymes found in GH5_28 are from Actinobacteria. The crystal structure of a Rhodococcus endo-glycoceramidase revealed an active site channel atypical for GH5 enzymes, which explains the unusual substrate for this type of enzymes . All the GH5_29 sequences reported here are from the genus Rhodococcus. The characterized bacterial enzymes of GH5_28 hydrolyze ganglio- and lacto-series glycosphingolipids. In contrast, the only GH5_29 enzyme investigated is not capable of hydrolyzing these substrates. Instead, this enzyme shows activity against 6-gala series glycosphingolipids and the designation oligogalactosyl-N-acylsphingosine 1,1’-β-galacto-hydrolase has been proposed .
GH5_52 : cellodextrinases (EC 126.96.36.199)
Two enzymes in the small intracellular bacterial GH5_52 subfamily exhibited cellodextrinase activity (EC 188.8.131.52). In addition two enzymes isolated from the cow rumen have been shown to have hydrolytic activity on carboxymethyl cellulose (CMC) agar .
Some GH5 subfamilies group together a panel of activities and are described as polyspecific subfamilies. The apparent plasticity of these subfamilies suggests that only few subtle changes could be sufficient to switch from one activity to the other. More likely, many of the enzymes present on the subfamily level are polyspecific to some extent and therefore should have more than a single EC attributed.
GH5_1 and GH5_2 : β-1,4-glucan cleaving enzymes
Extracellular enzymes from archaea, bacteria and uncultured symbiotic protists are represented in subfamily GH5_1 (Additional file 1: Figure S1: GH5_1). The activity observed for most characterized GH5_1 enzymes is endo-β-1,4-glucanase activity (EC 184.108.40.206). The apparent exception in subfamily GH5_1 is an exo-acting cellobiohydrolase activity (EC 220.127.116.11) from Clostridium thermocellum. However, one should note that the biochemical distinction of exo- versus endo-acting cellulases is particularly difficult to establish experimentally. Interestingly, an enzyme from Ruminococcus albus able to cleave CMC and glucomannan but particularly active on lichenin was recently described . Many proteins in this subfamily are modular (data not shown), a feature shared with many members of subfamily GH5_8 as described previously.
Subfamily GH5_2 is currently the largest in family GH5. This subfamily of extracellular enzymes, many of which are multimodular, contains a large number of characterized members that display endo-β-1,4-glucanase activity (EC 18.104.22.168). These endo-glucanases are distributed across the subfamily tree and are found in every major clade Additional file 1: Figure S1: GH5_2). One endo-glucanase from Fibrobacter succinogenes S85 in this subfamily was reported to be active both on CMC and oat spelt xylan , but xylanase activity has not been observed in other members thus far. Interestingly, one representative of this subfamily has been reported as a chitosanase (EC 22.214.171.124) with transglycosylation activity . A bifunctional cellulase/chitosanase has also been identified in Bacillus sp. NBL420  in a different clade. Significantly, a closely related N-terminal sequence of a bifunctional cellulase/chitosanase from Myxobacter sp. AL-1, suggests that a specific subgroup bearing both activities may exist. As for subfamily GH5_1, many members of this subfamily are multimodular, having both CBMs and cellulosomal-like dockerins.
GH5_4 : (xyloglucan-specific) endo-β-1,4-glucanases (EC 126.96.36.199 and EC 188.8.131.52), licheninases (EC 184.108.40.206), and xylanases (EC 220.127.116.11)
Subfamily GH5_4 members are typically extracellular bacterial enzymes, although some members come from ciliates and fungi, predominantly rumen organisms. In total, four enzyme activities have been reported for subfamily GH5_4. Thus far, GH5_4 is the only subfamily containing enzymes with reported xyloglucanase activity (EC 18.104.22.168) . Interestingly, the xyloglucanases are found in two different clades of the GH5_4 subfamily tree suggesting that the switch of enzyme activity inside this subfamily occurred at different times. (Additional file 1: Figure S1: GH5_4). Licheninases (EC 22.214.171.124) have also been described in this subfamily. Significantly, endo-β-1,4-xylanase activities (EC 126.96.36.199) have been reported for a few enzymes, but always in conjunction with other activities. For example, the xylan degrading specific activity of Clostridium cellulovorans EngB and EngD are both of approximately 14% of their respective specific activities on lichenan . Such features suggest an important degree of enzyme promiscuity given the structural similarity of the β-linked substrates. However, the most commonly reported EC number for representatives of GH5_4 is EC 188.8.131.52. Except for a few fungal pathogen members of the subfamily that bear a CBM1 module, no other known CBM is present, in sharp contrast to what has been found for extracellular subfamilies GH5_1 and GH5_2.
GH5_7 : β-1,4-mannan-cleaving enzymes (EC 184.108.40.206 and EC 220.127.116.11)
In subfamily GH5_7, previously named A7 (originally comprised of only eukaryotic mannanases but now also containing archaeal and bacterial members) , the three reported enzyme activities are associated to the degradation or to the modification of β-mannan-containing polysaccharides. Virtually all examined GH5_7 enzymes possess endo-β-1,4-mannanase activity (EC 18.104.22.168). One exception is the tomato protein LeMan4a, which in addition to hydrolytic activity, can act in vitro as a mannan transglycosylase (EC 2.4.1-) . Recently, mannan transglycosylase activity was also reported for two fungal GH5_7 enzymes . It is not impossible that further GH5_7 enzymes may reveal a transglycosylase activity in the future, since the distinction between hydrolytic and transglycosylase activity is dictated by the tendency of the glycosyl-enzyme intermediate to be intercepted by a water molecule or a saccharide molecule, respectively. Another interesting case among GH5_7 members is the β-mannosidase CmMan5A, which is able to release mannose from the non-reducing end of manno-oligosaccharides and -polysaccharides (EC 22.214.171.124) and is thus exo- and not endo-acting. This difference was explained by the length of three loops which modify the active center accessibility . Furthermore, this subfamily has both extracellular and intracellular enzymes. Many GH5_7 extracellular enzymes from fungi contain a CBM. In many bacteria, CBMs from different families are also found appended to GH5_7 catalytic modules (data not shown).
GH5_9 : fungal cell wall modifying enzymes
This subfamily contains only sequences of fungal origin putatively found in different cell locations: some are secreted, several present a GPI-anchor, others have single transmembrane segments and yet others appear to be intracellular. The activity exo-β-1,3-glucanase (EC 126.96.36.199) has been described for most of the characterized examples, all found among the apparently secreted enzymes. Structural investigation of the Candida albicans Exg protein provides a clue to the evolution of these exo-hydrolases; the structure reveals an active site pocket shaped for cleavage of β-1,3- but not β-1,4 glycosidic linkages . Surprisingly, the periplasmic Exg1 protein from Schizosaccharomyces pombe was demonstrated to be an endo-β-1,6-glucanase (EC 188.8.131.52), while the membrane-anchored protein Exg2 protein was shown to produce cell wall material when over-expressed . Two of the characterized GH5_9 exo-β-1,3-glucanases have been also described as β-glucosidases (EC 184.108.40.206) .
GH5_12 : β-glucosylceramidases (EC 220.127.116.11) and (flavonoid) β-glucosidase (EC 18.104.22.168)
Previously known to contain a flavonoid β-glucosidase (EC 22.214.171.124) from yeast , the panorama of specificities in this subfamily recently expanded to include several fungal β-glucosylceramidases (EC 126.96.36.199). The subfamily is grouped into a fungal clade, including three Cryptococcus sequences (one now characterized), and a clade dominated by bacterial enzymes (Additional file 1: Figure S1: GH5_12).
GH5_23 : fungal β-diglycosidases (EC 188.8.131.52 and EC 184.108.40.206)
Subfamily GH5_23 is composed of secreted fungal proteins, two of which have been characterized as β-diglycosidases that break down plant diglycoconjugated flavonoids. Deglycosylation of these compounds most often involves the sequential action of two β-glycosidases in contrast to the one-step hydrolytic release of the disaccharide moiety from the aglycone by β-diglycosidases . The characterized enzymes are a hesperidin 6-O-α-L-rhamnosyl-β-glucosidase (EC 220.127.116.11) from Stilbella fimetaria from which a partial sequence has been obtained  and a reported β-primeverosidase (EC 18.104.22.168) from Penicillium multicolor TS-5 . This subfamily is present in the genera Aspergillus and Penicillium known for interaction with plants.
GH5_25 : endo-β-1,4-glycanases (EC 22.214.171.124 and EC 126.96.36.199)
Most enzymes found in subfamily GH5_25 are derived from thermophiles. Interestingly, characterized enzymes in this subfamily represent examples of GH5 enzymes with multiple activities. For instance, Cel5A from Thermotoga maritima exhibits activity on both β-mannan-based and β-glucan-based polymers. Analyses of the TmCel5A structure have highlighted features important both for the nature of the duality and the thermostability [58, 59].
GH5_26 : endo-β-1,4-glycanases (EC 188.8.131.52 and EC 184.108.40.206)
GH5_26 is a small subfamily with a majority of sequences from uncultured microorganisms. The dominating activity found is endo-β-1,4-glucanase (EC 220.127.116.11), but one enzyme has also high activity against lichenan (EC 18.104.22.168) .
GH5_36 endo-b-1,4-glycanases (EC 22.214.171.124 and EC 126.96.36.199)
Several bacterial phyla are represented in subfamily GH5_36. One enzyme has a demonstrated endo-β-1,4-mannanase activity (EC 188.8.131.52), whereas a second enzyme exhibits licheninase activity (EC 184.108.40.206) in addition to β-mannanase activity [61, 62]; a 3-D structure is available for the latter enzyme.
GH5_37 : endo-β-1,3/4-glycanases (EC 220.127.116.11 and EC 18.104.22.168) + cellodextrinase (EC 22.214.171.124)
Three different activities are found in subfamily GH5_37, which consists of sequences of bacterial origin encoding intracellular proteins. The majority of the characterized enzymes are cellulases (EC 126.96.36.199), but there are also examples of licheninase activity (EC 188.8.131.52) , and cellodextrinase activity (EC 184.108.40.206) .
GH5_38 and GH5_46 : bacterial enzymes active on model plant cell wall (PCW) compounds
In subfamily GH5_38, three enzymes isolated from rumen metagenomic projects have been partially characterized [12, 14]. To these, we may add a Prevotella ruminicola 23 enzyme described as cellulase (PRU_1856) and shown to be active against CMC, Avicel, and lichenan . Another enzyme discovered in the microbial community of the cow rumen is active on CMC and belongs to subfamily GH5_46 . This is the single evidence of activity in this subfamily.
GH5_48 : bacterial enzymes active on chitin and chitosan derivatives
Several bacterial phyla are currently represented in subfamily GH5_48, mostly composed of extracellular and membrane-anchored proteins of unknown function. A single member, a partial sequence (not shown in trees because incomplete) coding for a protein from Pseudomonas putida P3(4)), over 98% identical to locus PPS_1333 (GenBank accession AEJ11908) from Pseudomonas putida S16, has been described as a bifunctional enzyme as it was active on preparative forms of chitin and colloidal chitosan and on the model compounds pNP-β-N-acetylglucosaminide and 4-methylumbelliferyl-N-acetyl β-D-glucosaminide . Given that chitosanases are already present in family GH5, related activities are not unexpected. Interestingly, one sequence in GH5_48 seems to lack the catalytic machinery (Additional file 1: Figure S1).
Subfamilies lacking experimental characterization
A large number of subfamilies still require evidence of enzyme activity (Table 1). A large characterization effort is still needed to identify the hidden activities of these subfamilies. For a number of subfamilies, particularly for those containing bacterial enzymes, direct hints of the activities to be tested can be obtained from the additional modules present, such as CBMs, but also from operon-like organizations. But the contribution of other more indirect approaches like transcriptomics and functional metagenomics using innovative sets of conditions and enlarged classes of substrates are equally potent to provide clues to enzyme function. The uncharacterized subfamilies are distributed across the GH5 phylogenetic tree. However, in one of the three major clades, a large block of uncharacterized subfamilies (GH5_13, GH5_18, GH5_19, GH5_30, GH5_41 and GH5_42) is located between characterized subfamilies containing β-mannan-acting enzymes (GH5_7, GH5_10, GH17, and GH5_31) (Figure 1). It is tempting to speculate that at least some of these uncharacterized subfamilies target substrates related to β-mannan.
Unclassified GH5 proteins
Roughly 20% of the analyzed GH5 sequences were not assigned to subfamilies, although some of these proteins have been characterized to different levels. The main reason for the inability to assign them to subfamilies was the lack of sufficiently related sequences to define a subfamily with at least five members. This is likely to change in the future due to a predicable increase in the number of available sequences and these unclassified sequences represent a pool out of which many more subfamilies will emerge. A total of 17 characterized enzymes were identified in this diverse set of non-classified sequences (see Table 2). Interestingly, two main groups arise: (i) post-genomic enzyme characterizations, and (ii) enzymes originating from functional metagenomics-based discovery and characterization. Here, the cytoplasmic endo-β-1,6-glucanase Exg3 from the model organism Schizosaccharomyces pombe is a representative of the former efforts focused on the identification of fundamental activities. On the other hand, the archael multidomain hyperthermophilic cellulase EBI-244 screened for its ability to degrade CMC at high temperature represents the latter category of efforts .
Non-catalytic GH5 subfamilies and GH5 modules
Interestingly, the analysis subjacent to the subfamily classification revealed proteins that are likely to be catalytically inactive (or perhaps possess alternate mechanisms and activities to the canonical GHs), due to incomplete catalytic machinery (see Figure 2d). GH5_30 is the only subfamily where all sequences are lacking the essential amino acids for GH activity. In addition, some members of other groups also appear to have lost their catalytic machinery, such as a Xanthomonas-specific subgroup that appears to rapidly emerge from subfamily GH5_1 (Additional file 1: Figure S1). All three members of this subgroup present an architecture where the inactive GH5 module is appended to an expansin module and an adjacent CBM63 module at the C-terminus. Although no catalytic chemical activity has been identified for expansin modules, it is significant that the knockout of CelXoB renders Xanthomonas oryzae pv. oryzae KACC 10331 avirulent . Another emerging non-catalytic GH5 subfamily is closely related to subfamily GH5_8. All three members of this new subgroup are lipoproteins that combine the apparently catalytically inactive GH5 module with a large C-terminal extension. These two subgroups share long branches in the common GH5 tree suggesting rapid evolution. Interestingly, a single member of subfamily GH5_48 has lost the catalytic machinery. Besides this loss, it is however still similar to other members of the subfamily. Whether this is a snapshot of the early steps of a new arising function or a sequencing error it is premature to say. Finally, the only distant GH5 member having lost its catalytic acid–base that has been functionally characterized is the putative carbohydrate biosensor Rsi24C-GH5 from Clostridium thermocellum ATCC 27405. Although its catalytic activity was lost, the extracellular GH5 module was shown to interact with crystalline cellulose so that a recognition signal could be conveyed by to its intracellular N-terminal anti-σ factor. The losses of the catalytic machinery here described convey that family GH5 sequences are also subject to recurring evolution that leads to novel functions. This type of evolutionary event has been described previously in other GH families. For instance, amino acid transporters derived from ancestral α-amylases are found in family GH13 , inactivated chitinases evolved into xylanase inhibitors in family GH18 , and mammalian lactalbumins, which are related to GH22 lysozymes, are all well-known examples of the recent evolution of glycosidases to acquire novel functionalities .
When the first five historical subfamilies of GH5 were established in 1990, the total number of GH5 protein sequences was 21 . More than 20 years later (August 2012), this number has increased approximately 150 times to exceed 3,200. The practical difficulties of handling such large datasets notwithstanding, this abundance of sequences is both a boon and a bane for phylogenetic analysis and functional prediction.
Assigning proteins to a large GH-family, like GH5, which harbors multiple specificities and activities, does not unlock the full potential of sequence-based classification. Thus, one aim of the present investigation was to obtain an improved correlation between protein sequences and catalytic specificity by refining a finer hierarchical level, the subfamily, for GH5 members. Up to 80 percent of the existing GH5 sequences were segregated into 51 subfamilies. Of these subfamilies, a total of 31 contained at least one member characterized to some degree, whereas 20 lacked enzymatically-characterized members altogether. Out of the 31 subfamilies characterized to some extent, 17 were monospecific and eleven were polyspecific (containing two or more enzyme activities). Nonetheless, one activity typically predominates within polyspecific subfamilies. Interestingly, both endo- and exo-acting enzymes have been observed in the same subfamily, e.g. GH5_1, GH5_7 and GH5_9, illustrating that (if real) the two types of activities reflect details of the three-dimensional structures. As a consequence of the canonical double displacement mechanism employed by GH5 enzymes, which involves the formation of a covalent glycosyl-enzyme intermediate, GH5 members can potentially catalyze transglycosylation in addition to, or instead of, hydrolysis . Although the amount of biochemical data is presently limited, we observed that subfamily classification in GH5 does not appear to correlate with transglycosylation activity, thus indicating that this property is also a consequence of subtle protein structural details.
This effort provides a first comprehensive view of the coverage and distribution of the curated set of 400 experimentally characterized enzymes in the GH5 family and couples this information with an extensively updated sequence-based GH5 subfamily division. In particular, it provides insights into the evolution of GH5 proteins, and the classification results can be used to assist in candidate protein selection for enzyme discovery and bioprospecting projects. For instance, both the members from the twenty defined subfamilies lacking functional characterization, as well as the numerous phylogenetic outliers, provide a vast number of interesting targets for future studies. In particular, although a significant amount of tertiary structural data is already available for GH5, the present work highlights that a large number subfamilies would benefit from a 3-D structure for at least one subfamily member. Moreover, the data presented here, and available at the CAZy database  as a community resource, will serve as a guide for protein engineering approaches exploiting the diverse activities found within the GH5 family.
Finally, in the present climate in which sequence data is literally flooding public databases, incorrect protein function annotations are too easily propagated by automated computer-based prediction methods, thereby jeopardizing the usefulness of these annotations. Increasingly rapid sequence accumulation is worsening the scenario. This problem is particularly illustrated by the GH5_11 subfamily of plant and fungal proteins that are annotated as cellulases in public databases, including the widely-used Arabidopsis Information Resource , despite a complete lack of experimental support for any one of its members. Such excesses of over-annotation equally affect the presumed non-catalytically active proteins and subfamilies. For example, subfamily GH5_30 further exemplifies the pitfalls of automated (mis-)annotation: several members are publicly annotated as mannanases, although they lack the conserved catalytic machinery of the family. To avoid such error propagation, we strongly advocate designating all predicted enzymes as “GH5_n” (where n is the subfamily number) until an activity has been rigorously demonstrated by biochemical experimentation.
The GH5 subfamily classification presented here provides a framework to sort family members into meaningful, predictive categories. By taking a conservative approach to protein annotation, this method offers a rigorous strategy to avoid misleading functional prediction in large-scale genomic sequencing projects. Whilst the subfamilies described herein generally act on a single substrate (seventeen monospecific subfamilies were identified), it is important to stress that precise details of glycoside hydrolase function, such as the extent of endo- vs. exo- modes of cleavage or the transglycosylation-to-hydrolysis ratio is unlikely to be predictable from sequence alone. We therefore recommend that such over-reaching predictions be altogether abandoned in genomic sequence annotation. To aid and advance global efforts in de novo sequence annotation, the GH5 subfamily classification scheme is now publicly available at the CAZy database .
The GH5 subfamilies were defined using the methods described for the subfamily classification of GH13 and all PL families [24, 25], which are briefly summarized here. After an initial removal of obviously incomplete and/or erroneous sequences, a total of 2347 full length GH5 catalytic module sequences were retrieved from the CAZy database (October 2011) and subdivided into two sets. One set of 414 sequences contained all GH5 modules from biochemically characterized and sequences positively tested in activity tests against a variety of substrates. The second set was composed of 1957 non-characterized sequences The latter subset was clustered at 75% identity using UCLUST4.0, a part of the USEARCH 4.0 package , and was reduced to 971 sequences. When combined with the 414 GH5 module sequences from a biochemically characterized and positively assayed subset we obtained a total of 1385 sequences. These sequences were aligned using MUSCLE 3.7  in two steps. An initial alignment was performed and its quality visually inspected so that the remaining incomplete and problem sequences were identified and edited or removed. This procedure ensured that the GH5 module boundaries were clear and that the sequences were trimmed if required, and that a majority of the residues constituting the catalytic site was present or that the alignment was not ambiguous. A final set of 1367 remaining sequences were then realigned using the same procedure. The eliminated, often fragmentary, sequences were used to complement biochemical activity information when relevant at later stages.
The resulting multisequence alignment of family GH5 catalytic domain sequences was used to infer an approximate-maximum-likelihood phylogenetic tree with FASTTREE 2.1 , a program adapted to the analysis of large sequence sets, using the Whelan Goldman model of amino acid evolution, the gamma option to rescale the branch lengths and compute a Gamma20-based likelihood, a total of four rounds of minimum-evolution moves, and options to make the maximum-likelihood nearest-neighbor interchanges more exhaustive. The identification and tagging of subfamilies followed a multi-step procedure. First, the tree was analyzed to tag the different sequences and nodes corresponding to the first 10 “historical” subfamilies (A1 to A10, described in the Introduction). These steps were performed to ensure continuity in subfamily definitions. For the remaining sequences in the tree, distinct nodes corresponding potential subfamilies were visually identified and their consistency checked using a procedure similar to that described . Within each of these groups, different starting sequences were selected and manual BLAST2 queries  performed against all the sequences found in CAZy in order to identify self-contained ensembles and establish subfamily limits. This analysis ensured that the each subfamily that was retained was singular and that the removal of sequences by the initial UCLUST filtering procedure and of fragmentary sequences did not introduce any bias. Finally, only subfamilies containing at least five sequences from different organisms were considered.
Subsequently, for each defined subfamily, maximum likelihood (ML) phylogenetic trees were built by using PhyML , and the reliability of the inferred relationships the trees was tested by bootstrap analysis using 100 resamplings of the data set.
aNo EC number is presently available for β-1,3-(gluco)mannanase activity.
bNo EC number is presently available for arabinoxylanase activity.
Varki A: Essentials of glycobiology. 2009, Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory Press, 2
Lichtenthaler FW: Carbohydrates as renewable raw materials: a major challenge of green chemistry. Methods and reagents for green chemistry: an introduction. Edited by: Tundo P, Perosa A, Zecchini F. 2007, Hoboken, NJ: J. Wiley, 23-63.
Laine RA: A calculation of all possible oligosaccharide isomers both branched and linear yields 1.05x10(12) structures for a reducing hexasaccharide: the isomer barrier to development of single-method saccharide sequencing or synthesis systems. Glycobiology. 1994, 4: 759-767. 10.1093/glycob/4.6.759.
Sharon N: The conquest of the last frontier of molecular and cell biology. Foreword Biochimie. 2001, 83: 555-555. 10.1016/S0300-9084(01)01310-4.
Anonymous: 10 Emerging Technologies That Will Change the World. Technology Review. 2003, 33-51.
Davies GJ, Gloster TM, Henrissat B: Recent structural insights into the expanding world of carbohydrate-active enzymes. Curr Opin Struct Biol. 2005, 15: 637-645. 10.1016/j.sbi.2005.10.008.
Nomenclature Committee of the International Union of Biochemistry and Molecular Biology (NC-IUBMB). http://www.chem.qmul.ac.uk/iubmb/enzyme/,
Henrissat B: A classification of glycosyl hydrolases based on amino-acid-sequence similarities. Biochem J. 1991, 280: 309-316.
Davies GJ, Sinnott ML: Sorting the diverse: the sequence-based classifications of carbohydrate-active enzymes. Biochem J. 2008, 1-5. 10.1042/BJ20080382. (online only)
Henrissat B, Claeyssens M, Tomme P, Lemesle L, Mornon JP: Cellulase families revealed by hydrophobic cluster-analysis. Gene. 1989, 81: 83-95. 10.1016/0378-1119(89)90339-9.
The Carbohydrate-Active enZYme (CAZy) database. http://www.cazy.org,
Duan CJ, Xian L, Zhao GC, Feng Y, Pang H, Bai XL, Tang JL, Ma QS, Feng JX: Isolation and partial characterization of novel genes encoding acidic cellulases from metagenomes of buffalo rumens. J Appl Microbiol. 2009, 107: 245-256. 10.1111/j.1365-2672.2009.04202.x.
Elifantz H, Waidner LA, Michelou VK, Cottrell MT, Kirchman DL: Diversity and abundance of glycosyl hydrolase family 5 in the North Atlantic Ocean. FEMS Microbiol Ecol. 2008, 63: 316-327. 10.1111/j.1574-6941.2007.00429.x.
Hess M, Sczyrba A, Egan R, Kim TW, Chokhawala H, Schroth G, Luo SJ, Clark DS, Chen F, Zhang T, et al: Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011, 331: 463-467. 10.1126/science.1200387.
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.
Béguin P: Molecular biology of cellulose degradation. Annu Rev Microbiol. 1990, 44: 219-248. 10.1146/annurev.mi.44.100190.001251.
Lo Leggio L, Parry NJ, VanBeeumen J, Claeyssens M, Bhat MK, Pickersgill RW: Crystallization and preliminary X-ray analysis of the major endoglucanase from Thermoascus aurantiacus. Acta Crystallogr D Biol Crystallogr. 1997, 53: 599-604. 10.1107/S0907444997005404.
Hilge M, Gloor SM, Rypniewski W, Sauer O, Heightman TD, Zimmermann W, Winterhalter K, Piontek K: High-resolution native and complex structures of thermostable beta-mannanase from Thermomonospora fusca - substrate specificity in glycosyl hydrolase family 5. Structure Fold Des. 1998, 6: 1433-1444. 10.1016/S0969-2126(98)00142-7.
Lo Leggio L, Larsen S: The 1.62 angstrom structure of Thermoascus aurantiacus endoglucanase: completing the structural picture of subfamilies in glycoside hydrolase family 5. FEBS Lett. 2002, 523: 103-108. 10.1016/S0014-5793(02)02954-X.
Larsson AM, Anderson L, Xu BZ, Munoz IG, Uson I, Janson JC, Stalbrand H, Stahlberg J: Three-dimensional crystal structure and enzymic characterization of beta-mannanase Man5A from blue mussel Mytilus edulis. J Mol Biol. 2006, 357: 1500-1510. 10.1016/j.jmb.2006.01.044.
Costanzo S, Ospina-Giraldo MD, Deahl KL, Baker CJ, Jones RW: Alternate intron processing of family 5 endoglucanase transcripts from the genus Phytophthora. Curr Genet. 2007, 52: 115-123. 10.1007/s00294-007-0144-z.
Opassiri R, Pomthong B, Akiyama T, Nakphaichit M, Onkoksoong T, Cairns MK, Cairns JRK: A stress-induced rice (Oryza sativa L.) beta-glucosidase represents a new subfamily of glycosyl hydrolase family 5 containing a fascin-like domain. Biochem J. 2007, 408: 241-249. 10.1042/BJ20070734.
St John FJ, Gonzalez JM, Pozharski E: Consolidation of glycosyl hydrolase family 30: A dual domain 4/7 hydrolase family consisting of two structurally distinct groups. FEBS Lett. 2010, 584: 4435-4441. 10.1016/j.febslet.2010.09.051.
Stam MR, Danchin EGJ, 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.
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.
Hancock SM, Rich JR, Caines MEC, Strynadka NCJ, Withers SG: Designer enzymes for glycosphingolipid synthesis by directed evolution. Nat Chem Biol. 2009, 5: 508-514. 10.1038/nchembio.191.
Danchin EGJ, Rosso MN, Vieira P, de Almeida-Engler J, Coutinho PM, Henrissat B, Abad P: Multiple lateral gene transfers and duplications have promoted plant parasitism ability in nematodes. Proc Natl Acad Sci USA. 2010, 107: 17651-17656. 10.1073/pnas.1008486107.
Acuña R, Padilla BE, Flórez-Ramos CP, Rubio JD, Herrera JC, Benavides P, Lee S-J, Yeats TH, Egan AN, Doyle JJ, et al: Adaptive horizontal transfer of a bacterial gene to an invasive insect pest of coffee. Proc Natl Acad U S A. 2012, 109: 4197-4202.
Boraston AB, Bolam DN, Gilbert HJ, Davies GJ: Carbohydrate-binding modules: fine-tuning polysaccharide recognition. Biochem J. 2004, 382: 769-781. 10.1042/BJ20040892.
Zhao YJ, Zhang YH, Cao Y, Qi JX, Mao LW, Xue YF, Gao F, Peng H, Wang XW, Gao GF, et al: Structural analysis of alkaline beta-mannanase from alkaliphilic Bacillus sp N16-5: implications for adaptation to alkaline conditions. PLoS One. 2011, 6: e14608-10.1371/journal.pone.0014608.
Suen G, Weimer PJ, Stevenson DM, Aylward FO, Boyum J, Deneke J, Drinkwater C, Ivanova NN, Mikhailova N, Chertkov O, et al: The complete genome sequence of Fibrobacter succinogenes S85 reveals a cellulolytic and metabolic specialist. PLoS One. 2011, 6: e18814-10.1371/journal.pone.0018814.
Sakamoto T, Taniguchi Y, Suzuki S, Ihara H, Kawasaki H: Characterization of Fusarium oxysporum beta-1,6-galactanase, an enzyme that hydrolyzes larch wood arabinogalactan. Appl Environ Microbiol. 2007, 73: 3109-3112. 10.1128/AEM.02101-06.
Luonteri E, Laine C, Uusitalo S, Teleman A, Siika-aho M, Tenkanen M: Purification and characterization of Aspergillus beta-D-galactanases acting on beta-1,4- and beta-1,3/6-linked arabinogalactans. Carbohyd Polym. 2003, 53: 155-168. 10.1016/S0144-8617(02)00303-X.
Dodd D, Moon YH, Swaminathan K, Mackie RI, Cann IKO: Transcriptomic analyses of xylan degradation by Prevotella bryantii and insights into energy acquisition by xylanolytic bacteroidetes. J Biol Chem. 2010, 285: 30261-30273. 10.1074/jbc.M110.141788.
Sugino H, Furuichi S, Murao S, Arai M, Fujii T: Molecular characterization of a Rhodotorula-lytic enzyme from Paecilomyces lilacinus having beta-1,3-mannanase activity. Biosci Biotechnol Biochem. 2004, 68: 757-760. 10.1271/bbb.68.757.
Correia MAS, Mazumder K, Bras JLA, Firbank SJ, Zhu YP, Lewis RJ, York WS, Fontes CMGA, Gilbert HJ: Structure and function of an arabinoxylan-specific xylanase. J Biol Chem. 2011, 286: 22510-22520. 10.1074/jbc.M110.217315.
Deboy RT, Mongodin EF, Fouts DE, Tailford LE, Khouri H, Emerson JB, Mohamoud Y, Watkins K, Henrissat B, Gilbert HJ, et al: Insights into plant cell wall degradation from the genome sequence of the soil bacterium Cellvibrio japonicus. J Bacteriol. 2008, 190: 5455-5463. 10.1128/JB.01701-07.
Caines MEC, Vaughan MD, Tarling CA, Hancock SM, Warren RAJ, Withers SG, Strynadka NCJ: Structural and mechanistic analyses of endo-glycoceramidase II, a membrane-associated family 5 glycosidase in the Apo and G(M3) ganglioside-bound forms. J Biol Chem. 2007, 282: 14300-14308. 10.1074/jbc.M611455200.
Ishibashi Y, Nakasone T, Kiyohara M, Horibata Y, Sakaguchi K, Hijikata A, Ichinose S, Omori A, Yasui Y, Imamura A, et al: A novel endoglycoceramidase hydrolyzes oligogalactosylceramides to produce galactooligosaccharides and ceramides. J Biol Chem. 2007, 282: 11386-11396. 10.1074/jbc.M608445200.
Zverlov VV, Velikodvorskaya GA, Schwarz WH: A newly described cellulosomal cellobiohydrolase, CelO, from Clostridium thermocellum: investigation of the exo-mode of hydrolysis, and binding capacity to crystalline cellulose. Microbiology. 2002, 148: 247-255.
Iakiviak M, Mackie RI, Cann IKO: Functional analyses of multiple lichenin-degrading enzymes from the rumen Bacterium Ruminococcus albus 8. Appl Environ Microbiol. 2011, 77: 7541-7550. 10.1128/AEM.06088-11.
Cho KK, Kim SC, Woo JH, Bok JD, Choi YJ: Molecular cloning and expression of a novel family A endoglucanase gene from Fibrobacter succinogenes S85 in Escherichia coli. Enzyme Microb Technol. 2000, 27: 475-481. 10.1016/S0141-0229(00)00256-8.
Tanabe T, Morinaga K, Fukamizo T, Mitsutomi M: Novel chitosanase from Streptomyces griseus HUT 6037 with transglycosylation activity. Biosci Biotechnol Biochem. 2003, 67: 354-364. 10.1271/bbb.67.354.
Hong I-P, Jang H-K, Lee S-Y, Choi S-G: Cloning and characterization of a bifunctional cellulase-chitosanase gene from Bacillus licheniformis NBL420. J Microbiol Biotechnol. 2003, 13: 35-42.
Pedraza-Reyes M, Gutierrez-Corona F: The bifunctional enzyme chitosanase-cellulase produced by the gram-negative microorganism Myxobacter sp. AL-1 is highly similar to Bacillus subtilis endoglucanases. Arch Microbiol. 1997, 168: 321-327. 10.1007/s002030050505.
Yaoi K, Nakai T, Kameda Y, Hiyoshi A, Mitsuishi Y: Cloning and characterization of two xyloglucanases from Paenibacillus sp strain KM21. Appl Environ Microbiol. 2005, 71: 7670-7678. 10.1128/AEM.71.12.7670-7678.2005.
Foong FCF, Doi RH: Characterization and comparison of Clostridium cellulovorans endoglucanases-xylanases EngB and EngD hyperexpressed in Escherichia coli. J Bacteriol. 1992, 174: 1403-1409.
Schroder R, Wegrzyn TF, Sharma NN, Atkinson RG: LeMAN4 endo-beta-mannanase from ripe tomato fruit can act as a mannan transglycosylase or hydrolase. Planta. 2006, 224: 1091-1102. 10.1007/s00425-006-0286-0.
Dilokpimol A, Nakai H, Gotfredsen CH, Baumann MJ, Nakai N, Abou Hachem M, Svensson B: Recombinant production and characterisation of two related GH5 endo-beta-1,4-mannanases from Aspergillus nidulans FGSC A4 showing distinctly different transglycosylation capacity. Biochim Biophys Acta. 2011, 1814: 1720-1729. 10.1016/j.bbapap.2011.08.003.
Dias FMV, Vincent F, Pell G, Prates JAM, Centeno MSJ, Tailford LE, Ferreira LMA, Fontes CMGA, Davies GJ, Gilbert HJ: Insights into the molecular determinants of substrate specificity in glycoside hydrolase family 5 revealed by the crystal structure and kinetics of Cellvibrio mixtus mannosidase 5A. J Biol Chem. 2004, 279: 25517-25526. 10.1074/jbc.M401647200.
Cutfield SM, Davies GJ, Murshudov G, Anderson BF, Moody PCE, Sullivan PA, Cutfield JF: The structure of the exo-beta-(1,3)-glucanase from Candida albicans in native and bound forms: Relationship between a pocket and groove in family 5 glycosyl hydrolases. J Mol Biol. 1999, 294: 771-783. 10.1006/jmbi.1999.3287.
Duenas-Santero E, Martin-Cuadrado AB, Fontaine T, Latge JP, del Rey F, de Aldana CV: Characterization of glycoside hydrolase family 5 proteins in Schizosaccharomyces pombe. Eukaryot Cell. 2010, 9: 1650-1660. 10.1128/EC.00187-10.
Schmidt S, Rainieri S, Witte S, Matern U, Martens S: Identification of a Saccharomyces cerevisiae glucosidase that hydrolyzes flavonoid glucosides. Appl Environ Microbiol. 2011, 77: 1751-1757. 10.1128/AEM.01125-10.
Ishibashi Y, Ikeda K, Sakaguchi K, Okino N, Taguchi R, Ito M: Quality control of fungus-specific glucosylceramide in Cryptococcus neoformans by endoglycoceramidase-related protein 1 (EGCrP1). J Biol Chem. 2012, 287: 368-381. 10.1074/jbc.M111.311340.
Mazzaferro LS, Breccia JD: Functional and biotechnological insights into diglycosidases. Biocatal Biotransfor. 2011, 29: 103-112. 10.3109/10242422.2011.594882.
Mazzaferro L, Pinuel L, Minig M, Breccia JD: Extracellular monoenzyme deglycosylation system of 7-O-linked flavonoid beta-rutinosides and its disaccharide transglycosylation activity from Stilbella fimetaria. Arch Microbiol. 2010, 192: 383-393. 10.1007/s00203-010-0567-7.
Tsuruhami K, Mori S, Amarume S, Saruwatari S, Murata T, Hirakake J, Sakata K, Usui T: Isolation and characterization of a beta-primeverosidase-like enzyme from Penicillium multicolor. Biosci Biotechnol Biochem. 2006, 70: 691-698. 10.1271/bbb.70.691.
Pereira JH, Chen ZW, McAndrew RP, Sapra R, Chhabra SR, Sale KL, Simmons BA, Adams PD: Biochemical characterization and crystal structure of endoglucanase Cel5A from the hyperthermophilic Thermotoga maritima. J Struct Biol. 2010, 172: 372-379. 10.1016/j.jsb.2010.06.018.
Wu TH, Huang CH, Ko TP, Lai HL, Ma YH, Chen CC, Cheng YS, Liu JR, Guo RT: Diverse substrate recognition mechanism revealed by Thermotoga maritima Cel5A structures in complex with cellotetraose, cellobiose and mannotriose. Biochim Biophys Acta. 2011, 1814: 1832-1840. 10.1016/j.bbapap.2011.07.020.
Voget S, Steele HL, Streit WR: Characterization of a metagenome-derived halotolerant cellulase. J Biotechnol. 2006, 126: 26-36. 10.1016/j.jbiotec.2006.02.011.
Khan MAS, Akbar M, Kitaoka M, Hayashi K: A unique thermostable lichenase from Thermotoga maritima MSB8 with divergent substrate specificity. Indian J Biotechnol. 2007, 6: 315-320.
Han YJ, Dodd D, Hespen CW, Ohene-Adjei S, Schroeder CM, Mackie RI, Cann IKO: Comparative analyses of two thermophilic enzymes exhibiting both beta-1,4 mannosidic and beta-1,4 glucosidic cleavage activities from Caldanaerobius polysaccharolyticus. J Bacteriol. 2010, 192: 4111-4121. 10.1128/JB.00257-10.
Wang WY, Thomson JA: Nucleotide sequence of the celA gene encoding a cellodextrinase of Ruminococcus flavefaciens Fd-1. Mol Gen Genet. 1990, 222: 265-269. 10.1007/BF00633827.
Matsui H, Ogata K, Tajima K, Nagamine T, Nakamura M, Aminov R, Benno Y: Cloning, expression, and characterization of a cellulase gene from Prevotella ruminicola. Anim Sci J. 2001, 72: 421-426.
Singh NA, Shanmugam V: Cloning and characterization of a bifunctional glycosyl hydrolase from an antagonistic Pseudomonas putida strain P3(4). J Basic Microbiol. 2012, 52: 340-349. 10.1002/jobm.201100232.
Graham JE, Clark ME, Nadler DC, Huffer S, Chokhawala HA, Rowland SE, Blanch HW, Clark DS, Robb FT: Identification and characterization of a multidomain hyperthermophilic cellulase from an archaeal enrichment. Nat Commun. 2011, 2: 375-
Temuujin U, Kim JW, Kim JK, Lee BM, Kang HW: Identification of novel pathogenicity-related cellulase genes in Xanthomonas oryzae pv. oryzae. Physiol Mol Plant Pathol. 2011, 76: 152-157. 10.1016/j.pmpp.2011.08.004.
Bahari L, Gilad Y, Borovok I, Kahel-Raifer H, Dassa B, Nataf Y, Shoham Y, Lamed R, Bayer EA: Glycoside hydrolases as components of putative carbohydrate biosensor proteins in Clostridium thermocellum. J Ind Microbiol Biotechnol. 2011, 38: 825-832. 10.1007/s10295-010-0848-9.
Durand A, Hughes R, Roussel A, Flatman R, Henrissat B, Juge N: Emergence of a subfamily of xylanase inhibitors within glycoside hydrolase family 18. FEBS J. 2005, 272: 1745-1755. 10.1111/j.1742-4658.2005.04606.x.
Qasba PK, Kumar S: Molecular divergence of lysozymes and alpha-lactalbumin. Crit Rev Biochem Mol Biol. 1997, 32: 255-306. 10.3109/10409239709082574.
Sinnott ML: Catalytic Mechanisms of Enzymatic Glycosyl Transfer. Chem Rev. 1990, 90: 1171-1202. 10.1021/cr00105a006.
The Arabidopsis Information Resource. http://www.arabidopsis.org,
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Price MN, Dehal PS, Arkin AP: FastTree 2-approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: 9490-10.1371/journal.pone.0009490.
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
We wish to thank Franz St John for his work on family GH30, which stimulated us to refine family GH5. Work at UMR 7257 of Aix-Marseille Université was funded by ANR grant BIP:BIP (ANR-10-BINF-03-04). Work in the Division of Glycoscience, KTH Biotechnology was funded by an Ingvar Carlsson Award to H.A. from the Swedish Foundation for Strategic Research. The Swedish Research Council Formas provided additional funding via CarboMat – The KTH Advanced Carbohydrate Materials Consortium (A Formas Strong Research Environment). H.B. gratefully acknowledges faculty support from The Michael Smith Laboratories and the University of British Columbia.
The authors declare they have no competing interests.
Pedro M. Coutinho performed the global phylogenetic analysis; Henrik Aspeborg and Yang Wang performed the phylogenetic analyses of each subfamily; Bernard Henrissat and Pedro M. Coutinho collected and curated all GH5 sequences used in the present analysis; Harry Brumer and Bernard Henrissat designed research and analyzed the data. The paper was collectively written. All authors read and approved the final manuscript.
Henrik Aspeborg, Pedro M Coutinho contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1. Rectangular phylogram view of the phylogenetic tree of family GH5. Branches corresponding to subfamilies 1–53 are shown in color and the individual subfamilies have their corresponding subfamily numbers as indicated in Figure 1. The branches corresponding to sequences not included into subfamilies are in black. Each individual protein module node is identified by a varying number of fields separated by “|” indicating: (i) the organism, with 3 letters for the genre and either 5 letters for the species or full strain code; (ii) the protein accession in public databases, typically GenBank; (iii) if attributed, the subfamily number or other information; (iv) if available, EC numbers (node in bold) or a “*” (node in bold and italic) to indicate precise enzyme characterizations or a simple activity tests, respectively. A suffix like “_2” may indicate the module position if more than one GH5 module is present on peptide. Lower confidence nodes with a SH-like local support below 0.7 (varying from low 0 to strong 1) are indicated with a black dot. Identified sequences without complete catalytic machinery are in red. Individual subfamily trees are also included in this file. (PDF 5 MB)
About this article
Cite this article
Aspeborg, H., Coutinho, P.M., Wang, Y. et al. Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5). BMC Evol Biol 12, 186 (2012). https://doi.org/10.1186/1471-2148-12-186
- Protein evolution
- Enzyme evolution
- Functional prediction
- Glycoside hydrolase family 5
- Phylogenetic analysis
- Subfamily classification