Large-scale genomic analyses with machine learning uncover predictive patterns associated with fungal phytopathogenic lifestyles and traits

Invasive plant pathogenic fungi have a global impact, with devastating economic and environmental effects on crops and forests. Biosurveillance, a critical component of threat mitigation, requires risk prediction based on fungal lifestyles and traits. Recent studies have revealed distinct genomic patterns associated with specific groups of plant pathogenic fungi. We sought to establish whether these phytopathogenic genomic patterns hold across diverse taxonomic and ecological groups from the Ascomycota and Basidiomycota, and furthermore, if those patterns can be used in a predictive capacity for biosurveillance. Using a supervised machine learning approach that integrates phylogenetic and genomic data, we analyzed 387 fungal genomes to test a proof-of-concept for the use of genomic signatures in predicting fungal phytopathogenic lifestyles and traits during biosurveillance activities. Our machine learning feature sets were derived from genome annotation data of carbohydrate-active enzymes (CAZymes), peptidases, secondary metabolite clusters (SMCs), transporters, and transcription factors. We found that machine learning could successfully predict fungal lifestyles and traits across taxonomic groups, with the best predictive performance coming from feature sets comprising CAZyme, peptidase, and SMC data. While phylogeny was an important component in most predictions, the inclusion of genomic data improved prediction performance for every lifestyle and trait tested. Plant pathogenicity was one of the best-predicted traits, showing the promise of predictive genomics for biosurveillance applications. Furthermore, our machine learning approach revealed expansions in the number of genes from specific CAZyme and peptidase families in the genomes of plant pathogens compared to non-phytopathogenic genomes (saprotrophs, endo- and ectomycorrhizal fungi). Such genomic feature profiles give insight into the evolution of fungal phytopathogenicity and could be useful to predict the risks of unknown fungi in future biosurveillance activities.

101 Dothideomycetes genomes: A test case for predicting lifestyles and emergence of pathogens

Dothideomycetes is the largest class of kingdom Fungi and comprises an incredible diversity of lifestyles, many of which have evolved multiple times. Plant pathogens represent a major ecological niche of the class Dothideomycetes and they are known to infect most major food crops and feedstocks for biomass and biofuel production. Studying the ecology and evolution of Dothideomycetes has significant implications for our fundamental understanding of fungal evolution, their adaptation to stress and host specificity, and practical implications with regard to the effects of climate change and on the food, feed, and livestock elements of the agro-economy. In this study, we present the first large-scale, whole-genome comparison of 101 Dothideomycetes introducing 55 newly sequenced species. The availability of whole-genome data produced a high-confidence phylogeny leading to reclassification of 25 organisms, provided a clearer picture of the relationships among the various families, and indicated that pathogenicity evolved multiple times within this class. We also identified gene family expansions and contractions across the Dothideomycetes phylogeny linked to ecological niches providing insights into genome evolution and adaptation across this group. Using machine-learning methods we classified fungi into lifestyle classes with >95 % accuracy and identified a small number of gene families that positively correlated with these distinctions. This can become a valuable tool for genome-based prediction of species lifestyle, especially for rarely seen and poorly studied species.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota)

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Comparative genomics of Mortierella elongata and its bacterial endosymbiont Mycoavidus cysteinexigens

Endosymbiosis of bacteria by eukaryotes is a defining feature of cellular evolution. In addition to well-known bacterial origins for mitochondria and chloroplasts, multiple origins of bacterial endosymbiosis are known within the cells of diverse animals, plants and fungi. Early-diverging lineages of terrestrial fungi harbor endosymbiotic bacteria belonging to the Burkholderiaceae. We sequenced the metagenome of the soil-inhabiting fungus Mortierella elongata and assembled the complete circular chromosome of its endosymbiont, Mycoavidus cysteinexigens, which we place within a lineage of endofungal symbionts that are sister clade to Burkholderia. The genome of M. elongata strain AG77 features a core set of primary metabolic pathways for degradation of simple carbohydrates and lipid biosynthesis, while the M. cysteinexigens (AG77) genome is reduced in size and function. Experiments using antibiotics to cure the endobacterium from the host demonstrate that the fungal host metabolism is highly modulated by presence/absence of M. cysteinexigens. Independent comparative phylogenomic analyses of fungal and bacterial genomes are consistent with an ancient origin for M. elongata - M. cysteinexigens symbiosis, most likely over 350 million years ago and concomitant with the terrestrialization of Earth and diversification of land fungi and plants.

Comparative genomics analysis of Lactobacillus species associated with weight gain or weight protection

BACKGROUND:

Some Lactobacillus species are associated with obesity and weight gain while others are associated with weight loss. Lactobacillus spp. and bifidobacteria represent a major bacterial population of the small intestine where lipids and simple carbohydrates are absorbed, particularly in the duodenum and jejunum. The objective of this study was to identify Lactobacillus spp. proteins involved in carbohydrate and lipid metabolism associated with weight modifications.

METHODS:

We examined a total of 13 complete genomes belonging to seven different Lactobacillus spp. previously associated with weight gain or weight protection. We combined the data obtained from the Rapid Annotation using Subsystem Technology, Batch CD-Search and Gene Ontology to classify gene function in each genome.

RESULTS:

We observed major differences between the two groups of genomes. Weight gain-associated Lactobacillus spp. appear to lack enzymes involved in the catabolism of fructose, defense against oxidative stress and the synthesis of dextrin, L-rhamnose and acetate. Weight protection-associated Lactobacillus spp. encoded a significant gene amount of glucose permease. Regarding lipid metabolism, thiolases were only encoded in the genome of weight gain-associated Lactobacillus spp. In addition, we identified 18 different types of bacteriocins in the studied genomes, and weight gain-associated Lactobacillus spp. encoded more bacteriocins than weight protection-associated Lactobacillus spp.

CONCLUSIONS:

The results of this study revealed that weight protection-associated Lactobacillus spp. have developed defense mechanisms for enhanced glycolysis and defense against oxidative stress. Weight gain-associated Lactobacillus spp. possess a limited ability to breakdown fructose or glucose and might reduce ileal brake effects.

Carbohydrate-active enzymes revealed in Coptotermes formosanus (Isoptera: Rhinotermitidae) transcriptome

Coptotermes formosanus is one of the most destructive wood-feeding termites. To understand the molecular mechanisms that regulate the development of the termite, a normalized C. formosanus cDNA library was constructed using mixed RNA isolated from workers, soldiers, nymphs and alates of both sexes. The sequencing of this library generated 131 636 expressed sequence tags (ESTs) and 25 939 assembled unigenes. The carbohydrate-active enzymes (CAZymes) revealed in this library were analysed in the present report. A total of 509 putative CAZymes were identified. Diverse cellulolytic enzymes were uncovered from both the host termite and from symbionts harboured by the termite, which were possibly the result of the high efficiency of cellulose utilization. CAZymes associated with trehalose biosynthetic and metabolic pathways were also identified, which are potential regulators of the physiological activities of trehalose, an important insect blood sugar. Representative CAZyme coding genes in glycoside hydrolase family 1 (GH1) were quantitatively analysed. The results showed that the five GH1 β-glucosidase genes were expressed differentially among different castes and one of them was female alate-specific. Overall, the normalized EST library provides a comprehensive genetic resource of C. formosanus and will serve a diverse range of research areas. The CAZymes represent one of the repositories of enzymes useful for physiological studies and applications in sugar-based biofuel production.

Glycosyltransferases: structures, functions, and mechanisms

Glycosyltransferases catalyze glycosidic bond formation using sugar donors containing a nucleoside phosphate or a lipid phosphate leaving group. Only two structural folds, GT-A and GT-B, have been identified for the nucleotide sugar-dependent enzymes, but other folds are now appearing for the soluble domains of lipid phosphosugar-dependent glycosyl transferases. Structural and kinetic studies have provided new insights. Inverting glycosyltransferases utilize a direct displacement S(N)2-like mechanism involving an enzymatic base catalyst. Leaving group departure in GT-A fold enzymes is typically facilitated via a coordinated divalent cation, whereas GT-B fold enzymes instead use positively charged side chains and/or hydroxyls and helix dipoles. The mechanism of retaining glycosyltransferases is less clear. The expected two-step double-displacement mechanism is rendered less likely by the lack of conserved architecture in the region where a catalytic nucleophile would be expected. A mechanism involving a short-lived oxocarbenium ion intermediate now seems the most likely, with the leaving phosphate serving as the base.

The genome of black cottonwood, Populus trichocarpa (Torr. & Gray)

We report the draft genome of the black cottonwood tree, Populus trichocarpa. Integration of shotgun sequence assembly with genetic mapping enabled chromosome-scale reconstruction of the genome. More than 45,000 putative protein-coding genes were identified. Analysis of the assembled genome revealed a whole-genome duplication event; about 8000 pairs of duplicated genes from that event survived in the Populus genome. A second, older duplication event is indistinguishably coincident with the divergence of the Populus and Arabidopsis lineages. Nucleotide substitution, tandem gene duplication, and gross chromosomal rearrangement appear to proceed substantially more slowly in Populus than in Arabidopsis. Populus has more protein-coding genes than Arabidopsis, ranging on average from 1.4 to 1.6 putative Populus homologs for each Arabidopsis gene. However, the relative frequency of protein domains in the two genomes is similar. Overrepresented exceptions in Populus include genes associated with lignocellulosic wall biosynthesis, meristem development, disease resistance, and metabolite transport.

Characterization of some mesogenic alkyl 1-thioglycosides

A series of dodecyl 1-thio-β-D-glycosides has been synthesized and characterized (DSC, NMR, CP MAS, X-ray diffraction) as possible new marking materials with liquid-crystalline properties. These compounds undergo solid to liquid crystal phase transitions at various temperatures, which depend on the nature of the carbohydrate part of the structure. Their liquid-crystalline phases show extreme shear thinning behaviour.Key words: liquid crystal, powder X-ray diffraction, phase transition, thioglycoside, solid-state NMR, marking material

Structural enzymology of carbohydrate-active enzymes: implications for the post-genomic era

Simple and complex carbohydrates have been described as "the last frontier of molecular and cell biology". The enzymes that are required for the synthesis and degradation of these compounds provide an enormous challenge in the post-genomic era. This reflects both the extreme chemical and functional diversity of sugars and the difficulties in characterizing both the substrates and the enzymes themselves. The vast myriad of enzymes involved in the synthesis, modification and degradation of oligosaccharides and polysaccharides is only just being unveiled by genomic sequencing. These so-called "carbohydrate-active enzymes" lend themselves to classification by sensitive sequence similarity detection methods. The modularity, often extremely complex, of these enzymes must first be dissected and annotated before high throughput characterization or "structural genomics" approaches may be employed. Once achieved, modular analysis also permits collation of a detailed "census" of carbohydrate-active enzymes for a whole organism or throughout an ecosystem. At the structural level, improvements in X-ray crystallography have opened up a three-dimensional understanding of the way these enzymes work. The mechanisms of many of the glycoside hydrolase families are becoming clearer, yet glycosyltransferases are only slowly revealing their secrets. What is clear from the genomic and structural data is that if we are to harness the latent power of glycogenomics, scientists must consider distant sequence relatives revealed by the sequence families or other sensitive detection methods.

The location of the ligand-binding site of carbohydrate-binding modules that have evolved from a common sequence is not conserved

Polysaccharide-degrading enzymes are generally modular proteins that contain non-catalytic carbohydrate-binding modules (CBMs), which potentiate the activity of the catalytic module. CBMs have been grouped into sequence-based families, and three-dimensional structural data are available for half of these families. Clostridium thermocellum xylanase 11A is a modular enzyme that contains a CBM from family 6 (CBM6), for which no structural data are available. We have determined the crystal structure of this module to a resolution of 2.1 A. The protein is a beta-sandwich that contains two potential ligand-binding clefts designated cleft A and B. The CBM interacts primarily with xylan, and NMR spectroscopy coupled with site-directed mutagenesis identified cleft A, containing Trp-92, Tyr-34, and Asn-120, as the ligand-binding site. The overall fold of CBM6 is similar to proteins in CBM families 4 and 22, although surprisingly the ligand-binding site in CBM4 and CBM22 is equivalent to cleft B in CBM6. These structural data define a superfamily of CBMs, comprising CBM4, CBM6, and CBM22, and demonstrate that, although CBMs have evolved from a relatively small number of ancestors, the structural elements involved in ligand recognition have been assembled at different locations on the ancestral scaffold.

The iota-carrageenase of Alteromonas fortis. A beta-helix fold-containing enzyme for the degradation of a highly polyanionic polysaccharide

Carrageenans are gel-forming hydrocolloids extracted from the cell walls of marine red algae. They consist of d-galactose residues bound by alternate alpha(1-->3) and beta(1-->4) linkages and substituted by one (kappa-carrageenan), two (iota-carrageenan), or three (lambda-carrageenan) sulfate-ester groups per disaccharide repeating unit. Both the kappa- and iota-carrageenan chains adopt ordered conformations leading to the formation of highly ordered aggregates of double-stranded helices. Several kappa-carrageenases and iota-carrageenases have been cloned from marine bacteria. Kappa-carrageenases belong to family 16 of the glycoside hydrolases, which essentially encompasses polysaccharidases specialized in the hydrolysis of the neutral polysaccharides such as agarose, laminarin, lichenan, and xyloglucan. In contrast, iota-carrageenases constitute a novel glycoside hydrolase structural family. We report here the crystal structure of Alteromonas fortis iota-carrageenase at 1.6 A resolution. The enzyme folds into a right-handed parallel beta-helix of 10 complete turns with two additional C-terminal domains. Glu(245), Asp(247), or Glu(310), in the cleft of the enzyme, are proposed as candidate catalytic residues. The protein contains one sodium and one chloride binding site and three calcium binding sites shown to be involved in stabilizing the enzyme structure.

Glycoside hydrolases and glycosyltransferases: families and functional modules

The past year has witnessed the expected increase in the number of solved structures of glycoside hydrolases and glycosyltransferases, and their constitutive modules. These structures show that, while glycoside hydrolases display an extraordinary variety of folds, glycosyltransferases and carbohydrate-binding modules appear to belong to a much smaller number of folding families.

A census of carbohydrate-active enzymes in the genome of Arabidopsis thaliana

The synthesis, modification, and breakdown of carbohydrates is one of the most fundamentally important reactions in nature. The structural and functional diversity of glycosides is mirrored by a vast array of enzymes involved in their synthesis (glycosyltransferases), modification (carbohydrate esterases) and breakdown (glycoside hydrolases and polysaccharide lyases). The importance of these processes is reflected in the dedication of 1-2% of an organism's genes to glycoside hydrolases and glycosyltransferases alone. In plants, these processes are of particular importance for cell-wall synthesis and expansion. starch metabolism, defence against pathogens, symbiosis and signalling. Here we present an analysis of over 730 open reading frames representing the two main classes of carbohydrate-active enzymes, glycoside hydrolases and glycosyltransferases, in the genome of Arabidopsis thaliana. The vast importance of these enzymes in cell-wall formation and degradation is revealed along with the unexpected dominance of pectin degradation in Arabidopsis, with at least 170 open-reading frames dedicated solely to this task.

Glycosidase active site mutations in human alpha-L-iduronidase

Mucopolysaccharidosis type I (MPS I; McKusick 25280) results from a deficiency in alpha-L-iduronidase activity. Using a bioinformatics approach, we have previously predicted the putative acid/base catalyst and nucleophile residues in the active site of this human lysosomal glycosidase to be Glu182 and Glu299, respectively. To obtain experimental evidence supporting these predictions, wild-type alpha-L-iduronidase and site-directed mutants E182A and E299A were individually expressed in Chinese hamster ovary-K1 cell lines. We have compared the synthesis, processing, and catalytic properties of the two mutant proteins with wild-type human alpha-L-iduronidase. Both E182A and E299A transfected cells produced catalytically inactive human alpha-L-iduronidase protein at levels comparable to the wild-type control. The E182A protein was synthesized, processed, targeted to the lysosome, and secreted in a similar fashion to wild-type alpha-L-iduronidase. The E299A mutant protein was also synthesized and secreted similarly to the wild-type enzyme, but there were alterations in its rate of traffic and proteolytic processing. These data indicate that the enzymatic inactivity of the E182A and E299A mutants is not due to problems of synthesis/folding, but to the removal of key catalytic residues. In addition, we have identified a MPS I patient with an E182K mutant allele. The E182K mutant protein was expressed in CHO-K1 cells and also found to be enzymatically inactive. Together, these results support the predicted role of E182 and E299 in the catalytic mechanism of alpha-L-iduronidase and we propose that the mutation of either of these residues would contribute to a very severe clinical phenotype in a MPS I patient.

Enzyme-catalyzed condensation reaction in a mammalian alpha-amylase. High-resolution structural analysis of an enzyme-inhibitor complex

Mammalian alpha-amylases catalyze the hydrolysis of alpha-linked glucose polymers according to a complex processive mechanism. We have determined the X-ray structures of porcine pancreatic alpha-amylase complexes with the smallest molecule of the trestatin family (acarviosine-glucose) which inhibits porcine pancreatic alpha-amylase and yet is not hydrolyzed by the enzyme. A structure analysis at 1.38 A resolution of this complex allowed for a clear identification of a genuine single hexasaccharide species composed of two alpha-1,4-linked original molecules bound to the active site of the enzyme. The structural results supported by mass spectrometry experiments provide evidence for an enzymatically catalyzed condensation reaction in the crystal.

Streptomyces matensis laminaripentaose hydrolase is an 'inverting' beta-1,3-glucanase

The laminaripentaose-producing beta-1,3-glucanase of Streptomyces matensis is a member of the glycoside hydrolase family GH-64. We have constructed and purified a recombinant hexahistidine-tagged form of the enzyme for characterisation. The enzyme, which exists as a monomer in solution, hydrolyses beta-1,3-glucan by a mechanism leading to overall inversion of the anomeric configuration. This is the first determination of the mechanism prevailing in glycoside hydrolase family GH-64 and this is the first characterisation of an 'inverting' beta-1,3-glucanase.

The kappa-carrageenase of P. carrageenovora features a tunnel-shaped active site: a novel insight in the evolution of Clan-B glycoside hydrolases

BACKGROUND

kappa-carrageenans are gel-forming, sulfated 1,3-alpha-1,4-beta-galactans from the cell walls of marine red algae. The kappa-carrageenase from the marine, gram-negative bacterium Pseudoalteromonas carrageenovora degrades kappa-carrageenan both in solution and in solid state by an endoprocessive mechanism. This beta-galactanase belongs to the clan-B of glycoside hydrolases.

RESULTS

The structure of P. carrageenovora kappa-carrageenase has been solved to 1.54 A resolution by the multiwavelength anomalous diffraction (MAD) method, using a seleno-methionine-substituted form of the enzyme. The enzyme folds into a curved beta sandwich, with a tunnel-like active site cavity. Another remarkable characteristic is the presence of an arginine residue at subsite -1.

CONCLUSIONS

The crystal structure of P. carrageenovora kappa-carrageenase is the first three-dimensional structure of a carrageenase. Its tunnel-shaped active site, the first to be reported for enzymes other than cellulases, suggests that such tunnels are associated with the degradation of solid polysaccharides. Clan-B glycoside hydrolases fall into two subgroups, one with catalytic machinery held by an ancestral beta bulge, and the other in which it is held by a regular beta strand. At subsite -1, all of these hydrolases exhibit an aromatic amino acid that interacts with the hexopyranose ring of the monosaccharide undergoing catalysis. In addition, in kappa-carrageenases, an arginine residue recognizes the sulfate-ester substituents of the beta-linked kappa-carrageenan monomers. It also appears that, in addition to the nucleophile and acid/base catalysts, two other amino acids are involved with the catalytic cycle, accelerating the deglycosylation step.

The celA gene, encoding a glycosyl hydrolase family 3 beta-glucosidase in Azospirillum irakense, is required for optimal growth on cellobiosides

The CelA beta-glucosidase of Azospirillum irakense, belonging to glycosyl hydrolase family 3 (GHF3), preferentially hydrolyzes cellobiose and releases glucose units from the C(3), C(4), and C(5) oligosaccharides. The growth of a DeltacelA mutant on these cellobiosides was affected. In A. irakense, the GHF3 beta-glucosidases appear to be functional alternatives for the GHF1 beta-glucosidases in the assimilation of beta-glucosides by other bacteria.

Crystal structure of a monocotyledon (maize ZMGlu1) beta-glucosidase and a model of its complex with p-nitrophenyl beta-D-thioglucoside

The maize beta-glucosidase isoenzymes ZMGlu1 and ZMGlu2 hydrolyse the abundant natural substrate DIMBOAGlc (2-O-beta-D-glucopyranosyl-4-hydroxy-7-methoxy-1,4-benzoxazin-3-one), whose aglycone DIMBOA (2,4-hydroxy-7-methoxy-1,4-benzoxazin-3-one) is the major defence chemical protecting seedlings and young plant parts against herbivores and other pests. The two isoenzymes hydrolyse DIMBOAGlc with similar kinetics but differ from each other and their sorghum homologues with respect to specificity towards other substrates. To gain insights into the mechanism of substrate (i.e. aglycone) specificity between the two maize isoenzymes and their sorghum homologues, ZMGlu1 was produced in Escherichia coli, purified, crystallized and its structure solved at 2.5 Angstrom resolution by X-ray crystallography. In addition, the complex of ZMGlu1 with the non-hydrolysable inhibitor p-nitrophenyl beta-D-thioglucoside was crystallized and, based on the partial electron density, a model for the inhibitor molecule within the active site is proposed. The inhibitor is located in a slot-like active site where its aromatic aglycone is held by stacking interactions with Trp-378. Whereas some of the atoms on the non-reducing end of the glucose moiety can be modelled on the basis of the electron density, most of the inhibitor atoms are highly disordered. This is attributed to the requirement of the enzyme to accommodate two different species, namely the substrate in its ground state and in its distorted conformation, for catalysis.

Optimized mixtures of recombinant Humicola insolens cellulases for the biodegradation of crystalline cellulose

The digestion of bacterial cellulose ribbons by ternary mixtures of enzymes consisting of recombinant cellulases (two cellobiohydrolases, Cel6A and Cel7A, and the endoglucanase Cel45A) from Humicola insolens was investigated over a wide range of mixture composition. The extent of digestion was followed by soluble sugar release (saccharification) analysis together with transmission electron microscopy (TEM) observations. It was found that the addition of minute quantities of Cel45A induced a spectacular increase in saccharification of the substrate with either Cel7A or the mixture of Cel6A and Cel7A. Conversely, only a moderate saccharification resulted from the mixing of Cel45A and Cel6A. This difference is believed to originate from (1) the occasional endo character of Cel6A and (2) the competition of Cel6A and Cel45A for the substrate sites that are sensitive to endo activity. Interestingly, the mixture of enzymes giving rise to the highest saccharification rate did not always correspond to mixtures of enzymes generating the highest synergy. TEM images revealed that the bacterial cellulose ribbons became at the same time cut and narrowed down under the action of an optimized mixture of the three enzymes.

High resolution X-ray crystallography shows that ascorbate is a cofactor for myrosinase and substitutes for the function of the catalytic base

Myrosinase, an S-glycosidase, hydrolyzes plant anionic 1-thio-beta-d-glucosides (glucosinolates) considered part of the plant defense system. Although O-glycosidases are ubiquitous, myrosinase is the only known S-glycosidase. Its active site is very similar to that of retaining O-glycosidases, but one of the catalytic residues in O-glycosidases, a carboxylate residue functioning as the general base, is replaced by a glutamine residue. Myrosinase is strongly activated by ascorbic acid. Several binary and ternary complexes of myrosinase with different transition state analogues and ascorbic acid have been analyzed at high resolution by x-ray crystallography along with a 2-deoxy-2-fluoro-glucosyl enzyme intermediate. One of the inhibitors, d-gluconhydroximo-1,5-lactam, binds simultaneously with a sulfate ion to form a mimic of the enzyme-substrate complex. Ascorbate binds to a site distinct from the glucose binding site but overlapping with the aglycon binding site, suggesting that activation occurs at the second step of catalysis, i.e. hydrolysis of the glycosyl enzyme. A water molecule is placed perfectly for activation by ascorbate and for nucleophilic attack on the covalently trapped 2-fluoro-glucosyl-moiety. Activation of the hydrolysis of the glucosyl enzyme intermediate is further evidenced by the observation that ascorbate enhances the rate of reactivation of the 2-fluoro-glycosyl enzyme, leading to the conclusion that ascorbic acid substitutes for the catalytic base in myrosinase.

The mechanism of substrate (aglycone) specificity in beta -glucosidases is revealed by crystal structures of mutant maize beta -glucosidase-DIMBOA, -DIMBOAGlc, and -dhurrin complexes

The mechanism and the site of substrate (i.e., aglycone) recognition and specificity were investigated in maize beta-glucosidase (Glu1) by x-ray crystallography by using crystals of a catalytically inactive mutant (Glu1E191D) in complex with the natural substrate 2-O-beta-d-glucopyranosyl-4-hydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOAGlc), the free aglycone DIMBOA, and competitive inhibitor para-hydroxy-S-mandelonitrile beta-glucoside (dhurrin). The structures of these complexes and of the free enzyme were solved at 2.1-, 2.1-, 2.0-, and 2.2-A resolution, respectively. The structural data from the complexes allowed us to visualize an intact substrate, free aglycone, or a competitive inhibitor in the slot-like active site of a beta-glucosidase. These data show that the aglycone moiety of the substrate is sandwiched between W378 on one side and F198, F205, and F466 on the other. Thus, specific conformations of these four hydrophobic amino acids and the shape of the aglycone-binding site they form determine aglycone recognition and substrate specificity in Glu1. In addition to these four residues, A467 interacts with the 7-methoxy group of DIMBOA. All residues but W378 are variable among beta-glucosidases that differ in substrate specificity, supporting the conclusion that these sites are the basis of aglycone recognition and binding (i.e., substrate specificity) in beta-glucosidases. The data also provide a plausible explanation for the competitive binding of dhurrin to maize beta-glucosidases with high affinity without being hydrolyzed.

Solution structure of the module X2 1 of unknown function of the cellulosomal scaffolding protein CipC of Clostridium cellulolyticum

Multidimensional, homo- and heteronuclear magnetic resonance spectroscopy combined with dynamical annealing has been used to determine the structure of a 94 residue module (X2 1) of the scaffolding protein CipC from the anaerobic bacterium Clostridium cellulolyticum. An experimental data set comprising 1647 nuclear Overhauser effect-derived restraints, 105 hydrogen bond restraints and 66 phi torsion angle restraints was used to calculate 20 converging final solutions. The calculated structures have an average rmsd about the mean structure of 0.55(+/-0.11) A for backbone atoms and 1.40(+/-0.11) A for all heavy atoms when fitted over the secondary structural elements. The X2 1 module has an immunoglobulin-like fold with two beta-sheets packed against each other. One sheet contains three strands, the second contains four strands. An additional strand is intercalated between the beta-sandwich, as well as two turns of a 3(.10) helix. X2 1 has a surprising conformational stability and may act as a conformational linker and solubility enhancer within the scaffolding protein. The fold of X2 1 is very similar to that of telokin, titin Ig domain, hemolin D2 domain, twitchin immunoglobulin domain and the first four domains of the IgSF portion of transmembrane cell adhesion molecule. As a consequence, the X2 1 module is the first prokaryotic member assigned to the I set of the immunoglobulin superfamily even though no sequence similarity with any member of this superfamily could be detected.

Crystal structure of a cohesin module from Clostridium cellulolyticum: implications for dockerin recognition

In the assembly of the Clostridium cellulolyticum cellulosome, the multiple cohesin modules of the scaffolding protein CipC serve as receptors for cellulolytic enzymes which bear a dockerin module. The X-ray structure of a type I C. cellulolyticum cohesin module (Cc-cohesin) has been solved using molecular replacement, and refined at 2.0 A resolution. Despite a rather low sequence identity of 32 %, this module has a fold close to those of the two Clostridium thermocellum cohesin (Ct-cohesin) modules whose 3D structures have been determined previously. Cc-cohesin forms a dimer in the crystal, as do the two Ct-cohesins. We show here that the dimer exists in solution and that addition of dockerin-containing proteins dissociates the dimer. This suggests that the dimerization interface and the cohesin/dockerin interface may overlap. The nature of the overall surface and of the dimer interface of Cc-cohesin differ notably from those of the Ct-cohesin modules, being much less polar, and this may explain the species specificity observed in the cohesin/dockerin interaction of C. cellulolyticum and C. thermocellum. We have produced a topology model of a C. cellulolyticum dockerin and of a Cc-cohesin/dockerin complex using homology modeling and available biochemical data. Our model suggests that a special residue pair, already identified in dockerin sequences, is located at the center of the cohesin surface putatively interacting with the dockerin.

iota-Carrageenases constitute a novel family of glycoside hydrolases, unrelated to that of kappa-carrageenases

iota-Carrageenases are polysaccharide hydrolases that cleave the beta-1,4 linkages between the d-galactose-4-sulfate and 3, 6-anhydro-d-galactose-2-sulfate residues in the red algal galactans known as iota-carrageenans. We report here on the purification of iota-carrageenase activity from the marine bacterium Zobellia galactanovorans and on the characterization of iota-carrageenase structural genes. Genomic libraries from this latter bacterium as well as from Alteromonas fortis were functionally screened for the presence of iota-carrageenase(+) clones. The Z. galactanovorans and A. fortis iota-carrageenase genes encode homologous proteins of 53.4 and 54.8 kDa, respectively. Based on hydrophobic cluster analysis and on the (1)H NMR monitoring of the products of the overexpressed A. fortis iota-carrageenase, these enzymes appear to form a new family of glycoside hydrolases, unrelated to that of kappa-carrageenases and with an inverting mechanism of hydrolysis. They both feature a 45-amino acid-long N-terminal segment with sequence similarity to the N-terminal region of several other polysaccharidases. In those for which a three-dimensional structure is available, this conspicuous segment, also deemed "glycanase motif" (Chua, J. E. H., Manning, P. A., and Morona, R. (1999) Microbiology (Reading) 145, 1649-1659), corresponds to a strand-helix-strand "cap" that covers the N-terminal end of a common, right-handed beta-helical fold.

Human glucocerebrosidase: heterologous expression of active site mutants in murine null cells

Using bioinformatics methods, we have previously identified Glu235 and Glu340 as the putative acid/base catalyst and nucleophile, respectively, in the active site of human glucocerebrosidase. Thus, we undertook site-directed mutagenesis studies to obtain experimental evidence supporting these predictions. Recombinant retroviruses were used to express wild-type and E235A and E340A mutant proteins in glucocerebrosidase-deficient murine cells. In contrast to wild-type enzyme, the mutants were found to be catalytically inactive. We also report the results of various studies (Western blotting, glycosylation analysis, subcellular fractionation, and confocal microscopy) indicating that the wild-type and mutant enzymes are identically processed and sorted to the lysosomes. Thus, enzymatic inactivity of the mutant proteins is not the result of incorrect folding/processing. These findings indicate that Glu235 plays a key role in the catalytic machinery of human glucocerebrosidase and may indeed be the acid/base catalyst. As concerns Glu340, the results both support our computer-based predictions and confirm, at the biological level, previous identification of Glu340 as the nucleophile by use of active site labeling techniques. Finally, our findings may help to better understand the molecular basis of Gaucher disease, the human lysosomal disease resulting from deficiency in glucocerebrosidase.

Classification of chitinases modules

Chitinases frequently display a modular structure featuring a catalytic domain attached to one or several ancillary noncatalytic domains whose function is often chitin binding. Gene cloning and DNA sequencing have allowed the determination of a massive number of amino acid sequences of chitinases during the last 10 years. This chapter presents a unifying classification system of the various chitinase modules that combines specific features of their sequences, three-dimensional structures and reaction mechanisms.

Identification of the catalytic residues of the first family of beta(1-3)glucanosyltransferases identified in fungi

A new family of glycosylphosphatidylinositol-anchored beta(1-3)glucanosyltransferases (Gelp), recently identified and characterized in the filamentous fungus Aspergillus fumigatus, showed functional similarity to the Gas/Phr/Epd protein families, which are involved in yeast morphogenesis. Sequence comparisons and hydrophobic cluster analysis (HCA) showed that all the Gas/Phr/Epd/Gel proteins belong to a new family of glycosylhydrolases, family 72. We confirmed by site-directed mutagenesis and biochemical analysis that the two conserved glutamate residues (the putative catalytic residues of this family, as determined by HCA) are involved in the active site of this family of glycosylhydrolases.

Structural features of normal and mutant human lysosomal glycoside hydrolases deduced from bioinformatics analysis

Lysosomal storage diseases are due to inherited deficiencies in various enzymes involved in basic metabolic processes. As with other genetic diseases, accurate structure data for these enzymatic proteins should help in better understanding the molecular effects of mutations identified in patients with the corresponding lysosomal diseases; however, no such three-dimensional (3D) structure data are available for many lysosomal enzymes. Thus, we herein intend to illustrate for an audience of molecular geneticists how structure information can nonetheless be obtained via a bioinformatics approach in the case of five human lysosomal glycoside hydrolases. Indeed, using the two-dimensional hydrophobic cluster analysis method to decipher the sequence information available in data banks for the large group of glycoside hydrolases (clan GH-A) to which these human lysosomal enzymes belong, we could deduce structure predictions for their catalytic domains and propose explanations for the molecular effects of mutations described in patients. In addition, in the case of human beta-glucuronidase for which experimental 3D data have been reported, we also show here that bioinformatics methods relying on the available 3D structure information can be used to obtain further insights into the effects of various mutations described in patients with Sly disease. In a broader perspective, our work stresses that, in the context of a rapid increase in protein sequence information through genome sequencing, bioinformatics approaches might be highly useful for generating structure-function predictions based on sequence-structure interrelationships.

Imaging the enzymatic digestion of bacterial cellulose ribbons reveals the endo character of the cellobiohydrolase Cel6A from Humicola insolens and its mode of synergy with cellobiohydrolase Cel7A

Dispersed cellulose ribbons from bacterial cellulose were subjected to digestion with cloned Cel7A (cellobiohydrolase [CBH] I) and Cel6A (CBH II) from Humicola insolens either alone or in a mixture and in the presence of an excess of beta-glucosidase. Both Cel7A and Cel6A were effective in partially converting the ribbons into soluble sugars, Cel7A being more active than Cel6A. In combination, these enzymes showed substantial synergy culminating with a molar ratio of approximately two-thirds Cel6A and one-third Cel7A. Ultrastructural transmission electron microscopy (TEM) observations indicated that Cel7A induced a thinning of the cellulose ribbons, whereas Cel6A cut the ribbons into shorter elements, indicating an endo type of action. These observations, together with the examination of the digestion kinetics, indicate that Cel6A can be classified as an endo-processive enzyme, whereas Cel7A is essentially a processive enzyme. Thus, the synergy resulting from the mixing of Cel6A and Cel7A can be explained by the partial endo character of Cel6A. A preparation of bacterial cellulose ribbons appears to be an appropriate substrate, superior to Valonia or bacterial cellulose microcrystals, to visualize directly by TEM the endo-processivity of an enzyme such as Cel6A.

Crystal structures of human pancreatic alpha-amylase in complex with carbohydrate and proteinaceous inhibitors

Crystal structures of human pancreatic alpha-amylase (HPA) in complex with naturally occurring inhibitors have been solved. The tetrasaccharide acarbose and a pseudo-pentasaccharide of the trestatin family produced identical continuous electron densities corresponding to a pentasaccharide species, spanning the -3 to +2 subsites of the enzyme, presumably resulting from transglycosylation. Binding of the acarviosine core linked to a glucose residue at subsites -1 to +2 appears to be a critical part of the interaction process between alpha-amylases and trestatin-derived inhibitors. Two crystal forms, obtained at different values of pH, for the complex of HPA with the protein inhibitor from Phaseolus vulgaris (alpha-amylase inhibitor) have been solved. The flexible loop typical of the mammalian alpha-amylases was shown to exist in two different conformations, suggesting that loop closure is pH-sensitive. Structural information is provided for the important inhibitor residue, Arg-74, which has not been observed previously in structural analyses.

Biochemical analysis of recombinant fungal mutanases. A new family of alpha1,3-glucanases with novel carbohydrate-binding domains

Nucleotide sequence analysis shows that Trichoderma harzianum and Penicillium purpurogenum alpha1,3-glucanases (mutanases) have homologous primary structures (53% amino acid sequence identity), and are composed of two distinct domains: a NH(2)-terminal catalytic domain and a putative COOH-terminal polysaccharide-binding domain separated by a O-glycosylated Pro-Ser-Thr-rich linker peptide. Each mutanase was expressed in Aspergillus oryzae host under the transcriptional control of a strong alpha-amylase gene promoter. The purified recombinant mutanases show a pH optimum in the range from pH 3.5 to 4.5 and a temperature optimum around 50-55 degrees C at pH 5.5. Also, they exhibit strong binding to insoluble mutan with K(D) around 0.11 and 0.13 microM at pH 7 for the P. purpurogenum and T. harzianum mutanases, respectively. Partial hydrolysis showed that the COOH-terminal domain of the T. harzianum mutanase binds to mutan. The catalytic domains and the binding domains were assigned to a new family of glycoside hydrolases and to a new family of carbohydrate-binding domains, respectively.

Cellulosome-like sequences in Archaeoglobus fulgidus: an enigmatic vestige of cohesin and dockerin domains

The distribution of cellulosomal cohesin domains among the sequences currently compiled in various sequence databases was investigated. Two cohesin domains were detected in two consecutive open reading frames (ORFs) of the recently sequenced genome of the archaeon Archaeoglobus fulgidus. Otherwise, no cohesin-like sequence could be detected in organisms other than those of the Eubacteria. One of the A. fulgidus cohesin-containing ORFs also harbored a dockerin domain, but the additional modular portions of both genes are undefined, both with respect to sequence homology and function. It is currently unclear what function(s) the putative cohesin and dockerin-containing proteins play in the life cycle of this organism. In particular, since A. fulgidus contains no known glycosyl hydrolase gene, the presence of a cellulosome can be excluded. The results suggest that cohesin and dockerin signature sequences cannot be used alone for the definitive identification of cellulosomes in genomes.

Localization of synthesis of beta1,6-glucan in Saccharomyces cerevisiae

Beta1,6-Glucan is a key component of the yeast cell wall, interconnecting cell wall proteins, beta1,3-glucan, and chitin. It has been postulated that the synthesis of beta1,6-glucan begins in the endoplasmic reticulum with the formation of protein-bound primer structures and that these primer structures are extended in the Golgi complex by two putative glucosyltransferases that are functionally redundant, Kre6 and Skn1. This is followed by maturation steps at the cell surface and by coupling to other cell wall macromolecules. We have reinvestigated the role of Kre6 and Skn1 in the biogenesis of beta1,6-glucan. Using hydrophobic cluster analysis, we found that Kre6 and Skn1 show significant similarities to family 16 glycoside hydrolases but not to nucleotide diphospho-sugar glycosyltransferases, indicating that they are glucosyl hydrolases or transglucosylases instead of genuine glucosyltransferases. Next, using immunogold labeling, we tried to visualize intracellular beta1,6-glucan in cryofixed sec1-1 cells which had accumulated secretory vesicles at the restrictive temperature. No intracellular labeling was observed, but the cell surface was heavily labeled. Consistent with this, we could detect substantial amounts of beta1,6-glucan in isolated plasma membrane-derived microsomes but not in post-Golgi secretory vesicles. Taken together, our data indicate that the synthesis of beta1, 6-glucan takes place largely at the cell surface. An alternative function for Kre6 and Skn1 is discussed.

The third chitinase gene (chiC) of Serratia marcescens 2170 and the relationship of its product to other bacterial chitinases

The third chitinase gene (chiC) of Serratia marcescens 2170, specifying chitinases C1 and C2, was identified. Chitinase C1 lacks a signal sequence and consists of a catalytic domain belonging to glycoside hydrolase family 18, a fibronectin type III-like domain (Fn3 domain) and a C-terminal chitin-binding domain (ChBD). Chitinase C2 corresponds to the catalytic domain of C1 and is probably generated by proteolytic removal of the Fn3 and ChBDs. The loss of the C-terminal portion reduced the hydrolytic activity towards powdered chitin and regenerated chitin, but not towards colloidal chitin and glycol chitin, illustrating the importance of the ChBD for the efficient hydrolysis of crystalline chitin. Phylogenetic analysis showed that bacterial family 18 chitinases can be clustered in three subfamilies which have diverged at an early stage of bacterial chitinase evolution. Ser. marcescens chitinase C1 is found in one subfamily, whereas chitinases A and B of the same bacterium belong to another subfamily. Chitinase C1 is the only Ser. marcescens chitinase that has an Fn3 domain. The presence of multiple, divergent, chitinases in a single chitinolytic bacterium is perhaps necessary for efficient synergistic degradation of chitin.

Digestion of crystalline cellulose substrates by the clostridium thermocellum cellulosome: structural and morphological aspects

The action of cellulosomes from Clostridium thermocellum on model cellulose microfibrils from Acetobacter xylinum and cellulose microcrystals from Valonia ventricosa was investigated. The biodegradation of these substrates was followed by transmission electron microscopy, Fourier-transform IR spectroscopy and X-ray diffraction analysis, as a function of the extent of degradation. The cellulosomes were very effective in catalysing the complete digestion of bacterial cellulose, but the total degradation of Valonia microcrystals was achieved more slowly. Ultrastructural observations during the digestion process suggested that the rapid degradation of bacterial cellulose was the result of a very efficient synergistic action of the various enzymic components that are attached to the scaffolding protein of the cellulosomes. The degraded Valonia sample assumed various shapes, ranging from thinned-down microcrystals to crystals where one end was pointed and the other intact. This complexity may be correlated with the multi-enzyme content of the cellulosomes and possibly to a diversity of the cellulosome composition within a given batch. Another aspect of the digestion of model celluloses by cellulosomes is the relative invariability of their crystallinity, together with their Ialpha/Ibeta composition throughout the degradation process. Comparison of the action of cellulosomes with that of fungal enzymes indicated that the degradation of cellulose crystals by cellulosomes occurred with only limited levels of processivity, in contrast with the observations reported for fungal enzymes. The findings were consistent with a mechanism whereby initial attack by a cellulosome of an individual cellulose crystal results in its 'commitment' towards complete degradation.

Growth of Azospirillum irakense KBC1 on the aryl beta-glucoside salicin requires either salA or salB

The rhizosphere nitrogen-fixing bacterium Azospirillum irakense KBC1 is able to grow on pectin and beta-glucosides such as cellobiose, arbutin, and salicin. Two adjacent genes, salA and salB, conferring beta-glucosidase activity to Escherichia coli, have been identified in a cosmid library of A. irakense DNA. The SalA and SalB enzymes preferentially hydrolyzed aryl beta-glucosides. A Delta(salA-salB) A. irakense mutant was not able to grow on salicin but could still utilize arbutin, cellobiose, and glucose for growth. This mutant could be complemented by either salA or salB, suggesting functional redundancy of these genes in salicin utilization. In contrast to this functional homology, the SalA and SalB proteins, members of family 3 of the glycosyl hydrolases, show a low degree of amino acid similarity. Unlike SalA, the SalB protein exhibits an atypical truncated C-terminal region. We propose that SalA and SalB are representatives of the AB and AB' subfamilies, respectively, in glycosyl hydrolase family 3. This is the first genetic implication of this beta-glucosidase family in the utilization of beta-glucosides for microbial growth.

The gene glvA of Bacillus subtilis 168 encodes a metal-requiring, NAD(H)-dependent 6-phospho-alpha-glucosidase. Assignment to family 4 of the glycosylhydrolase superfamily

The gene glvA (formerly glv-1) from Bacillus subtilis has been cloned and expressed in Escherichia coli. The purified protein GlvA (449 residues, Mr = 50,513) is a unique 6-phosphoryl-O-alpha-D-glucopyranosyl:phosphoglucohydrolase (6-phospho-alpha-glucosidase) that requires both NAD(H) and divalent metal (Mn2+, Fe2+, Co2+, or Ni2+) for activity. 6-Phospho-alpha-glucosidase (EC 3.2.1.122) from B. subtilis cross-reacts with polyclonal antibody to maltose 6-phosphate hydrolase from Fusobacterium mortiferum, and the two proteins exhibit amino acid sequence identity of 73%. Estimates for the Mr of GlvA determined by SDS-polyacrylamide gel electrophoresis (51,000) and electrospray-mass spectroscopy (50,510) were in excellent agreement with the molecular weight of 50,513 deduced from the amino acid sequence. The sequence of the first 37 residues from the N terminus determined by automated analysis agreed precisely with that predicted by translation of glvA. The chromogenic and fluorogenic substrates, p-nitrophenyl-alpha-D-glucopyranoside 6-phosphate and 4-methylumbelliferyl-alpha-D-glucopyranoside 6-phosphate were used for the discontinuous assay and in situ detection of enzyme activity, respectively. Site-directed mutagenesis shows that three acidic residues, Asp41, Glu111, and Glu359, are required for GlvA activity. Asp41 is located at the C terminus of a betaalphabeta fold that may constitute the dinucleotide binding domain of the protein. Glu111 and Glu359 may function as the catalytic acid (proton donor) and nucleophile (base), respectively, during hydrolysis of 6-phospho-alpha-glucoside substrates including maltose 6-phosphate and trehalose 6-phosphate. In metal-free buffer, GlvA exists as an inactive dimer, but in the presence of Mn2+ ion, these species associate to form the NAD(H)-dependent catalytically active tetramer. By comparative sequence alignment with its homologs, the novel 6-phospho-alpha-glucosidase from B. subtilis can be assigned to the nine-member family 4 of the glycosylhydrolase superfamily.

The kappa-carrageenase of the marine bacterium Cytophaga drobachiensis. Structural and phylogenetic relationships within family-16 glycoside hydrolases

We report here cloning from the marine gliding bacterium Cytophaga drobachiensis of kappa-carrageenase, a glycoside hydrolase involved in the degradation of kappa-carrageenan. Structural features in the nucleotide sequence are pointed out, including the presence of an octameric omega sequence similar to the ribosome-binding sites of various eukaryotes and prokaryotes. The cgkA gene codes for a protein of 545 aa, with a signal peptide of 35 aa and a 229-aa-long posttranslationaly processed C-terminal domain. The enzyme displays the overall folding and catalytic domain characteristics of family 16 of glycoside hydrolases, which comprises other beta-1,4-alpha-1,3-D/L-galactan hydrolases, beta-1,3-D-glucan hydrolases (laminarinases), beta-1,4-1,3-D-glucan hydrolases (lichenases), and beta-1,4-D-xyloglucan endotransglycosylases. In order to address the origin and evolution of CgkA, a comprehensive phylogenetic tree of family 16 was built using parsimony analysis. Family-16 glycoside hydrolases cluster according to their substrate specificity, regardless of their phylogenetic distribution over eubacteria and eukaryotes. Such a topology suggests that the general homology between laminarinases, agarases, kappa-carrageenases, lichenases, and xyloglucan endotransglycosylases has arisen through gene duplication, likely from an ancestral protein with laminarinase activity.