Skip to main content

A network-based approach to identify substrate classes of bacterial glycosyltransferases

Abstract

Background

Bacterial interactions with the environment- and/or host largely depend on the bacterial glycome. The specificities of a bacterial glycome are largely determined by glycosyltransferases (GTs), the enzymes involved in transferring sugar moieties from an activated donor to a specific substrate. Of these GTs their coding regions, but mainly also their substrate specificity are still largely unannotated as most sequence-based annotation flows suffer from the lack of characterized sequence motifs that can aid in the prediction of the substrate specificity.

Results

In this work, we developed an analysis flow that uses sequence-based strategies to predict novel GTs, but also exploits a network-based approach to infer the putative substrate classes of these predicted GTs. Our analysis flow was benchmarked with the well-documented GT-repertoire of Campylobacter jejuni NCTC 11168 and applied to the probiotic model Lactobacillus rhamnosus GG to expand our insights in the glycosylation potential of this bacterium. In L. rhamnosus GG we could predict 48 GTs of which eight were not previously reported. For at least 20 of these GTs a substrate relation was inferred.

Conclusions

We confirmed through experimental validation our prediction of WelI acting upstream of WelE in the biosynthesis of exopolysaccharides. We further hypothesize to have identified in L. rhamnosus GG the yet undiscovered genes involved in the biosynthesis of glucose-rich glycans and novel GTs involved in the glycosylation of proteins. Interestingly, we also predict GTs with well-known functions in peptidoglycan synthesis to also play a role in protein glycosylation.

Background

The glycome, playing a crucial role in allowing bacteria to establish environment- and host-specific interactions [1, 2] consists of a wide variety of glycoconjugates, i.e. glycans being covalently linked to other macromolecules. In Gram-negatives, these glycoconjugates occur mainly in the outer membrane as a thin layer of peptidoglycan (PG) and lipopolysaccharides (LPS) or lipo-oligosaccharides (LOS). Across the outer membrane, exopolysaccharides (EPS) or capsular polysaccharides (CPS), glycoproteins and glycolipids can further decorate the cell surface [2]. In Gram-positives, which in contrast to Gram-negatives lack an outer membrane, complex polymers such as teichoic acids in Firmicutes and lipoglycans in Actinobacteria strengthen a thick layer of PG. CPS or EPS are also often found as most external layer in Gram-positive bacteria. Bacteria can also produce intracellular glycoconjugates, such as glycosylated secondary metabolites and storage polysaccharides like glycogen [2].

Glycosyltransferases (GTs), transferring sugar moieties from an activated donor to a specific substrate [3], are key enzymes in the biosynthesis of glycoconjugates. Depending on their specificity, the substrates of GTs range from lipids, proteins, saccharides, nucleic acids to small molecules [3]. In bacteria, two different glycosylation mechanisms have been described: sequential glycosylation, in which either soluble or membrane-associated GTs transfer glycan monomers directly to the final substrate and en bloc glycosylation, in which the sugar moiety is first assembled and only then transferred to the final substrate by an specialized GT (oligosaccharyltransferase (OST) or polymerase) [4, 5]. The latter mechanism is by far the best documented, and is involved in the biosynthesis of heteropolymeric EPS/CPS, O-antigens in LPS, and even PG biosynthesis, highlighting the commonalities in the biosynthesis of these glycoconjugates [5]. Apart from their general role in glycosylation, the specificities of most of the GTs and the cellular role of their end products are still largely unknown. In addition, most of the substrate specificities of GTs involved in LPS, PG and glycoproteins have been described in Gram-negatives [6, 7], while glycosylation in Gram-positives is much less studied.

Whereas sequence-based predictions have shown useful to identify potential GTs [810], predicting the specificity of those identified GTs is less trivial, definitely for prokaryotes for which no clear sequence motifs determining substrate specificity have been described [11]. In addition, many GTs and OSTs show substrate promiscuity [12, 13], hampering the identification of clear substrate motifs.

To improve the annotation of GTs in prokaryotes, we developed an analysis flow that uses a sequence-based strategy to predict GTs and a network-based approach [14] to identify links between these predicted GTs and other genes/proteins. Although such links do not give insights into the precise biochemical mechanisms of a GT with its substrate, they aid in relating the GT to possible classes of molecules that could accept the sugar moieties from these GTs (referred to as substrate classes).

We tested our analysis flow on the genome of C. jejuni NCTC 11168, in which the important classes of glycoconjugates (N- and O-glycoproteins, PG, LOS, and CPS) are well characterized [4].

Further applying our analysis flow on the probiotic bacterium Lactobacillus rhamnosus GG provided a comprehensive re-annotation of putative GTs in this species, the possible substrate classes of these GTs and their mode of action. These predictions are a very useful resource for experimentalists, predominantly because the study of (protein) glycosylation in lactobacilli and related organisms is not straightforward [15]. Our predictions unveil putative novel mechanisms of (protein) glycosylation, involving the potential, promiscuous role of GTs with known function in PG biosynthesis.

Methods

Bacterial proteomes

The proteomes and current genome annotations of Lactobacillus rhamnosus GG (NC_013198.1) and Campylobacter jejuni NCTC 11168 (NC_002163.1) were obtained from GenBank (http://www.ncbi.nlm.nih.gov/genbank/).

Hidden Markov Model profile searches

Hidden Markov Models (HMMs) describing known GT signatures were collected from SUPERFAMILY (http://supfam.cs.bris.ac.uk/SUPERFAMILY/), CAZy (http://www.cazy.org/) and Pfam (http://pfam.sanger.ac.uk/) and subdivided into three groups depending on their expected specificity for GTs (Table 1). For CAZy, a thorough search of this database was performed, and all the HMMs covering GT classes that had bacterial representatives were included in our analysis (see below).

Table 1 Summary of the Hidden Markov Models (HMMs) used to screen for glycosyltransferases in the proteomes of Campylobacter jejuni NCTC 11168 and Lactobacillus rhamnosus GG

The first and least specific group contains the HMM representing ‘Rossmann-fold domains’, which are known to resemble the GT-A and GT-B folds typical for GTs using sugar nucleotides as donor [3, 8, 16]. A second group comprises the HMMs for ‘Sugar transferases’ and ‘UDP-Glycosyltransferases’ respectively, both HMMs of intermediate specificity covering a broad class of GTs [8, 10]. A last group combines a set of more GT-specific HMMs (10 in total), all of which are based on a small number of family-specific sequences [1726]. This group combines HMMs extracted from CAZy [27], representative for enzymes that catalyze glycosidic bonds (strictu-sensu GTs) with HMMs extracted from Pfam [28] that are representative for non-Leloir GTs that use non-nucleotide sugar donors or oligo/polysaccharides. Enzymes involved in the transfer of the sugar moiety to the final substrate (such as OTases and priming GTs) are examples of this latter class of non-Leloir GTs.

The collected HMMs were used to screen entire proteomes (C. jejuni NCTC 11168 and Lactobacillus rhamnosus GG) with hmmsearch from the HMMER package version 2.2 [29]. Hits were filtered using an E-value cut-off of 0.1.

Protein fold recognition

The profile based fold recognition method pGenTHREADER [30], accessible via the PSIRED server (http://bioinf.cs.ucl.ac.uk/psipred/) was used to detect known GT-A/GT-B folds in proteins predicted to be GTs by the HMM search. Each of the input sequences was aligned against a library of 3D folds based on CATH v3.3 (the Protein Structure Database, available at http://www.cathdb.info/) by pGenTHREADER. The library of 3D folds contains a total of 684 PDB structures of known GTs. Putative GTs were only retained if they predicted fold showed significant homology (net score > 46) to the one of a resolved 3D structure with known GT activity present in the library (refined set). We selected a cutoff > 46 on the net score of pGenTHREADER since any values higher than this threshold are categorized as HIGH to CERTIFIED confidence predictions (default conservative setting of the tool).

Detecting functional partners of glycosyltransferases

The STRING database (http://string-db.org/) was used as the source of functional networks [14, 31]. We interrogated STRING using as queries our predicted GTs from both L. rhamnosus GG and C. jejuni NCTC 11168 to retrieve the network of functional partners associated to each query (query-based subnetwork). We only considered functional interactions with a score higher than 0.7, which is the default value in STRING for high confidence interactions. A total of 1112 functional interactions were retrieved for L. rhamnosus GG, supported by 2338 independent evidences distributed as follows: 1682 evidences based on the genomic context of the interacting partners (e.g. physical closeness, co-occurrence in closely related species, gene fusion events); 153 evidences based on the co-expression of the interacting partners; 28 evidences derived from high-throughput experiments (e.g. protein-protein interaction data); 465 evidences derived from the literature (text-mining). For C. jejuni NCTC 11168 a total of 1727 functional interactions were retrieved supported by 3190 independent evidences from the following data sources: 2520 evidences based on the genomic context of the interacting partners; 47 evidences based on co-expression; 37 evidences from high-throughput experiments; 584 evidences derived from the literature.

Gene Ontology annotation files for L. rhamnosus GG and C. jejuni NCTC 11168 were obtained from http://www.ebi.ac.uk/GOA/proteomes.html. To calculate which functional GO classes were enriched amongst interacting partners of a certain GT, we used the hypergeometric test, corrected for multiple testing using False Discovery Rate [32].

We then created ‘consensus networks’ that combine the local network neighborhood of all GTs, predicted to belong to the same specificity class and of which the local subnetworks are enriched in the same GO terms. GT-specific subnetworks were merged in a consensus network by retaining the edges from all the composing subnetworks that either reflect GT-GT interactions, interactions between a GT and one or more transmembrane proteins (membrane associations) or interactions between GTs and proteins with predicted glycosylation signals (predicted protein substrate relation).

Detection of putative protein glycosylation sites

Glycosylation sites were predicted in the proteomes of C. jejuni NCTC 11168 and L. rhamnosus GG using the GlycoPP webserver (http://www.imtech.res.in/raghava/glycopp/), specially developed for the analysis of prokaryotic protein sequences. Predictions were made using the hybrid approaches: BPP + ASA (for N-glycosites predictions) and PPP + ASA (for O-glycosites prediction) as suggested by the developers. A SVM threshold of 0.5 was used to reduce the probability of false positive predictions.

Prediction of transmembrane helices

Transmembrane helices were predicted using the TMHMM server version 2.0 (http://www.cbs.dtu.dk/services/TMHMM/).

Benchmark

The available data on glycosylation in the paradigm organism C. jejuni NCTC 11168 was used for benchmarking purposes and helped us to fine-tune and evaluate our workflow. C. jejuni is considered as a model for bacterial glycosylation, since it can not only N- and O- glycosylate proteins by both sequential and en bloc transfer [33, 34], but also produces a wide variety of glycoconjugates, including PG, LOS and CPS. Because glycosylation is extensively studied in C. jejuni NCTC 11168 we used this model system to compile a literature benchmark dataset. We obtained information on 10 proteins with experimentally verified glycosyltransferase activity and known substrate specificity in C. jejuni (Cj1124c, Cj1125c, Cj1126c, Cj1127c, Cj1128c and Cj1129c involved in protein N-glycosylation and Cj1133, Cj1136, Cj1139c and Cj1148 involved in LOS biosynthesis). Proteins annotated in C. jejuni NCTC 11168 as GTs based on indirect evidence (e.g. through homology assignment) were omitted from the benchmark dataset.

Reannotation of GTs in C. jejuni and L. rhamnosus GGbased on our predictions and literature

For the GTs that were previously annotated with a GT-related function, a simplified annotation is proposed when the evidence on the exact GT activity is not available for L. rhamnosus GG (such as for LGG_00279, LGG_00280 and LGG_00281). In addition, gene names inferred from non-strong homology searches (i.e. BLASTn E-value > 0.01) were removed (e.g. LGG_00348). For GTs putatively involved in polysaccharide biosynthesis (LGG_00279-LGG_00283, see below), gene names were corrected in agreement with the correct gene nomenclature [35].

Experimental work

L. rhamnosus GG and its mutant derivatives were grown in MRS without agitation. A new ΔwelI::TcR gene deletion mutant, lacking the LGG_02047 gene, termed CMPG 10811, was constructed as described earlier [36], using the pro-7946 (5′-ATACTAGTTCTTATCATAGTTTCCAGACC-3′) and pro-7947 (5′-ATCCCGGGGTGGGGAACTTGCTG-3′) primers. As this is a gene deletion mutant in an operon, polar effects can not completely be ruled out. Total EPS determination, monomer analysis and adhesion assays were performed as previously described [37]. Statistical analysis (One-way ANOVA) was performed using GraphPad Prism 6 on data corresponding to three technical repeats of three independent biological samples.

Results

Annotating putative glycosyltransferases

To predict additional GTs, we used an HMM based screening (Figure 1A). To maximize the sensitivity of our screening, the heterogeneous functional family of GTs was represented by a collection of 12 different HMMs, each of which captures a different characteristic of known GTs (Table 1). These 12 HMMs were subdivided into three groups depending on their expected specificity for GTs, referred to as respectively I) ‘Rossmann-fold domains’, II) ‘Sugar transferase’ and ‘UDP-Glycosyltranferase’ and III) a set of nine more GT-specific HMMs.

Figure 1
figure 1

Glycosyltransferase annotation flow. A: Genome-wide annotation of glycosyltransferases (GTs). Glycosyltransferases are predicted by scanning the proteomes of the studied species for GT-specific signatures using Hidden Markov Models (HMM) from SUPERFAMILY, CAZy and Pfam. An additional fold recognition filtering step is applied to only retain those genes containing a three-dimensional fold (inferred by the PGenTHREADER algorithm) with significant homology to folds present in experimentally confirmed GTs (deposited in the SCOP database). B: Predicting GT substrate class and putative mode of action (bottom panel). The local network neighborhood of each query GT (black node) in a functional interaction network (STRING) is used to extract a GT-specific local subnetwork for each query GT. The local subnetwork of a GT comprises predicted functional partners (proteins being functionally related to the query GT). Based on the GO enrichment analysis of these genes in this local subnetwork, the substrate class of the query GT is derived. To gain information on the mode of glycosylation, the GT specific local subnetwork is further annotated with either membrane associations between a query GT and a predicted transmembrane protein (blue edge) and with relations indicative for protein glycosylation (yellow edge).

As HMM-based screenings, definitely those performed with the least GT-specific HMMs, tend to also find many non-specific hits (false positives), predictions were further filtered using a protein fold recognition step: GTs predicted by the HMM profiling were only retained if they contained a three-dimensional fold with significant homology to folds present in experimentally confirmed GTs from any species (referred to as the refined set in Figure 1) (see Methods).

The results of the HMM based screening in both L. rhamnosus GG and C. jejuni NCTC 11168 before and after filtering with the fold based predictions are shown in Figure 2, together with the most abundant GO categories present amongst the predicted GTs. Filtering successfully reduced potential false positive predictions, for instance, a large fraction of oxidoreductases (all binding the cofactor NAD) obtained by screening with the least specific ‘Rossmann-fold domain’ HMM were removed after the fold recognition based filtering (Figure 2A). The three predictions in C. jejuni (Additional file 1: Table S1) and the five in L. rhamnosus GG (Table 2) made by the ‘Rossmann-fold domain’ HMM and retained after the fold recognition could not be retrieved by any of the other HMM models, showing the added value of also using this least specific class of HMMs. Screening with the ‘Sugar transferases’ and ‘UDP-glycosyltransferases’ HMMs in contrast resulted in predictions that were quite GT-specific, as indeed approximately 50% of the originally obtained predictions also contain a GT-like fold (Figure 2B and C). Fold-based filtering here removed mainly predicted DNA-binding proteins, as their mechanism of binding DNA is also based on recognizing the sugar moieties of the nucleotides. As expected, screening with the HMMs obtained from Pfam and CAZy resulted both in C. jejuni NCTC 11168 and L. rhamnosus GG in the highest fraction of hits that also displayed a GT-like fold (Figure 2D).

Figure 2
figure 2

Annotated glycosyltransferases. Results for the model system Campylobacter jejuni are shown on the left panel and for L. rhamnosus GG on the right panel. Putative GTs were predicted using an HMM based screening. A: results obtained with an HMM recognizing ‘Rossmann-fold domains’, expected to be the HMM with the lowest specificity towards GTs (Table 1, class I). B and C: results obtained with a family of HMMs of intermediate specificity for GTs (Table 1, class II). D: results obtained with the class of HMMs, most specific for GTs (Table 1, class III). Pie charts indicate the extent to which different functional classes were enriched amongst the predictions obtained with the respective classes of HMMs. Slices indicated in red on the pie chart correspond to the functional classes of the predictions that were retained after the fold recognition filtering step. For each group of HMMs, the total number of predictions is denoted in black on top of every pie chart and the number of predictions retained after applying the fold recognition step is denoted in red.

Table 2 Updated annotation of glycosyltransferases predicted in the genome of Lactobacillus rhamnosus GG

The performance of our GT prediction flow with and without the fold recognition filtering step was also evaluated in terms of the true-positive rate on the C. jejuni benchmark (containing 10 proteins with experimentally validated GT activity in C. jejuni NCTC 11168, see Methods). To obtain a full recall of 100% (that is retrieving all 10 positives), we had to make 184 predictions before the filtering. After the filtering the true positive rate increased from 10/184 to 10/44 (Additional file 1: Table S1). In addition to recovering all benchmark GTs (those indicated with experimental validation in Additional file 1: Table S1), most other predictions corresponded to previously made GT related annotations in C. jejuni NCTC 11168 that were based on indirect evidence (e.g. through experimental validation in other closely related species), such as the loci comprising the GT genes responsible for the synthesis of LOS (CJ1133 – CJ1148) [38], the GTs for N- (CJ1121c– CJ1129c) [33] and O-glycoprotein biosynthesis (CJ1311CJ1333) [34] and the CPS biosynthesis cluster (CJ1416cCJ1442c) [39, 40]. In addition, we made a total of 17 new predictions for yet unannotated genes in C. jejuni NCTC 11168 (Additional file 1: Table S1). Finally, we also retrieved four potential false positives (Additional file 1: Table S1).

The good agreement between our predictions and known information on glycosylation in C. jejuni NCTC 11168 [33], suggests that also for L. rhamnosus GG, the predictions summarized in Table 2 reflect true GTs. In addition, Table 2 provides a curated annotation update of GTs in L. rhamnosus GG: besides adding novel predictions, we removed potential erroneous annotations that originated through homology-based associations (indicated by conservation in Additional file 1: Table S1) as especially for GTs it is difficult to extrapolate the functional annotation without further experimental evidence (e.g. for LGG_00279). For GTs putatively involved in polysaccharide biosynthesis (LGG_00279-LGG_00283, see below), gene names were corrected in agreement with the conventional gene nomenclature [35].

Of the total number of 48 final predictions in L. rhamnosus GG (Table 2), five correspond to the experimentally documented locus encoding the enzymes involved in the synthesis of the complex galactose-rich EPS of L. rhamnosus GG [37, 41]. We also recovered the conserved cluster of GTs involved in the production of the intracellular storage glycogen-like polysaccharides [42] and the GTs necessary for the biosynthesis of PG [17]. In 33 cases, our predictions were consistent with previously annotated GTs (supported either by sequence conservation or by experimental evidence in related species. In five cases, indicated in Table 2 with a hash, our predictions are likely false positives. Eight of the 48 predicted GTs in L. rhamnosus GG were completely novel (indicated with a star in Table 2).

Among the novel predictions, two resulted from the screening with the ‘Rossmann-fold domain’ (class I) (LGG_01412 and LGG_00928, see Table 2). The other novel predictions LGG_01195 (previously annotated as ‘ABC transporter’), LGG_00985 (previously annotated as ‘integral membrane protein’) and LGG_02347 (previously annotated as ‘hypothetical protein’ were all detected by screening with the dedicated HMMs of class III (Table 1), further confirming the added value of these HMMs to find additional GTs. The screening with the HMMs of class II predicted as potential GTs LGG_00283 (a yet unannotated protein), LGG_01991 and LGG_01992. Both latter enzymes exhibit a high similarity with experimentally validated GTs in E. coli of the UDP-glycosyltransferase/Glycogen phosphorylase superfamily [42], further confirming their GT activity. However, they also show high sequence homology with UDP-N-acetylglucosamine 2-epimerases. This would be in agreement with the work of Campbell et al. (2000) showing that UDP-N-acetylglucosamine 2-epimerase has homology to phosphoglycosyl transferases and shares the same catalytic mechanism [43].

Despite the similar number of predicted GTs, the genomic organization of these predicted GTs is very different in C. jejuni NCTC 11168 and L. rhamnosus GG. In C. jejuni NCTC 11168, about 82% of the predicted GTs (corresponding to 36 GTs) are clustered into seven genomic regions, each of which contains at least two and on average five GTs that are physically located next to each other. The remaining eight predicted C. jejuni NCTC 11168 GTs are scattered in the genome (i.e. with no other GT present immediately up- or downstream). For L. rhamnosus GG, a smaller fraction of the predicted GTs is organized in clusters: about 56% of the predicted GTs (corresponding to 28 GTs) are located in 9 clusters, that are on average slightly smaller (with a mean size of three GTs) than those found in C. jejuni NCTC 11168. The remaining 20 predicted GTs in L. rhamnosus GG are isolated in the genome. For both species, most of the well-studied experimentally verified GTs are localized in these clusters, e.g. in C. jejuni NCTC 11168 these clusters correspond to the genomic regions involved in the synthesis of LOS, CPS and N- and in O-protein glycosylation [4], whereas in L. rhamnosus GG one of the predicted clusters correspond to the known region for galactose-rich EPS [37, 41] and one to the cluster for the biosynthesis of intracellular storage glycogen-like polysaccharides [42, 44]. The function of the remaining seven clusters in L. rhamnosus GG is yet unknown.

Compared to the ones organized in clusters in both genomes, most of the GTs found in isolation appear to be much less studied. A closer inspection of these isolated GTs showed that in L. rhamnosus GG (in 7 of the 20 cases (LGG_01057, LGG_01069, LGG_01147, LGG_01412, LGG_01487, LGG_01538, LGG_02004)), but not in C. jejuni NCTC 11168, these isolated GTs are flanked by DNA topoisomerases, tyrosine recombinases, Holliday junction-specific endonucleases, phage-related resolvases and transposases (according to the current genome annotation of L. rhamnosus GG (NC_013198.1)). In addition, overlaying our predictions with the results of a previous comparative analysis between L. rhamnosus GG and its close relative L. rhamnosus LC705 [44], indicates that many of the isolated GTs we identified are specific for L. rhamnosus GG (such as LGG_02004). These observations, together with the lower fraction of GTs occurring in large genomic clusters, indicates that in L. rhamnosus GG, much more than in C. jejuni NCTC 11168, the glycosylation potential has been shaped by horizontal gene transfer and intra-genomic rearrangements, similarly to what has been observed for GTs belonging to family 6 of GTs in bacteria and vertebrates (CAZy database) [45, 46].

Network-based strategy relating GTs to their substrate classes

To relate the predicted GTs to their potential substrates, we exploit the ‘local neighborhood’ of these GTs in a functional network, hereby assuming that GTs should be connected to their substrates, either directly or indirectly, via other GTs or enzymes. For the network, we relied on STRING, of which the functional interactions are inferred from physical (genome-wide protein-protein interactions, literature) and functional data (genomic co-localization, co-expression, co-occurrences, gene fusion-fission events) [14, 31]. The local neighborhood of a predicted GT (or local subnetwork) is here defined as the nodes that directly connect to the predicted GT (the latter of which is also referred to as the query GT) in the STRING network. We could derive 44 subnetworks for C. jejuni NCTC 11168, and 48 for L. rhamnosus GG. For each GT-specific subnetwork, the GO categories that were most overrepresented amongst the members of the subnetwork were used to infer for the query GT of each subnetwork a putative substrate class. As such we could predict a substrate class for 30/44 GTs in C. jejuni NCTC 11168 and for 20/48 GTs in L. rhamnosus GG which related to either saccharides, PG, proteins and lipids (see Additional file 2: Table S2 for C. jejuni NCTC 11168 and Table 3 for L. rhamnosus GG).

Table 3 Proposed substrate classess of predicted glycosyltransferases in Lactobacillus rhamnosus GG

The relation of the predicted GTs with their network neighbours was further specified using information on putative membrane associations or presence of glycosylation sites in the network members (Methods): a query-GT being connected to a transmembrane protein is referred to as a ‘membrane association’ and is indicative for soluble GTs that exert their action by interacting with transmembrane proteins, e.g. a transporter of glycoconjugates [4751]. A query-GT being connected to proteins with putative glycosylation sites hints towards the glycosylation of those proteins by the query-GT (substrate relation).

To gain insight in the mutual interactions between GTs and of these GTs with other proteins involved in the same process, we created ‘consensus networks’ that combine the local network neighbourhood of all GTs, predicted to belong to the same specificity class and of which the local subnetworks are enriched in the same GO terms (Figure 3).

Figure 3
figure 3

Consensus networks derived for each of the predicted substrate classes of putative GTs in L. rhamnosus GG. Consensus networks show all GTs, having the same substrate class, together with their protein neighbors that are hypothesized to contribute to the same common glycosylation mechanism as the one the GTs are involved in. On the consensus networks, nodes are proteins than can either be GTs (green nodes), transmembrane proteins (orange nodes) or proteins containing glycosylation signals (violet nodes). Membrane associations established between GTs and transmembrane proteins are represented by blue edges while predicted substrate relations between GT and proteins containing glycosylation signals are represented by yellow edges. Black edges refer to interactions between predicted GTs. If the local network neighborhood of GTs (local subnetwork) belonging to the same substrate class shows enrichment in more than one GO category (e.g. both the GO terms of EPS and glycogen biosynthesis), the consensus network is shown for each of the enriched GO categories. A: consensus networks involving GTs, predicted to glycosylate saccharides. Note that here two independent consensus networks were derived corresponding to respectively extracellular and intracellular PS biosynthesis. B: consensus network involving GTs, predicted to glycosylate peptidoglycan (PG). C: consensus network involving GTs, predicted to glycosylate lipids. D: consensus networks involving GTs, predicted to glycosylate proteins. Three independent consensus networks were derived corresponding to respectively cell cycle regulation, protein translation and DNA metabolic processes. Our analysis suggests substrate promiscuity for MurG, PBP1A, PBP1B and PBPA, all of which were predicted to be involved in the glycosylation of both peptidoglycan and proteins.

Inferred substrate classes of predicted GTs in the benchmark

To assess the extent to which our network-based approach was able to correctly infer substrate classes, we used as benchmark again the 10 GTs in C. jejuni NCTC 11168 for which also the substrate specificity is known (see Methods). Our strategy was able to recover the known substrate class of all 10 GTs (sensitivity of 100%) on a total of 31 predicted substrate classes for GTs in C. jejuni (true positive rate of 10/31).

Inferred substrate classes of predicted GTs in L. rhamnosusGG

The 20 GTs in L. rhamnosus GG for which we could predict their putative substrate class are summarized in Table 3.

GTs predicted to glycosylate saccharides

In L. rhamnosus GG, the substrate class saccharides (Figure 3A) comprises the largest number of GTs, which is to be expected as saccharides are the most common substrates for GTs [3]. The group of GTs that could be related to saccharides comprises two consensus networks: the first consensus network consists of GTs that, according to their GO annotation are involved in the biosynthesis of extracellular polysaccharides (WclC, WclB, WelE, WelG, WelH, WelI, RmlA2, LGG_00295) [37, 41]. The topology of this consensus network is indicative for en bloc glycosylation [4, 5] because it contains several interconnected soluble GTs, all linked to a membrane-bound priming GT together with Wzx flippases that transfer the subunits en bloc (see below).

This consensus network (Figure 3A) can be further subdivided into two cliques of interconnected GTs. The first clique (welI-welG-welH-rmlA1-rmlA2) contains genes involved in the synthesis of galactose-rich EPS, such as amongst others WelE (LGG_02043), the priming GT, with an experimentally verified substrate [37]. From the previously annotated gene cluster for galactose-rich EPS [37, 44], our analysis only missed welJ, annotated as alpha-1,3-galactosyltransferase (LGG_02048), as this gene was not predicted as a GT in our analysis. This gene does not appear to contain any signatures of the currently known HMMs for GTs and might represent a false negative of our analysis or an erroneous annotation in the current release of the L. rhamnosus GG genome NC_013198.1. This last hypothesis is supported by the small gene size of welJ, which would be atypical for a GT.

Regarding the second clique (wclC-LGG_00295-wclB), it contains genes for which the substrate specificity towards saccharides is known from homology-based extrapolation only. As we know from previous work that L. rhamnosus GG contains, besides its galactose-rich EPS also shorter, glucose-rich polysaccharides structures, we would hypothesize that this clique contains the missing genes for those glucose-rich polysaccharides structures [52]. The prediction of an independent Wzx flippase for each of the sets of interconnected GTs (cliques) (i.e. LGG_02049 for the galactose-rich clique and WclC and WclB for the clique putatively responsible for glucose-rich EPS synthesis), together with the known exquisite substrate specificity of Wzx flippases [53] further supports the hypothesis of each clique being responsible for the biosynthesis of another glycan type. Assuming that indeed the upper clique is involved in the synthesis of glucose-rich saccharide structures implies that the predicted link between WelE and this second clique (WclC, LGG_00295 and WclB) must be mere functional (i.e. not invoking a direct interaction), since knock-out experiments indicate that WelE is not the direct priming GT of the glucose-rich EPS structures [37].

The second consensus network (Figure 3A lower part, GlgA, GlgC, GlgD, GlgP, GalU) recapitulates all known members of the glycosylation system involved in glycogen synthesis except GlgB (LGG_02027), a conserved glycogen branching enzyme with transglycosylase activity, i.e. an enzyme that has both hydrolase and GT characteristics [54], which was not picked up by our HMM-based search step. From the predicted GTs in this network only GlgA, previously already known as a glycogen synthase, seems to be a genuine GT [42, 55]. For the other proteins GlgC, GlgD and GlgP, GalU -though related to glycan biosynthesis- enzyme activities other than GT activity have been documented [56]. The consensus network of the glycogen enzymes is composed solely of soluble proteins, which is in agreement with the intracellular nature of the glycogen-like polysaccharides. The connectivity between only soluble GTs points towards a sequential glycosylation mechanism in which sugar monomers are directly transferred from activated sugar-nucleotide donors (probably produced by GalU) to the respective substrates.

GTs predicted to glycosylate peptidoglycans

Five GTs could be related to PG precursors (PBP1A, PBP1B, PBP2A, MurG and LGG_01538), an annotation that has previously been suggested based on sequence conservation of these GTs across species (Figure 3B). GO enrichment analysis of their functional subnetworks suggests, both in L. rhamnosus GG (Table 3) and C. jejuni NCTC 11168 (Additional file 2: Table S2), a link between PG biosynthesis and a diverse set of processes, such as the regulation of cell shape, cell cycle and response to antibiotics, in agreement with the well-known functions of PG. Compared to the genes involved in EPS biosynthesis, it is remarkable that the GT genes involved in PG biosynthesis and remodelling do not occur in genomic clusters. The diversity of the processes in which these PG GTs are involved, might imply their necessity to be expressed under different environmental stimuli, which in turn can explain their organization in individual transcriptional units rather than in operons.

The consensus network of this class of GTs (Figure 3B) shows that all of these GTs are predicted to have transmembrane domains except for the soluble protein encoded by murG. The network organization is consistent with the known two-stage mechanism of bacterial PG biosynthesis consisting of cytoplasmic glycosylation reactions mediated by soluble GTs, followed by membrane-bound transglycosylation activities [57, 58].

GTs predicted to glycosylate lipids

The group of GTs that could be related to lipids contains three predicted GTs (LGG_00998, LGG_00999, LGG_01057) (Figure 3C). For these three GTs, their respective functional subnetworks showed enrichment for the terms ‘carbohydrate’ and ‘lipid metabolism’, suggesting that they are involved in the synthesis of lipoglycans present on the cell wall of the Gram-positive bacterium L. rhamnosus GG. This predicted role is more plausible than their homology based annotated role as ‘LPS biosynthesis glycosyltransferases’, as LPS molecules are absent in Gram-positives. The sparsity of the consensus network of these three GTs might be due to the incompleteness of the STRING network. So far, the existence of lipoglycans in L. rhamnosus GG has not yet been shown by biochemical studies.

GTs predicted to glycosylate proteins

A final group of seven GTs could be related to protein substrates and contains both predicted transmembrane (PBP1A, PBP1B, PBP2A) and predicted soluble GTs (LGG_00825, LGG_00826, LGG_01147, MurG). The GTs in this class were classified as protein GTs because the putative protein substrates in their subnetworks carry glycosylation signals. The GTs fall apart in three consensus subnetworks related to respectively cell cycle regulation, protein translation and DNA metabolic processes (Figure 3D).

A first consensus network comprises three transmembrane GTs (PBP1A, PBP2A, PBP1B) and MurG all predicted to be involved in ‘cell cycle regulation’ (according to the GO enrichment analysis of their respective subnetworks). Their consensus network points towards a substrate relation between each of the four GTs MurG, PBP1A, PBP2A and PBP1B, and cell division proteins (between MurG, PBP1A, PBP2A and PBP1B and the cell division protein FtsI on the one hand and between PBP1A, LGG_01706 and LGG_00254 on the other hand). Two previous studies further support our predictions: in Bacteroides fragilis FtsI, and other cell cycle related proteins such as FtsX and FtsQ, have been shown to be glycosylated [59]. In addition, a very recent study in L. plantarum WCFS1 [60] provides experimental evidence for the glycosylation of the cell division proteins FtsY, FtsZ, and FtsK 1 [60]. Our results – on the other hand- indicate that the three transmembrane GTs and MurG, known to be involved in PG biosynthesis show substrate promiscuity and would also have relations with protein substrates in L. rhamnosus GG (Table 3). A link between PG biosynthesis and protein glycosylation is not completely impossible given the fact that these predicted ‘promiscuous’ GTs co-occur with their predicted protein substrates including FtsI in cell division multi-enzyme complexes (Figure 4).

Figure 4
figure 4

Protein glycosylation of the cell division machinery. Schematic overview of the cell division machinery of L. rhamnosus. PBP1A, PBP1B, PBPB2A and MurG are predicted to be putative GTs. Our network-based analysis predicted PBP3, FtsI and PBP2B as putative substrates of the indicated GTs. The Msp1 cell wall hydrolase is the experimentally validated glycoprotein in L. rhamnosus GG [36].

This link between PG biosynthesis and protein glycosylation is further supported by the fact that the other predicted protein substrate of PBP1A (the D-alanyl-D-alanine carboxypeptidase (LGG_00254)), is also known to be directly involved in PG biosynthesis by introducing interpeptide cross-links. Although not yet reported for D,D trans-peptidases, other PG remodeling enzymes such as the PG hydrolases Msp1 in L. rhamnosus GG [36] and Acm2 in L. plantarum WCFS1 [61] were recently shown to be glycosylated [62].

A second consensus cluster is composed of two soluble GTs predicted to be involved in ‘protein translation’ (LGG_00825-LGG_00826). Both of these GTs were predicted to participate in the glycosylation of YkuJ, a protein co-translated with CcpC, a repressor of the tricarboxylic acid cycle in Bacillus subtillis (Figure 3C) [63]. LGG_00825 and LGG_00826 also exhibit a membrane association mediated by LGG_00751, annotated in L. rhamnosus GG as a hypothetical protein with a pfam09335 domain typical for SNARE associated Golgi proteins in eukaryotes. The membrane association of both GTs via a protein involved in translation, together with the fact that the subnetwork of LGG_00825 is enriched in the function ‘protein translation’ is consistent with the existence of an eukaryotic counterpart of sequential co-translational glycosylation in bacteria [51].

A last consensus network comprises only one GT, LGG_01147, predicted to be involved in ‘DNA metabolic processes’. LGG_01147 shows a substrate relation with LGG_01145, encoding a putative DNA entry nuclease, while establishing a membrane association mediated by LGG_01146 (Figure 3C). Little is known about these interacting partners, but nucleases are often glycosylated in eukaryotes [64]. Although not specifically related to nucleases, glycosylation of extracellular enzymes has been reported in prokaryotes [36, 61, 6567] and is thought to promote their stability [36]. Whether this is also the case in LGG_01146 needs to be further substantiated.

Experimental analysis of the GT network for EPS biosynthesis

We experimentally validated the GT network hierarchy within the clique for galactose-rich EPS (Figure 3A) by constructing a gene deletion mutant in the welI gene and comparing its phenotype to the phenotypes of the wild type (WT) and the gene deletion mutant of the priming GT WelE. As phenotypes, we tested the amount and monomer composition of EPS, and the adhesion capacity to the intestinal epithelial cell line Caco-2 as an indirect measurement of the EPS level [37]. According to our predictions, WelI would be one of the GTs that transfer sugar moieties to the sugar subunit initiated by the priming GT WelE. Based on these predictions, a gene deletion mutant of WelI would be expected to affect the amount of EPS, as in the absence of WelI less sugar moieties will be transferred to the subunit initiated by the WelE, but the effect of the WelI deletion on the phenotype should be less severe than the effect observed when deleting the priming GT WelE. A phenotype for the welI mutant intermediate between the WT and the welE gene deletion mutant is indeed observed for both assays confirming the predicted role of WelI upstream of WelE: the ΔwelI::TcR mutant displays a lower galactose-rich EPS content than the WT, but a higher content and more galactose than the gene deletion mutant of the priming GT WelE (Figure 5A and B). In agreement with EPS having a negative effect on adhesion, the adherence capacity is the highest for the welE mutant, intermediate for the welI mutant and lowest for the WT (Figure 5C).

Figure 5
figure 5

Experimental validation of the EPS network hierarchy. A: Total cell wall polysaccharides were extracted from respectively LGG wild-type, a ΔwelE::TcR gene deletion mutant (CMPG5351) and ΔwelI::TcR gene deletion mutant (CMPG10811). The total amount of EPS was measured. Error bars indicate standard deviations (of three repeats). One-way ANOVA statistical analysis rendered a p-value smaller than 0.05 for the variation of EPS across strains. B: Sugar monomer composition. The data are expressed as relative amounts, taking the total amount of detected monomeric sugars as 100%. Error bars indicate standard deviations (of three repeats). One-way ANOVA analyses (performed independently on each of the three datasets) rendered significant p-values (<0.05) for the variation of each sugar monomer across strains. C: Adhesion capacity. The adhesion capacity of wild type and mutants to Caco-2 cells is compared. Error bars indicate standard deviations (of three repeats). A One-way ANOVA analysis rendered a significant p-value (<0.05) for the variation of the adhesion capacity of the strains.

Discussion

In this work we developed an analysis flow that uses sequence-based strategies to predict novel GTs, but also exploits a network-based approach to infer the substrate classes of these putative GT. Using a broad definition of GT activity, including also HMMs for OSTs and other non-typical GTs, allowed covering a large part of the glycosylation potential. Applying our flow resulted in a careful revision of GTs in the current genome annotation of L. rhamnosus GG (NC_013198.1). We confirmed the identity of 33 GTs and predicted 8 novel ones. In contrast to what is observed in C. jejuni NCTC 11168, GTs appear to be much less clustered in genomic regions, but rather occur as isolated genes flanked by transposable elements. This points towards a key role of horizontal gene transfer in the acquisition of the glycosylation potential of L. rhamnosus GG.

Complementing the sequence-based with a network based-approach allowed us to also relate some of those GTs to their potential substrates. Most prior experimental studies focused on analyzing the specificity of GTs organized in clusters together with their auxiliary enzymes, as this allows for the straightforward extrapolation of known specificities of some members to all members in the cluster. By considering, next to the genomic organization, also links in a functional network, we could predict the substrate classes for the numerous, isolated GTs in L. rhamnosus GG. Exploiting membrane associations and substrate relations for the nodes in the GT-centered networks helped predicting the mutual relations between the GTs and between the GTs and their substrates.

Our analysis contributed to the annotation of GTs in L. rhamnosus GG. For instance, we hypothesize that one of the genomic regions that was previously annotated to be involved in EPS biosynthesis in general would contain the missing genes involved in the biosynthesis of short glucose-rich polysaccharides that are known to decorate the surface of L. rhamnosus GG [68]. In addition, we uncovered several novel interactions. For instance, for the isolated GTs known to be involved in PG biosynthesis (PBP1B, PBP2A, PBP1A and MurG), our network-based approach suggests an additional role in the glycosylation of proteins that are either involved in the biosynthesis of the PG (LGG_00254) or in cell division (LGG_01280 or FtsI). Substrate promiscuity of GTs is not uncommon in bacteria as for instance in Gram-negative pathogens, enzymes with relaxed specificity are shared between different processes, such as LPS and glycoprotein biosynthesis [4, 38]. Validating the activity of GTs that were predicted to glycosylate proteins- is cumbersome, as in vitro enzymatic assays do not represent the cellular conditions that are relevant for the assembly of these GTs in multi-enzyme membrane-associated complexes [69]. However, because PG biosynthesis is a process involving multi-enzyme complexes for which the assembly is tightly regulated [69], it is not unlikely that also protein glycosylation would act as an additional regulatory layer in this structural complex formation. Provided our hypothesis on their substrate specificity towards both proteins and PG would be true, these promiscuous GTs (PBP1B, PBP2A, PBP1A and MurG) are unlikely to be the priming GTs of their putative protein substrates, given their well characterized specificities towards PG precursors in both Gram-positives and negatives [70]. We hypothesize that the priming GTs predicted to be involved in protein glycosylation must be (Lactobacillus) species- or strain-specific rather than generally conserved in prokaryotes. This is supported by the observation that the best documented glycoprotein in L. rhamnosus GG, i.e. Msp1, another protein associated to the divisome [36] (see Figure 4), was no longer glycosylated after transfer to the Gram-negative E. coli[15] despite the fact that E. coli also has PBP1A, PBP1B, PBP2A and MurG homologs. In addition, the sugar monomers added on Msp1 [36] and related PG hydrolases such as Acm2 [71] show different sugar lectin specificities in L. rhamnosus, L. casei and L. plantarum.

Conclusions

Our results show how combining sequence- and network-based computational predictions can unveil insights in the bacterial glycosylation potential, thereby providing novel links and interesting hypotheses for further investigation.

Authors’ information

The authors wish it to be known that, in their opinion, Aminael Sánchez-Rodríguez and Hanne L.P. Tytgat should be regarded as joint first authors.

References

  1. Kay E, Lesk VI, Tamaddoni-Nezhad A, Hitchen PG, Dell A, Sternberg MJ, Muggleton S, Wren BW: Systems analysis of bacterial glycomes. Biochem Soc Trans. 2010, 38 (5): 1290-1293. 10.1042/BST0381290.

    Article  CAS  PubMed  Google Scholar 

  2. Upreti RK, Kumar M, Shankar V: Bacterial glycoproteins: functions, biosynthesis and applications. Proteomics. 2003, 3 (4): 363-379. 10.1002/pmic.200390052.

    Article  CAS  PubMed  Google Scholar 

  3. Lairson LL, Henrissat B, Davies GJ, Withers SG: Glycosyltransferases: structures, functions, and mechanisms. Annu Rev Biochem. 2008, 77: 521-555. 10.1146/annurev.biochem.76.061005.092322.

    Article  CAS  PubMed  Google Scholar 

  4. Guerry P, Szymanski CM: Campylobacter sugars sticking out. Trends Microbiol. 2008, 16 (9): 428-435. 10.1016/j.tim.2008.07.002.

    Article  CAS  PubMed  Google Scholar 

  5. Hug I, Feldman MF: Analogies and homologies in lipopolysaccharide and glycoprotein biosynthesis in bacteria. Glycobiology. 2011, 21 (2): 138-151. 10.1093/glycob/cwq148.

    Article  CAS  PubMed  Google Scholar 

  6. Typas A, Banzhaf M, Gross CA, Vollmer W: From the regulation of peptidoglycan synthesis to bacterial growth and morphology. Nat Rev Microbiol. 2012, 10 (2): 123-136.

    CAS  Google Scholar 

  7. Lerouge I, Vanderleyden J: O-antigen structural variation: mechanisms and possible roles in animal/plant-microbe interactions. FEMS Microbiol Rev. 2002, 26 (1): 17-47. 10.1111/j.1574-6976.2002.tb00597.x.

    Article  CAS  PubMed  Google Scholar 

  8. Hansen SF, Bettler E, Rinnan A, Engelsen SB, Breton C: Exploring genomes for glycosyltransferases. Mol Biosyst. 2010, 6 (10): 1773-1781. 10.1039/c000238k.

    Article  CAS  PubMed  Google Scholar 

  9. Hansen SF, Bettler E, Wimmerova M, Imberty A, Lerouxel O, Breton C: Combination of several bioinformatics approaches for the identification of new putative glycosyltransferases in Arabidopsis. J Proteome Res. 2009, 8 (2): 743-753. 10.1021/pr800808m.

    Article  PubMed  Google Scholar 

  10. Egelund J, Skjot M, Geshi N, Ulvskov P, Petersen BL: A complementary bioinformatics approach to identify potential plant cell wall glycosyltransferase-encoding genes. Plant Physiol. 2004, 136 (1): 2609-2620. 10.1104/pp.104.042978.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Weerapana E, Imperiali B: Asparagine-linked protein glycosylation: from eukaryotic to prokaryotic systems. Glycobiology. 2006, 16 (6): 91R-101R. 10.1093/glycob/cwj099.

    Article  CAS  PubMed  Google Scholar 

  12. Faridmoayer A, Fentabil MA, Haurat MF, Yi W, Woodward R, Wang PG, Feldman MF: Extreme substrate promiscuity of the Neisseria oligosaccharyl transferase involved in protein O-glycosylation. J Biol Chem. 2008, 283 (50): 34596-34604. 10.1074/jbc.M807113200.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Wacker M, Feldman MF, Callewaert N, Kowarik M, Clarke BR, Pohl NL, Hernandez M, Vines ED, Valvano MA, Whitfield C, Aebi M: Substrate specificity of bacterial oligosaccharyltransferase suggests a common transfer mechanism for the bacterial and eukaryotic systems. Proc Natl Acad Sci U S A. 2006, 103 (18): 7088-7093. 10.1073/pnas.0509207103.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, von Mering C: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011, 39 (Database issue): D561-D568.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Claes IJ, Schoofs G, Regulski K, Courtin P, Chapot-Chartier MP, Rolain T, Hols P, von Ossowski I, Reunanen J, de Vos WM, Palva A, Vanderleyden J, De Keersmaecker SC, Lebeer S: Genetic and biochemical characterization of the cell wall hydrolase activity of the major secreted protein of Lactobacillus rhamnosus GG. PLoS One. 2012, 7 (2): e31588-10.1371/journal.pone.0031588.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Ha S, Gross B, Walker S: E. Coli MurG: a paradigm for a superfamily of glycosyltransferases. Curr Drug Targets Infect Disord. 2001, 1 (2): 201-213. 10.2174/1568005014606116.

    Article  CAS  PubMed  Google Scholar 

  17. Di Guilmi AM, Dessen A, Dideberg O, Vernet T: The glycosyltransferase domain of penicillin-binding protein 2a from Streptococcus pneumoniae catalyzes the polymerization of murein glycan chains. J Bacteriol. 2003, 185 (15): 4418-4423. 10.1128/JB.185.15.4418-4423.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Maldonado-Barragan A, Caballero-Guerrero B, Lucena-Padros H, Ruiz-Barba JL: Genome sequence of Lactobacillus pentosus IG1, a strain isolated from Spanish-style green olive fermentations. J Bacteriol. 2011, 193 (19): 5605-10.1128/JB.05736-11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Yoshida Y, Nakano Y, Yamashita Y, Koga T: Identification of a genetic locus essential for serotype b-specific antigen synthesis in Actinobacillus actinomycetemcomitans. Infect Immun. 1998, 66 (1): 107-114.

    CAS  PubMed Central  PubMed  Google Scholar 

  20. Sun Y, Wang M, Wang Q, Cao B, He X, Li K, Feng L, Wang L: Genetic analysis of the Cronobacter sakazakii O4 to O7 O-antigen gene clusters and development of a PCR assay for identification of all C. sakazakii O serotypes. Appl Environ Microbiol. 2012, 78 (11): 3966-3974. 10.1128/AEM.07825-11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Provencher C, LaPointe G, Sirois S, Van Calsteren MR, Roy D: Consensus-degenerate hybrid oligonucleotide primers for amplification of priming glycosyltransferase genes of the exopolysaccharide locus in strains of the Lactobacillus casei group. Appl Environ Microbiol. 2003, 69 (6): 3299-3307. 10.1128/AEM.69.6.3299-3307.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Baiet B, Burel C, Saint-Jean B, Louvet R, Menu-Bouaouiche L, Kiefer-Meyer MC, Mathieu-Rivet E, Lefebvre T, Castel H, Carlier A, Cadoret JP, Lerouge P, Bardor M: N-glycans of Phaeodactylum tricornutum diatom and functional characterization of its N-acetylglucosaminyltransferase I enzyme. J Biol Chem. 2011, 286 (8): 6152-6164. 10.1074/jbc.M110.175711.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  23. Knauer R, Lehle L: The oligosaccharyltransferase complex from Saccharomyces cerevisiae. Isolation of the OST6 gene, its synthetic interaction with OST3, and analysis of the native complex. J Biol Chem. 1999, 274 (24): 17249-17256. 10.1074/jbc.274.24.17249.

    Article  CAS  PubMed  Google Scholar 

  24. Silberstein S, Collins PG, Kelleher DJ, Gilmore R: The essential OST2 gene encodes the 16-kD subunit of the yeast oligosaccharyltransferase, a highly conserved protein expressed in diverse eukaryotic organisms. J Cell Biol. 1995, 131 (2): 371-383. 10.1083/jcb.131.2.371.

    Article  CAS  PubMed  Google Scholar 

  25. Campbell JA, Davies GJ, Bulone V, Henrissat B: A classification of nucleotide-diphospho-sugar glycosyltransferases based on amino acid sequence similarities. Biochem J. 1997, 326 (Pt 3): 929-939.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Mengin-Lecreulx D, Texier L, Rousseau M, van Heijenoort J: The murG gene of Escherichia coli codes for the UDP-N-acetylglucosamine: N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase involved in the membrane steps of peptidoglycan synthesis. J Bacteriol. 1991, 173 (15): 4625-4636.

    CAS  PubMed Central  PubMed  Google Scholar 

  27. Coutinho PM, Deleury E, Davies GJ, Henrissat B: An evolving hierarchical family classification for glycosyltransferases. J Mol Biol. 2003, 328 (2): 307-317. 10.1016/S0022-2836(03)00307-3.

    Article  CAS  PubMed  Google Scholar 

  28. Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (Database issue): D211-D222.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39 (Web Server issue): W29-W37.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Lobley A, Sadowski MI, Jones DT: pGenTHREADER and pDomTHREADER: new methods for improved protein fold recognition and superfamily discrimination. Bioinformatics. 2009, 25 (14): 1761-1767. 10.1093/bioinformatics/btp302.

    Article  CAS  PubMed  Google Scholar 

  31. von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA, Bork P: STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005, 33 (Database issue): D433-D437.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Nothaft H, Szymanski CM: Protein glycosylation in bacteria: sweeter than ever. Nat Rev Microbiol. 2010, 8 (11): 765-778. 10.1038/nrmicro2383.

    Article  CAS  PubMed  Google Scholar 

  34. Szymanski CM, Logan SM, Linton D, Wren BW: Campylobacter–a tale of two protein glycosylation systems. Trends Microbiol. 2003, 11 (5): 233-238. 10.1016/S0966-842X(03)00079-9.

    Article  CAS  PubMed  Google Scholar 

  35. Reeves PR, Hobbs M, Valvano MA, Skurnik M, Whitfield C, Coplin D, Kido N, Klena J, Maskell D, Raetz CR, Rick PD: Bacterial polysaccharide synthesis and gene nomenclature. Trends Microbiol. 1996, 4 (12): 495-503. 10.1016/S0966-842X(97)82912-5.

    Article  CAS  PubMed  Google Scholar 

  36. Lebeer S, Claes IJ, Balog CI, Schoofs G, Verhoeven TL, Nys K, von Ossowski I, de Vos WM, Tytgat HL, Agostinis P, Palva A, Van Damme EJ, Deelder AM, De Keersmaecker SC, Wuhrer M, Vanderleyden J: The major secreted protein Msp1/p75 is O-glycosylated in Lactobacillus rhamnosus GG. Microb Cell Fact. 2012, 11: 15-10.1186/1475-2859-11-15.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Lebeer S, Verhoeven TL, Francius G, Schoofs G, Lambrichts I, Dufrene Y, Vanderleyden J, De Keersmaecker SC: Identification of a gene cluster for the Biosynthesis of a long, Galactose-Rich Exopolysaccharide in Lactobacillus rhamnosus GG and functional analysis of the priming Glycosyltransferase. Appl Environ Microbiol. 2009, 75 (11): 3554-3563. 10.1128/AEM.02919-08.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Karlyshev AV, Ketley JM, Wren BW: The Campylobacter jejuni glycome. FEMS Microbiol Rev. 2005, 29 (2): 377-390.

    CAS  PubMed  Google Scholar 

  39. Linton D, Karlyshev AV, Wren BW: Deciphering Campylobacter jejuni cell surface interactions from the genome sequence. Curr Opin Microbiol. 2001, 4 (1): 35-40. 10.1016/S1369-5274(00)00161-2.

    Article  CAS  PubMed  Google Scholar 

  40. Gundogdu O, Bentley SD, Holden MT, Parkhill J, Dorrell N, Wren BW: Re-annotation and re-analysis of the Campylobacter jejuni NCTC11168 genome sequence. BMC Genomics. 2007, 8: 162-10.1186/1471-2164-8-162.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Lebeer S, Claes IJ, Verhoeven TL, Vanderleyden J, De Keersmaecker SC: Exopolysaccharides of Lactobacillus rhamnosus GG form a protective shield against innate immune factors in the intestine. Microb Biotechnol. 2011, 4 (3): 368-374. 10.1111/j.1751-7915.2010.00199.x.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Kiel JA, Boels JM, Beldman G, Venema G: Glycogen in Bacillus subtilis: molecular characterization of an operon encoding enzymes involved in glycogen biosynthesis and degradation. Mol Microbiol. 1994, 11 (1): 203-218. 10.1111/j.1365-2958.1994.tb00301.x.

    Article  CAS  PubMed  Google Scholar 

  43. Campbell RE, Mosimann SC, Tanner ME, Strynadka NC: The structure of UDP-N-acetylglucosamine 2-epimerase reveals homology to phosphoglycosyl transferases. Biochemistry. 2000, 39 (49): 14993-15001. 10.1021/bi001627x.

    Article  CAS  PubMed  Google Scholar 

  44. Kankainen M, Paulin L, Tynkkynen S, von Ossowski I, Reunanen J, Partanen P, Satokari R, Vesterlund S, Hendrickx APA, Lebeer S, De Keersmaecker SCJ, Vanderleyden J, Hamalainen T, Laukkanen S, Salovuori N, Ritari J, Alatalo E, Korpela R, Mattila-Sandholm T, Lassig A, Hatakka K, Kinnunen KT, Karjalainen H, Saxelin M, Laakso K, Surakka A, Palva A, Salusjarvi T, Auvinen P, de Vos WM: Comparative genomic analysis of Lactobacillus rhamnosus GG reveals pili containing a human-mucus binding protein. Proc Natl Acad Sci USA. 2009, 106 (40): 17193-17198. 10.1073/pnas.0908876106.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 2009, 37 (Database issue): D233-D238.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Brew K, Tumbale P, Acharya KR: Family 6 glycosyltransferases in vertebrates and bacteria: inactivation and horizontal gene transfer may enhance mutualism between vertebrates and bacteria. J Biol Chem. 2010, 285 (48): 37121-37127. 10.1074/jbc.R110.176248.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  47. Mohammadi T, Karczmarek A, Crouvoisier M, Bouhss A, Mengin-Lecreulx D, den Blaauwen T: The essential peptidoglycan glycosyltransferase MurG forms a complex with proteins involved in lateral envelope growth as well as with proteins involved in cell division in Escherichia coli. Mol Microbiol. 2007, 65 (4): 1106-1121. 10.1111/j.1365-2958.2007.05851.x.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Charbonneau ME, Cote JP, Haurat MF, Reiz B, Crepin S, Berthiaume F, Dozois CM, Feldman MF, Mourez M: A structural motif is the recognition site for a new family of bacterial protein O-glycosyltransferases. Mol Microbiol. 2012, 83 (5): 894-907. 10.1111/j.1365-2958.2012.07973.x.

    Article  CAS  PubMed  Google Scholar 

  49. Grass S, Lichti CF, Townsend RR, Gross J, St Geme JW: The Haemophilus influenzae HMW1C protein is a glycosyltransferase that transfers hexose residues to asparagine sites in the HMW1 adhesin. PLoS Pathog. 2010, 6 (5): e1000919-10.1371/journal.ppat.1000919.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Wu R, Wu H: A molecular chaperone mediates a two-protein enzyme complex and glycosylation of serine-rich streptococcal adhesins. J Biol Chem. 2011, 286 (40): 34923-34931. 10.1074/jbc.M111.239350.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Dell A, Galadari A, Sastre F, Hitchen P: Similarities and differences in the glycosylation mechanisms in prokaryotes and eukaryotes. Int J Microbiol. 2010, 2010: 148178-

    Article  PubMed Central  PubMed  Google Scholar 

  52. Francius G, Lebeer S, Alsteens D, Wildling L, Gruber HJ, Hols P, Keersmaecker SD, Vanderleyden J, Dufrêne YF: Detection, localization, and conformational analysis of single polysaccharide molecules on live bacteria. ACS Nano. 2008, 2 (9): 1921-1929. 10.1021/nn800341b.

    Article  CAS  PubMed  Google Scholar 

  53. Islam ST, Lam JS: Wzx flippase-mediated membrane translocation of sugar polymer precursors in bacteria. Environ Microbiol. 2013, 15 (4): 1001-1015. 10.1111/j.1462-2920.2012.02890.x.

    Article  CAS  PubMed  Google Scholar 

  54. Lim WJ, Park SR, Kim MK, An CL, Yun HJ, Hong SY, Kim EJ, Shin EC, Lee SW, Lim YP, Yun HD: Cloning and characterization of the glycogen branching enzyme gene existing in tandem with the glycogen debranching enzyme from Pectobacterium chrysanthemi PY35. Biochem Biophys Res Commun. 2003, 300 (1): 93-101. 10.1016/S0006-291X(02)02763-8.

    Article  CAS  PubMed  Google Scholar 

  55. Wilson WA, Roach PJ, Montero M, Baroja-Fernandez E, Munoz FJ, Eydallin G, Viale AM, Pozueta-Romero J: Regulation of glycogen metabolism in yeast and bacteria. FEMS Microbiol Rev. 2010, 34 (6): 952-985.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Ballicora MA, Iglesias AA, Preiss J: ADP-glucose pyrophosphorylase, a regulatory enzyme for bacterial glycogen synthesis. Microbiol Mol Biol Rev. 2003, 67 (2): 213-225. 10.1128/MMBR.67.2.213-225.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. van Heijenoort J: Formation of the glycan chains in the synthesis of bacterial peptidoglycan. Glycobiology. 2001, 11 (3): 25R-36R. 10.1093/glycob/11.3.25R.

    Article  CAS  PubMed  Google Scholar 

  58. Peregrin-Alvarez JM, Xiong X, Su C, Parkinson J: The modular organization of protein interactions in Escherichia coli. PLoS Comput Biol. 2009, 5 (10): e1000523-10.1371/journal.pcbi.1000523.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Fletcher CM, Coyne MJ, Comstock LE: Theoretical and experimental characterization of the scope of protein O-glycosylation in Bacteroides fragilis. J Biol Chem. 2011, 286 (5): 3219-3226. 10.1074/jbc.M110.194506.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  60. Fredriksen L, Moen A, Adzhubei AA, Mathiesen G, Eijsink VG, Egge-Jacobsen W: Lactobacillus plantarum WCFS1 O-linked protein glycosylation: An extended spectrum of target proteins and modification sites detected by mass spectrometry. Glycobiology. 2013, 23 (12): 1439-1451. 10.1093/glycob/cwt071.

    Article  CAS  PubMed  Google Scholar 

  61. Fredriksen L, Mathiesen G, Moen A, Bron PA, Kleerebezem M, Eijsink VG, Egge-Jacobsen W: The major autolysin Acm2 from Lactobacillus plantarum undergoes cytoplasmic O-glycosylation. J Bacteriol. 2012, 194 (2): 325-333. 10.1128/JB.06314-11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  62. Barinka C, Sacha P, Sklenar J, Man P, Bezouska K, Slusher BS, Konvalinka J: Identification of the N-glycosylation sites on glutamate carboxypeptidase II necessary for proteolytic activity. Protein Sci. 2004, 13 (6): 1627-1635. 10.1110/ps.04622104.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  63. Commichau FM, Forchhammer K, Stulke J: Regulatory links between carbon and nitrogen metabolism. Curr Opin Microbiol. 2006, 9 (2): 167-172. 10.1016/j.mib.2006.01.001.

    Article  CAS  PubMed  Google Scholar 

  64. Pimkin M, Miller CG, Blakesley L, Oleykowski CA, Kodali NS, Yeung AT: Characterization of a periplasmic S1-like nuclease coded by the Mesorhizobium loti symbiosis island. Biochem Biophys Res Commun. 2006, 343 (1): 77-84. 10.1016/j.bbrc.2006.02.117.

    Article  CAS  PubMed  Google Scholar 

  65. Brechtel E, Matuschek M, Hellberg A, Egelseer EM, Schmid R, Bahl H: Cell wall of Thermoanaerobacterium thermosulfurigenes EM1: isolation of its components and attachment of the xylanase XynA. Arch Microbiol. 1999, 171 (3): 159-165. 10.1007/s002030050694.

    Article  CAS  PubMed  Google Scholar 

  66. Huang L, Forsberg C, Thomas D: Purification and characterization of a chloride-stimulated cellobiosidase from Bacteroides succinogenes S85. J Bacteriol. 1988, 170 (7): 2923-2932.

    CAS  PubMed Central  PubMed  Google Scholar 

  67. Olsen O, Thomsen KK: Improvement of bacterial β-glucanase thermostability by glycosylation. J Gen Microbiol. 1991, 137 (3): 579-585. 10.1099/00221287-137-3-579.

    Article  CAS  Google Scholar 

  68. Prachi T, Beaussart A, Andre G, Rolain T, Lebeer S, Vanderleyden J, Hols P, Dufrene YF: Towards a nanoscale view of lactic acid bacteria. Micron. 2012, 43 (12): 1323-1330. 10.1016/j.micron.2012.01.001.

    Article  Google Scholar 

  69. Zapun A, Noirclerc-Savoye M, Helassa N, Vernet T: Peptidoglycan assembly machines: the biochemical evidence. Microb Drug Resist. 2012, 18 (3): 256-260. 10.1089/mdr.2011.0236.

    Article  CAS  PubMed  Google Scholar 

  70. Sauvage E, Kerff F, Terrak M, Ayala JA, Charlier P: The penicillin‒binding proteins: structure and role in peptidoglycan biosynthesis. FEMS Microbiol Rev. 2008, 32 (2): 234-258. 10.1111/j.1574-6976.2008.00105.x.

    Article  CAS  PubMed  Google Scholar 

  71. Rolain T, Bernard E, Beaussart A, Degand H, Courtin P, Egge-Jacobsen W, Bron PA, Morsomme P, Kleerebezem M, Chapot-Chartier M-P: O-glycosylation as a novel control mechanism of peptidoglycan hydrolase activity. J Biol Chem. 2013, 288 (31): 22233-22247. 10.1074/jbc.M113.470716.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Tine Verhoeven for her dedicated work in the genetic analysis and mutant construction; Geert Schoofs, Prof. C. Courtin and Jeroen Snelders for their help with respectively the EPS extraction, concentration determination and monomer analysis.

Funding

This work is supported by: 1) Katholieke Universiteit Leuven funding: GOA/08/011, project NATAR and PF 10/018; 2) Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT Vlaanderen): SBO-BioFrame, SBO-NEMOA, SB-Hanne Tytgat; 3) Fonds Wetenschappelijk Onderzoek-Vlaanderen (FWO) G.0428.13 N, Postdoc grant Sarah Lebeer and Research Grant 152011 N 4) Ghent University [Multidisciplinary Research Partnership “N2N”].

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Sarah Lebeer or Kathleen Marchal.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ASR conceived the study, performed the computational analysis and interpreted the results. HT performed the experimental work and provided biological input for the computational analysis. ASR and HT drafted the manuscript. KM participated in the design of the computational analysis and in the writing of the manuscript. SL participated in the design of the experimental work, provided biological input for the computational analyses and helped drafting the manuscript. JVDL, JW, SL and KM coordinated and managed the research. All authors contributed to writing the manuscript and approved its final version.

Aminael Sánchez-Rodríguez, Hanne LP Tytgat contributed equally to this work.

Electronic supplementary material

12864_2013_6075_MOESM1_ESM.pdf

Additional file 1: Table S1: List of glycosyltransferases predicted in the genome of Campylobacter jejuni NCTC 11168. Locus tag: gene identifier of the predicted GT. Genes for which a GT activity was predicted in this study that was not present in the current annotation are marked with a star (*). Potential false positive results are indicated with a hash (#). Current annotation: functional annotation as in the current genome release of GenBank (NC_002163.1). Proposed annotation: new annotation based on the results of our analysis. HMM: Description of the Hidden Markov Model (HMM) with which the indicated GT was identified. Note that all predicted GTs also passed the fold based filtering. Evidence: Type of evidence for the GT activity. Conservation: shows significant sequence conservation with an experimentally validated GT in a closely related species. Experimental validation: the GT activity has been experimentally validated in Campylobacter jejuni NCTC 11168. Reference: reference to the publication(s) supporting the prediction. (PDF 51 KB)

12864_2013_6075_MOESM2_ESM.pdf

Additional file 2: Table S2: Proposed substrate classes of glycosyltransferases in Campylobacter jejuni NCTC 11168. Locus tag: gene identifier of the predicted GT used as query in STRING to obtain a query-dependent subnetwork. Localization: indicates whether the query-GT was predicted to be cytoplasmic (C) or a transmembrane protein (TM). Enriched GO categories: GO categories enriched amongst the members of the query-dependent subnetwork of the indicated query-GT. Only categories showing an enrichment value p < 0.05 are shown (according to a hypergeometric test corrected for multiple testing using False Discovery Rate). Membrane association: it refers to edges between the query-GT and members of its subnetwork predicted to be transmembrane proteins. Partner GTs: predicted/experimentally validated GTs that belong to the subnetwork of the query-GT. Predicted substrate class of a query-GT: inferred from the GO enrichment analysis of the query-dependent subnetwork of the indicated query-GT derived from STRING. Potential protein substrate: it refers to edges between the query-GT and members of its subnetwork predicted to have N- or O-glycosylation signals. Such proteins are therefore suggested to be potential substrates of the query-GT in the cases where proteins are the proposed substrate. Evidence: level of evidence for the substrate class prediction. Conservation: shows a significant sequence conservation with a GT for which a susbtrate specificity has been experimentally validated in a closely related species. Experimental validation: the substrate specificity of the GT has been experimentally validated in Campylobacter jejuni NCTC 11168. Reference: publication(s) supporting the predicted substrate class of the query-GT. (PDF 131 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sánchez-Rodríguez, A., Tytgat, H.L., Winderickx, J. et al. A network-based approach to identify substrate classes of bacterial glycosyltransferases. BMC Genomics 15, 349 (2014). https://doi.org/10.1186/1471-2164-15-349

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-349

Keywords