- Research article
- Open Access
Phylogenomic analysis of glycogen branching and debranching enzymatic duo
BMC Evolutionary Biology volume 14, Article number: 183 (2014)
Branched polymers of glucose are universally used for energy storage in cells, taking the form of glycogen in animals, fungi, Bacteria, and Archaea, and of amylopectin in plants. Some enzymes involved in glycogen and amylopectin metabolism are similarly conserved in all forms of life, but some, interestingly, are not. In this paper we focus on the phylogeny of glycogen branching and debranching enzymes, respectively involved in introducing and removing of the α(1–6) bonds in glucose polymers, bonds that provide the unique branching structure to glucose polymers.
We performed a large-scale phylogenomic analysis of branching and debranching enzymes in over 400 completely sequenced genomes, including more than 200 from eukaryotes. We show that branching and debranching enzymes can be found in all kingdoms of life, including all major groups of eukaryotes, and thus were likely to have been present in the last universal common ancestor (LUCA) but have been lost in seemingly random fashion in numerous single-celled eukaryotes. We also show how animal branching and debranching enzymes evolved from their LUCA ancestors by acquiring additional domains. Furthermore, we show that enzymes commonly perceived as orthologous, such as human branching enzyme GBE1 and E. coli branching enzyme GlgB, are in fact related by a gene duplication and consequently paralogous.
Despite being usually associated with animal liver glycogen and plant starch, energy storage in the form of branched glucose polymers is clearly an ancient process and has probably been present in the last universal common ancestor of all present life. The evolution of the enzymes enabling this form of energy storage is more complex than previously thought and illustrates the need for explicit phylogenomic analysis in the study of even seemingly “simple” metabolic enzymes. Patterns of conservation in the evolution of the glycogen/starch branching and debranching enzymes hint at some as yet unknown mechanisms, as mutations disrupting these patterns lead to a variety of genetic diseases in humans and other mammals.
In animals, glucose is stored as glycogen, whereas plants store glucose as starch. Starch is a mixture of α-amylose, a linear polysaccharide made of α(1–4) linked glucose molecules and amylopectin, a branched polysaccharide that varies from α-amylose by the presence of α(1–6) linked branches every 24 to 30 residues. Glycogen differs from amylopectin in that its α(1–6) branches occur more frequently, typically every 8 to 14 residues . In animals, glycogen forms 100 to 400 Å diameter cytoplasmic granules, which in mammals are especially noticeable in cells that have the greatest need of glycogen—liver and muscle cells, but it is also produced in other types of cells, including neurons where it can have deleterious effects . The branching is important for fast response to metabolic needs, because synthesis and degradation of the glycogen polymer can only occur from the non-reducing ends of the α-1,4 chains; therefore, highly branched glycogen has a higher number of “ends” per volume. Additionally, branching increases the water solubility of glycogen –. While the glycogen role in mammals is best known, it has also been shown to be used as a metabolic reserve in yeast and various bacteria .
Glycogen synthesis and breakdown involves a number of enzymes, such as glycogen synthase (EC 184.108.40.206), which adds glucose to the growing glycogen chain, and glycogen phosphorylase (EC 220.127.116.11), which cleaves linear α(1–4) linked glycogen chains to produce monomers of glucose-1-phosphate. The activity of glycogen phosphorylase, however, comes to a halt when it approaches an α(1–6) linked branch point four units away. In this situation, the action of a debranching enzyme, which removes α(1–6) linkages, becomes necessary for continued glycogen breakdown ,. Such debranching enzymes—together with their enzymatic opposites, branching enzymes that introduce α(1–6) linkages—are the focus of this work. Humans, and other mammals, possess one branching enzyme and one debranching enzyme, occurring in various isoforms  (see Figure 1 for an overview).
The human glycogen branching enzyme (EC 18.104.22.168), also referred to as Amylo-(1,4–1,6)-transglycosylase, is encoded by the gene GBE1. This enzyme is involved in glycogen synthesis by transferring α(1–4) linked glucosyl blocks from the outer end of a growing glycogen chain to an α(1–6) position on the same or on an adjacent chain . Also, this enzyme (together with other glycogen and starch branching enzymes) has been characterized in the CAZy database  as a member of the Glycoside Hydrolase Family 57 ,. The human glycogen branching enzyme is a large, multidomain enzyme composed of three domains. The N-terminal domain of this enzyme is classified in Pfam  and CAZy as Carbohydrate-Binding Module 48 (also called Isoamylase N-terminal domain; and abbreviated as CBM_48) . The central domain is a TIM barrel glycosyl hydrolase superfamily member (Pfam: Alpha-amylase) ; and localized at the C-terminus is an all-beta domain (Pfam: Alpha-amylase_C). The N- and C-terminal domains, CBM_48 and Alpha-amylase_C, respectively, are distantly homologous and structurally similar. Both are classified as members of the Glycosyl hydrolase domain (GHD) superfamily (Pfam clan CL0369), which contains substrate binding domains of many carbohydrate hydrolases. Branching enzymes with this domain architecture are well conserved throughout all kingdoms of life, with homologs possessing all three domains having been found in plants (as starch branching enzymes) , yeast , and various Bacteria, including E. coli (gene name glgB) . The three dimensional structures of human (PDB identifier: 4BZY) as well as a variety of bacterial (for example, E. coli: 1M7X , Mycobacterium tuberculosis: 3K1D ) and archaeal (Thermococcus kodakarensis: 3N8T, 3 N92, 3 N98 ) glycogen branching enzymes, together with plant starch branching enzymes (rice Oryza sativa: 3AMK ), were determined experimentally, allowing for precise domain boundary definitions and detailed comparisons. For the human glycogen branching enzyme, these data are shown in Figure 2 which depicts the three-dimensional structure of the GBE1 gene product, in concert with domain boundaries as defined by Pfam HMMs and by the three-dimensional structure itself.
The human glycogen debranching enzyme, is encoded by the gene GDE (also called AGL). This enzyme, similar to its homologs from some other species (such as other mammals, yeast, and TreX from Sulfolobus acidocaldarium,), has two biochemical functions—that of amylo-alpha-1,6-glucosidase (EC 22.214.171.124) and of 4-alpha-glucanotransferase (EC 126.96.36.199) —. 4-alpha-glucanotransferase transfers a segment of three glucose units from α(1—6) branched four-unit chains (the result of glycogen phosphorylase activity) to an adjacent branch of the glycogen chain. Amylo-alpha-1,6-glucosidase then cleaves the α(1—6) linkage to release the remaining glucose . In other species, such as plants and E. coli (GlgX ), the glucosidase and glucanotransferase activities are carried out by two distinct enzymes (despite the high structural similarity between the E. coli glucosidase GlgX and Sulfolobus acidocaldarium TreX ), in which case only the glucosidase is referred to as a glycogen debranching enzyme —.
The human glycogen debranching enzyme is almost twice as large as the GBE1 enzyme and is composed of at least four domains (using Pfam classification): hGDE_N—hGDE_amylase—hGDE_central—GDE_C. The hGDE_amylase domain and the central domain of the branching enzyme, Alpha-amylase, are distantly related, as both are members of the TIM barrel fold containing the glycosyl hydrolase superfamily (Pfam clan CL0058) ,. On the other hand, the N_terminal DGE_C domain is predicted to have an alpha/alpha toroidal structure, consisting of several alpha hairpins arranged in a closed circular array, similar to bacterial glucoamylases. As of this writing, there are no experimentally determined three-dimensional structures of human debranching enzymes or any of its close homologs, albeit, as shown in Figure 3, reliable predictions can be made for all its domains, except the hGDE_central domain.
In contrast to the universally conserved branching enzyme, some bacterial—for instance E. coli (GlgX) , and plant debranching enzymes —are not homologous to the human debranching enzyme. In fact, they are related to the branching enzymes containing an N-terminal CBM_48 domain followed by Alpha-amylase domains. On the other hand, many Bacteria and Archaea do not have E. coli–type debranching enzymes; instead, they have a homolog of the human debranching enzyme, consisting of the prokaryote-specific C-terminal GDE domain preceded by a GDE_N domain. No direct experimental evidence for the function of the bacterial proteins with the GDE domain currently exists; however, because of their distant homology to eukaryotic debranching enzymes and their genomic distribution, where they are often found in species that lack the E. coli–type of a debranching enzyme, it is often assumed that they indeed function as debranching enzymes. As mentioned above, another difference between human and plant/E. coli debranching enzymes is that human glycogen debranching enzyme possesses a second enzymatic activity, that of a 4-alpha-glucanotransferase –.
Glycogen storage diseases
Enzymes involved in glycogen metabolism and its regulation are of great medical interest because mutations in these enzymes have been shown to lead to a wide variety of genetic diseases, collectively called glycogen storage diseases (GSDs). At least ten different types (identified with numerical sub-type designations) of GSDs have been described, with the phenotypes depending on the enzyme affected and the specific positions of the mutations within a given enzyme. Some of the mutations result in the abnormal accumulation of glycogen and/or abnormal glycogen structure (or both). For example, Cori’s disease (Type III GSD), a rare autosomal recessive genetic disorder, is caused by mutations resulting in deficiencies in the glycogen debranching enzyme, preventing depolymerization of glycogen at the α-1,6 branching points and result in the accumulation, in the liver and muscle, of abnormal glycogen with very short outer chains that cannot be broken down further. The symptoms, such as enlarged liver and hypoglycemia, are similar, but tend to be less severe, than those of Type I GSD. Interestingly, the liver symptoms usually disappear after puberty ,. One of the most severe glycogen storage diseases is Anderson’s disease (Type IV GSD) –, a very rare autosomal recessive genetic disorder caused by a defective glycogen branching enzyme (EC 188.8.131.52), leading to the formation and accumulation of abnormal glycogen with long, unbranched chains. GSDIV is also called amylopectinosis, since the glycogen in the affected cells resembles plant amylopectin. The abnormal, elongated glycogen particles lead to cell degeneration and eventually death. The exact mechanism of the cell death in GSDIV is still unknown. The phenotype of this disease is variable, involving the liver, skeletal muscle, heart, and central nervous system, alone or in combinations. The typical GSD Type IV disease presents in the first 18 months of life with an enlarged liver and cirrhosis, leading to liver failure and death by 5 years of age . Another manifestation of the same disease, in the late-onset variant, is as a neurogenerative disease called adult polyglucosan body disease (APBD) . GSDIV is also found in horses and cats ,. Table 1 lists some of the point mutations causing GSDIII, GSDIV and APBD.
Although phylogenetic studies have been performed on subsets of branching and debranching enzymes from a limited range of species (e.g., ,,), and experimental work (biochemistry, gene expression, protein structures, and biotechnological applications ) has been done on individual branching and debranching enzymes, especially those from plants, Bacteria, and Archaea, not much is known about the deep evolutionary histories of these enzyme families. The objective of the analysis presented here is to elucidate the evolution of glycogen and starch branching and debranching enzymes from a wide range of species, covering Bacteria, Archaea, and all major groups of eukaryotes. In particular, due to their involvement in glycogen storage diseases, we are interested in the evolutionary relationships of the human glycogen branching and debranching enzymes to their well-studied bacterial counterparts GlgB and GlgX/TreX. Due to the availability of more than 200 completely sequenced eukaryotic genomes, we were able to perform a large scale, protein domain-centric, comparative genomics analysis to assess the linage specific distributions, domain compositions and patterns of sequence conservation of these two important enzymes.
Results and discussion
We extracted protein sequences of glycogen branching and debranching enzyme homologs with the characteristic combinations of CBM48 and GDE_C Pfam domains (using a per domain cutoff E-value of 10–3) from 276 completely sequenced eukaryotic genomes, covering most major eukaryotic groups, as well as from select archaeal and bacterial genomes (listed in Additional file 1). Proteins with these domains were then analyzed for their overall domain architectures and for their phylogenetic relationships (listed in Additional files 2 and 3).
CBM_48—Alpha–amylase containing branching and debranching enzymes
Phylogenetic analysis of enzymes with a CBM_48—Alpha-amylase architecture (see Figure 4) shows that these enzymes can be divided into two well separated groups (100% support based on Bayesian, ML, and distance based methods): branching enzymes with a CBM_48–Alpha—amylase–Alpha-amylase_C architecture and debranching enzymes with mostly a CBM_48—Alpha-amylase architecture with either a very divergent form of the Alpha-amylase_C domain, or, in some cases, containing additional domains at the N- (such as the Bacterial pullanase-associated domain, PUD ) and C-termini (such as DUF 3372).
Proteins from the first group, with the well-defined CBM_48—Alpha-amylase—Alpha-amylase_C architecture are present in species from all kingdoms of life and perform branching functions in glycogen and starch biosynthesis pathways (EC 184.108.40.206). This group includes human GBE1, yeast GLC3, Dictyostelium GlgB, Arabidopsis SBEs, and E. coli GlgB. Most organisms have just one representative of this group, with the exception of land plants and green algae (Viridiplantae, “green plants”) and the photosynthetic cyanobacteria that tend to contain multiple paralogs from this group. For instance, in land plants three sub-groups of starch branching enzymes exist, usually named SBE1, SBE2, and SBE3. However, not all plants possess one member of each subgroup; for example, Arabidopsis underwent a recent duplication of SBE2, resulting in SBE2-1 and SBE2-2, and also has one SBE3 member, but lacks a representative of SBE1. Careful phylogenetic analysis shows that this group, due to at least two ancient gene duplications, one of which occurred pre-LUCA (last universal common ancestor) and one pre-LECA (last eukaryotic common ancestor), followed by lineage specific gene losses, has to be divided into a minimum of three sub-groups of orthologous proteins (labeled A, B, and C in Figure 4). While the existence of a separate, plant-specific subgroup containing SBE3 has been reported previously ,, our results show that the well-studied E. coli branching enzyme GlgB, together with GlgB from the cyanobacterium Nostoc punctiforme, are clearly not orthologous to human GBE1, yeast GLC3, Dictyostelium GlgB, and plant SBE1, SBE 2, and SBE 3, but instead represent a branch that emerged by an ancient duplication and was lost in most eukaryotes. On the other hand, other bacterial branching enzymes, such as the one of Bacteroides thetaiotaomicron (UniProt: Q8A9P4), are indeed orthologous to human GBE1. We employed the RIO approach (Resampled Inference of Orthologs) , on MrBayes  output gene trees to confirm these findings. In short, this approach allows to calculate the probability of orthology relationships by integrating orthology assignments over a distribution of gene trees (produced by MrBayes, in this case) . According to this, the posterior probability of E. coli GlgB being orthologous to human GBE1 is 0.0, whereas the posterior probability of Bacteroides thetaiotaomicron Q8A9P4 of being orthologous to human GBE1 is 1.0. These results are further supported by analysis of conserved residues, as described below (see Table 1). Finally, it is likely that this group was affected by even more basal gene duplications, but due to relatively poor phylogenetic resolution at the base of this sub-tree, this remains speculative at this moment.
Enzymes from the second main group, with CBM_48—Alpha-amylase architectures, are only found in land plants and green algae, red algae (Rhodophyta), and some bacterial and archaeal species. Phylogenetic analysis further subdivides this group into two sub-groups of orthologous proteins, correlating with their annotated functions (100% support based on Bayesian, ML, and distance based methods). Enzymes in one sub-group perform debranching functions in glycogen and starch catabolic pathways (EC 220.127.116.11). Similar to branching enzymes, these enzymes underwent expansion in land plants. For example, Arabidopsis and Oryza sativa japonica (rice) contain three paralogs—Isoamylase 1, 2, and 3. The E. coli glycogen debranching enzyme GlgX is a member of this group as well. Pullanases (EC 18.104.22.168) form the second sub-group. These enzymes are found in land plants and green algae, red algae, and Bacteria (we were unable to detect any likely archaeal orthologs in our set of complete genomes) and have additional domains (Bacterial pullanase-associated domain, PUD , or DUF3372). The significantly different lengths of the Alpha_amylase domain of different species depicted in Figure 4 are likely artifacts of the Pfam HMM used to identify them. The three dimensional structure of the human glycogen branching enzyme shows that Alpha_amylase occupies most of the space between CBM_48 and Alpha_amylase_C (Figure 2). This is likely to be the case in all species.
GDE_C domain containing debranching enzymes
Phylogenetic analysis of the debranching enzymes paints a very different picture from that of branching enzymes (see Figures 1 and 5). These proteins are found in animals and fungi and their relatives (members of the Unikonta) and in certain single-celled eukaryotes form the Bikonta group (such as Paramecium tetraurelia), as well as in Bacteria and Archaea, but are not present in land plants, green algae, and Rhodophyta (which are all members of the group Archaeplastida). They exhibit diverse domain architectures, especially between eukaryotes and Bacteria/Archaea. In eukaryotes, the Pfam domain architecture is generally hGDE_N—hGDE-amylase—hGDE_central—GDE_C, whereas bacterial and archaeal enzymes are much shorter and have a different domain architecture, with a GDE_N domain substituting for the three N-terminal domains of the eukaryotic enzymes, resulting in a GDE_N—GDE_C arrangement. As mentioned above, no direct experimental evidence for the function of these bacterial and archaeal proteins currently exists, although because of their distant homology to eukaryotic debranching enzymes and their genomic distribution, we speculate that they function as debranching enzymes. Furthermore, no three-dimensional structure of any protein from this group of bacterial and archaeal GDE_N—GDE_C enzymes is available as of this writing, and no reliable predictions can be made about possible relationships between hGDE_central and GDE_N domains and any other protein domains. A preliminary of analysis of the draft genome of Cyanophora paradoxa indicates that this representative of the Glaucophyta (a small group of freshwater algae  which, together with Rhodophyta, are estimated to be the earliest branching members of Archaeplastida ) contains a putative debranching enzyme with a hGDE_central—GDE_C architecture (see Additional file 3) and thus is conceivable to have a pattern of branching/debranching enzymes dissimilar to that of other Archaeplastida. More genomic data from Glaucophyta will be needed to precisely determine where and when during Archaeplastida evolution the loss of the hGDE-amylase, hGDE_central, GDE_C, and hGDE_N domains occurred.
Distribution of branching and debranching enzymes in major groups of eukaryotes
We also investigated the distribution of branching and debranching enzymes over all major groups of eukaryotes with at least one completely sequenced genome (see Figure 6 for percentages, and Figure 1 for a simplified overview). The result is that branching and debranching enzymes can be found in all major groups of eukaryotes. The only possible exception to this is Rhizaria (a large group of mostly unicellular eukaryotes ), even though with only two completely sequenced genomes in this group, a conclusive answer is impossible at this point. Animals (and their closest relatives, the single-celled choanoflagelates), land plants, and green algae have the highest percentage of genomes with both branching and debranching enzymes (we suspect that the real number is close to 100%; missing enzymes in either category in some animal and plant genomes are mostly likely due to sequencing, assembly, and gene prediction errors and do not represent actual gene losses, as the “losses” appear randomly). For fungi and Amoebozoa, these percentages are lower but are still above 60% and 80%, respectively. On the other hand, the majority of the Alveolata, stramenopiles, and Excavata lack both enzymes but still contain some species with both enzymes (see Additional 4). For other groups, due to limited genomes sequenced, a reliable percentage cannot yet be calculated.
Human disease mutations in GSD3 and GSD4 and their counterparts
Finally (see Table 1), we investigated the amino-acid conservation in the positions mutated in the human glycogen storage diseases GSD3, GSD4, and APBD in orthologs of the human proteins (i.e. from sub-tree “A” for CBM_48—Alpha-amylase branching enzymes; see Figure 4) from a wide variety of species, as well as in the paralogous GlgB enzymes from E. coli and Nostoc punctiforme. Our results show that all disease mutations occur in highly conserved positions, even when compared to species as distantly related as plants and bacteria, stressing the importance of these positions/residues (strong conservation of the glycine at position 1448 in the human glycogen debranching enzyme, mutated in GSD3, has been noted previously ). On the other hand, in the case of CBM_48—Alpha-amylase branching enzymes, this conservation is not maintained in the paralogous enzymes from E. coli and Nostoc punctiforme. We have no explanation for unexpectedly low conversion in these positions in the putative enzyme from Tetrahymena thermophila, especially since our phylogenetic analysis reveals nothing unusual in its sequence.
Branching enzymes (EC 22.214.171.124) with a CBM 48—Alpha-amylase—Alpha-amylase C architecture are present in all the major group of eukaryotes, as well as in Archaea and Bacteria, and are therefore likely to have been present in the last universal common ancestor (LUCA). While they are found in the vast majority of all animal and plant genomes (including green algae) sequenced so far, and are fairly common in fungi, many individual species of single-celled eukaryotes lack identifiable homologs of these enzymes, likely due to gene loss.
For debranching enzymes (EC 126.96.36.199 and EC 188.8.131.52) with a CBM 48—Alpha-amylase architecture, the distribution is very different. On the eukaryotic side, these enzymes are limited to plants and green algae (for which they are found in the vast majority of all sequenced genomes). They are also fairly widespread in Bacteria and Archaea. In contrast, non-homologous debranching enzymes (EC 184.108.40.206) containing the GDE_C domain can be found in species from all kingdoms of life, except green plants and algae.
Comparing these two families allows us to conclude, that in plants, GDE_C-containing debranching enzymes have been replaced by CBM48-containing enzymes. In certain Bacteria (e.g. Nostoc punctiforme) both types of debranching enzymes exist in parallel.
The only major eukaryotic groups for which we are unable to make reliable conclusions are Rhizaria and Glaucophyta. Rhizaria is the only major group of eukaryotes with at least two fully sequenced genomes for which we were unable to detect any glycogen/starch branching or debranching enzymes (such enzymes could be found neither in the two completely sequenced species from Rhizaria, Bigelowiella natans and Reticulomyxa filosa, nor by searching for homologs in Rhizaria in UniProt and Genbank). Despite their importance in the study of plant evolution (due to their placement at the root of the Archaeplastida sub-tree), only one draft genome for Glaucophyta has been released as of this writing; preventing us from conclusively determining whether the pattern of branching/debranching enzymes in this group of alga is indeed not like that of plants and more like that of the rest of eukaryotes (as our preliminary results indicate).
As for human glycogen branching enzyme GBE1, the evolutionary history of this protein can be traced back to Bacteria, for a putative ortholog of human branching enzyme exists in several bacteria, for instance, in a dominant human gut symbiont—Bacteroides thetaiotaomicron. In contrast, our study shows that the well-known E. coli branching enzyme GlgB is not an ortholog of its human homolog, but a member of a separate branch that has been lost in many eukaryotes (including mammals). This obeservation, combined with the low conservation of residues mutated in human diseases in E.coli GlgB, has implications against its use as a model for studying human GSD4/APBD. On the other hand, the Bacteroides thetaiotaomicron branching enzyme is an attractive target for modeling human glycogen storage diseases in Bacteria.
Finally, these results show that not only regulatory proteins, such as those involved in apoptosis regulation , but also basic metabolic enzymes may have a complex evolutionary history, rich in ancient and recent gene duplications, combined with lineage specific gene losses and dynamic domain architectures, with frequent and surreptitious addition and loss of individual domains. Such a history can only be revealed by explicit phylogenetic and comparative domain architecture analysis.
Protein predictions for organisms with a completely sequenced genome were obtained from the sources listed in Additional file 1, covering the following species: 93 animals, 2 choanoflagellates, Capsaspora owczarzaki, Sphaeroforma arctica, 78 fungi, Fonticula alba, 7 amoebozoans, 30 land plants, 10 green algae, Cyanidioschyzon merolae, Cyanophora paradoxa (draft), Emiliania huxleyi, Guillardia theta, 12 Alveolata, 16 Stramenopiles, Bigelowiella natans, Reticulomyxa filose, Thecamonas trahens, 8 Excavata, 49 Archaea, and 133 Bacteria. Hmmscan from HMMER 3.0  together with HMMs for Carbohydrate-binding module 48 (Isoamylase N-terminal domain, CBM_40, PF02922) and Amylo-alpha-1,6-glucosidase (GDE_C, PF06202) from Pfam 27.0  were used to extract putative branching and debranching enzyme sequences. We experimented with different per-domain E-value thresholds to ensure that the results presented here are robust and not simply artifacts of an arbitrarily chosen threshold. For the phylogenetic analyses we generally used a per-domain E-value threshold of 10-3 (unless noted otherwise).
Multiple sequence alignments
Multiple sequence alignments were calculated using MAFFT 7.017b (with “localpair” and “maxiterate 1000” options) and ProbCons 1.12 (default options) . Prior to phylogenetic inference, multiple sequence alignment columns with more than 50% gaps were deleted; for comparison we also performed the analyses based on alignments for which we only deleted columns with more than 90% gaps.
Distance-based minimal evolution trees were inferred by FastME 2.0  (with balanced tree swapping and “GME” initial tree options) based on pairwise distances calculated by TREE-PUZZLE 5.2  (using the WAG substitution model  as recommended by PROTTEST 1.4 , a uniform model of rate heterogeneity, estimation of amino acid frequencies from the dataset, and approximate parameter estimation using a Neighbor-joining tree). For maximum likelihood and Bayesian approaches we employed PhyML 2.4.4  (using 100 bootstrapped data sets, the WAG substitution model, 4 substitution rate categories, estimated proportion of invariable sites, estimated Gamma distribution parameter, and an initial tree calculated by the BIONJ algorithm) and MrBayes 3.2.2  (with 106 generations, a sample frequency of 100, a mixture of amino-acid models with fixed rate matrices and equal rates, and 25% burn-in). For the calculations of typed support values from different sources, confadd 1.01 was used . Tree and domain composition diagrams were drawn using Archaeopteryx . All conclusions presented in this work are robust relative to the alignment methods, the alignment processing, the phylogeny reconstruction methods, and the parameters used. All sequence, alignment, and phylogeny files are available upon request.
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad repository, doi:10.5061/dryad.34vq1, http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.34vq1, in phyloXML format .
CZ participated in the conception and design of the study, performed the data collection, sequence analyses, phylogenetic calculations, and contributed to the interpretation of the results. AG participated in the conception and design of the study, performed the structural analyses, and contributed to the interpretation of the results. Both authors contributed to the writing of the manuscript and read and approved the final text.
Adult polyglucosan body disease
Carbohydrate-active enzymes database
Glycogen branching enzyme
Glycogen debranching enzyme
Glycosyl hydrolase domain
Glycogen storage disease
Last eukaryotic common ancestor
Last universal common ancestor
Starch branching enzyme
Voet D, Voet JG: Biochemistry. In 4th edition. New York: Wiley; 2010.
Magistretti PJ, Allaman I: Glycogen: a Trojan horse for neurons. Nat Neurosci. 2007, 10 (11): 1341-1342. 10.1038/nn1107-1341.
Ball SG, Morell MK: From bacterial glycogen to starch: understanding the biogenesis of the plant starch granule. Annu Rev Plant Biol. 2003, 54: 207-233. 10.1146/annurev.arplant.54.031902.134927.
Meléndez R, Meléndez-Hevia E, Canela EI: The fractal structure of glycogen: a clever solution to optimize cell metabolism. Biophys J. 1999, 77 (3): 1327-1332. 10.1016/S0006-3495(99)76982-1.
Roach PJ: Glycogen and its metabolism. Curr Mol Med. 2002, 2: 101-120. 10.2174/1566524024605761.
Roach PJ, Depaoli-Roach A, Hurley TD, Tagliabracci VS: Glycogen and its metabolism: some new developments and old themes. Biochem J. 2012, 441: 763-787. 10.1042/BJ20111416.
Wilson W, Roach PJ, Montero M, Baroja-Fernández E, Muñoz FJ, Eydallin G, Viale AM, Pozueta-Romero J: Regulation of glycogen metabolism in yeast and bacteria. FEMS Microbiol Rev. 2010, 34: 952-985.
Greenberg CC, Jurczak MJ, Danos AM, Brady MJ: Glycogen branches out: new perspectives on the role of glycogen metabolism in the integration of metabolic pathways. Am J Physiol Endocrinol Metab. 2006, 291: E1-E8. 10.1152/ajpendo.00652.2005.
Bao Y, Yang BZ, Dawson TL, Chen YT: Isolation and nucleotide sequence of human liver glycogen debranching enzyme mRNA: identification of multiple tissue-specific isoforms. Gene. 1997, 197: 389-398. 10.1016/S0378-1119(97)00291-6.
Burki F, Pawlowski J: Monophyly of Rhizaria and Multigene Phylogeny of unicellular Bikonts. Mol Biol Evol. 2006, 23 (10): 1922-1930. 10.1093/molbev/msl055.
Roger AJ, Simpson AGB: Evolution: revisiting the root of the eukaryote tree. Curr Biol. 2009, 19: R165-R167. 10.1016/j.cub.2008.12.032.
Shalchian-Tabrizi K, Minge M, Espelund M, Orr R, Ruden T, Jakobsen KS, Cavalier-Smith T: Multigene phylogeny of choanozoa and the origin of animals. PLoS One. 2008, 3: e2098-10.1371/journal.pone.0002098.
Ward TL, Valberg SJ, Adelson DL, Abbey CA, Binns MM, Mickelson JR: Glycogen branching enzyme (GBE1) mutation causing equine glycogen storage disease IV. Mamm Genome. 2004, 15 (7): 570-577. 10.1007/s00335-004-2369-1.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B: The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42 (D1): D490-D495. 10.1093/nar/gkt1178.
Santos CR, Tonoli CCC, Trindade DM, Betzel C, Takata H, Kuriki T, Kanai T, Imanaka T, Arni RK, Murakami MT: Structural basis for branching-enzyme activity of glycoside hydrolase family 57: Structure and stability studies of a novel branching enzyme from the hyperthermophilic archaeon Thermococcus Kodakaraensis KOD1. Proteins. 2011, 79 (2): 547-557. 10.1002/prot.22902.
Palomo M, Pijning T, Booiman T, Dobruchowska JM, van der Vlist J, Kralj S, Planas A, Loos K, Kamerling JP, Dijkstra BW, van der Maarel MJ, Dijkhuizen L, Leemhuis H: Thermus thermophilus glycoside hydrolase family 57 branching enzyme. J Biol Chem. 2011, 286 (5): 3520-3530. 10.1074/jbc.M110.179515.
Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer EL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (Database issue): D290-D301. 10.1093/nar/gkr1065.
Katsuya Y, Mezaki Y, Kubota M, Matsuura Y: Three-dimensional structure of Pseudomonas isoamylase at 2.2 Å resolution. J Mol Biol. 1998, 281 (5): 885-897. 10.1006/jmbi.1998.1992.
Larson SB, Greenwood A, Cascio D, Day J, McPherson A: Refined molecular structure of pig pancreatic α-amylase at 2 • 1 Å resolution. J Mol Biol. 1994, 235 (5): 1560-1584. 10.1006/jmbi.1994.1107.
Han Y, Sun F-J, Rosales-Mendoza S, Korban SS: Three orthologs in rice, Arabidopsis, and Populus encoding starch branching enzymes (SBEs) are different from other SBE gene families in plants. Gene. 2007, 401: 123-130. 10.1016/j.gene.2007.06.026.
Thon VJ, Vigneron-Lesens C, Marianne-Pepin T, Montreuil J, Decq A, Rachez C, Ball SG, Cannon JF: Coordinate regulation of glycogen metabolism in the yeast Saccharomyces cerevisiae. Induction of glycogen branching enzyme. J Biol Chem. 1992, 267 (21): 15224-15228.
Baecker PA, Greenberg E, Preiss J: Biosynthesis of bacterial glycogen. Primary structure of Escherichia coli 1,4-alpha-D-glucan:1,4-alpha-D-glucan 6-alpha-D-(1, 4-alpha-D-glucano)-transferase as deduced from the nucleotide sequence of the glg B gene. J Biol Chem. 1986, 261 (19): 8738-8743.
Abad MC, Binderup K, Rios-Steiner J, Arni RK, Preiss J, Geiger JH: The X-ray crystallographic structure of Escherichia coli branching enzyme. J Biol Chem. 2002, 277 (44): 42164-42170. 10.1074/jbc.M205746200.
Pal K, Kumar S, Sharma S, Garg SK, Alam MS, Xu HE, Agrawal P, Swaminathan K: Crystal structure of full-length Mycobacterium tuberculosis H37Rv glycogen branching enzyme: insights of N-terminal beta-sandwich in substrate specificity and enzymatic activity. J Biol Chem. 2010, 285 (27): 20897-20903. 10.1074/jbc.M110.121707.
Noguchi J, Chaen K, Vu NT, Akasaka T, Shimada H, Nakashima T, Nishi A, Satoh H, Omori T, Kakuta Y, Kimura M: Crystal structure of the branching enzyme I (BEI) from Oryza sativa L with implications for catalysis and substrate binding. Glycobiology. 2011, 21 (8): 1108-1116. 10.1093/glycob/cwr049.
Maruta K, Mitsuzumi H, Nakada T, Kubota M, Chaen H, Fukuda S, Sugimoto T, Kurimoto M: Cloning and sequencing of a cluster of genes encoding novel enzymes of trehalose biosynthesis from thermophilic archaebacterium Sulfolobus acidocaldarius. Biochim Biophys Acta. 1996, 1291 (3): 177-181. 10.1016/S0304-4165(96)00082-7.
Nguyen DHD, Park J-T, Shim J-H, Tran PL, Oktavina EF, Nguyen TLH, Lee S-J, Park C-S, Li D, Park S-H, Stapleton D, Lee JS, Park KH: Reaction kinetics of substrate transglycosylation catalyzed by TreX of sulfolobus solfataricus and effects on glycogen breakdown. J Bacteriol. 2014, 196 (11): 1941-1949. 10.1128/JB.01442-13.
Park H-S, Park J-T, Kang H-K, Cha H, Kim D-S, Kim J-W, Park K-H: TreX from sulfolobus solfataricus ATCC 35092 displays isoamylase and 4-alpha-glucanotransferase activities. Biosci Biotechnol Biochem. 2007, 71 (5): 1348-1352. 10.1271/bbb.70016.
Woo E-J, Lee S, Cha H, Park J-T, Yoon S-M, Song H-N, Park K-H: Structural insight into the bifunctional mechanism of the glycogen-debranching enzyme TreX from the archaeon sulfolobus solfataricus. J Biol Chem. 2008, 283 (42): 28641-28648. 10.1074/jbc.M802560200.
Nakayama A, Yamamoto K, Tabata S: Identification of the catalytic residues of bifunctional glycogen debranching enzyme. J Biol Chem. 2001, 276 (31): 28824-28828. 10.1074/jbc.M102192200.
Romeo T, Kumar A, Preiss J: Analysis of the Escherichia coli glycogen gene cluster suggests that catabolic enzymes are encoded among the biosynthetic genes. Gene. 1988, 70 (2): 363-376. 10.1016/0378-1119(88)90208-9.
Song H-N, Jung T-Y, Park J-T, Park B-C, Myung PK, Boos W, Woo E-J, Park K-H: Structural rationale for the short branched substrate specificity of the glycogen debranching enzyme GlgX. Proteins. 2010, 78 (8): 1847-1855.
Anantharaman V, Aravind L, Koonin EV: Emergence of diverse biochemical activities in evolutionarily conserved structural scaffolds of proteins. Curr Opin Chem Biol. 2003, 7: 12-20. 10.1016/S1367-5931(02)00018-2.
Nagano N, Orengo C, Thornton JM: One fold with many functions: the evolutionary relationships between TIM barrel families based on their sequences, structures and functions. J Mol Biol. 2002, 321: 741-765. 10.1016/S0022-2836(02)00649-6.
Yang H, Liu MY, Romeo T: Coordinate genetic regulation of glycogen catabolism and biosynthesis in Escherichia coli via the CsrA gene product. J Bacteriol. 1996, 178 (4): 1012-1017.
Streb S, Delatte T, Fau-Umhang M, Umhang M, Fau-Eicke S, Eicke S, Fau-Schorderet M, Schorderet M, Fau-Reinhardt D, Reinhardt D, Fau-Zeeman SC, Zeeman SC: Starch granule biosynthesis in Arabidopsis is abolished by removal of all debranching enzymes but restored by the subsequent removal of an endoamylase. Plant Cell. 2008, 20: 1040-4651. 10.1105/tpc.108.063487.
Lucchiari S, Fogh I, Prelle A, Parini R, Bresolin N, Melis D, Fiori L, Scarlato G, Comi GP: Clinical and genetic variability of glycogen storage disease type IIIa: seven novel AGL gene mutations in the Mediterranean area. Am J Med Genet. 2002, 109 (3): 183-190. 10.1002/ajmg.10347.
Glycogen Storage Disease Type III: Glycogen Storage Disease Type III.., [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/books/NBK26372/]
Andersen DH: Familial cirrhosis of the liver with storage of abnormal glycogen. Lab Invest. 1956, 5 (1): 11-20.
Glycogen Storage Disease Type IV: Glycogen Storage Disease Type IV.., [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/books/NBK115333/]
Moses SW, Parvari R: The variable presentations of glycogen storage disease type IV: a review of clinical, enzymatic and molecular studies. Curr Mol Med. 2002, 2 (2): 177-188. 10.2174/1566524024605815.
Burrow TA, Hopkin RJ, Bove KE, Miles L, Wong BL, Choudhary A, Bali D, Li SC, Chen Y-T: Non-lethal congenital hypotonia due to glycogen storage disease type IV. Am J Med Genet A. 2006, 140A (8): 878-882. 10.1002/ajmg.a.31166.
Bruno C, Servidei S, Shanske S, Karpati G, Carpenter S, McKee D, Barohn RJ, Hirano M, Rifai Z, Dimauro S: Glycogen branching enzyme deficiency in adult polyglucosan body disease. Ann Neurol. 1993, 33 (1): 88-93. 10.1002/ana.410330114.
Fyfe JC, Kurzhals RL, Hawkins MG, Wang P, Yuhki N, Giger U, Van Winkle TJ, Haskins ME, Patterson DF, Henthorn PS: A complex rearrangement in GBE1 causes both perinatal hypoglycemic collapse and late-juvenile-onset neuromuscular degeneration in glycogen storage disease type IV of Norwegian forest cats. Mol Genet Metab. 2007, 90 (4): 383-392. 10.1016/j.ymgme.2006.12.003.
Rahman S, Regina A, Li Z, Mukai Y, Yamamoto M, Kosar-Hashemi B, Abrahams S, Morell MK: Comparison of starch-branching enzyme genes reveals evolutionary relationships among isoforms. Characterization of a gene for starch-branching enzyme IIa from the wheat genome donor Aegilops tauschii. Plant Physiol. 2001, 125 (3): 1314-1324. 10.1104/pp.125.3.1314.
Nougue O, Corbi J, Ball S, Manicacci D, Tenaillon M: Molecular evolution accompanying functional divergence of duplicated genes along the plant starch biosynthesis pathway. BMC Evol Biol. 2014, 14 (1): 103-10.1186/1471-2148-14-103.
Zeeman SC, Kossmann J, Smith AM: Starch: its metabolism, evolution, and biotechnological modification in plants. Annu Rev Plant Biol. 2010, 61 (1): 209-234. 10.1146/annurev-arplant-042809-112301.
Yeats C, Bentley S, Bateman A: New knowledge from old: in silico discovery of novel protein domains in Streptomyces coelicolor. BMC Microbiol. 2003, 3: 3-10.1186/1471-2180-3-3.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP: MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012, 61 (3): 539-542. 10.1093/sysbio/sys029.
Katoh K, Standley DM: MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013, 30 (4): 772-780. 10.1093/molbev/mst010.
Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18 (3): 502-504. 10.1093/bioinformatics/18.3.502.
Desper R, Gascuel O: Fast and accurate phylogeny reconstruction algorithms based on the minimum-evolution principle. J Comput Biol. 2002, 9 (5): 687-705. 10.1089/106652702761034136.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Zmasek CM, Godzik A: Data from: Phylogenomic Analysis of Glycogen Branching and Debranching Enzymatic Duo. Dryad Data Repository. 2014
Zmasek CM, Eddy SR: RIO: analyzing proteomes by automated phylogenomics using resampled inference of orthologs. BMC Bioinformatics. 2002, 3: 14-10.1186/1471-2105-3-14.
forester: software libraries for evolutionary biology and comparative genomics research.., [https://sites.google.com/site/cmzmasek/home/software/forester]
Zmasek CM, Eddy SR: A simple algorithm to infer gene duplication and speciation events on a gene tree. Bioinformatics. 2001, 17 (9): 821-828. 10.1093/bioinformatics/17.9.821.
Price DC, Chan CX, Yoon HS, Yang EC, Qiu H, Weber APM, Schwacke R, Gross J, Blouin NA, Lane C, Reyes-Prieto A, Durnford DG, Neilson JA, Lang BF, Burger G, Steiner JM, Löffelhardt W, Meuser JE, Posewitz MC, Ball S, Arias MC, Henrissat B, Coutinho PM, Rensing SA, Symeonidi A, Doddapaneni H, Green BR, Rajah VD, Boore J, Bhattacharya D: Cyanophora paradoxa genome elucidates origin of photosynthesis in algae and plants. Science. 2012, 335 (6070): 843-847. 10.1126/science.1213561.
Keeling PJ: Diversity and evolutionary history of plastids and their hosts. Am J Bot. 2004, 91 (10): 1481-1493. 10.3732/ajb.91.10.1481.
Deschamps P, Moreira D: Signal conflicts in the phylogeny of the primary photosynthetic eukaryotes. Mol Biol Evol. 2009, 26 (12): 2745-2753. 10.1093/molbev/msp189.
Brown MW, Kolisko M, Silberman JD, Roger AJ: Aggregative multicellularity evolved independently in the eukaryotic supergroup Rhizaria. Current Biol: CB. 2012, 22 (12): 1123-1127. 10.1016/j.cub.2012.04.021.
Cheng A, Zhang M, Okubo M, Omichi K, Saltiel AR: Distinct mutations in the glycogen debranching enzyme found in glycogen storage disease type III lead to impairment in diverse cellular functions. Hum Mol Genet. 2009, 18 (11): 2045-2052. 10.1093/hmg/ddp128.
Zmasek CM, Zhang Q, Ye Y, Godzik A: Surprising complexity of the ancestral apoptosis network. Genome Biol. 2007, 8 (10): R226-10.1186/gb-2007-8-10-r226.
Eddy SR: Accelerated Profile HMM Searches. PLoS Comput Biol. 2011, 7 (10): e1002195-10.1371/journal.pcbi.1002195.
Do CB, Mahabhashyam MS, Brudno M, Batzoglou S: ProbCons: probabilistic consistency-based multiple sequence alignment. Genome Res. 2005, 15 (2): 330-340. 10.1101/gr.2821705.
Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18 (5): 691-699. 10.1093/oxfordjournals.molbev.a003851.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.
Han M, Zmasek CM: phyloXML: XML for evolutionary biology and comparative genomics. BMC Bioinformatics. 2009, 10 (1): 356-10.1186/1471-2105-10-356.
Nozaki H, Takano H, Misumi O, Terasawa K, Matsuzaki M, Maruyama S, Nishida K, Yagisawa F, Yoshida Y, Fujiwara T, Takio S, Tamura K, Chung SJ, Nakamura S, Kuroiwa H, Tanaka K, Sato N, Kuroiwa T: A 100%-complete sequence reveals unusually simple genomic features in the hot-spring red alga Cyanidioschyzon merolae. BMC Biol. 2007, 5: 28-10.1186/1741-7007-5-28.
This work was supported by the Human Frontier Science Project (grant number RGP0027/2011). The authors acknowledge the sequencing centers listed in Additional file 1 for their efforts in sequencing, assembling, and annotating the genomes analyzed in this study.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 2: Protein identifiers (mostly from UniProt; for others, see legend for Additional file 4) for branching and debranching enzymes with CBM_48 and Alpha-amylase domains are listed (per-domain E-value cutoff: 10-3). For eukaryotic enzymes, taxonomic groups (such as Alveolata) are indicated. Simplified domain architecture overviews are given for each enzyme (“~” is used to indicate linkers between domains shorter than 11aa, whereas “----”stands for linkers longer than 10aa). Individual E-values for CBM_48 and Alpha-amylase domains are shown as well. (XLSX 73 KB)
Additional file 3: Protein identifiers (mostly from UniProt; for others, see legend for Additional file 4) for debranching enzymes with hGDE-amylase and GDE_C domains are listed (per-domain E-value cutoff: 10-3). Taxonomic groups (such as Alveolata) are indicated. Simplified domain architecture overviews are given for each enzyme (“~” is used to indicate linkers between domains shorter than 11aa, whereas “----” stands for linkers longer than 10aa). Individual E-values for hGDE-amylase and GDE_C domains are shown as well. (XLSX 31 KB)
Additional file 4: Protein identifiers are generally from the UniProt database, except for those marked with an asterisk, which are from GenBank, and those from Sphaeroforma arctica, Thecamonas trahens, and Fonticula alba which originate from the Origins of Multicellularity Sequencing Project (Broad Institute of Harvard and MIT: http://www.broadinstitute.org), and those from Cyanidioschyzon merolae which are from the National Institute of Genetics, Japan . All examples are from completely sequenced genomes (see Additional file 1). (XLSX 18 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zmasek, C.M., Godzik, A. Phylogenomic analysis of glycogen branching and debranching enzymatic duo. BMC Evol Biol 14, 183 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-014-0183-2
- Glycogen storage disease