Comparative genomic and phylogenetic analyses of Gammaproteobacterial glg genes traced the origin of the Escherichia coli glycogen glgBXCAP operon to the last common ancestor of the sister orders Enterobacteriales and Pasteurellales

• Published in: PloS one 2015-01, Vol.10 (1), p.e0115516-e0115516
• Main Authors: Almagro, Goizeder, Viale, Alejandro M, Montero, Manuel, Rahimpour, Mehdi, Muñoz, Francisco José, Baroja-Fernández, Edurne, Bahaji, Abdellatif, Zúñiga, Manuel, González-Candelas, Fernando, Pozueta-Romero, Javier
• Format: Article
• Language: English
• Subjects: Adenosine diphosphate; Analysis; Bacteria; Biodegradation; Bioinformatics; Biological evolution; Comparative analysis; Copy number; E coli; Enzymes; Escherichia coli; Escherichia coli - genetics; Evolution (Biology); Evolution, Molecular; Evolutionary genetics; Gammaproteobacteria; Gene transfer; Genes; Genetic aspects; Genomes; Genomics; Glucan; Glucose; Glycogen; Glycogen - genetics; Glycogen phosphorylase; Glycogen synthase; Glycogen synthesis; Metabolism; Operon; Pasteurellaceae - genetics; Phosphorylase; Phylogenetics; Phylogeny; Polymers; Proteins; Species; Trends
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0115516

Production of branched α-glucan, glycogen-like polymers is widely spread in the Bacteria domain. The glycogen pathway of synthesis and degradation has been fairly well characterized in the model enterobacterial species Escherichia coli (order Enterobacteriales, class Gammaproteobacteria), in which the cognate genes (branching enzyme glgB, debranching enzyme glgX, ADP-glucose pyrophosphorylase glgC, glycogen synthase glgA, and glycogen phosphorylase glgP) are clustered in a glgBXCAP operon arrangement. However, the evolutionary origin of this particular arrangement and of its constituent genes is unknown. Here, by using 265 complete gammaproteobacterial genomes we have carried out a comparative analysis of the presence, copy number and arrangement of glg genes in all lineages of the Gammaproteobacteria. These analyses revealed large variations in glg gene presence, copy number and arrangements among different gammaproteobacterial lineages. However, the glgBXCAP arrangement was remarkably conserved in all glg-possessing species of the orders Enterobacteriales and Pasteurellales (the E/P group). Subsequent phylogenetic analyses of glg genes present in the Gammaproteobacteria and in other main bacterial groups indicated that glg genes have undergone a complex evolutionary history in which horizontal gene transfer may have played an important role. These analyses also revealed that the E/P glgBXCAP genes (a) share a common evolutionary origin, (b) were vertically transmitted within the E/P group, and (c) are closely related to glg genes of some phylogenetically distant betaproteobacterial species. The overall data allowed tracing the origin of the E. coli glgBXCAP operon to the last common ancestor of the E/P group, and also to uncover a likely glgBXCAP transfer event from the E/P group to particular lineages of the Betaproteobacteria.