Complete genome sequence of the cellulolytic planctomycete Telmatocola sphagniphila SP2T and characterization of the first cellulolytic enzyme from planctomycetes

Planctomycetes of the family Gemmataceae are strictly aerobic chemo-organotrophs that display a number of hydrolytic capabilities. A member of this family, Telmatocola sphagniphila SP2^T, is the first described planctomycete with experimentally proven ability for growth on cellulose. In this study, the complete genome sequence of strain SP2^T was obtained and the genome-encoded determinants of its cellulolytic potential were analyzed. The T. sphagniphila SP2^T genome was 6.59 Mb in size and contained over 5200 potential protein-coding genes. The search for enzymes that could be potentially involved in cellulose degradation identified a putative cellulase that contained a domain from the GH44 family of glycoside hydrolases. Homologous enzymes were also revealed in the genomes of two other Gemmataceae planctomycetes, Zavarzinella formosa A10^T and Tuwongella immobilis MBLW1^T. The gene encoding this predicted cellulase in strain SP2^T was expressed in E. coli and the hydrolytic activity of the recombinant enzyme was confirmed in tests with carboxymethyl cellulose but not with crystalline cellulose, xylan, mannan or laminarin. This is the first experimentally characterized cellulolytic enzyme from planctomycetes.

Frigoriglobus tundricola gen. nov., sp. nov., a psychrotolerant cellulolytic planctomycete of the family Gemmataceae from a littoral tundra wetland

Planctomycetes of the family Gemmataceae are characterized by large genome sizes and cosmopolitan distribution in freshwater and terrestrial environments but their ecological functions remain poorly understood. In this study, we characterized a novel representative of this family, strain PL17T, which was isolated from a littoral tundra wetland and was capable of growth on xylan and cellulose. Cells of this isolate were represented by pink-pigmented spheres that multiplied by budding and occurred singly or in short chains and aggregates. Strain PL17T was obligately aerobic, mildly acidophilic chemoorganotrophic bacterium, which displayed good tolerance of low temperatures. The major fatty acids were C18:0, C16:1ω5, and βOH-C16:1; the major polar lipid was trimethylornithine. The genome of strain PL17T consisted of a 9.83 Mb chromosome and a 24.69kb plasmid. The G+C contents of the chromosomal and plasmid DNA were 67.4 and 62.3mol%, respectively. Over 8900 potential protein-coding genes were identified in the genome including a putative cellulase that contains a domain from the GH5 family of glycoside hydrolases. The genome of strain PL17T contained one linked and one unlinked rRNA operons with 16S rRNA gene sequences displaying 94.5% similarity to that in Gemmata obscuriglobus UQM2246T. Based on the results of comparative phenotypic, chemotaxonomic and phylogenomic analyses, we propose to classify strain PL17T (= CECT 9407T=VKM B-3467T) as representing a novel genus and species of the family Gemmataceae, Frigoriglobus tundricola gen. nov., sp. nov.

Limnoglobus roseus gen. nov., sp. nov., a novel freshwater planctomycete with a giant genome from the family Gemmataceae

The family Gemmataceae accommodates aerobic, chemoorganotrophic planctomycetes, which inhabit various freshwater ecosystems, wetlands and soils. Here, we describe a novel member of this family, strain PX52T, which was isolated from a boreal eutrophic lake in Northern Russia. This isolate formed pink-pigmented colonies and was represented by spherical cells that occurred singly, in pairs or aggregates and multiplied by budding. Daughter cells were highly motile. PX52T was an obligate aerobic chemoorganotroph, which utilized various sugars and some heteropolysaccharides. Growth occurred at pH 5.0-7.5 (optimum pH 6.5) and at temperatures between 10 and 30 °C (optimum 20-25 °C). The major fatty acids were C18 : 1ɷ7c, C18 : 0 and βOH-C16:0; the major intact polar lipid was trimethylornithine, and the quinone was MK-6. The complete genome of PX52T was 9.38 Mb in size and contained nearly 8000 potential protein-coding genes. Among those were genes encoding a wide repertoire of carbohydrate-active enzymes (CAZymes) including 33 glycoside hydrolases (GH) and 87 glycosyltransferases (GT) affiliated with 17 and 12 CAZy families, respectively. DNA G+C content was 65.6 mol%. PX52T displayed only 86.0-89.8 % 16S rRNA gene sequence similarity to taxonomically described Gemmataceae planctomycetes and differed from them by a number of phenotypic characteristics and by fatty acid composition. We, therefore, propose to classify it as representing a novel genus and species, Limnoglobus roseus gen. nov., sp. nov. The type strain is strain PX52T (=KCTC 72397T=VKM B-3275T).

[GH10 Family of Glycoside Hydrolases: Structure and Evolutionary Connections]

Evolutionary connections were analyzed for endo-β-xylanases, which possess the GH10 family catalytic domains. A homology search yielded thrice as many proteins as are available from the Carbohydrate-Active Enzymes (CAZy) database. Lateral gene transfer was shown to play an important role in evolution of bacterial proteins of the family, especially in the phyla Acidobacteria, Cyanobacteria, Planctomycetes, Spirochaetes, and Verrucomicrobia. In the case of Verrucomicrobia, 23 lateral transfers from organisms of other phyla were detected. Evolutionary relationships were observed between the GH10 family domains and domains with the TIM-barrel tertiary structure from several other glycosidase families. The GH39 family of glycoside hydrolases showed the closest relationship. Unclassified homologs were grouped into 12 novel families of putative glycoside hydrolases (GHL51-GHL62).

Comparative Genomics of Four Isosphaeraceae Planctomycetes: A Common Pool of Plasmids and Glycoside Hydrolase Genes Shared by Paludisphaera borealis PX4T, Isosphaera pallida IS1BT, Singulisphaera acidiphila DSM 18658T, and Strain SH-PL62

The family Isosphaeraceae accommodates stalk-free planctomycetes with spherical cells, which can be assembled in short chains, long filaments, or aggregates. These bacteria inhabit a wide variety of terrestrial environments, among those the recently described Paludisphaera borealis PX4T that was isolated from acidic boreal wetlands. Here, we analyzed its finished genome in comparison to those of three other members of the Isosphaeraceae: Isosphaera pallida IS1BT, Singulisphaera acidiphila DSM 18658T, and the uncharacterized planctomycete strain SH-PL62. The complete genome of P. borealis PX4T consists of a 7.5 Mb chromosome and two plasmids, 112 and 43 kb in size. Annotation of the genome sequence revealed 5802 potential protein-coding genes of which 2775 could be functionally assigned. The genes encoding metabolic pathways common for chemo-organotrophic bacteria, such as glycolysis, citrate cycle, pentose-phosphate pathway, and oxidative phosphorylation were identified. Several genes involved in the synthesis of peptidoglycan as well as N-methylated ornithine lipids were present in the genome of P. borealis PX4T. A total of 26 giant genes with a size >5 kb were detected. The genome encodes a wide repertoire of carbohydrate-active enzymes (CAZymes) including 44 glycoside hydrolases (GH) and 83 glycosyltransferases (GT) affiliated with 21 and 13 CAZy families, respectively. The most-represented families are GH5, GH13, GH57, GT2, GT4, and GT83. The experimentally determined carbohydrate utilization pattern agrees well with the genome-predicted capabilities. The CAZyme repertoire in P. borealis PX4T is highly similar to that in the uncharacterized planctomycete SH-PL62 and S. acidiphila DSM 18658T, but different to that in the thermophile I. pallida IS1BT. The latter strain has a strongly reduced CAZyme content. In P. borealis PX4T, many of its CAZyme genes are organized in clusters. Contrary to most other members of the order Planctomycetales, all four analyzed Isosphaeraceae planctomycetes have plasmids in numbers varying from one to four. The plasmids from P. borealis PX4T display synteny to plasmids from other family members, providing evidence for their common evolutionary origin.

[Phylogeny of beta-xylanases from Planctomycetes]

Here, we present the results of a computational analysis of a group of hypothetical GH10 endo-beta-xylanases from the Planctomycetes, a bacterial phylum with poorly characterized functional capabilities. These proteins are encoded in all analyzed genomes of heterotrophic Planctomycetes and form a phylogenetically distinct and tight cluster. In addition, we determined nucleotide sequences for endo-beta-xylanase genes from five strains of Isosphaera-Singulisphaera group of the Planctomycetes. The trees constructed for the 16S rRNA genes and the inferred amino acid sequences of endo-beta-xylanases were highly congruent, thus suggesting the vertical transfer of endo-beta-xylanase genes and their functional importance in Planctomycetes.

Hierarchical classification of glycoside hydrolases

This review deals with structural and functional features of glycoside hydrolases, a widespread group of enzymes present in almost all living organisms. Their catalytic domains are grouped into 120 amino acid sequence-based families in the international classification of the carbohydrate-active enzymes (CAZy database). At a higher hierarchical level some of these families are combined in 14 clans. Enzymes of the same clan have common evolutionary origin of their genes and share the most important functional characteristics such as composition of the active center, anomeric configuration of cleaved glycosidic bonds, and molecular mechanism of the catalyzed reaction (either inverting, or retaining). There are now extensive data in the literature concerning the relationship between glycoside hydrolase families belonging to different clans and/or included in none of them, as well as information on phylogenetic protein relationship within particular families. Summarizing these data allows us to propose a multilevel hierarchical classification of glycoside hydrolases and their homologs. It is shown that almost the whole variety of the enzyme catalytic domains can be brought into six main folds, large groups of proteins having the same three-dimensional structure and the supposed common evolutionary origin.

GH101 family of glycoside hydrolases: subfamily structure and evolutionary connections with other families

The GH101 family is composed of endo-alpha-N-acetylgalactosaminidases and their homologues. Pairwise sequence comparison and phylogenetic analysis allowed us to distinguish five to six subfamilies in this family. Diverse domain structures were found among the family members. Usually they have five irreplaceable and some optional domains. Iterative screening of the protein database revealed an evolutionary relationship of the GH101 catalytic domain with glycoside hydrolase domains from GH13, GH31, and GH70 families. Among other homologous proteins we have found representatives of COG1649, as well as members of four new families of predicted glycoside hydrolases (GHL1-GHL4).

GH97 is a new family of glycoside hydrolases, which is related to the alpha-galactosidase superfamily

Background

As a rule, about 1% of genes in a given genome encode glycoside hydrolases and their homologues. On the basis of sequence similarity they have been grouped into more than ninety GH families during the last 15 years. The GH97 family has been established very recently and initially included only 18 bacterial proteins. However, the evolutionary relationship of the genes encoding proteins of this family remains unclear, as well as their distribution among main groups of the living organisms.

Results

The extensive search of the current databases allowed us to double the number of GH97 family proteins. Five subfamilies were distinguished on the basis of pairwise sequence comparison and phylogenetic analysis. Iterative sequence analysis revealed the relationship of the GH97 family with the GH27, GH31, and GH36 families of glycosidases, which belong to the α-galactosidase superfamily, as well as a more distant relationship with some other glycosidase families (GH13 and GH20).

Conclusion

The results of this study show an unexpected sequence similarity of GH97 family proteins with glycoside hydrolases from several other families, that have (β/α)₈-barrel fold of the catalytic domain and a retaining mechanism of the glycoside bond hydrolysis. These data suggest a common evolutionary origin of glycosidases representing different families and clans.

Retrieving sequences of enzymes experimentally characterized but erroneously annotated : the case of the putrescine carbamoyltransferase

Background

Annotating genomes remains an hazardous task. Mistakes or gaps in such a complex process may occur when relevant knowledge is ignored, whether lost, forgotten or overlooked. This paper exemplifies an approach which could help to ressucitate such meaningful data.

Results

We show that a set of closely related sequences which have been annotated as ornithine carbamoyltransferases are actually putrescine carbamoyltransferases. This demonstration is based on the following points : (i) use of enzymatic data which had been overlooked, (ii) rediscovery of a short NH₂-terminal sequence allowing to reannotate a wrongly annotated ornithine carbamoyltransferase as a putrescine carbamoyltransferase, (iii) identification of conserved motifs allowing to distinguish unambiguously between the two kinds of carbamoyltransferases, and (iv) comparative study of the gene context of these different sequences.

Conclusions

We explain why this specific case of misannotation had not yet been described and draw attention to the fact that analogous instances must be rather frequent. We urge to be especially cautious when high sequence similarity is coupled with an apparent lack of biochemical information. Moreover, from the point of view of genome annotation, proteins which have been studied experimentally but are not correlated with sequence data in current databases qualify as "orphans", just as unassigned genomic open reading frames do. The strategy we used in this paper to bridge such gaps in knowledge could work whenever it is possible to collect a body of facts about experimental data, homology, unnoticed sequence data, and accurate informations about gene context.

Three acidic residues are at the active site of a beta-propeller architecture in glycoside hydrolase families 32, 43, 62, and 68

Multiple-sequence alignment of glycoside hydrolase (GH) families 32, 43, 62, and 68 revealed three conserved blocks, each containing an acidic residue at an equivalent position in all the enzymes. A detailed analysis of the site-directed mutations so far performed on invertases (GH32), arabinanases (GH43), and bacterial fructosyltransferases (GH68) indicated a direct implication of the conserved residues Asp/Glu (block I), Asp (block II), and Glu (block III) in substrate binding and hydrolysis. These residues are close in space in the 5-bladed beta-propeller fold determined for Cellvibrio japonicus alpha-L-arabinanase Arb43A [Nurizzo et al., Nat Struct Biol 2002;9:665-668] and Bacillus subtilis endo-1,5-alpha-L-arabinanase. A sequence-structure compatibility search using 3D-PSSM, mGenTHREADER, INBGU, and SAM-T02 programs predicted indistinctly the 5-bladed beta-propeller fold of Arb43A and the 6-bladed beta-propeller fold of sialidase/neuraminidase (GH33, GH34, and GH83) as the most reliable topologies for GH families 32, 62, and 68. We conclude that the identified acidic residues are located at the active site of a beta-propeller architecture in GH32, GH43, GH62, and GH68, operating with a canonical reaction mechanism of either inversion (GH43 and likely GH62) or retention (GH32 and GH68) of the anomeric configuration. Also, we propose that the beta-propeller architecture accommodates distinct binding sites for the acceptor saccharide in glycosyl transfer reaction.

Using Quaternary Structures to Assess the Evolutionary History of Proteins: The Case of the Aspartate Carbamoyltransferase

Many evolutionary scenarios describing the history of proteins are based solely on phylogenetic studies. We have designed a new approach that allows ascertainment of such questionable scenarios by taking into account quaternary structures: we used aspartate carbamoyltransferase (ATCase) as a case study. Prokaryotic ATCases correspond to different classes of quaternary structures according to the mode of association of the catalytic PyrB subunit with other polypeptides, either the PyrI regulatory subunit (class B) or a dihydroorotase (class A), which may be active (PyrC, subclass A1) or inactive (PyrC', subclass A2). Class C is uniquely made up of trimers of PyrB. The PyrB phylogenetic tree is not congruent with the tree of life, but it became coherent when we recognized the existence of two families of ATCases, ATC I and ATC II. Remarkably, a very strong correlation was found between the pattern of PyrB phylogenetic clustering and the different classes of quaternary structures of ATCases. All class B ATCases form a clade in family ATC II, which also contains all eukaryotic sequences. In contrast, family ATC I is made up of classes A and C. These results suggest unexpected common ancestry for prokaryotic B and eukaryotic ATCases on the one hand, and for A and C on the other. Thus, the emergence of specific quaternary structures appears to have been a more recent event than the separation into the ATC I and ATC II families. We propose that different evolutionary constraints, depending on the identity of the partners interacting in the different kinds of holoenzymes, operated in a concerted way on the ancestral pyrB genes and the respective associated genes pyrI or pyrC, so as to maintain appropriate inter-polypeptides interactions at the level of quaternary structure. The process of coevolution of genes encoding proteins interacting in various holoenzymes has been assessed by calculating the correlation coefficient between their respective phylogenetic trees. Our approach integrating data obtained from the separate fields of structural biology and molecular evolution could be useful in other cases where pure statistical data need to receive independent confirmation.

beta-fructosidase superfamily: homology with some alpha-L-arabinases and beta-D-xylosidases

Comparison of the amino acid sequences of four families of glycosyl hydrolases reveals that they are homologous and have several common conserved regions. Two of these families contain beta-fructosidases (glycosyl hydrolase families GH32 and GH68) and the other two include alpha-L-arabinases and beta-xylosidases (families GH43 and GH62). The latter two families are proposed to be grouped together with the former two into the beta-fructosidase (furanosidase) superfamily. Several ORFs can be considered as a fifth family of the superfamily on the basis of sequence similarity. It is shown for the first time that a glycosyl hydrolase superfamily can include enzymes with both inversion and retention mechanism of action. Composition of the active center for enzymes of the superfamily is discussed.

[Changes in the composition of populations of Phytophthora infestans (Mont.) de Bary in Russia in the 1990s based on the results of mitochondrial DNA analysis]

Changes in the composition of the populations of Phytophthora infestans, the agent of late blight of potato and open-ground tomato, occurred in Russia in the early 1990s. These changes were detected via analysis of mitochondrial DNA (mtDNA) and other characteristics in 219 isolates of P. infestans. It was demonstrated that these changes were similar to those that occurred in European populations of the fungus in the early 1980s, when "old" genotypes were rapidly displaced by "new" ones. Possible causes of this phenomenon, which is related to the characteristics of the parasite's biology, and of host specialization are discussed.

Conserved sequence motifs in levansucrases and bifunctional beta-xylosidases and alpha-L-arabinases

Comparison of the amino acid sequences of two families of glycosyl hydrolases reveals that they are related in a region in the central part of the sequences. One of these families (GH family 68) includes levansucrases and the other one (glycosyl hydrolase family 43) includes bifunctional beta-xylosidases and alpha-L-arabinofuranosidases. The similarity of the primary structure of proteins from these families allows us to consider the invariant glutamate residue as a component of their active center. It is shown for the first time that glycosyl hydrolases recognizing different glycofuranoside residues can have a common sequence motif.